On the relevance of numerical simulations to booming sand

We have performed a simulation study of 3D cohesionless granular flows down an inclined chute. We find that the oscillations observed in [L.E. Silbert, Phys. Rev. Lett., 94, 098002 (2005)] near the angle of repose are harmonic vibrations of the lowest normal mode. Their frequencies depend on the contact stiffness as well as on the depth of the flow. Could these oscillations account for the phenomena of "booming sand"? We estimate an effective contact stiffness from the Hertz law, but this leads to frequencies several times higher than observed. However, the Hertz law also predicts interpenetrations of a few nanometers, indicating that the oscillations frequencies are governed by the surface stiffness, which can be much lower than the bulk one. This is in agreement with previous studies ascribing the ability to sing to the presence of a soft coating on the grain surface.Comment: accepted for publication in Physical Review E http://pre.aps.org; Physical Review E (2012) to be publishe

Homogeneous cooling of rough, dissipative particles: Theory and simulations

We investigate freely cooling systems of rough spheres in two and three dimensions. Simulations using an event driven algorithm are compared with results of an approximate kinetic theory, based on the assumption of a generalized homogeneous cooling state. For short times $t$, translational and rotational energy are found to change linearly with $t$. For large times both energies decay like $t^{-2}$ with a ratio independent of time, but not corresponding to equipartition. Good agreement is found between theory and simulations, as long as no clustering instability is observed. System parameters, i.e. density, particle size, and particle mass can be absorbed in a rescaled time, so that the decay of translational and rotational energy is solely determined by normal restitution and surface roughness.Comment: 10 pages, 10 eps-figure

A quasi-elastic regime for vibrated granular gases

Using simple scaling arguments and two-dimensional numerical simulations of a granular gas excited by vibrating one of the container boundaries, we study a double limit of small $1-r$ and large $L$, where $r$ is the restitution coefficient and $L$ the size of the container. We show that if the particle density $n_0$ and $(1-r^2)(n_0 Ld)$ where $d$ is the particle diameter, are kept constant and small enough, the granular temperature, i.e. the mean value of the kinetic energy per particle, $/N$, tends to a constant whereas the mean dissipated power per particle, $/N$, decreases like $1/\sqrt{N}$ when $N$ increases, provided that $(1-r^2)(n_0 Ld)^2 < 1$. The relative fluctuations of $E$, $D$ and the power injected by the moving boundary, $I$, have simple properties in that regime. In addition, the granular temperature can be determined from the fluctuations of the power $I(t)$ injected by the moving boundary.

Energy flows in vibrated granular media

We study vibrated granular media, investigating each of the three components of the energy flow: particle-particle dissipation, energy input at the vibrating wall, and particle-wall dissipation. Energy dissipated by interparticle collisions is well estimated by existing theories when the granular material is dilute, and these theories are extended to include rotational kinetic energy. When the granular material is dense, the observed particle-particle dissipation rate decreases to as little as 2/5 of the theoretical prediction. We observe that the rate of energy input is the weight of the granular material times an average vibration velocity times a function of the ratio of particle to vibration velocity. `Particle-wall' dissipation has been neglected in all theories up to now, but can play an important role when the granular material is dilute. The ratio between gravitational potential energy and kinetic energy can vary by as much as a factor of 3. Previous simulations and experiments have shown that E ~ V^delta, with delta=2 for dilute granular material, and delta ~ 1.5 for dense granular material. We relate this change in exponent to the departure of particle-particle dissipation from its theoretical value.Comment: 19 pages revtex, 10 embedded eps figures, accepted by PR

Active galactic nucleus feedback in clusters of galaxies

Observations made during the last ten years with the Chandra X-ray Observatory have shed much light on the cooling gas in the centers of clusters of galaxies and the role of active galactic nucleus (AGN) heating. Cooling of the hot intracluster medium in cluster centers can feed the supermassive black holes found in the nuclei of the dominant cluster galaxies leading to AGN outbursts which can reheat the gas, suppressing cooling and large amounts of star formation. AGN heating can come in the form of shocks, buoyantly rising bubbles that have been inflated by radio lobes, and the dissipation of sound waves.Comment: Refereed review article published in Chandra's First Decade of Discovery Special Feature edition of the Proceedings of the National Academy of Science

Greater data science at baccalaureate institutions

Donoho's JCGS (in press) paper is a spirited call to action for statisticians, who he points out are losing ground in the field of data science by refusing to accept that data science is its own domain. (Or, at least, a domain that is becoming distinctly defined.) He calls on writings by John Tukey, Bill Cleveland, and Leo Breiman, among others, to remind us that statisticians have been dealing with data science for years, and encourages acceptance of the direction of the field while also ensuring that statistics is tightly integrated. As faculty at baccalaureate institutions (where the growth of undergraduate statistics programs has been dramatic), we are keen to ensure statistics has a place in data science and data science education. In his paper, Donoho is primarily focused on graduate education. At our undergraduate institutions, we are considering many of the same questions.Comment: in press response to Donoho paper in Journal of Computational Graphics and Statistic

Neural activation and functional connectivity during motor imagery of bimanual everyday actions

© 2012 Szameitat et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.Bimanual actions impose intermanual coordination demands not present during unimanual actions. We investigated the functional neuroanatomical correlates of these coordination demands in motor imagery (MI) of everyday actions using functional magnetic resonance imaging (fMRI). For this, 17 participants imagined unimanual actions with the left and right hand as well as bimanual actions while undergoing fMRI. A univariate fMRI analysis showed no reliable cortical activations specific to bimanual MI, indicating that intermanual coordination demands in MI are not associated with increased neural processing. A functional connectivity analysis based on psychophysiological interactions (PPI), however, revealed marked increases in connectivity between parietal and premotor areas within and between hemispheres. We conclude that in MI of everyday actions intermanual coordination demands are primarily met by changes in connectivity between areas and only moderately, if at all, by changes in the amount of neural activity. These results are the first characterization of the neuroanatomical correlates of bimanual coordination demands in MI. Our findings support the assumed equivalence of overt and imagined actions and highlight the differences between uni- and bimanual actions. The findings extent our understanding of the motor system and may aid the development of clinical neurorehabilitation approaches based on mental practice.This study was funded by the Medical Research Council, UK (CEG 61501; Dr Sterr)

The energy flux into a fluidized granular medium at a vibrating wall

We study the power input of a vibrating wall into a fluidized granular medium, using event driven simulations of a model granular system. The system consists of inelastic hard disks contained between a stationary and a vibrating elastic wall, in the absence of gravity. Two scaling relations for the power input are found, both involving the pressure. The transition between the two occurs when waves generated at the moving wall can propagate across the system. Choosing an appropriate waveform for the vibrating wall removes one of these scalings and renders the second very simple.Comment: 5 pages, revtex, 7 postscript figure