BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Vicente Garzó

Vicente Garzó appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softphysics.flu-dyn

Catalog footprint

What is connected

38works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

38 published item(s)

preprint2025arXiv

Nonlinear transport of tracer particles immersed in a strongly sheared dilute gas with inelastic collisions

The diffusion of tracer particles immersed in a granular gas under uniform shear flow (USF) is analyzed within the framework of the inelastic Boltzmann equation. Two different but complementary approaches are followed to achieve exact results. First, we maintain the structure of the Boltzmann collision operator but consider inelastic Maxwell models (IMM). Using IMM allows us to compute the collisional moments of the Boltzmann operator without knowing the velocity distribution functions of the granular binary mixture explicitly. Second, we consider a kinetic model of the Boltzmann equation for inelastic hard spheres (IHS). This kinetic model is based on the equivalence between a gas of elastic hard spheres subjected to a drag force proportional to the particle velocity and a gas of IHS. We solve the Boltzmann--Lorentz kinetic equation for tracer particles using a generalized Chapman--Enskog--like expansion around the shear flow distribution. This reference distribution retains all hydrodynamic orders in the shear rate. The mass flux is obtained to first order in the deviations of the concentration, pressure, and temperature from their values in the reference state. Due to the velocity space anisotropy induced by the shear flow, the mass flux is expressed in terms of tensorial quantities rather than the conventional scalar diffusion coefficients. The exact results derived here are compared with those previously obtained for IHS by using different approximations [JSTAT P02012 (2007)]. The comparison generally shows reasonable quantitative agreement, especially for IMM results. Finally, we study segregation by thermal diffusion as an application of the theory. The phase diagrams illustrating segregation are shown and compared with IHS results, demonstrating qualitative agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Enskog kinetic theory of binary granular suspensions: heat flux and stability analysis of the homogeneous steady state

The Enskog kinetic theory of multicomponent granular suspensions employed previously [Gómez González, Khalil, and Garzó, Phys. Rev. E \textbf{101}, 012904 (2020)] is considered further to determine the four transport coefficients associated with the heat flux. These transport coefficients are obtained by solving the Enskog equation by means of the application of the Chapman--Enskog method around the local version of the homogeneous state. Explicit forms of the heat flux transport coefficients are provided in steady-state conditions by considering the so-called second Sonine approximation to the distribution function of each species. Their quantitative variation on the control parameters of the mixture (masses and diameters, coefficients of restitution, concentration, volume fraction, and the background temperature) is demonstrated and the results show that in general the dependence of the heat flux transport coefficients on inelasticity is clearly different from that found in the absence of the gas phase (\emph{dry} granular mixtures). As an application of the general results, the stability of the homogeneous steady state is analyzed by solving the linearized Navier--Stokes hydrodynamic equations. The linear stability analysis (which holds for wavelengths long compared with the mean free path) shows that the transversal and longitudinal modes are always stable with respect to long-enough wavelength excitations. This conclusion agrees with previous results derived for monocomponent and (dilute) bidisperse granular suspensions but contrasts with the instabilities found in previous works in dry (no gas phase) granular mixtures.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Kinetic theory of polydisperse granular mixtures: influence of the partial temperatures on transport properties. A review

It is well-recognized that granular media under rapid flow conditions can be modeled as a gas of hard spheres with inelastic collisions. At moderate densities, a fundamental basis for the determination of the granular hydrodynamics is provided by the Enskog kinetic equation conveniently adapted to account for inelastic collisions. A surprising result (compared to its molecular gas counterpart) for granular mixtures is the failure of the energy equipartition, even in homogeneous states. This means that the partial temperatures $T_i$ (measuring the mean kinetic energy of each species) are different to the (total) granular temperature $T$. The goal of this paper is to provide an overview on the effect of different partial temperatures on the transport properties of the mixture. Our analysis addresses first the impact of energy nonequipartition on transport which is only due to the inelastic character of collisions. This effect (which is absent for elastic collisions) is shown to be significant in important problems in granular mixtures such as thermal diffusion segregation. Then, an independent source of energy nonequipartition due to the existence of a divergence of the flow velocity is studied. This effect (which was already analyzed in several pioneering works on dense hard-sphere molecular mixtures) affects to the bulk viscosity coefficient. Analytical (approximate) results are compared against Monte Carlo and molecular dynamics simulations, showing the reliability of kinetic theory for describing granular flows.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the mean square displacement of intruders in freely cooling granular gases

We compute the mean square displacement (MSD) of intruders immersed in a freely cooling granular gas made up of smooth inelastic hard spheres. In general, intruders and particles of the granular gas are assumed to have different mechanical properties, implying that non-equipartition of energy must be accounted for in the computation of the diffusion coefficient $D$. In the hydrodynamic regime, the time decay of the granular temperature $T$ of the cooling granular gas is known to be dictated by Haff's law; the corresponding decay of the intruder's collision frequency entails a time decrease of the diffusion coefficient $D$. Explicit knowledge of this time dependence allows us to determine the MSD by integrating the corresponding diffusion equation. As in previous studies of self-diffusion (intruders mechanically equivalent to gas particles) and the Brownian limit (intruder's mass much larger than the grain's mass), we find a logarithmic time dependence of the MSD as a consequence of Haff's law. Beyond the logarithmic time growth, we find that the MSD depends on the mechanical system parameters in a highly complex way. To explain the observed behaviour, we analyze in detail the intruder's random walk, consisting of ballistic displacements interrupted by anisotropic deflections caused by the collisions with the hard spheres. We also show that the MSD can be thought of as arising from an equivalent random walk with isotropic, uncorrelated steps.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Kinetic theory of granular particles immersed in a molecular gas

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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. As expected, they are given in terms of the solutions of a set of coupled linear integral equations. In addition, in contrast to 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 the 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Enskog kinetic theory of rheology for a moderately dense inertial suspension

The Enskog kinetic theory for moderately dense inertial suspensions under simple shear flow is considered as a model to analyze the rheological properties of the system. The influence of the background fluid on suspended particles is modeled via a viscous drag force plus a Langevin-like term defined in terms of the background temperature. In a previous paper [Hayakawa et al., Phys. Rev. E 96, 042903 (2017)], Grad's moment method with the aid of a linear shear-rate expansion was employed to obtain a theory which gave good agreement with the results of event-driven Langevin simulations of hard spheres for low densities and/or small shear rates. Nevertheless, the previous approach had a limitation of applicability to the high shear-rate and high density regime. Thus, in the present paper, we extend the previous work and develop Grad's theory including higher order terms in the shear rate. This improves significantly the theoretical predictions, a quantitative agreement between theory and simulation being found in the high-density region (volume fractions smaller than or equal to $0.4$).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-Newtonian rheology in inertial suspensions of inelastic rough hard spheres under simple shear flow

Non-Newtonian transport properties of an inertial suspension of inelastic rough hard spheres under simple shear flow are determined from the Boltzmann kinetic equation. The influence of the interstitial gas on rough hard spheres is modeled via a Fokker-Planck generalized equation for rotating spheres accounting for the coupling of both the translational and rotational degrees of freedom of grains with the background viscous gas. The generalized Fokker-Planck term is the sum of two ordinary Fokker-Planck differential operators in linear $\mathbf{v}$ and angular $\boldsymbolω$ velocity space. As usual, each Fokker-Planck operator is constituted by a drag force term (proportional to $\mathbf{v}$ and/or $\boldsymbolω$) plus a stochastic Langevin term defined in terms of the background temperature $T_\text{ex}$. The Boltzmann equation is solved by two different but complementary approaches: (i) by means of Grad's moment method, and (ii) by using a Bhatnagar-Gross-Krook (BGK)-type kinetic model adapted to inelastic rough hard spheres. As occurs in the case of \emph{smooth} inelastic hard spheres, our results show that both the temperature and the non-Newtonian viscosity increase drastically with increasing the shear rate (discontinuous shear thickening effect) while the fourth-degree velocity moments also exhibit an $S$-shape. In particular, while high levels of roughness may slightly attenuate the jump of the viscosity in comparison to the smooth case, the opposite happens for the rotational temperature. As an application of these results, a linear stability analysis of the steady simple shear flow solution is also carried out showing that there are regions of the parameter space where the steady solution becomes linearly unstable.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unified hydrodynamic description for driven and undriven inelastic Maxwell mixtures at low density

A hydrodynamic description for inelastic Maxwell mixtures driven by a stochastic bath with friction is derived. Contrary to previous works where constitutive relations for the fluxes were restricted to states near the homogeneous steady state, here the set of Boltzmann kinetic equations is solved by means of the Chapman--Enskog method by considering a more general time-dependent reference state. Due to this choice, the transport coefficients are given in terms of the solutions of a set of nonlinear differential equations which must be in general numerically solved. The solution to these equations gives the transport coefficients in terms of the parameters of the mixture (masses, diameters, concentration, and coefficients of restitution) and the time-dependent (scaled) parameter $ξ^*$ which determines the influence of the thermostat on the system. The Navier--Stokes transport coefficients are exactly obtained in the special cases of undriven mixtures ($ξ^*=0$) and driven mixtures under steady conditions ($ξ^*=ξ_\text{st}^*$, where $ξ_\text{st}^*$ is the value of the reduced noise strength at the steady state). As a complement, the results for inelastic Maxwell models (IMM) in both undriven and driven steady states are compared against approximate results for inelastic hard spheres (IHS) [Khalil and Garzó, Phys. Rev. E \textbf{88}, 052201 (2013)]. While the IMM predictions for the diffusion transport coefficients show an excellent agreement with those derived for IHS, significant quantitative differences are specially found in the case of the heat flux transport coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Enskog kinetic theory for multicomponent granular suspensions

The Navier--Stokes transport coefficients of multicomponent granular suspensions at moderate densities are obtained in the context of the (inelastic) Enskog kinetic theory. The suspension is modeled as an ensemble of solid particles where the influence of the interstitial gas on grains is via a viscous drag force plus a stochastic Langevin-like term defined in terms of a background temperature. In the absence of spatial gradients, it is shown first that the system reaches a homogeneous steady state where the energy lost by inelastic collisions and viscous friction is compensated for by the energy injected by the stochastic force. Once the homogeneous steady state is characterized, a \emph{normal} solution to the set of Enskog equations is obtained by means of the Chapman--Enskog expansion around the \emph{local} version of the homogeneous state. To first-order in spatial gradients, the Chapman--Enskog solution allows us to identify the Navier--Stokes transport coefficients associated with the mass, momentum, and heat fluxes. In addition, the first-order contributions to the partial temperatures and the cooling rate are also calculated. Explicit forms for the diffusion coefficients, the shear and bulk viscosities, and the first-order contributions to the partial temperatures and the cooling rate are obtained in steady-state conditions by retaining the leading terms in a Sonine polynomial expansion. The results show that the dependence of the transport coefficients on inelasticity is clearly different from that found in its granular counterpart (no gas phase). The present work extends previous theoretical results for \emph{dilute} multicomponent granular suspensions [Khalil and Garzó, Phys. Rev. E \textbf{88}, 052201 (2013)] to higher densities.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tracer diffusion coefficients in a sheared inelastic Maxwell gas

We study the transport properties of an impurity in a sheared granular gas, in the framework of the Boltzmann equation for inelastic Maxwell models. We investigate here the impact of a nonequilibrium phase transition found in such systems, where the tracer species carries a finite fraction of the total kinetic energy (ordered phase). To this end, the diffusion coefficients are first obtained for a granular binary mixture in spatially inhomogeneous states close to the simple shear flow. In this situation, the set of coupled Boltzmann equations are solved by means of a Chapman-Enskog-like expansion around the (local) shear flow distributions for each species, thereby retaining all the hydrodynamic orders in the shear rate $a$. Due to the anisotropy induced by the shear flow, three tensorial quantities $D_{ij}$, $D_{p,ij}$, and $D_{T,ij}$ are required to describe the mass transport process instead of the conventional scalar coefficients. These tensors are given in terms of the solutions of a set of coupled algebraic equations, which can be \emph{exactly} solved as functions of the shear rate $a$, the coefficients of restitution $α_{sr}$ and the parameters of the mixture (masses and composition). Once the forms of $D_{ij}$, $D_{p,ij}$, and $D_{T,ij}$ are obtained for arbitrary mole fraction $x_1=n_1/(n_1+n_2)$ (where $n_r$ is the number density of species $r$), the tracer limit ($x_1\to 0$) is carefully considered for the above three diffusion tensors. Explicit forms for these coefficients are derived showing that their shear rate dependence is significantly affected by the order-disorder transition.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Transport coefficients of solid particles immersed in a viscous gas

Transport properties of a suspension of solid particles in a viscous gas are studied. The dissipation in such systems arises from two sources: inelasticity in particle collisions and viscous dissipation due to the effect of the gas phase on the particles. Here, we consider a simplified case in which the mean relative velocity between the gas and solid phases is taken to be zero, such that "thermal drag" is the only remaining gas-solid interaction. Unlike the previous more general treatment of the drag force [Garzó \emph{et al.}, J. Fluid Mech. \textbf{712}, 129 (2012)], here we take into account contributions to the (scaled) transport coefficients $η^*$ (shear viscosity), $κ^*$ (thermal conductivity) and $μ^*$ (Dufour-like coefficient) coming from the temperature-dependence of the (dimensionless) friction coefficient $γ^*$ characterizing the amplitude of the drag force. At moderate densities, the thermal drag model (which is based on the Enskog kinetic equation) is solved by means of the Chapman-Enskog method and the Navier-Stokes transport coefficients are determined in terms of the coefficient of restitution, the solid volume fraction and the friction coefficient. The results indicate that the effect of the gas phase on $η^*$ and $μ^*$ is non-negligible (especially in the case of relatively dilute systems) while the form of $κ^*$ is the same as the one obtained in the dry granular limit. Finally, as an application of these results, a linear stability analysis of the hydrodynamic equations is carried out to analyze the conditions for stability of the homogeneous cooling state. A comparison with direct numerical simulations shows a good agreement for conditions of practical interest.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Anomalous transport of impurities in inelastic Maxwell gases

A mixture of dissipative hard grains generically exhibits a breakdown of kinetic energy equipartition. The undriven and thus freely cooling binary problem, in the tracer limit where the density of one species becomes minute, may exhibit an extreme form of this breakdown, with the minority species carrying a finite fraction of the total kinetic energy of the system. We investigate the fingerprint of this non-equilibrium phase transition, akin to an ordering process, on transport properties. The analysis, performed by solving the Boltzmann kinetic equation from a combination of analytical and Monte Carlo techniques, hints at the possible failure of hydrodynamics in the ordered region. As a relevant byproduct of the study, the behaviour of the second and fourth-degree velocity moments is also worked out.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generalized transport coefficients for inelastic Maxwell mixtures under shear flow

The Boltzmann equation framework for inelastic Maxwell models is considered to determine the transport coefficients associated with the mass, momentum and heat fluxes of a granular binary mixture in spatially inhomogeneous states close to the simple shear flow. The Boltzmann equation is solved by means of a Chapman-Enskog-like expansion around the (local) shear flow distributions $f_r^{(0)}$ for each species that retain all the hydrodynamic orders in the shear rate. Due to the anisotropy induced by the shear flow, tensorial quantities are required to describe the transport processes instead of the conventional scalar coefficients. These tensors are given in terms of the solutions of a set of coupled equations, which can be analytically solved as functions of the shear rate $a$, the coefficients of restitution $α_{rs}$ and the parameters of the mixture (masses, diameters and composition). Since the reference distribution functions $f_r^{(0)}$ apply for arbitrary values of the shear rate and are not restricted to weak dissipation, the corresponding generalized coefficients turn out to be nonlinear functions of both $a$ and $α_{rs}$. The dependence of the relevant elements of the three diffusion tensors on both the shear rate and dissipation is illustrated in the tracer limit case, the results showing that the deviation of the generalized transport coefficients from their forms for vanishing shear rates is in general significant. A comparison with the previous results obtained analytically for inelastic hard spheres by using Grad's moment method is carried out showing a good agreement over a wide range of values for the coefficients of restitution. Finally, as an application of the theoretical expressions derived here for the transport coefficients, thermal diffusion segregation of an intruder immersed in a granular gas is also studied.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inelastic Maxwell models for monodisperse gas-solid flows

The Boltzmann equation for $d$-dimensional inelastic Maxwell models is considered to analyze transport properties for monodisperse gas-solid suspensions. The influence of the interstitial gas phase on the dynamics of solid particles is modeled via a viscous drag force. The Chapman-Enskog method is applied to solve the inelastic Boltzmann equation to first order in the deviations of the hydrodynamic fields from their values in the homogeneous cooling state. Explicit expressions for the Navier-Stokes transport coefficients are \emph{exactly} obtained in terms of both the coefficient of restitution and the friction coefficient characterizing the amplitude of the external force. The conditions under which a hydrodynamic regime independent of the initial conditions is reached are widely discussed. Finally, the results derived here are compared with those previously obtained for inelastic hard spheres in steady state conditions by using the so-called first Sonine approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stability of freely cooling granular mixtures at moderate densities

The formation of velocity vortices and density clusters is an intriguing phenomenon of freely cooling granular flows. In this work, the critical length scale $L_c$ for the onset of instability is determined via stability analysis of the linearized Navier-Stokes hydrodynamic equations of $d$-dimensional granular binary mixtures at moderate densities. In contrast to previous attempts, the analysis is not restricted to nearly elastic systems since it takes into account the nonlinear dependence of the transport coefficients and the cooling rate on the collisional dissipation. As expected from previous results obtained in the very dilute regime, linear stability shows $d-1$ transversal (shear) modes and a longitudinal ("heat") mode to be unstable with respect to long enough wavelength excitations. The theoretical predictions also show that the origin of the instability is driven by the transversal component of the velocity field that becomes \emph{unstable} when the system length $L>L_c$. An explicit expression of $L_c$ is obtained in terms of the masses and diameters of the mixture, the composition, the volume fraction and the coefficients of restitution. Previous results derived in the limit of both mechanically equivalent particles and low-density mixtures are consistently recovered. Finally, a comparison with previous theoretical works which neglect the influence of dissipation on the transport coefficients shows quantitative discrepancies for strong dissipation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Homogeneous states in driven granular mixtures: Enskog kinetic theory versus molecular dynamics simulations

The homogeneous state of a binary mixture of smooth inelastic hard disks or spheres is analyzed. The mixture is driven by a thermostat composed by two terms: a stochastic force and a drag force proportional to the particle velocity. The combined action of both forces attempts to model the interaction of the mixture with a bath or surrounding fluid. The problem is studied by means of two independent and complementary routes. First, the Enskog kinetic equation with a Fokker-Planck term describing interactions of particles with thermostat is derived. The ratio of kinetic temperatures $T_1/T_2$ and the fourth-degree velocity moments $λ_1$ and $λ_2$ (which measure non-Gaussian properties of $φ_i$) are explicitly determined as a function of the mass ratio, size ratio, composition, density and coefficients of restitution. Secondly, to assess the reliability of the theoretical results, molecular dynamics simulations of a binary granular mixture of spheres are performed for two values of the coefficient of restitution ($α=0.9$ and 0.8) and three different solid volume fractions ($ϕ=0.00785$, 0.1 and 0.2). Comparison between kinetic theory and computer simulations for the temperature ratio shows excellent agreement, even for moderate densities and strong dissipation. In the case of the cumulants $λ_1$ and $λ_2$, good agreement is found for the lower densities although significant discrepancies between theory and simulation are observed with increasing density.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hydrodynamic Burnett equations for inelastic Maxwell models of granular gases

The hydrodynamic Burnett equations and the associated transport coefficients are exactly evaluated for generalized inelastic Maxwell models. In those models, the one-particle distribution function obeys the inelastic Boltzmann equation, with a velocity-independent collision rate proportional to the $γ$ power of the temperature. The pressure tensor and the heat flux are obtained to second order in the spatial gradients of the hydrodynamic fields with explicit expressions for all the Burnett transport coefficients as functions of $γ$, the coefficient of normal restitution, and the dimensionality of the system. Some transport coefficients that are related in a simple way in the elastic limit become decoupled in the inelastic case. As a byproduct, existing results in the literature for three-dimensional elastic systems are recovered, and a generalization to any dimension of the system is given. The structure of the present results is used to estimate the Burnett coefficients for inelastic hard spheres.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hydrodynamic granular segregation induced by boundary heating and shear

Segregation induced by a thermal gradient of an impurity in a driven low-density granular gas is studied. The system is enclosed between two parallel walls from which we input thermal energy to the gas. We study here steady states occurring when the inelastic cooling is exactly balanced by some external energy input (stochastic force or viscous heating), resulting in a uniform heat flux. A segregation criterion based on Navier-Stokes granular hydrodynamics is written in terms of the tracer diffusion transport coefficients, whose dependence on the parameters of the system (masses, sizes and coefficients of restitution) is explicitly determined from a solution of the inelastic Boltzmann equation. The theoretical predictions are validated by means of Monte Carlo and molecular dynamics simulations, showing that Navier-Stokes hydrodynamics produces accurate segregation criteria even under strong shearing and/or inelasticity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Influence of a drag force on linear transport in low-density gases. Stability analysis

The transport coefficients of a dilute classical gas in the presence of a drag force proportional to the velocity of the particle are determined from the Boltzmann equation. The viscous drag force could model the friction of solid particles with a surrounding fluid (interstitial gas phase). First, when the drag force is the only external action on the state of the system, the Boltzmann equation admits a Maxwellian solution $f_0(\mathbf{v},t)$ with a time-dependent temperature. Then, the Boltzmann equation is solved by means of the Chapman-Enskog expansion around the local version of the distribution $f_0$ to obtain the relevant transport coefficients of the system: the shear viscosity $η$, the thermal conductivity $κ$, and a new transport coefficient $μ$ (which is also present in granular gases) relating the heat flux with the density gradient. The results indicate that while $η$ is not affected by the drag force, the impact of this force on the transport coefficients $κ$ and $μ$ may be significant. Finally, a stability analysis of the linear hydrodynamic equations with respect to the time-dependent equilibrium state is performed, showing that the onset of instability is associated with the transversal shear mode that could be unstable for wave numbers smaller than a certain critical wave number.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Instabilities in moderately dense granular binary mixtures

A linear stability analysis of the Navier-Stokes (NS) granular hydrodynamic equations is performed to determine the critical length scale for the onset of vortices and clusters instabilities in granular dense binary mixtures. In contrast to previous attempts, our results (which are based on the solution to the inelastic Enskog equation to NS order) are not restricted to nearly elastic systems since they take into account the complete nonlinear dependence of the NS transport coefficients on the coefficients of restitution $α_{ij}$. The theoretical predictions for the critical length scales are compared to molecular dynamics (MD) simulations in flows of strong dissipation ($α_{ij}\geq 0.7$) and moderate solid volume fractions ($ϕ\leq 0.2$). We find excellent agreement between MD and kinetic theory for the onset of velocity vortices, indicating the applicability of NS hydrodynamics to polydisperse flows even for strong inelasticity, finite density, and particle dissimilarity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Navier-Stokes transport coefficients for driven inelastic Maxwell models

We calculate in this work the Navier-Stokes transport coefficients from the Boltzmann equation for $d$-dimensional inelastic Maxwell models. By granular gas we mean here a low density system of identical spheres that lose a fraction of their kinetic energy after collisions. In the present work, the granular gas is fluidized by the presence of a thermostat that aides the system to reach a steady state. The thermostat is composed by two terms: a random force and a drag force. The combined action of both forces, that act homogeneously on the granular gas, tries to mimic the interaction of the set of particles with a surrounding fluid. The Chapman-Enskog method is applied to solve the inelastic Boltzmann equation to first order in the deviations of the hydrodynamic fields from their values in the homogeneous steady state. Since the collisional cooling cannot be compensated locally for by the heat produced by the driving forces, the reference (zeroth-order) distribution function $f^{(0)}$ depends on time through its dependence on the granular temperature. To simplify the analysis and obtain explicit forms for the transport coefficients, the steady state conditions are considered. A comparison with previous results obtained for inelastic hard spheres is also carried out.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Thermal diffusion segregation of an impurity in a driven granular fluid

We study segregation of an impurity in a driven granular fluid under two types of \emph{steady} states. In the first state, the granular gas is driven by a stochastic volume force field with a Fourier-type profile while in the second state, the granular gas is sheared in such a way that inelastic cooling is balanced by viscous heating. We compare theoretical results derived from a solution of the (inelastic) Boltzmann equation at Navier-Stokes (NS) order with those obtained from the Direct Monte Carlo simulation (DSMC) method and molecular dynamics (MD) simulations. Good agreement is found between theory and simulation, which provides strong evidence of the reliability of NS granular hydrodynamics for these steady states (including the dynamics of the impurity), even at high inelasticities. In addition, preliminary results for thermal diffusion in granular fluids at moderate densitis are also presented. As for dilute gases \cite{VGK14}, excellent agreement is also found in this more general case.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Transport coefficients of a granular gas of inelastic rough hard spheres

The Boltzmann equation for inelastic and rough hard spheres is considered as a model of a dilute granular gas. In this model, the collisions are characterized by constant coefficients of normal and tangential restitution and hence the translational and rotational degrees of freedom are coupled. A normal solution to the Boltzmann equation is obtained by means of the Chapman-Enskog method for states near the homogeneous cooling state. The analysis is carried out to first order in the spatial gradients of the number density, the flow velocity, and the granular temperature. The constitutive equations for the momentum and heat fluxes and for the cooling rate are derived and the associated transport coefficients are expressed in terms of the solutions of linear integral equations. For practical purposes, a first Sonine approximation is used to obtain explicit expressions of the transport coefficients as nonlinear functions of both coefficients of restitution and the moment of inertia. Known results for purely smooth inelastic spheres and perfectly elastic and rough spheres are recovered in the appropriate limits.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A numerical study of the Navier-Stokes transport coefficients for 2D granular hydrodynamics

A numerical study is presented to analyze the thermal mechanisms of unsteady, supersonic granular flow, by means of hydrodynamic simulations of the Navier-Stokes granular equations. For this purpose a paradigmatic problem in granular dynamics such as the Faraday instability is selected. Two different approaches for the Navier-Stokes transport coefficients for granular materials are considered, namely the traditional Jenkins-Richman theory for moderately dense quasi-elastic grains, and the improved Garzó-Dufty-Lutsko theory for arbitrary inelasticity, which we also present here. Both solutions are compared with event-driven simulations of the same system under the same conditions, by analyzing the density, the temperature and the velocity field. Important differences are found between the two approaches leading to interesting implications. In particular, the heat transfer mechanism coupled to the density gradient which is a distinctive feature of inelastic granular gases, is responsible for a major discrepancy in the temperature field and hence in the diffusion mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Diffusion transport coefficients for granular binary mixtures at low density. Thermal diffusion segregation

The mass flux of a low-density granular binary mixture obtained previously by solving the Boltzmann equation by means of the Chapman-Enskog method is considered further. As in the elastic case, the associated transport coefficients $D$, $D_p$ and $D'$ are given in terms of the solutions of a set of coupled linear integral equations which are approximately solved by considering the first and second Sonine approximations. The diffusion coefficients are explicitly obtained as functions of the coefficients of restitution and the parameters of the mixture (masses, diameters and concentration) and their expressions hold for an arbitrary number of dimensions. In order to check the accuracy of the second Sonine correction for highly inelastic collisions, the Boltzmann equation is also numerically solved by means of the direct simulation Monte Carlo (DSMC) method to determine the mutual diffusion coefficient $D$ in some special situations (self-diffusion problem and tracer limit). The comparison with DSMC results reveals that the second Sonine approximation to $D$ improves the predictions made from the first Sonine approximation. We also study the granular segregation driven by a uni-directional thermal gradient. The segregation criterion is obtained from the so-called thermal diffusion factor $Λ$, which measures the amount of segregation parallel to the temperature gradient. The factor $Λ$ is determined here by considering the second-order Sonine forms of the diffusion coefficients and its dependence on the coefficients of restitution is widely analyzed across the parameter space of the system. The results obtained in this paper extend previous works carried out in the tracer limit (vanishing mole fraction of one of the species) by some of the authors of the present paper.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Grad's moment method for a granular fluid at moderate densities. Navier-Stokes transport coefficients

The Navier-Stokes transport coefficients of a granular dense fluid of smooth inelastic hard disks or spheres are explicitly determined by solving the inelastic Enskog equation by means of Grad's moment method. The transport coefficients are explicitly determined as functions of the (constant) coefficient of restitution and the solid volume fraction. In addition, the cooling rate is also calculated to first order in the spatial gradients. The calculations are performed for an arbitrary number of dimensions. The results are not limited to small dissipation and are expected to apply at moderate densities. It is found that the expressions of the Navier-Stokes transport coefficients and the cooling rate agree with those previously obtained from the Chapman-Enskog method by using the leading terms in a Sonine polynomial expansion. This shows the equivalence between both methods for granular fluids in the Navier-Stokes approximation. A comparison with previous results derived from Grad's moment method for inelastic disks and spheres is also carried out.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Homogeneous steady states in a granular fluid driven by a stochastic bath with friction

The homogeneous state of a granular flow of smooth inelastic hard spheres or disks described by the Enskog-Boltzmann kinetic equation is analyzed. The granular gas is fluidized by the presence of a random force and a drag force. The combined action of both forces, that act homogeneously on the granular gas, tries to mimic the interaction of the set of particles with a surrounding fluid. The first stochastic force thermalizes the system, providing for the necessary energy recovery to keep the system in its gas state at all times, whereas the second force allows us to mimic the action of the surrounding fluid viscosity. After a transient regime, the gas reaches a steady state characterized by a \emph{scaled} distribution function $φ$ that does not only depend on the dimensionless velocity $\mathbf{c}\equiv \mathbf{v}/v_0$ ($v_0$ being the thermal velocity) but also on the dimensionless driving force parameters. The dependence of $φ$ and its first relevant velocity moments $a_2$ and $a_3$ (which measure non-Gaussian properties of $φ$) on both the coefficient of restitution $\al$ and the driven parameters is widely investigated by means of the direct simulation Monte Carlo method. In addition, approximate forms for $a_2$ and $a_3$ are also derived from an expansion of $φ$ in Sonine polynomials. The theoretical expressions of the above Sonine coefficients agree well with simulation data, even for quite small values of $α$. Moreover, the third order expansion of the distribution function makes a significant accuracy improvement for larger velocities and inelasticities. Results also show that the non-Gaussian corrections to the distribution function $φ$ are smaller than those observed for undriven granular gases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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 locally for 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 [V. Garzó, M. G. Chamorro, and F. Vega Reyes, Phys. Rev. E \textbf{87}, 032201 (2013)] to granular mixtures at low density.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Assessing a Hydrodynamic Description for Instabilities in Highly Dissipative, Freely Cooling Granular Gases

An intriguing phenomenon displayed by granular flows and predicted by kinetic-theory-based models is the instability known as particle "clustering," which refers to the tendency of dissipative grains to form transient, loose regions of relatively high concentration. In this work, we assess a modified-Sonine approximation recently proposed [Garzó et al., Physica A 376, 94 (2007)] for a granular gas via an examination of system stability. In particular, we determine the critical length scale associated with the onset of two types of instabilities -vortices and clusters- via stability analyses of the Navier-Stokes-order hydrodynamic equations by using the expressions of the transport coefficients obtained from both the standard and the modified-Sonine approximations. We examine the impact of both Sonine approximations over a range of solids fraction ϕ<0.2 for small restitution coefficients e=0.25--0.4, where the standard and modified theories exhibit discrepancies. The theoretical predictions for the critical length scales are compared to molecular dynamics (MD) simulations, of which a small percentage were not considered due to inelastic collapse. Results show excellent quantitative agreement between MD and the modified-Sonine theory, while the standard theory loses accuracy for this highly dissipative parameter space. The modified theory also remedies a (highdissipation) qualitative mismatch between the standard theory and MD for the instability that forms more readily. Furthermore, the evolution of cluster size is briefly examined via MD, indicating that domain-size clusters may remain stable or halve in size, depending on system parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Collisional rates for the inelastic Maxwell model: application to the divergence of anisotropic high-order velocity moments in the homogeneous cooling state

The collisional rates associated with the isotropic velocity moments $<V^{2r}>$ and the anisotropic moments $<V^{2r}V_i>$ and $<V^{2r}(V_iV_j-d^{-1}V^2δ_{ij})>$ are exactly derived in the case of the inelastic Maxwell model as functions of the exponent $r$, the coefficient of restitution $α$, and the dimensionality $d$. The results are applied to the evolution of the moments in the homogeneous free cooling state. It is found that, at a given value of $α$, not only the isotropic moments of a degree higher than a certain value diverge but also the anisotropic moments do. This implies that, while the scaled distribution function has been proven in the literature to converge to the isotropic self-similar solution in well-defined mathematical terms, nonzero initial anisotropic moments do not decay with time. On the other hand, our results show that the ratio between an anisotropic moment and the isotropic moment of the same degree tends to zero.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Enskog kinetic theory for $d$-dimensional dense granular gases

The goal of this note is to provide most of the technical details involved in the application of the Chapman-Enskog method to solve the revised Enskog equation to Navier-Stokes order. Explicit expressions for the transport coefficients and the cooling rate are obtained in terms of the coefficient of restitution and the solid volume fraction by using a new Sonine approach. This new approach consists of replacing, where appropriate in the Chapman-Enskog procedure, the local equilibrium distribution (used in the standard first Sonine approximation) by the homogeneous cooling state distribution. The calculations are performed in an arbitrary number of dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Segregation of an intruder in a heated granular dense gas

A recent segregation criterion [V. Garzó, Phys. Rev. E \textbf{78}, 020301(R) (2008)] based on the thermal diffusion factor $Λ$ of an intruder in a heated granular gas described by the inelastic Enskog equation is revisited. The sign of $Λ$ provides a criterion for the transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE). The present theory incorporates two extra ingredients not accounted for by the previous theoretical attempt. First, the theory is based upon the second Sonine approximation to the transport coefficients of the mass flux of intruder. Second, the dependence of the temperature ratio (intruder temperature over that of the host granular gas) on the solid volume fraction is taken into account in the first and second Sonine approximations. In order to check the accuracy of the Sonine approximation considered, the Enskog equation is also numerically solved by means of the direct simulation Monte Carlo (DSMC) method to get the kinetic diffusion coefficient $D_0$. The comparison between theory and simulation shows that the second Sonine approximation to $D_0$ yields an improvement over the first Sonine approximation when the intruder is lighter than the gas particles in the range of large inelasticity. With respect to the form of the phase diagrams for the BNE/RBNE transition, the kinetic theory results for the factor $Λ$ indicate that while the form of these diagrams depends sensitively on the order of the Sonine approximation considered when gravity is absent, no significant differences between both Sonine solutions appear in the opposite limit (gravity dominates the thermal gradient). In the former case (no gravity), the first Sonine approximation overestimates both the RBNE region and the influence of dissipation on thermal diffusion segregation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermal diffusion segregation in granular binary mixtures described by the Enskog equation

Diffusion induced by a thermal gradient in a granular binary mixture is analyzed in the context of the (inelastic) Enskog equation. Although the Enskog equation neglects velocity correlations among particles which are about to collide, it retains spatial correlations arising from volume exclusion effects and thus it is expected to apply to moderate densities. In the steady state with gradients only along a given direction, a segregation criterion is obtained from the thermal diffusion factor $Λ$ measuring the amount of segregation parallel to the thermal gradient. As expected, the sign of the factor $Λ$ provides a criterion for the transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE) by varying the parameters of the mixture (masses, sizes, concentration, solid volume fraction, and coefficients of restitution). The form of the phase diagrams for the BNE/RBNE transition is illustrated in detail for several systems, with special emphasis on the significant role played by the inelasticity of collisions. In particular, an effect already found in dilute gases (segregation in a binary mixture of identical masses and sizes {\em but} different coefficients of restitution) is extended to dense systems. A comparison with recent computer simulation results shows a good qualitative agreement at the level of the thermal diffusion factor. The present analysis generalizes to arbitrary concentration previous theoretical results derived in the tracer limit case.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Transport coefficients of driven granular fluids at moderate volume fractions

In a recent publication [Phys. Rev. E \textbf{83}, 011301 (2011)], Vollmayr--Lee \emph{et al.} have determined by computer simulations the thermal diffusivity and the longitudinal viscosity coefficients of a driven granular fluid of hard spheres at intermediate volume fractions. Although they compare their simulation results with the predictions of kinetic theory, they use the dilute expressions for the driven system and the modified Sonine approximations for the undriven system. The goal here is to carry out this comparison by proposing a modified Sonine approximation to the Enskog equation for driven granular fluids that leads to a better quantitative agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Energy production rates in fluid mixtures of inelastic rough hard spheres

The aim of this work is to explore the combined effect of polydispersity and roughness on the partial energy production rates and on the total cooling rate of a granular fluid mixture. We consider a mixture of inelastic rough hard spheres of different number densities, masses, diameters, moments of inertia, and mutual coefficients of normal and tangential restitution. Starting from the first equation of the BBGKY hierarchy, the collisional energy production rates associated with the translational and rotational temperatures ($T_i^\text{tr}$ and $T_i^\text{rot}$) are expressed in terms of two-body average values. Next, those average values are estimated by assuming a velocity distribution function based on maximum-entropy arguments, allowing us to express the energy production rates and the total cooling rate in terms of the partial temperatures and the parameters of the mixture. Finally, the results are applied to the homogeneous cooling state of a binary mixture and the influence of inelasticity and roughness on the temperature ratios $T_1^\text{tr}/T_1^\text{rot}$, $T_2^\text{tr}/T_1^\text{tr}$, and $T_2^\text{rot}/T_1^\text{rot}$ is analyzed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Non-Newtonian granular hydrodynamics. What do the inelastic simple shear flow and the elastic Fourier flow have in common?

We describe a special class of steady Couette flows in dilute granular gases admitting a non-Newtonian hydrodynamic description for strong dissipation. The class occurs when viscous heating exactly balances inelastic cooling, resulting in a uniform heat flux. It includes the Fourier flow of ordinary gases and the simple or uniform shear flow (USF) of granular gases as special cases. The rheological functions have the same values as in the USF and generalized thermal conductivity coefficients can be identified. These points are confirmed by molecular dynamics simulations, Monte Carlo simulations of the Boltzmann equation, and analytical results from Grad's 13-moment method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Segregation by thermal diffusion in granular shear flows

Segregation by thermal diffusion of an intruder immersed in a sheared granular gas is analyzed from the (inelastic) Boltzmann equation. Segregation is induced by the presence of a temperature gradient orthogonal to the shear flow plane and parallel to gravity. We show that, like in analogous systems without shear, the segregation criterion yields a transition between upwards segregation and downwards segregation. The form of the phase diagrams is illustrated in detail showing that they depend sensitively on the value of gravity relative to the thermal gradient. Two specific situations are considered: i) absence of gravity, and ii) homogeneous temperature. We find that both mechanisms (upwards and downwards segregation) are stronger and more clearly separated when compared with segregation criteria in systems without shear.

0 citations0 reviewsNo code linkedOpen access