Next Article in Journal
The Case of Small Prime Numbers Versus the Joye–Libert Cryptosystem
Next Article in Special Issue
A Fast Galerkin Approach for Solving the Fractional Rayleigh–Stokes Problem via Sixth-Kind Chebyshev Polynomials
Previous Article in Journal
The Extendability of Cayley Graphs Generated by Transposition Trees
Previous Article in Special Issue
Solitary Wave Solutions for the Higher Dimensional Jimo-Miwa Dynamical Equation via New Mathematical Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Non-Linear Fokker–Planck Equation in Low-Regularity Space

1
School of Mathematics and Statistics, Nanyang Normal University, Nanyang 473061, China
2
School of Mathematics and Statistics, Hubei University of Arts and Science, Xiangyang 441053, China
3
School of Mathematics and Computational Science, Xiangtan University, Xiangtan 411105, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2022, 10(9), 1576; https://doi.org/10.3390/math10091576
Submission received: 16 March 2022 / Revised: 5 April 2022 / Accepted: 7 April 2022 / Published: 7 May 2022

Abstract

:
We construct the global existence and exponential time decay rates of mild solutions to the non-linear Fokker–Planck equation near a global Maxwellians with small-amplitude initial data in the low regularity function space L k 1 L T L v 2 where the regularity assumption on the initial data is weaker.

1. Introduction and Main Results

In this paper, we are concerned with the Cauchy problem to the nonlinear Fokker–Planck equation as follows
t F + v · x F = ρ v · ( v F + v F ) , F ( 0 , x , v ) = F 0 ( x , v ) ,
where the non-negative unknown function F ( t , x , v ) is the spatially periodic density distribution function of particles with position x = ( x 1 , x 2 , x 3 ) T 3 : = [ π , π ] 3 and velocity v = ( v 1 , v 2 , v 3 ) R 3 at time t 0 , and the density ρ ( t , x ) is defined as ρ = R 3 F d v .
In statistical mechanics, nonlinear Fokker–Planck equation is a partial differential equation which describes the Brownian motion of particles. This equation illustrates the evolution of particle probability density function with velocity, time and space position under the influence of resistance or random force. This equation is also widely used in various fields such as plasma physics, astrophysics, nonlinear hydrodynamics, theory of electronic circuitry and laser arrays, population dynamics, human movement sciences and marketing.
The global equilibria for the nonlinear Fokker–Planck Equation (1) is the normalized global Maxwellian
μ = μ ( v ) = ( 2 π ) 3 2 e | v | 2 2 .
Therefore, we can define the perturbation f = f ( t , x , v ) by
F ( t , x , v ) = μ + μ 1 2 f ( t , x , v ) ,
then the Cauchy problem (1) of the nonlinear Fokker–Planck is reformulated as [1]
t f + v · x f = ρ L f , f ( 0 , x , v ) = f 0 ( x , v ) = μ 1 2 ( F 0 ( x , v ) μ ) ,
where the density ρ ( t , x ) and the linear Fokker–Planck operator L are given by
ρ = 1 + R 3 μ 1 2 f d v
and
L f = μ 1 2 v · μ v ( μ 1 2 f ) = Δ v f + 1 4 ( 6 | v | 2 ) f .
Defining the velocity orthogonal projection
P : L 2 ( R v 3 ) Span μ 1 2 , v i μ 1 2 ( 1 i 3 ) ,
then for any given function f ( t , x , v ) L 2 ( R v 3 ) , one has
P f = a ( t , x ) μ 1 2 + b ( t , x ) · v μ 1 2 ,
with
a = R 3 μ 1 2 f d v , b = R 3 v · μ 1 2 f d v .
Therefore, we have the following macro–micro decomposition of solutions f ( t , x , v ) of the nonlinear Fokker–Planck Equation (1) with respect to the given global Maxwellian μ which was introduced in [2]
f ( t , x , v ) = P f ( t , x , v ) + { I P } f ( t , x , v ) ,
where I denotes the identity operator, P f and { I P } f are called the macroscopic and the microscopic component of f ( t , x , v ) , respectively.
Furthermore, multiplying (2) by μ 1 2 and integrating with respect to v over R 3 to obtain
t a + x · b = 0 ,
then integrating the equality with respect to x over T 3 , we get the conservation of mass
T 3 R 3 μ 1 2 f ( t , x , v ) d v d x = T 3 R 3 μ 1 2 f ( 0 , x , v ) d v d x = 0 .
Let ν ( v ) = 1 + | v | 2 and denote the norm | · | ν by
| f | ν 2 = R v 3 ν ( v ) | f | 2 + | v f | 2 d v .
As is known [3,4,5], the Fokker–Planck operator L is coercive in the sense that there is a positive constant λ 0 such that
( L f , f ) L v 2 λ 0 | { I P } f | ν 2 + | b | 2 .
Motivated by [6], we use the low-regularity function space L k 1 L T L v 2 equipped with norm
f L k 1 L T L v 2 : = Z k 3 sup 0 t T f ^ ( t , k , · ) L v 2 d Σ ( k ) < ,
where the Fourier transformation of f ( t , x , v ) with respect to x T 3 is defined by
f ^ ( t , k , v ) = F x f ( t , k , v ) = T 3 e i k · x f ( t , x , v ) d x , k Z 3 .
In this paper, we will use d Σ ( k ) to denote the discrete measure in Z 3 , i.e.,
Z 3 g ( k ) d Σ ( k ) = k Z 3 g ( k ) .
Notations
  • A B means that there is a constant C > 0 such that A C B . A B means A B and B A .
  • Denoting the dot product ( f , g ) = f · g ¯ for any complex functions.
  • Denoting ( · , · ) L v 2 the complex inner product over L v 2 , i.e.,
    ( f , g ) L v 2 = R v 3 f ( v ) g ( v ) ¯ d v .
  • The convolution of f and g is defined as
    ( f ^ g ^ ) ( k ) = Z l 3 f ^ ( k l ) g ^ ( l ) d Σ ( l ) .
  • R denotes the real part of a complex number.
Based on this preparing work, our main result can be stated as follows.
Theorem 1.
Assume that f 0 ( x , v ) satisfies the conservation of mass
T 3 R 3 μ 1 2 f ( t , x , v ) d v d x = T 3 R 3 μ 1 2 f ( 0 , x , v ) d v d x = 0 .
and F 0 ( x , v ) = μ + μ 1 2 f 0 ( x , v ) 0 , there is a small sufficiently ϵ 0 such that if
f 0 ( x , v ) L k 1 L v 2 ϵ 0 ,
then the Cauchy problem (2) admits a unique global mild solution f ( t , x , v ) satisfying F ( t , x , v ) = μ + μ 1 2 f ( t , x , v ) 0 , and it holds that
f ( t , x , v ) L k 1 L T L v 2 + f ( t , x , v ) L k 1 L T 2 L ν 2 f 0 ( x , v ) L k 1 L v 2 ,
for any t > 0 . Moreover, there is a constant λ > 0 such that the solution also admits the time decay estimate
f ( t , x , v ) L k 1 L v 2 e λ t f 0 ( x , v ) L k 1 L v 2 ,
for any t 0 .
Remark 1.
Compared with the integer Sobolev space H x 4 used in [1], the regularity assumption on the initial data is weaker due to H x 4 L k 1 .
There are a lot of results about the global existence and large time behavior of solutions to Fokker–Planck-type equations, such as, for the Fokker–Planck–Boltzmann equation, the global existence and temporal decay estimates of classical solutions are established based on the nonlinear energy method developed in [2], under Grad’s angular cut-off in [5,7] and without cut-off in [8,9], respectively. As for the Vlasov–Poisson–Fokker–Planck equation, Duan and Liu [3] obtained the time-periodic small-amplitude solution in the three dimensional whole space by Serrin’s method. Hwang and Jang [10], Wang [11] obtained the global existence and the time decay of the solution. For the problem (1), the global existence is proved by combining uniform-in-time energy estimates and the decay rates of the solution is obtained by using the precise spectral analysis of the linearized Fokker–Planck operator as well as the energy method in [1]. Interested readers can refer to the references [4,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28] for more related details.
We note that the previous related results that are obtained in Sobolev space involved the v-derivatives or x-derivatives, which required high regularity of the initial data. In order to obtain the global in time solutions in low-regularity function space, Duan–Liu–Sakamoto–Strain [6] introduced the space L k 1 L T L v 2 to deal with the Landau and non-cutoff Boltzmann equation, where L 1 corresponds to the Weiner algebra over a torus satisfying f g L k 1 f L k 1 g L k 1 . Motivated by this method, we are desired to obtain the global existence of solutions to the Fokker–Planck–Boltzmann equation in low regularity function space.
The rest of this paper is organized as follows. In Section 2, we list some basic lemmas which will be used in the later proof. Section 3 and Section 4 are devoted to deducing global existence and exponential time decay rates for the solution to the Cauchy problem of the Fokker–Planck–Boltzmann equation respectively, where the proofs of Theorem 1 is complete.

2. Basic Lemmas

In this section, we give some results concerning the linear Fokker–Planck operator L and the nonlinear term.
Lemma 1
([3,4,5]). There is a constant λ 0 > 0 such that
R ( L f ^ , f ^ ) L v 2 λ 0 | { I P } f ^ ( t , k , · ) | ν 2 + | b ^ ( t , k ) | 2 .
Lemma 2.
It holds that
Z k 3 0 T a ^ L f ^ , f ^ L v 2 d t 1 / 2 d Σ ( k ) η + f L k 1 L T L v 2 4 η Z k 3 0 T | f ^ ( t , k , · ) | ν 2 d t 1 / 2 d Σ ( k )
for any T > 0 , where the constant η > 0 can be arbitrarily small. It also holds that
Z k 3 0 T | ( | a ^ L f ^ | , μ 1 4 ) L v 2 | 2 d t 1 / 2 d Σ ( k ) f L k 1 L T L v 2 f L k 1 L T 2 L ν 2 .
Proof. 
By Fubini’s theorem, one can obtain
a ^ Δ v f ^ , f ^ L v 2 = R v 3 Z l 3 a ^ ( t , k l ) Δ v f ^ ( t , l , v ) d Σ ( l ) f ^ ( t , k , v ) d v = Z l 3 a ^ ( t , k l ) R v 3 Δ v f ^ ( t , l , v ) f ^ ( t , k , v ) d v d Σ ( l ) = Z l 3 a ^ ( t , k l ) R v 3 v f ^ ( t , l , v ) v f ^ ( t , k , v ) d v d Σ ( l ) Z l 3 | a ^ ( t , k l ) | v f ^ ( t , l , · ) L v 2 v f ^ ( t , k , · ) L v 2 d Σ ( l ) .
By applying Cauchy–Schwarz’s inequality with respect to 0 T ( · ) d t and using Young’s inequality, we have
Z k 3 0 T a ^ Δ v f ^ , f ^ L v 2 d t 1 / 2 d Σ ( k ) Z k 3 0 T Z l 3 | a ^ ( t , k l ) | v f ^ ( t , l , · ) L v 2 v f ^ ( t , k , · ) L v 2 d Σ ( l ) d t 1 / 2 d Σ ( k ) Z k 3 0 T Z l 3 | a ^ ( t , k l ) | v f ^ ( t , l , · ) L v 2 d Σ ( l ) 2 d t 1 / 4 × 0 T v f ^ ( t , k , · ) L v 2 2 d t 1 / 4 d Σ ( k ) η Z k 3 0 T v f ^ ( t , k , · ) L v 2 2 d t 1 / 2 d Σ ( k ) + 1 4 η Z k 3 0 T Z l 3 | a ^ ( t , k l ) | v f ^ ( t , l , · ) L v 2 d Σ ( l ) 2 d t 1 / 2 d Σ ( k ) ,
where η > 0 is a sufficiently small universal constant. For the second term in the above inequality, we can obtain
0 T Z l 3 | a ^ ( t , k l ) | v f ^ ( t , l , · ) L v 2 d Σ ( l ) 2 d t 1 / 2 Z l 3 0 T | a ^ ( t , k l ) | 2 v f ^ ( t , l , · ) L v 2 2 d t 1 / 2 d Σ ( l ) ,
by the Minkowski’s inequality · L l 1 L t 2 · L t 2 L l 1 . By Fubini’s theorem and translation invariance with (16), we obtain
Z k 3 0 T Z l 3 | a ^ ( t , k l ) | v f ^ ( t , l , · ) L v 2 d Σ ( l ) 2 d t 1 / 2 d Σ ( k ) Z k 3 Z l 3 0 T | a ^ ( t , k l ) | 2 v f ^ ( t , l , · ) L v 2 2 d t 1 / 2 d Σ ( l ) d Σ ( k ) Z k 3 Z l 3 sup 0 t T | a ^ ( t , k l ) | 0 T v f ^ ( t , l , · ) L v 2 2 d t 1 / 2 d Σ ( l ) d Σ ( k ) = Z l 3 Z k 3 sup 0 t T | a ^ ( t , k l ) | 0 T v f ^ ( t , l , · ) L v 2 2 d t 1 / 2 d Σ ( k ) d Σ ( l ) = Z l 3 Z k 3 sup 0 t T | a ^ ( t , k l ) | d Σ ( k ) 0 T v f ^ ( t , l , · ) L v 2 2 d t 1 / 2 d Σ ( l ) = Z k 3 sup 0 t T | a ^ ( t , k ) | d Σ ( k ) Z l 3 0 T v f ^ ( t , l , · ) L v 2 2 d t 1 / 2 d Σ ( l ) .
Applying the Hölder inequality, we have
Z k 3 sup 0 t T | a ^ ( t , k ) | d Σ ( k ) = Z k 3 sup 0 t T | R 3 μ 1 2 f ^ ( t , k , v ) d v | d Σ ( k ) Z k 3 sup 0 t T f ^ ( t , k , · ) L v 2 d Σ ( k ) .
Combining (15), (17) and (18), we obtain
Z k 3 0 T a ^ Δ v f ^ , f ^ L v 2 d t 1 / 2 d Σ ( k ) η + f L k 1 L T L v 2 4 η Z k 3 0 T | f ^ ( t , k , · ) | ν 2 d t 1 / 2 d Σ ( k ) ,
with the fact
Z k 3 0 T v f ^ ( t , k , · ) L v 2 2 d t 1 / 2 d Σ ( k ) Z k 3 0 T | f ^ ( t , k , · ) | ν 2 d t 1 / 2 d Σ ( k ) .
Similarly, we can obtain
Z k 3 0 T a ^ 1 4 ( 6 | v | 2 ) f ^ , f ^ L v 2 d t 1 / 2 d Σ ( k ) η + f L k 1 L T L v 2 4 η Z k 3 0 T | f ^ ( t , k , · ) | ν 2 d t 1 / 2 d Σ ( k ) .
By (19) and (20) and the linear Fokker–Planck operator L f ^ = Δ v f ^ + 1 4 ( 6 | v | 2 ) f ^ , the desired result (13) can be obtained. The proof of (14) is also can be deduced similarly. □

3. Global Existence

Firstly, we need to obtain the estimates of the microscopic dissipation for the solution f in (2).
Proposition 1.
Under the assumptions in Theorem 1, it holds that
Z k 3 sup 0 t T f ^ ( t , k , · ) L v 2 d Σ ( k ) + Z k 3 0 T | { I P } f ^ ( t , k , · ) | ν 2 d t 1 / 2 d Σ ( k ) + Z k 3 0 T | b ^ ( t , k ) | 2 d t 1 / 2 d Σ ( k ) f 0 L k 1 L v 2 + η + f L k 1 L T L v 2 4 η Z k 3 0 T | f ^ ( t , k , · ) | ν 2 d t 1 / 2 d Σ ( k ) ,
for any T > 0 , where the constant η > 0 can be arbitrarily small.
Proof. 
Taking Fourier transform of (2) with respect x, we have
t f ^ ( t , k , v ) + i v · f ^ ( t , k , v ) L f ^ ( t , k , v ) = ( a ^ L f ^ ) ( t , k , v ) ,
where the convolutions are taken with respect to k Z 3 :
( a ^ L f ^ ) ( t , k , v ) = Z l 3 a ^ ( t , k l ) L f ^ ( t , l , v ) d Σ ( l ) .
Taking product with the complex conjugate of f ^ ( t , k , v ) and further taking the real part of the resulting equation, we have
1 2 d d t | f ^ ( t , k , v ) | 2 R ( L f ^ , f ^ ) = R ( a ^ L f ^ , f ^ ) ,
integrating the above identity with respect to v and then t, we obtain
1 2 f ^ ( t , k , · ) L v 2 2 0 t R ( L f ^ , f ^ ) L v 2 d τ = 1 2 f ^ 0 ( k , · ) L v 2 2 + 0 t R ( a ^ L f ^ , f ^ ) L v 2 d τ .
Recalling the coercivity estimates of L in Lemma 1, we can obtain
1 2 f ^ ( t , k , · ) L v 2 2 + λ 0 0 t | { I P } f ^ ( τ , k , · ) | ν 2 d τ + λ 0 0 t | b ^ ( τ , k ) | 2 d τ 1 2 f ^ 0 ( k , · ) L v 2 2 + 0 t R ( a ^ L f ^ , f ^ ) L v 2 d τ .
Taking the square root on both sides and using the inequality
1 3 ( A + B + C ) A 2 + B 2 + C 2 , A 2 + B 2 A + B , A , B , C 0 ,
we further have
1 6 f ^ ( t , k , · ) L v 2 + λ 0 3 0 t | { I P } f ^ ( τ , k , · ) | ν 2 d τ 1 / 2 + λ 0 3 0 t | b ^ ( τ , k ) | 2 d τ 1 / 2 1 2 f ^ 0 ( k , · ) L v 2 + 0 t R ( a ^ L f ^ , f ^ ) L v 2 d τ 1 / 2 .
Moreover, taking sup 0 t T on both sides of the above inequality and integrating the resulting inequality with respect to d Σ ( k ) over Z 3 , we have
Z k 3 sup 0 t T f ^ ( t , k , · ) L v 2 d Σ ( k ) + Z k 3 0 T | { I P } f ^ ( t , k , · ) | ν 2 d t 1 / 2 d Σ ( k ) + Z k 3 0 T | b ^ ( t , k ) | 2 d t 1 / 2 d Σ ( k ) f 0 L k 1 L v 2 + Z k 3 0 T a ^ L f ^ , f ^ L v 2 d t 1 / 2 d Σ ( k ) f 0 L k 1 L v 2 + η + f L k 1 L T L v 2 4 η Z k 3 0 T | f ^ ( t , k , · ) | ν 2 d t 1 / 2 d Σ ( k ) ,
by (13) in Lemma 2. Thus the desired estimate (21) are obtained. □
Now we give the estimate of the macroscopic component a ^ ( t , k ) by the dual argument.
Proposition 2.
Under the assumptions of Theorem 1, it holds that
Z k 3 0 T | a ^ ( t , k ) | 2 d t 1 / 2 d Σ ( k ) f L k 1 L T L v 2 + f 0 L k 1 L v 2 + f L k 1 L T L v 2 f L k 1 L T 2 L ν 2 + Z k 3 0 T | b ^ ( t , k ) | 2 d t 1 / 2 d Σ ( k ) + Z k 3 0 T | { I P } f ^ ( t , k , · ) | ν 2 d t 1 / 2 d Σ ( k ) .
Proof. 
In order to obtain the estimate of a, we take a test function as
Φ ^ ( t , k , v ) C 1 ( 0 , ) × Z 3 × R 3 ,
which will be fixed later. Applying the Fourier transform to (2), taking the inner product of it and Φ ^ in L v 2 , then integrating the resultant over [ 0 , T ] , one can obtain
( f ^ , Φ ^ ) | t = T ( f ^ , Φ ^ ) | t = 0 0 T ( f ^ , t Φ ^ ) d t 0 T ( f ^ , v · i k Φ ^ ) d t = 0 T ( a ^ L f ^ , Φ ^ ) d t + 0 T ( L f ^ , Φ ^ ) d t ,
where we have used the notation ( · , · ) = ( · , · ) L v 2 and set T > 0 be an arbitrary fixed constant. Denoting ( f ^ , Φ ^ ) ( T ) = ( f ^ , Φ ^ ) | t = T , ( f ^ , Φ ^ ) ( 0 ) = ( f ^ , Φ ^ ) | t = 0 and plugging in the macro–micro decomposition, we have
0 T ( P f ^ , v · i k Φ ^ ) d t = ( f ^ , Φ ^ ) ( 0 ) ( f ^ , Φ ^ ) ( T ) + 0 T ( f ^ , t Φ ^ ) d t + 0 T ( { I P } f ^ , v · i k Φ ^ ) d t J 1 + 0 T ( a ^ L f ^ , Φ ^ ) d t J 2 + 0 T ( L f ^ , Φ ^ ) d t J 3 .
Thanks to (6), we can get the conservation law for mass, i.e.,
T 3 a ( t , x ) d x = T 3 a ( 0 , x ) d x = 0 ,
so that a ^ ( t , 0 ) = 0 . Now we choose the test function as
Φ ^ ( t , k , v ) = ( | v | 2 10 ) v · i k ϕ ^ a ( t , k ) μ 1 2 ,
where ϕ ^ a ( t , k ) is a solution to
| k | 2 ϕ ^ a ( t , k ) = a ^ ( t , k ) .
Since a ^ ( t , 0 ) = 0 , we can formally write ϕ ^ a ( t , k ) = a ^ ( t , k ) / | k | 2 for any k Z 3 where ϕ ^ a ( t , 0 ) = 0 . Using the estimate of a in [6], we can obtain
0 T ( P f ^ , v · i k Φ ^ ) d t = 5 0 T | a ^ ( t , k ) | 2 d t , | J 1 | f ^ ( T , k , · ) L v 2 2 + f ^ 0 ( k ) L v 2 2 + η 0 T | a ^ ( t , k ) | 2 d t + 0 T | b ^ ( t , k ) | 2 d t + 0 T | { I P } f ^ ( t , k ) | ν 2 d t , | J 2 | η 0 T | a ^ ( t , k ) | 2 d t + 0 T | ( | a ^ L f ^ | , μ 1 4 ) L v 2 | 2 d t ,
where η > 0 is a sufficiently small universal constant. Now we concentrate on the estimate of J 3 , since
| J 3 | = | 0 T ( Δ v f ^ + 1 4 ( 6 | v | 2 ) f ^ , Φ ^ ) d t | | 0 T ( Δ v f ^ , Φ ^ ) d t | J 3 , 1 + | 0 T ( 1 4 ( 6 | v | 2 ) f ^ , Φ ^ ) d t | J 3 , 2 .
By virtue of the macro–micro decomposition, (5) gives
J 3 , 1 | 0 T ( Δ v P f ^ , Φ ^ ) d t | + | 0 T ( Δ v { I P } f ^ , Φ ^ ) d t | | 0 T ( a ^ Δ v μ 1 2 , Φ ^ ) d t | J 3 , 1 1 + | 0 T ( b ^ Δ v ( v μ 1 2 ) , Φ ^ ) d t | J 3 , 1 2 + | 0 T ( Δ v { I P } f ^ , Φ ^ ) d t | J 3 , 1 3 .
Owing to Δ v μ 1 2 = ( 1 4 v 2 3 2 ) μ 1 2 and (26), we obtain
J 3 , 1 1 = | 0 T a ^ Δ v μ 1 2 , ( | v | 2 10 ) v · i k ϕ ^ a ( t , k ) μ 1 2 d t | = | 0 T a ^ ( 1 4 v 2 3 2 ) μ 1 2 , ( | v | 2 10 ) v · i k ϕ ^ a ( t , k ) μ 1 2 d t | = 0 ,
as the integrand function is odd for v. Due to k Z 3 and | k | 2 ϕ ^ a ( t , k ) = a ^ ( t , k ) , it holds that with the Young’s inequality
J 3 , 1 2 0 T | b ^ ( t , k ) | | k | | ϕ ^ a ( t , k ) | d t η 0 T | a ^ ( t , k ) | 2 d t + C η 0 T | b ^ ( t , k ) | 2 d t .
Similarly, we can obtain
J 3 , 1 3 = | 0 T ( { I P } f ^ , Δ v Φ ^ ) d t | d t η 0 T | a ^ ( t , k ) | 2 d t + C η 0 T | { I P } f ^ ( t , k ) | ν 2 d t .
Regarding the estimate of J 3 , 2 , we can also deduce that
J 3 , 2 | 0 T ( 1 4 ( 6 | v | 2 ) P f ^ , Φ ^ ) d t | + | 0 T ( 1 4 ( 6 | v | 2 ) { I P } f ^ , Φ ^ ) d t | η 0 T | a ^ ( t , k ) | 2 d t + C η 0 T | b ^ ( t , k ) | 2 d t + C η 0 T | { I P } f ^ ( t , k ) | ν 2 d t ,
where
0 T ( 1 4 ( 6 | v | 2 ) a ^ μ 1 2 , Φ ^ ) d t = 0 T ( 1 4 ( 6 | v | 2 ) a ^ μ 1 2 , ( | v | 2 10 ) v · i k ϕ ^ a ( t , k ) μ 1 2 ) d t = 0 ,
and
0 T ( 1 4 ( 6 | v | 2 ) ( | v | 2 3 ) c ^ μ 1 2 , Φ ^ ) d t = 0 T ( 1 4 ( 6 | v | 2 ) ( | v | 2 3 ) c ^ μ 1 2 , ( | v | 2 10 ) v · i k ϕ ^ a ( t , k ) μ 1 2 ) d t = 0 ,
since the integrand function is odd for v. By collecting the above estimates and taking η > 0 is small enough, we obtain
0 T | a ^ ( t , k ) | 2 d t f ^ ( T , k , · ) L v 2 2 + f ^ 0 ( k ) L v 2 2 + 0 T | b ^ ( t , k ) | 2 d t + 0 T | { I P } f ^ ( t , k ) | ν 2 d t + 0 T | ( | a ^ L f ^ | , μ 1 4 ) L v 2 | 2 d t .
Thus, we have
Z k 3 0 T | a ^ ( t , k ) | 2 d t 1 / 2 d Σ ( k ) f L k 1 L T L v 2 + f 0 L k 1 L v 2 + f L k 1 L T L v 2 f L k 1 L T 2 L ν 2 + Z k 3 0 T | b ^ ( t , k ) | 2 d t 1 / 2 d Σ ( k ) + Z k 3 0 T | { I P } f ^ ( t , k , · ) | ν 2 d t 1 / 2 d Σ ( k ) ,
where we have used
Z k 3 f ^ ( T , k , · ) L v 2 d Σ ( k ) Z k 3 sup 0 t T f ^ ( t , k , · ) L v 2 d Σ ( k ) = f L k 1 L T L v 2 ,
and
Z k 3 0 T | ( | a ^ L f ^ | , μ 1 4 ) L v 2 | 2 d t 1 / 2 d Σ ( k ) f L k 1 L T L v 2 f L k 1 L T 2 L ν 2 ,
in Lemma 2. Thus the proof the proposition is complete. □
Proposition 3.
Under the assumptions of Theorem 1, it holds that
f L k 1 L T L v 2 + f L k 1 L T 2 L ν 2 f 0 L k 1 L v 2 + f L k 1 L T L v 2 f L k 1 L T 2 L ν 2 .
Proof. 
Taking the linear combination as (21) + M × (25) with M > 0 being small enough and noticing that
P f L k 1 L T 2 L ν 2 [ a , b , c ] L k 1 L T 2 , f L k 1 L T 2 L ν 2 { I P } f L k 1 L T 2 L ν 2 + P f L k 1 L T 2 L ν 2 ,
the desired estimate (27) are obtained. □
With the above preparation in hand, we are ready to deduce the global-in-time existence of the solution. Firstly, the local-in-time existence and uniqueness of the solutions to the Cauchy problem (2) can be established by performing the standard arguments as in [6], where we omit its proof for simplicity. To extend the local solution into the global one, we can deduce that
f L k 1 L T L v 2 + f L k 1 L T 2 L ν 2 f 0 L k 1 L v 2 ,
from (27) in Proposition 3 by virtue of the smallness assumption on f 0 L k 1 L v 2 . Combining this with the local existence, the global mild solution and uniqueness follows immediately from the standard continuity argument. This completes the proof of the global existence and the uniform estimate (10).

4. Large Time Behavior

To deduce the exponential time decay rates of the solution f ( t , x , v ) , we take
h ^ = e λ t f ^
with λ > 0 which will be chosen later. Since f satisfies the system (2), then h ^ satisfies
t h ^ + i k · v h ^ = e λ t ( a ^ L f ^ ) + λ h ^ + L h ^ ,
with initial data
h ^ ( 0 , k , v ) = h ^ 0 ( k , v ) .
Using the same method to deduce Proposition 3, we can obtain
h L k 1 L T L v 2 + h L k 1 L T 2 L ν 2 h 0 L k 1 L v 2 + λ Z k 3 0 T | | h ^ ( t , k , · ) | | L v 2 2 d t 1 / 2 d Σ ( k ) h 0 L k 1 L v 2 + λ h L k 1 L T 2 L ν 2 ,
for any T > 0 . Then, take λ > 0 to be small enough, yields that
h L k 1 L T L v 2 + h L k 1 L T 2 L ν 2 h 0 L k 1 L v 2 .
By the Minkowski’s inequality · L k 1 L T · L T L k 1 , we can deduce
Z k 3 h ^ ( t , k , · ) L v 2 d Σ ( k ) h 0 L k 1 L v 2 .
Since that h ^ = e λ t f ^ , then one can obtain
f ( t ) L k 1 L v 2 e λ t f 0 L k 1 L v 2 ,
Thus, we have completed the proof of Theorem 1.

5. Conclusions

In this paper, the global existence and exponential time decay rates of mild solutions to the nonlinear Fokker–Planck equation are obtained in the low regularity space L k 1 L T L v 2 by the nonlinear energy estimates. Compared with the previous results, the regularity assumption on the initial data is weaker.

Author Contributions

Y.F.: formal analysis, writing—original draft preparation. B.T.: proof of theorem, review and supervision. All authors have read and agreed to the published version of the manuscript.

Funding

The first author is supported by the National Natural Science Foundation of China under contracts 11801285 and Foundation of Nanyang Normal University under contracts 2018ZX014. The corresponding author is supported by the National Natural Science Foundation of China under contracts 12026263.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liao, J.; Wang, Q.R.; Yang, X.F. Global existence and decay rates of the solutions near Maxwellian for non-linear Fokker-Planck equations. J. Stat. Phys. 2018, 173, 222–241. [Google Scholar] [CrossRef]
  2. Guo, Y. The Boltzmann equation in the whole space. Indiana Univ. Math. J. 2004, 53, 1081–1094. [Google Scholar] [CrossRef]
  3. Duan, R.J.; Liu, S.Q. Time-periodic solutions of the Vlasov-Poisson-Fokker-Planck system. Acta Math. Sci. 2015, 35, 876–886. [Google Scholar] [CrossRef]
  4. Duan, R.J.; Fornasier, M.; Toscani, G. A kinetic flocking model with diffusion. Comm. Math. Phys. 2010, 300, 95–145. [Google Scholar] [CrossRef]
  5. Xiong, L.J.; Wang, T.; Wang, L.S. Global existence and decay of solutions to the Fokker-Planck-Boltzmann equation. Kinet Relat. Model. 2014, 1, 169–194. [Google Scholar] [CrossRef]
  6. Duan, R.J.; Liu, S.Q.; Sakamoto, S.; Strain, R.M. Global Mild Solutions of the Landau and Non-Cutoff Boltzmann Equations. Comm. Pure Appl. Math. 2021, 74, 932–1020. [Google Scholar] [CrossRef]
  7. Wang, X.L.; Shi, H.P. Decay and stability of solutions to the Fokker-Planck-Boltzmann equation in. Appl. Anal. 2018, 11, 1933–1959. [Google Scholar] [CrossRef]
  8. Liu, L.Q.; Wang, H. Global existence and decay of solutions for hard potentials to the Fokker-Planck-Boltzmann equation without cut-off. Commun. Pure Appl. Anal. 2020, 19, 3113–3136. [Google Scholar] [CrossRef] [Green Version]
  9. Wang, H. Global existence and decay of solutions for soft potentials to the Fokker-Planck-Boltzmann equation without cut-off. J. Math. Anal. Appl. 2020, 486, 123947. [Google Scholar] [CrossRef]
  10. Hwang, H.J.; Jang, J. On the Vlasov-Poisson-Fokker-Planck equation near Maxwellian. Discret. Contin. Dyn. Syst. Ser. B 2013, 18, 681–691. [Google Scholar] [CrossRef]
  11. Wang, X.L. Global Existence and Long-Time Behavior of Solutions to the Vlasov-Poisson-Fokker-Planck System. Acta Appl. Math. 2020, 170, 853–881. [Google Scholar] [CrossRef]
  12. Barbu, V.; Röckner, M. Nonlinear Fokker–Planck equations with timedependent coefficients. arXiv 2021, arXiv:2110.12460. [Google Scholar]
  13. Barbu, V.; Röckner, M. Solutions for nonlinear Fokker–Planck equations with measures as initial data and McKean-Vlasov equations. J. Funct. Anal. 2021, 280, 108926. [Google Scholar] [CrossRef]
  14. Bouchut, F. Existence and uniqueness of a global smooth solution for the Vlasov-Poisson-Fokker-Planck system in three dimensions. J. Funct. Anal. 1993, 111, 239–258. [Google Scholar] [CrossRef] [Green Version]
  15. Carrillo, J.A.; Soler, J. On the initial value problem for the Vlasov-Poisson-Fokker-Planck system with initial data in Lp spaces. Math. Methods Appl. Sci. 1995, 18, 825–839. [Google Scholar] [CrossRef]
  16. Cercignani, C. The Boltzmann Equation and Its Applications; Springer: New York, NY, USA, 1988. [Google Scholar]
  17. Frank, T.D. Nonlinear Fokker-Planck Equations: Fundamentals and Applications; Springer: New York, NY, USA, 2005. [Google Scholar]
  18. Golse, F.; Imbert, C.; Mouhot, C.; Vasseur, A. Harnack inequality for kinetic Fokker-Planck equations with rough coefficients and application to the Landau equation. arXiv 2016, arXiv:1607.08068. [Google Scholar] [CrossRef]
  19. He, L.B. Regularities of the solutions to the Fokker-Planck-Boltzmann equation. J. Differ. Equ. 2008, 244, 3060–3079. [Google Scholar] [CrossRef] [Green Version]
  20. Hassan, I. Application to differential transformation method for solving systems of differential equations. Appl. Math. Model. 2008, 32, 2552–2559. [Google Scholar] [CrossRef]
  21. Hatami, M.; Ganji, D.D.; Sheikholeslami, M. Differential Transformation Method for Mechanical Engineering Problems; Academic Press: Cambridge, MA, USA, 2016. [Google Scholar]
  22. Ji, M.; Qi, W.; Shen, Z.; Yi, Y. Existence of periodic probability solutions to Fokker-Planck equations with applications. J. Funct. Anal. 2019, 277, 108281. [Google Scholar] [CrossRef]
  23. Kanth, A.; Aruna, K. Differential transform method for solving linear and non-linear systems of partial differential equations. Phys. Lett. A 2008, 372, 6896–6898. [Google Scholar] [CrossRef]
  24. Li, H.L.; Matsumura, A. Behaviour of the Fokker-Planck-Boltzmann equation near a Maxwellian. Arch. Ration. Mech. Anal. 2008, 189, 1–44. [Google Scholar] [CrossRef]
  25. Pulvirenti, M.; Simeoni, C. L estimates for the Vlasov-Poisson-Fokker-Planck equation. Math. Methods Appl. Sci. 2000, 23, 923–935. [Google Scholar] [CrossRef]
  26. Santos, L. Microscopic dynamics of nonlinear Fokker-Planck equations. Phys. Rev. E 2021, 103, 032106. [Google Scholar] [CrossRef] [PubMed]
  27. Villani, C. A review of mathematical topics in collisional kinetic theory. In Handbook of Mathematical Fluid Dynamics; Elsevier: Amsterdam, The Netherlands, 2002; Volume 1, pp. 71–305. [Google Scholar]
  28. Yang, T.; Yu, H.J. Global classical solutions for the Vlasov-Maxwell-Fokker-Planck system. SIAM J. Math. Anal. 2010, 42, 459–488. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fan, Y.; Tang, B. The Non-Linear Fokker–Planck Equation in Low-Regularity Space. Mathematics 2022, 10, 1576. https://doi.org/10.3390/math10091576

AMA Style

Fan Y, Tang B. The Non-Linear Fokker–Planck Equation in Low-Regularity Space. Mathematics. 2022; 10(9):1576. https://doi.org/10.3390/math10091576

Chicago/Turabian Style

Fan, Yingzhe, and Bo Tang. 2022. "The Non-Linear Fokker–Planck Equation in Low-Regularity Space" Mathematics 10, no. 9: 1576. https://doi.org/10.3390/math10091576

APA Style

Fan, Y., & Tang, B. (2022). The Non-Linear Fokker–Planck Equation in Low-Regularity Space. Mathematics, 10(9), 1576. https://doi.org/10.3390/math10091576

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop