Showing papers by "Vicente Garzó published in 2021"

Mpemba-like effect in driven binary mixtures

Abstract: The Mpemba effect occurs when two samples at different initial temperatures evolve in such a way that the temperatures cross each other during the relaxation toward equilibrium. In this paper, we show the emergence of a Mpemba-like effect in a molecular binary mixture in contact with a thermal reservoir (bath). The interaction between the gaseous particles of the mixture and the thermal reservoir is modeled via a viscous drag force plus a stochastic Langevin-like term. The presence of the external bath couples the time evolution of the total and partial temperatures of each component allowing the appearance of the Mpemba phenomenon, even when the initial temperature differences are of the same order of the temperatures themselves. Analytical results are obtained by considering multitemperature Maxwellian approximations for the velocity distribution functions of each component. The theoretical analysis is carried out for initial states close to and far away (large Mpemba-like effect) from equilibrium. The former situation allows us to develop a simple theory where the time evolution equation for the temperature is linearized around its asymptotic equilibrium solution. This linear theory provides an expression for the crossover time. We also provide a qualitative description of the large Mpemba effect. Our theoretical results agree very well with computer simulations obtained by numerically solving the Enskog kinetic equation by means of the direct simulation Monte Carlo method and by performing molecular dynamics simulations. Finally, preliminary results for driven granular mixtures also show the occurrence of a Mpemba-like effect for inelastic collisions.

Time-dependent homogeneous states of binary granular suspensions

Abstract: The time evolution of a homogeneous bidisperse granular suspension is studied in the context of the Enskog kinetic equation. The influence of the surrounding viscous gas on the solid particles is modeled via a deterministic viscous drag force plus a stochastic Langevin-like term. It is found first that, regardless of the initial conditions, the system reaches (after a transient period lasting a few collisions per particle) a universal unsteady hydrodynamic regime where the distribution function of each species not only depends on the dimensionless velocity (as in the homogeneous cooling state) but also on the instantaneous temperature scaled with respect to the background temperature. To confirm this result, theoretical predictions for the time-dependent partial temperatures are compared against direct simulation Monte Carlo (DSMC) results; the comparison shows an excellent agreement confirming the applicability of hydrodynamics in granular suspensions. Also, in the transient regime, the so-called Mpemba-like effect (namely, when an initially hotter sample cools sooner than the colder one) is analyzed for inelastic collisions. The theoretical analysis of the Mpemba effect is performed for initial states close to and far away from the asymptotic steady state. In both cases, good agreement is found again between theory and DSMC results. As a complement to the previous studies, we determine in this paper the dependence of the steady values of the dynamic properties of the suspension on the parameter space of the system. More specifically, we focus our attention on the temperature ratio T1/T2 and the fourth degree cumulants c1 and c2 (measuring the departure of the velocity distributions f1 and f2 from their Maxwellian forms). While our approximate theoretical expression for T1/T2 agrees very well with computer simulations, some discrepancies are found for the cumulants. Finally, a linear stability analysis of the steady state solution is also carried out showing that the steady state is always linearly stable.

Time-dependent homogeneous states of binary granular suspensions

Abstract: The time evolution of a homogeneous bidisperse granular suspension is studied in the context of the Enskog kinetic equation. The influence of the surrounding viscous gas on the solid particles is modeled via a deterministic viscous drag force plus a stochastic Langevin-like term. It is found first that, regardless of the initial conditions, the system reaches (after a transient period lasting a few collisions per particle) a universal unsteady hydrodynamic regime where the distribution function of each species not only depends on the dimensionless velocity (as in the homogeneous cooling state) but also on the instantaneous temperature scaled with respect to the background temperature. To confirm this result, theoretical predictions for the time-dependent partial temperatures are compared against direct simulation Monte Carlo (DSMC) results; the comparison shows an excellent agreement confirming the applicability of hydrodynamics in granular suspensions. Also, in the transient regime, the so-called Mpemba-like effect (namely, when an initially hotter sample cools sooner than the colder one) is analyzed for inelastic collisions. The theoretical analysis of the Mpemba effect is performed for initial states close to and far away from the asymptotic steady state. In both cases, a good agreement is found again between theory and DSMC results. As a complement of the previous studies, we determine in this paper the dependence of the steady values of the dynamic properties of the suspension on the parameter space of the system. More specifically, we focus on our attention in the temperature ratio $T_1/T_2$ and the fourth degree cumulants $c_1$ and $c_2$ (measuring the departure of the velocity distributions $f_1$ and $f_2$ from their Maxwellian forms). Finally, a linear stability analysis of the steady state solution is also carried out showing that the steady state is always linearly stable.

Navier–Stokes transport coefficients for a model of a confined quasi-two-dimensional granular binary mixture

Abstract: The Navier–Stokes transport coefficients for a model of a confined quasi-two-dimensional granular binary mixture of inelastic hard spheres are determined from the Boltzmann kinetic equation. A normal or hydrodynamic solution to the Boltzmann equation is obtained via the Chapman–Enskog method for states near the local version of the homogeneous time-dependent state. 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. In addition, in contrast to the previous results obtained for low-density granular mixtures, there are also nonzero contributions to the first-order approximations to the partial temperatures T i ( 1 ) and the cooling rate ζ(1). Explicit forms for the diffusion transport coefficients, the shear viscosity coefficient, and the quantities T i ( 1 ) and ζ(1) are obtained by assuming steady state conditions and by considering the leading terms in a Sonine polynomial expansion. The above transport coefficients are given in terms of the coefficients of restitution, concentration, and the masses and diameters of the components of the mixture. The results apply, in principle, for arbitrary degree of inelasticity and are not limited to specific values of concentration, mass, and/or size ratios. As a simple application of these results, the violation of the Onsager reciprocal relations for a confined granular mixture is quantified in terms of the parameter space of the problem.

Comment on "Kinetic theory models for granular mixtures with unequal granular temperature: Hydrodynamic velocity" [Phys. Fluids 33, 043321 (2021)]

Abstract: Comment on the paper J. Solsvik and E. Manger, "Kinetic theory models for granular mixtures with unequal granular temperature: Hydrodynamic velocity," Phys. Fluids \textbf{33}, 043321 (2021).

Comment on "Kinetic theory models for granular mixtures with unequal granular temperature: Hydrodynamic velocity" [Phys. Fluids 33, 043321 (2021)]

Abstract: Comment on the paper J. Solsvik and E. Manger, "Kinetic theory models for granular mixtures with unequal granular temperature: Hydrodynamic velocity," Phys. Fluids \textbf{33}, 043321 (2021).

Non-monotonic Mpemba effect in binary molecular suspensions

Abstract: The Mpemba effect is a phenomenon in which an initially hotter sample cools sooner. In this paper, we show the emergence of a non-monotonic Mpemba-like effect in a molecular binary mixture immersed in a viscous gas. Namely, a crossover in the temperature evolution when at least one of the samples presents non-monotonic relaxation. The influence of the bath on the dynamics of the particles is modeled via a viscous drag force plus a stochastic Langevin-like term. Each component of the mixture interchanges energy with the bath depending on the mechanical properties of its particles. This discrimination causes the coupling between the time evolution of temperature with that of the partial temperatures of each component. The non-monotonic Mpemba effect—and its inverse and mixed counterparts—stems from this coupling. In order to obtain analytical results, the velocity distribution functions of each component are approximated by considering multitemperature Maxwellian distributions. The theoretical results derived from the Enskog kinetic theory show an excellent agreement with direct simulation Monte Carlo (DMSC) data.

Kinetic theory of granular particles immersed in a molecular gas

Abstract: The transport coefficients of a dilute gas of inelastic hard spheres immersed in a molecular gas are determined. We assume that the number density of the granular gas is much smaller than that of the surrounding molecular gas, so that the latter is not affected by the presence of solid particles. In this situation, the molecular gas may be treated as a thermostat (or bath) of elastic hard spheres at a fixed temperature. This system (granular gas thermostated by a bath of elastic hard spheres) can be considered as a reliable model for describing the dynamic properties of particle-laden suspensions. The Boltzmann kinetic equation is the starting point of the present work. First step is to characterise the reference state in the perturbation scheme, namely the homogeneous state. Theoretical results for the granular temperature and kurtosis obtained in the homogeneous steady state are compared against Monte Carlo simulations showing a good agreement. Then, the Chapman-Enskog method is employed to solve the Boltzmann equation to first order in spatial gradients. As expected, the Navier-Stokes-Fourier transport coefficients of the granular gas are given in terms of the solutions of a coupled set of linear integral equations which are approximately solved by considering the leading terms in a Sonine polynomial expansion. Our results show that the dependence of the transport coefficients on the coefficient of restitution is quite different from that found when the influence of the interstitial molecular gas is neglected (dry granular gas). When the granular particles are much more heavier than the gas particles (Brownian limit) the expressions of the transport coefficients are consistent with those previously derived from the Fokker-Planck equation. Finally, a linear stability analysis of the homogeneous steady state is performed showing this state is always linearly stable.

Stability of the homogeneous steady state for a model of a confined quasi-two-dimensional granular fluid

Abstract: A linear stability analysis of the hydrodynamic equations of a model for confined quasi-two-dimensional granular gases is carried out. The stability analysis is performed around the homogeneous steady state (HSS) reached eventually by the system after a transient regime. In contrast to previous studies (which considered dilute or quasielastic systems), our analysis is based on the results obtained from the inelastic Enskog kinetic equation, which takes into account the (nonlinear) dependence of the transport coefficients and the cooling rate on dissipation and applies to moderate densities. As in earlier studies, the analysis shows that the HSS is linearly stable with respect to long enough wavelength excitations.