1. 引言
近年来,由于非齐次Navier-Stokes方程在模拟许多物理现象中的突出作用,人们对非齐次Navier-Stokes方程进行了广泛的数学研究,并且取得了大量的成果。由于物理重要性、复杂性、丰富的现象和数学挑战,非齐次不可压缩Navier-Stokes方程的初边值问题拥有大量的研究文献。
在本文中,考虑以下粘性系数依赖于密度的不可压缩热传导Navier-Stokes方程的全局强解,区域
上的三维不可压缩Navier-Stokes方程组可由下面的方程组来描述:
(1.1)
其中,初始条件
(1.2)
边界条件:
(1.3)
其中
是单位外法线,
是时间,函数
分别是流体的密度、速度、绝对温度和压力,且
形变张量满足
。
粘性系数
和
是依赖于
的连续可微函数,满足
(1.4)
其中,
是正常数。
早期的关于Navier-Stokes方程解的存在性研究主要分为两类,即弱解的存在性和强解的存在性。1974年,Kazhikov [1]得到了,在初始密度
远离真空,非齐次Navier-Stokes方程在能量空间中至少有一个全局弱解。他还得到了该系统对于三个空间维度的小数据和二维的所有数据都存在全局强解。2003年,为了克服真空存在带来的困难,Choe和Kim [2]提出了以下兼容性条件:
并且研究了强解的局部存在性,后来2013年由Huang-Wang [3]推广到了全局强解。2019年,Danchin和Mucha [4]建立了三维有界域中具有非负密度的强解的全局存在性。2015年,Zhang和Tan [5]得到了强解的全局存在性和唯一性。最近,2021年,Guo和Li [6]将文献[5]中的结果推广到了粘性系数随温度变化的情况。
如果粘性系数
为常数,1974年,Kazhikov [1]证明了强解的存在性,但是并未证得唯一性。1978年,Ladyzenskaya和Solonnikov [7]证明了二维情况下强解的全局存在性,并且得到了三维情况下强解的局部存在性以及当初值
小时,强解的全局存在性。后来,诸如1999年Itoh和Tani [8],1990年Padula [9] [10],1991年,Salvi [11]等也得出了一些相关结果。2017年,Zhong [12]得到了在三维情况下具有非负密度不可压缩热传导Navier-Stokes方程强解的全局存在性。2020年,在此基础上,Zhong [13]将这个结果推广到了三维问题强解的全局存在性和强解大时间行为。
(1.5)
目前为止,只要粘性系数
取决于密度
,大多数结果都集中在二维情况下。正如文献[14]所指出的,密度和速度之间的强相互作用会使密度依赖粘性系数情况的全局理论以及恒定粘性系数情况不能直接应用。全局弱解是由1988年,DiPerna和Lions [15] [16]导出的。后来,1997年,在粘性系数
是
范数中正常数的小扰动的情况下,Desjardins [17]得到了二维情况下具有更高正则性的全局弱解。2015年,Abidi和Zhang [18]将这个二维结果推广到了强解的情况。对于三维情况,2004年,Cho和Kim [19]通过施加一些初始兼容性条件构建了一个独特的局部强解。1996年,Lions [16]建立了在允许真空的初始密度的任何空间维度上非齐次Navier-Stokes方程弱解的全局存在性。之后,2015年,Liang [20]和2018年,Lü-Shi-Zhong [21]分别建立了二维Cauchy问题解的全局存在性。后来,2015年Huang-Wang [22]研究了有界域上强解的全局存在性。2018年,Wang [23]研究了具有一般外力的三维初边值问题,在初始密度足够小的假设下得到了强解的全局存在性。2022年,Zhong [24]得到了具有大初值和允许真空的二维边值问题强解的全局存在性和强解的大时间行为。
受Zhong [12]-[14]启发,本文研究三维粘性系数依赖于密度的不可压缩热传导Navier-Stokes方程的全局强解。
2. 主要定理
定理2.1 当常数
,假设初始数据
满足正则性条件
(2.1)
和相容性条件
(2.2)
其中
且
。则对于问题(1.1),存在时间T和一个唯一强解,使得当
(2.3)
此外,如果
是局部强解
的最大存在时间,则对任意的
,有
或
(2.4)
定理2.2 假设初始数据
当
时,满足(2.1),(2.2)和定理1.1的条件。则存在一个
仅依赖于
的小正常数
,使得如果
(2.5)
则问题(1.1)有一个唯一的全局强解。
标注2.3 为了克服密度相关粘度带来的困难,先验估计由
开始,然后利用
Stokes方程的正则性结果,完成先验估计的证明。其数学意义在于假设粘性系数梯度的范数有界,确保粘性项与
的相互作用不会破坏解的光滑性。
3. 预备知识
3.1. 符号说明
3.2. 引理
引理3.1 假设当
时,有
和
。令
是边值问题的唯一弱解:
(3.1)
其中在
上有
和
,则有以下结果:
(1) 假设
,
和
(3.2)
(2) 假设当
时
,则
且
(3.3)
其中,常数
仅依赖于和
。
4. 定理2.1证明
设
是系统(1.1)~(1.3)强解
的最大存在时间。假设(2.4)不成立,即
(4.1)
其中
。
4.1. 先验估计
首先,由于传输方程(1.1)1和不可压缩条件
,容易获得下面的引理。
引理4.1.1 假设
是问题(1.1)~(1.3)在
上的一个强解,则对于任意的
,有
(4.2)
引理4.1.2 假设
是问题(1)在
上的一个强解,则对于任意的
,有
(4.3)
证明 将(1.1)1和(1.1)3分别乘以
和
,将两式相加,并在
上关于
积分。将所得方程,对时间t分部积分并利用(1.1)1和(1.1)4,易得(4.3)。
证毕。
引理4.1.3 假设
是问题(1)在
上的一个强解,则对于任意的
,有
(4.4)
证明 因为
是连续可微函数,可从(1.1)1中得
(4.5)
将(1.1)2乘以
并将所得方程在
上积分,可得
(4.6)
因为
结合(4.5)和(4.6)可得
(4.7)
利用Holder和Gagliardo-Nirenberg不等式,可得
(4.8)
利用Sobloev’s嵌入定理和(4.3),可得
(4.9)
将(4.8)和(4.9)代入至(4.7)中,并利用Young不等式,可得
(4.10)
应用引理3.1,(4.1)和(4.2)可得
(4.11)
因此,
(4.12)
将(4.12)代入(4.11),并应用Cauchy不等式,可得
(4.13)
因此,结合(4.3)可得
(4.14)
结合(4.3)与(4.14),并利用Gronwall不等式即可得(4.4)。
证毕。
注4.1.1 在条件(4.1)下,对于任意的
,有
(4.15)
引理4.1.4 当常数
,在条件(4.1)下,对于任意的
,有
(4.16)
证明 将(1.1)2关于时间
微分,可得
(4.17)
将(4.17)乘以
并在
上分部积分,利用(1.1)1得
(4.18)
下面分别估计
。
因为
(4.19)
将
代入(4.18),并结合(4.19)可得
(4.20)
由(4.15)和Sobolev不等式可得
(4.21)
因此,由(4.20),(4.21)和Gronwall不等式可得
(4.22)
由(4.15)和(4.22)可得
(4.23)
证毕。
引理4.1.5 当常数
,在条件(4.1)下,对于任意的
,有
(4.24)
证明 由引理3.1,(4.2),(4.23),Sobolev嵌入定理和Holder不等式可得
结合(4.22)可得
(4.25)
对输运方程(1.1)1去空间导数
,有
(4.26)
当
,用
乘以(4.26)可得
对上式分部积分并结合
可得
因此可得
即
结合上式和(4.25),并利用Gronwall不等式可得
(4.27)
由(4.2)和(4.27)可得
(4.28)
证毕。
引理4.1.6 当常数
,在条件(4.1)下,对于任意的
,有
(4.29)
证明 令是
的平均值,由(4.2),(4.3)和Poincare不等式可得
因此
(4.30)
相似地可得
(4.31)
将(1.1)3乘以
并将所得方程在
上积分,可得
(4.32)
由(4.2),(4.23)并利用Holder不等式和Sobolev不等式可得
(4.33)
结合(4.1),(4.5),(4.23),(4.30)和Sobolev不等式可得
(4.34)
因
(4.35)
对(4.35)分部积分并结合和(4.23)可得
(4.36)
将(4.33),(4.34)和(4.35)代入(4.32)可得
(4.37)
由(1.1)3,(1.3)和椭圆方程标准
估计可得
(4.38)
结合(4.23),(4.28),(4.30)和(4.38)可得
(4.39)
将(4.39)代入(4.37)可得
(4.40)
因为
结合上式与(4.22),(4.30),(4.40)并利用Gronwall不等式可得
(4.41)
对(1.1)3关于
求导并结合(1.1)1和(4.5)可得
(4.42)
将(4.42)乘以
并在
上分部积分
(4.43)
下面将分别估计
将
代入(4.43),可得
(4.44)
由(4.39),(4.40),(4.41)可得
(4.45)
结合(4.44)和(4.45)可得
(4.46)
结合(4.22)和(4.46)并利用Gronwall不等式可得
(4.47)
因此,由(4.46)和(4.47)可得
(4.48)
4.2. 定理2.1的证明
下面将利用反证法来证明。假设(2.4)是错的,即(4.1)成立。引理4.1.1~4.1.6中的常数C均与
无关,即引理4.1.1~4.1.6中得到的所有先验估计对于任意
都是一直有界的。因此函数
在
满足初始条件(1.5)。其中,标准论证所得
,而
,意味着
因此,
有
和
满足
。因此,
也满足(2.2),所以取
作为初始数据,可以将局部强解拓展到
之外,这与
的最大性相矛盾。因此完成了定理1.1的证明。
5. 定理2.2的证明
在本节中,定义
。通常,利用
表示
对的依赖,在本节中,所有的正常数C都与时间t无关。
5.1. 先验估计
引理5.1.1 假设是问题(1.1)在
上的唯一局部强解,且满足初始数据,有
(5.1)
和
(5.2)
此外,有
(5.3)
其中,
,
是
的直径。
证明 因为(1.1)1是输运方程,对其应用标准极大值原理,易得(5.1)。
将(3.1)2乘以
并分部积分,再关于t积分可得
(5.4)
由于
(5.5)
并结合(5.4)可得(5.2)。
由(4.2)和Poincaré不等式可得
(5.6)
因此,
(5.7)
易得
(5.8)
在
上积分,即可得(5.3)。
证毕。
引理5.1.2 假设
是问题(1.1)在
上的唯一局部强解,且满足初始数据,假设
(5.9)
则存在一个仅依赖于
和
的正常数
,使得对于
,有
(5.10)
证明 证明类似引理4.1.3。由(4.7),(4.8),(4.9)易得
(5.11)
其中,
仅依赖于
和
。
因此,结合(5.2)和(5.11)可得
(5.12)
将上式与(5.2)结合可得(5.10)。
将(5.12)乘以
并利用(4.4)可得
(5.13)
由(5.2),(5.3),(5.12)和Gronwall不等式可得(5.10)。
证毕。
引理5.1.3 假设
是问题(1.1)~(1.3)在
上的唯一局部强解,且(5.4)成立,则存在一个仅依赖于
和
的正常数
,使得对于
,有
(5.14)
和
(5.15)
证明 利用(4.18)和引理5.1.4中的
估计,可得
(5.16)
将(5.16)乘以
可得
(5.17)
由(5.10)可得
(5.18)
结合(4.15)和Sobolev不等式,可得
(5.19)
由(5.9),(5.10),(5.18)以及Cauchy-Schwarz不等式,可得
(5.20)
对于
,将(5.16)乘以
,可得
(5.21)
引理5.1.4 假设
是问题(1.1)~(1.3)在
上的唯一局部强解,且(5.4)成立,则存在一个仅依
赖于
和
的正常数
,使得
(5.22)
证明 由Sobolev嵌入定理,仅需估计
。令
并取
,其中
。
则由引理4.1.4和Sobolev不等式,可得
(5.23)
因此,当
时,可得
(5.24)
当
时,由(5.10),(5.11)和(5.16)可得
(5.25)
因此,由(5.24)和(5.25)可得
(5.26)
相似地,如果
,可得
(5.27)
如果
,可得
(5.28)
由(5.27)和(5.28)可得
(5.29)
结合(5.26)和(5.29)可得(5.22)。
引理5.1.5 假设
是问题(1.1)~(1.3)在
上的唯一局部强解,且(5.4)成立,则存在一个仅依
赖于
和
的正常数
,如果
(5.30)
则
(5.31)
证明 由标准能量估计可得
(5.32)
其中,
是一个仅依赖于
的正常数。因此由Gronwall不等式和(5.22)可得
(5.33)
令
,由(5.30)和(5.33)可得(5.31)。
5.2. 定理2.2的证明
下面计划将局部强解推广至全局强解。
令
为引理5.1.5中所述常数,并假设初始数据满足(2.1)和(2.2)且
(5.34)
依据([25],定理1.1),存在时间
使得问题(1.1)~(1.3)在
有一个唯一的局部强解
,
且
依赖于
,
,
和
。因为
和
在
空间的连续性,则存在
使得
。
令
则
。由定理5.1.5易得
(5.35)
显然,
。否则,如果
,由引理4.1.4,4.1.5,4.1.6,对于任意的
,存在一个一致常数
的,使得
应用定理1.1的证明,取
作为初始数据可将局部强解推广至
之外。至此完成定理2.2的证明。
致 谢
作者非常感谢审稿人提出的建设性的意见和友好的建议。
基金项目
作者由国家自然科学基金资助项目——青年科学基金项目(管道中的定常不可压缩与可压缩空泡流的适定性理论,项目编号:12001071)资助。