可压缩流体Navier-Slip边界条件问题解的存在性研究

王雪娜, 雍 燕

(上海理工大学 理学院,上海 200093)

可压缩流体Navier-Slip边界条件问题解的存在性研究

王雪娜, 雍 燕

(上海理工大学 理学院,上海 200093)

证明了在有界区域Ω⊂3中带Navier-Slip边界条件的可压缩Navier-Stokes方程的解的局部存在性.在证明过程中,首先利用线性化方法将方程转化为线性方程,再利用Galerkin逼近方法得到线性方程组的弱解,通过能量估计方法,得到关于逼近解的一致先验估计,取极限得到方程的解的局部存在性.

可压缩Navier-Stokes方程; Navier-Slip边界条件; 解的存在性

1 问题的提出

研究可压缩Navier-Stokes方程在t>0上的初边值问题的解的存在性已经有了很长的历史,并且划分了三类解,现将逐一介绍这三类解.

第一类解是光滑小解,最早由Matsumura等[1]提出,这种解的初始值和解在H3(Ω)中都近似等于一个常数.后来,Danchin[2-3]将这个解推广到齐次Besov空间.具体做法是将方程组线性化,然后对线性方程组进行迭代,得到一组逼近解,证明逼近解的渐近衰减速率,最后将逼近解取极限得到原方程组解的存在性.这类光滑小解是连续可微的,缺点是不能体现出Navier-Stokes方程组自身的奇性.

第二类解是一般的“大能量”弱解,该解由文献[4-5]引入,并由Feireisl推广.这类解的好处是对任意大的初始能量和非负密度都容易证明解的存在性.但是,这一类解的集合比较大,通常还包含了非物理解.由于这类解的光滑性不够,因此,对解的性质分析起来比较困难.

第三类解是Hoff解.对于Hoff解而言,优点是它保持了解的非线性性和物理上关心的一些有意义的特性.另一方面,这类解当初值具有足够的正则性时,解一般也具有较好的正则性,一般可以根据相关定理来证明解的唯一性和连续性.关于全空间的Cauchy问题,Hoff解的研究已经非常完备了.但是,对于有界区域,由于多了边界,处理起来要相对困难一些.对有界区域,边界条件有很多种提法,如Direchlet边界、No-Slip边界和Navier-Slip边界等.

前人关于Navier-Stokes方程组Direchlet边值问题解的存在性的研究结果比较丰富.例如,对于Direchlet边界条件,文献[6]给出了等熵Navier-Stokes方程

(1)

的解的局部存在性.

但是,关于Slip边值问题的解的局部存在性研究起来更加困难.最近文献[7]证明了等熵的Navier-Stokes方程组带有Slip边值问题的解的局部存在性.而本文研究了非等熵的Navier-Stokes方程组带有Slip边值问题的解的局部存在性.因为，非等熵Navier-Stokes方程组比等熵的Navier-Stokes方程组多了一个能量方程,在研究时除了要估计密度和速度,还要估计温度,而温度和速度是耦合在一起的,方程组更加复杂,因此,研究起来更加困难.

本文研究了带有Navier-Slip边界条件的非等熵Navier-Stokes方程

(2)

初始值条件为

(3)

Navier-Slip边界条件为

(4)

式中,n(x)为∂Ω上的单位外法向量;(·)τ为(·)的切向分量;α为正常数;B为对称矩阵,并且区域Ω和矩阵B满足下列条件：

a.Ω是三维空间中的有界开集,且∂Ω∈C3;

b. B在(ρ,u,θ)的邻域内是C1光滑的.

在本文中,如果不加特别说明,用字母C表示常数.为了方便后面的叙述,构造一个函数空间W.

其中,W是具有H1内积的Hilbert空间,参见文献[7].

现给出解的定义.

a. ρ∈C([0,T]; W1,q(Ω)),(u,θ)∈C([0,T];

H2(Ω))∩L2([0,T];W2,q(Ω))

ρt∈C([0,T]; Lq(Ω)),(ut,θt)∈C([0,T];

L2(Ω))∩L2([0,T];H1(Ω))

其中,φ,ψ,φ在Ω×[0,T]上是Lipschitz连续的.则对于任意t,φ,ψ∈W,有(ρ,u,θ)是式(2)～(4)的解.

现介绍本文的主要结果(定理1).

定理1 令q∈(3,6],且q是常数,如果初始值(ρ0,u0,θ0)以及h,f满足下列条件:

和相容性条件:

则存在时间T*>0和唯一解(ρ,u,θ)满足带有条件式(3)和式(4)的初边值问题,使得

u,θ满足边界条件(4).

2 几个重要引理[8]

对于u∈H1(Ω),有

(6)

引理2 对于p∈[2,6],u∈H1(Ω),存在一个常数C=C(Ω),使得

(7)

引理3 对于u∈H1(Ω),存在常数C=C(Ω),使得

(8)

3 Lame算子和它的谱理论[9]

讨论Lame算子L和特征函数的正则性.设式(2)和式(4)成立,

(9)

式中:β为常数;g为已知函数.

问题等价为：存在u∈W,使得对所有w∈W,有

S(w)+λdivudivw+βuw]dx

由Lax-Milgram定理,存在有界解算子Sβ∶L2(Ω)→W,使得对所有w∈W,满足

(10)

引理4 设Ω⊂3是有界开集,∂Ω∈C3(Ω),μ>0,2μ+3λ>0,B∈L(∂Ω),且B是对称矩阵,则对足够大的β,存在紧自伴算子Sβ∶L2(Ω)→L2(Ω),Sβ∈W,使得式(10)成立.对所有g∈L2(Ω)和w∈W,存在一组标准正交基{wk}k∈L2(Ω),且L2(Ω)⊂W∩C(Ω),使得Sβwk=λk-1wk,λk→.

引理5 设μ>0,2μ+3λ≥0,m≥0,Ω是有界开集,∂Ω∈Cm+1(Ω),B∈Cm+1(Ω),其中,Cm=C(μ,λ,Ω,B,m)依赖于μ,λ,Ω,B,m.若u是式(9)的解,β,g∈Hm(Ω),则u∈Hm+2(Ω),且

(11)

(12)

对所有φ∈H1(Ω),γ足够大且与κ,α,Ω有关.

(13)

4 先验估计

现推导式(2)的解的先验估计,在本文中假设存在常数c0,使得

(14)

然后作如下先验假设：存在c1和T*,使得

(15)

其中,1≤c0≤c1,0

本文中如果不加说明,C是依赖于μ,λ,κ,α,B,Ω以及f和h的范数的正常数.

4.1 密度的估计

关于密度ρ的存在唯一性,Diperna等[10]得到式(16).

(17)

的解.

由式(16)可知,

首先在方程(2)的第2个等式两边同乘以ρq-1,并在Ω上积分,有

(18)

联合式(18)和式(19),利用Gronwall不等式,有

(21)

4.2 热能估计

用θ乘以方程(2)的第3个等式,在Ω上积分,有

由引理1可知,

用θt乘以方程(2)的第3个等式,在Ω上积分,有

根据引理1,可得

对方程(2)的第3个等式关于t求偏导,乘以θt,在Ω上积分,有

根据引理1,可得

由上面2个式子,并且利用Gronwall不等式,有

4.3 速度估计

用u来乘以方程(2)的第2个等式,在Ω上积分,有

用ut乘以方程(2)的第2个等式,在Ω上积分,有

根据引理1,可得

对方程(2)的第2个等式关于t求导,乘以ut,在Ω上积分,有

根据引理1,可得

由式(54)和式(55)可知,对0≤t≤min (T*,T5),利用Gronwall不等式,有

由Cauchy不等式可得

同时,有

由前面的推异过程可得先验估计.

5 线性化方程组的解的存在性

考虑方程(2)的线性化方程组

(56)

现证明该线性化方程组在Navier-Slip边界条件下解的局部存在性.

引理8 令q∈(3,6],设v∈L([0,T];H2(Ω))∩L2([0,T];W2,q(Ω)),Ω是R3中的有界区域,vt∈L([0,T];L2(Ω))∩L2([0,T];H1(Ω)),(ρ0,u0,θ0),f,h是已知函数,且满足正则条件:

边界条件为

则带有条件式(3)和式(4)的初边值问题(56)存在唯一解(ρ,u,θ),使得

其中,(u,θ)满足边界条件式(4).

证明 可由文献[8]中的定理4得到方程(56)的第1个等式的唯一解的存在性和正则性,现在用v代u,类似于式(20),可得到唯一解ρ，并且

线性动量方程和能量方程可以写成如下形式：

现用Galerkin方法得到式(57)中(θ,u)的存在性.令

则对θ,φ∈H1(Ω),u,ψ∈W,式(57)和式(58)的弱形式为

取um(·,0)∈Wm,θm(·,0)∈Vm,则有引理9.

其中,M1与B,L和Ω有关.

类似于引理9的证明,可以证明引理10.

(62)

其中,常数M2与κ,α和Ω有关.

(63)

其中,(,)是L2(Ω)的内积,将θm(t)和um(t)代入式(59),则式(63)转为如下线性常微分方程组：

(64)

现通过迭代方法构造初边值问题(2)～(4)的近似解.

a. 定义u0=0;

b. 设uk-1,k≥1,则(ρk,uk,θk)是下面的初边值问题：

的解,且满足初始条件和边界条件：

6 定理1的证明

由引理8可知强解(ρk,uk,θk)的整体解的存在性.

由前面估计可知,

其中,C是仅依赖于α,μ,κ,B,Ω和f,h的范数的正常数.由Aubin-Lions引理可得

因此,可得

类似于文献[8]中的连续性讨论,可证u∈C([0,T*];H1(Ω))∩C([0,T*];H2(Ω)),并由标准嵌入定理可得,ρ∈C([0,T*];W1,q(Ω)).定理1证毕.

[1]MATSUMURAA,NISHIDAT.Initialboundaryvalueproblemsfortheequationsofmotionofcompressibleviscousandheat-conductivefluids[J].CommunicationsinMathematicalPhysics,1983,89(4):445-446.

[2]DANCHINR.GlobalexistenceincriticalspacesforcompressibleNavier-Stokesequations[J].InventionesMathematicae,2000,141(3):579-614.

[3]DANCHINR.Globalexistenceincriticalspacesforflowsofcompressibleviscousandheat-conductivegases[J].ArchiveforRationalMechanicsandAnalysis,2001,160(1):1-39.

[4]HOFFD.Existenceofsolutionstoamodelforsparse,one-dimensionalfluids[J].JournalofDifferentialEquations,2011,250(132):1083-1113.

[5]HOFFD,SERRED.ThefailureofcontinuousdependenceoninitialdatafortheNavier-Stokesequationsofcompressibleflow[J].SIAMJournalonAppliedMathematics,1991,51(4):887-898.

[6]HOFFD.LocalsolutionsofacompressibleflowproblemwithNavierboundaryconditionsingeneralthree-dimensionaldomains[J].SIAMJournalonMathematicalAnalysis,2012,44(2):633-650.

[7]CHOYG,CHOEHJ,KIMH.Uniquesolvabilityoftheinitialboundaryvalueproblemsforcompressibleviscousfluids[J].JournaldeMathématiquesPuresetAppliquées,2004,83(2):243-275.[8]EVANSLC.Partialdifferentialequations[M].2nded.Providence,RI:AmericanMathematicalSociety,2010.

[9] 刘炳初.泛函分析[M].2版.北京:科学出版社,2004.

[10]DIPERNARJ,LIONSPL.Ordinarydifferentialequations,transporttheoryandSobolevspaces[J].InventionesMathematicae,1989,98(3):511-547.

(编辑:石 瑛)

Local Solutions of a Compressible Flow Problem with Navier-Slip Boundary Conditions

WANG Xuena, YONG Yan

(CollegeofScience,UniversityofShanghaiforScienceandTechnology,Shanghai200093,China)

The study shows that there exist local-in-time solutions of the compressible Navier-Stoke equations with Navier-Slip boundary conditions in a bounded domainΩ⊂3.To this end,the linearization approach,Galerkin method,standard energy method,aprior estimates and uniform convergence theorem have been used.

compressibleNavier-Stokesequations;Navier-Slipboundaryconditions;solutionsexistence

1007-6735(2017)01-0015-10

10.13255/j.cnki.jusst.2017.01.004

2016-12-10

王雪娜(1990-),女,硕士研究生.研究方向：偏微分方程.E-mail：2423948351@qq.com

雍 燕(1982-),女,讲师.研究方向：偏微分方程.E-mail：yongyan_math@126.com

O 175.25

A