Uniform shear flow in dissipative gases. Computer simulations of inelastic hard spheres and (frictional) elastic hard spheres

In the preceding paper (cond-mat/0405252), we have conjectured that the main transport properties of a dilute gas of inelastic hard spheres (IHS) can be satisfactorily captured by an equivalent gas of elastic hard spheres (EHS), provided that the latter are under the action of an effective drag force and their collision rate is reduced by a factor $(1+\alpha)/2$ (where $\alpha$ is the constant coefficient of normal restitution). In this paper we test the above expectation in a paradigmatic nonequilibrium state, namely the simple or uniform shear flow, by performing Monte Carlo computer simulations of the Boltzmann equation for both classes of dissipative gases with a dissipation range $0.5\leq \alpha\leq 0.95$ and two values of the imposed shear rate $a$. The distortion of the steady-state velocity distribution from the local equilibrium state is measured by the shear stress, the normal stress differences, the cooling rate, the fourth and sixth cumulants, and the shape of the distribution itself. In particular, the simulation results seem to be consistent with an exponential overpopulation of the high-velocity tail. The EHS results are in general hardly distinguishable from the IHS ones if $\alpha\gtrsim 0.7$, so that the distinct signature of the IHS gas (higher anisotropy and overpopulation) only manifests itself at relatively high dissipationsComment: 23 pages; 18 figures; Figs. 2 and 9 include new simulations; two new figures added; few minor changes; accepted for publication in PR

Prediction of strong shock structure using the bimodal distribution function

A modified Mott-Smith method for predicting the one-dimensional shock wave solution at very high Mach numbers is constructed by developing a system of fluid dynamic equations. The predicted shock solutions in a gas of Maxwell molecules, a hard sphere gas and in argon using the newly proposed formalism are compared with the experimental data, direct-simulation Monte Carlo (DSMC) solution and other solutions computed from some existing theories for Mach numbers M<50. In the limit of an infinitely large Mach number, the predicted shock profiles are also compared with the DSMC solution. The density, temperature and heat flux profiles calculated at different Mach numbers have been shown to have good agreement with the experimental and DSMC solutionsComment: 22 pages, 9 figures, Accepted for publication in Physical Review

On the role of the Knudsen layer in rapid granular flows

A combination of molecular-dynamics simulations, theoretical predictions, and previous experiments are used in a two-part study to determine the role of the Knudsen layer in rapid granular flows. First, a robust criterion for the identification of the thickness of the Knudsen layer is established: a rapid deterioration in Navier-Stokes-order prediction of the heat flux is found to occur in the Knudsen layer. For (experimental) systems in which heat flux measurements are not easily obtained, a rule-of-thumb for estimating the Knudsen layer thickness follows, namely that such effects are evident within 2.5 (local) mean free paths of a given boundary. Second, comparisons of simulation and experimental data with Navier-Stokes order theory are used to provide a measure as to when Knudsen layer effects become non-negligible. Specifically, predictions that do not account for the presence of a Knudsen layer appear reliable for Knudsen layers collectively composing up to 20% of the domain, whereas deterioration of such predictions becomes apparent when the domain is fully comprised of the Knudsen layer.Comment: 9 figures, accepted to Journal of Fluid Mechanic

Computer simulations of an impurity in a granular gas under planar Couette flow

We present in this work results from numerical solutions, obtained by means of the direct simulation Monte Carlo (DSMC) method, of the Boltzmann and Boltzmann--Lorentz equations for an impurity immersed in a granular gas under planar Couette flow. The DSMC results are compared with the exact solution of a recent kinetic model for the same problem. The results confirm that, in steady states and over a wide range of parameter values, the state of the impurity is enslaved to that of the host gas: it follows the same flow velocity profile, its concentration (relative to that of the granular gas) is constant in the bulk region, and the impurity/gas temperature ratio is also constant. We determine also the rheological properties and nonlinear hydrodynamic transport coefficients for the impurity, finding a good semi-quantitative agreement between the DSMC results and the theoretical predictions.Comment: 23 pages, 11 figures; v2: minor change

Class of dilute granular Couette flows with uniform heat flux

In a recent paper [F. Vega Reyes et al., Phys. Rev. Lett. 104, 028001 (2010)] we presented a preliminary description of a special class of steady Couette flows in dilute granular gases. In all flows of this class the viscous heating is exactly balanced by inelastic cooling. This yields a uniform heat flux and a linear relationship between the local temperature and flow velocity. The class (referred to as the LTu class) includes the Fourier flow of ordinary gases and the simple shear flow of granular gases as special cases. In the present paper we provide further support for this class of Couette flows by following four different routes, two of them being theoretical (Grad's moment method of the Boltzmann equation and exact solution of a kinetic model) and the other two being computational (molecular dynamics and Monte Carlo simulations of the Boltzmann equation). Comparison between theory and simulations shows a very good agreement for the non-Newtonian rheological properties, even for quite strong inelasticity, and a good agreement for the heat flux coefficients in the case of Grad's method, the agreement being only qualitative in the case of the kinetic model.Comment: 15 pages, 10 figures; v2: change of title plus some other minor change

Segregation of an intruder in a heated granular dense gas

A recent segregation criterion [V. Garz\'o, Phys. Rev. E \textbf{78}, 020301(R) (2008)] based on the thermal diffusion factor $\Lambda$ of an intruder in a heated granular gas described by the inelastic Enskog equation is revisited. The sign of $\Lambda$ provides a criterion for the transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE). The present theory incorporates two extra ingredients not accounted for by the previous theoretical attempt. First, the theory is based upon the second Sonine approximation to the transport coefficients of the mass flux of intruder. Second, the dependence of the temperature ratio (intruder temperature over that of the host granular gas) on the solid volume fraction is taken into account in the first and second Sonine approximations. In order to check the accuracy of the Sonine approximation considered, the Enskog equation is also numerically solved by means of the direct simulation Monte Carlo (DSMC) method to get the kinetic diffusion coefficient $D_0$. The comparison between theory and simulation shows that the second Sonine approximation to $D_0$ yields an improvement over the first Sonine approximation when the intruder is lighter than the gas particles in the range of large inelasticity. With respect to the form of the phase diagrams for the BNE/RBNE transition, the kinetic theory results for the factor $\Lambda$ indicate that while the form of these diagrams depends sensitively on the order of the Sonine approximation considered when gravity is absent, no significant differences between both Sonine solutions appear in the opposite limit (gravity dominates the thermal gradient). In the former case (no gravity), the first Sonine approximation overestimates both the RBNE region and the influence of dissipation on thermal diffusion segregation.Comment: 9 figures; to be published in Phys. Rev.

Thermal diffusion segregation in granular binary mixtures described by the Enskog equation

Diffusion induced by a thermal gradient in a granular binary mixture is analyzed in the context of the (inelastic) Enskog equation. Although the Enskog equation neglects velocity correlations among particles which are about to collide, it retains spatial correlations arising from volume exclusion effects and thus it is expected to apply to moderate densities. In the steady state with gradients only along a given direction, a segregation criterion is obtained from the thermal diffusion factor $\Lambda$ measuring the amount of segregation parallel to the thermal gradient. As expected, the sign of the factor $\Lambda$ provides a criterion for the transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE) by varying the parameters of the mixture (masses, sizes, concentration, solid volume fraction, and coefficients of restitution). The form of the phase diagrams for the BNE/RBNE transition is illustrated in detail for several systems, with special emphasis on the significant role played by the inelasticity of collisions. In particular, an effect already found in dilute gases (segregation in a binary mixture of identical masses and sizes {\em but} different coefficients of restitution) is extended to dense systems. A comparison with recent computer simulation results shows a good qualitative agreement at the level of the thermal diffusion factor. The present analysis generalizes to arbitrary concentration previous theoretical results derived in the tracer limit case.Comment: 7 figures, 1 table. To appear in New J. Phys., special issue on "Granular Segregation

Shocks in supersonic sand

We measure time-averaged velocity, density, and temperature fields for steady granular flow past a wedge and calculate a speed of granular pressure disturbances (sound speed) equal to 10% of the flow speed. The flow is supersonic, forming shocks nearly identical to those in a supersonic gas. Molecular dynamics simulations of Newton's laws and Monte Carlo simulations of the Boltzmann equation yield fields in quantitative agreement with experiment. A numerical solution of Navier-Stokes-like equations agrees with a molecular dynamics simulation for experimental conditions excluding wall friction.Comment: 4 pages, 5 figure

Relevance of initial and final conditions for the Fluctuation Relation in Markov processes

Numerical observations on a Markov chain and on the continuous Markov process performed by a granular tracer show that the ``usual'' fluctuation relation for a given observable is not verified for finite (but arbitrarily large) times. This suggests that some terms which are usually expected to be negligible, i.e. ``border terms'' dependent only on initial and final states, in fact cannot be neglected. Furthermore, the Markov chain and the granular tracer behave in a quite similar fashion.Comment: 23 pages, 5 figures, submitted to JSTA

Diffusion in a multi-component Lattice Boltzmann Equation model

Diffusion phenomena in a multiple component lattice Boltzmann Equation (LBE) model are discussed in detail. The mass fluxes associated with different mechanical driving forces are obtained using a Chapman-Enskog analysis. This model is found to have correct diffusion behavior and the multiple diffusion coefficients are obtained analytically. The analytical results are further confirmed by numerical simulations in a few solvable limiting cases. The LBE model is established as a useful computational tool for the simulation of mass transfer in fluid systems with external forces.Comment: To appear in Aug 1 issue of PR