7 results
Higher-order moment theories for dilute granular gases of smooth hard spheres
- Vinay Kumar Gupta, Priyanka Shukla, Manuel Torrilhon
-
- Journal:
- Journal of Fluid Mechanics / Volume 836 / 10 February 2018
- Published online by Cambridge University Press:
- 12 December 2017, pp. 451-501
-
- Article
- Export citation
-
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.
Nonlinear instability and convection in a vertically vibrated granular bed
- Priyanka Shukla, Istafaul H. Ansari, Devaraj van der Meer, Detlef Lohse, Meheboob Alam
-
- Journal:
- Journal of Fluid Mechanics / Volume 761 / 25 December 2014
- Published online by Cambridge University Press:
- 17 November 2014, pp. 123-167
-
- Article
- Export citation
-
The nonlinear instability of the density-inverted granular Leidenfrost state and the resulting convective motion in strongly shaken granular matter are analysed via a weakly nonlinear analysis of the hydrodynamic equations. The base state is assumed to be quasi-steady and the effect of harmonic shaking is incorporated by specifying a constant granular temperature at the vibrating plate. Under these mean-field assumptions, the base-state temperature decreases with increasing height away from the vibrating plate, but the density profile consists of three distinct regions: (i) a collisional dilute layer at the bottom, (ii) a levitated dense layer at some intermediate height and (iii) a ballistic dilute layer at the top of the granular bed. For the nonlinear stability analysis (Shukla & Alam, J. Fluid Mech., vol. 672, 2011b, pp. 147–195), the nonlinearities up to cubic order in the perturbation amplitude are retained, leading to the Landau equation, and the related adjoint stability problem is formulated taking into account appropriate boundary conditions. The first Landau coefficient and the related modal eigenfunctions (the fundamental mode and its adjoint, the second harmonic and the base-flow distortion, and the third harmonic and the cubic-order distortion to the fundamental mode) are calculated using a spectral-based numerical method. The genesis of granular convection is shown to be tied to a supercritical pitchfork bifurcation from the density-inverted Leidenfrost state. Near the bifurcation point the equilibrium amplitude ($A_{e}$) is found to follow a square-root scaling law, $A_{e}\sim \sqrt{{\it\Delta}}$, with the distance ${\it\Delta}$ from the bifurcation point. We show that the strength of convection (measured in terms of velocity circulation) is maximal at some intermediate value of the shaking strength, with weaker convection at both weaker and stronger shaking. Our theory predicts that at very strong shaking the convective motion remains concentrated only near the top surface, with the bulk of the expanded granular bed resembling the conduction state of a granular gas, dubbed as a floating-convection state. The linear and nonlinear patterns of the density and velocity fields are analysed and compared with experiments qualitatively. Evidence of 2:1 resonance is shown for certain parameter combinations. The influences of bulk viscosity, effective Prandtl number, shear work and free-surface boundary conditions on nonlinear equilibrium states are critically assessed.
Nonlinear vorticity-banding instability in granular plane Couette flow: higher-order Landau coefficients, bistability and the bifurcation scenario
- Priyanka Shukla, Meheboob Alam
-
- Journal:
- Journal of Fluid Mechanics / Volume 718 / 10 March 2013
- Published online by Cambridge University Press:
- 08 February 2013, pp. 131-180
-
- Article
- Export citation
-
The rapid granular plane Couette flow is known to be unstable to pure spanwise perturbations (i.e. perturbations having variations only along the mean vorticity direction) below some critical density (volume fraction of particles), resulting in the banding of particles along the mean vorticity direction: this is dubbed ‘vorticity banding’ instability. The nonlinear state of this instability is analysed using quintic-order Landau equation that has been derived from the pertinent hydrodynamic equations of rapid granular fluid. We have found analytical solutions for related modal/harmonic equations of finite-size perturbations up to quintic order in perturbation amplitude, leading to an exact calculation of both first and second Landau coefficients. This helped to identify the bistable nature of nonlinear vorticity-banding instability for a range of densities spanning from moderately dense to dense flows. For perturbations with small spanwise wavenumbers, the bifurcation scenario for vorticity banding unfolds, with increasing density from the dilute limit, as supercritical pitchfork $\rightarrow $ subcritical pitchfork $\rightarrow $ subcritical Hopf bifurcations. The transition from supercritical to subcritical pitchfork bifurcations is found to occur via the appearance of a degenerate/bicritical point (at which both the linear growth rate and the first Landau coefficient are simultaneously zero) that divides the critical line into two parts: one representing the first-order and the other the second-order phase transitions. Both subcritical oscillatory and stationary solutions have also been uncovered for dilute and dense flows, respectively, when the spanwise wavenumber is large. In all cases, the nonlinear solutions correspond to inhomogeneous states of shear stress and pressure along the vorticity direction, and hence are analogues of vorticity banding in other complex fluids. The quartic-order mean-flow resonance is evidenced in the parameter space for which the second Landau coefficient undergoes a jump discontinuity of infinite order. The importance of retaining higher-order terms to calculate the second Landau coefficient and their possible effects on the nature of bifurcations are elucidated.
Nonlinear stability, bifurcation and vortical patterns in three-dimensional granular plane Couette flow
- Meheboob Alam, Priyanka Shukla
-
- Journal:
- Journal of Fluid Mechanics / Volume 716 / 10 February 2013
- Published online by Cambridge University Press:
- 25 January 2013, pp. 349-413
-
- Article
- Export citation
-
The effects of three-dimensional (3D) perturbations, having wave-like modulations along both the streamwise and spanwise/vorticity directions, on the nonlinear states of five types of linear instability modes, the nature of their bifurcations and the resulting nonlinear patterns are analysed for granular plane Couette flow using an order-parameter theory which is an extension of our previous work on two-dimensional (2D) perturbations (Shukla & Alam, J. Fluid Mech., vol. 672, 2011b, pp. 147–195). The differential equations for modal amplitudes (the fundamental mode, the mean-flow distortion, the second harmonic and the distortion of the fundamental mode), up to cubic-order in perturbation amplitude, are solved using a spectral-based numerical technique, yielding an estimate of the first Landau coefficient that accounts for the leading-order nonlinear effect on finite-amplitude perturbations. In the near-critical regime of flows, we found evidence of mean-flow resonance, characterized by the divergence of the first Landau coefficient, that occurs due to the interaction/resonance between a linear instability mode and a mean-flow mode. The nonlinear solutions are found to appear via both pitchfork and Hopf bifurcations from the underlying linear instability modes, leading to supercritical nonlinear states of stationary and travelling wave solutions. The subcritical travelling wave solutions have also been uncovered in the linearly stable regimes of flow. It is shown that multiple nonlinear states of both stationary and travelling waves can coexist for a given parameter combination of mean density and Couette gap. The 3D nonlinear solutions persist for a range of spanwise wavenumbers up to ${k}_{z} = O(1)$ that originate from 2D instabilities which occur beyond a moderate value of the mean density. For purely 3D instabilities in dilute flows (having no analogue in 2D flows), the supercritical finite-amplitude solutions persist for a much larger range of spanwise wavenumber up to ${k}_{z} = O(10)$. For all instabilities, the vortical motion on the cross-stream plane has been characterized in terms of the fixed/critical points of the underlying flow field: saddles, nodes (sources and sinks) and vortices have been identified. While the cross-stream velocity field for supercritical solutions in dilute flows contains nodes and saddles, the subcritical solutions are dominated by large-scale vortices in the background of saddle-node-type motions. The latter type of flow pattern also persists at moderate densities in the form of supercritical nonlinear solutions that originate from the dominant 2D instability modes for which the vortex appears to be driven by two nearby saddles. The location of this vortex is found to be correlated with the local maxima of the streamwise vorticity.
Nonlinear stability and patterns in granular plane Couette flow: Hopf and pitchfork bifurcations, and evidence for resonance
- PRIYANKA SHUKLA, MEHEBOOB ALAM
-
- Journal:
- Journal of Fluid Mechanics / Volume 672 / 10 April 2011
- Published online by Cambridge University Press:
- 18 February 2011, pp. 147-195
-
- Article
- Export citation
-
The first evidence of a variety of nonlinear equilibrium states of travelling and stationary waves is provided in a two-dimensional granular plane Couette flow via nonlinear stability analysis. The relevant order-parameter equation, the Landau equation, has been derived for the most unstable two-dimensional perturbation of finite size. Along with the linear eigenvalue problem, the mean-flow distortion, the second harmonic, the distortion to the fundamental mode and the first Landau coefficient are calculated using a spectral-based numerical method. Two types of bifurcations, Hopf and pitchfork, that result from travelling and stationary instabilities, respectively, are analysed using the first Landau coefficient. The present bifurcation theory shows that the flow is subcritically unstable to stationary finite-amplitude perturbations of long wavelengths (kx ~ 0, where kx is the streamwise wavenumber) in the dilute limit that evolve from subcritical shear-banding modes (kx = 0), but at large enough Couette gaps there are stationary instabilities with kx = O(1) that lead to supercritical pitchfork bifurcations. At moderate-to-large densities, in addition to supercritical shear-banding modes, there are long-wave travelling instabilities that lead to Hopf bifurcations. It is shown that both supercritical and subcritical nonlinear states exist at moderate-to-large densities that originate from the dominant stationary and travelling instabilities for which kx = O(1). Nonlinear patterns of density, velocity and granular temperature for all types of instabilities are contrasted with their linear eigenfunctions. While the supercritical solutions appear to be modulated forms of the fundamental mode, the structural features of unstable subcritical solutions are found to be significantly different from their linear counterparts. It is shown that the granular plane Couette flow is prone to nonlinear resonances in both stable and unstable regimes, the signature of which is implicated as a discontinuity in the first Landau coefficient. Our analysis identified two types of modal resonances that appear at the quadratic order in perturbation amplitude: (i) a ‘mean-flow resonance’ which occurs due to the interaction between a streamwise-independent shear-banding mode (kx = 0) and a linear/fundamental mode kx ≠ 0, and (ii) an exact ‘1 : 2 resonance’ that results from the interaction between two waves with their wavenumber ratio being 1 : 2.
Weakly nonlinear theory of shear-banding instability in a granular plane Couette flow: analytical solution, comparison with numerics and bifurcation
- PRIYANKA SHUKLA, MEHEBOOB ALAM
-
- Journal:
- Journal of Fluid Mechanics / Volume 666 / 10 January 2011
- Published online by Cambridge University Press:
- 16 November 2010, pp. 204-253
-
- Article
- Export citation
-
A weakly nonlinear theory, in terms of the well-known Landau equation, has been developed to describe the nonlinear saturation of the shear-banding instability in a rapid granular plane Couette flow using the amplitude expansion method. The nonlinear modes are found to follow certain symmetries of the base flow and the fundamental mode, which helped to identify analytical solutions for the base-flow distortion and the second harmonic, leading to an exact calculation of the first Landau coefficient. The present analytical solutions are used to validate a spectral-based numerical method for the nonlinear stability calculation. The regimes of supercritical and subcritical bifurcations for the shear-banding instability have been identified, leading to the prediction that the lower branch of the neutral stability contour in the (H, φ0)-plane, where H is the scaled Couette gap (the ratio between the Couette gap and the particle diameter) and φ0 is the mean density or the volume fraction of particles, is subcritically unstable. The predicted finite-amplitude solutions represent shear localization and density segregation along the gradient direction. Our analysis suggests that there is a sequence of transitions among three types of pitchfork bifurcations with increasing mean density: from (i) the bifurcation from infinity in the Boltzmann limit to (ii) subcritical bifurcation at moderate densities to (iii) supercritical bifurcation at larger densities to (iv) subcritical bifurcation in the dense limit and finally again to (v) supercritical bifurcation near the close packing density. It has been shown that the appearance of subcritical bifurcation in the dense limit depends on the choice of the contact radial distribution function and the constitutive relations. The scalings of the first Landau coefficient, the equilibrium amplitude and the phase diagram, in terms of mode number and inelasticity, have been demonstrated. The granular plane Couette flow serves as a paradigm that supports all three possible types of pitchfork bifurcations, with the mean density (φ0) being the single control parameter that dictates the nature of the bifurcation. The predicted bifurcation scenario for the shear-band formation is in qualitative agreement with particle dynamics simulations and the experiment in the rapid shear regime of the granular plane Couette flow.
Universality of shear-banding instability and crystallization in sheared granular fluid
- MEHEBOOB ALAM, PRIYANKA SHUKLA, STEFAN LUDING
-
- Journal:
- Journal of Fluid Mechanics / Volume 615 / 25 November 2008
- Published online by Cambridge University Press:
- 25 November 2008, pp. 293-321
-
- Article
- Export citation
-
The linear stability analysis of an uniform shear flow of granular materials is revisited using several cases of a Navier–Stokes-level constitutive model in which we incorporate the global equation of states for pressure and thermal conductivity (which are accurate up to the maximum packing density νm) and the shear viscosity is allowed to diverge at a density νμ (<νm), with all other transport coefficients diverging at νm. It is shown that the emergence of shear-banding instabilities (for perturbations having no variation along the streamwise direction), that lead to shear-band formation along the gradient direction, depends crucially on the choice of the constitutive model. In the framework of a dense constitutive model that incorporates only collisional transport mechanism, it is shown that an accurate global equation of state for pressure or a viscosity divergence at a lower density or a stronger viscosity divergence (with other transport coefficients being given by respective Enskog values that diverge at νm) can induce shear-banding instabilities, even though the original dense Enskog model is stable to such shear-banding instabilities. For any constitutive model, the onset of this shear-banding instability is tied to a universal criterion in terms of constitutive relations for viscosity and pressure, and the sheared granular flow evolves toward a state of lower ‘dynamic’ friction, leading to the shear-induced band formation, as it cannot sustain increasing dynamic friction with increasing density to stay in the homogeneous state. A similar criterion of a lower viscosity or a lower viscous-dissipation is responsible for the shear-banding state in many complex fluids.