Theory of thermostatted inhomogeneous granular fluids: A self-consistent density functional description

Dynamical density functional theory with hydrodynamic interactions in confined geometries

We study the dynamics of colloidal fluids in both unconfined geometries and when confined by a hard wall. Under minimal assumptions, we derive a dynamical density functional theory (DDFT) which includes hydrodynamic interactions (HI; bath-mediated forces). By using an efficient numerical scheme based on pseudospectral methods for integro-differential equations, we demonstrate its excellent agreement with the full underlying Langevin equations for systems of hard disks in partial confinement. We further use the derived DDFT formalism to elucidate the crucial effects of HI in confined systems.

Heat flux of driven granular mixtures at low density: Stability analysis of the homogeneous steady state

The Navier-Stokes order hydrodynamic equations for a low-density driven granular mixture obtained previously [Khalil and Garzó, Phys. Rev. E 88, 052201 (2013)PLEEE81539-375510.1103/PhysRevE.88.052201] from the Chapman-Enskog solution to the Boltzmann equation are considered further. The four transport coefficients associated with the heat flux are obtained in terms of the mass ratio, the size ratio, composition, coefficients of restitution, and the driven parameters of the model. Their quantitative variation on the control parameters of the system is demonstrated by considering the leading terms in a Sonine polynomial expansion to solve the exact integral equations. As an application of these results, the stability of the homogeneous steady state is studied. In contrast to the results obtained in undriven granular mixtures, the stability analysis of the linearized Navier-Stokes hydrodynamic equations shows that the transversal and longitudinal modes are (linearly) stable with respect to long enough wavelength excitations. This conclusion agrees with a previous analysis made for single granular gases.

Towards a statistical mechanical theory of active fluids

We present a stochastic description of a model of N mutually interacting active particles in the presence of external fields and characterize its steady state behavior in the absence of currents. To reproduce the effects of the experimentally observed persistence of the trajectories of the active particles we consider a Gaussian force having a non-vanishing correlation time τ, whose finiteness is a measure of the activity of the system. With these ingredients we show that it is possible to develop a statistical mechanical approach similar to the one employed in the study of equilibrium liquids and to obtain the explicit form of the many-particle distribution function by means of the multidimensional unified colored noise approximation. Such a distribution plays a role analogous to the Gibbs distribution in equilibrium statistical mechanics and provides complete information about the microscopic state of the system. From here we develop a method to determine the one- and two-particle distribution functions in the spirit of the Born-Green-Yvon (BGY) equations of equilibrium statistical mechanics. The resulting equations which contain extra-correlations induced by the activity allow us to determine the stationary density profiles in the presence of external fields, the pair correlations and the pressure of active fluids. In the low density regime we obtained the effective pair potential ϕ(r) acting between two isolated particles separated by a distance, r, showing the existence of an effective attraction between them induced by activity. Based on these results, in the second half of the paper we propose a mean field theory as an approach simpler than the BGY hierarchy and use it to derive a van der Waals expression of the equation of state.

Transport coefficients for driven granular mixtures at low density

The transport coefficients of a granular binary mixture driven by a stochastic bath with friction are determined from the inelastic Boltzmann kinetic equation. A normal solution is obtained via the Chapman-Enskog method for states near homogeneous steady states. The mass, momentum, and heat fluxes are determined to first order in the spatial gradients of the hydrodynamic fields, and the associated transport coefficients are identified. They are given in terms of the solutions of a set of coupled linear integral equations. As in the monocomponent case, since the collisional cooling cannot be compensated for locally by the heat produced by the external driving, the reference distributions (zeroth-order approximations) f(i)((0)) (i=1,2) for each species depend on time through their dependence on the pressure and the temperature. Explicit forms for the diffusion transport coefficients and the shear viscosity coefficient are obtained by assuming the steady-state conditions and by considering the leading terms in a Sonine polynomial expansion. A comparison with previous results obtained for granular Brownian motion and by using a (local) stochastic thermostat is also carried out. The present work extends previous theoretical results derived for monocomponent dense gases [Garzó, Chamorro, and Vega Reyes, Phys. Rev. E 87, 032201 (2013)] to granular mixtures at low density.

Unification of dynamic density functional theory for colloidal fluids to include inertia and hydrodynamic interactions: derivation and numerical experiments

Starting from the Kramers equation for the phase-space dynamics of the N-body probability distribution, we derive a dynamical density functional theory (DDFT) for colloidal fluids including the effects of inertia and hydrodynamic interactions (HI). We compare the resulting theory to extensive Langevin dynamics simulations for both hard rod systems and three-dimensional hard sphere systems with radially symmetric external potentials. As well as demonstrating the accuracy of the new DDFT, by comparing with previous DDFTs which neglect inertia, HI, or both, we also scrutinize the significance of including these effects. Close to local equilibrium we derive a continuum equation from the microscopic dynamics which is a generalized Navier-Stokes-like equation with additional non-local terms governing the effects of HI. For the overdamped limit we recover analogues of existing configuration-space DDFTs but with a novel diffusion tensor.

Charge transport in nanochannels: a molecular theory

We introduce a theoretical and numerical method to investigate the flow of charged fluid mixtures under extreme confinement. We model the electrolyte solution as a ternary mixture comprising two ionic species of opposite charge and a third uncharged component. The microscopy approach is based on kinetic theory and is fully self-consistent. It allows us to determine configurational properties, such as layering near the confining walls and the flow properties. We show that, under the appropriate assumptions, the approach reproduces the phenomenological equations used to describe electrokinetic phenomena, without requiring the introduction of constitutive equations to determine the fluxes. Moreover, we model channels of arbitrary shape and nanometric roughness, features that have important repercussions on the transport properties of these systems. Numerical simulations are obtained by solving the evolution dynamics of the one-particle phase-space distributions of each species by means of a lattice Boltzmann method for flows in straight and wedged channels. Results are presented for the microscopic density, the velocity profiles, and the volumetric and charge flow rates. Strong departures from electroneutrality are shown to appear at the molecular level.

Dynamic density functional theory versus kinetic theory of simple fluids

By combining methods of kinetic and density functional theory, we present a description of molecular fluids which accounts for their microscopic structure and thermodynamic properties as well as their hydrodynamic behavior. We focus on the evolution of the one-particle phase space distribution, rather than on the evolution of the average particle density which features in dynamic density functional theory. The resulting equation can be studied in two different physical limits: diffusive dynamics, typical of colloidal fluids without hydrodynamic interaction where particles are subject to overdamped motion resulting from coupling with a solvent at rest, and inertial dynamics, typical of molecular fluids. Finally, we propose an algorithm to solve numerically and efficiently the resulting kinetic equation by employing a discretization procedure analogous to the one used in the lattice Boltzmann method.

Kinetic theory of correlated fluids: from dynamic density functional to Lattice Boltzmann methods

Using methods of kinetic theory and liquid state theory we propose a description of the nonequilibrium behavior of molecular fluids, which takes into account their microscopic structure and thermodynamic properties. The present work represents an alternative to the recent dynamic density functional theory, which can only deal with colloidal fluids and is not apt to describe the hydrodynamic behavior of a molecular fluid. The method is based on a suitable modification of the Boltzmann transport equation for the phase space distribution and provides a detailed description of the local structure of the fluid and its transport coefficients. Finally, we propose a practical scheme to solve numerically and efficiently the resulting kinetic equation by employing a discretization procedure analogous to the one used in the Lattice Boltzmann method.

Beyond the dynamic density functional theory for steady currents: application to driven colloidal particles in a channel

Motivated by recent studies on the dynamics of colloidal solutions in narrow channels, we consider the steady state properties of an assembly of noninteracting particles subject to the action of a traveling potential moving at a constant speed, while the solvent is modeled by a heat bath at rest in the laboratory frame. Here, since the description we propose takes into account the inertia of the colloidal particles, it is necessary to consider the evolution of both positions and momenta and study the governing equation for the one-particle phase-space distribution. First, we derive the asymptotic form of its solutions as an expansion in Hermite polynomials and their generic properties, such as the force and energy balance, and then we particularize our study to the case of an inverted parabolic potential barrier. We numerically obtain the steady state density and temperature profile and show that the expansion is rapidly convergent for large values of the friction constant and small drifting velocities. On the one hand, the present results confirm the previous studies based on the dynamic density functional theory (DDFT): On the other hand, when the friction constant is large, it display effects such as the presence of a wake behind the barrier and a strong inhomogeneity in the temperature field which are beyond the DDFT description.