Research Article | Open Access

# Boltzmann’s Six-Moment One-Dimensional Nonlinear System Equations with the Maxwell-Auzhan Boundary Conditions

**Academic Editor:**Peter G. L. Leach

#### Abstract

We prove existence and uniqueness of the solution of the problem with initial and Maxwell-Auzhan boundary conditions for nonstationary nonlinear one-dimensional Boltzmann’s six-moment system equations in space of functions continuous in time and summable in square by a spatial variable. In order to obtain a priori estimation of the initial and boundary value problem for nonstationary nonlinear one-dimensional Boltzmann’s six-moment system equations we get the integral equality and then use the spherical representation of vector. Then we obtain the initial value problem for Riccati equation. We have managed to obtain a particular solution of this equation in an explicit form.

#### 1. Introduction

In case of one-atom gas any macroscopic system during process of its evolution to an equilibrium state passes 3 stages: initial transition period, described in terms of a full function of system distribution; the kinetic period, by means of a one-partial distribution function; and the hydrodynamic period, by means of five first moments of the distribution function. In kinetic regime the behavior of a rarefied gas in space of time and velocity is described by Boltzmann’s equation. It is known from gas dynamics that in most of the encountered problems there is no need to use detailed microscopic gas description with help of the distribution function. Therefore it is natural to look for a less detailed description using macroscopic hydrodynamic variables (density, hydrodynamic velocity, temperature, etc.). As these variables are defined in terms of moments of the distribution function, we are faced with the problem of analyzing the various moments of Boltzmann’s equation.

Note that Boltzmann’s moment equations are intermediate between Boltzmann (kinetic theory) and hydrodynamic levels of description of state of the rarefied gas and form class of nonlinear partial differential equations. Existence of such class of equations was noticed by Grad [1, 2] in 1949. He obtained the moment system by expanding the particle distribution function in Hermite polynomials near the local Maxwell distribution. Grad used Cartesian coordinates of velocities and Grad’s moment system contained as coefficients such unknown hydrodynamic characteristics like density, temperature, average speed, and so forth. Formulation of boundary conditions for Grad’s system is almost impossible, as the characteristic equations for various approximations of Grad’s hyperbolic system contain unknown parameters like density, temperature, and average speed. However, 13- and 20-moment Grad equations are widely used in solving many problems of the kinetic theory of gases and plasma.

Boltzmann’s equation is equivalent to an infinite system of differential equations relative to the moments of the particle distribution function in the complete system of eigenfunctions of linearized operator. As a rule we limit study by finite moment system of equations because solving infinite system of equations does not seem to be possible.

Finite system of moment equations for a specific task with a certain degree of accuracy replaces Boltzmann’s equation. It is necessary, also roughly, to replace the boundary conditions for the particle distribution function by a number of macroscopic conditions for the moments; that is, there arises the problem of boundary conditions for a finite system of equations that approximate the microscopic boundary conditions for Boltzmann’s equation. The question of boundary conditions for a finite system of moment equations can be divided into two parts: how many conditions must be imposed and how they should be prepared. From microscopic boundary conditions for Boltzmann’s equation there can be obtained an infinite set of boundary conditions for each type of decomposition. However, the number of boundary conditions is determined not by the number of moment equations; that is, it is impossible, for example, to take as many boundary conditions as equations, although the number of moment equations affects the number of boundary conditions. In addition, the boundary conditions must be consistent with the moment equations and the resulting problem must be correct. Boundary condition problems arise in the following tasks: (1) moment boundary conditions in rarefied gas slip-flow problems; (2) definition of boundary conditions on the surfaces of streamlined rarefied gas; (3) prediction aerospace aerodynamic characteristics of aircraft at very high speeds and at high altitudes and so forth.

In work [3] we have obtained the moment system which differs from Grad’s system of equations. We used spherical velocity coordinates and decomposed the distribution function into a series of eigenfunctions of the linearized collision operator [4, 5], which is the product of Sonine polynomials and spherical functions. The resulting system of equations, which correspond to the partial sum of series and which we called Boltzmann’s moment system of equations, is a nonlinear hyperbolic system in relation to the moments of the particles distribution function.

The structure of Boltzmann’s moment system of equations corresponds to the structure of Boltzmann’s equation; namely, the differential part of the resulting system is linear in relation to the moments of the distribution function and nonlinearity is included as moments of collision integral [6].

The linearity of differential part of Boltzmann’s moment system of equations simplifies the task of formulation of the boundary conditions. In work [3] a homogeneous boundary condition for particles distribution function was approximated and proved the correctness of initial and boundary value problem for nonlinear nonstationary Boltzmann’s moment system of equations in three-dimensional space. In work [7] the initial and boundary value problem for one-dimensional nonstationary Boltzmann’s equation with boundary conditions of Maxwell was approximated by a corresponding problem for Boltzmann’s moment system of equations. The boundary conditions for Boltzmann’s moment system of equations were called Maxwell-Auzhan conditions.

In work [8] a systematic nonperturbative derivation of a hierarchy of closed systems of moment equations corresponding to any classical theory has been presented. This paper is a fundamental work where closed systems of moment equations describe a transition regime. Moreover, hydrodynamical model is used to describe charge transport in a generic compound semiconductor. Compound semiconductors have found wide use in the microelectronic industry. The evolution equation for macroscopic variables is obtained by taking moments of the transport equation [9–13].

The study of various problems for Boltzmann’s moment system of equations is an important and actual task in the theory of a rarefied gas and other applications of the moment system equations. The correctness of initial and boundary problems for Boltzmann’s moment system of equations with Maxwell-Auzhan boundary conditions is being studied for the first time.

#### 2. Existence and Uniqueness of the Solutions of Initial and Boundary Value Problem for Six-Moment One-Dimensional Boltzmann’s System of Equations with Maxwell-Auzhan Boundary Conditions

In this section we prove the existence and uniqueness of solutions of the initial and boundary value problem for six-moment one-dimensional Boltzmann’s system of equations with Maxwell-Auzhan boundary conditions in space of functions, continuous in time and summable in square by spatial variable. Note that the theorem of the existence of a global solution in time of the initial and boundary value problem for 3-dimensional nonlinear Boltzmann’s equation with boundary conditions of Maxwell is proved in work [14].

We write in an expanded form system of one-dimensional Boltzmann’s moment equations in the* k*th approximation, which corresponds to decomposition of the particle distribution function by eigenfunctions of the linearized collision operator where the moments of a nonlinear collision operator are expressed through coefficients of Talmi and Klebsh-Gordon [15, 16] and have the form where are generalized coefficients of Talmi, are Klebsh-Gordon coefficients, is the constant, is Boltzmann’s constant, and is the ideal gas temperature.

For generalized coefficients of Talmi there exists a table [15] for each value of quantum number from 0 to 6. Moreover, there is a program on IBM for calculation of generalized coefficients of Talmi.

If in (1) takes values from 0 to 3, we get Boltzmann’s moment system equations in the third approximation. We write it in an expanded form: where and are the moments of the particle distribution function; , , and are moments of the collision integral, where , , , and are the Fourier coefficients of the cross section expansion by the Legendre polynomials. The first, third, and fourth equations of system (3) correspond to mass conservation law, momentum conservation law, and energy conservation law, respectively.

We study the correctness of initial and boundary value problem for six-moment one-dimensional Boltzmann’s system equations where is the transpose matrix; is the positive definite matrix; and are given initial vector-functions; are vector moments of falling to the boundary particle distribution function; are vector moments of reflecting from the boundary particle distribution function. Equation (4) is vector matrix form of system equations (3).

It is possible to check through direct calculations thatand matrix has three positive and the same number of negative nonzero eigenvalues. From (6)-(7) it follows that the number of boundary conditions on the left and right ends of interval is equal to the number of positive and negative eigenvalues of the matrix .

Thus system (4) is a symmetric hyperbolic nonlinear partial differential equations system. Let us show that is a sign-nondefined square form. It is easy to check that , where andEigenvalues of matrix are −1/2, , 0, 0, 0, and 1/2. Therefore is a sign-nondefined square form. Similarly, we can show that and are also sign-nondefined square forms.

For problem (4)–(7) the following theorem takes place.

Theorem 1. *If , then problem (4)–(7) has a unique solution in domain belonging to the space ; moreoverwhere is a constant independent from and is the particular solution of Riccati equation (15) in an explicit form.*

*Proof. *Let . Let us prove estimation (11). We multiply the first equation in system (4) by and the second equation by and perform integration from to : After integration by parts we getTaking into account boundary conditions (6)-(7) we rewrite equality (13) in the following form:where .

Let us use the spherical representation [17] of vector , where , , and . Furthermore assume that , where the value of we define below.

Substituting the values and into (14) we have thatwhereLet us study (15) with an initial condition Let be a particular solution of Riccati equation (15). Then the general solution of (15) has the formHence, taking into consideration initial condition (17) we obtain If , then the second term of is bounded under the condition that difference is positive. Further we suppose a difference is a positive number. Let . We denote by the moment of time at which Then is bounded , where and . Analyzing structure of (see (31)) from equality (20) we obtain for any .

Now we found one particular solution of (15). The nonlinear Riccati equation can always be reduced to a second-order linear ordinary differential equation [18]whereSolution of this equation will lead to a solution of of the original Riccati equation. By substitution we write (21) in the following form [18]: where .

If instead of and we substitute values from (22), then we getIn work [19] it was proven that function is the solution of (24), where and the series on the right hand side of (26) uniformly converge.

Thenis a particular solution of Riccati equation (15).

Since , value ofNow instead of we substitute its value and get .

The value of with issince for any

An expansion form of the function It is not difficult to prove using immediate substitution that function satisfies Riccati equation (15).

From equality (19) we obtainHenceWe take as interval of the solution of problem (15)–(17) segment , since the integrand is bounded. Hence (where ) a priori estimation (11) takes place.

Now we prove the existence of solution for (4)–(7) using Galerkin method. Let be a basis in space , where the dimension of vector is equal to the dimension of vector . For each we define an approximate solution of (4)–(7) as follows: where is the orthogonal projection in of the function on the subspace, spanned by , We represent in form , where The coefficients are determined from the equationswhere is the th component of .

We multiply (35) by and sum over from to :With help of the above shown arguments we now prove that is bounded in some time interval , where , , , and where is constant and independent from and is a particular solution of Riccati equation about .

Problem (40)-(41) represents a Cauchy problem for the ordinary system of differential equations. Existence of the solution of problem (40)-(41) follows from theory of ordinary system of differential equations (Picard’s existence and uniqueness theorem) [20]. Hence the existence of the solution of problem (34)–(37) follows.

Thus, according to estimation (43) the sequence of approximate solutions of problem (5)–(7) is uniformly bounded in function space . Moreover, the homogeneous system of equations with respect to has only a trivial solution. Then it follows from results in [21] that is week in and is week in as . Further it can be shown with standard methods that limit element is a weak solution of problem (5)–(7).

The uniqueness of solution of problem (5)–(7) is proved by contradiction. Let problem (5)–(7) have two different solutions and . We denote them again by and . Then with respect to new values of and we obtain the following problem: We prove that solution of problem (44)–(46) is trivial. Hence the uniqueness of the solution of problem (5)–(7) follows.

Using the method which was mentioned before we get (see equality (14))Let us transform the integrand on the right side of the equality asOnce again we use spherical representation of and , where . Then concerning we obtain new initial value problemwhere has the same value as in (15) andThe general solution of (49) isSolution of (49), which is satisfying homogeneous condition (50), is trivial; that is, . Hence, and and .

The theorem is proved.

#### Competing Interests

The authors declare that they have no competing interests.

#### References

- H. Grad, “On the kinetic theory of rarefied gases,”
*Communications on Pure and Applied Mathematics*, vol. 2, no. 4, pp. 331–407, 1949. View at: Publisher Site | Google Scholar - H. Grad, “Principle of the kinetic theory of gases,” in
*Thermodynamics of Gases*, vol. 12 of*Handbuch der Physik*, pp. 205–294, Springer, Berlin, Germany, 1958. View at: Google Scholar - A. Sakabekov,
*Initial-Boundary Value Problems for the Boltzmann's Moment System Equations Gylym, Almaty*, 2002. - C. Cercignani,
*Theory and Application of the Boltzmann Equation*, Instituto di Matematica, Milano, Italy, 1975. - M. N. Kogan,
*Dynamic of Rarefied Gas*, Nauka, Moscow, Russia, 1967. - K. Kumar, “Polynomial expansions in kinetic theory of gases,”
*Annals of Physics*, vol. 37, no. 1, pp. 113–141, 1966. View at: Publisher Site | Google Scholar - A. Sakabekov and Y. Auzhani, “Boundary conditions for the onedimensional nonlinear nonstationary Boltzmann's moment system equations,”
*Journal of Mathematical Physics*, vol. 55, Article ID 123507, 2014. View at: Google Scholar - C. D. Levermore, “Moment closure hierarchies for kinetic theories,”
*Journal of Statistical Physics*, vol. 83, no. 5-6, pp. 1021–1065, 1996. View at: Publisher Site | Google Scholar | Zentralblatt MATH - G. Mascali and V. Romano, “A hydrodynamical model for holes in silicon semiconductors: the case of non-parabolic warped bands,”
*Mathematical and Computer Modelling*, vol. 53, no. 1-2, pp. 213–229, 2011. View at: Publisher Site | Google Scholar - G. Mascali and V. Romano, “A non parabolic hydrodynamical subband model for semiconductors based on the maximum entropy principle,”
*Mathematical and Computer Modelling*, vol. 55, no. 3-4, pp. 1003–1020, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH - V. D. Camiola and V. Romano, “2DEG-3DEG charge transport model for MOSFET based on the maximum entropy principle,”
*SIAM Journal on Applied Mathematics*, vol. 73, no. 4, pp. 1439–1459, 2013. View at: Publisher Site | Google Scholar - G. Alì, G. Mascali, V. Romano, and R. C. Torcasio, “A hydrodynamical model for covalent semiconductors with a generalized energy dispersion relation,”
*European Journal of Applied Mathematics*, vol. 25, no. 2, pp. 255–276, 2014. View at: Publisher Site | Google Scholar | Zentralblatt MATH - V. D. Camiola and V. Romano, “Hydrodynamical model for charge transport in graphene,”
*Journal of Statistical Physics*, vol. 157, no. 6, pp. 1114–1137, 2014. View at: Publisher Site | Google Scholar - S. Mischler, “Kinetic equations with Maxwell boundary conditions,”
*Annales Scientifiques de l'Ecole Normale Superieure*, vol. 43, no. 5, pp. 719–760, 2010. View at: Google Scholar - V. G. Neudachin and U. F. Smirnov,
*Nucleon Association of Easy Kernel*, Nauka, Moscow, Russia, 1969. - M. Moshinsky,
*The Harmonic Oscillator in Modern Physics: from Atoms to Quarks*, 1960. - S. I. Pokhozhaev, “On an approach to nonlinear equation,”
*Doklady Akademii Nauk SSSR*, vol. 247, pp. 1327–1331, 1979. View at: Google Scholar - E. Kamke,
*Differentialgleichungen Lösungsmethoden und Lösungen*, I. Gewohuliche Differentialgleichungen, B.G.Teubner, Leipzig, Germany, 1977. View at: Publisher Site - A. Tungatarov and D. K. Akhmed-Zaki, “Cauchy problem for one class of ordinary differential equations,”
*International Journal of Mathematical Analysis*, vol. 6, no. 14, pp. 695–699, 2012. View at: Google Scholar - M. Tenenbaum and H. Pollard,
*Ordinary Differential Equations*, Harper and Row, New York, NY, USA, 1963. - L. Tartar, “Compensated compactness and applications to partial differential equations,” in
*Proceedings of the Non-Linear Analysis and Mechanics, Heriot-Watt Symposium*, R. J. Knops, Ed., vol. 4 of*Research Notes in Math*, pp. 136–212, 1979. View at: Google Scholar

#### Copyright

Copyright © 2016 A. Sakabekov and Y. Auzhani. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.