三维粘性系数依赖于密度的不可压缩热传导Navier-Stokes方程的全局强解
Global Strong Solution for 3D Viscous Incompressible Heat Conducting Navier-Stokes Equations with Density-Dependent Viscosity
摘要: 本文研究了三维粘性系数依赖于密度的非齐次不可压缩热传导Navier-Stokes方程。首先,当粘性系数的梯度的范数满足 μ( ρ ) L ( 0,T; L p ) < 时,存在一个整体强解,此外,如果初始能量适当小,证明了三维粘性非齐次热传导变粘性Navier-Stokes方程整体强解的唯一性。
Abstract: In this paper, we investigate an 3D viscosity incompressible heat conducting Navier-Stokes equations with density-dependent viscosity. First, we obtain that there exists a global strong solution provided the norm of the gradient of viscosity satisfies μ( ρ ) L ( 0,T; L p ) < . Moreover, if energy is suitably small, we show the uniqueness of the global strong solution to the three-dimensional viscous non-homogeneous heat conducting Navier-Stokes equations with variable viscosity.
文章引用:王智辉. 三维粘性系数依赖于密度的不可压缩热传导Navier-Stokes方程的全局强解[J]. 应用数学进展, 2025, 14(4): 826-842. https://doi.org/10.12677/aam.2025.144210

1. 引言

近年来,由于非齐次Navier-Stokes方程在模拟许多物理现象中的突出作用,人们对非齐次Navier-Stokes方程进行了广泛的数学研究,并且取得了大量的成果。由于物理重要性、复杂性、丰富的现象和数学挑战,非齐次不可压缩Navier-Stokes方程的初边值问题拥有大量的研究文献。

在本文中,考虑以下粘性系数依赖于密度的不可压缩热传导Navier-Stokes方程的全局强解,区域 Ω R 3 上的三维不可压缩Navier-Stokes方程组可由下面的方程组来描述:

{ ρ t +div( ρu )=0, ( ρu ) t +div( ρuu )+P=div( 2μ( ρ )D( u ) ), c ν [ ( ρθ ) t +div( ρuθ ) ]div( κ( ρ )θ )=2μ( ρ ) | Du | 2 , div u=0, (1.1)

其中,初始条件

( ρ,u,θ )( 0,x )=( ρ 0 , u 0 , θ 0 )( x ),xΩ, (1.2)

边界条件:

u=0, θ ν =0,onΩ. (1.3)

其中 ν 是单位外法线, t0 是时间,函数 ρ,u,θ,P 分别是流体的密度、速度、绝对温度和压力,且 D(u)

形变张量满足 D( u )= 1 2 ( u+ ( u ) tr )

粘性系数 μ( ρ ) κ( ρ ) 是依赖于 ρ 的连续可微函数,满足

μ ¯ μ( ρ ) μ _ >0, κ ¯ κ( ρ ) κ _ >0, (1.4)

其中, μ _ , μ ¯ , κ _ , κ ¯ 是正常数。

早期的关于Navier-Stokes方程解的存在性研究主要分为两类,即弱解的存在性和强解的存在性。1974年,Kazhikov [1]得到了,在初始密度 ρ 0 远离真空,非齐次Navier-Stokes方程在能量空间中至少有一个全局弱解。他还得到了该系统对于三个空间维度的小数据和二维的所有数据都存在全局强解。2003年,为了克服真空存在带来的困难,Choe和Kim [2]提出了以下兼容性条件:

{ μΔ u 0 + P 0 = ρ 0 g 1 κΔ θ 0 2μ | u 0 | 2 = ρ 0 g 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]证明了二维情况下强解的全局存在性,并且得到了三维情况下强解的局部存在性以及当初值 u 0 小时,强解的全局存在性。后来,诸如1999年Itoh和Tani [8],1990年Padula [9] [10],1991年,Salvi [11]等也得出了一些相关结果。2017年,Zhong [12]得到了在三维情况下具有非负密度不可压缩热传导Navier-Stokes方程强解的全局存在性。2020年,在此基础上,Zhong [13]将这个结果推广到了三维问题强解的全局存在性和强解大时间行为。

div( μ( ρ 0 ) u 0 )+P= ρ 0 g, (1.5)

目前为止,只要粘性系数 μ( ρ ) 取决于密度 ρ ,大多数结果都集中在二维情况下。正如文献[14]所指出的,密度和速度之间的强相互作用会使密度依赖粘性系数情况的全局理论以及恒定粘性系数情况不能直接应用。全局弱解是由1988年,DiPerna和Lions [15] [16]导出的。后来,1997年,在粘性系数 μ( ρ ) L 范数中正常数的小扰动的情况下,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 当常数 q( 3, ) ,假设初始数据 ( ρ 0 , u 0 , θ 0 ) 满足正则性条件

0 ρ 0 W 1,q , u 0 H 0 1 H 2 , θ 0 H n 2 ,div  u 0 =0, (2.1)

和相容性条件

{ div( 2μ( ρ )D( u 0 ) )+ P 0 = ρ 0 g 1 , div( κ( ρ 0 ) θ 0 )+2μ( ρ 0 ) | D( u 0 ) | 2 = ρ 0 g 2 , (2.2)

其中 P 0 H 1 g 1 , g 2 L 2 。则对于问题(1.1),存在时间T和一个唯一强解,使得当 q 0 ( 3, )

{ ρ0,ρC( [ 0,T ]; W 1, q 0 ), ρ t C( [ 0,T ]; L q 0 ), uC( [ 0,T ]; H 0 1 H 2 ) L 2 ( 0,T; W 2, q 0 ), θ0,θC( [ 0,T ]; H n 2 ) L 2 ( 0,T; W 2, q 0 ), ( u t , θ t ) L 2 ( 0,T; H 1 ),( ρ u t , ρ θ t ) L ( 0,T; L 2 ), (2.3)

此外,如果 T * 是局部强解 ( ρ,u,θ ) 的最大存在时间,则对任意的 3<pq ,有 T * =

lim T T μ( ρ ) L ( 0,T; L p ) =. (2.4)

定理2.2 假设初始数据 ( ρ 0 , u 0 , θ 0 ) 0 ρ 0 ρ ¯ 时,满足(2.1),(2.2)和定理1.1的条件。则存在一个

仅依赖于 Ω, ρ 0 u 0 L 2 , ρ 0 L , u 0 L 2 的小正常数 ε 0 ,使得如果

μ( ρ 0 ) L p ε 0 , (2.5)

则问题(1.1)有一个唯一的全局强解。

标注2.3 为了克服密度相关粘度带来的困难,先验估计由 sup 0t< T μ( ρ ) L p =M< 开始,然后利用

Stokes方程的正则性结果,完成先验估计的证明。其数学意义在于假设粘性系数梯度的范数有界,确保粘性项与 ( u )u 的相互作用不会破坏解的光滑性。

3. 预备知识

3.1. 符号说明

{ L p = L p ( Ω ), H k = H k,2 ( Ω ), W k,p = W k,p ( Ω ), H 0 1 ={ u H 1 |u=0onΩ }, H ν 2 ={ u H 2 |uv=0onΩ }.

3.2. 引理

引理3.1 假设当 3<q< 时,有 0ρ ρ ¯ ρ W 1,q 。令 (u,P) 是边值问题的唯一弱解:

div( μ( ρ )u )+P=F,divu=0inΩ, P μ( ρ ) dx =0, (3.1)

其中在 [ 0, ρ ¯ ] 上有 μ C 1 [ 0, ) μ _ μ( ρ ) μ ¯ ,则有以下结果:

(1) 假设 f L 2 ( u,P ) H 2 × H 1

u H 2 + P μ( ρ ) H 1 C f L 2 ( 1+ μ( ρ ) L q ) q q3 . (3.2)

(2) 假设当 r( 3,q ) F L r ,则 ( u,Π ) W 2,r × W 1,r

u W 2,r + P μ( ρ ) W 1,r C F L r ( 1+ μ( ρ ) L q ) qr 2( qr ) . (3.3)

其中,常数 C 仅依赖于 Ω,q,r, μ ¯ μ

4. 定理2.1证明

T 是系统(1.1)~(1.3)强解 ( ρ,u,θ ) 的最大存在时间。假设(2.4)不成立,即

sup 0t< T μ( ρ ) L p =M<, (4.1)

其中 3<pq

4.1. 先验估计

首先,由于传输方程(1.1)1和不可压缩条件 div u=0 ,容易获得下面的引理。

引理4.1.1 假设 ( ρ,u,θ,P ) 是问题(1.1)~(1.3)在 [ 0, T * ) 上的一个强解,则对于任意的 t[ 0, T * ) ,有

ρ( t ) L = ρ 0 L ρ ¯ . (4.2)

引理4.1.2 假设 ( ρ,u,θ,P ) 是问题(1)在 [ 0, T * ) 上的一个强解,则对于任意的 T[ 0, T * ) ,有

sup 0tT ( ρθ L 1 + ρ u L 2 2 )+ μ _ 0 T u L 2 2 dt C. (4.3)

证明 将(1.1)1和(1.1)3分别乘以 u θ ,将两式相加,并在 Ω 上关于 x 积分。将所得方程,对时间t分部积分并利用(1.1)1和(1.1)4,易得(4.3)。

证毕

引理4.1.3 假设 ( ρ,u,θ,P ) 是问题(1)在 [ 0, T * ) 上的一个强解,则对于任意的 T[ 0, T * ) ,有

sup 0tT u L 2 2 + 0 T ( ρ u t L 2 2 + u H 1 )dt C. (4.4)

证明 因为 μ( ρ ) 是连续可微函数,可从(1.1)1中得

[ μ( ρ ) ] t +uμ( ρ )=0. (4.5)

将(1.1)2乘以 u t 并将所得方程在 Ω 上积分,可得

2 μ ( ρ )D( u ) u t dx+ ρ | u t | 2 dx= ρ uu u t dx. (4.6)

因为

2 μ ( ρ )D( u ) u t dx= d dt μ( ρ ) | D( u ) | 2 dx [ μ( ρ ) ] t | D( u ) | 2 dx

结合(4.5)和(4.6)可得

d dt μ ( ρ ) | D( u ) | 2 dx+ ρ | u t | 2 dx= ρ uu u t dx ( uμ( ρ ) ) | D( u ) | 2 dx i=1 2 I k . (4.7)

利用Holder和Gagliardo-Nirenberg不等式,可得

| I 1 | | ρ u t || ρ u || u |dx C u L 6 u L 3 ρ u t L 2 C ρ u t L 2 u L 2 3 2 u H 1 1 2 1 4 ρ u t L 2 2 +C u L 2 3 u H 1 , (4.8)

利用Sobloev’s嵌入定理和(4.3),可得

| I 2 |=| ( uμ( ρ ) ) | D( u ) | 2 dx |C | μ( ρ ) || u || u |dx C μ( ρ ) L 3 u L 6 u L 4 2 C u L 2 3 2 u H 1 3 2 , (4.9)

将(4.8)和(4.9)代入至(4.7)中,并利用Young不等式,可得

d dt μ ( ρ ) | D( u ) | 2 dx+ ρ u t L 2 2 1 2 ρ u t L 2 2 +C u L 2 3 u H 1 +C u L 2 3 2 u H 1 3 2 1 2 ρ u t L 2 2 +C u L 2 6 +C u H 1 2 , (4.10)

应用引理3.1,(4.1)和(4.2)可得

u H 1 C( ρ u t L 2 + ρuu L 2 ) ( 1+ μ( ρ ) L p ) p p3 C ρ u t L 2 +C ρuu L 2 C ρ u t L 2 +C ρ L u L 6 u L 3 C ρ u t L 2 + u L 2 3 2 u H 1 1 2 1 2 u H 1 +C ρ u t L 2 +C u L 2 3 , (4.11)

因此,

u H 1 C ρ u t L 2 +C u L 2 3 . (4.12)

将(4.12)代入(4.11),并应用Cauchy不等式,可得

d dt μ ( ρ ) | D( u ) | 2 dx+ ρ u t L 2 2 3 4 ρ u t L 2 2 + u L 2 6 + u H 1 2 , (4.13)

因此,结合(4.3)可得

d dt μ ( ρ ) | D( u ) | 2 dx+ ρ u t L 2 2 C, (4.14)

结合(4.3)与(4.14),并利用Gronwall不等式即可得(4.4)。

证毕

4.1.1 在条件(4.1)下,对于任意的 T[ 0, T * ) ,有

u H 1 C ρ u t L 2 . (4.15)

引理4.1.4 当常数 q( 3, ) ,在条件(4.1)下,对于任意的 T[ 0, T * ) ,有

sup 0tT u H 2 C. (4.16)

证明 将(1.1)2关于时间 t 微分,可得

ρ u tt +ρu u t div( 2μ( ρ )D( u t ) )= P t + ρ t ( u t +uu )ρ u t u+div( 2 μ t D( u ) ). (4.17)

将(4.17)乘以 u t 并在 Ω 上分部积分,利用(1.1)1

1 2 d dt ρ | u t | 2 dx +2 μ ( ρ )D( u t ) u t dx = div ( ρu ) | u t | 2 dx+ div ( ρu )uu u t dx ρ u t u u t dx 2 μ t D( u ) u t dx k=1 4 M k . (4.18)

下面分别估计 M k

| M 1 |=| ρ u | u t | 2 dx | 2 ρ L 1 2 u L 6 ρ u t L 3 u t L 2 C ρ L 1 2 u L 2 ρ u t L 2 1 2 ρ u t L 2 1 2 u t L 2 μ _ 8 u t L 2 2 +C ρ u t L 2 2 ,

| M 2 |=| ρ u( uu u t )dx | ρ L u L 6 u L 2 u L 6 u t L 6 + ρ L u L 6 2 2 u L 2 u t L 6 + ρ L u L 6 2 u L 6 u t L 2 C u H 2 u t L 2 μ _ 8 u t L 2 2 +C ρ u t L 2 2 ,

| M 3 |=| ρ u t u u t dx | ρ L 1 2 ρ u t L 3 u t L 6 u L 2 C u L 2 ρ u t L 2 1 2 u t L 2 3 2 μ _ 8 u t L 2 2 +C ρ u t L 2 2 ,

| M 4 |=| 2 μ t D( u ) u t dx | C | u || μ( ρ ) || u || u t |dx C u L μ( ρ ) L p u L 2p p2 u t L 2 μ _ 8 u t L 2 2 +C u L 2 u H 2 2 μ _ 8 u t L 2 2 +C u L 2 ρ u t L 2 2 +C u L 2 .

因为

2 μ ( ρ )D( u t ) u t dx=2 μ ( ρ ) | D( u t ) | 2 dx2 μ _ | D( u t ) | 2 dx = μ _ u t L 2 2 , (4.19)

M k ( k=1,2,3,4 ) 代入(4.18),并结合(4.19)可得

d dt ρ u t L 2 2 + μ _ u t L 2 2 C( 1+ u L 2 ) ρ u t L 2 2 +C u L 2 . (4.20)

由(4.15)和Sobolev不等式可得

0 T u L 2 dt C 0 T u L 2 u H 1 dt C 0 T u H 2 2 dt C 0 T ρ u t L 2 2 dt C, (4.21)

因此,由(4.20),(4.21)和Gronwall不等式可得

sup 0tT ρ u t L 2 2 + 0 T u t L 2 2 dt C, (4.22)

由(4.15)和(4.22)可得

sup 0tT u H 2 C sup 0tT ρ u t L 2 +CC. (4.23)

证毕

引理4.1.5 当常数 q( 3,6 ] ,在条件(4.1)下,对于任意的 T[ 0, T * ) ,有

sup 0tT ρ W 1,q C. (4.24)

证明 由引理3.1,(4.2),(4.23),Sobolev嵌入定理和Holder不等式可得

0 T u L dt 0 T C u W 1,r dt 0 T C ( ρ u t L r + ρuu L r ) ( 1+ μ( ρ ) L p ) pr 2( pr ) dt C 0 T ( ρ u t L r + ρuu L r )dt C 0 T ρ u t L 2 6r 2r u t L 2 3r6 2r dt+C 0 T u L 2 3 r u H 1 2r3 r dt,

结合(4.22)可得

0 T u L dt C. (4.25)

对输运方程(1.1)1去空间导数 ,有

t ρ+u 2 ρ+uρ=0. (4.26)

q( 3,6 ] ,用 q | p | q2 ρ 乘以(4.26)可得

d dt | ρ | q dx +q u 2 ρ | ρ | q2 ρdx=q uρ | ρ | q2 ρdx,

对上式分部积分并结合 div u=0 可得

q u 2 ρ | ρ | q2 ρdx= u ( | ρ | q )dx= | ρ | q divudx =0,

因此可得

d dt ρ L q q q | u | | ρ | q dx q u L ρ L q q ,

d dt ρ L q u L ρ L q ,

结合上式和(4.25),并利用Gronwall不等式可得

sup 0tT ρ L q C, (4.27)

由(4.2)和(4.27)可得

sup 0tT ρ W 1,q C. (4.28)

证毕

引理4.1.6 当常数 q( 3,6 ] ,在条件(4.1)下,对于任意的 T[ 0, T * ) ,有

sup 0tT θ H 2 C. (4.29)

证明 θ ¯ 1 | Ω | θ dx θ 的平均值,由(4.2),(4.3)和Poincare不等式可得

| θ ¯ | ρdx | ρ θdx |+| ρ ( θ θ ¯ )dx |C+C θ L 2 ,

因此

θ H 1 C+C θ L 2 . (4.30)

相似地可得

θ t H 1 C ρ θ t L 2 +C θ t L 2 . (4.31)

将(1.1)3乘以 θ t 并将所得方程在 Ω 上积分,可得

1 2 d dt κ ( ρ ) | θ | 2 dx+ c v ρ | θ t | 2 dx = c v ρ ( uθ ) θ t dx+2 μ ( ρ ) | D( u ) | 2 θ t dx+ 1 2 κ ( ρ ) t | θ | 2 dx I 1 + I 2 + I 3 . (4.32)

由(4.2),(4.23)并利用Holder不等式和Sobolev不等式可得

| I 1 |=| c v ρ ( uθ ) θ t dx | c v 2 ρ θ t L 2 2 +C θ L 2 2 . (4.33)

结合(4.1),(4.5),(4.23),(4.30)和Sobolev不等式可得

I 2 2 d dt μ ( ρ ) | D( u ) | 2 θdx+C u L μ( ρ ) L 6 u L 6 2 θ L 2 +C μ ¯ θ L 6 u L 3 u t L 2 2 d dt μ ( ρ ) | D( u ) | 2 θdx+C u H 2 3 θ L 2 +C θ L 2 u H 1 u t L 2 2 d dt μ ( ρ ) | D( u ) | 2 θdx+C u t L 2 2 +C θ L 2 2 +C. (4.34)

κ ( ρ ) t +uκ( ρ )=0, (4.35)

对(4.35)分部积分并结合和(4.23)可得

I 3 = 1 2 u κ( ρ ) | θ | 2 dx= 1 2 κ ( ρ )u | θ | 2 dx C κ( ρ ) L u L θ L 2 2 θ L 2 C θ L 2 2 +C 2 θ L 2 2 . (4.36)

将(4.33),(4.34)和(4.35)代入(4.32)可得

d dt ( κ( ρ ) | θ | 2 4μ( ρ ) | D( u ) | 2 θ )dx + c v ρ θ t L 2 2 c v 2 ρ θ t L 2 2 +C u t L 2 2 +C 2 θ L 2 2 +C θ L 2 2 +C C u t L 2 2 +C 2 θ L 2 2 +C θ L 2 2 +C. (4.37)

由(1.1)3,(1.3)和椭圆方程标准 H 2 估计可得

θ H 2 2 C κ 1 L ( μ( ρ ) | u | 2 L 2 2 + ρ θ t L 2 2 + ρuθ L 2 2 + κ( ρ )θ L 2 2 + θ L 2 2 ) C u L 4 4 +C ρ L ρ θ t L 2 2 +C ρ L u L 2 θ L 2 2 +C κ L ρ L 4 2 θ L 4 2 +C θ L 2 2 C u H 2 4 +C ρ θ t L 2 2 +C θ H 1 2 +C θ L 2 θ H 2 c v 2 ρ θ t L 2 2 +C θ L 2 2 +C+ 1 2 θ H 2 2 , (4.38)

结合(4.23),(4.28),(4.30)和(4.38)可得

θ H 2 2 c v 2 ρ θ t L 2 2 +C θ L 2 2 +C, (4.39)

将(4.39)代入(4.37)可得

d dt ( κ( ρ ) | θ | 2 4μ( ρ ) | D( u ) | 2 θ )dx + c v ρ θ t L 2 2 C u t L 2 2 +C θ L 2 2 +C. (4.40)

因为

4 μ ( ρ ) | D( u ) | 2 θdxC μ ¯ θ L 6 u L 12 5 2 C θ L 2 u H 2 2 k _ 2 θ L 2 2 +C,

结合上式与(4.22),(4.30),(4.40)并利用Gronwall不等式可得

sup 0tT θ H 1 2 + 0 T ρ θ t L 2 2 dt C. (4.41)

对(1.1)3关于 t 求导并结合(1.1)1和(4.5)可得

c v [ ρ θ tt +ρu θ t ]div( κ( ρ ) θ t ) = c v div( ρu )( θ t +uθ ) c v ρ u t θ+2μ( ρ ) ( | D( u ) | 2 ) t 2( uμ( ρ ) ) | D( u ) | 2 +div( κ ( ρ ) t θ ), (4.42)

将(4.42)乘以 θ t 并在 Ω 上分部积分

c v 2 d dt ρ | θ t | 2 dx+ κ ( ρ ) | θ t | 2 dx = c v div ( ρu ) | θ t | 2 dx+ c v div ( ρu )( uθ ) θ t dx c v ρ ( u t θ ) θ t dx +2 μ ( ρ ) ( | D( u ) | 2 ) t θ t dx2 ( uμ( ρ ) ) | D( u ) | 2 θ t dx κ ( ρ ) t θ θ t dx k=1 6 M ¯ k . (4.43)

下面将分别估计 M ¯ k ( k=1,2,,6 )

| M ¯ 1 |=| c v ρ u | θ t | 2 dx |2 c v ρ L 1 2 u L ρ θ t L 2 θ t L 2 C u L 6 1 2 u L 6 1 2 ρ θ t L 2 θ t L 2 κ _ 12 θ t L 2 2 +C ρ θ t L 2 2 ,

M ¯ k ( k=1,2,,6 ) 代入(4.43),可得

c v d dt ρ θ t L 2 2 + κ _ θ L 2 2 C ρ θ t L 2 2 +C θ H 2 2 +C u t L 2 2 +C. (4.44)

由(4.39),(4.40),(4.41)可得

θ H 2 2 C ρ θ t L 2 2 +C, (4.45)

结合(4.44)和(4.45)可得

c v d dt ρ θ t L 2 2 + κ _ θ t L 2 2 C ρ θ t L 2 2 +C u t L 2 2 +C, (4.46)

结合(4.22)和(4.46)并利用Gronwall不等式可得

sup 0tT ρ θ t L 2 2 + 0 T θ t L 2 2 C, (4.47)

因此,由(4.46)和(4.47)可得

sup 0tT θ H 2 2 C sup 0tT ρ θ t L 2 2 +CC. (4.48)

4.2. 定理2.1的证明

下面将利用反证法来证明。假设(2.4)是错的,即(4.1)成立。引理4.1.1~4.1.6中的常数C均与 t< T * 无关,即引理4.1.1~4.1.6中得到的所有先验估计对于任意 t< T * 都是一直有界的。因此函数

( ρ,u,θ )( T * ,x ) lim t T * ( ρ,u,θ )( t,x ),

t= T * 满足初始条件(1.5)。其中,标准论证所得 ρ, u ˙ , θ ˙ C( [ 0,T ]; L 2 ) ,而 f ˙ f t +uf ,意味着

( ρ u ˙ ,ρ θ ˙ )( T * ,x ) lim t T * ( ρ u ˙ ,ρ θ ˙ )( t,x ) L 2 .

因此,

{ div( 2μ( ρ )D( u ) )+ P| t= T * = ρ ( T * ,x ) g 1 ( x ), div( κ( ρ )θ )+ 2μ( ρ ) | D( u ) | 2 | t= T * = ρ ( T * ,x ) g 2 ( x ),

g 1 ( x ){ ρ 1 2 ( T * ,x )( ρ u ˙ )( T * ,x ),       forx{ x|ρ( T * ,x )>0 }, 0,                                          forx{ x|ρ( T * ,x )=0 },

g 2 ( x ){ c ν ρ 1 2 ( T * ,x )( ρ θ ˙ )( T * ,x ),       forx{ x|ρ( T * ,x )>0 }, 0,                                             forx{ x|ρ( T * ,x )=0 },

满足 g 1 , g 2 L 2 。因此, ( ρ,u,θ )( T * ,x ) 也满足(2.2),所以取 ( ρ,u,θ )( T * ,x ) 作为初始数据,可以将局部强解拓展到 T * 之外,这与 T * 的最大性相矛盾。因此完成了定理1.1的证明。

5. 定理2.2的证明

在本节中,定义 C 0 ρ 0 u 0 L 2 2 。通常,利用 C( f ) 表示 C 对的依赖,在本节中,所有的正常数C都与时间t无关。

5.1. 先验估计

引理5.1.1 假设是问题(1.1)在 [ 0,T ] 上的唯一局部强解,且满足初始数据,有

ρ( t ) L = ρ 0 L , (5.1)

ρ u( t ) L 2 2 + μ _ 0 T u L 2 2 dtdt C 0 . (5.2)

此外,有

sup 0tT e σt ρ u( t ) L 2 2 + μ _ e σt 0 T u L 2 2 dt C 0 . (5.3)

其中, σ= μ ρ ¯ d 2 d Ω 的直径。

证明 因为(1.1)1是输运方程,对其应用标准极大值原理,易得(5.1)。

将(3.1)2乘以 u t 并分部积分,再关于t积分可得

ρ u( t ) L 2 2 +4 0 t μ ( ρ ) | D( u ) | 2 dxds= ρ 0 u 0 L 2 2 . (5.4)

由于

2 μ ( ρ ) | D( u ) | 2 dx2 μ _ | D( u ) | 2 dx = μ _ u L 2 2 , (5.5)

并结合(5.4)可得(5.2)。

由(4.2)和Poincaré不等式可得

ρ u L 2 2 ρ L u L 2 2 ρ ¯ d 2 u L 2 2 , (5.6)

因此,

d dt ρ u L 2 2 +σ ρ u L 2 2 + μ _ u L 2 2 0, (5.7)

易得

d dt ( e σt ρ u L 2 2 )+ μ _ e σt u L 2 2 0. (5.8)

[ 0,T ] 上积分,即可得(5.3)。

证毕

引理5.1.2 假设 ( ρ,u,θ ) 是问题(1.1)在 [ 0,T ] 上的唯一局部强解,且满足初始数据,假设

sup 0tT μ( ρ )( t ) L q 1, (5.9)

则存在一个仅依赖于 Ω, C 0 ,q, μ _ , μ ¯ , ρ 0 L u 0 L 2 的正常数 C ,使得对于 i{ 1,2 } ,有

sup 0tT t i u L 2 2 + 0 T t i ρ u t L 2 2 dtC. (5.10)

证明 证明类似引理4.1.3。由(4.7),(4.8),(4.9)易得

d dt μ ( ρ ) | D( u ) | 2 dx+ 1 2 ρ | u t | 2 dxC u L 2 3 u H 1 +C u L 2 3 2 u H 1 , (5.11)

其中, C 仅依赖于 Ω, ρ 0 L , μ _ , μ ¯ q

因此,结合(5.2)和(5.11)可得

d dt μ( ρ ) | D( u ) | 2 dx+ 1 2 ρ u t L 2 2 C 1 ( C 0 ,C ) u L 2 3 + C 1 ( C 0 ,C ) u L 2 6 + C 1 ( C 0 ,C ) u L 2 9 2 , (5.12)

将上式与(5.2)结合可得(5.10)。

将(5.12)乘以 t 2 并利用(4.4)可得

d dt t μ( ρ ) | D( u ) | 2 dx+ 1 2 t ρ u t L 2 2 μ ( ρ ) | D( u ) | 2 dx+ C 1 ( C 0 ,C )t u L 2 3 + C 1 ( C 0 ,C )t u L 2 3 u L 2 3 2 + u L 2 3 C 1 ( C 0 ,C ) μ ¯ u 2 2 , (5.13)

由(5.2),(5.3),(5.12)和Gronwall不等式可得(5.10)。

证毕

引理5.1.3 假设 ( ρ,u,θ ) 是问题(1.1)~(1.3)在 [ 0,T ] 上的唯一局部强解,且(5.4)成立,则存在一个仅依赖于 Ω, C 0 , ρ 0 L , u 0 L 2 , μ _ , μ ¯ q 的正常数 C ,使得对于 i{ 1,2 } ,有

sup 0tT t ρ u t L 2 2 + 0 T t u t L 2 2 dtC, (5.14)

sup 0tT t 2 ρ u t L 2 2 + 0 T t 2 u t L 2 2 dtC. (5.15)

证明 利用(4.18)和引理5.1.4中的 J k 估计,可得

d dt ρ u t L 2 2 + μ _ u t L 2 2 C( 1+ u L 2 ) ρ u t L 2 2 +C u L 2 u L 2 3 , (5.16)

将(5.16)乘以 t 可得

d dt ( t ρ u t L 2 2 )+ μ _ t u t L 2 2 ρ u t L 2 2 +Ct ρ u t L 2 2 +Ct u L ρ u t L 2 2 +Ct u L 2 u L 2 3 . (5.17)

由(5.10)可得

0 T t ρ u t L 2 2 dt ( 0 T ρ u t L 2 2 dt ) 1 2 ( 0 T t 2 ρ u t L 2 2 dt ) 1 2 C. (5.18)

结合(4.15)和Sobolev不等式,可得

0 T u L 2 dt C 0 T u L 2 u H 1 dt C 0 T u H 2 2 dt C 0 T ρ u t L 2 2 dt C. (5.19)

由(5.9),(5.10),(5.18)以及Cauchy-Schwarz不等式,可得

0 T t u L 2 u L 2 3 dt sup 0tT ( t u L 2 2 ) sup 0tT u L 2 0 T u L 2 dtC sup 0tT t u L 2 2 C. (5.20)

对于 i=2 ,将(5.16)乘以 t 2 ,可得

d dt ( t 2 ρ u t L 2 2 )+ μ _ t 2 u t L 2 2 2t ρ u t L 2 2 +C t 2 ρ u t L 2 2 +C t 2 u L ρ u t L 2 2 +C t 2 u L 2 u L 2 3 . (5.21)

引理5.1.4 假设 ( ρ,u,θ ) 是问题(1.1)~(1.3)在 [ 0,T ] 上的唯一局部强解,且(5.4)成立,则存在一个仅依

赖于 Ω, C 0 , ρ 0 L , u 0 L 2 , μ _ , μ ¯ q 的正常数 C ,使得

0 T u L dt C (5.22)

证明 由Sobolev嵌入定理,仅需估计 u W 1,r 。令 Fρ u t +ρuu 并取 r ,其中 3<r<min{ 6,p }

则由引理4.1.4和Sobolev不等式,可得

0 T u W 1,r C 0 T F L r ( 1+ μ( ρ ) L p ) pr 2( pr ) dt C 0 T ( ρ u t L r + ρuu L r )dt C 0 T ρ u t L 2 6r 2r u t L 2 3r6 2r dt+C 0 T u L 2 3 r u H 1 2r3 r dt. (5.23)

因此,当 0<T1 时,可得

1 T ρ u t L 2 6r 2r u t L 2 3r6 2r dt sup 1tT ( t 2 ρ u t L 2 2 ) 6r 4r ( 1 T t 2 u t L 2 2 dt ) 3r6 4r ( 1 T t 4r r+6 dt ) r+6 4r C. (5.24)

T>1 时,由(5.10),(5.11)和(5.16)可得

1 T ρ u t L 2 6r 2r u t L 2 3r6 2r dt sup 1tT ( t 2 ρ u t L 2 2 ) 6r 4r ( 1 T t 2 u t L 2 2 dt ) 3r6 4r ( 1 T t 4r r+6 dt ) r+6 4r C. (5.25)

因此,由(5.24)和(5.25)可得

0 T ρ u t L 2 6r 2r u t L 2 3r6 2r dtC. (5.26)

相似地,如果 0<T1 ,可得

0 T u L 2 3 r u H 1 2r3 r dtC sup 0tT ( t u L 2 2 ) 3 2r ( 0 T t u H 1 2 dt ) 2r3 2r ( 0 T t 4r r+6 dt ) r+6 4r C. (5.27)

如果 T>1 ,可得

1 T u L 2 3 r u H 1 2r3 r dtC 1 T u H 1 2 dt C sup 1tT ( t 2 u H 1 2 ) 3 2r ( 1 T t 2 dt ) 3 2r C. (5.28)

由(5.27)和(5.28)可得

0 T u L 2 3 r u H 1 2r3 r dtC. (5.29)

结合(5.26)和(5.29)可得(5.22)。

引理5.1.5 假设 ( ρ,u,θ ) 是问题(1.1)~(1.3)在 [ 0,T ] 上的唯一局部强解,且(5.4)成立,则存在一个仅依

赖于 Ω, C 0 , ρ 0 L , u 0 L 2 , μ _ , μ ¯ q 的正常数 C ,如果

μ( ρ 0 ) L q ε 0 , (5.30)

sup 0tT μ( ρ )( t ) L q 1 2 . (5.31)

证明 由标准能量估计可得

d dt μ( ρ ) L p C( p ) u L μ( ρ ) L p , (5.32)

其中, C 是一个仅依赖于 p 的正常数。因此由Gronwall不等式和(5.22)可得

sup 0tT μ( ρ ) L p μ( ρ 0 ) L p exp( C( p ) 0 T u L dt )L μ( ρ 0 ) L p . (5.33)

ε 0 1 2L ,由(5.30)和(5.33)可得(5.31)。

5.2. 定理2.2的证明

下面计划将局部强解推广至全局强解。

ε 0 为引理5.1.5中所述常数,并假设初始数据满足(2.1)和(2.2)且

μ( ρ 0 ) L q ε 0 1 2 . (5.34)

依据([25],定理1.1),存在时间 T * >0 使得问题(1.1)~(1.3)在 Ω×( 0, T * ) 有一个唯一的局部强解 ( ρ,u,θ )

T * 依赖于 ρ 0 W 1,q u 0 H 1 g L 2 μ 。因为 μ( ρ 0 ) L q ε 0 1 2 μ( ρ ) L q 空间的连续性,则存在 T 1 ( 0, T * ) 使得 sup 0tT μ( ρ )( t ) L q 1

T * sup{ T|( ρ,u,θ )( 1.1 )Ω×[ 0,T ] },

T 1 * sup{ T|( ρ,u,θ )( 1.1 )Ω×( 0,T ] sup 0tT μ( ρ ) L q 1 }.

T 1 * T 1 >0 。由定理5.1.5易得

T * = T 1 * . (5.35)

显然, T * = 。否则,如果 T * < ,由引理4.1.4,4.1.5,4.1.6,对于任意的 t( 0, T * ) ,存在一个一致常数 C( T * ) 的,使得

ρ( t ) W 1,q + u( t ) H 2 + θ( t ) H 2 C( T * ),

应用定理1.1的证明,取 ( ρ,u,θ )( T * ,x ) 作为初始数据可将局部强解推广至 T * 之外。至此完成定理2.2的证明。

致 谢

作者非常感谢审稿人提出的建设性的意见和友好的建议。

基金项目

作者由国家自然科学基金资助项目——青年科学基金项目(管道中的定常不可压缩与可压缩空泡流的适定性理论,项目编号:12001071)资助。

参考文献

[1] Kazhikov, A.V. (1974) Resolution of Boundary Value Problems for Nonhomogeneous Viscous Fluids. Doklady Natsionalnoi Akademii Nauk Belarusi, 216, 1008-1010.
[2] Jun Choe, H. and Kim, H. (2003) Strong Solutions of the Navier-Stokes Equations for Nonhomogeneous Incompressible Fluids. Communications in Partial Differential Equations, 28, 1183-1201.
https://doi.org/10.1081/pde-120021191
[3] Huang, X. and Wang, Y. (2013) Global Strong Solution to the 2D Nonhomogeneous Incompressible MHD System. Journal of Differential Equations, 254, 511-527.
https://doi.org/10.1016/j.jde.2012.08.029
[4] Danchin, R. and Mucha, P.B. (2018) The Incompressible Navier-Stokes Equations in Vacuum. Communications on Pure and Applied Mathematics, 72, 1351-1385.
https://doi.org/10.1002/cpa.21806
[5] Zhang, X. and Tan, Z. (2015) The Global Wellposedness of the 3D Heat-Conducting Viscous Incompressible Fluids with Bounded Density. Nonlinear Analysis: Real World Applications, 22, 129-147.
https://doi.org/10.1016/j.nonrwa.2014.08.001
[6] Guo, Z. and Li, Q. (2021) Global Existence and Large Time Behaviors of the Solutions to the Full Incompressible Navier-Stokes Equations with Temperature-Dependent Coefficients. Journal of Differential Equations, 274, 876-923.
https://doi.org/10.1016/j.jde.2020.10.031
[7] Ladyzhenskaya, O.A. and Solonnikov, V.A. (1978) Unique Solvability of an Initial and Boundary-Value Problem for Viscous Incompressible Nonhomogeneous Fluids. Journal of Soviet Mathematics, 9, 697-749.
https://doi.org/10.1007/bf01085325
[8] Itoh, S. and Tani, A. (1999) Solvability of Nonstationary Problems for Nonhomogeneous Incompressible Fluids and the Convergence with Vanishing Viscosity. Tokyo Journal of Mathematics, 22, 17-42.
https://doi.org/10.3836/tjm/1270041610
[9] Padula, M. (1982) An Existence Theorem for Non-Homogeneous Incompressible Fluids. Rendiconti del Circolo Matematico di Palermo, 31, 119-124.
https://doi.org/10.1007/bf02849542
[10] Padula, M. (1990) On the Existence and Uniqueness of Non-Homogeneous Motions in Exterior Domains. Mathematische Zeitschrift, 203, 581-604.
https://doi.org/10.1007/bf02570758
[11] Salvi, R. (1991) The Equations of Viscous Incompressible Non-Homogeneous Fluids: On the Existence and Regularity. The Journal of the Australian Mathematical Society, Series B, Applied Mathematics, 33, 94-110.
https://doi.org/10.1017/s0334270000008651
[12] Zhong, X. (2017) Global Strong Solution for 3D Viscous Incompressible Heat Conducting Navier-Stokes Flows with Non-Negative Density. Journal of Differential Equations, 263, 4978-4996.
https://doi.org/10.1016/j.jde.2017.06.004
[13] Zhong, X. (2020) Global Existence and Large Time Behavior of Strong Solutions for 3D Nonhomogeneous Heat Conducting Navier-Stokes Equations. Journal of Mathematical Physics, 61, Article 111503.
https://doi.org/10.1063/5.0012871
[14] Zhong, X. (2018) Global Strong Solution for Viscous Incompressible Heat Conducting Navier-Stokes Flows with Density-Dependent Viscosity. Analysis and Applications, 16, 623-647.
https://doi.org/10.1142/s0219530518500069
[15] DiPerna, R.J. and Lions, P.L. (1988) Equations différentielles ordinaires et équations de transport avec descoefficients irréguliers. In: Séminaire Équations aux dérivées partielles (Polytechnique) dit aussiSéminaire Goulaouic-Schwartz”, 1-9.
[16] Lions, P.L. (1996) Mathematical Topics in Fluid Mechanics: Volume I: Incompressible Models. Oxford University Press.
[17] Desjardins, B. (1997) Regularity Results for Two-Dimensional Flows of Multiphase Viscous Fluids. Archive for Rational Mechanics and Analysis, 137, 135-158.
https://doi.org/10.1007/s002050050025
[18] Abidi, H. and Zhang, P. (2015) On the Global Well-Posedness of 2-D Inhomogeneous Incompressible Navier-Stokes System with Variable Viscous Coefficient. Journal of Differential Equations, 259, 3755-3802.
https://doi.org/10.1016/j.jde.2015.05.002
[19] Cho, Y. and Kim, H. (2004) Unique Solvability for the Density-Dependent Navier-Stokes Equations. Nonlinear Analysis: Theory, Methods & Applications, 59, 465-489.
https://doi.org/10.1016/j.na.2004.07.020
[20] Liang, Z. (2015) Local Strong Solution and Blow-Up Criterion for the 2D Nonhomogeneous Incompressible Fluids. Journal of Differential Equations, 258, 2633-2654.
https://doi.org/10.1016/j.jde.2014.12.015
[21] Lü, B., Shi, X. and Zhong, X. (2018) Global Existence and Large Time Asymptotic Behavior of Strong Solutions to the Cauchy Problem of 2D Density-Dependent Navier-Stokes Equations with Vacuum. Nonlinearity, 31, 2617-2632.
https://doi.org/10.1088/1361-6544/aab31f
[22] Huang, X. and Wang, Y. (2014) Global Strong Solution with Vacuum to the Two Dimensional Density-Dependent Navier-Stokes System. SIAM Journal on Mathematical Analysis, 46, 1771-1788.
https://doi.org/10.1137/120894865
[23] Wang, W., Yu, H. and Zhang, P. (2018) Global Strong Solutions for 3D Viscous Incompressible Heat Conducting Navier-Stokes Flows with the General External Force. Mathematical Methods in the Applied Sciences, 41, 4589-4601.
https://doi.org/10.1002/mma.4915
[24] Zhong, X. (2022) Global Existence and Large Time Behavior of Strong Solutions for Nonhomogeneous Heat Conducting Navier-Stokes Equations with Large Initial Data and Vacuum. Communications in Mathematical Sciences, 20, 1193-1209.
https://doi.org/10.4310/cms.2022.v20.n5.a1
[25] Cho, Y. and Kim, H. (2008) Existence Result for Heat-Conducting Viscous Incompressible Fluids with Vacuum. Journal of the Korean Mathematical Society, 45, 645-681.
https://doi.org/10.4134/jkms.2008.45.3.645