No Arabic abstract
Transport coefficients associated with the mass flux of a binary mixture of Maxwell molecules under uniform shear flow are exactly determined from the Boltzmann kinetic equation. A normal solution is obtained via a Chapman--Enskog-like expansion around a local shear flow distribution that retains all the hydrodynamics orders in the shear rate. In the first order of the expansion the mass flux is proportional to the gradients of mole fraction, pressure, and temperature but, due to the anisotropy induced in the system by the shear flow, mutual diffusion, pressure diffusion and thermal diffusion tensors are identified instead of the conventional scalar coefficients. These tensors are obtained in terms of the shear rate and the parameters of the mixture (particle masses, concentrations, and force constants). The description is made both in the absence and in the presence of an external thermostat introduced in computer simulations to compensate for the viscous heating. As expected, the analysis shows that there is not a simple relationship between the results with and without the thermostat. The dependence of the three diffusion tensors on the shear rate is illustrated in the tracer limit case, the results showing that the deviation of the generalized transport coefficients from their equilibrium forms is in general quite important. Finally, the generalized transport coefficients associated with the momentum and heat transport are evaluated from a model kinetic equation of the Boltzmann equation.
We simulate the motion of spherical particles in a phase-separating binary mixture. By combining cell dynamical equations with Langevin dynamics for particles, we show that the addition of hard particles significantly changes both the speed and the morphology of the phase separation. At the late stage of the spinodal decomposition process, particles significantly slow down the domain growth, in qualitative agreement with earlier experimental data.
Population annealing is a sequential Monte Carlo scheme well-suited to simulating equilibrium states of systems with rough free energy landscapes. Here we use population annealing to study a binary mixture of hard spheres. Population annealing is a parallel version of simulated annealing with an extra resampling step that ensures that a population of replicas of the system represents the equilibrium ensemble at every packing fraction in an annealing schedule. The algorithm and its equilibration properties are described and results are presented for a glass-forming fluid composed of a 50/50 mixture of hard spheres with diameter ratio of 1.4:1. For this system, we obtain precise results for the equation of state in the glassy regime up to packing fractions $varphi approx 0.60$ and study deviations from the BMCSL equation of state. For higher packing fractions, the algorithm falls out of equilibrium and a free volume fit predicts jamming at packing fraction $varphi approx 0.667$. We conclude that population annealing is an effective tool for studying equilibrium glassy fluids and the jamming transition.
We consider the steady states of a driven inelastic Maxwell gas consisting of two types of particles with scalar velocities. Motivated by experiments on bilayers where only one layer is driven, we focus on the case when only one of the two types of particles are driven externally, with the other species receiving energy only through inter-particle collision. The velocity $v$ of a particle that is driven is modified to $-r_w v+eta$, where $r_w$ parameterises the dissipation upon the driving and the noise $eta$ is taken from a fixed distribution. We characterize the statistics for small velocities by computing exactly the mean energies of the two species, based on the simplifying feature that the correlation functions are seen to form a closed set of equations. The asymptotic behaviour of the velocity distribution for large speeds is determined for both components through a combination of exact analysis for a range of parameters or obtained numerically to a high degree of accuracy from an analysis of the large moments of velocity. We show that the tails of the velocity distribution for both types of particles have similar behaviour, even though they are driven differently. For dissipative driving ($r_w<1$), the tails of the steady state velocity distribution show non-universal features and depend strongly on the noise distribution. On the other hand, the tails of the velocity distribution are exponential for diffusive driving ($r_w=1$) when the noise distribution decays faster than exponential.
The objective of this study is to assess the impact of a dense-phase treatment on the hydrodynamic description of granular, binary mixtures relative to a previous dilute-phase treatment. Two theories were considered for this purpose. The first, proposed by Garzo and Dufty (GD) [Phys. Fluids {bf 14}, 146 (2002)], is based on the Boltzmann equation which does not incorporate finite-volume effects, thereby limiting its use to dilute flows. The second, proposed by Garzo, Hrenya and Dufty (GHD) [Phys. Rev. E {bf 76}, 31303 and 031304 (2007)], is derived from the Enskog equation which does account for finite-volume effects; accordingly this theory can be applied to moderately dense systems as well. To demonstrate the significance of the dense-phase treatment relative to its dilute counterpart, the ratio of dense (GHD) to dilute (GD) predictions of all relevant transport coefficients and equations of state are plotted over a range of physical parameters (volume fraction, coefficients of restitution, material density ratio, diameter ratio, and mixture composition). These plots reveal the deviation between the two treatments, which can become quite large ($>$100%) even at moderate values of the physical parameters. Such information will be useful when choosing which theory is most applicable to a given situation, since the dilute theory offers relative simplicity and the dense theory offers improved accuracy. It is also important to note that several corrections to original GHD expressions are presented here in the form of a complete, self-contained set of relevant equations.
The Boltzmann equation for d-dimensional inelastic Maxwell models is considered to analyze transport properties in spatially inhomogeneous states close to the simple shear flow. A normal solution is obtained via a Chapman--Enskog--like expansion around a local shear flow distribution f^{(0)} that retains all the hydrodynamic orders in the shear rate. The constitutive equations for the heat and momentum fluxes are obtained to first order in the deviations of the hydrodynamic field gradients from their values in the reference state and the corresponding generalized transport coefficients are {em exactly} determined in terms of the coefficient of restitution alpha and the shear rate a. Since f^{(0)} applies for arbitrary values of the shear rate and is not restricted to weak dissipation, the transport coefficients turn out to be nonlinear functions of both parameters a and alpha. A comparison with previous results obtained for inelastic hard spheres from a kinetic model of the Boltzmann equation is also carried out.