A nonlinear theory of non-stationary low Mach number channel flows of freely cooling nearly elastic granular gases

We use hydrodynamics to investigate non-stationary channel flows of freely cooling dilute granular gases. We focus on the regime where the sound travel time through the channel is much shorter than the characteristic cooling time of the gas. As a result, the gas pressure rapidly becomes almost homogeneous, while the typical Mach number of the flow drops well below unity. Eliminating the acoustic modes, we reduce the hydrodynamic equations to a single nonlinear and nonlocal equation of a reaction-diffusion type in Lagrangian coordinates. This equation describes a broad class of channel flows and, in particular, can follow the development of the clustering instability from a weakly perturbed homogeneous cooling state to strongly nonlinear states. If the heat diffusion is neglected, the reduced equation is exactly soluble, and the solution develops a finite-time density blowup. The heat diffusion, however, becomes important near the attempted singularity. It arrests the density blowup and brings about novel inhomogeneous cooling states (ICSs) of the gas, where the pressure continues to decay with time, while the density profile becomes time-independent. Both the density profile of an ICS, and the characteristic relaxation time towards it are determined by a single dimensionless parameter that describes the relative role of the inelastic energy loss and heat diffusion. At large values of this parameter, the intermediate cooling dynamics proceeds as a competition between low-density regions of the gas. This competition resembles Ostwald ripening: only one hole survives at the end.Comment: 20 pages, 15 figures, final versio

Simulation of an avalanche in a fluid with a soft-sphere / immersed boundary method including a lubrication force.

The present work aims at reproducing the local dynamics of a dense granular media evolving in a viscous fluid from the grain scale to that of thousands of grains, encountered in environmental multiphase flows. To this end a soft-sphere collision / immersed-boundary method is developed. The methods are validated alone through various standard configurations including static and dynamical situations. Then, simulations of binary wall-particle collisions in a fluid are performed for a wide range of Stokes number ranging in [10-1, 104]. Good agreement with available experimental data is found provided that a local lubrication model is used. Finally, three-dimensional simulations of gravity/shear-driven dense granular flows in a viscous fluid are presented. The results open the way for a parametric study in the parameter space initial aspect ratio - initial packing