Kinetic theory and hydrodynamics of dense, reacting fluids far from equilibrium

Energy nonequipartition in a collisional model of a confined quasi-two-dimensional granular mixture

A collisional model of a confined quasi-two-dimensional granular mixture is considered to analyze homogeneous steady states. The model includes an effective mechanism to transfer the kinetic energy injected by vibration in the vertical direction to the horizontal degrees of freedom of grains. The set of Enskog kinetic equations for the velocity distribution functions of each component is derived first to analyze the homogeneous state. As in the one-component case, an exact scaling solution is found where the time dependence of the distribution functions occurs entirely through the granular temperature T. As expected, the kinetic partial temperatures T_{i} of each component are different and, hence, energy equipartition is broken down. In the steady state, explicit expressions for the temperature T and the ratio of partial kinetic temperatures T_{i}/T_{j} are obtained by considering Maxwellian distributions defined at the partial temperatures T_{i}. The (scaled) granular temperature and the temperature ratios are given in terms of the coefficients of restitution, the solid volume fraction, the (scaled) parameters of the collisional model, and the ratios of mass, concentration, and diameters. In the case of a binary mixture, the theoretical predictions are exhaustively compared with both direct simulation Monte Carlo and molecular dynamics simulations with a good agreement. The deviations are identified to be originated in the non-Gaussianity of the velocity distributions and on microsegregation patterns, which induce spatial correlations not captured in the Enskog theory.

Homogeneous steady state of a confined granular gas

The nonequilibrium statistical mechanics and kinetic theory for a model of a confined quasi-two-dimensional gas of inelastic hard spheres is presented. The dynamics of the particles includes an effective mechanism to transfer the energy injected in the vertical direction to the horizontal degrees of freedom. The Enskog approximation is formulated and used as the basis to investigate the temperature and the distribution function of the steady state eventually reached by the system. An exact scaling of the distribution function of the system having implications on the form of its moments is pointed out. The theoretical predictions are compared with numerical results obtained by a particle simulation method, and a good agreement is found.

Hydrodynamic modes in a confined granular fluid

Confined granular fluids, placed in a shallow box that is vibrated vertically, can achieve homogeneous stationary states due to energy injection mechanisms that take place throughout the system. These states can be stable even at high densities and inelasticities allowing for a detailed analysis of the hydrodynamic modes that govern the dynamics of granular fluids. By analyzing the decay of the time correlation functions it is shown that there is a crossover from a quasielastic regime in which energy evolves as a slow mode to an inelastic regime with energy slaved to the other conserved fields. The two regimes have well differentiated transport properties and in the inelastic regime the dynamics can be described by a reduced hydrodynamics with modified longitudinal viscosity and sound speed. The crossover between the two regimes takes place at a wave vector that is proportional to the inelasticity. A two-dimensional granular model, with collisions that mimic the energy transfers that take place in a confined system, is studied by means of microscopic simulations. The results show excellent agreement with the theoretical framework and allow validation of hydrodynamiclike models.

On the form of the kinetic energy balance equation in the kinetic variational theory

An alternative balance equation for the kinetic energy density used in the kinetic variational theory (KVT) is proposed. The new equation is consistent with the well-known standard form interpreted in terms of a flux and a source for both a one-component fluid and a mixture. Within the proposed version, the KVT mean-field collision term produces contributions to the heat flux and the source of kinetic energy being absent in the original formulation. It is shown that the introduced modification can affect the KVT thermal conductivity of the mixture while for the one-component fluid it becomes important only in the second and higher orders in gradients.

Transport coefficients for the hard-sphere granular fluid

In the preceding paper, linear response methods have been applied to obtain formally exact expressions for the parameters of Navier-Stokes order hydrodynamics. The analysis there is general, applying to both normal and granular fluids with a wide range of collision rules. Those results are specialized here to the case of smooth, inelastic, hard spheres with constant coefficient of normal restitution, for further elaboration. Explicit expressions for the cooling rate, pressure, and transport coefficients are given and compared with the corresponding expressions for a system of elastic hard spheres. The scope of the results for further analytical explorations and possible numerical evaluation is discussed.

Phase-space approach to dynamical density functional theory

The authors consider a system of interacting particles subjected to Langevin inertial dynamics and derive the governing time-dependent equation for the one-body density. They show that, after suitable truncations of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy, and a multiple time scale analysis, they obtain a self-consistent equation involving only the one-body density. This study extends to arbitrary dimensions previous work on a one-dimensional fluid and highlights the subtleties of kinetic theory in the derivation of dynamical density functional theory.

Hydrodynamics of an endothermic gas with application to bubble cavitation

The hydrodynamics for a gas of hard spheres which sometimes experience inelastic collisions resulting in the loss of a fixed, velocity-independent, amount of energy Delta is investigated with the goal of understanding the coupling between hydrodynamics and endothermic chemistry. The homogeneous cooling state of a uniform system and the modified Navier-Stokes equations are discussed and explicit expressions given for the pressure, cooling rates, and all transport coefficients for D dimensions. The Navier-Stokes equations are solved numerically for the case of a two-dimensional gas subject to a circular piston so as to illustrate the effects of the energy loss on the structure of shocks found in cavitating bubbles. It is found that the maximal temperature achieved is a sensitive function of Delta with a minimum occurring near the physically important value of Delta approximately 12,000 K approximately 1 eV.

Nonequilibrium inertial dynamics of colloidal systems

We consider the properties of a one-dimensional fluid of Brownian inertial hard-core particles, whose microscopic dynamics is partially damped by a heat bath. Direct interactions among the particles are represented as binary, instantaneous elastic collisions. Collisions with the heat bath are accounted for by a Fokker-Planck collision operator, whereas direct collisions among the particles are treated by a well known method of kinetic theory, the revised Enskog theory. By means of a time multiple time-scale method we derive the evolution equation for the average density. Remarkably, for large values of the friction parameter and/or of the mass of the particles we obtain the same equation as the one derived within the dynamic density functional theory (DDF). In addition, at moderate values of the friction constant, the present method allows to study the inertial effects not accounted for by DDF method. Finally, a numerical test of these corrections is provided.

Chapman-Enskog expansion about nonequilibrium states with application to the sheared granular fluid

The Chapman-Enskog method of solution of kinetic equations, such as the Boltzmann equation, is based on an expansion in gradients of the deviations of the hydrodynamic fields from a uniform reference state (e.g., local equilibrium). This paper presents an extension of the method so as to allow for expansions about arbitrary, far-from-equilibrium reference states. The primary result is a set of hydrodynamic equations for studying variations from the arbitrary reference state which, unlike the usual Navier-Stokes hydrodynamics, does not restrict the reference state in any way. The method is illustrated by application to a sheared granular gas which cannot be studied using the usual Navier-Stokes hydrodynamics.

Transport properties of dense dissipative hard-sphere fluids for arbitrary energy loss models

The revised Enskog approximation for a fluid of hard spheres which lose energy upon collision is discussed for the case that the energy is lost from the normal component of the velocity at collision but is otherwise arbitrary. Granular fluids with a velocity-dependent coefficient of restitution are an important special case covered by this model. A normal solution to the Enskog equation is developed using the Chapman-Enskog expansion. The lowest order solution describes the general homogeneous cooling state and a generating function formalism is introduced for the determination of the distribution function. The first order solution, evaluated in the lowest Sonine approximation, provides estimates for the transport coefficients for the Navier-Stokes hydrodynamic description. All calculations are performed in an arbitrary number of dimensions.

Hard sphere dynamics for normal and granular fluids

A fluid of N smooth, hard spheres is considered as a model for normal (elastic collision) and granular (inelastic collision) fluids. The potential energy is discontinuous for hard spheres so that the pairwise forces are singular and the usual forms of Newtonian and Hamiltonian mechanics do not apply. Nevertheless, particle trajectories in the N particle phase space are well defined and the generators for these trajectories can be identified. The first part of this presentation is a review of the generators for the dynamics of observables and probability densities. The new results presented in the second part refer to applications of these generators to the Liouville dynamics for granular fluids. A set of eigenvalues and eigenfunctions of the generator for this Liouville dynamics system is identified in a special stationary representation. This provides a class of exact solutions to the Liouville equation that are closely related to hydrodynamics for granular fluids.

Rheology of dense polydisperse granular fluids under shear

The solution of the Enskog equation for the one-body velocity distribution of a moderately dense arbitrary mixture of inelastic hard spheres undergoing planar shear flow is described. A generalization of the Grad moment method, implemented by means of a novel generating function technique, is used so as to avoid any assumptions concerning the size of the shear rate. The result is illustrated by using it to calculate the pressure, normal stresses, and shear viscosity of a model polydisperse granular fluid in which grain size, mass, and coefficient of restitution vary among the grains. The results are compared to a numerical solution of the Enskog equation as well as molecular-dynamics simulations. Most bulk properties are well described by the Enskog theory and it is shown that the generalized moment method is more accurate than the simple (Grad) moment method. However, the description of the distribution of temperatures in the mixture predicted by Enskog theory does not compare well to simulation, even at relatively modest densities.

Imploding shock wave in a fluid of hard-core particles

We report the study of a fluid of hard-disk particles in a contracting cavity. Under supersonic contraction speed, a shock wave converges to the center of the cavity where it implodes, creating a central peak in temperature. The dynamics of the fluid is studied by solving the Euler and Navier-Stokes equations, as well as by molecular dynamics simulations and the Enskog direct simulation Monte Carlo method. The value of the maximum temperature reached at the center of the cavity is systematically investigated with the different methods which give consistent results. Moreover, we develop a scaling theory for the maximum temperature based on the self-similar solutions of Euler's equations and mean-free-path considerations. This scaling theory provides a comprehensive scheme for the interpretation of the numerical results. In addition, the effects of the imploding shock wave on an passively driven isomerization reaction A <= => B are also studied.