Asymptotic Solutions of the Kinetic Boltzmann Equation and Multicomponent Non-Equilibrium Gas Dynamics

Affiliation(s)

^{1}
Fundamental Researches Department, Russian Federal Nuclear Centre, All-Russian Scientific Research
Institute of Experimental Physics, Sarov, Russia.

^{2}
St. Petersburg State University, St. Petersburg, Russia.

ABSTRACT

In the article correct method for the kinetic Boltzmann equation asymptotic solution is formulated, the Hilbert’s and Enskog’s methods are discussed. The equations system of multicomponent non- equilibrium gas dynamics is derived, that corresponds to the first order in the approximate (asym- ptotic) method for solution of the system of kinetic Boltzmann equations.

1. Introduction

In 1912 Hilbert considered the kinetic Boltzmann equation for one-component gas as an example of integral equation and proposed a “recipe” for its approximate (asymptotic) solution (see [1], Chapter~XXII). Hilbert’s “recipe” was inconvenient for practical use, because the five arbitrary functional parameters of the first and the following approximations of the velocity distribution function had to be found by solving the differential equations in partial derivatives (equations of gas dynamics of the first and higher orders). Five years later Enskog proposed to use zero conditions, conditions with zero right-hand sides, to determine the five arbitrary functional parameters of the first and following approximations of the velocity distribution function. The imposition of the zero conditions leads, in fact, to using different comparison scales in the asymptotic expansion of the velocity distribution function and in the asymptotic expansion of the particle number density, the mean (mass) velocity and the temperature, that are derived from the asymptotic expansion of the velocity distribution function by integration over velocities with different weighting functions. As a result of paralogism of the method of successive approximations (one has to set variable coefficients of the same terms of the unified comparison scale equal to each other) partial time derivatives vanish in the necessary conditions of solutions existence of integral equations of higher orders (see below) and with them terms of gas-dynamic equations, corresponding to viscosity, heat conduction, … vanish. Enskog “improved” the situation by the introducing (see, for example, [2], Chapter 7, § 1, Section 5) of the unsubstantiated expansion of partial time derivative:

(1)

The approach of Struminskii, who had proposed in 1974 in [3] his approximate (asymptotic) method of solution of the system of kinetic Boltzmann equations for multicomponent gas, differs from the approach of Enskog to asymptotic solution of the Boltzmann equations system for gas mixture in that, how the infinitesimal parameter is introduced in the Boltzmann equations system for gas mixture, i.e. the solution is constructing in another asymptotic limit. In substance, Struminskii’s method of solution of kinetic equations system is the same as Enskog’s method (Struminskii used the partial time derivative expansion, as Enskog did).

In section 2 below will be proposed the correct method of asymptotic solution of the kinetic Boltzmann equations system for multicomponent gas mixture for the approach, that combines Enskog’s and Struminskii’s approaches; in particular, it will be shown, how one has to modify Enskog’s method: in addition to asymptotic expansion of the velocity distribution function i-component particles of gas mixture it is necessary to determine and to use the expansion of the particle number density of i-component, mean mass velocity and temperature T of the gas mixture.

Further, in the Section 3 the system of infinitesimal first order equations of multicomponent non-equilibrium gas dynamics, appearing during the process of the solution of the system of Boltzmann equations by successive approximations method in the Section 2 as necessary condition of the existence of approximate (asymptotic) solution of the integral equations system, is considered in more detail.

This article is condensed version of our article arXiv:1303.6275. Notations, used below, are close to notations in [2]; it is assumed, that all regarded functions are continuous and continuously differentiable so many times as it is necessary, if their derivatives are considered, and all regarded integrals converge.

2. Correct Method of Solution of the Kinetic Boltzmann Equations System

The Boltzmann equations system, that describes change of dependent on t and spatial coordinates, prescribed by radius-vector, the velocity distribution functions due to collision with particles of other components of mixture of rarefied monatomic gases, where are the velocities of particles of i-component of the mixture {see [2], Chapter 8, Equation (1.1); discussion of the derivation of the Boltzmann equations system and its applicability range see, for example, in [2], Chapters 3 and 18, [4], Chapter 7, § 1; below the central interaction of molecules are considered only, when the force, with which each molecule acts on the other, is directed along the line, connecting the centers of the molecules}, could be written as:

(2)

in (2) N is a set of indexes, that are numbering components of the mixture; is an external force, which acts on the molecule of the i-component; is the mass of the molecule of the i-component; is the modulus of the relative velocity of colliding particles; b is the impact distance, is the azimuth angle, is the unit vector, directed to the center of mass of the colliding particles from the point of closest approach―see [2], Chapter 3, Figure 3; the scalar function is determined by equality

(3)

by prime in (20) and below the velocities and the functions of velocities after the collision are denoted.

Let us introduce following notations:

(4)

(5)

to differ velocities of colliding molecules of the same kind in (22) the one velocity is denoted by and the other is denoted by (without any index) and the index of the corresponding velocity distribution function f is omitted.

In Enskog’s approach the differential parts of the Boltzmann Equations (2), that are denoted by below, are considered to be small as compared with the right-hand sides of Equations (2)―see [2], Chapter 7, § 1, Section 5―therefore the indicator of infinite smallness is formally introduced in the Boltzmann equations system in the following way:

(6)

In Struminskii’s approach to the asymptotic solution of the Boltzmann equations system the differential parts of the Boltzmann Equations (2) and the collision integrals of the particles of i-component with the particles of the other components are considered to be small as compared with the collision integral of the particles of i-component between each other, therefore the indicator of infinitesimality is introduced in the Boltzmann equations system in another way:

(7)

It is possible to combine Enskog’s approach with Struminskii’s approach. For this purpose we divide the set of mixture components N into two subsets: the subset of components, that we call formally inner components (we could consider the case, when there are some subsets of inner components, but this case does not fundamentally differ from the one, considered below, the only difference is that the notation become more complicated) and the subset of components, that we call external components. To differ the two groups of mixture components we

denote the subset of indexes of inner components as well as the indexes of inner components and the subset of indexes of external components as well as the indexes of external components; the

intersection of the sets and is the empty set―and the union of these sets is the set of indexes of all mixture components; if an assertion concerns both kinds of components the special symbols will be omitted. In new notations the Boltzmann equations system can be rewritten as:

(8)

(9)

Let us write the asymptotic expansion of the velocity distribution function of particles of i-component as formal series of successive approximations in powers of:

(10)

The differential parts of the Equations (3) are written as:

(11)

where

(12)

―cf. with [2], Chapter 7, § 1, Sections 4, 5 and [3]. In (11)-(12) the partial time derivative expansion (1) is not used in contrast to that, how it was made by Enskog and further by Struminskii. As result, described below method for solution of the system of kinetic Boltzmann equations differ fundamentally from Enskog’s method and Struminskii’s method.

Substituting (10) and (11) in (8) and equating coefficients at the same powers of to each other, we obtain the equations system of the method of successive approximations for finding the velocity distribution functions of inner components particles of gas mixture; taking introduced notations (4), (5) and (12) into account, the system can be rewritten as:

(13)

(14)

Similarly substituting (10) and (11) in (9) and equating coefficients at the same powers of to each other, we obtain the equations system of the method of successive approximations for finding the velocity distribution

functions of particles of external components of gas mixture:

(15)

(16)

Speaking about an order of approximation below, we assume the order to be equal to the value of index r in (14), (16). According to (5), (13), in zero order approximation we have the following system of integral equations to find the velocity distribution functions of particles of inner components of gas mixture:

(17)

The general solution of the equations system (17) can be written as a set of the Maxwell functions:

(18)

where k is the Boltzmann constant.

Particle number density of the i-component, mean mass velocity and temperature T of inner com- ponents of mixture are introduced by definitions:

(19)

(20)

(21)

in (21) is the Boltzmann constant. From (19)-(21) the equality is obtained:

(22)

that is convenient to use below instead of definition (21).

According to definitions (19), (20), (21), in addition to the asymptotic expansion (10) it is necessary to determine asymptotic expansions for particle number density of the i-component

(23)

mean mass velocity

(24)

and temperature of inner components of mixture

(25)

Substituting (10) and (23)-(25) in (19), (20), (22) and equating terms of the same infinitesimal order we obtain scalar relations, that connect asymptotic expansions (10) and (23)-(25):

(26)

(27)

(28)

In (27), (28) the notations are introduced

(29)

(30)

In particular, for from (26)-(28) we obtain expressions for arbitrary functions, and in (18) through the zero order approximations to local values of the -component number

density, the mean mass velocity and the temperature of inner components of the mixture:

(31)

(32)

(33)

According to (4), (15), zero order integral equations, from which the velocity distribution functions of particles of outer components of the mixture are found:

(34)

―are simpler than Equations (17) and differ actually from (17) only by lack of summation over components. Therefore, similarly (18), the general solution of the equations system (34) can be written as a set of the Maxwell functions:

(35)

where and are some, independent of, scalar functions of spatial coordinates, defined by the radius vector, and time t, and is a vector function of and t.

Let’s add to the definition of the number density of particles of i-component definitions of mean velocity and temperature of outer component of mixture:

(36)

(37)

from (19), (36), (37) the equality is obtained:

(38)

that is convenient to use below instead of definition (37).

Let’s enter similar (24)-(25) asymptotic expansions of outer -component mean velocity

(39)

and outer -component temperature

(40)

Substituting (10), (23), (39), (40) in (19), (36), (38) and equating terms of the same infinitesimal order we obtain for each 5 (scalar) relations, that connect asymptotic expansions (10), (23), (39), (40):

(41)

(42)

(43)

cf. with (26)-(28). In (42), (43) the notation is used

(44)

For from (41)-(43) we obtain expressions for arbitrary functions, and in (35) through the zero order approximations to local values of the number density, the mean

velocity and the temperature of outer -component of the mixture:

(45)

(46)

(47)

For the velocity distribution functions of inner components of gas mixture are found from the integral equations system (14), which, taking (5) and equality

(48)

into account, can be rewritten in the form

(49)

in (49) functions are written as, where are new unknown functions.

The left-hand sides of Equations (49) involves functions, that are known from the previous step of the successive approximations method. Unknown functions enter, linearly, only into the right-hand sides of Equations (49). Therefore the general solution of the system of Equations (14) is a family of functions of a form

, where, , a family of functions is a particular solution of the system of inhomogeneous Equations (49) and a family of functions is the general solution of the system of homogeneous equations

(50)

Multiplying Equations (50) by, integrating over all values of, summing over and transforming integrals, we obtain

(51)

From (51) we conclude, that are linear combinations of the summational invariants of the collision :

(52)

where and are some, independent of, scalar functions of spatial coordinates, defined by the radius vector, and time t, and is a vector function of and t (as well as above, arbitrary functions and are identical for all inner components of the mixture), and, hence,

(53)

To simplify further evaluations according to the expression for, see (18) and (31)-(33), let us rewrite (53) as

(54)

where, and are new functions of and t. Family of functions is a solution of the system of inhomogeneous equations

(55)

where denote left-hand sides of the Equations (49), taken with opposite sign.

Multiplying Equations (55) by , integrating over all values of and transforming integrals as above, we obtain as necessary condition for the existence of solutions of the system of integral Equations (55), the necessity of the fulfillment of equalities:

(56)

(57)

Among (infinitesimal) set of particular solutions of the system of Equations (55), different from each other on some solution of the system of homogeneous Equations (50), unique solution may be chosen such that

(58)

(59)

Having substituted expression for

(60)

in (26)-(28), taking (18), (29)-(33) and (58)-(59) into account, we obtain a system of algebraic equations [constraint equations for asymptotic expansions (10) and (23)-(25)]:

(61)

(62)

(63)

from which we find expressions for functions, and through (variable) coe-

fficients of asymptotic expansions of the particle number density of -component, of the mean mass velocity and of the temperature of inner components of the mixture

(64)

(65)

(66)

Then the fulfillment of equalities (56)-(57) can be considered as the differential equations, the r-order equations of gas dynamics, for finding, ,.

The partial solution of the system of inhomogeneous Equations (55), satisfying (58)-(59), may be constructed, for example, using expansion of in series in terms of Sonine polynomials with expansion

coefficients, depending on and t (see [2] or [4]); such construction proves existence of solutions of the system of integral Equations (49).

For the velocity distribution functions of outer components of gas mixture may be similarly found from the integral equations system (16):

(67)

where

(68)

(69)

(70)

The fulfillment of analogous (56)-(57) equalities

(71)

can be considered as the differential equations, the r-order equations of gas dynamics, for finding, ,.

3. The System of First Order Equations of Multicomponent Non-Equilibrium Gas Dynamics

Let us consider in more detail the system of infinitesimal first order Equations (56)-(57), (71), derived above as the necessary (and sufficient) condition of the solution existence of the first order integral equations system (14), (16).

To simplify transformations, according to the expressions for velocity distribution functions of particles of infinitesimal zero order (18), (35), functions, may be used in (56)-(57), (71) rather than functions, , respectively:

(72)

(73)

(74)

for inner components, for outer components.

At transformation of differential parts of the Equations (56)-(57) and (71) we use equalities:

(75)

(76)

(77)

In (75)-(77) the bar above symbol with index denotes the average of the value:

(78)

and are considered as independent variables. At averaging in (77) it is assumed, that external force, acting on the particle of species i, is independent of the particle velocity, it is assumed also, that integrals,

depending on external forces, are convergent, and product tends to zero, when tends to infinity.

After simple transformations from (56)-(57) and (71) we obtain following system of infinitesimal first order equations of multicomponent non-equilibrium gas dynamics:

(79)

(80)

(81)

(82)

(83)

(84)

In accordance with the general definition of pressure tensor of i-component of gas mixture

(85)

and with the general definition of i-component heat flux vector

(86)

(cf. with [2], Chapter 2, §§ 3, 4) in (79)-(84)

(87)

is inner components pressure tensor of zero order, is inner components hydrostatic pressure of zero order, is the unit tensor, double product of two second rank tensors and ([2], Chapter 1, § 3) is the scalar,

(88)

is inner components heat flux vector of zero order,

(89)

is zero order internal energy of particles of inner components per unit volume, which is equal, in this case, to energy of their translational chaotic motion, however, the energy transfer equations, written in form (81) and (84) can be used in more general cases as well (cf. with [4], Chapter 7, § 6), in (87)-(89) averaging (78) is performed with Maxwell function from (18);

(90)

is -component pressure tensor of zero order, is -component hydrostatic pressure of zero order,

(91)

is -component heat flux vector of zero order,

(92)

is zero order internal energy of particles of -component per unit volume, in (90)-(92) averaging (78) is performed with Maxwell function from (35).

General analytic expressions for integrals, from (80), (81) and (83), (84), that depend on the interaction cross-section, can be derived in general case, when separate components (with Maxwell velocity distribution function of particles) have different mean velocities and temperatures.

System of infinitesimal first order equations of multicomponent non-equilibrium gas dynamics (79)-(84) is proposed to use for describing turbulent flows

Cite this paper

Serov, S. and Serova, S. (2016) Asymptotic Solutions of the Kinetic Boltzmann Equation and Multicomponent Non-Equilibrium Gas Dynamics.*Journal of Applied Mathematics and Physics*, **4**, 1687-1697. doi: 10.4236/jamp.2016.48177.

Serov, S. and Serova, S. (2016) Asymptotic Solutions of the Kinetic Boltzmann Equation and Multicomponent Non-Equilibrium Gas Dynamics.

References

[1] Hilbert, D. (1912) Grundzüge einer Allgemeinen Theorie der Linearen Integralgleichungen. Teubner, Leipzig and Berlin. (In Ger-man)

[2] Chapman, S. and Cowling, T.G. (1952) The Mathematical Theory of Non-uniform Gases. Cambridge University Press, Cambridge.

[3] Struminskii, V.V. (1974) Influence of Diffusion Velocity on Flow of Gas Mixtures. Prikladnaya Mathematica i Me-chanica [Applied Mathematics and Mechanics], 38, 203-210. (In Russian)

[4] Hirschfelder, J.O., Curtiss, Ch.F. and Bird, R.B. (1954) Molecular Theory of Gases and Liquids. Wiley, New York.

[1] Hilbert, D. (1912) Grundzüge einer Allgemeinen Theorie der Linearen Integralgleichungen. Teubner, Leipzig and Berlin. (In Ger-man)

[2] Chapman, S. and Cowling, T.G. (1952) The Mathematical Theory of Non-uniform Gases. Cambridge University Press, Cambridge.

[3] Struminskii, V.V. (1974) Influence of Diffusion Velocity on Flow of Gas Mixtures. Prikladnaya Mathematica i Me-chanica [Applied Mathematics and Mechanics], 38, 203-210. (In Russian)

[4] Hirschfelder, J.O., Curtiss, Ch.F. and Bird, R.B. (1954) Molecular Theory of Gases and Liquids. Wiley, New York.