不可压缩微极流方程在临界Besov空间中的局部适定性
Local Well-Posedness of Incompressible Micropolar Equations in Critical Besov Spac-es
DOI: 10.12677/AAM.2023.125257, PDF, HTML, XML, 下载: 135  浏览: 206 
作者: 邓迎伊:浙江师范大学数学科学学院,浙江 金华
关键词: 微极流系统Besov空间局部适定性Micropolar System Besov Space Local Well-Posedness
摘要: 本文主要研究d (d ≥ 2)维不可压缩微极流系统在L2框架下临界Besov空间中的局部适定性。证明了当方程组(1.1)的初值满足时,(1.1)存在一个唯一的局部解。
Abstract: This paper mainly studies the local well-posedness of d (d ≥ 2)-dimensional incompressible mi-cropolar system in the critical Besov space under the L2 framework. It is proved that (1.1) has a unique local solution when the initial data of (1.1) satisfies .
文章引用:邓迎伊. 不可压缩微极流方程在临界Besov空间中的局部适定性[J]. 应用数学进展, 2023, 12(5): 2561-2568. https://doi.org/10.12677/AAM.2023.125257

1. 引言

本文我们考虑了在 d ( d 2 ) 空间中微极流体方程的局部适定性问题:

{ t u ( χ + ν ) Δ u + u u + π 2 χ curl ω = 0 , t ω μ Δ ω + u ω + 4 χ ω κ div ω 2 χ curl u = 0 , div u = 0 , ( u , ω ) | t = 0 ( u 0 , ω 0 ) , (1.1)

其中 u = u ( x , t ) 表示流体的速度, ω ( x , t ) 表示流体的为旋转速度, π = π ( x , t ) 表示压力, χ , ν , μ , κ > 0

常粘性系数。在本文中,为方便起见,令 κ = μ = 1 ν = χ = 1 2

微极流方程于1966年首次由C. A. Eringen [1] 引入,用于模拟微极流体。它可以看作是一种非牛顿流体,具有微旋转效应和微转动惯量。微极流体是具有微观结构的流体,涵盖了许多现象(如由悬浮在粘性介质中的颗粒所组成的流体等),是Navier-Stokes方程的重要推广。

ω = 0 χ = 0 时,系统(1.1)退化为经典的Navier-Stokes方程。Navier-Stokes方程有以下的尺度不变性:

( u λ ( t , x ) , π λ ( t , x ) ) = ( λ u ( λ 2 t , λ x ) , λ 2 π ( λ 2 t , λ x ) ) , (1.2)

而临界空间是指对于尺度变换保持空间范数不变的空间。

对于Navier-Stokes方程适定性的研究已有大量的文献。对于临界空间,Fujita和Kato [2] 得到了

Navier-Stokes方程在临界齐次Sobolev空间 H ˙ 1 2 中的全局适定性结果。Cannone [3] 证明了Navier-Stokes

方程当 p > 3 时初值在 B ˙ p , 3 / p 1 空间中的适定性;Chemin [4] 推广了Cannone的结果,研究了初值在 B ˙ p , q 3 / p 1 ( 1 p < , 1 q ) 中Navier-Stokes 方程的适定性。

对于不可压缩微极流方程,Chen和Miao [5] 证明了三维不可压缩微极系统对于小初值在 B ˙ p , 3 / p 1 ( 1 p < 6 ) 中的整体适定性。Song受到Chemin [4] 的启发,在 [6] 中推广了Chen [5] 的结果,将 p , q 的范围从 1 p < 6 q = 分别扩展到 1 p < 1 q ,得到了更一般的全局适定性。另外,关于具有全粘性和部分粘性的二维不可压缩微极系统的适定性可以参考 [7] [8] 。

对(1.1)1作用Leray投影算子 P ,可得

{ t u Δ u + P ( u u ) curl ω = 0 , t ω Δ ω + u ω + 2 ω div ω curl u = 0 , div u = 0 , ( u , ω ) | t = 0 ( u 0 , ω 0 ) , (1.3)

其中 P I Q Q Δ 1 div 且由 div u = 0 易知 P u = u

显然,这个系统相较于Navier-Stokes方程已不满足尺度不变性。而处理微极流方程的困难主要在于处理方程中的耦合项 curl u , curl ω ,我们利用Qian等人在 [9] 中的技巧,可在 L 2 能量方法中将(1.3)中的 2 ω curl u , curl ω 相抵消。

在本文中,对于给定的Banach空间X,我们记 ( a , b ) X = a X + b X 。对于 q [ 1 , + ] ,记号 L q ( I ; X ) 表示I上取值于X的可测函数的集合,使得 t f ( t ) X 属于 L q ( I ) 。当 I = [ 0 , T ] 时,记 L q ( 0 , T ; X ) L T q ( X )

主要结论如下。

2. 主要定理

定理2.1令 ( u 0 , ω 0 ) B ˙ 2 , 1 d 2 1 ( d ) × B ˙ 2 , 1 d 2 1 ( d ) ,则存在 T > 0 使得系统(1.1)有一个局部且唯一的解

( u , ω ) ( L ˜ T ( B ˙ 2 , 1 d 2 1 ) L T 1 ( B ˙ 2 , 1 d 2 + 1 ) ) × ( L ˜ T ( B ˙ 2 , 1 d 2 1 ) L T 1 ( B ˙ 2 , 1 d 2 + 1 ) )

3. 预备知识

首先回顾齐次Besov空间的定义。

定义3.1 [10] 令 ( p , r ) [ 1 , ] 2 s ,齐次Besov空间是由分布 u S h ( d ) (即 u S ( d ) 且满足

lim j S ˙ j u L = 0 )所组成的,且有

u B ˙ p , r s ( 2 j s Δ ˙ j u L p ) j l r ( ) .

引理3.1 [11] 令 σ T > 0 1 p 1 q 2 q 1 。令u是热方程(3.1)的解

{ t u Δ u = f , u | t = 0 = u 0 . (3.1)

q 1 = ( 1 + 1 / q 1 1 / q 2 ) 1 ,则存在两个依赖于维数d的正常数c和C使得

u L ˜ T q 1 ( B ˙ p , 1 σ + 2 q 1 ) C ( j 2 j σ Δ ˙ j u 0 L p ( 1 e c T 2 2 j q 1 c q 1 ) 1 / q 1 + j 2 j ( σ 2 + 2 / q 2 ) Δ ˙ j f L T q 2 ( L p ) ( 1 e c T 2 2 j q 1 c q 1 ) 1 / q 1 ) . (3.2)

特别地,有以下估计成立

u L ˜ T q 1 ( B ˙ p , 1 σ + 2 q 1 ) C ( u 0 B ˙ p , 1 σ + f L ˜ T q 2 ( B ˙ p , 1 σ 2 + 2 q 2 ) ) . (3.3)

注3.2通过类似于引理3.1的证明过程,对于线性热方程

{ t u Δ u + 2 u = f , u | t = 0 = u 0 . (3.4)

的解也有估计式(3.2)成立.

引理3.3 [12] 令 1 p , q s 1 d q s 2 d min { 1 p , 1 q } s 1 + s 2 > d max { 0 , 1 p + 1 q 1 } 。对于

( f , g ) B ˙ q , 1 s 1 ( d ) × B ˙ p , 1 s 2 ( d ) ,我们有

f g B ˙ p , 1 s 1 + s 2 d q C f B ˙ q , 1 s 1 g B ˙ p , 1 s 2 . (3.5)

引理3.4 [10] 对于 0 < s 1 < s 2 0 θ 1 ,且 1 p , q 1 , q 2 ,有以下插值不等式成立:

u B ˙ p , 1 s C u B ˙ p , 1 s 1 θ u B ˙ p , 1 s 2 1 θ , u L ˜ T q ( B ˙ p , 1 s ) C u L ˜ T q 1 ( B ˙ p , 1 s 1 ) θ u L ˜ T q 2 ( B ˙ p , 1 s 2 ) 1 θ , (3.6)

其中 s = θ s 1 + ( 1 θ ) s 2 1 q = θ q 1 + 1 θ q 2

4. 定理证明

4.1. 构造逼近解列

首先,构造(1.3)的逼近解列 ( u n , ω n )

( u L n , ω L n ) 是以下方程组的解:

{ t u L n Δ u L n curl ω L n = 0 , t ω L n Δ ω L n div ω L n + 2 ω L n curl u L n = 0 , div u L n = 0 , ( u L n , ω L n ) | t = 0 ( u n , ω n ) , (4.1)

其中 u n j n Δ ˙ j u 0 ω n j n Δ ˙ j ω 0

( u n , ω n ) ( u L n , ω L n ) + ( u ¯ n , ω ¯ n ) ( u ¯ 0 , ω ¯ 0 ) = ( 0 , 0 ) ( u ¯ n + 1 , ω ¯ n + 1 ) 是以下线性系统的解:

{ t u ¯ n + 1 Δ u ¯ n + 1 = P ( u n u n ) + curl ω ¯ n + 1 , t ω ¯ n + 1 Δ ω ¯ n + 1 div ω ¯ n + 1 + 2 ω ¯ n + 1 = u n ω n + curl u ¯ n + 1 , div u ¯ n + 1 = 0 , ( u ¯ n + 1 , ω ¯ n + 1 ) | t = 0 ( 0 , 0 ) . (4.2)

4.2. 逼近解列的一致有界性

先对 ( u L n , ω L n ) 进行估计,对(4.1)1作用 Δ ˙ j

t Δ ˙ j u L n Δ Δ ˙ j u L n = curl Δ ˙ j ω L n ,

再对上式两端乘以 Δ ˙ j u L n 并在 d 上积分,利用Young不等式可得

1 2 d d t Δ ˙ j u L n L 2 2 + Δ ˙ j u L n L 2 2 = d curl Δ ˙ j ω L n Δ ˙ j u L n d x = d Δ ˙ j ω L n curl Δ ˙ j u L n d x Δ ˙ j ω L n L 2 2 + 1 4 Δ ˙ j u L n L 2 2 . (4.3)

同理有

1 2 d d t Δ ˙ j ω L n L 2 2 + Δ ˙ j ω L n L 2 2 + div Δ ˙ j ω L n L 2 2 + 2 Δ ˙ j ω L n L 2 2 = d curl Δ ˙ j u L n Δ ˙ j ω L n d x Δ ˙ j ω L n L 2 2 + 1 4 Δ ˙ j u L n L 2 2 . (4.4)

将(4.3)和(4.4)相加,利用Bernstein不等式得到

1 2 d d t ( Δ ˙ j u L n L 2 2 + Δ ˙ j ω L n L 2 2 ) + c 2 2 j ( Δ ˙ j u L n L 2 2 + Δ ˙ j ω L n L 2 2 ) 0 ,

则有

d d t ( Δ ˙ j u L n L 2 + Δ ˙ j ω L n L 2 ) + c 2 2 j ( Δ ˙ j u L n L 2 + Δ ˙ j ω L n L 2 ) 0.

对上式在 [ 0 , T ] 上积分,再乘以 2 j ( d 2 1 ) 并对 j 求和即得

( u L n , ω L n ) L ˜ T ( B ˙ 2 , 1 d 2 1 ) L T 1 ( B ˙ 2 , 1 d 2 + 1 ) C ( u 0 , ω 0 ) B ˙ 2 , 1 d 2 1 . (4.5)

又由引理3.1

u L n L T 1 ( B ˙ 2 , 1 d 2 + 1 ) C j 2 j ( d 2 1 ) Δ ˙ j u 0 L 2 ( 1 e c T 2 2 j ) + C j 2 j ( d 2 1 ) Δ ˙ j curl ω L n L T 1 ( L 2 ) ( 1 e c T 2 2 j ) . (4.6)

由于 div Q ω L n = Δ Q ω L n ,将 ω L n 的方程(4.1)2分解为

{ t P ω L n Δ P ω L n + 2 P ω L n = curl P u L n , t Q ω L n 2 Δ Q ω L n + 2 Q ω L n = 0 , (4.7)

易知 curl P = curl ,则对(4.7)1和(4.7)2分别运用引理3.1、注3.2,再由 P Q 的性质可得

ω L n L T 1 ( B ˙ 2 , 1 d 2 + 1 ) C ( P ω L n L T 1 ( B ˙ 2 , 1 d 2 + 1 ) + Q ω L n L T 1 ( B ˙ 2 , 1 d 2 + 1 ) ) C j 2 j ( d 2 1 ) Δ ˙ j ω 0 L 2 ( 1 e c T 2 2 j ) + C j 2 j ( d 2 1 ) Δ ˙ j curl u L n L T 1 ( L 2 ) ( 1 e c T 2 2 j ) . (4.8)

根据 ( u L n , ω L n ) 的估计(4.5) (4.6)和(4.8),我们可选取一个时间 T > 0 使得对,有以下不等式成立:

T η 3 , ( u L n , ω L n ) L ˜ T ( B ˙ 2 , 1 d 2 1 ) C U 0 , ( u L n , ω L n ) L T 1 ( B ˙ 2 , 1 d 2 + 1 ) C η 2 , (4.9)

其中 η 1 U 0 ( u 0 , ω 0 ) B ˙ 2 , 1 d 2 1

对于 ( u ¯ n + 1 , ω ¯ n + 1 ) ,我们运用归纳法证明以下估计成立

( u ¯ n , ω ¯ n ) L ˜ T ( B ˙ 2 , 1 d 2 1 ) L T 1 ( B ˙ 2 , 1 d 2 + 1 ) η . (4.10)

假设(4.10)成立,则下证(4.10)对于 n + 1 成立。

通过类似于(4.5)式的估计过程可得

( u ¯ n + 1 , ω ¯ n + 1 ) L ˜ T ( B ˙ 2 , 1 d 2 1 ) L T 1 ( B ˙ 2 , 1 d 2 + 1 ) C ( P ( u n u n ) , u n ω n ) L T 1 ( B ˙ 2 , 1 d 2 1 ) . (4.11)

对于上式右端两项,利用引理3.3、引理3.4以及(4.9)和(4.10)有

P ( u n u n ) L T 1 ( B ˙ 2 , 1 d 2 1 ) C u n L T 2 ( B ˙ 2 , 1 d 2 ) 2 C u L n L T 2 ( B ˙ 2 , 1 d 2 ) 2 + C u ¯ n L T 2 ( B ˙ 2 , 1 d 2 ) 2 C u L n L ˜ T ( B ˙ 2 , 1 d 2 1 ) u L n L T 1 ( B ˙ 2 , 1 d 2 + 1 ) + C u ¯ n L ˜ T ( B ˙ 2 , 1 d 2 1 ) u ¯ n L T 1 ( B ˙ 2 , 1 d 2 + 1 ) C ( U 0 + 1 ) η 2 , (4.12)

u n ω n L T 1 ( B ˙ 2 , 1 d 2 1 ) C u n L T 2 ( B ˙ 2 , 1 d 2 ) ω n L T 2 ( B ˙ 2 , 1 d 2 ) C u n L T 2 ( B ˙ 2 , 1 d 2 ) 2 + C ω n L T 2 ( B ˙ 2 , 1 d 2 ) 2 C ( u L n L T 2 ( B ˙ 2 , 1 d 2 ) 2 + u ¯ n L T 2 ( B ˙ 2 , 1 d 2 ) 2 + ω L n L T 2 ( B ˙ 2 , 1 d 2 ) 2 + ω ¯ n L T 2 ( B ˙ 2 , 1 d 2 ) 2 ) C ( u L n L ˜ T ( B ˙ 2 , 1 d 2 1 ) u L n L T 1 ( B ˙ 2 , 1 d 2 + 1 ) + u ¯ n L ˜ T ( B ˙ 2 , 1 d 2 1 ) u ¯ n L T 1 ( B ˙ 2 , 1 d 2 + 1 ) ) + C ( ω L n L ˜ T ( B ˙ 2 , 1 d 2 1 ) ω L n L T 1 ( B ˙ 2 , 1 d 2 + 1 ) + ω ¯ n L ˜ T ( B ˙ 2 , 1 d 2 1 ) ω ¯ n L T 1 ( B ˙ 2 , 1 d 2 + 1 ) ) C ( U 0 + 1 ) η 2 . (4.13)

将(4.12)和(4.13)代入(4.11),并取 η 足够小,使得 C ( U 0 + 1 ) η 1 ,即得

( u ¯ n + 1 , ω ¯ n + 1 ) L ˜ T ( B ˙ 2 , 1 d 2 1 ) L T 1 ( B ˙ 2 , 1 d 2 + 1 ) η . (4.14)

因此,(4.10)对于 n + 1 成立。

4.3. 解的存在性和唯一性

下证:对于所有的 n 1 T * > 0 足够小, ( u n , ω n ) 为空间 L 2 ( 0 , T * ; B ˙ 2 , 1 d 2 ) 中的柯西序列。

δ u n + 1 = u n + 1 u n δ ω n + 1 = ω n + 1 ω n n 1

(4.15)

对任意的 T ] 0 , T * [ ,记 V n + 1 = ( δ u n + 1 , δ ω n + 1 ) L T 2 ( B ˙ 2 , 1 d 2 )

分别将算子 P , Q 作用于方程(4.15)2,对(4.15)2分解后的两个方程以及(4.15)1分别运用引理3.1、注3.2,再相加得到

(4.16)

下面对上式右端项进行估计:

由Bernstein不等式、Holder不等式以及引理3.3,可得

curl ( δ u n + 1 , δ ω n + 1 ) L T 1 ( B ˙ 2 , 1 d 2 1 ) C ( δ u n + 1 , δ ω n + 1 ) L T 1 ( B ˙ 2 , 1 d 2 ) C T 1 2 ( δ u n + 1 , δ ω n + 1 ) L T 2 ( B ˙ 2 , 1 d 2 ) , (4.17)

(4.18)

u n ω n u n 1 ω n 1 L T 1 ( B ˙ 2 , 1 d 2 1 ) C u n ω n u n ω n 1 + u n ω n 1 u n 1 ω n 1 L T 1 ( B ˙ 2 , 1 d 2 1 ) C u n δ ω n δ u n ω n 1 L T 1 ( B ˙ 2 , 1 d 2 1 ) C ( u n 1 L T 2 ( B ˙ 2 , 1 d 2 ) δ ω n L T 2 ( B ˙ 2 , 1 d 2 ) + ω n 1 L T 2 ( B ˙ 2 , 1 d 2 ) δ u n L T 2 ( B ˙ 2 , 1 d 2 ) ) . (4.19)

Δ ˙ n + 1 ( u 0 , ω 0 ) B ˙ 2 , 1 d 2 1 C ( u 0 , ω 0 ) B ˙ 2 , 1 d 2 1 . (4.20)

δ n Δ ˙ n + 1 ( u 0 , ω 0 ) B ˙ 2 , 1 d 2 1 ,我们知道初值 ( u 0 , ω 0 ) B ˙ 2 , 1 d 2 1 ( d ) × B ˙ 2 , 1 d 2 1 ( d ) ,由此可得当 n 时, δ n 0 。通过(4.9)和(4.10)以及引理3.4,易知 ( u n , ω n ) L T 2 ( B ˙ 2 , 1 d 2 ) 是有界的。结合(4.16)~(4.20),我们得到对任

意的 T ] 0 , T * [

V n + 1 C * δ n + C * T 1 2 V n + 1 + C * D ( T ) V n , (4.21)

其中 C * > 0 是一个常数。取 T * 使得 C * T * 1 2 1 4 C * D ( T * ) 1 4 ,则

V n + 1 4 C * 3 δ n + 1 3 V n . (4.22)

这表明序列 { ( u n , ω n ) } L 2 ( 0 , T * ; B ˙ 2 , 1 d 2 ) 中收敛。结合(4.1)和(4.2),可得极限 ( u , ω ) 是系统(1.3)对于初值 ( u 0 , ω 0 ) B ˙ 2 , 1 d 2 1 × B ˙ 2 , 1 d 2 1 的解且

( u , ω ) L ˜ T ( B ˙ 2 , 1 d 2 1 ) L T 1 ( B ˙ 2 , 1 d 2 + 1 ) . (4.23)

类似于解局部存在的证明过程,可以得到解 ( u , ω ) 的唯一性。

5. 总结

微极流方程又称非对称流方程,由于微极流系统中耦合项的存在,使得对该系统的处理存在一定困难,对于微极流系统适定性的研究也有待完善。本文研究了 d ( d 2 ) 维不可压缩微极流方程在 L 2 框架下的临界Besov空间中的局部适定性,主要借鉴 [13] 的方法。首先构造逼近解列,证明逼近解列一致有界后,再通过证明逼近解列在解空间中为柯西序列,得到极限存在,随后证明该极限在分布的意义下为方程(1.3)的解,从而完成了存在性的证明。唯一性可由类似于柯西列的证明过程得到。本文的创新点在于在证明逼近解列一致有界的过程中,利用Qian等人在 [9] 中的技巧,可在使用 L 2 能量方法时将(1.3)中的 2 ω curl u , curl ω 相抵消,相较于 [5] 中运用Green矩阵的方法更为简便,当然本文这种方法只适用于 L 2 框架。

参考文献

[1] Eringen, A.C. (1966) Theory of Micropolar Fluid. Journal of Mathematics and Mechanics, 16, 1-18.
https://doi.org/10.1512/iumj.1967.16.16001
[2] Fujita, H. and Kato, T. (1964) On the Navier-Stokes Initial Value Problem. Archive for Rational Mechanics and Analysis, 16, 269-315.
https://doi.org/10.1007/BF00276188
[3] Cannone, M. (1995) Ondelettes, Paraproduits et Navier-Stokes. Diderot Editeur, Paris.
[4] Chemin, J.Y. (1999) Théorèmesd’unicité pour le système de Navier-Stokes tridimensionnel. Journal d’Analyse Mathématique, 77, 27-50.
https://doi.org/10.1007/BF02791256
[5] Chen, Q. and Miao, C. (2012) Global Well-Posedness for the Micropolar Fluid System in Critical Besov Spaces. Journal of Differential Equations, 252, 2698-2724.
https://doi.org/10.1016/j.jde.2011.09.035
[6] Song, Z. (2021) The Gevrey Analyticity and Decay for the Micropolar System in the Critical Besov Space. Journal of Evolution Equations, 21, 4751-4771.
https://doi.org/10.1007/s00028-021-00731-0
[7] Chen, M. (2013) Global Well-Posedness of the 2D Incompressi-ble Micropolar Fluid Flows with Partial Viscosity and Angular Viscosity. Acta Mathematica Scientia, 33, 929-935.
https://doi.org/10.1016/S0252-9602(13)60051-X
[8] Dong, B., Li, J. and Wu, J. (2017) Global Well-Posedness and Large-Time Decay for the 2D Micropolar Equations. Journal of Differential Equations, 262, 3488-3523.
https://doi.org/10.1016/j.jde.2016.11.029
[9] Qian, C., Chen, H. and Zhang, T. (2022) Global Existence of Weak Solutions for 3D Incompressible Inhomogeneous Asymmetric Fluids. Mathematische Annalen, 1-39.
https://doi.org/10.1007/s00208-022-02427-3
[10] Bahouri, H., Chemin, J.Y. and Danchin, R. (2011) Fourier Analysis and Nonlinear Partial Differential Equations. Grundlehren der Mathematischen Wissenschaften, Vol. 343, Springer.
https://doi.org/10.1007/978-3-642-16830-7
[11] Danchin, R. (2001) Local Theory in Critical Spaces for Compressible Viscous and Heat-Conductive Gases. Commu- nications in Partial Differential Equations, 26, 1183-1233.
https://doi.org/10.1081/PDE-100106132
[12] Danchin, R. and Xu, J. (2017) Optimal Time-Decay Estimates for the Compressible Navier-Stokes Equations in the Critical Lp Framework. Archive for Rational Mechanics and Analysis, 224, 53-90.
https://doi.org/10.1007/s00205-016-1067-y
[13] Danchin, R. (2006) The Inviscid Limit for Density-Dependent Incompressible Fluids. Annales de la Faculté des Sciences de Toulouse: Mathématiques, 15, 637-688.
https://doi.org/10.5802/afst.1133