Onset of patterns in an oscillated granular layer: continuum and molecular dynamics simulations

Time dependence and density inversion in simulations of vertically oscillated granular layers

We study a layer of grains atop a plate which oscillates sinusoidally in the direction of gravity, using three-dimensional, time-dependent numerical solutions of continuum equations to Navier-Stokes order as well as hard-sphere molecular dynamics simulations. For high accelerational amplitudes of the plate, the layer exhibits a steady-state "density inversion" in which a high-density portion of the layer is supported by a lower-density portion. At low accelerational amplitudes, the layer exhibits oscillatory time dependence that is strongly correlated to the motion of the plate. We show that continuum simulations yield results consistent with molecular dynamics results in both regimes.

Continuum simulations of shocks and patterns in vertically oscillated granular layers

We study interactions between shocks and standing-wave patterns in vertically oscillated layers of granular media using three-dimensional, time-dependent numerical solutions of continuum equations to Navier-Stokes order. We simulate a layer of grains atop a plate that oscillates sinusoidally in the direction of gravity. Standing waves form stripe patterns when the accelerational amplitude of the plate's oscillation exceeds a critical value. Shocks also form with each collision between the layer and the plate; we show that pressure gradients formed by these shocks cause the flow to reverse direction within the layer. This reversal leads to an oscillatory state of the pattern that is subharmonic with respect to the plate's oscillation. Finally, we study the relationship between shocks and patterns in layers oscillated at various frequencies and show that the pattern wavelength increases monotonically as the shock strength increases.

Subharmonic wave transition in a quasi-one-dimensional noisy fluidized shallow granular bed

We present an experimental and theoretical study of the pattern formation process of standing subharmonic waves in a fluidized quasi-one-dimensional shallow granular bed. The fluidization process is driven by means of a time-periodic air flow, analogous to a tapping type of forcing. Measurements of the amplitude of the critical mode close to the transition are in quite good agreement with those inferred from a universal stochastic amplitude equation. This allows us to determine both the bifurcation point of the deterministic system and the corresponding noise intensity. We also show that the probability density distribution is well described by a generalized Rayleigh distribution, which is the stationary solution of the corresponding Fokker-Planck equation of the universal stochastic amplitude equation that describes our system.

Effects of thermal noise on pattern onset in continuum simulations of shaken granular layers

The author investigates the onset of patterns in vertically oscillated layers of dissipative particles using numerical solutions of continuum equations to Navier-Stokes order. Above a critical accelerational amplitude of the cell, standing waves form stripe patterns which oscillate subharmonically with respect to the cell. Continuum simulations neglecting interparticle friction yield pattern wavelengths consistent with experiments using frictional particles. However, the critical acceleration for standing-wave formation is approximately 10% lower in continuum simulations without added noise than in molecular-dynamics simulations. This Brief Report incorporates fluctuating hydrodynamics theory into continuum simulations by adding noise terms with no fit parameters; this modification yields a critical acceleration in agreement with molecular-dynamics simulations.