Lattice-gas and lattice-Boltzmann models of miscible fluids

We introduce new lattice-gas and lattice-Boltzmann models for simulating miscible fluids in two dimensions. The inclusion of a nonlocal interaction produces a lattice gas with lower diffusivity than achieved before. To overcome some observed unphysical properties of this lattice gas, we introduce a lattice-Boltzmann analogue of the model. We first formulate a miscible two-component lattice-Boltzmann model with local interactions only, and show that its diffusivity is determined by an eigenvalue of the linearized collision operator. Diffusivity is then reduced by including nonlocal interactions. The utility of the model is demonstrated by a simulation of two-dimensional viscous fingering.     

Shear Stress in a Couette Flow of Liquid-Particle Suspensions

The mechanisms of momentum transfer and shear stress of liquid-particle suspensions in two-dimensional Couette flow are studied using direct numerical simulation by lattice-Boltzmann techniques. The results obtained display complex flow phenomena that arise from the two-phase nature of the fluid including a nonlinear velocity profile, layering of particles, and apparent slip near the solid walls. The general rheological behaviour of the suspension is dilatant. A detailed study of the various momentum transfer mechanisms that contribute to the total shear stress indicates that the observed shear thickening is related to enhanced relative solid phase stress for increasing shear rates.     

Dense colloidal suspensions under time-dependent shear

We consider the nonlinear rheology of dense colloidal suspensions under a time-dependent simple shear flow. Starting from the Smoluchowski equation for interacting Brownian particles advected by shearing (ignoring fluctuations in fluid velocity), we develop a formalism which enables the calculation of time-dependent, far-from-equilibrium averages. Taking shear stress as an example, we derive exactly a generalized Green-Kubo relation and an equation of motion for the transient density correlator, involving a three-time memory function. Mode coupling approximations give a closed constitutive equation yielding the time-dependent stress for arbitrary shear rate history. We solve this equation numerically for the special case of a hard sphere glass subject to step strain.     

Kubo formulas for relativistic fluids in strong magnetic fields

Magnetohydrodynamics of strongly magnetized relativistic fluids is derived in the ideal and dissipative cases, taking into account the breaking of spatial symmetries by a quantizing magnetic field. A complete set of transport coefficients, consistent with the Curie and Onsager principles, is derived for thermal conduction, as well as shear and bulk viscosities. It is shown that in the most general case the dissipative function contains five shear viscosities, two bulk viscosities, and three thermal conductivity coefficients. We use Zubarev’s non-equilibrium statistical operator method to relate these transport coefficients to correlation functions of the equilibrium theory. The desired relations emerge at linear order in the expansion of the non-equilibrium statistical operator with respect to the gradients of relevant statistical parameters (temperature, chemical potential, and velocity.) The transport coefficients are cast in a form that can be conveniently computed using equilibrium (imaginary-time) infrared Green’s functions defined with respect to the equilibrium statistical operator.     

Lattice-Boltzmann Simulation of Particle Suspensions in Shear Flow

Inclusion of short-range particle–particle interactions for increased numerical stability in a lattice-Boltzmann code for particle-fluid suspensions, and handling of the particle phase for an effective implementation of the code for parallel computing, are discussed and formulated. In order to better understand the origin of the shear-thickening behavior observed in real suspensions, two simplified cases are considered with the code thus developed. A chain-like cluster of suspended particles is shown to increase the momentum transfer in a shear flow between channel walls, and thereby the effective viscosity of the suspension in comparison with random configurations of particles. A single suspended particle is also shown to increase the effective viscosity under shear flow of this simple suspension for particle Reynolds numbers above unity, due to inertial effects that change the flow configuration around the particle. These mechanisms are expected to carry over to large-scale particle-fluid suspensions.     

Viscoelastic model of phase separation in colloidal suspensions and emulsions

We propose a simple physical model of phase separation of colloidal suspensions and emulsions, which we call the "viscoelastic model." On the basis of this model, we consider two poorly understood phenomena: (i) phase separation accompanying the formation of a transient gel, and its collapse, and (ii) shear effects on composition fluctuations and phase separation. These phenomena can be explained by "asymmetric stress division" between the components of a mixture due to their size difference; the interaction network of particles can store elastic energy, while a fluid component cannot. The importance of the bulk stress stemming from an interaction network is discussed, using a concept of self-induced elastic constraint due to connectivity. We argue that there are common features to polymer solutions, colloidal suspensions, emulsions, and possibly protein solutions. They originate from dynamic asymmetry between the components and the resulting interaction network of the slower component of a mixture, which leads to the formation of a transient gel.     

A Conservative,Entropic Multispecies BGK Model

We derive a conservative multispecies BGK model that follows the spirit of the original, single species BGK model by making the specific choice to conserve species masses, total momentum, and total kinetic energy and to satisfy Boltzmann’s \(\mathcal {H}\)-Theorem. The derivation emphasizes the connection to the Boltzmann operator which allows for direct inclusion of information from higher-fidelity collision physics models. We also develop a complete hydrodynamic closure via the Chapman-Enskog expansion, including a general procedure to generate symmetric diffusion coefficients based on this model. We numerically investigate velocity and temperature relaxation in dense plasmas and compare the model with previous multispecies BGK models and discuss the trade-offs that are made in defining and using them. In particular, we demonstrate that the BGK model in the NRL plasma formulary does not conserve momentum or energy in general.     

Normal modes of a relativistic quantum plasma; The one-component plasma

The normal modes of a relativistic electron gas are studied on the basis of the Boltzmann-Vlasov kinetic equation via a projection operator formalism. A general framework is constructed in which the fully relativistic Vlasov self-consistent force term appears as a symmetric operator acting in the Hilbert space of one-particle states. The plasma-dynamical equations are obtained by projecting onto the subspace consisting of the charge, energy and momentum densities, plus the nonconserved current density. The eigenmodes of these equations include two transverse and two longitudinal plasma modes, and one damped heat mode. They are explicitly calculated up to second order in the wave vector and to first order in the collision frequency.     

Self-diffusion in sheared colloidal suspensions: violation of fluctuation-dissipation relation

Using memory-function formalism we show that in sheared colloidal suspensions the fluctuation-dissipation theorem for self-diffusion, i.e., Einstein's relation between self-diffusion and mobility tensors, is violated and propose a new way to measure this violation in Brownian dynamics simulations. We derive mode-coupling expressions for the tagged particle friction tensor and for an effective, shear-rate dependent temperature.     

Gyrokinetic Fokker-Planck collision operator

The gyrokinetic linearized exact Fokker-Planck collision operator is obtained in a form suitable for plasma gyrokinetic equations, for arbitrary mass ratio. The linearized Fokker-Planck operator includes both the test-particle and field-particle contributions, and automatically conserves particles, momentum, and energy, while ensuring non-negative entropy production. Finite gyroradius effects in both field-particle and test-particle terms are evaluated. When implemented in gyrokinetic simulations, these effects can be precomputed. The field-particle operator at each time step requires the evaluation of a single two-dimensional integral, and is not only more accurate, but appears to be less expensive to evaluate than conserving model operators.