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

Authors: A. Sakabekov; Y. Auzhani
Journal: Journal of Applied Mathematics
 (2016)

---

## 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.


---

## Body

## 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 thekth approximation, which corresponds to decomposition of the particle distribution function by eigenfunctions of the linearized collision operator (1)
∂
f
n
l
∂
t
+
1
α
∂
∂
x
l
2
n
+
l
+
1
/
2
2
l
-
1
2
l
+
1
f
n
,
l
-
1
-
2
n
+
1
2
l
-
1
2
l
+
1
f
n
+
1
,
l
-
1
+
l
+
1
2
n
+
l
+
1
/
2
2
l
+
1
2
l
+
3
f
n
,
l
+
1
-
2
n
2
l
+
1
2
l
+
3
f
n
-
1
,
l
+
1
=
I
n
l
,
2
n
+
l
=
0,1
,
…
,
k
,where the moments I
n
l of a nonlinear collision operator are expressed through coefficients of Talmi and Klebsh-Gordon [15, 16] and have the form (2)
I
n
l
=
∑
N
3
L
3
n
3
l
3
:
l
∣
n
l
00
:
l
N
3
L
3
n
3
l
3
:
l
∣
n
1
l
1
n
2
l
2
:
l
l
1
0
l
2
0
l
0
σ
l
3
-
σ
0
f
n
1
l
1
f
n
2
l
2
,where N
3
L
3
n
3
l
3
:
l
∣
n
1
l
1
n
2
l
2
:
l are generalized coefficients of Talmi, l
1
0
l
2
0
/
l
0 are Klebsh-Gordon coefficients, α
=
1
/
R
θ is the constant, R 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 ξ
=
2
n
+
l from 0 to 6. Moreover, there is a program on IBM for calculation of generalized coefficients of Talmi.If in (1) 2
n
+
l takes values from 0 to 3, we get Boltzmann’s moment system equations in the third approximation. We write it in an expanded form: (3)
∂
f
00
∂
t
+
1
α
∂
f
01
∂
x
=
0
,
∂
f
02
∂
t
+
1
α
∂
∂
x
2
3
f
01
+
3
5
f
03
-
2
2
15
f
11
=
J
02
,
∂
f
10
∂
t
+
1
α
∂
∂
x
-
2
3
f
01
+
5
3
f
11
=
0
,
∂
f
01
∂
t
+
1
α
∂
∂
x
f
00
+
2
3
f
02
-
2
3
f
10
=
0
,
∂
f
03
∂
t
+
1
α
∂
∂
x
3
5
f
02
=
J
03
,
∂
f
11
∂
t
+
1
α
∂
∂
x
-
2
2
15
f
02
+
5
3
f
10
=
J
11
,
x
∈
-
a
,
a
,
t
>
0
,where f
00
=
f
00
t
,
x and f
01
=
f
01
t
,
x
,
…
,
f
11
=
f
11
t
,
x are the moments of the particle distribution function; J
02
=
σ
2
-
σ
0
f
00
f
02
-
f
01
2
/
3
/
2, J
03
=
(
1
/
4
)
σ
3
+
3
σ
1
-
4
σ
0
f
00
f
03
+
(
1
/
4
5
)
(
2
σ
1
+
σ
0
-
3
σ
3
)
f
01
f
02, and J
11
=
σ
1
-
σ
0
(
f
00
f
11
+
(
1
/
2
)
5
/
3
f
10
f
01
-
(
2
/
15
)
f
01
f
02
) are moments of the collision integral, where σ
0, σ
1, σ
2, and σ
3 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(4)
∂
u
∂
t
+
A
∂
w
∂
x
=
J
1
u
,
w
∂
w
∂
t
+
A
′
∂
u
∂
x
=
J
2
u
,
w
,
t
∈
0
,
T
,
x
∈
-
a
,
a
,
(5)
u
t
=
0
=
u
0
x
,
w
t
=
0
=
w
0
x
,
x
∈
-
a
,
a
,
(6)
A
w
-
+
B
u
-
x
=
-
a
=
1
β
A
w
+
-
B
u
+
x
=
-
a
-
1
-
β
η
α
β
π
F
,
t
∈
0
,
T
,
(7)
A
w
-
-
B
u
-
x
=
a
=
1
β
A
w
+
+
B
u
+
x
=
a
+
1
-
β
η
α
β
π
F
,
t
∈
0
,
T
,where(8)
A
=
1
α
1
0
0
2
3
3
5
-
2
2
15
-
2
3
0
5
3
,
B
=
1
α
π
2
2
3
-
1
3
2
3
2
2
-
1
-
1
3
-
1
3
2
J
1
u
,
w
=
0
,
J
02
,
0
′
,
J
2
u
,
w
=
0
,
J
03
,
J
11
′
,
u
=
f
00
,
f
02
,
f
10
′
,
w
=
f
01
,
f
03
,
f
11
′
,
F
=
1
4
2
,
1
8
6
,
1
8
3
′
.
A
′ is the transpose matrix; B is the positive definite matrix; u
0
x
=
f
00
0
x
,
f
02
0
x
,
f
10
0
x
′ and w
0
x
=
f
01
0
x
,
f
03
0
x
,
f
11
0
x
′ are given initial vector-functions; w
+
,
u
+ are vector moments of falling to the boundary particle distribution function; w
-
,
u
-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 that(9)
det
⁡
A
1
=
det
⁡
0
A
A
′
0
≠
0
,and matrix A
1 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 -
a
,
a is equal to the number of positive and negative eigenvalues of the matrix A
1.Thus system (4) is a symmetric hyperbolic nonlinear partial differential equations system. Let us show that J
02 is a sign-nondefined square form. It is easy to check that f
00
f
02
-
(
f
01
2
/
3
)
=
(
C
U
,
U
), where U
=
u
,
w
′ and(10)
C
=
0
1
2
0
0
0
0
1
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-
1
3
0
0
0
0
0
0
0
0
0
0
0
0
0
0
.Eigenvalues of matrix C are −1/2, -
1
/
3, 0, 0, 0, and 1/2. Therefore J
02 is a sign-nondefined square form. Similarly, we can show that J
03 and J
11 are also sign-nondefined square forms.For problem (4)–(7) the following theorem takes place.Theorem 1.
IfU
0
=
(
u
0
x
,
w
0
x
)
∈
L
2
[
-
a
,
a
], then problem (4)–(7) has a unique solution in domain -
a
,
a
×
[
0
,
T
] belonging to the space C
0
,
T
;
L
2
-
a
,
a; moreover(11)
U
L
2
-
a
,
a
-
r
1
C
0
,
T
≤
C
1
U
0
L
2
-
a
,
a
-
r
1
0
,where C
1 is a constant independent from U and T
~
Ο
U
0
L
2
-
a
,
a
-
r
1
(
0
)
)
-
1
,
r
1
t is the particular solution of Riccati equation (15) in an explicit form.Proof.
LetU
0
∈
L
2
[
-
a
,
a
]. Let us prove estimation (11). We multiply the first equation in system (4) by u and the second equation by w and perform integration from -
a to a: (12)
1
2
d
d
t
∫
-
a
a
u
,
u
+
w
,
w
d
x
+
∫
-
a
a
A
∂
w
∂
x
,
u
+
A
′
∂
u
∂
x
,
w
d
x
=
∫
-
a
a
J
1
,
u
+
J
2
,
w
d
x
.
After integration by parts we get(13)
1
2
d
d
t
∫
-
a
a
u
,
u
+
w
,
w
d
x
+
u
-
,
A
w
-
x
=
a
-
u
-
,
A
w
-
x
=
-
a
=
∫
-
a
a
J
1
,
u
+
J
2
,
w
d
x
.Taking into account boundary conditions (6)-(7) we rewrite equality (13) in the following form:(14)
1
2
d
d
t
∫
-
a
a
u
,
u
+
w
,
w
d
x
+
B
u
-
,
u
-
x
=
a
+
B
u
-
,
u
-
x
=
-
a
-
1
β
A
w
+
-
B
u
+
,
u
-
x
=
-
a
+
1
β
A
w
+
+
B
u
+
,
u
-
x
=
a
+
F
1
,
u
-
x
=
a
+
F
1
,
u
-
x
=
-
a
=
∫
-
a
a
J
1
u
,
w
,
u
+
J
2
u
,
w
,
w
d
x
,where F
1
=
(
(
1
-
β
)
η
/
α
β
π
)
F.
Let us use the spherical representation [17] of vector U
t
,
x
=
r
t
ω
t
,
x, where ω
t
,
x
=
ω
1
t
,
x
,
ω
2
t
,
x
′, r
t
=
U
t
,
·
L
2
-
a
,
a, and ω
L
2
-
a
,
a
=
1. Furthermore assume that ω
t
,
x
∈
C
2
[
0
,
T
], where the value of T we define below.
Substituting the valuesu
=
r
t
ω
1
t
,
x and w
=
r
t
ω
2
t
,
x into (14) we have that(15)
d
r
d
t
+
r
P
t
=
r
2
Q
t
-
f
t
,where(16)
P
t
=
B
ω
1
-
,
ω
1
-
x
=
a
+
B
ω
1
-
,
ω
1
-
x
=
-
a
+
1
β
A
ω
2
+
,
ω
1
-
x
=
a
+
B
ω
1
+
,
ω
1
-
x
=
a
+
B
ω
1
+
,
ω
1
-
x
=
-
a
-
A
ω
2
+
,
ω
1
-
x
=
-
a
,
Q
t
=
∫
-
a
a
J
1
ω
1
,
ω
2
,
ω
1
+
J
2
ω
1
,
ω
2
,
ω
2
d
x
,
f
t
=
F
1
,
ω
1
-
x
=
a
+
F
1
,
ω
1
-
x
=
-
a
.Let us study (15) with an initial condition (17)
r
0
=
U
0
=
U
0
L
2
-
a
,
a
.Let r
1
t be a particular solution of Riccati equation (15). Then the general solution of (15) has the form(18)
r
t
=
r
1
t
+
exp
⁡
∫
0
t
2
Q
τ
r
1
τ
-
P
τ
d
τ
C
-
∫
0
t
Q
τ
exp
⁡
∫
0
τ
2
Q
ξ
r
1
ξ
-
P
ξ
d
ξ
d
τ
-
1
.Hence, taking into consideration initial condition (17) we obtain (19)
r
t
=
r
1
t
+
exp
⁡
∫
0
t
2
Q
τ
r
1
τ
-
P
τ
d
τ
1
U
0
-
r
1
0
-
∫
0
t
Q
τ
exp
⁡
∫
0
τ
2
Q
ξ
r
1
ξ
-
P
ξ
d
ξ
d
τ
-
1
.If R
t
≡
∫
0
t
Q
τ
exp
⁡
⁡
∫
0
τ
2
Q
ξ
r
1
ξ
-
P
ξ
d
ξ
d
τ
≤
0
∀
t, then the second term of r
(
t
) is bounded ∀
t
∈
[
0
,
+
∞
) under the condition that difference U
0
-
r
1
(
0
) is positive. Further we suppose a difference U
0
-
r
1
0 is a positive number. Let R
t
>
0. We denote by T
1 the moment of time at which (20)
1
U
0
-
r
1
0
-
∫
0
T
1
Q
τ
exp
⁡
∫
0
τ
2
Q
ξ
r
1
ξ
-
P
ξ
d
ξ
d
τ
=
0
.Then r
t
-
r
1
(
t
) is bounded ∀
t
∈
[
0
,
T
], where T
<
T
1 and T
1
~
Ο
1
/
U
0
-
r
1
0. Analyzing structure of r
1
(
t
) (see (31)) from equality (20) we obtain Q
t
>
0 for any t
∈
0
,
T.
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](21)
y
′
′
+
f
1
t
y
′
+
f
2
t
y
=
0
,where(22)
y
t
=
exp
⁡
-
∫
0
t
Q
τ
r
τ
d
τ
,
f
1
t
=
Q
t
P
t
-
Q
′
t
Q
t
,
f
2
t
=
Q
t
f
t
.
Solution of this equation will lead to a solution ofr
t
=
-
y
′
/
y
Q
(
t
) of the original Riccati equation. By substitution (23)
y
t
=
z
t
exp
⁡
-
1
2
∫
0
t
f
1
τ
d
τwe write (21) in the following form [18]: (24)
z
′
′
t
+
a
t
z
t
=
0
,where a
t
=
f
2
t
-
(
1
/
2
)
f
1
′
(
t
)
-
(
1
/
4
)
f
1
2
(
t
).
If instead off
1
t and f
2
t we substitute values from (22), then we get(25)
a
t
=
Q
f
t
-
2
Q
2
P
′
-
Q
′
′
Q
+
Q
′
2
+
Q
P
-
Q
′
2
4
Q
2
.
In work [19] it was proven that function (26)
z
t
=
1
+
∑
k
=
1
∞
b
k
tis the solution of (24), where (27)
b
1
t
=
-
∫
0
t
∫
0
ξ
a
τ
d
τ
d
ξ
,
b
k
t
=
-
∫
0
t
∫
0
ξ
a
τ
b
k
-
1
τ
d
τ
d
ξand the series on the right hand side of (26) uniformly converge.
Then(28)
y
=
1
+
∑
k
=
1
∞
b
k
t
exp
⁡
-
1
2
∫
0
t
f
1
τ
d
τ
,
r
1
t
=
-
y
′
y
Q
=
-
z
′
t
-
1
/
2
z
t
f
1
t
z
t
Q
tis a particular solution of Riccati equation (15).
Sincez
′
0
=
0, value of(29)
r
1
0
=
-
z
′
0
-
1
/
2
z
0
f
1
0
z
0
Q
0
=
f
1
0
2
Q
0
.
Now instead off
1
(
0
) we substitute its value P
0
-
Q
′
(
0
)
/
Q
(
0
) and get r
1
0
=
(
P
0
Q
0
-
Q
′
(
0
)
)
/
2
Q
2
(
0
).
The value ofQ
(
t
) with t
=
0 is(30)
Q
0
=
∫
-
a
a
J
1
ω
1
,
ω
2
,
ω
1
+
J
2
ω
1
,
ω
2
,
ω
2
t
=
0
d
x
>
0
,since Q
t
>
0 for any t
∈
0
,
T
.
An expansion form of the function(31)
r
1
t
=
1
Q
t
1
z
t
∫
0
t
a
τ
z
τ
d
τ
+
Q
t
P
t
-
Q
′
t
2
Q
t
.
It is not difficult to prove using immediate substitution that functionr
1
t satisfies Riccati equation (15).
From equality (19) we obtain(32)
r
t
-
r
1
t
≤
C
U
0
-
r
1
0
.
Hence(33)
U
t
,
·
L
2
-
a
,
a
-
r
1
t
≤
C
U
0
-
r
1
0
,
U
L
2
-
a
,
a
-
r
1
C
0
,
T
≤
C
U
0
-
r
1
0
.
We take as interval of the solution of problem (15)–(17) segment 0
;
1
/
(
U
0
-
r
1
0
), since the integrand Q
τ
exp
⁡
∫
0
τ
2
Q
ξ
r
1
ξ
-
P
ξ
d
ξ is bounded. Hence ∀
t
∈
[
0
,
T
] (where T
≈
T
1) a priori estimation (11) takes place.
Now we prove the existence of solution for (4)–(7) using Galerkin method. Let v
l
(
x
)
l
=
1
∞ be a basis in space L
2
[
-
a
,
a
], where the dimension of vector v
l
(
x
) is equal to the dimension of vector U. For each m we define an approximate solution U
m of (4)–(7) as follows: (34)
U
m
=
∑
j
=
1
m
c
j
m
t
v
j
x
,
(35)
∫
-
a
a
∂
U
m
∂
t
+
A
1
∂
U
m
∂
x
,
v
i
x
d
x
=
∫
-
a
a
J
U
m
,
v
i
x
d
x
,
i
=
1
,
m
¯
,
t
∈
0
,
T
,
(36)
U
m
t
=
0
=
U
0
m
x
,
x
∈
R
,
(37)
A
w
m
-
∓
B
u
m
-
x
=
∓
a
=
1
β
A
w
m
+
±
B
u
m
+
x
=
±
a
±
1
-
β
η
α
β
π
F
,where U
0
m is the orthogonal projection in L
2 of the function U
0 on the subspace, spanned by v
1
,
…
,
v
m, (38)
J
U
m
=
J
1
u
m
,
w
m
,
J
2
u
m
,
w
m
′
.
We representv
j
(
x
) in form v
j
x
=
v
j
1
,
v
j
2
′, where (39)
v
j
1
=
v
j
1
,
v
j
2
,
v
j
3
′
,
v
j
2
=
v
j
4
,
v
j
5
,
v
j
6
′
.
The coefficientsc
j
m
(
t
) are determined from the equations(40)
∑
j
=
1
m
d
c
j
m
d
t
∫
-
a
a
v
j
,
v
i
d
x
+
c
j
m
B
v
j
-
1
,
v
i
-
1
x
=
a
+
B
v
j
-
1
,
v
i
-
1
x
=
-
a
+
B
v
i
-
1
,
v
j
-
1
x
=
a
+
B
v
i
-
1
,
v
j
-
1
x
=
-
a
+
1
β
A
v
i
+
2
+
B
v
i
+
1
,
v
j
-
1
x
=
a
-
1
β
A
v
i
+
2
-
B
v
i
+
1
,
v
j
-
1
x
=
-
a
+
1
β
A
v
j
+
2
+
B
v
j
+
1
,
v
i
-
1
x
=
a
-
1
β
A
v
j
+
2
-
B
v
j
+
1
,
v
i
-
1
x
=
-
a
+
F
1
,
v
i
-
1
x
=
a
+
F
1
,
v
i
-
1
x
=
-
a
+
F
1
,
v
j
-
1
x
=
a
+
F
1
,
v
j
-
1
x
=
-
a
-
∫
-
a
a
A
∂
v
i
2
∂
x
,
v
j
1
+
A
′
∂
v
i
1
∂
x
,
v
j
2
d
x
=
∫
-
a
a
J
∑
j
=
1
m
c
j
m
v
j
,
v
i
d
x
,
i
=
1
,
m
¯
,
t
∈
0
,
T
(41)
c
i
m
0
=
d
i
m
,
i
=
1
,
m
¯
,where d
i
m is the ith component of U
0
m.
We multiply (35) by c
i
m
t and sum over i from 1 to m:(42)
∫
-
a
a
∂
U
m
∂
t
+
A
1
∂
U
m
∂
x
,
U
m
d
x
=
∫
-
a
a
J
U
m
,
U
m
d
x
.With help of the above shown arguments we now prove that r
m
(
t
) is bounded in some time interval [
0
,
T
m
], where U
m
t
,
x
=
r
m
(
t
)
ω
m
t
,
x, T
m
≈
Ο
U
0
m
-
r
m
1
(
0
)
-
1, T
m
≥
T
∀
m, and (43)
U
m
L
2
-
a
,
a
-
r
m
1
C
0
,
T
≤
C
2
U
0
L
2
-
a
,
a
-
r
1
0
,where C
2 is constant and independent from m and r
m
1
(
t
) is a particular solution of Riccati equation about r
m
t.
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 U
m
-
r
m
1 of approximate solutions of problem (5)–(7) is uniformly bounded in function space C
0
,
T
;
L
2
[
-
a
,
a
]. Moreover, the homogeneous system of equations τ
E
+
(
1
/
α
)
A
ξ with respect to τ
,
ξ has only a trivial solution. Then it follows from results in [21] that U
m
-
r
m
1
→
U
-
r
1 is week in C
0
,
T
;
L
2
[
-
a
,
a
] and J
(
U
m
,
)
→
J
U is week in C
0
,
T
;
L
2
[
-
a
,
a
] as m
→
∞. 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 u
1
,
w
1 and u
2
,
w
2. We denote them again by u
=
u
1
-
u
2 and w
=
w
1
-
w
2. Then with respect to new values of u and w we obtain the following problem: (44)
∂
u
∂
t
+
A
∂
w
∂
x
=
J
1
u
1
,
w
1
-
J
1
u
2
,
w
2
∂
w
∂
t
+
A
′
∂
u
∂
x
=
J
2
u
1
,
w
1
-
J
2
u
2
,
w
2
,
x
∈
-
a
,
a
(45)
u
t
=
0
=
0
,
w
t
=
0
=
0
,
(46)
A
w
-
+
B
u
-
x
=
-
a
=
1
β
A
w
+
-
B
u
+
x
=
-
a
,
A
w
-
-
B
u
-
x
=
a
=
1
β
A
w
+
+
B
u
+
x
=
a
.
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))(47)
1
2
d
d
t
∫
-
a
a
u
,
u
+
w
,
w
d
x
+
B
u
-
,
u
-
x
=
a
+
B
u
-
,
u
-
x
=
-
a
-
1
β
A
w
+
-
B
u
+
,
u
-
x
=
-
a
+
1
β
A
w
+
+
B
u
+
,
u
-
x
=
a
=
∫
-
a
a
J
1
u
1
,
w
1
-
J
1
u
2
,
w
2
,
u
+
J
2
u
1
,
w
1
-
J
2
u
2
,
w
2
,
w
d
x
.
Let us transform the integrand on the right side of the equality as(48)
J
1
u
1
,
w
1
-
J
1
u
2
,
w
2
,
u
=
σ
2
-
σ
0
2
f
00,1
f
02,1
-
f
01,1
2
3
-
f
00,2
f
02,2
-
f
01,2
2
3
f
02,1
-
f
02,2
=
σ
2
-
σ
0
2
f
00,1
f
02,1
-
f
00,2
f
02,2
-
1
3
f
01,1
2
-
f
01,2
2
f
02,1
-
f
02,2
=
σ
2
-
σ
0
2
f
00,1
-
f
00,2
f
02,1
+
f
00,2
f
02,1
-
f
02,2
-
1
3
f
01,1
-
f
01,2
f
01,1
+
f
01,2
f
02,1
-
f
02,2
;
J
2
u
1
,
w
1
-
J
2
u
2
,
w
2
,
u
=
J
03
u
1
,
w
1
-
J
03
u
2
,
w
2
,
f
03,1
-
f
03,2
+
J
11
u
1
,
w
1
-
J
11
u
2
,
w
2
,
f
11,1
-
f
11,2
=
1
4
σ
3
+
3
σ
1
-
4
σ
0
f
00,1
-
f
00,2
f
03,1
+
f
00,2
f
03,1
-
f
03,2
+
1
4
5
2
σ
1
+
σ
0
-
3
σ
3
f
01,1
-
f
01,2
f
02,1
+
f
01,2
f
02,1
-
f
02,2
f
03,1
-
f
03,2
+
σ
1
-
σ
0
f
00,1
-
f
00,2
f
11,1
+
f
00,2
f
11,1
-
f
11,2
+
0.5
5
2
f
10,1
-
f
10,2
f
01,1
+
f
10,2
f
01,1
-
f
01,2
-
2
15
f
01,1
-
f
01,2
f
02,1
+
f
01,2
f
02,1
-
f
02,2
f
11,1
-
f
11,2
.Once again we use spherical representation of u
=
r
t
w
1
t
,
x and w
=
r
t
w
2
t
,
x, where r
t
=
U
t
,
·
L
2
-
a
,
a. Then concerning r
t we obtain new initial value problem(49)
d
r
d
t
+
r
P
t
=
r
2
Q
1
t
,
(50)
r
0
=
0
,where P
t has the same value as in (15) and(51)
Q
1
t
=
σ
2
-
σ
0
2
∫
-
a
a
ω
00,1
-
ω
00,2
ω
02,1
+
ω
00,2
ω
02,1
-
ω
02,2
-
1
3
ω
01,1
-
ω
01,2
ω
01,1
+
ω
01,2
ω
02,1
-
ω
02,2
d
x
+
1
4
σ
3
+
3
σ
1
-
4
σ
0
∫
-
a
a
ω
00,1
-
ω
00,2
ω
03,1
+
+
ω
00,2
ω
03,1
-
ω
03,2
ω
03,1
-
ω
03,2
d
x
+
1
4
5
2
σ
1
+
σ
0
-
3
σ
3
∫
-
a
a
ω
01,1
-
ω
01,2
ω
02,1
+
ω
01,2
ω
02,1
-
ω
02,2
ω
03,1
-
ω
03,2
d
x
+
σ
1
-
σ
0
∫
-
a
a
ω
00,1
-
ω
00,2
ω
11,1
+
ω
00,2
ω
11,1
-
ω
11,2
+
1
2
5
2
ω
10,1
-
ω
10,2
ω
01,1
+
ω
10,2
ω
01,1
-
ω
01,2
-
2
15
ω
01,1
-
ω
01,2
ω
02,1
+
ω
01,2
ω
02,1
-
ω
02,2
ω
11,1
-
ω
11,2
d
x
.
The general solution of (49) is(52)
r
t
=
exp
⁡
-
∫
0
t
P
τ
d
τ
C
-
∫
0
t
Q
1
τ
exp
⁡
-
∫
0
τ
P
ξ
d
ξ
d
τ
-
1
.Solution of (49), which is satisfying homogeneous condition (50), is trivial; that is, r
t
=
0. Hence, U
C
(
0
,
T
;
L
2
-
a
,
a
)
=
0 and u
1
=
u
2 and w
1
=
w
2.
The theorem is proved.



---

*Source: 5834620-2016-07-10.xml*