2.1. Liouville Equation and Reduced Forms
In Cartesian vector notation, the Liouville equation for particle number preserving systems reads [
11]
where
is the N-particle density function,
and
are the position and momentum coordinates for the ith particle of mass
, and we have assumed that the force acting on particle
i is only a function of the position coordinates
).
Now consider a reduced form of the Liouville equation for a set of molecules
that can be obtained by integrating the Liouville equation over the phase-space of the other
set of molecules. Following standard procedures [
11], we can integrate Equation (
1) over
space to obtain the
reduced Liouville equation
The force acting on the
i-th molecule can be written as the gradient of a potential as
and for the purposes of specific applications below, the potential will be approximated by a sum over the pair interaction potential as (“pairwise additivity”)
where the pair potential
is the interaction potential between any two molecules in the system, i.e., the effects of three or more body interactions on the pair potential expression are neglected. Please note that this assumption does not limit the generality of our results given below.
Substituting Equation (
4) into Equation (
2) and integration on the right-hand side over all the (
) space except the (
) molecule gives
Equation (
6) is the reduced Liouville equation for pairwise additive interaction forces, and it is an integro-differential equation where the evolution of the
density depends on the next higher-order,
, density. This nonhomogeneous dimensionality feature is the so-called BBGKY hierarchy named after its originators: Bogoliubov, Born, Green, Kirkwood, and Yvon [
11].
2.2. Entropy Conservation
Irving and Kirkwood derived the mass, momentum, and energy conservation equations from the Liouville equation, and their paradigm will be followed here unless otherwise noted. For more details the reader is referred to the original manuscript [
2]. To obtain an entropy conservation equation following the IK (Irving and Kirkwood) approach, we define a general quantity
as [
3]
where, following Green [
12],
is the natural logarithm of an
s-th ordered, normalized density function that depends on the multiparticle expansion method [
12,
13]. Please note that
is an arbitrary locator vector in the system. The phase-space average of
,
, is computed by integrating
over all s phase-space coordinates as
We call this average the Boltzmann–Gibbs entropy, where
is the entropy per molecule and
n is the local molecular number density. For
,
and we recover Boltzmann’s definition [
3,
4]
where
ℏ is Planck’s constant, and in the limit
, it can be shown that we recover Gibbs’ entropy [
12]
We note that for systems at global equilibrium
is independent of the locator vector
and, using the fact that
at global equilibrium, Equation (
9) becomes equal to Gibbs’ equilibrium entropy function.
Extensions to other orders can also be shown. Specifically, the nonequilibrium, two-space Green’s entropy would be defined according to Equation (
7) as [
14]
where
, which again under global equilibrium conditions can readily be shown to reduce to Green’s two-particle equilibrium expression [
14]. Please note that various multiparticle expansion or closure methods for the s-order density functions have been recently studied and reviewed by Singer [
13].
Now, the dynamical variable for entropy given above depends explicitly on time through the s-order density functions, whereas the dynamical variables for mass, momentum, and energy introduced by IK have no explicit time dependence. So, to obtain an entropy conservation equation following the IK approach, we can work with the Liouville equation and dynamic variable modified to include the explicit time dependence in
, or it is somewhat easier and equivalent to work directly with the reduced Liouville equation, Equation (
6); we choose the latter approach.
Multiplying Equation (
6) by
and integrating over all
space gives the following
entropy conservation equation [
3]
where
n is the local number density and the bulk or mass average velocity
is defined by
and the entropy flux vector
also follows as
which represents the flux of entropy relative to the bulk average velocity, where
.
Now, the last term on the right-hand side represents the entropy generation, specifically
Please note that for
, we obtain the well-known Boltzmann’s entropy generation term [
11],
where we have used Newton’s Third Law in writing this last expression. Please note that further simplification of Equation (
14) is possible, but our focus here will be on the
generation term and its particular forms.
2.3. Asymptotic Expansions
We now look at the specific expressions for the entropy flux and entropy generation terms based on asymptotic expansions. The truncation of these expansions will be clearly shown to result in finite entropy generation; consequently, the truncations must be consistent across the transport equations for any given system. For the sake of generality, we introduce the following dimensionless variables denoted by an asterisk into the reduced Liouville equation: ; ; ; ; ; where and are characteristic values of length, momentum, and intermolecular interaction force, respectively, where the characteristic force is expressed as a characteristic interaction energy, , divided by a characteristic interaction length scale, .
Substituting these dimensionless quantities into Equation (
5) leads to the non-dimensionless, reduced Liouville equation [
14]
where
and
are dimensionless groups.
For most problems, the characteristic velocity would be the molecular velocity and, thus,
, where
is a characteristic temperature. Also, the characteristic intermolecular potential,
, is usually on the order of magnitude of
, so we can select
and
giving
Specifically in Equation (
16) above,
is the length scale over which an appreciable change in
with respect to
occurs and
is the length scale over which an appreciable change in
with respect to
occurs. For example, in most gas and liquid systems
, and this disparity is used to obtain approximate asymptotic solutions to the reduced Liouville equation. Since appreciable changes in
are concomitant with changes in its moments (or the transport properties),
is also referred to as characteristic macroscopic length scale. We will also assume
to be on the order of one for a moderate density gas. The characteristic time is normally associated with the macroscopic scale and
may be taken to be unity for our purposes here. Please note that the scaling analysis given above is a generalization of that given by Frieman [
15].
In both gases and liquids
is a smallness parameter. Under these conditions, the solution to Equation (
16) is written in the following regular or Poincare’ perturbation expansion
and it follows that
Using the Chapman–Enskog expansion method, it was previously shown that the entropy flux and generation expressions can be obtained from their definitions given above [
3]. The expressions include only the so-called kinetic contributions, which is consistent with the
level of analysis. Here we demonstrate a more general method for obtaining the entropy flux and generation expressions. For completeness, we also note that in the case of rarefied and dilute gases the expansion is not “regular”, but can still be readily treated as originally shown by Frieman [
14,
15]. Also, note that a complete analysis for gases requires both and
and
analysis, the latter of which includes two body interactions and leads to the “potential part” of the flux relationships [
11,
14]. However, the paradigm given below is readily extended to a complete
analysis and will not be considered here in order to simplify the presentation and results.
For
, substitution of Equation (
20) into Equation (
16) leads to a hierarchy of order equations [
14]:
Now, through straightforward analysis of the
expressions, it can be readily shown that the leading order solution to Equation (
22) follows the
local equilibrium density [
14], given in
dimensionless terms as
where
is the dimensionless momentum of molecule j relative to the bulk motion and
is the dimensionless two particle density function [
14]. The relations given above are all that are needed to derive general forms of the entropy flux and entropy generation in gases, as we will now show.