Open AccessArticle

Exact Results for Non-Newtonian Transport Properties in Sheared Granular Suspensions: Inelastic Maxwell Models and BGK-Type Kinetic Model

Rubén Gómez González

^1,†

and

Vicente Garzó

^2,*,†,‡

Departamento de Física, Universidad de Extremadura, Avda. de Elvas s/n, E-06006 Badajoz, Spain

Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, Avda. de Elvas s/n, E-06006 Badajoz, Spain

Author to whom correspondence should be addressed.

^†

These authors contributed equally to this work.

^‡

https://fisteor.cms.unex.es/investigadores/vicente-garzo-puertos/.

Entropy 2024, 26(3), 265; https://doi.org/10.3390/e26030265

Submission received: 21 February 2024 / Revised: 12 March 2024 / Accepted: 13 March 2024 / Published: 15 March 2024

(This article belongs to the Special Issue 180th Anniversary of Ludwig Boltzmann)

Download

Browse Figures

Versions Notes

Abstract

The Boltzmann kinetic equation for dilute granular suspensions under simple (or uniform) shear flow (USF) is considered to determine the non-Newtonian transport properties of the system. In contrast to previous attempts based on a coarse-grained description, our suspension model accounts for the real collisions between grains and particles of the surrounding molecular gas. The latter is modeled as a bath (or thermostat) of elastic hard spheres at a given temperature. Two independent but complementary approaches are followed to reach exact expressions for the rheological properties. First, the Boltzmann equation for the so-called inelastic Maxwell models (IMM) is considered. The fact that the collision rate of IMM is independent of the relative velocity of the colliding spheres allows us to exactly compute the collisional moments of the Boltzmann operator without the knowledge of the distribution function. Thanks to this property, the transport properties of the sheared granular suspension can be exactly determined. As a second approach, a Bhatnagar–Gross–Krook (BGK)-type kinetic model adapted to granular suspensions is solved to compute the velocity moments and the velocity distribution function of the system. The theoretical results (which are given in terms of the coefficient of restitution, the reduced shear rate, the reduced background temperature, and the diameter and mass ratios) show, in general, a good agreement with the approximate analytical results derived for inelastic hard spheres (IHS) by means of Grad’s moment method and with computer simulations performed in the Brownian limiting case (

m / m_{g} \to \infty

, where

m_{g}

and m are the masses of the particles of the molecular and granular gases, respectively). In addition, as expected, the IMM and BGK results show that the temperature and non-Newtonian viscosity exhibit an S shape in a plane of stress–strain rate (discontinuous shear thickening, DST). The DST effect becomes more pronounced as the mass ratio

m / m_{g}

increases.

Keywords:

granular suspensions; Boltzmann kinetic equation; inelastic Maxwell models; BGK-type kinetic model; non-Newtonian transport properties

1. Introduction

A very usual way of assessing the effect of the surrounding fluid on the dynamics properties of solid particles is through an effective fluid–solid force [1,2,3,4]. In some models, this force is simply proportional to the velocity particle (Stokes linear drag law) [5,6,7,8,9,10,11,12]. This type of force attempts to mimic the energy dissipated by grains due to their friction on the interstitial viscous gas. A more sophisticated model [13] also incorporates a Langevin-like stochastic term that accounts for the energy transferred to grains due to their “collisions” with particles of the background gas. However, although this coarse-grained approach has provided reliable results in the past, it would be desirable to consider a suspension model that takes into account the real collisions between grains and particles of the surrounding (molecular) gas. This sort of suspension model (which was inspired by a previous work of Biben et al. [14]) has been recently proposed [15]. In this model, granular particles are assumed to be sufficiently rarefied so that they do not disturb the state of the molecular (background) gas. As a consequence, the interstitial gas may be treated as a thermostat at the temperature

T_{g}

. Moreover, although the concentration (mole fraction) of grains is quite small, apart from the elastic collisions between solid and molecular gas particles, one has to consider the inelastic collisions between grains themselves. This model can be useful to analyze transport properties in particle-laden suspensions [16] where very dilute particles (like aerosols) are immersed in a fluid (like air).

The rheological properties of a granular suspension under simple (or uniform) shear flow (USF) have been recently determined [17]. In contrast to previous attempts [5,6,7,8,10,11,12,18], the results obtained in Ref. [17] were derived from the collisional model proposed in Ref. [15]. On the other hand, a limitation of these results is that they were approximately obtained by employing Grad’s moment method [19], namely, a method based on the truncation of a series expansion of the velocity distribution function in (orthogonal) Sonine polynomials. The use of this approximate method is essentially motivated by the form of the collision rate for inelastic hard spheres (IHS) appearing inside the Boltzmann collision operator. The collision rate for IHS is proportional to the magnitude of the normal component of the relative velocity of the two spheres that are about to collide. This velocity dependence of the collision rate for IHS prevents the possibility of deriving exact expressions for the transport properties in the USF problem, even in the case of elastic collisions.

A possible way to overcome the technical difficulty of the hard-sphere kernel is to consider the so-called inelastic Maxwell models (IMM). As for the conventional Maxwell molecules [20,21,22], the collision rate of IMM is independent of the relative velocity of the two colliding spheres [23]. The use of IMM instead of IHS opens up the possibility of obtaining exact analytical results of the Boltzmann equation in some specific nonequilibrium situations, like the USF. In particular, the knowledge of the collisional moments of the Boltzmann equation for IMM enables a clear exploration of the impact of inelasticity on the non-Newtonian transport properties of the granular suspension without introducing uncontrolled approximations.

Another possible alternative for obtaining exact results is to consider a kinetic model that retains the relevant physical properties of the Boltzmann collision operator but turns out to be more tractable than the true kinetic equation. This route has been widely employed in the past in the case of molecular dilute gases where it has been shown that several exact solutions in far from equilibrium states agree very well with Monte Carlo simulations of the Boltzmann equation. Here, we will consider a kinetic model for granular suspensions [24] to complement the theoretical expressions obtained from the Boltzmann equation for IMM. Since this kinetic model is based on the well known Bhatnagar–Gross–Krook (BGK) model [22] for molecular gases, we will refer to it as a BGK-type kinetic model.

The objective of this paper is to determine the rheological properties of granular particles immersed in a bath of elastic hard spheres under USF. At a macroscopic level, the USF is characterized by constant number densities for solid and gas particles, a uniform temperature, and a (common) linear velocity profile

U_{g, x} = U_{x} = a y

, where a is the constant shear rate. Here,

U_{g}

and

U

denote the mean flow velocities of the molecular and granular gases, respectively. Since we are interested here in the steady state where the system admits a non-Newtonian hydrodynamic description, an external thermostat force (proportional to the peculiar velocity) must be introduced to keep the temperature

T_{g}

of the molecular gas constant.

The use of an IMM as well as a BGK-type kinetic model allows us to exactly compute the rheological properties of the granular suspension. These properties are expressed as nonlinear functions of the (reduced) shear rate

a^{*} = a / γ

(where

γ

is a drift coefficient characterizing the friction of solid particles on the viscous gas), the coefficient of restitution, the (reduced) background temperature

T_{g}^{*}

, and the diameter

σ / σ_{g}

and mass

m / m_{g}

ratios. Here,

σ_{g}

and

m_{g}

are the diameter and mass of the particles of the molecular gas, respectively, while

σ

and m are the diameter and mass of the solid particles, respectively.

A very interesting phenomenon appearing in gas–solid suspensions is the so-called discontinuous shear thickening (DST), where the non-Newtonian shear viscosity of the granular suspension drastically increases with the shear rate. In most of the cases, the DST phenomena occurs in highly concentrated suspensions of particles such as mixtures of cornstarch in water [25]. Although several mechanisms have been proposed in the granular literature [25] to explain this nonequilibrium discontinuous transition, the problem is still not understood.

As occurs for IHS [17], our results show that the kinetic granular temperature and the non-Newtonian viscosity exhibit a DST effect for sufficiently large values of the mass ratio

m / m_{g}

. In fact, in the Brownian limiting case (

m / m_{g} \to \infty

), the expressions of the rheological properties derived here reduce to those previously obtained [26] from a coarse-grained description based on the Fokker–Planck operator. This agreement justifies the use of this latter approach to analyze the DST effect in dilute granular suspensions [27]. It is important to remark that the DST effect found here occurs in a system with a simple structure (low-density granular gas immersed in a dilute molecular gas). This means that the origin of this phenomenon is essentially associated here with the nonlinear response of the granular suspension to the presence of large shear rates rather than to the complex structure of the system.

Apart from the transport properties (which are related to the second-degree velocity moments), the explicit forms of the higher-degree velocity moments as well as the velocity distribution function of the granular gas were also obtained from the BGK model. This is one of the main advantages of using a kinetic model instead of the true Boltzmann equation. Our results show, in particular, that the fourth-degree moments of the distribution function also exhibit a DST effect. With respect to the velocity distribution function, as expected, we find that its distortion from equilibrium is more significant as both the mass

m / m_{g}

and diameter

σ / σ_{g}

ratios depart from 1. In addition, a comparison between the BGK results and numerical solutions of the Boltzmann equation from the direct simulation Monte Carlo (DSMC) method [28] for IHS shows a generally good qualitative agreement between both approaches.

The plan of the paper is as follows. The Boltzmann kinetic equation for a granular gas immersed in a bath of elastic hard spheres under USF is presented in Section 2. The balance equations for the temperatures of the molecular and granular gases are also displayed. Section 3 deals with the calculations carried out for IMM of the rheological properties of the granular suspension. While a shear thinning effect is always found for the nonlinear shear viscosity of the molecular gas, the corresponding shear viscosity of the granular gas exhibits a DST effect for sufficiently large values of the mass ratio

m / m_{g}

. The results derived from the BGK-type kinetic model are provided in Section 4, while a comparison between the theoretical results obtained for IHS, IMM, and BGK model is displayed in Section 5 for several systems. Our results highlight a good agreement for the rheology between the three different approaches. Moreover, theoretical results obtained from the IMM and BGK model are also compared against computer simulations in the Brownian limit (

m / m_{g} \to \infty

), showing a good agreement. The paper is closed in Section 6 with some concluding remarks.

2. Boltzmann Kinetic Equation for Sheared Granular Suspensions

We consider a set of solid particles (granular gas) of mass m and diameter

σ

which are immersed in a solvent (molecular gas) constituted by particles of mass

m_{g}

and diameter

σ_{g}

. As usual, the granular gas is modeled as a gas of hard disks (

d = 2

) or spheres (

d = 3

) with inelastic collisions. In the simplest model, the inelasticity of collisions is characterized by a constant (positive) coefficient of normal restitution

α \leq 1

, where

α = 1

refers to elastic collisions. On the other hand, collisions between solid particles and particles of the molecular gas are elastic. We also assume that the number density of grains is much smaller than that of solvent so that the state of the latter is not perturbed by the presence of the former. In these conditions, we can treat the molecular gas as a bath or thermostat at the temperature

T_{g}

(once the parameters of the system, specifically the shear rate, have been set). Moreover, although the granular gas is sufficiently rarefied, we take into account the collisions among grains in its corresponding Boltzmann kinetic equation.

We assume that the system (granular particles plus solvent) is subjected to USF. As said in the Introduction section, this state is characterized by constant densities

n_{g}

and n, uniform temperatures

T_{g}

and T, and by a (common) linear profile of the x component of the flow velocities along the y axis:

n_{g} \equiv const, n \equiv const,

(1)

\nabla T_{g} = \nabla T = 0,

(2)

U_{g, i} = U_{i} = a_{i j} r_{j}, a_{i j} = a δ_{i x} δ_{j y},

(3)

with a being the constant shear rate. Here,

n_{g}

U_{g}

, and

T_{g}

are the number density, the mean flow velocity, and the temperature, respectively, of the molecular gas. In terms of its one-particle velocity distribution function

f_{g} (r, v; t)

, these hydrodynamic fields are defined as

\{n_{g}, n_{g} U_{g}, d n_{g} T_{g}\} = \int d v \{1, v, m_{g} V^{2}\} f_{g} (v),

(4)

where

V = v - U

is the peculiar velocity. Note that in Equation (4), the Boltzmann constant

k_{B} = 1

. We will take this value throughout the paper for the sake of simplicity. In addition, in Equations (1)–(3), n,

U

, and T denote the number density, the mean flow velocity, and the (granular) temperature, respectively, of the granular gas. They are defined as

\{n, n U, d n T\} = \int d v \{1, v, m V^{2}\} f (v) .

(5)

Since the only spatial gradient present in the USF problem is the shear rate, the pressure tensor

P_{g} = m_{g} \int d v V V f_{g} (v)

(6)

of the molecular gas, and the pressure tensor

P = m \int d v V V f (v)

(7)

of the granular gas are the relevant fluxes in the problem. They provide information on the transport of momentum across the system. Our main target is to determine

P_{g}

and

P

for arbitrary shear rates.

One of the main advantages of the USF at a microscopic level is that it becomes a spatially homogeneous state when the velocities of the particles are referred to a Lagrangian frame moving with the linear velocity

U_{i} = a_{i j} r_{j}

. In this new frame and in the steady state, the distribution functions of the molecular and granular gases adopt the form

f_{g} (r, v) = f_{g} (V), f (r, v) = f (V) .

(8)

In addition, as the state of the solvent is not perturbed by the solid particles, the temperature

T_{g}

in the USF state increases in time due to the viscous heating term

- a P_{x y} > 0

. Thus, as usual in nonequilibrium molecular dynamics simulations [29], an external nonconservative force (thermostat) is introduced in the molecular gas to achieve a stationary state. Among the different possibilities, for simplicity, a force proportional to the particle velocity (Gaussian thermostat) of the form

F_{g} = - m_{g} ξ V

is considered in this paper. The parameter

ξ

is chosen to be a function of the shear rate by the condition that

T_{g}

reaches a constant value in the long time limit. Analogously, the granular gas is also subjected to this kind of Gaussian thermostat (i.e.,

F = - m ξ V

), where

ξ

is the same quantity for the solvent and the solid particles.

Under the above conditions, in the low-density regime, the distribution function

f_{g} (V)

of the molecular gas obeys the nonlinear (closed) Boltzmann equation

- a V_{y} \frac{\partial f_{g}}{\partial V_{x}} - ξ \frac{\partial}{\partial V} \cdot V f_{g} = J_{g} [V | f_{g}, f_{g}],

(9)

while the distribution function

f (V)

of the granular gas obeys the kinetic equation

- a V_{y} \frac{\partial f}{\partial V_{x}} - ξ \frac{\partial}{\partial V} \cdot V f = J [V | f, f] + J_{BL} [V | f, f_{g}] .

(10)

Here,

J_{g} [f_{g}, f_{g}]

and

J [f, f]

are the nonlinear Boltzmann collision operators for the molecular and granular gases, respectively, and

J_{BL} [f, f_{g}]

is the linear Boltzmann–Lorentz collision operator [30,31]. The balance equations for

T_{g}

and T can easily obtained by multiplying both sides of Equations (9) and (10) by

m_{g} V^{2}

and

m V^{2}

, respectively, and integrating over velocity. The results are

- a P_{g, x y} = d ξ p_{g},

(11)

- a P_{x y} - \frac{d}{2} p ζ_{g} = d ξ p + \frac{d}{2} p ζ,

(12)

where

p_{g} = n_{g} T_{g}

and

p = n T

are the hydrostatic pressures of the molecular and granular gases, respectively, and the partial production rates

ζ

and

ζ_{g}

are defined, respectively, as

ζ = - \frac{m}{d n T} \int d v V^{2} J [v | f, f], ζ_{g} = - \frac{m}{d n T} \int d v V^{2} J_{BL} [v | f, f_{g}] .

(13)

The cooling rate

ζ

gives the rate of kinetic energy loss due to inelastic collisions between particles of the granular gas. It vanishes for elastic collisions

(α = 1)

. The term

ζ_{g}

gives the transfer of kinetic energy between the particles of the granular gas and the solvent. This term vanishes when the granular and molecular gases are at the same temperature (

T_{g} = T

). Equation (11) implies that, in the steady state, the viscous heating term (

- a P_{g, x y} > 0

) is exactly balanced by the heat extracted in the gas by the external thermostat. On the other hand, since

ζ_{g}

can be positive or negative, Equation (12) implies that, in the steady state, the term

- a P_{x y} - (d / 2) p ζ_{g}

is exactly compensated for the cooling terms arising from collisional dissipation (

ζ p

) and the thermostat term (

ξ p

The USF state is, in general, a non-Newtonian state characterized by shear-rate-dependent transport coefficients. In particular, one can define the non-Newtonian shear viscosity of the molecular gas as

η_{g} = - \frac{P_{g, x y}}{a} .

(14)

Analogously, the non-Newtonian shear viscosity of the granular gas is given by

η = - \frac{P_{x y}}{a} .

(15)

In addition, beyond the Navier–Stokes domain, normal stress differences are expected in the USF. This means that

P_{g, x x} \neq P_{g, y y} \neq P_{g, z z}

and

P_{x x} \neq P_{y y} \neq P_{z z}

It is quite evident that the evaluation of the rheological properties of the molecular and granular gases requires the knowledge of the pressure tensors

P_{g}

and

P

. The nonzero elements of these tensors can be obtained by multiplying by

m_{g} V V

and

m V V

both sides of Equations (9) and (10), respectively, and integrating over

V

. However, to achieve explicit forms for

P_{g}

and

P

, one has to compute the collisional moments

A_{g} = m_{g} \int d V V V J_{g} [f_{g}, f_{g}],

(16)

B = m \int d V V V J [f, f], C = m \int d V V V J_{BL} [f, f_{g}] .

(17)

In the case of IHS, the collisional moments

A_{g}

B

, and

C

cannot be exactly computed. As said in the Introduction section, a good estimate of them for IHS was made in Ref. [17] by means of Grad’s moment method [19]. This method is based on the expansion of the distributions

f_{g} (V)

and

f (V)

in a complete set of orthogonal polynomials, with the coefficients being the corresponding velocity moments of those distributions. The above expansion generates an infinite hierarchy of moment equations that must be truncated at a given order. This truncation allows one to arrive at a closed set of coupled equations for the velocity moments that can be recursively solved. Thus, since the results derived in Ref. [17] for the rheological properties of molecular and granular gases are approximated, it is convenient to revisit the problem and determine the exact expressions for the non-Newtonian transport properties of the granular suspension. To achieve such exact forms, two independent approaches will be considered in this paper: (i) the Boltzmann kinetic equation for IMM; and (ii) a BGK-type kinetic model for IHS. This task will be carried out in the next two sections.

3. Rheology from Inelastic Maxwell Models

In this section, we will consider an IMM, namely, a collisional model where the collision rate of the two colliding spheres are independent of their relative velocity. In this case, the Boltzmann collision operator

J_{g} [f_{g}, f_{g}]

(this is a simple version of the Boltzmann collision operator for Maxwell molecules) of the molecular gas can be written as [22]

J_{g} [v_{1} | f_{g}, f_{g}] = \frac{ν_{g}^{M}}{n_{g} S_{d}} \int d v_{2} \int d \hat{σ} [f_{g} (V_{1}^{″}) f_{g} (V_{2}^{″}) - f_{g} (V_{1}) f_{g} (V_{2})],

(18)

where

S_{d} = 2 π^{d / 2} / Γ (\frac{d}{2})

is the total solid angle in d dimensions and

ν_{g}^{M}

is an independent-velocity collision frequency. In Equation (18), the primes on the velocities denote the initial values

{V_{1}^{″}, V_{2}^{″}}

that lead to

{V_{1}, V_{2}}

following a binary collision:

V_{1}^{″} = V_{1} - (\hat{σ} \cdot g_{12}) \hat{σ}, V_{2}^{″} = V_{2} + (\hat{σ} \cdot g_{12}) \hat{σ} .

(19)

The effective collision frequency

ν_{g}^{M}

can be seen as a free parameter of the model to be chosen to attain agreement with the properties of interest of the original Boltzmann equation for IHS. For instance, to correctly capture the velocity dependence of the original IHS collision rate, we can assume that the IMM collision rate is proportional to

{\sqrt{T}}_{g}

In the context of IMM, the inelastic Boltzmann collision operator

J [f, f]

is [23,31]

J [V | f, f] = \frac{ν^{M}}{n S_{d}} \int d V_{2} \int d \hat{σ} [α^{- 1} f (V_{1}^{″}) f (V_{2}^{″}) - f (V_{1}) f (V_{2})],

(20)

while the Boltzmann–Lorentz collision operator

J_{BL} [f, f_{g}]

is defined as [22]

J_{BL} [V | f, f_{g}] = \frac{ν_{0}^{M}}{n S_{d}} \int d V_{2} \int d \hat{σ} [f (V_{1}^{″}) f_{g} (V_{2}^{″}) - f (V_{1}) f_{g} (V_{2})] .

(21)

The relationship between

(V_{1}^{″}, V_{2}^{″})

and

(V_{1}, V_{2})

in Equation (20) is

V_{1}^{″} = V_{1} - \frac{1 + α}{2 α} (\hat{σ} \cdot g_{12}) \hat{σ}, V_{2}^{″} = V_{2} + \frac{1 + α}{2 α} (\hat{σ} \cdot g_{12}) \hat{σ},

(22)

while in Equation (21), it is

V_{1}^{″} = V_{1} - 2 μ_{g} (\hat{σ} \cdot g_{12}) \hat{σ}, V_{2}^{″} = V_{2} + 2 μ (\hat{σ} \cdot g_{12}) \hat{σ}

(23)

where

μ_{g} = \frac{m_{g}}{m + m_{g}}, μ = \frac{m}{m + m_{g}} .

(24)

In addition, as in the case of the collision frequency

ν_{g}^{M}

, the collision frequencies

ν^{M}

and

ν_{0}^{M}

for granular–granular and granular–molecular collisions, respectively, can be chosen to optimize the agreement with the results derived from IHS. We will choose them later.

As mentioned in previous works on IMM [26], the main advantage of computing the collisional moments of the Boltzmann operator for Maxwell models (both elastic and inelastic models) is that they can be exactly provided in terms of the velocity moments of the distribution functions without the explicit knowledge of the latter. This property has been exploited to compute the second-, third-, and fourth-degree collisional moments of IMM for monocomponent and multicomponent granular gases [31]. The exact knowledge of the second-degree collisional moments allow us to obtain exact expressions for the rheological properties of the molecular and granular gases. Let us separately evaluate the rheology of both gases.

3.1. Rheological Properties of the Molecular Gas

The pressure tensor of the molecular gas is defined by Equation (6). To obtain the nonzero elements of this tensor, one multiplies both sides of the Boltzmann equation (9) by

m_{g} V_{i} V_{j}

and integrates over velocity. The result is

a_{i k} P_{g, k j} + a_{j k} P_{g, k i} + 2 ξ P_{g, i j} = - \frac{2 ν_{g}^{M}}{d + 2} (P_{g, i j} - p_{g} δ_{i j}),

(25)

where

p_{g} = n_{g} T_{g}

is the hydrostatic pressure of the molecular gas. Upon obtaining Equation (25), use is made of the result [31]:

A_{g, i j} = m_{g} \int d V V_{i} V_{j} J_{g} [f_{g}, f_{g}] = - \frac{2 ν_{g}^{M}}{d + 2} (P_{g, i j} - p_{g} δ_{i j}) .

(26)

The (reduced) elements of the pressure tensor

P_{g}^{*} = P_{g} / (n_{g} T_{g})

can be easily obtained from Equation (25). They are given by

P_{g, x x}^{*} = \frac{1}{1 + 2 \tilde{ξ}} [1 + \frac{2 {\tilde{a}}^{2}}{{(1 + 2 \tilde{ξ})}^{2}}], P_{g, y y}^{*} = P_{g, z z}^{*} = \frac{1}{1 + 2 \tilde{ξ}}, P_{g, x y}^{*} = - \frac{\tilde{a}}{{(1 + 2 \tilde{ξ})}^{2}} .

(27)

Here, we introduce the quantities

\tilde{a} = \frac{a}{{\tilde{ν}}_{g}^{M}}, \tilde{ξ} = \frac{ξ}{{\tilde{ν}}_{g}^{M}}, {\tilde{ν}}_{g}^{M} = \frac{2}{d + 2} ν_{g}^{M} .

(28)

The constraint

P_{g, x x}^{*} + (d - 1) P_{g, y y}^{*} = d

leads to a cubic equation relating

\tilde{ξ}

and

\tilde{a}

{\tilde{a}}^{2} = d \tilde{ξ} {(1 + 2 \tilde{ξ})}^{2} .

(29)

The real root of Equation (29) gives the shear-rate dependence of

\tilde{ξ} (\tilde{a})

. It is given by

\tilde{ξ} (\tilde{a}) = \frac{2}{3} {sinh}^{2} [\frac{1}{6} {cosh}^{- 1} (1 + \frac{27}{d} {\tilde{a}}^{2})] .

(30)

Comparison between the results derived here for

P_{g, i j}^{*}

with those recently [17] obtained for IHS by means of Grad’s moment method [19] shows that both results are identical if the effective collision frequency

ν_{g}^{M}

is given by

ν_{g}^{M} = \frac{4 π^{(d - 1) / 2}}{Γ (\frac{d}{2})} n_{g} σ_{g}^{d - 1} \sqrt{\frac{T_{g}}{m_{g}}} .

(31)

Henceforth, we will take the choice (31) for

ν_{g}^{M}

From Equations (27), one can identify the (dimensionless) non-Newtonian shear viscosity

η_{g}^{*} = ν_{g}^{M} η_{g} / p_{g} = - P_{g, x y}^{*} / \tilde{a}

and the (dimensionless) normal stress difference

Ψ_{g}^{*} = P_{g, x x}^{*} - P_{g, y y}^{*}

η_{g}^{*} = \frac{1}{{(1 + 2 \tilde{ξ})}^{2}}, Ψ_{g}^{*} = \frac{2 {\tilde{a}}^{2}}{{(1 + 2 \tilde{ξ})}^{3}} .

(32)

Note that the results derived here for Maxwell molecules yield

P_{g, y y} = P_{g, z z}

. This result contrasts with the one obtained for hard spheres by numerically solving the Boltzmann equation by means of the DSMC method [28], where it has been shown that

P_{g, y y} \neq P_{g, z z}

. However, the difference (

P_{g, y y} - P_{g, z z}

) found in the simulations is, in general, quite small [9].

It is also important to note that in the case of Maxwell molecules there is an exact equivalence between the description with and without the drag force

- m ξ V

. Nevertheless, for non-Maxwell molecules, this type of force does not play a neutral role in the transport properties of the system [32].

The shear-rate dependence of

η_{g}^{*}

and

Ψ_{g}^{*}

is plotted in Figure 1 for a three-dimensional system (

d = 3

). As expected, the nonlinear viscosity

η_{g}^{*}

decreases with increasing (reduced) shear rate

\tilde{a}

(shear thinning effect). The opposite effect is observed for the normal stress difference function

Ψ_{g}^{*}

since it increases with the shear rate. Figure 1 also highlights the excellent agreement found between the theoretical results for the Maxwell molecules with those obtained by numerically solving the Boltzmann equation for hard spheres from the DSMC method [28].

3.2. Rheological Properties of the Granular Gas

As in the case of the molecular gas, the rheology of the granular gas can be also determined by multiplying both sides of Equation (10) by

m V_{i} V_{j}

and integrating over

V

. After some algebra, one achieves the result

\begin{matrix} a_{i k} P_{k j} + a_{j k} P_{k i} + 2 ξ P_{i j} = - ν^{M} [ν_{η}^{*} P_{i j} + p (ζ^{*} - ν_{η}^{*}) δ_{i j}] - \frac{4}{d} ν_{0}^{M} μ_{g} \\ \times [(1 - \frac{2}{d + 2} μ_{g}) P_{i j} - \frac{2}{d + 2} μ \frac{n}{n_{g}} P_{g, i j} - \frac{d}{d + 2} μ_{g} n T_{g} (χ + \frac{m}{m_{g}}) δ_{i j}], \end{matrix}

(33)

where use is made of the results [31]:

B_{i j} = m \int d V V_{i} V_{j} J [f, f] = - ν^{M} [ν_{η}^{*} P_{i j} + p (ζ^{*} - ν_{η}^{*}) δ_{i j}],

(34)

\begin{matrix} C_{i j} & = & m \int d V V_{i} V_{j} J_{BL} [f, f_{g}] = - \frac{4}{d} ν_{0}^{M} μ_{g} [(1 - \frac{2}{d + 2} μ_{g}) P_{i j} - \frac{2}{d + 2} μ \frac{n}{n_{g}} P_{g, i j} \\ - \frac{d}{d + 2} μ_{g} n T_{g} (χ + \frac{m}{m_{g}}) δ_{i j}] . \end{matrix}

(35)

In Equation (33),

ζ^{*} = \frac{ζ}{ν^{M}} = \frac{1 - α^{2}}{2 d}

(36)

is the (reduced) cooling rate for the granular gas,

χ = T / T_{g}

is the temperature ratio, and

ν_{η}^{*} = \frac{(d + 1 - α) (1 + α)}{d (d + 2)} .

(37)

The partial cooling rate

ζ_{g}

can be exactly obtained from Equation (34) as

ζ_{g} = \frac{4 ν_{0}^{M}}{d} μ_{g} [1 - μ_{g} (1 + θ)],

(38)

where

θ = \frac{m T_{g}}{m_{g} T}

(39)

is the ratio of the mean square velocities of granular and molecular gas particles. The forms (36) and (38) can be employed to fix the values of the free parameters

ν^{M}

and

ν_{0}^{M}

. They are chosen under the criterion that

ζ

and

ζ_{g}

of IMM are the same as that of IHS of diameters

σ

and

σ_{0}

. In this latter case, the above cooling rates are estimated by using Grad’s approximation [17]. In this approximation,

ζ^{IHS} = \frac{2 π^{(d - 1) / 2}}{d Γ (\frac{d}{2})} n σ^{d - 1} \sqrt{\frac{T}{m}} (1 - α^{2}),

(40)

ζ_{g}^{IHS} = \frac{8 π^{(d - 1) / 2}}{d Γ (\frac{d}{2})} n_{g} {\bar{σ}}^{d - 1} μ_{g} {(\frac{1 + θ}{θ})}^{1 / 2} \sqrt{\frac{2 T_{g}}{m_{g}}} [1 - μ_{g} (1 + θ)],

(41)

where

\bar{σ} = (σ + σ_{g}) / 2

. Equations (38), (40), and (41) yield the identities

ν^{M} = \frac{4 π^{(d - 1) / 2}}{Γ (\frac{d}{2})} n σ^{d - 1} \sqrt{\frac{T}{m}}, ν_{0}^{M} = \frac{2 π^{(d - 1) / 2}}{Γ (\frac{d}{2})} n_{g} {\bar{σ}}^{d - 1} {(\frac{1 + θ}{θ})}^{1 / 2} \sqrt{\frac{2 T_{g}}{m_{g}}} .

(42)

To compare with the rheological properties of IHS [17], it is convenient at this level of the description to identify the friction (or drift) coefficient

γ

appearing in the Brownian limiting case (

m / m_{g} \to \infty

) when the molecular gas is at equilibrium. In fact, this limiting case is the situation considered when one employs a coarse-grained approach [1,2,3,4] to assess the impact of the interstitial gas on the dynamics properties of grains. In this limiting case, the expression (35) of the collisional moment

C_{i j}

reduces to

C_{i j}^{Br} = - \frac{8 π^{(d - 1) / 2}}{d Γ (\frac{d}{2})} n T_{g} n_{g} {\bar{σ}}^{d - 1} {(\frac{m_{g}}{m})}^{1 / 2} \sqrt{\frac{2 T_{g}}{m}} (P_{k ℓ}^{*} - δ_{k ℓ}),

(43)

where

P_{i j}^{*} = P_{i j} / (n T_{g})

, and we take into account that in the Brownian limit

μ_{g} \to m_{g} / m

and

(1 + θ) / θ \to 1

in the expression (42) of

ν_{0}^{M}

. The form of

C_{i j}

derived in Ref. [26] by replacing the Boltzmann–Lorentz collisional operator (21) by the Fokker–Planck operator

γ \frac{\partial}{\partial V} \cdot V f + γ \frac{T_{g}}{m} \frac{\partial^{2} f}{\partial V^{2}}

(44)

C_{i j} = - 2 γ n T_{g} (P_{k ℓ}^{*} - δ_{k ℓ}) .

(45)

Comparison between Equations (43) and (45) allows us to identify

γ

for IMM as

γ = \frac{4 π^{(d - 1) / 2}}{d Γ (\frac{d}{2})} n_{g} {\bar{σ}}^{d - 1} {(\frac{m_{g}}{m})}^{1 / 2} {(\frac{2 T_{g}}{m})}^{1 / 2} .

(46)

The expression (46) for the friction coefficient

γ

for IMM is the same as the one obtained for IHS [15,17].

We are now in a position to determine the nonzero elements of the (reduced) pressure tensor

P_{i j}^{*}

. From Equations (33) and (46), one obtains the equation

\begin{matrix} a_{i k}^{*} P_{k j}^{*} + a_{j k}^{*} P_{k i}^{*} + 2 ξ^{*} P_{i j}^{*} & = & - ν^{* M} [ν_{η}^{*} P_{i j}^{*} + χ (ζ^{*} - ν_{η}^{*}) δ_{i j}] - 2 μ {(\frac{1 + θ}{θ})}^{1 / 2} \\ \times (X δ_{i j} + Y P_{i j}^{*} + Z P_{g, i j}^{*}), \end{matrix}

(47)

where

a_{i j}^{*} = a_{i j} / γ

ξ^{*} = ξ / γ

ν^{* M} = \frac{ν^{M}}{γ} = \frac{2^{d + 1} d}{\sqrt{π}} ϕ \sqrt{χ T_{g}^{*}},

(48)

and

X = - \frac{d}{d + 2} μ (\frac{1 + θ}{θ}), Y = 1 - \frac{2 μ_{g}}{d + 2}, Z = - \frac{2 μ}{d + 2} .

(49)

Here,

ϕ = \frac{π^{d / 2}}{2^{d - 1} d Γ (\frac{d}{2})} n σ^{d}

(50)

is the solid volume fraction of the granular gas,

T_{g}^{*} = T_{g} / m σ^{2} γ^{2}

, and upon deriving Equation (48), use is made of the identity

\frac{n σ^{d - 1}}{n_{g} {\bar{σ}}^{d - 1}} = \frac{2^{d + \frac{3}{2}}}{\sqrt{π}} {(\frac{m_{g}}{m})}^{1 / 2} ϕ \sqrt{T_{g}^{*}} .

(51)

As occurs for the rheology of the molecular gas, Equation (47) shows that the diagonal elements of the pressure tensor

P_{i j}^{*}

orthogonal to the shear plane

x y

are equal to

P_{y y}^{*}

(i.e.,

P_{y y}^{*} = P_{z z}^{*} = \dots = P_{d d}^{*}

). This implies that the

x x

element is given by

P_{x x}^{*} = d χ - (d - 1) P_{y y}^{*}

. The

y y

and

x y

elements of the (reduced) pressure tensor can be written as

P_{y y}^{*} = \frac{Ω_{y y}}{ν_{y y}}, P_{x y}^{*} = \frac{Ω_{x y} - a^{*} P_{y y}^{*}}{ν_{y y}},

(52)

where

ν_{y y} = 2 ξ^{*} + ν^{* M} ν_{η}^{*} + 2 μ {(\frac{1 + θ}{θ})}^{1 / 2} Y,

(53)

Ω_{y y} = - ν^{* M} χ (ζ^{*} - ν_{η}^{*}) - 2 μ {(\frac{1 + θ}{θ})}^{1 / 2} (X + Z P_{g, y y}^{*}),

(54)

Ω_{x y} = - 2 μ {(\frac{1 + θ}{θ})}^{1 / 2} Z P_{g, x y}^{*} .

(55)

Note that the elements of the pressure tensor

P_{g, y y}^{*}

and

P_{g, x y}^{*}

of the molecular gas must be expressed in terms of the (reduced) shear rate

a^{*}

and the (reduced) thermostat parameter

ξ^{*}

. For this, one has to take into account the relationships between

\tilde{a}

and

\tilde{ξ}

with

a^{*}

and

ξ^{*}

, respectively. They are given by

\tilde{a} = (γ / {\tilde{ν}}_{g}^{M}) a^{*}

and

\tilde{ξ} = (γ / {\tilde{ν}}_{g}^{M}) ξ^{*}

, where

\frac{γ}{{\tilde{ν}}_{g}^{M}} = \frac{d + 2}{\sqrt{2} d} {(\frac{\bar{σ}}{σ_{g}})}^{d - 1} \frac{m_{g}}{m} .

(56)

The equation defining the temperature ratio

χ

can be easily derived from Equation (47) as

\frac{2}{d} a^{*} P_{x y}^{*} + 2 ξ^{*} χ = - ν^{* M} χ ζ^{*} + 2 μ^{2} {(\frac{1 + θ}{θ})}^{1 / 2} (1 - χ) .

(57)

From Equations (52) and (57), one finally obtains

a^{*}

in terms of the parameter space of the system:

a^{*} = \sqrt{\frac{d}{2} \frac{2 μ^{2} {(1 + θ^{- 1})}^{1 / 2} (1 - χ) - (ν^{* M} ζ^{*} + 2 ξ^{*}) χ}{\frac{Ω_{x y} / a^{*}}{ν_{y y}} - \frac{Ω_{y y}}{ν_{y y}^{2}}}} .

(58)

As happens in the case of IHS [17], the temperature ratio

χ

cannot be expressed in Equation (58) as an explicit function of the (reduced) shear rate and the remaining parameters of the system. On the other hand, for given values of the parameter space

Ξ \equiv (χ, α, σ / σ_{g}, m / m_{g}, ϕ, T_{g}^{*})

, the temperature ratio can be implicitly determined from the physical solution to Equation (58).

3.3. Brownian Limit

Before illustrating the shear-rate dependence of the rheological properties of the molecular gas for arbitrary values of the mass ratio

m / m_{g}

, it is convenient to check the consistency of the present results with those derived in Ref. [26] for IMM by using the Fokker–Planck operator (44). This consistency is expected to apply in the Brownian limit

m / m_{g} \to \infty

. In this limiting case, at a given value of the (reduced) shear rate

a^{*}

θ \to \infty

γ / ν_{g}^{M} \propto m_{g} / m \to 0

\tilde{a} \propto m_{g} / m \to 0

\tilde{ξ} \propto {\tilde{a}}^{2} \propto {(m_{g} / m)}^{2} \to 0

, and

ξ^{*} = \tilde{ξ} ({\tilde{ν}}_{g}^{M} / γ) \propto m_{g} / m \to 0

. Consequently,

P_{g, i j}^{*} = δ_{i j}

and the expressions of

P_{y y}^{*}

P_{x y}^{*}

, and

a^{*}

are

P_{y y}^{*} = \frac{2 - ν^{* M} χ (ζ^{*} - ν_{η}^{*})}{2 + ν^{* M} ν_{η}^{*}}, P_{x y}^{*} = - \frac{2 - ν^{* M} χ (ζ^{*} - ν_{η}^{*})}{{(2 + ν^{* M} ν_{η}^{*})}^{2}} a^{*},

(59)

a^{*} = \sqrt{\frac{d}{2} \frac{ν^{* M} ζ^{*} + 2 (1 - χ^{- 1})}{ν^{* M} (ν_{η}^{*} - ζ^{*}) + 2 χ^{- 1}}} (2 + ν^{* M} ν_{η}^{*}) .

(60)

Equations (59) and (60) are consistent with Equations (32), (33), and (35) of Ref. [26]. It is important to note that, to assess consistency with the Fokker–Planck results, the size ratio has been kept constant or proportional to the mass ratio so that

ξ^{*} \to 0

4. Rheology from a BGK-Type Kinetic Model of the Boltzmann Equation

We consider in this section the results derived for the USF from a BGK-type kinetic model of the Boltzmann equation [24,33]. In the problem for the granular suspension considered here, one has to replace the true Boltzmann operators

J_{g} [f_{g}, f_{g}]

J [f, f]

, and

J_{BL} [f, f_{g}]

by simpler relaxation terms that retain the relevant physical properties of those operators but are more tractable than the true kinetic equations. As in the case of IMM, let us separately determine the rheological properties of the molecular and granular gases by starting from these kinetic models.

4.1. Rheological Properties of the Molecular Gas

In the case of the molecular gas, the Boltzmann collision operator

J_{g} [f_{g}, f_{g}]

is replaced by the conventional BGK kinetic model [22,34]:

J_{g} [f_{g}, f_{g}] \to - ν_{g} (f_{g} - f_{g}^{M}),

(61)

where

ν_{g}

is an effective velocity-independent collision frequency and

f_{g}^{M}

is the Maxwellian distribution

f_{g}^{M} (V) = n_{g} {(\frac{m_{g}}{2 π T_{g}})}^{d / 2} exp (- \frac{m_{g} V^{2}}{2 T_{g}}) .

(62)

Thus, according to Equation (9), the velocity distribution function

f_{g} (V)

obeys the BGK kinetic equation

- a V_{y} \frac{\partial f_{g}}{\partial V_{x}} - ξ \frac{\partial}{\partial V} \cdot V f_{g} = - ν_{g} (f_{g} - f_{g}^{M}) .

(63)

The nonzero elements of the pressure tensor

P_{g, i j}

can be easily obtained from Equation (63) by multiplying both sides of this equation by

m_{g} V_{i} V_{j}

and integrating over

V

. The BGK expressions of the (reduced) elements of the pressure tensor

P_{g, i j}^{*}

are given by Equations (27) with the replacement

{\tilde{ν}}_{g}^{M} \to ν_{g}

. As a consequence, the results derived from the BGK equation for the rheological properties agree with those obtained from the Boltzmann equation for hard spheres when

ν_{g} = {\tilde{ν}}_{g}^{M} = \frac{8 π^{(d - 1) / 2}}{(d + 2) Γ (\frac{d}{2})} n_{g} σ_{g}^{d - 1} \sqrt{\frac{T_{g}}{m_{g}}} .

(64)

4.2. Rheological Properties of the Granular Gas

In the case of the granular gas, we consider the kinetic model proposed by Vega Reyes et al. [24] for granular mixtures. This kinetic model is based on the equivalence between a system of elastic spheres subject to a drag force proportional to the peculiar velocity

V

with a gas of IHS [35]. Thus, the true inelastic Boltzmann collision operator is replaced by a relaxation term plus a drag force term. This (approximate) mapping between elastic hard spheres plus drag force with IHS allows us to extend known kinetic models of molecular gases to inelastic gases. Here, we consider the well-known Gross and Krook (GK) model for molecular (elastic) mixtures [36] for the corresponding relaxation terms. With this approach, the collision operators

J [f, f]

and

J_{BL} [f, f_{g}]

are defined, respectively, as

J [f, f] \to - ν^{'} (α) (f - f^{M}) + \frac{ϵ}{2} \frac{\partial}{\partial V} \cdot V f,

(65)

J_{BL} [f, f_{g}] \to - ν (f - {\tilde{f}}_{g}) .

(66)

In Equations (65) and (66),

f^{M} (V) = n {(\frac{m}{2 π T})}^{d / 2} exp (- \frac{m V^{2}}{2 T})

(67)

is the Maxwellian distribution of the granular gas,

ϵ = \frac{2 π^{(d - 1) / 2}}{d Γ (\frac{d}{2})} n σ^{d - 1} \sqrt{\frac{T}{m}} (1 - α^{2}),

(68)

and the reference distribution function

{\tilde{f}}_{g} (V)

is [36]

{\tilde{f}}_{g} (V) = n {(\frac{m}{2 π \tilde{T}})}^{d / 2} exp (- \frac{m V^{2}}{2 \tilde{T}}) .

(69)

Note that in the definition (68) of

ϵ

, the cooling rate

ζ

was evaluated by using the Maxwellian approximation for the distribution f. In this case,

ζ

is given by Equation (40).

In Equations (65) and (66), the quantities

ν^{'}

ν

, and

\tilde{T}

are chosen to optimize the agreement with some properties of interest of the Boltzmann equation for IHS. The usual method of obtaining the above parameters is to ensure that the kinetic model reproduces the collisional transfer equations of momentum and energy for elastic collisions (

α = 1

). However, since

U = U_{g}

in the USF, we only have one constraint (the one associated with the transfer of energy) instead of two, so that

\tilde{T}

and

ν

admit more than one form. Here, we propose a choice (see Appendix A for more technical details) that leads to an excellent agreement with the results obtained for IHS from Grad’s moment method [17]. More specifically, we take the following values of

\tilde{T}

and

ν

\tilde{T} = T_{g}, ν = \frac{8 π^{(d - 1) / 2}}{d Γ (\frac{d}{2})} n_{g} {\bar{σ}}^{d - 1} \frac{m m_{g}}{{(m + m_{g})}^{2}} {(\frac{2 T_{g}}{m_{g}} + \frac{2 T}{m})}^{1 / 2} .

(70)

It still remains to completely define the model to chose the effective collision frequency

ν^{'}

. It is defined here to reproduce the collisional moment

\int d v m V_{i} V_{j} J [f, f]

(71)

of the Boltzmann equation for IHS when one takes Grad’s trial distribution for f [17]. This leads to the expression

ν^{'} = \frac{2 π^{(d - 1) / 2}}{d (d + 2) Γ (\frac{d}{2})} n σ^{d - 1} \sqrt{\frac{T}{m}} (1 + α) [d + 1 + (d - 1) α] .

(72)

Therefore, the BGK kinetic equation for the sheared granular gas is given by

- a V_{y} \frac{\partial f}{\partial V_{x}} - ξ \frac{\partial}{\partial V} \cdot V f = - ν^{'} (f - f^{M}) + \frac{ϵ}{2} \frac{\partial}{\partial V} \cdot V f - ν (f - {\tilde{f}}_{g}),

(73)

where

ν

and

ν^{'}

are defined by Equations (70) and (72), respectively, and the Maxwellian distribution

{\tilde{f}}_{g}

is given by Equation (69) with

\tilde{T} = T_{g}

The possibility of determining all the velocity moments of the distribution function is likely one of the main advantages of employing a kinetic model instead of the true Boltzmann equation. In the USF problem, it is convenient to define the general velocity moments

M_{k_{1}, k_{2}, k_{3}} = \int d V V_{x}^{k_{1}} V_{y}^{k_{2}} V_{z}^{k_{3}} f (V) .

(74)

Although, here, we are essentially interested in the three-dimensional case, we will compute the velocity moments for

d = 3

and

d = 2

. Note that for hard disks (

d = 2

k_{3} = 0

in Equation (74), since the z-axis is meaningless. The hierarchy of moment equations can be obtained by multiplying Equation (73) by

V_{x}^{k_{1}} V_{y}^{k_{2}} V_{z}^{k_{3}}

and integrating over

V

. The result is

a k_{1} M_{k_{1} - 1, k_{2} + 1, k_{3}} + (ν^{'} + ν + k λ) M_{k_{1}, k_{2}, k_{3}} = N_{k_{1}, k_{2}, k_{3}},

(75)

where

λ = ξ + ϵ / 2

k = k_{1} + k_{2} + k_{3}

, and

N_{k_{1}, k_{2}, k_{3}} = n {(\frac{2 T_{g}}{m})}^{k / 2} (ν + χ^{k / 2} ν^{'}) M_{k_{1}, k_{2}, k_{3}}^{L} .

(76)

In Equation (75), for hard spheres (

d = 3

)

M_{k_{1}, k_{2}, k_{3}}^{L} = \int d c c_{x}^{k_{1}} c_{y}^{k_{2}} c_{z}^{k_{3}} e^{- c^{2}} = π^{- 3 / 2} Γ (\frac{k_{1} + 1}{2}) Γ (\frac{k_{2} + 1}{2}) Γ (\frac{k_{3} + 1}{2}),

(77)

k_{1}

k_{2}

, and

k_{3}

are even, being zero otherwise. For hard disks (

d = 2

M_{k_{1}, k_{2}, 0}^{L} = π^{- 1} Γ (\frac{k_{1} + 1}{2}) Γ (\frac{k_{2} + 1}{2})

(78)

k_{1}

and

k_{2}

are even, being zero otherwise.

The solution to Equation (75) can be written as (see Appendix B of Ref. [26] for some details)

M_{k_{1}, k_{2}, k_{3}} = \sum_{q = 0}^{k_{1}} \frac{k_{1}!}{(k_{1} - q)!} \frac{{(- a)}^{q}}{{(ν^{'} + ν + k λ)}^{1 + q}} N_{k_{1} - q, k_{2} + q, k_{3}} .

(79)

The nonzero elements of the pressure tensor

P_{i j}

can be easily obtained from Equation (79). In dimensionless form, the BGK expressions of the elements of

P_{i j}^{*} = P_{i j} / n T_{g}

are

P_{y y}^{*} = P_{z z}^{*} = \frac{ν + ν^{'} χ}{ν^{'} + ν + 2 λ}, P_{x y}^{*} = - \frac{ν + ν^{'} χ}{{(ν^{'} + ν + 2 λ)}^{2}} a,

(80)

P_{x x}^{*} = \frac{ν + ν^{'} χ}{ν^{'} + ν + 2 λ} [1 + \frac{2 a^{2}}{{(ν^{'} + ν + 2 λ)}^{2}}] .

(81)

The (steady) temperature ratio

χ = T / T_{g}

can be obtained from the constraint

P_{x x}^{*} + (d - 1) P_{y y}^{*} = d χ

. This yields the implicit equation

a^{*} = (ν^{' *} + ν^{*} + 2 λ^{*}) \sqrt{\frac{d}{2} χ \frac{ν^{*} (1 - χ^{- 1}) + 2 λ^{*}}{ν^{*} + ν^{' *} χ}},

(82)

where

a^{*} = a / γ

ν^{' *} = \frac{ν^{'}}{γ} = \frac{2^{d}}{d + 2} ϕ \sqrt{\frac{χ T_{g}^{*}}{π}} (1 + α) [d + 1 + (d - 1) α],

(83)

ν^{*} = \frac{ν}{γ} = 2 μ^{2} {(\frac{1 + θ}{θ})}^{1 / 2}, λ^{*} = \frac{λ}{γ} = ξ^{*} + \frac{ϵ^{*}}{2},

(84)

ξ^{*} = \frac{\sqrt{2} d}{d + 2} {(\frac{σ_{g}}{\bar{σ}})}^{d - 1} \frac{m}{m_{g}} \tilde{ξ}, ϵ^{*} = 2^{d} ϕ \sqrt{\frac{χ T_{g}^{*}}{π}} (1 - α^{2}) .

(85)

Here,

\tilde{ξ}

is given by Equation (30) and we recall that

T_{g}^{*} = T_{g} / m σ^{2} γ^{2}

4.3. Brownian Limit

As in the case of IMM, it is quite interesting to consider the Brownian limiting case

m / m_{g} \to \infty

. In this case,

P_{g, i j}^{*} = δ_{i j}

μ \to 1

θ \to \infty

, and

ξ^{*} \to 0

. Thus, following similar steps to those made for IMM, one obtains the expressions

P_{y y}^{*} = \frac{2 + ν^{' *} χ}{2 + ν^{' *} + ϵ^{*}}, P_{x y}^{*} = - \frac{2 + ν^{' *} χ}{{(2 + ν^{' *} + ϵ^{*})}^{2}} a^{*},

(86)

a^{*} = (2 + ν^{' *} + ϵ^{*}) \sqrt{\frac{d}{2} χ \frac{2 (1 - χ^{- 1}) + ϵ^{*}}{2 + ν^{' *} χ}} .

(87)

4.4. Velocity Distribution of the Granular Gas

Apart from obtaining all the velocity moments, the use of kinetic models allows us, in some cases, to explicitly determine the velocity distribution function

f (V)

. The BGK-type Equation (73) reads

(1 - d \hat{λ} - \hat{a} V_{y} \frac{\partial}{\partial V_{x}} - \hat{λ} V \cdot \frac{\partial}{\partial V}) f (V) = f_{R} (V),

(88)

where

\hat{a} = a / (ν^{'} + ν)

\hat{λ} = λ / (ν^{'} + ν)

, and

f_{R} (V) = \frac{ν^{'}}{ν^{'} + ν} f^{M} (V) + \frac{ν}{ν^{'} + ν} {\tilde{f}}_{g} (V) .

(89)

The hydrodynamic solution to Equation (88) can be formally written as

\begin{matrix} f (V) & = & {(1 - d \hat{λ} - \hat{a} V_{y} \frac{\partial}{\partial V_{x}} - \hat{λ} V \cdot \frac{\partial}{\partial V})}^{- 1} f_{R} (V) \\ = & \int_{0}^{\infty} d s e^{- (1 - d \hat{λ}) s} e^{\hat{a} s V_{y} \frac{\partial}{\partial V_{x}}} e^{\hat{λ} s V \cdot \frac{\partial}{\partial V}} f_{R} (V) . \end{matrix}

(90)

The action of the velocity operators in Equation (90) on an arbitrary function

g (V)

is [26]

e^{\hat{a} s V_{y} \frac{\partial}{\partial V_{x}}} g (V_{x}, V_{y}, V_{z}) = g (V_{x} + \hat{a} s V_{y}, V_{y}, V_{z}),

(91)

e^{\hat{λ} s V \cdot \frac{\partial}{\partial V}} g (V_{x}, V_{y}, V_{z}) = g (e^{\hat{λ} s} V_{x}, e^{\hat{λ} s} V_{y}, e^{\hat{λ} s} V_{z}) .

(92)

The explicit form of the one-particle velocity distribution function can be explicitly obtained when one takes into account in Equation (90) the action of the velocity operators given by Equations (91) and (92). The result can be written as

f (V) = n {(\frac{m}{2 T})}^{d / 2} φ (c), c = {(\frac{m}{2 T})}^{1 / 2} V,

(93)

where

\begin{matrix} φ (c) & = & π^{- d / 2} \int_{0}^{\infty} d s e^{- (1 - d \hat{λ}) s} {\frac{ν^{'}}{ν^{'} + ν} exp [- e^{2 \hat{λ} s} {(c + s \hat{a} \cdot c)}^{2}] \\ + \frac{ν}{ν^{'} + ν} χ^{d / 2} exp [- χ e^{2 \hat{λ} s} {(c + s \hat{a} \cdot c)}^{2}]} . \end{matrix}

(94)

Here, we introduce the tensor

{\hat{a}}_{i j} = \hat{a} δ_{i x} δ_{j y}

In order to illustrate the shear-rate dependence of the distribution function, it is convenient to consider the marginal distribution for

d = 3

\begin{matrix} φ_{x} (c_{x}) & = & \int_{- \infty}^{\infty} d c_{y} \int_{- \infty}^{\infty} d c_{z} φ (c) \\ = & \frac{1}{\sqrt{π}} \int_{0}^{\infty} d s \frac{e^{- (1 - \hat{λ}) s}}{\sqrt{1 + {\hat{a}}^{2} s^{2}}} {\frac{ν^{'}}{ν^{'} + ν} exp (- e^{2 \hat{λ} s} \frac{c_{x}^{2}}{1 + {\hat{a}}^{2} s^{2}}) \\ + \frac{ν}{ν^{'} + ν} χ^{1 / 2} exp (- χ e^{2 \hat{λ} s} \frac{c_{x}^{2}}{1 + {\hat{a}}^{2} s^{2}})} . \end{matrix}

(95)

In the Brownian limit,

ξ^{*} \to 0

μ \to 1

θ \to \infty

, and so

ν^{*} \to 2

λ^{*} \to ϵ^{*} / 2

, and

\hat{λ} \to \frac{ϵ^{*} / 2}{2 + ν^{' *}}, \hat{a} \to \frac{a^{*}}{2 + ν^{' *}} .

(96)

Thus, when

m / m_{g} \to \infty

, Equation (95) becomes

\begin{matrix} φ_{x} (c_{x}) & = & \frac{1}{\sqrt{π}} \int_{0}^{\infty} d s \frac{e^{- (1 - \hat{λ}) s}}{\sqrt{1 + {\hat{a}}^{2} s^{2}}} {\frac{ν^{' *}}{2 + ν^{' *}} exp (- e^{2 \hat{λ} s} \frac{c_{x}^{2}}{1 + {\hat{a}}^{2} s^{2}}) \\ + \frac{2}{2 + ν^{' *}} χ^{1 / 2} exp (- χ e^{2 \hat{λ} s} \frac{c_{x}^{2}}{1 + {\hat{a}}^{2} s^{2}})}, \end{matrix}

(97)

where

\hat{λ}

and

\hat{a}

are given by Equations (96).

5. Comparison between IMM and BGK Results

In Section 3 and Section 4, we made use of the Boltzmann equation for IMM and the BGK-type kinetic model to investigate the shear-rate dependence of rheological properties in a sheared granular suspension. These properties are expressed in terms of the coefficient of restitution

α

, the reduced background temperature

T_{g}^{*}

, and the diameter

σ / σ_{g}

and mass

m / m_{g}

ratios. Additionally, there exists a residual dependence on density through the volume fraction

ϕ

. To avoid that, one could, for instance, reduce the shear rate using the effective collision frequencies

ν^{M} (T)

ν (T)

. However, for consistency with simulations and considering the background temperature

T_{g}

as a known quantity, we opted to employ

γ (T_{g})

as the reference frequency.

Given that in this section the second-degree moments of the distribution function are compared with molecular dynamics (MD) simulations for IHS in the Brownian limiting case [18], we set fixed values of

T_{g}^{*} = 1

and

ϕ = 0.0052

for subsequent analysis. The selection of

T_{g}^{*}

as a free parameter imposes a constraint between the diameter

σ / σ_{g}

and mass

m / m_{g}

ratios [17]:

\frac{σ}{σ_{g}} = {[{(\frac{\sqrt{π}}{4 \sqrt{2}} \frac{n}{n_{g}} \sqrt{\frac{m}{m_{g}}} \frac{1}{ϕ \sqrt{T_{g}^{*}}})}^{1 / (d - 1)} - 1]}^{- 1} .

(98)

This relation ensures convergence of results to those obtained via the Fokker–Planck equation as

m / m_{g} \to \infty

, since

ξ^{*} \to 0

. Furthermore, since we want to recover the results obtained in Ref. [17] derived from Grad’s method, we take

n / n_{g} = 10^{- 3}

(rarefied granular gas).

The second-degree moments expressed through the reduced temperature

χ

, non-Newtonian shear viscosity

η^{*}

, and the normal stress difference

Ψ^{*}

are plotted in Figure 2 for

α = 0.9

and 1. Here,

Ψ^{*} = P_{x x}^{*} - P_{y y}^{*} = d χ - d P_{y y}^{*}

. Equations (52) and (80) provide analytical expressions for IMM and BGK-type kinetic model, respectively, from which rheological properties are illustrated. Notably, there is nearly perfect agreement between Grad’s solution for IHS, as obtained in Ref. [17], and both IMM and BGK-type results for any mass ratio, highlighting the ability of relatively simple models to capture essential properties of granular suspensions.

In particular, a DST transition characterized by S-shaped curves becomes more pronounced as the mass ratio

m / m_{g}

increases. Specifically, the non-Newtonian shear viscosity

η^{*}

exhibits a discontinuous transition (at a certain value of

a^{*}

) which intensifies as the particles of the granular gas become heavier. The theoretical results are validated with MD simulations [18] in the Brownian limiting case (

m / m_{g} \to 0

), showing generally good agreement despite slight discrepancies in the transition zone. Simulations suggest a sharper transition, likely due to the absence of molecular chaos in highly nonequilibrium situations. To address this, DSMC simulations for IHS are performed in the same limit, showing good agreement with theoretical results and further emphasizing a more pronounced transition. This phenomenon is likely attributable to a sudden growth of higher-order moments, resulting in a proportional increase in the deviation from theoretical predictions. As a result, the lower moments are also affected. Some technical details of the application of the DSMC method are available in the supplementary material of Ref. [15].

The simplicity of the BGK and IMM models enables exploration beyond second-degree moments. Accordingly, we utilize the BGK-type kinetic equation to compute the fourth-degree moments. Although similar calculations could be performed in the case of IMM, we opt to omit them due to their extensive analytical effort. Additionally, drawing insights from the Fokker–Planck model [26] and dry granular gases, we anticipate potential divergences of the moments derived from IMM under certain shear rate conditions. We focus our efforts on calculating the following fourth-degree moments:

M_{4 | 0}^{*} = \frac{m^{2}}{n T_{g}^{2}} \int d V V^{4} f (V),

(99)

M_{2 | x y}^{*} = \frac{m^{2}}{n T_{g}^{2}} \int d V V^{2} V_{x} V_{y} f (V) .

(100)

Thus, in terms of the generic moments

M_{k_{1}, k_{2}, k_{3}}

and according to the expression (79), the moments

M_{4 | 0}^{*}

and

M_{2 | x y}^{*}

are given by

M_{4 | 0}^{*} = \frac{m^{2}}{n T_{g}^{2}} (M_{4, 0, 0} + M_{0, 4, 0} + M_{0, 0, 4} + 2 M_{2, 2, 0} + 2 M_{2, 0, 2} + 2 M_{0, 2, 2}),

(101)

M_{2 | x y}^{*} = \frac{m^{2}}{n T_{g}^{2}} (M_{3, 1, 0} + M_{1, 3, 0} + M_{1, 1, 2}) .

(102)

Figure 3a illustrates the ratio

M_{4 | 0}^{*} (a^{*}) / M_{4 | 0}^{*} (0)

as a function of

a^{*}

for

α = 0.9

and 1 and three values of the mass ratio. We observe that variations in the mass ratio do not significantly alter the trends observed in the Brownian limiting case [26]. An abrupt transition in the higher-order moments is evident within a small region of

a^{*}

. Specifically, the kurtosis

M_{4 | 0}^{*}

increases with the mass ratio

m / m_{g}

until it converges to the value obtained in the Brownian limit. Consistent with the conclusions drawn in Ref. [17], an increase in the mass of the particles of the granular gas results in an elevation of the granular temperature. Consequently, energy nonequipartition accentuates and moves the suspension away from equilibrium, leading to an increase in kurtosis as the distribution function deviates from its Maxwellian form. Regarding the influence of collisional dissipation, we observe that the effect of

α

M_{4 | 0}^{*}

remains relatively discrete. Figure 3b illustrates the shear-rate dependence of the (reduced) moment

M_{2 | x y}^{*}

. This moment vanishes in the absence of shear rate (

a^{*} = 0

). Similar conclusions to those made for the moment

M_{4 | 0}^{*}

can be drawn.

Theoretical predictions for the fourth-degree moments are compared against DSMC simulations for IHS conducted in this paper in the Brownian limiting case. A qualitative agreement is observed, although simulations suggest a sharper transition. Some quantitative discrepancies are noticeable, which are mainly disguised by the scale. To assess the reliability of the BGK-type results, we focus on the region

0 < a^{*} < 1

, where all the fourth-degree velocity moments of IMM are well-defined functions of the shear rate. In addition, non-Newtonian effects are still significant in the range of values of the (reduced) shear rate

a^{*} \leq 1

. To this purpose, Figure 4 shows the (reduced) fourth-degree moments

M_{4 | 0}^{*} (a^{*}) / M_{4 | 0}^{*} (0)

and

M_{2 | x y}^{*} (a^{*})

for

α = 0.7

and 1. These moments are also illustrated as obtained for IMM in the Brownian limit [26]. It is worth noting that the results derived in Ref. [26] stem from considering an effective force modeling the interstitial gas, diverging from the limit of a Boltzmann–Lorentz operator modeled by a BGK-type equation as

m / m_{g} \to \infty

. Consequently, since DSMC simulations employ the exact Fokker–Planck operator, they perfectly align with the IMM results, while discrepancies emerge when compared with the BGK-type results. It is noteworthy that the BGK-type model slightly overestimates the deviation from the Newtonian situation (

a^{*} = 0

), a phenomenon also observed for molecular gases [37]. Moreover, non-Newtonian effects become apparent even at low values of

a^{*}

Finally, in Figure 5, the ratio

R_{x} (c_{x}) = φ_{x} (c_{x}) / (π^{- 1 / 2} e^{- c_{x}^{2}})

is plotted for

a^{*} = 1

and four different values of the mass ratio. Here, the marginal distribution

φ_{x} (c_{x})

is given by Equation (95). It is evident that the deviation from equilibrium (

R_{x} \neq 1

) becomes more significant as the mass ratio

m / m_{g}

increases. Moreover, a comparison between theory and DSMC simulations reveals some disagreement in the BGK-type solution. Although the relative difference of these discrepancies is relatively small (it is about 8%), this contradicts what was observed in Ref. [9], where good agreement between the BGK solution and DSMC data was shown in the region of thermal velocities.

6. Concluding Remarks

In our study, we explored the non-Newtonian transport properties of a dilute granular suspension subjected to USF using the Boltzmann kinetic equation. The particles are represented as d-dimensional hard spheres with mass m and diameter

σ

, immersed in an interstitial gas acting as a thermostat at temperature

T_{g}

. Various models for granular suspensions incorporate a gas–solid force to represent the influence of the external fluid. While some models consider only isolated body resistance via a linear drag law [5,6,8,10,11,12,38,39], others [13,40] include an additional Langevin-type stochastic term. In this paper, we consider a suspension model where the collisions between grains and particles of the interstitial (molecular) gas are taken into account. Thus, based on previous studies [14,41], we discretize the surrounding molecular gas, assigning individual particles with mass

m_{g}

and diameter

σ_{g}

, thereby accounting for elastic collisions between grains and background gas particles in the starting kinetic equation.

Under USF conditions, the system is characterized by constant density profiles n and

n_{g}

, uniform temperatures T and

T_{g}

, and a (common) flow velocity

U_{x} = U_{g, x} = a y

, where a denotes the shear rate. In agreement with previous investigations on uniform sheared suspensions, the mean flow velocity

U

is coupled to that of the gas phase

U_{g}

. Consequently, the viscous heating term due to shear and the energy transferred by the grains from collisions with the molecular gas is compensated by the cooling terms derived from collisional dissipation, allowing the achievement of a steady state. A distinctive feature of the USF is that the one-particle velocity distribution function

f (r, v)

depends on space only through its dependence on the peculiar velocity

V = v - U

. Consequently, the velocity distribution function becomes uniform in the Lagrangian reference frame, moving with

V

. This means that

f (r, v) \equiv f (V)

. Based on symmetry considerations, the heat flux

q

vanishes, making the pressure tensor

P

the relevant flux. Therefore, to understand the intricate dynamics of granular suspensions under shear flow, it is imperative to focus on their non-Newtonian properties—derived from the pressure tensor. These include the (reduced) temperature

χ = T / T_{g}

, the (reduced) nonlinear shear viscosity

η^{*}

, and the (reduced) normal stress difference

Ψ^{*}

Given that the most challenging aspect of dealing with the Boltzmann equation lies in the collision operator, it is reasonable to explore alternatives that render this operator more analytically tractable than in the case of IHS. Among the most sophisticated techniques in this regard is to consider the Boltzmann equation for IMM. As in the case of elastic collisions [20,22], the collision rate for IMM is independent of the relative velocity of the colliding particles. As a consequence, the collisional moments of degree k of the Boltzmann collisional operator can be expressed as a linear combination of velocity moments of degree less than or equal to k. To complement the results derived for IMM, we also considered in this paper the use of a BGK-type kinetic model where the true Boltzmann operator is replaced by a simple relaxation term. Here, we employed both approaches to compute the rheological properties of the sheared granular suspension. Thus, our objective was twofold. Firstly, we aimed to assess the reliability and compatibility of the proposed models with previous results [17] obtained for IHS using Grad’s moment method. Additionally, DSMC simulations for IHS were performed as an alternative method to validate any potential discrepancies identified. Secondly, taking advantage of the capabilities provided by the BGK model, we endeavored to calculate the velocity distribution function and the higher-order moments that offer insights into its characteristics.

Before proceeding with the computation of rheological properties, it is necessary to understand the response of the molecular fluid to shear stress. This assessment was also conducted using both the Maxwell molecules and BGK-type kinetic model that were later used to model the granular gas. A novelty here is the application of a force (Gaussian thermostat) of the form

F = - m ξ V

to compensate for the energy gained through viscous shear stresses. This force, by consistency, also applies to the granular gas, maintaining convergence to a steady state. As anticipated, the results agree well with those obtained through Grad’s moment method [17] (see Equations (32) and (64)). Consequently, once the problem conditions (including the shear rate a) are defined, the molecular temperature

T_{g}

is determined, effectively serving as a thermostat for the granular gas.

After accurately describing the rheology of the molecular gas, we focused on modeling the granular gas. Using both the IMM and BGK-type model separately, we calculated the nonzero elements of the the pressure tensor. The knowledge of these elements allowed us to identify the relevant rheological properties of the granular suspension. As shown in Figure 2, these quantities are represented as functions of the coefficient of restitution

α

and the mass ratio

m / m_{g}

. In particular, we find that the theoretical results obtained from the Grad’s method for IHS, IMM, and BGK-type model show remarkable agreement, with almost indistinguishable curves. This underlines the effectiveness of structurally simple models in capturing the complexities of sheared granular suspensions. We observe a DST-type transition starting at a certain value of

a^{*}

, which increases with the mass ratio

m / m_{g}

. Interestingly, similar to the MD simulations performed for IHS [18], the DSMC data suggest a more abrupt transition than predicted by theory. Given that the main divergences between Grad’s (for IHS) and DSMC results arise from the form of the distribution function, the significance of investigating higher-order moments to assess the deviation of the distribution function from its Maxwellian reference is then justified.

Based on the previous literature where discrepancies in fourth-degree moments have been observed [26], and acknowledging the potential lengthiness of calculations, for the sake of simplicity, we decided to employ only the BGK-type model to compute higher-order moments. Specifically, we concentrated on the (symmetric) fourth-degree moments

M_{4 | 0}

and

M_{2 | x y}

. The shear-rate dependence of these moments is illustrated in Figure 3 for the same parameter values of

α

and

m / m_{g}

as those employed for the rheological quantities. Initially, we note that the fourth-degree moments also exhibit an abrupt transition at a value of

a^{*}

that increases with

m / m_{g}

until reaching the Brownian limit. We think that the DST behavior will also appear in all higher moments.

Furthermore, DSMC simulations in the Brownian limit qualitatively capture the profile of these fourth-degree moments, although some quantitative disparities are apparent. To ascertain the extent of these discrepancies, we narrowed our focus to the interval

0 < a^{*} < 1

, where non-Newtonian effects are apparent. Additionally, we included IMM results directly as obtained from an effective Fokker–Planck-type model [26]. Figure 4 illustrates that BGK results overestimate the deviation from the moments computed when no shear stress is applied compared to DSMC simulations and the results obtained using an effective force to model the interstitial gas. These disparities are also observed in the marginal distribution function

φ_{x}

The theoretical findings presented here motivate the comparison with computer simulations. Although the observed agreement in the Brownian limit is encouraging, there is scope to extend this agreement to scenarios with finite mass ratios. Our plan is to carry out simulations of this type in the near future, which we expect will further validate and improve our theoretical framework. In addition, we plan to extend our current findings to finite densities by exploring the Enskog kinetic equation, which will allow us to evaluate the involvement of density in the occurrence of these phenomena. Recent results [42] derived in the context of the Enskog equation by using the Fokker–Planck operator have shown that there is a transition from the discontinuous shear thickening (observed in dilute gases) to the continuous shear thickening for denser systems. We want to see if this behavior is also present for large but finite mass ratios. The above lines of research will be some of the main objectives of our upcoming projects.

Author Contributions

Conceptualization, R.G.G. and V.G.; software, validation, R.G.G. and V.G.; formal analysis, R.G.G. and V.G.; writing—original draft preparation, V.G.; writing—review and editing, R.G.G. and V.G. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge financial support from Grant PID2020-112936GB-I00 funded by MCIN/AEI/10.13039/501100011033, and from Grant IB20079 funded by Junta de Extremadura (Spain) and by ERDF “A way of making Europe”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors report no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:

IMM	Inelastic Maxwell model
IHS	Inelastic hard spheres
USF	Uniform shear flow
DST	Discontinuous shear thickening
DSMC	Direct simulation Monte Carlo
BGK	Bhatnagar–Gross–Krook kinetic model
GK	Gross–Krook kinetic model

Appendix A. Some Technical Details in the BGK-Type Kinetic Model

In this Appendix, we give some details on the determination of the parameters

\tilde{T}

and

ν

appearing in the relaxation term (66). To obtain them, we require the collisional transfer of energy of grains due to their elastic collisions with particles of the molecular gas to be the same as that obtained from the true Boltzmann–Lorentz collision operator. This implies that

\int d v V^{2} J_{BL} [V | f, f_{g}] = - ν \int d v V^{2} [f (V) - {\tilde{f}}_{g} (V)] .

(A1)

Given that the collisional moment involving the operator

J_{BL} [f, f_{g}]

cannot be exactly computed for IHS, one estimates this moment by replacing f and

f_{g}

with their Grad’s solutions [19]. In this approximation, one achieves the results [17]

\int d v V^{2} J_{BL} [V | f, f_{g}] = - \frac{8 π^{(d - 1) / 2}}{Γ (\frac{d}{2})} n n_{g} {\bar{σ}}^{d - 1} \frac{m_{g}}{{(m + m_{g})}^{2}} {(\frac{2 T_{g}}{m_{g}} + \frac{2 T}{m})}^{1 / 2} (T - T_{g}) .

(A2)

Moreover,

\int d v V^{2} [f (V) - {\tilde{f}}_{g} (V)] = - d \frac{n ν}{m} (T - \tilde{T}) .

(A3)

From Equations (A2) and (A3), one obtains the identity

\frac{8 π^{(d - 1) / 2}}{d Γ (\frac{d}{2})} n_{g} {\bar{σ}}^{d - 1} \frac{m m_{g}}{{(m + m_{g})}^{2}} {(\frac{2 T_{g}}{m_{g}} + \frac{2 T}{m})}^{1 / 2} (T - T_{g}) = ν (T - \tilde{T}) .

(A4)

Equation (A4) allows us to make the identifications (70).

Finally, from Grad’s moment method [19], the collisional moment (71) can be written as [17]

\int d v m V_{i} V_{j} J [f, f] = - ν^{M} n T_{g} [ν_{η}^{* IHS} P_{i j}^{*} - χ (ν_{η}^{* IHS} - ζ^{* IHS}) δ_{i j}],

(A5)

where

ν^{M}

is defined in Equation (42) and

ν_{η}^{* IHS} = \frac{(1 + α) [3 (1 - α) + 2 d]}{2 d (d + 2)}, ζ^{* IHS} = \frac{1 - α^{2}}{2 d} .

(A6)

The BGK-type kinetic model (65) yields the result

\int d v m V_{i} V_{j} [- ν^{'} (f - f^{M}) + \frac{ϵ}{2} \frac{\partial}{\partial V} \cdot V f] = - n T_{g} [(ν^{'} + ϵ) P_{i j}^{*} - ν^{'} χ δ_{i j}],

(A7)

where we recall that

P_{i j}^{*} = P_{i j} / n T_{g}

. Comparison between Equations (A5) and (A7) yields Equation (72) for

ν^{'} (α)

References

Koch, D.L. Kinetic theory for a monodisperse gas-solid suspension. Phys. Fluids A 1990, 2, 1711–1722. [Google Scholar] [CrossRef]
Gidaspow, D. Multiphase Flow and Fluidization; Academic Press: Cambridge, MA, USA, 1994. [Google Scholar]
Fullmer, W.D.; Hrenya, C.M. The clustering instability in rapid granular and gas-solid flows. Annu. Rev. Fluid Mech. 2017, 49, 485–510. [Google Scholar] [CrossRef]
Jackson, R. The Dynamics of Fluidized Particles; Cambridge University Press: New York, NY, USA, 2000. [Google Scholar]
Tsao, H.K.; Koch, D.L. Simple shear flows of dilute gas-solid suspensions. J. Fluid Mech. 1995, 296, 211–245. [Google Scholar] [CrossRef]
Sangani, A.S.; Mo, G.; Tsao, H.K.; Koch, D.L. Simple shear flows of dense gas-solid suspensions at finite Stokes numbers. J. Fluid Mech. 1996, 313, 309–341. [Google Scholar] [CrossRef]
Wylie, J.J.; Zhang, Q.; Li, Y.; Hengyi, X. Driven inelastic-particle systems with drag. Phys. Rev. E 2009, 79, 031301. [Google Scholar] [CrossRef] [PubMed]
Heussinger, C. Shear thickening in granular suspensions: Interparticle friction and dynamically correlated clusters. Phys. Rev. E 2013, 88, 050201. [Google Scholar] [CrossRef] [PubMed]
Chamorro, M.G.; Vega Reyes, F.; Garzó, V. Non-Newtonian hydrodynamics for a dilute granular suspension under uniform shear flow. Phys. Rev. E 2015, 92, 052205. [Google Scholar] [CrossRef]
Saha, S.; Alam, M. Revisiting ignited-quenched transition and the non-Newtonian rheology of a sheared dilute gas-solid suspension. J. Fluid Mech. 2017, 833, 206–246. [Google Scholar] [CrossRef]
Alam, M.; Saha, S.; Gupta, R. Unified theory for a sheared gas-solid suspension: From rapid granular suspension to its small-Stokes-number limit. J. Fluid Mech. 2019, 870, 1175–1193. [Google Scholar] [CrossRef]
Saha, S.; Alam, M. Burnett-order constitutive relations, second moment anisotropy and co-existing states in sheared dense gas-solid suspensions. J. Fluid Mech. 2020, 887, A9. [Google Scholar] [CrossRef]
Garzó, V.; Tenneti, S.; Subramaniam, S.; Hrenya, C.M. Enskog kinetic theory for monodisperse gas-solid flows. J. Fluid Mech. 2012, 712, 129–168. [Google Scholar] [CrossRef]
Biben, T.; Martin, P.; Piasecki, J. Stationary state of thermostated inelastic hard spheres. Physica. A 2002, 310, 308–324. [Google Scholar] [CrossRef]
Gómez González, R.; Garzó, V. Kinetic theory of granular particles immersed in a molecular gas. J. Fluid Mech. 2022, 943, A9. [Google Scholar] [CrossRef]
Subramaniam, S. Multiphase flows: Rich physics, challenging theory, and big simulations. Phys. Rev. Fluids 2020, 5, 110520. [Google Scholar] [CrossRef]
Gómez González, R.; García Chamorro, M.; Garzó, V. Rheology of granular particles immersed in a molecular gas under uniform shear flow. arXiv 2023, arXiv:2311.16717. [Google Scholar]
Hayakawa, H.; Takada, S. Kinetic theory of discontinuous rheological phase transition for a dilute inertial suspension. Prog. Theor. Exp. Phys 2019, 2019, 083J01. [Google Scholar] [CrossRef]
Grad, H. On the kinetic theory of rarefied gases. Commun. Pure Appl. Math. 1949, 2, 331–407. [Google Scholar] [CrossRef]
Ikenberry, E.; Truesdell, C. On the pressures and the flux of energy in a gas according to Maxwell’s kinetic theory. J. Rat. Mech. Anal. 1956, 5, 1–54. [Google Scholar] [CrossRef]
Chapman, S.; Cowling, T.G. The Mathematical Theory of Nonuniform Gases; Cambridge University Press: Cambridge, UK, 1970. [Google Scholar]
Cercignani, C. The Boltzmann Equation and Its Applications; Springer: New York, NY, USA, 1988. [Google Scholar]
Ben-Naim, E.; Krapivsky, P.L. The Inelastic Maxwell Model. In Granular Gas Dynamics; Pöschel, T., Luding, S., Eds.; Springer: Berlin/Heidelberg, Germany, 2003; Volume 624, pp. 65–94. [Google Scholar]
Vega Reyes, F.; Garzó, V.; Santos, A. Granular mixtures modeled as elastic hard spheres subject to a drag force. Phys. Rev. E 2007, 75, 061306. [Google Scholar] [CrossRef] [PubMed]
Brown, E.; Jeager, H.M. Dynamic jamming point for shear thickening suspensions. Rep. Prog. Phys. 2014, 77, 046602. [Google Scholar] [CrossRef]
Gómez González, R.; Garzó, V. Simple shear flow in granular suspensiones: Inelastic Maxwell models and BGK-type kinetic model. J. Stat. Mech. 2019, 2019, 013206. [Google Scholar] [CrossRef]
Pelargonio, S.; Zaccone, A. Generalized Langevin equation with shear flow and its fluctuation-dissipation theorems derived from a Caldeira-Leggett Hamiltonian. Phys. Rev. E 2023, 107, 064102. [Google Scholar] [CrossRef]
Bird, G.A. Molecular Gas Dynamics and the Direct Simulation Monte Carlo of Gas Flows; Clarendon: Oxford, UK, 1994. [Google Scholar]
Evans, D.J.; Morriss, G.P. Statistical Mechanics of Nonequilibrium Liquids; Academic Press: London, UK, 1990. [Google Scholar]
Brilliantov, N.; Pöschel, T. Kinetic Theory of Granular Gases; Oxford University Press: Oxford, UK, 2004. [Google Scholar]
Garzó, V. Granular Gaseous Flows; Springer Nature: Cham, Switzerland, 2019. [Google Scholar]
Dufty, J.W.; Santos, A.; Brey, J.J.; Rodríguez, R.F. Model for nonequilibrium computer simulation methods. Phys. Rev. A 1986, 33, 459–466. [Google Scholar] [CrossRef]
Brey, J.J.; Dufty, J.W.; Santos, A. Kinetic models for granular flow. J. Stat. Phys. 1999, 97, 281–322. [Google Scholar] [CrossRef]
Bhatnagar, P.L.; Gross, E.P.; Krook, M. A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component system. Phys. Rev. 1954, 94, 511–525. [Google Scholar] [CrossRef]
Santos, A.; Astillero, A. System of elastic hard spheres which mimics the transport properties of a granular gas. Phys. Rev. E 2005, 72, 031308. [Google Scholar] [CrossRef] [PubMed]
Gross, E.P.; Krook, M. Model for collision processes in gases. Small amplitude oscillations of charged two-component systems. Phys. Rev. 1956, 102, 593–604. [Google Scholar] [CrossRef]
Garzó, V.; Santos, A. Kinetic Theory of Gases in Shear Flows. Nonlinear Transport; Springer: Dordrecht, The Netherlands, 2003. [Google Scholar]
Parmentier, J.F.; Simonin, O. Transition models from the quenched to ignited states for flows of inertial particles suspended in a simple sheared viscous fluid. J. Fluid Mech. 2012, 711, 147–160. [Google Scholar] [CrossRef]
Wang, T.; Grob, M.; Zippelius, A.; Sperl, M. Active microrhelogy of driven granular particles. Phys. Rev. E 2014, 89, 042209. [Google Scholar] [CrossRef]
Puglisi, A.; Loreto, V.; Marconi, U.M.B.; Petri, A.; Vulpiani, A. Clustering and non-Gaussian behavior in granular matter. Phys. Rev. Lett. 1998, 81, 3848–3851. [Google Scholar] [CrossRef]
Santos, A. Granular fluid thermostated by a bath of elastic hard spheres. Phys. Rev. E 2003, 67, 051101. [Google Scholar] [CrossRef] [PubMed]
Takada, S.; Hayakawa, H.; Santos, A.; Garzó, V. Enskog kinetic theory of rheology for a moderately dense inertial suspension. Phys. Rev. E 2020, 102, 022907. [Google Scholar] [CrossRef] [PubMed]

Figure 1. Plot of the nonlinear shear viscosity

η_{g}^{*}

and the normal stress difference

Ψ_{g}^{*}

for hard spheres (

d = 3

) as functions of the (reduced) shear rate

\tilde{a}

. Symbols refer to the DSMC results for hard spheres.

Figure 1. Plot of the nonlinear shear viscosity

η_{g}^{*}

and the normal stress difference

Ψ_{g}^{*}

for hard spheres (

d = 3

) as functions of the (reduced) shear rate

\tilde{a}

. Symbols refer to the DSMC results for hard spheres.

Figure 2. Plots of the (steady) granular temperature

χ

(a), the non-Newtonian shear viscosity

η^{*}

(b), and the normal stress difference

Ψ^{*}

a^{*}

for two different values of the coefficient of restitution

α

: 0.9 (left panel) and 1 (right panel). The graphs represent four distinct mass ratio values

m / m_{g}

10^{3}

(yellow lines),

10^{4}

(blue lines),

10^{6}

(red lines), and the Brownian limit (black lines). Here,

T_{g}^{*} = 1

d = 3

, and

ϕ = 0.0052

. The solid lines correspond to the IMM results, the dashed lines are the BGK-like results, and the dotted lines refer to Grad’s solution for IHS. Symbols denote computer simulation results performed in the Brownian limit: circles refer to the DSMC data obtained in this paper for IHS, while squares are MD results obtained in Ref. [18] for IHS.

Figure 2. Plots of the (steady) granular temperature

χ

(a), the non-Newtonian shear viscosity

η^{*}

(b), and the normal stress difference

Ψ^{*}

a^{*}

for two different values of the coefficient of restitution

α

: 0.9 (left panel) and 1 (right panel). The graphs represent four distinct mass ratio values

m / m_{g}

10^{3}

(yellow lines),

10^{4}

(blue lines),

10^{6}

(red lines), and the Brownian limit (black lines). Here,

T_{g}^{*} = 1

d = 3

, and

ϕ = 0.0052

Figure 3. Plots of the (reduced) fourth-degree moments

M_{4 | 0}^{*} (a^{*}) / M_{4 | 0}^{*} (0)

(a) and

- M_{2 | x y}^{*} (a^{*})

(b) as functions of the (reduced) shear rate

a^{*}

obtained from the BGK-type equation for two different values of the coefficient of restitution

α

: 0.9 (left panel) and 1 (right panel). The graphs represent four distinct mass ratio values

m / m_{g}

10^{3}

(yellow lines),

10^{4}

(blue lines),

10^{6}

(red lines), and the Brownian limit (black lines). Here,

T_{g}^{*} = 1

d = 3

, and

ϕ = 0.0052

. Symbols refer to the DSMC results obtained for IHS in this paper in the Brownian limit.

Figure 3. Plots of the (reduced) fourth-degree moments

M_{4 | 0}^{*} (a^{*}) / M_{4 | 0}^{*} (0)

(a) and

- M_{2 | x y}^{*} (a^{*})

(b) as functions of the (reduced) shear rate

a^{*}

obtained from the BGK-type equation for two different values of the coefficient of restitution

α

: 0.9 (left panel) and 1 (right panel). The graphs represent four distinct mass ratio values

m / m_{g}

10^{3}

(yellow lines),

10^{4}

(blue lines),

10^{6}

(red lines), and the Brownian limit (black lines). Here,

T_{g}^{*} = 1

d = 3

, and

ϕ = 0.0052

. Symbols refer to the DSMC results obtained for IHS in this paper in the Brownian limit.

Figure 4. Plots of the (reduced) fourth-degree moments

M_{4 | 0}^{*} (a^{*}) / M_{4 | 0}^{*} (0)

(a) and

- M_{2 | x y}^{*} (a^{*})

(b) as functions of the (reduced) shear rate

a^{*}

obtained from the BGK-type equation for two different values of the coefficient of restitution

α

: 0.7 (solid lines) and 1 (dashed lines). The graphs represent four distinct mass ratio values

m / m_{g}

: 5 (yellow lines), 10 (blue lines), 50 (red lines), and the Brownian limit (black lines). Here,

T_{g}^{*} = 1

d = 3

, and

ϕ = 0.0052

. Symbols refer to the DSMC results for IHS in the Brownian limit (squares for

α = 1

and circles for

α = 0.7

). The green lines are the IMM results as obtained in Ref. [26] in the Brownian limit.

Figure 4. Plots of the (reduced) fourth-degree moments

M_{4 | 0}^{*} (a^{*}) / M_{4 | 0}^{*} (0)

(a) and

- M_{2 | x y}^{*} (a^{*})

(b) as functions of the (reduced) shear rate

a^{*}

obtained from the BGK-type equation for two different values of the coefficient of restitution

α

: 0.7 (solid lines) and 1 (dashed lines). The graphs represent four distinct mass ratio values

m / m_{g}

: 5 (yellow lines), 10 (blue lines), 50 (red lines), and the Brownian limit (black lines). Here,

T_{g}^{*} = 1

d = 3

, and

ϕ = 0.0052

. Symbols refer to the DSMC results for IHS in the Brownian limit (squares for

α = 1

and circles for

α = 0.7

). The green lines are the IMM results as obtained in Ref. [26] in the Brownian limit.

Figure 5. Plot of the ratio

R_{x} (c_{x}) = φ_{x} (c_{x}) / (π^{- 1 / 2} e^{- c_{x}^{2}})

for

a^{*} = 1

as a function of the (reduced) velocity

c_{x}

for

α = 0.9

and five different values of the mass ratio

m / m_{g}

: 1 (yellow line), 5 (blue line), 10 (red line), 50 (green line), and the Brownian limit (black lines). Here,

T_{g}^{*} = 1

d = 3

, and

ϕ = 0.0052

. The dashed line refers to the DSMC results for IHS in the Brownian limit.

Figure 5. Plot of the ratio

R_{x} (c_{x}) = φ_{x} (c_{x}) / (π^{- 1 / 2} e^{- c_{x}^{2}})

for

a^{*} = 1

as a function of the (reduced) velocity

c_{x}

for

α = 0.9

and five different values of the mass ratio

m / m_{g}

: 1 (yellow line), 5 (blue line), 10 (red line), 50 (green line), and the Brownian limit (black lines). Here,

T_{g}^{*} = 1

d = 3

, and

ϕ = 0.0052

. The dashed line refers to the DSMC results for IHS in the Brownian limit.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Gómez González, R.; Garzó, V. Exact Results for Non-Newtonian Transport Properties in Sheared Granular Suspensions: Inelastic Maxwell Models and BGK-Type Kinetic Model. Entropy 2024, 26, 265. https://doi.org/10.3390/e26030265

AMA Style

Gómez González R, Garzó V. Exact Results for Non-Newtonian Transport Properties in Sheared Granular Suspensions: Inelastic Maxwell Models and BGK-Type Kinetic Model. Entropy. 2024; 26(3):265. https://doi.org/10.3390/e26030265

Chicago/Turabian Style

Gómez González, Rubén, and Vicente Garzó. 2024. "Exact Results for Non-Newtonian Transport Properties in Sheared Granular Suspensions: Inelastic Maxwell Models and BGK-Type Kinetic Model" Entropy 26, no. 3: 265. https://doi.org/10.3390/e26030265

APA Style

Gómez González, R., & Garzó, V. (2024). Exact Results for Non-Newtonian Transport Properties in Sheared Granular Suspensions: Inelastic Maxwell Models and BGK-Type Kinetic Model. Entropy, 26(3), 265. https://doi.org/10.3390/e26030265

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Exact Results for Non-Newtonian Transport Properties in Sheared Granular Suspensions: Inelastic Maxwell Models and BGK-Type Kinetic Model

Abstract

1. Introduction

2. Boltzmann Kinetic Equation for Sheared Granular Suspensions

3. Rheology from Inelastic Maxwell Models

3.1. Rheological Properties of the Molecular Gas

3.2. Rheological Properties of the Granular Gas

3.3. Brownian Limit

4. Rheology from a BGK-Type Kinetic Model of the Boltzmann Equation

4.1. Rheological Properties of the Molecular Gas

4.2. Rheological Properties of the Granular Gas

4.3. Brownian Limit

4.4. Velocity Distribution of the Granular Gas

5. Comparison between IMM and BGK Results

6. Concluding Remarks

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Conflicts of Interest

Abbreviations

Appendix A. Some Technical Details in the BGK-Type Kinetic Model

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI