Dynamics and Stability of Shock Waves in Granular Gases Undergoing Activated Inelastic Collisions

Abstract

The present work investigates the dynamics and stability of shock waves in granular gases. The problem was modelled for a piston propagating into a system of disks that can undergo inelastic collisions if an impact threshold is exceeded. The model was addressed numerically at the microscopic and macroscopic levels. The molecular dynamics methodology employed the Event-Driven Molecular Dynamics method, and the continuum model was formulated using the Navier-Stokes equations for granular gases with the transport terms of Jenkins and Richman and a modified cooling rate term.

The inviscid steady state shock structure was derived and analyzed. The results indicated that a relaxing shock structure is expected for sufficiently strong shock waves. Beyond this limit the structure was shown to be independent of the initial energy, a finding similar to the strong shock approximation in molecular gases.

One-dimensional simulations demonstrated that the molecular dynamics and continuum models yield similar evolutions and structures of the shock wave, validating the continuum description of this study. Two-dimensional results showed that sufficiently strong shock waves can exhibit multi-dimensional instability with high density non-uniformities and convective rolls within the structure, with the size of instabilities shown to scale with the relaxation length of the shock structure. Instabilities were observed with the continuum description only with the inclusion of statistical fluctuations to density mimicking the molecular model. The cases that were unstable were shown to be in a regime whereby statistical fluctuations can become important, following the description for this regime by Bird.

Based on these findings, it is proposed that unstable shock behaviour can be observed for highly dissipative shock waves that yield short relaxation length scales, where fluctuations become important. The current work may shed light on unstable shock behaviour observed in dissipative gases, having implications for both granular media and molecular gases.