Priyanka Shukla

Grad's method of moments is employed to develop higher-order Grad moment equations-up to the first 26 moments-for dilute granular gases within the framework of the (inelastic) Boltzmann equation. The homogeneous cooling state of a freely cooling granular gas is investigated with the Grad 26-moment equations in a semi-linearized setting and it is shown that the granular temperature in the homogeneous cooling state still decays according to Haff's law while the other higher-order moments decay on a faster time scale. The nonlinear terms of the fully contracted fourth moment are also considered and, by exploiting the stability analysis of fixed points, it is shown that these nonlinear terms have a negligible effect on Haff's law. Furthermore, an even larger Grad moment system, which includes the fully contracted sixth moment, is also scrutinized and the stability analysis of fixed points is again exploited to conclude that even the inclusion of the scalar sixth-order moment into the Grad moment system has a negligible effect on Haff's law. The constitutive relations for the stress and heat flux (i.e. the Navier-Stokes and Fourier relations) are derived through the Grad 26-moment equations and compared with those obtained via the Chapman-Enskog expansion and via computer simulations. The linear stability of the homogeneous cooling state is analysed through the Grad 26-moment system and various subsystems by decomposing them into longitudinal and transverse systems. It is found that one eigenmode in both longitudinal and transverse systems in the case of inelastic gases is unstable. By comparing the eigenmodes from various theories, it is established that the 13-moment eigenmode theory predicts that the unstable heat mode of the longitudinal system remains unstable for all wavenumbers below a certain coefficient of restitution, while any other higher-order moment theory shows that this mode becomes stable above some critical wavenumber for all values of the coefficient of restitution. In particular, the Grad 26-moment theory leads to a smooth profile for the critical wavenumber, in contrast to the other considered theories. Furthermore, the critical system size obtained through the Grad 26-moment theory is in excellent agreement with that obtained through existing theories.

One prototypical instability in granular flows is the shear-banding instability, in which a uniform granular shear flow breaks into alternating bands of dense and dilute clusters of particles having low and high shear (shear stress or shear rate), respectively. In this work, the shear-banding instability in an arbitrarily inelastic granular shear flow is analyzed through the linear stability analysis of granular hydrodynamic equations closed with Navier-Stokes-level constitutive relations. It is shown that the choice of appropriate constitutive relations plays an important role in predicting the shear-banding instability. A parametric study is carried out to study the effect of the restitution coefficient, channel width, and mean density. Two global criteria relating the control parameters are found for the onset of the shear-banding instability.

A fingering instability can develop at the interface between two fluids when the more mobile fluid is injected into the less-mobile one. For example, viscous fingering appears when a less viscous (i.e., more mobile) fluid displaces a more viscous (and hence less mobile) one in a porous medium. Fingering can also be due to a local change in mobility arising when a precipitation reaction locally decreases the permeability. We numerically analyze the properties of the related precipitation fingering patterns occurring when an A+B→C chemical reaction takes place, where A and B are reactants in solution and C is a solid product. We show that, similarly to reactive viscous fingering patterns, the precipitation fingering structures differ depending on whether A invades B or vice versa. This asymmetry can be related to underlying asymmetric concentration profiles developing when diffusion coefficients or initial concentrations of the reactants differ. In contrast to reactive viscous fingering, however, precipitation fingering patterns appear at shorter time scales than viscous fingers because the solid product C has a diffusivity tending to zero which destabilizes the displacement. Moreover, contrary to reactive viscous fingering, the system is more unstable with regard to precipitation fingering when the high-concentrated solution is injected into the low-concentrated one or when the faster diffusing reactant displaces the slower diffusing one.