Search

The hydrodynamics and rheology of a sheared dilute gas-solid suspension, consisting of inelastic hard-spheres suspended in a gas, are analysed using anisotropic Maxwellian as the single particle distribution function. The closed-form solutions for granular temperature and three invariants of the second-moment tensor are obtained as functions of the Stokes number ($St$), the mean density ($\nu$) and the restitution coefficient ($e$). Multiple states of high and low temperatures are found when the Stokes number is small, thus recovering the "ignited" and "quenched" states, respectively, of Tsao \& Koch (J. Fluid Mech.,1995). The phase diagram is constructed in the three-dimensional ($\nu, St, e$)-space that delineates the regions of ignited and quenched states and their coexistence. Analytical expressions for the particle-phase shear viscosity and the normal stress differences are obtained, along with related scaling relations on the quenched and ignited states. At any $e$, the shear-viscosity undergoes a discontinuous jump with increasing shear rate (i.e.~ discontinuous shear-thickening) at the "quenched-ignited" transition. The first (${\mathcal N}_1$) and second (${\mathcal N}_2$) normal-stress differences also undergo similar first-order transitions: (i) ${\mathcal N}_1$ jumps from large to small positive values and (ii) ${\mathcal N}_2$ from positive to negative values with increasing $St$, with the sign-change of ${\mathcal N}_2$ identified with the system making a transition from the quenched to ignited states. The superior prediction of the present theory over the standard Grad's method and the Chapman-Enskog solution is demonstrated via comparisons of transport coefficients with simulation data for a range of Stokes number and restitution coefficient.

A new method is introduced to simulate jamming of polyhedral grains under controlled stress that incorporates global degrees of freedom through the metric tensor of a periodic cell containing grains. Jamming under hydrostatic/isotropic stress and athermal conditions leads to a precise definition of the ideal jamming point at zero shear stress. The structures of tetrahedra jammed hydrostatically exhibit less translational order and lower jamming-point density than previously described `maximally random jammed' hard tetrahedra. Under the same conditions, cubes jam with negligible nematic order. Grains with octahedral symmetry jam in the large-system limit with an abundance of face-face contacts in the absence of nematic order. For sufficiently large face-face contact number, percolating clusters form that span the entire simulation box. The response of hydrostatically jammed tetrahedra and cubes to shear-stress perturbation is also demonstrated with the variable-cell method., Comment: 10 pages, 8 figures

The average number of constraints per particle $< C_{total} >$ in mechanically stable systems of Platonic solids (except cubes) approaches the isostatic limit at the jamming point ($< C_{total} > \rightarrow 12$), though average number of contacts are hypostatic. By introducing angular alignment metrics to classify the degree of constraint imposed by each contact, constraints are shown to arise as a direct result of local orientational order reflected in edge-face and face-face alignment angle distributions. With approximately one face-face contact per particle at jamming chain-like face-face clusters with finite extent form in these systems., Comment: 4 pages, 3 figures, 4 table

A mechanically-based structural optimization method is utilized to explore the phenomena of jamming for assemblies of frictionless Platonic solids. Systems of these regular convex polyhedra exhibit mechanically stable phases with density substantially less than optimal for a given shape, revealing that thermal motion is necessary to access high density phases. We confirm that the large system jamming threshold of $0.623 \pm 0.003$ for tetrahedra is consistent with experiments on tetrahedral dice. Also, the extremely short-ranged translational correlations of packed tetrahedra observed in experiments are confirmed here, in contrast with those of thermally simulated glasses. Though highly ordered phases are observed to form for small numbers of cubes and dodecahedra, the short correlation length scale suppresses ordering in large systems, resulting in packings that are mechanically consistent with `orientationally disordered' contacts (point-face and edge-edge contacts). Mild nematic ordering is observed for large systems of cubes, whereas angular correlations for the remaining shapes are ultra short-ranged. Power-law scaling exponents for energy with respect to distance from the jamming threshold exhibit a clear dependence on the `highest order' percolating contact topology. These nominal exponents are 6, 4, and 2 for configurations having percolating point-face (or edge-edge), edge-face, and face-face contacts, respectively. These Platonic solids exhibit hypostatic behavior, with average jamming contact number between the isostatic value for spheres and that of asymmetric particles. These shapes violate the isostatic conjecture, displaying contact number that decreases monotonically with sphericity. The common symmetry of dual polyhedra results in local translational structural similarity. (Note abstract abridged due to length constraint), Comment: 49 pages, 13 figures, 2 tables

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 $\nu_{m}$) and the shear viscosity is allowed to diverge at a density $\nu_\mu$ ($< \nu_{m}$), with all other transport coefficients diverging at $\nu_{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 $\nu_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 {\it 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., Comment: 26 pages

Linear stability and the non-modal transient energy growth in compressible plane Couette flow are investigated for two prototype mean flows: (a) the {\it uniform shear} flow with constant viscosity, and (b) the {\it non-uniform shear} flow with {\it stratified} viscosity. Both mean flows are linearly unstable for a range of supersonic Mach numbers ($M$). For a given $M$, the critical Reynolds number ($Re$) is significantly smaller for the uniform shear flow than its non-uniform shear counterpart. An analysis of perturbation energy reveals that the instability is primarily caused by an excess transfer of energy from mean-flow to perturbations. It is shown that the energy-transfer from mean-flow occurs close to the moving top-wall for ``mode I'' instability, whereas it occurs in the bulk of the flow domain for ``mode II''. For the non-modal analysis, it is shown that the maximum amplification of perturbation energy, $G_{\max}$, is significantly larger for the uniform shear case compared to its non-uniform counterpart. For $\alpha=0$, the linear stability operator can be partitioned into ${\cal L}\sim \bar{\cal L} + Re^2{\cal L}_p$, and the $Re$-dependent operator ${\cal L}_p$ is shown to have a negligibly small contribution to perturbation energy which is responsible for the validity of the well-known quadratic-scaling law in uniform shear flow: $G(t/{\it Re}) \sim {\it Re}^2$. A reduced inviscid model has been shown to capture all salient features of transient energy growth of full viscous problem. For both modal and non-modal instability, it is shown that the {\it viscosity-stratification} of the underlying mean flow would lead to a delayed transition in compressible Couette flow.

Nonmodal transient growth studies and estimation of optimal perturbations have been made for the compressible plane Couette flow with three-dimensional disturbances. The maximum amplification of perturbation energy over time, $G_{\max}$, is found to increase with increasing Reynolds number ${\it Re}$, but decreases with increasing Mach number $M$. More specifically, the optimal energy amplification $G_{\rm opt}$ (the supremum of $G_{\max}$ over both the streamwise and spanwise wavenumbers) is maximum in the incompressible limit and decreases monotonically as $M$ increases. The corresponding optimal streamwise wavenumber, $\alpha_{\rm opt}$, is non-zero at M=0, increases with increasing $M$, reaching a maximum for some value of $M$ and then decreases, eventually becoming zero at high Mach numbers. While the pure streamwise vortices are the optimal patterns at high Mach numbers, the modulated streamwise vortices are the optimal patterns for low-to-moderate values of the Mach number. Unlike in incompressible shear flows, the streamwise-independent modes in the present flow do not follow the scaling law $G(t/{\it Re}) \sim {\it Re}^2$, the reasons for which are shown to be tied to the dominance of some terms in the linear stability operator. Based on a detailed nonmodal energy analysis, we show that the transient energy growth occurs due to the transfer of energy from the mean flow to perturbations via an inviscid {\it algebraic} instability. The decrease of transient growth with increasing Mach number is also shown to be tied to the decrease in the energy transferred from the mean flow ($\dot{\mathcal E}_1$) in the same limit.

Using particle simulations of the uniform shear flow of a rough dilute granular gas, we show that the translational and rotational velocities are strongly correlated in direction, but there is no orientational correlation-induced singularity at perfectly smooth ($\beta=-1$) and rough ($\beta=1$) limits for elastic collisions ($e=1$); both the translational and rotational velocity distribution functions remain close to a Gaussian for these two limiting cases. Away from these two limits, the orientational as well as spatial velocity correlations are responsible for the emergence of non-Gaussian high velocity tails. The tails of both distribution functions follow stretched exponentials, with the exponents depending on normal ($e$) and tangential ($\beta$) restitution coefficients., Comment: Physical Review Letters (accepted)

The effect of Prandtl number on the linear stability of a plane thermal plume is analyzed under quasi-parallel approximation. At large Prandtl numbers ($Pr>100$), we found that there is an additional unstable loop whose size increases with increasing $Pr$. The origin of this new instability mode is shown to be tied to the coupling of the momentum and thermal perturbation equations. Analyses of the perturbation kinetic energy and thermal energy suggest that the buoyancy force is the main source of perturbation energy at high Prandtl numbers that drives this instability.

Based on the Boltzmann-Enskog kinetic theory, we develop a hydrodynamic theory for the well known (reverse) Brazil nut segregation in a vibro-fluidized granular mixture. Using an analogy with standard fluid mechanics, we have recently suggested a novel mechanism of segregation in granular mixtures based on a {\it competition between buoyancy and geometric forces}: the Archimedean buoyancy force, a pseudo-thermal buoyancy force due to the difference between the energies of two granular species, and two geometric forces, one compressive and the other-one tensile in nature, due to the size-difference. For a mixture of perfectly hard-particles with elastic collisions, the pseudo-thermal buoyancy force is zero but the intruder has to overcome the net compressive geometric force to rise. For this case, the geometric force competes with the standard Archimedean buoyancy force to yield a threshold density-ratio, $R_{\rho 1}=\rho_l/\rho_s < 1$, above which the {\it lighter intruder sinks}, thereby signalling the {\it onset} of the {\it reverse buoyancy} effect. For a mixture of dissipative particles, the non-zero pseudo-thermal buoyancy force gives rise to another threshold density-ratio, $R_{\rho 2}$ ($> R_{\rho 1}$), above which the intruder rises again. Focussing on the {\it tracer} limit of intruders in a dense binary mixture, we find that the rise-time of the intruder could vary {\it non-monotonically} with the density-ratio. For a given size-ratio, there is a threshold density-ratio for the intruder at which it takes the maximum time to rise, and above(/below) which it rises faster, implying that {\it the heavier (and larger) the intruder, the faster it ascends}. Our theory offers a unified description for the (reverse) Brazil-nut segregation and the non-monotonic ascension dynamics of Brazil-nuts.

From event-driven simulations of a gravity-driven channel flow of inelastic hard-disks, we show that the velocity distribution function remains close to a Gaussian for a wide range densities (even when the Knudsen number is of order one) if the walls are smooth and the particle collisions are nearly elastic. For dense flows, a transition from a Gaussian to a power-law distribution for the high velocity tails occurs with increasing dissipation in the center of the channel, irrespective of wall-roughness. For a rough wall, the near-wall distribution functions are distinctly different from those in the bulk even in the quasielastic limit.

A numerical bifurcation analysis is carried out to understand the role of gravity on the shear-band formation in rapid granular Couette flow. At {\it low} shear rates, there is a unique solution with a {\it plug} near the bottom wall and a {\it shear-layer} near the top-wall; this solution mirrors typical shearbanding-type profiles in earth-bound shear-cell experiments. Interestingly, a {\it stable} plug near the top-wall is also a solution of these equations at {\it high} shear rates; there is a multitude of other plugged states, with the plugs being located in an ordered fashion within the Couette gap. The origin of such shearbanding solutions is tied to the spontaneous symmetry-breaking {\it shearbanding} instabilities of the gravity-free uniform shear flow, leading to both subcritical and supercritical pitchfork bifurcations. In the language of singularity theory, we have established that this bifurcation problem admits {\it universal unfolding} of pitchfork bifurcations., Comment: 11 pages and 6 figures; Proceedings of the Symposium on Trends in Applications of Mathematics to Mechanics, 22-28 August, 2004, Seeheim-Jugenheim, Germnay (Editors: K. Hutter and Y. Wang)

The first normal stress difference (${\mathcal N}_1$) and the microstructure in a dense sheared granular fluid of smooth inelastic hard-disks are probed using event-driven simulations. While the anisotropy in the second moment of fluctuation velocity, which is a Burnett-order effect, is known to be the progenitor of normal stress differences in {\it dilute} granular fluids, we show here that the collisional anisotropies are responsible for the normal stress behaviour in the {\it dense} limit. As in the elastic hard-sphere fluids, ${\mathcal N}_1$ remains {\it positive} (if the stress is defined in the {\it compressive} sense) for dilute and moderately dense flows, but becomes {\it negative} above a critical density, depending on the restitution coefficient. This sign-reversal of ${\mathcal N}_1$ occurs due to the {\it microstructural} reorganization of the particles, which can be correlated with a preferred value of the {\it average} collision angle $\theta_{av}=\pi/4 \pm \pi/2$ in the direction opposing the shear. We also report on the shear-induced {\it crystal}-formation, signalling the onset of fluid-solid coexistence in dense granular fluids. Different approaches to take into account the normal stress differences are discussed in the framework of the relaxation-type rheological models., Comment: 21 pages, 13 figures