Convergence to the Self-similar Solutions to the Homogeneous Boltzmann Equation
Yoshinori Morimoto, Tong Yang, Huijiang Zhao

TL;DR
This paper rigorously proves that solutions to the spatially homogeneous Boltzmann equation with infinite initial energy converge to self-similar solutions, extending understanding beyond finite energy cases.
Contribution
It provides a rigorous mathematical justification for convergence to self-similar solutions in the infinite energy setting for the Boltzmann equation without angular cutoff.
Findings
Solutions tend to self-similar profiles over time
Extension of convergence results to infinite energy initial data
Analysis specific to Maxwellian molecule type cross sections
Abstract
The Boltzmann H-theorem implies that the solution to the Boltzmann equation tends to an equilibrium, that is, a Maxwellian when time tends to infinity. This has been proved in varies settings when the initial energy is finite. However, when the initial energy is infinite, the time asymptotic state is no longer described by a Maxwellian, but a self-similar solution obtained by Bobylev-Cercignani. The purpose of this paper is to rigorously justify this for the spatially homogeneous problem with Maxwellian molecule type cross section without angular cutoff.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsGas Dynamics and Kinetic Theory · Advanced Mathematical Modeling in Engineering · Numerical methods in inverse problems
Convergence to the Self-similar Solutions to the Homogeneous Boltzmann Equation
**Yoshinori Morimoto
**Graduate School of Human and Environmental Studies, Kyoto University
Kyoto 606-8501, Japan
E-mail address: [email protected]
**Tong Yang
**Department of Mathematics, City University of Hong Kong
Kowloon Tang, Hong Kong, China
E-mail address: [email protected]
**Huijiang Zhao
**School of Mathematics and Statistics, Wuhan University
Wuhan 430072, China
E-mail address: [email protected]
Abstract
The Boltzmann H-theorem implies that the solution to the Boltzmann equation tends to an equilibrium, that is, a Maxwellian when time tends to infinity. This has been proved in varies settings when the initial energy is finite. However, when the initial energy is infinite, the time asymptotic state is no longer described by a Maxwellian, but a self-similar solution obtained by Bobylev-Cercignani. The purpose of this paper is to rigorously justify this for the spatially homogeneous problem with Maxwellian molecule type cross section without angular cutoff.
AMS subject classifications: 82C40; 76P05.
Keywords: Measure valued solution, infinite energy, self-similar solutions, time asymptotic states.
1 Introduction
Consider the homogeneous Boltzmann equation
[TABLE]
with initial data
[TABLE]
where the non-negative unknown function is the distribution density function of particles with velocity at time . The right hand side of (1.1) is the Boltzmann bilinear collision operator corresponding to the Maxwellian molecule type cross section
[TABLE]
Here for
[TABLE]
from the conservation of momentum and energy,
[TABLE]
The Maxwellian molecule type cross section in (1.3) is a non-negative function depending only on the deviation angle . As usual, is restricted to by replacing by its “symmetrized” version . Moreover, motivated by inverse power laws, throughout this paper, we assume
[TABLE]
for positive constants and .
As in [7, 8, 9, 11, 16], the Cauchy problem (1.1) and (1.2) is considered in the set of probability measures on . For presentation, we first introduce some function spaces defined in the previous literatures. For , denotes the probability density function on such that
[TABLE]
and moreover when , it requires that
[TABLE]
Following [7], a characteristic function is the Fourier transform of with respect to :
[TABLE]
For each , set with and
[TABLE]
Here the distance between two suitable functions and with is defined by
[TABLE]
Then the set endowed with the distance is a complete metric space. It follows from Lemma 3.12 of [7], that for all and the embeddings hold for .
The advantage of considering the Maxwellian molecule cross section is that the Bobylev formula is in a simple form. That is, by taking the Fourier transform (1.5) of the equation (1.1) leads to the following equation for the new unknown :
[TABLE]
where we have used
[TABLE]
Here,
[TABLE]
satisfying
[TABLE]
From now on, we consider the Cauchy problem for (1.6) with initial condition
[TABLE]
For , it is shown in [7, 11, 12] that this Cauchy problem admits a unique global solution for every . Moreover, for any if is not a single Dirac mass, cf. [12, 13].
To study the large time behavior of the solution, it depends on whether the initial energy is finite or not, and in the above setting, it depends on the parameter , cf. [2, 7, 8, 9, 11, 12, 14, 15, 16] and the references cited therein:
- •
When , the initial datum has finite energy so that the solution tends to the Maxwellian defined by the initial datum. This was indeed proved in the early work by Tanaka [15] using probability theory in the weak convergence in probability. And it was also proved later in [9, 14, 16] by using analytic methods about convergence in Toscani metrics. Moreover, if some moment higher than the second order is assumed to be bounded, the convergence in the distance with is shown to be exponentially decay in time, cf. [9];
- •
When , the initial energy is infinite so that the solution will no longer tend to an equilibrium, but to a self-similar solution
[TABLE]
[TABLE]
Here, is any given constant and is a radially symmetric non-negative function satisfying
[TABLE]
The regularity of the self-similar solution in was proved in [12, 13]. However, the convergence to the self-similar solution is not well understood even though there are some works, cf. [6, 7, 8] about pointwise convergence in radially symmetric setting or in weak topology with scaling. In fact, even how to show convergence in distribution sense has been a problem.
The main difficulties in studying the convergence to the self-similar solutions come from the fact that the self-similar solution has infinite energy and it decays to zero exponentially in time except in norm. The purpose of this paper is to show strong convergence holds when under some conditions on the initial perturbation.
For this, we first consider the distance between two solutions. For and , as in [9, 10], set
[TABLE]
where
[TABLE]
is the Kronecker delta and is a smooth radially symmetric function satisfying and for and for .
The first result in this paper on the time asymptotic stability of the solutions is given by
Theorem 1.1**.**
Suppose for . Let and be the corresponding two global solutions of the Cauchy problem (1.6) with initial data and respectively. Assume for some , the initial data satisfy
[TABLE]
[TABLE]
Then there exists some positive constant independent of and such that
[TABLE]
Here, and
[TABLE]
Note that for the convergence to the self-similar solution, one can simply take . Based on this, in order to obtain a convergence in strong topology, such as in the Sobolev norms, we will give a uniform in time estimate on the solution in -norm that is given in
Theorem 1.2**.**
For , assume that satisfies (1.14)-(1.15) and is not a single Dirac mass, . Then for any given positive constant and any , there exists a positive constant independent of such that
[TABLE]
Consequently, there exists a positive constant independent of such that
[TABLE]
holds for any . Since
[TABLE]
(1.19) and (1.20) imply that when
[TABLE]
the convergence rate given in (1.19) is faster than the decay rate of the self-similar solution itself. Hence in this case, the infinite energy solution converges to the self-similar solution exponentially in time.
Remark 1.1**.**
Since as , the condition (1.21) holds when is close to .
For the case with finite energy, the above stability estimates give a better convergence description on the solution obtained in the previous literatures, which extends the exponential convergence result in the Toscani metrics with , cf. [9], to the Sobolev space for any . In fact, we have
Corollary 1.1**.**
Suppose that is not a single Dirac mass and satisfies
[TABLE]
for some positive constant with . Then for any , there exist positive constants independent of such that
[TABLE]
and
[TABLE]
Here is any given positive constant and .
A direct consequence of (1.23) and (1.24) gives
[TABLE]
for some positive constant depending only on and .
Remark 1.2**.**
Two comments on the above two theorems:
- •
By Lemma 2.6, sufficient conditions for the requirements (1.15) and (1.22) are
[TABLE]
and
[TABLE]
respectively.
- •
The convergence rate in Corollary 1.1 is faster than the corresponding rates in Theorem 1.1 and Theorem 1.2.
Before the end of the introduction, we list some notations used throughout the paper. Firstly, , with , and are used for some generic large positive constants and stand for some generic small positive constants. When the dependence needs to be explicitly pointed out, then the notations like are used. For multi-index , . And means that there is a constant such that , and means and .
The rest of this paper will be organized as follows: Some known results concerning the global solvability, stability, and regularity of solutions to the Cauchy problem (1.6) and (1.9) in are recalled in Section 2. Moreover, some properties of the approximations of the initial data in will also be given in this section. And then the proofs of Theorem 1.1, Theorem 1.2, and Corollary 1.1 will be given in the next three sections respectively.
2 Preliminaries
In this section, we wil first recall the global solvability, stability and regularity results on the Cauchy problem (1.6) and (1.9) obtained in [5, 6, 7, 8, 11, 12, 13]. And then we will study the properties of the approximation on the initial data defined in (2.3) for later stability estimates.
For the Cauchy problem (1.6) and (1.9), the following estimates are proved in [7, 8, 11, 12, 13].
Lemma 2.1**.**
For , if , then the Cauchy problem (1.6) and (1.9) admits a unique global classical solution satisfying
[TABLE]
If is another solution with initial data , then
[TABLE]
Furthermore, if is not a single Dirac mass, then for and .
And for self-similar solution constructed in [5, 6], by [12, 13], we have
Lemma 2.2**.**
For and a constant , there exists a radially symmetric function satisfying (1.11) such that
[TABLE]
is a solution of the Cauchy problem (1.1) with initial datum . Moreover, for .
The relation between and was given in [7] and [12] and it can be stated as follows.
Lemma 2.3**.**
It holds that
- (i).
For , if , then ;
- (ii).
For , if , then for any , .
Since the energy of the initial data is infinite, for analysis, we will first approximate it by a cutoff on the large velocity so that the moment of any order is bounded. And then it remains to show that the solution with this kind of approximation has uniform bound independent of the cutoff paremeter. On the other hand, the approximate solution can not be arbitrary because it has to be in the function space .
For and , let be the smooth function defined in the construction of and set , define
[TABLE]
with
[TABLE]
The properties of the approximation function are given in
Lemma 2.4**.**
For , if we choose sufficiently large, then
- (i).
, and for sufficiently large it holds
[TABLE]
Here the positive constant depends only on ;
- (ii).
*For and sufficiently large , with *norm being uniformly bounded, precisely,
[TABLE]
Thus
[TABLE]
and
[TABLE]
Proof.
We first prove (2.6)-(2.8). Since it is straightforward to verify (2.6), (2.7) is a direct consequence of (2.6) and Lemma 2.3. We only prove (2.8) as follows: For this, note that
[TABLE]
Choose sufficiently large, we have
[TABLE]
Thus
[TABLE]
Similarly,
[TABLE]
From (2.9), (2), and the fact that
[TABLE]
we obtain
[TABLE]
On the other hand, implies that . Consequently, for , it holds that
[TABLE]
so that for each , there exists a such that
[TABLE]
holds for any .
Choose so that is bounded by a constant independent of because of (2.7). Then
[TABLE]
provided that for some sufficiently small .
(2.13) together with (2.14) imply that for any and , we have
[TABLE]
And (2.8) follows directly from (2.12) and (2.15).
Now it remains to prove (2.5). Set
[TABLE]
then
[TABLE]
Firstly, from (2.9), (2.11) and the fact
[TABLE]
we have
[TABLE]
For , by noticing
[TABLE]
we have for that
[TABLE]
For , it holds that
[TABLE]
Thus, (2.19), (2.20) and (2) imply that
[TABLE]
and consequently
[TABLE]
Inserting (2) and (2.23) into (2) yields (2.5) and this completes the proof of Lemma 2.4. ∎
Now let
[TABLE]
be the approximation of defined by (1)3. The following lemma gives the convergence of to as .
Lemma 2.5**.**
Assume
[TABLE]
then
[TABLE]
Proof.
Notice that
[TABLE]
We have from
[TABLE]
and (2.9) that for
[TABLE]
This together with (2)-(2.11) imply that
[TABLE]
For , if is sufficiently large, we have from (2.9), (2.11) and the assumption (2.25) that
[TABLE]
Thus
[TABLE]
Finally for , from (2), (2.11), the assumption (2.25) and , the dominated convergence theorem yields
[TABLE]
Inserting (2.28), (2.29) and (2.30) into (2.28) gives (2.26). This completes the proof of the lemma. ∎
In the last lemma of this section, a sufficient condition for (1.15) on used in Theorem 1.1 is given.
Lemma 2.6**.**
Let , it holds that
[TABLE]
Proof.
In fact, by the assumption (1.14), we have
[TABLE]
The Taylor expansion of to the second order implies that
[TABLE]
and the Taylor expansion to to the third order gives
[TABLE]
Thus interpolation yields
[TABLE]
for . With this, (2.31) follows. And this completes the proof of the lemma. ∎
3 Proof of Theorem 1.1
To prove Theorem 1.1, as in [7], we first approximate the cross section by a sequence of bounded cross sections defined by
[TABLE]
Then consider
[TABLE]
Here
[TABLE]
For and , let and be the approximation of and constructed in the previous section. Since , it follows from Lemma 2.1 that the Cauchy problem (3.2)-(3.3) with ( ) admits a unique non-negative global solution () satisfying (). Moreover, for , (2.2), Lemma 2.3 and Lemma 2.4 imply that
[TABLE]
Here and denote the unique non-negative solutions of the Cauchy problem (3.2)-(3.3) with initial data and respectively, and
[TABLE]
Furthermore,
[TABLE]
Consequently, Lemma 2.1 yields
[TABLE]
Noticing that
[TABLE]
where (3.5) has been used, from (2.8), we have
Lemma 3.1**.**
The limit
[TABLE]
holds uniformly, locally with respect to and .
Putting
[TABLE]
with
[TABLE]
we now deduce the equation for . Set
[TABLE]
Since and satisfy
[TABLE]
we have
[TABLE]
Let
[TABLE]
By taking , we have from Lemma 3.1 and Lemma 2.5 that uniformly, locally with respect to and as . To derive the equation for , we firstly study
[TABLE]
In fact, for , we have
Lemma 3.2**.**
It holds that
[TABLE]
uniformly, locally with respect to and . And satisfies
[TABLE]
for any and some positive constant independent of and .
Proof.
Since
[TABLE]
it follows from Lemma 3.1 and Lemma 2.5 that
[TABLE]
uniformly, locally with respect to and . Here
[TABLE]
and
[TABLE]
To estimate the bounds on and , firstly note that for ,
[TABLE]
imply that
[TABLE]
Here we have used the fact that has a uniform upper bound for any .
If , then so that . Hence, as obtained in [9], we have
[TABLE]
and
[TABLE]
Then
[TABLE]
Thus for , it holds that
[TABLE]
For when , we have from the assumption and Lemma 2.1 that
[TABLE]
The above estimates together with and imply that
[TABLE]
for any . Consequently, for ,
[TABLE]
(3.20) together with (3.22) imply that
[TABLE]
With (3.19) and (3.23), let , the estimate (3.15) follows immediately. This completes the proof of the lemma. ∎
Now by letting in (3.11), we get from Lemma 2.5, Lemma 3.1 and Lemma 3.2 that solves
[TABLE]
Here satisfies (3.15). By Lemmas 2.5, 3.1 and 3.2, , , and are continuous functions of and satisfy (3.24) in the sense of distribution. Since , , , and are uniformly bounded, is also uniformly bounded so that is globally Lipschitz continuous with respect to . Hence (3.24) holds almost everywhere. Furthermore, by the continuity of , , and , we have that is a continuous function of and consequently satisfies (3.24) everywhere.
The next lemma is about the upper bound on for some .
Lemma 3.3**.**
If with , then
[TABLE]
Here with
[TABLE]
Proof.
The proof is divided into two steps, the first step is to show that . Indeed, for , we have from (3.24) that
[TABLE]
On the other hand, by letting in (3.7), we have from Lemma 3.1 that for
[TABLE]
(3.29) together with the definition of imply that
[TABLE]
for any . Hence, by (3), Lemma 3.2 and the fact that , we can deduce by using the Gronwall inequality that there exists a positive constant independent of and such that
[TABLE]
holds for . Here is any given positive constant.
Since the positive constant in (3.30) is independent of , we have from (3.30) by letting that
[TABLE]
With (3.31), set
[TABLE]
we can get from (3.24) and the fact that
[TABLE]
A direct consequence of (3) yields
[TABLE]
Since , we can apply the argument used in [9] to have
[TABLE]
so that (3.26) follows. This completes the proof of the lemma. ∎
We now turn to prove Theorem 1.1. Let and be the unique solutions of the Cauchy problem (3.2)-(3.3) with initial data and respectively, then
[TABLE]
solve
[TABLE]
and
[TABLE]
respectively.
The estimate (3.26) in Lemma 3.3 gives
[TABLE]
Putting (3.36) and (3.39) together yields
[TABLE]
Noticing that
[TABLE]
we have
[TABLE]
On the other hand, it is shown in [7, 11] that uniformly as , locally in with respect to . By (3.40), we obtain from (3.41) that
[TABLE]
(3.42) is exactly (1.16) and thus the proof of Theorem 1.1 is completed.
4 Proof of Theorem 1.2
To prove Theorem 1.2, compared with Theorem 1.1, we only need to obtain the uniform estimate (1.18) on and the key point is to deduce the following coercivity estimate.
Lemma 4.1**.**
There exists a sufficiently large positive constant such that
[TABLE]
holds for any and some positive constant which depends only on .
Proof.
Notice that
[TABLE]
we have from Theorem 1.1 of [12] that for that , and consequently Lemma 3 of [1] shows that there exists a positive constant independent of and such that
[TABLE]
Hence, we have
[TABLE]
On the other hand, we have from the stability estimate given in Theorem 1.1 that
[TABLE]
holds for any with a constant independent of and .
A direct consequence of (4.2) and (4.3) is
[TABLE]
Thus
[TABLE]
With (4.5), we now turn to prove (4.1). Firstly, note that . If we choose sufficiently large such that
[TABLE]
then for , and sufficiently small such that
[TABLE]
we have
[TABLE]
Thus for the case when , and satisfies (4.7), one can deduce from the assumption (1.4), the estimates (4.5), (4.6) and (4) that
[TABLE]
Here we have used the fact that for . This completes the proof of the lemma. ∎
With Lemma 4.1, we now deduce the uniform estimate on . Let be the Fourier transform of with respect to . For any , let with and defined as in Theorem 1.1. Multiplying (1.6) by with being the complex conjugate of gives
[TABLE]
We estimate and term by term as follows. Since , it follows firstly from Lemma 4.1 that
[TABLE]
because
[TABLE]
For , if we use the change of variable for the term , the cancelation lemma (Lemma 1 of [1]) implies that
[TABLE]
For , note that
[TABLE]
We estimate and separately.
For , since for and , we have
[TABLE]
and consequently
[TABLE]
Here we have used the fact that .
For , since
[TABLE]
and
[TABLE]
we obtain
[TABLE]
Hence, there exists a constant depending on such that
[TABLE]
because of . (4) together with (4) shows that there exists a such that
[TABLE]
Inserting (4.10), (4) and (4.14) into (4), we have, for another ,
[TABLE]
which gives
[TABLE]
Noting again, by means of (4.15) we see that for any there exists a such that
[TABLE]
This and (1.20) give
[TABLE]
Moreover,
[TABLE]
[TABLE]
Here we have used the fact that
[TABLE]
(4) is exactly (1.19) and the proof of Theorem 1.2 is completed.
5 Proof of Corollary 1.1
We prove Corollary 1.1 in this last section. Firstly of all, note that Theorem 1.1 and Theorem 1.2 hold for . The purpose of Corollary 1.1 is to have a better convergence rate in the case of finite energy.
In fact, compared with Theorem 1.1 and Theorem 1.2, the main difference is that now the initial data is of finite energy and consequently the corresponding global solution of the Cauchy problem (3.2)-(3.3) with also has finite energy, i.e.
[TABLE]
With (5.1), it is straightforward to show that
[TABLE]
Consequently, the term in (3) can be estimated by
[TABLE]
Here, .
Having (5.3), the proof of Corollary 1.1 is the same as the ones for Theorem 1.1 and Theorem 1.2. Thus, we omit the detail for brevity.
Acknowledgment
The research of the first author was supported by Grant-in-Aid for Scientific Research No. 25400160, Japan Society of the Promotion of Science. The research of the second author was supported by the NSFC-RGC Grant, N-CityU102/12. The research of the third author was supported by the Fundamental Research Funds for the Central Universities of China and three grants from the National Natural Science Foundation of China under contracts 10925103, 11271160 and 11261160485.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] R. Alexandre, L. Desvillettes, C. Villani, and B. Wennberg, Entropy dissipation and longrange interactions. Arch. Rational Mech. Anal. 152 (2000), 327-355.
- 2[2] L. Arkeryd, Asymptotic behaviour of the Boltzmann equation with infinite range forces. Comm. Math. Phys. 86 (1982), 475-484.
- 3[3] A. V. Bobylev, The method of the Fourier transform in the theory of the Boltzmann equation for Maxwell molecules. Dokl. Akad. Nauk SSSR 225 (6) (1975), 1041-1044.
- 4[4] A. V. Bobylev, The theory of the nonlinear spatially uniform Boltzmann equation for Maxwell molecules. In Mathematical physics reviews, Vol. 7, volume 7 of Soviet Sci. Rev. Sect. C Math. Phys. Rev. , pages 111 C 233. Harwood Academic Publ., Chur, 1988.
- 5[5] A. V. Bobylev and C. Cercignani, Exact eternal solutions of the Boltzmann equation. J. Statist. Phys. 106 (5-6) (2002), 1019-1038.
- 6[6] A. V. Bobylev and C. Cercignani, Self-similar solutions of the Boltzmann equation and their applications. J. Statist. Phys. 106 (5-6) (2002), 1039-1071.
- 7[7] M. Cannone and G. Karch, Infinite energy solutions to the homogeneous Boltzmann equation. Comm. Pure Appl. Math. 63 (6) (2010), 747-778.
- 8[8] M. Cannone and G. Karch, On self-similar solutions to the homogeneous Boltzmann equation. Kinetic and Related Models 6 (4) (2013), 801-808.
