張敏心, 陳 敏, 羅 宏

(四川師范大學(xué) 數(shù)學(xué)科學(xué)學(xué)院, 四川 成都 610066)

Boussinesq方程由流體速度場(chǎng)與溫度場(chǎng)耦合而成,是地球物理動(dòng)力學(xué)中的重要模型,該方程被廣泛應(yīng)用于大氣科學(xué)和海洋湍流研究[1-3],并且在Rayleigh-Benard對(duì)流研究中也發(fā)揮了重要作用[4].

本文主要研究如下的二維不可壓縮Boussinesq方程

(1)

其中,Ω?R2是有界開(kāi)集且有光滑邊界?Ω,u=(u1,u2)為速度場(chǎng),T(x,t)為溫度,p為壓強(qiáng),已知函數(shù)f和j分別表示外力和熱源,正常數(shù)γ和κ分別表示黏性系數(shù)和熱擴(kuò)散系數(shù),e2=(0,1).

Boussinesq方程的穩(wěn)態(tài)解對(duì)于理解其動(dòng)力學(xué)行為具有重要意義,關(guān)于二維Boussinesq系統(tǒng)的穩(wěn)態(tài)解問(wèn)題已經(jīng)有了一些研究[5-12].1991年,Morimoto[5]利用Galerkin方法證明了在有界域中Boussinesq近似方程穩(wěn)態(tài)解的存在唯一性.1996年,Lorca等[6]利用譜Galerkin方法結(jié)合不動(dòng)點(diǎn)變?cè)玫搅藦V義Boussinesq方程解的存在性、正則性和唯一性條件.1998年,文獻(xiàn)[7]證明了在外區(qū)域中的廣義Boussinesq模型穩(wěn)態(tài)弱解的存在性.2000年,Lee等[8]研究了二維穩(wěn)態(tài)Boussinesq方程的最優(yōu)控制問(wèn)題,不僅證明了該方程光滑解和最優(yōu)解的存在性,而且利用拉格朗日乘數(shù)法推導(dǎo)出一個(gè)最優(yōu)系統(tǒng).2011年,Alekseev等[9]考慮了在非齊次Dirichlet速度邊界條件和混合溫度邊界條件下的黏性穩(wěn)態(tài)Boussinesq方程最優(yōu)控制的穩(wěn)定性.2012年,Kim[10]利用經(jīng)典的Banach不動(dòng)點(diǎn)定理得到了在有界光滑區(qū)域上具有任意大邊界條件的穩(wěn)態(tài)Boussinesq系統(tǒng)極弱解的存在性,也證明了具有小邊界條件的穩(wěn)態(tài)Boussinesq系統(tǒng)極弱解的唯一性.2021年,黃娟麗等[11]采用分離變量法得到在水平與垂直條帶區(qū)域上Boussinesq方程一些顯式穩(wěn)態(tài)解的結(jié)果.2022年,Chen等[12]用時(shí)間-權(quán)重能量估計(jì),得到了在平面帶區(qū)域上具有部分耗散的二維Boussinesq方程穩(wěn)態(tài)解的整體存在性和穩(wěn)定性,同時(shí)也建立了二維Boussinesq方程解的顯式衰減率.

本文借助文獻(xiàn)[13-14]的方法,研究二維Boussinesq方程穩(wěn)態(tài)解的存在性與正則性.首先,利用弱連續(xù)算子銳角原理證明二維穩(wěn)態(tài)Boussinesq方程弱解的存在性;其次,利用ADN理論以及線性橢圓方程理論證明方程強(qiáng)解的存在性,得到強(qiáng)解的范數(shù)估計(jì)式,并通過(guò)Sard-Smale定理得到解的有限性;最后,再次利用ADN理論以及線性橢圓方程理論得到方程古典解的存在性,實(shí)現(xiàn)了解的正則性的提升.

1 預(yù)備知識(shí)

本節(jié)介紹一些預(yù)備知識(shí),包括弱連續(xù)算子銳角原理、ADN理論、線性橢圓方程理論和Sard-Smale定理等知識(shí).

引理1.1[15]對(duì)于如下的橢圓方程:

(2)

0<λ1≤λ2≤…,

λk→∞,k→∞.

引理 1.2[15-16](Sobolev嵌入定理) 令Ω?R2是一個(gè)有界開(kāi)集,對(duì)任意區(qū)域Ω有:

Wk,p0(Ω)

特別地,這個(gè)嵌入是連續(xù)的,并且有以下不等式成立

其中,C=C(n,p)是依賴于n,p的常數(shù).

設(shè)X是一個(gè)線性空間,X1、X2是Banach空間,X1是可分的,X2是自反的.設(shè)X?X2,于是存在一個(gè)線性映射L:X→X1是一一稠密的.

定義 1.1[15](弱連續(xù)算子的定義) 映射G:X2→X1*是弱連續(xù)的,如果對(duì)所有的{Φn}?X2,Φn?Φ0有

引理 1.3[15](弱連續(xù)算子銳角原理) 設(shè)G:X2→X1*是弱連續(xù)的,U?X2是有界開(kāi)集,0∈U.若

〈G(Φ),LΦ〉≥0, ?Φ∈?U∩X,

則方程G(Φ)=0在X2中存在一個(gè)解.

引理 1.4[17-18](ADN理論) 考慮如下Stokes方程

(3)

則有如下結(jié)論:

1) 令(u,p)∈C2,α(Ω,Rn)×C1,α(Ω)是問(wèn)題(3)的解.若f∈Ck,α(Ω,Rn),k≥0,則

(u,p)∈Ck+2,α(Ω,Rn)×Ck+1,α(Ω),

且滿足如下估計(jì)

‖u‖Ck+2,α+‖p‖Ck+1,α≤

C(‖f‖Ck,α+‖(u,p)‖C0),

其中,C>0僅依賴于μ、n、k、α、Ω.

2) 令(u,p)∈W2,p(Ω,Rn)×W1,p(Ω)(1

(u,p)∈Wk+2,p(Ω,Rn)×Wk+1,p(Ω),

且有

‖u‖Wk+2,p+‖p‖Wk+1,p≤

C(‖f‖Wk,p+‖(u,p)‖Lp),

這里C>0只依賴于μ、n、k、p、Ω.

引理 1.5[16](線性橢圓方程理論) 考慮如下橢圓方程

(4)

其中,aij,bi,c∈L∞(Ω),f∈L2(Ω),aij=aji是一致橢圓的,即存在常數(shù)0<λ1≤λ2使得

λ1|ξ|2≤aijξiξj≤λ2|ξ|2,

?ξ∈Rn,x∈Ω,

則有如下結(jié)論成立:

1) 令Ω?R2是C2,α的,aij,bi,c,f∈C0,α(Ω).如果u∈C2,α(Ω,Rn)是問(wèn)題(4)的解,那么

‖u‖C2,α≤C(‖f‖C0,α+‖u‖C0),

其中,C>0僅依賴于n、α、λ、Ω和aij、bi、c的C0,α范數(shù).

2) 令Ω?R2是C2的,aij∈C0(Ω),bi,c∈L∞(Ω),f∈Lp(Ω).如果u∈W2,p(Ω,Rn)是問(wèn)題(4)的解,那么

‖u‖W2,p≤C(‖f‖Lp+‖u‖Lp),

其中,C>0僅依賴于n、p、λ、Ω和aij、bi、c的C0范數(shù)或L∞范數(shù).

引理 1.6[19-20](Sard-Smale定理) 令G:X→Y是C1的零指標(biāo)的Fredholm算子,則G的正則值在Y中稠密.進(jìn)一步,若p∈Y是G的臨界點(diǎn),則G-1(p)是一個(gè)離散集.

2 主要結(jié)論

引入空間

H={Φ=(u,T)∈L2(Ω,R3)|Φ滿足

H1=H1(Ω,R3)∩H.

證明對(duì)于?Ψ=(v,s)∈H1,定義G:H1→H1*,

fv-Tve2+κ▽T·▽s+(u·▽)T·s-js]dx.

首先,驗(yàn)證G算子滿足

〈GΦ,Φ〉≥0.

(5)

根據(jù)Young不等式以及Poincaré不等式,可得

κ▽T·▽T-jT)dx≥

‖T‖‖u‖-‖j‖‖T‖≥

‖j‖‖T‖=

‖f‖‖u‖-‖j‖‖T‖≥

則存在恰當(dāng)大的常數(shù)N使得

〈GΦ,Φ〉≥0, ?Φ∈?BN∩X,

于是(5)式成立.

再驗(yàn)證G的弱連續(xù)性.假設(shè)在空間H1中有Φk?Φ,根據(jù)引理1.2可得

uk→u0在Lp(Ω,R2)中,Tk→T0在

Lp(Ω)中,其中1≤p<∞.

(6)

dx+

dx=

(7)

作下面分解

dx+

由于在H01(Ω,R2)中有un?u0,即在L2(Ω,R2)中有Dun?Du0.同理,在H01(Ω)中有Tn?T0,即在L2(Ω)中有DTn?DT0,并且u0·v∈L2,u0·s∈L2,因此

dx=0,

由H?lder不等式和(6)式可得

dx|≤

于是(7)式成立.由引理1.3可得,問(wèn)題(1)存在一個(gè)弱解Φ∈H1.

證明由H?lder不等式,可得

(8)

考慮如下的Stokes方程

由f∈Lq(Ω,R2),T∈L2(Ω)和(8)式可知

因此,由引理1.4可得

對(duì)于橢圓方程:

(10)

由j∈Lq(Ω)和(8)式可知

因此,由引理1.5可得

由引理1.2可得

1≤q≤6,

因此

于是可推出

(u·▽)u∈L2(Ω,R2),

(u·▽)T∈L2(Ω).

(11)

由(11)式、引理1.4和引理1.5可得

(u,T,p)∈W2,2(Ω,R2)×W2,2(Ω)×W1,2(Ω).

由引理1.2可得

1≤q<∞.

可以推出

(u·▽)u∈L3(Ω,R2),

(u·▽)T∈L3(Ω).

利用引理1.4和引理1.5可得

(u,T,p)∈W2,3(Ω,R2)×W2,3(Ω)×W1,3(Ω).

如此重復(fù)以上步驟,便證得(1)式的解(u,T,p)∈W2,q(Ω,R2)×W2,q(Ω)×W1,q(Ω),q≥2.

證明對(duì)于(1)式的解(u,T,p)需要證明如下估計(jì):

‖u‖W2,q+‖T‖W2,q+‖p‖W1,q≤

C(‖f‖Lq+‖j‖

q≥2.

(12)

由引理1.4和(9)式可得

‖u‖W2,q+‖p‖W1,q≤C‖g‖Lq≤

C(‖f‖Lq+‖T‖Lq+‖u·▽u‖Lq).

由引理1.5和(10)式可得

‖T‖W2,q≤C‖h‖Lq≤

C(‖j‖Lq+‖u·▽T‖Lq),

因此

‖u‖W2,q+‖T‖W2,q+‖p‖W1,q≤

C(‖f‖Lq+‖T‖Lq+‖u·▽u‖Lq+

‖j‖Lq+‖u·▽T‖Lq).

(13)

由H?lder不等式可得

‖u·▽u‖Lq≤‖u‖L2q‖▽u‖L2q,

‖u·▽T‖Lq≤‖u‖L2q‖▽T‖L2q.

(14)

由Gagliardo-Nirenberg不等式以及Young不等式可得

‖▽u‖

ε‖Δu‖Lq+Cε-1‖▽u‖,

(15)

‖▽T‖

ε‖ΔT‖Lq+Cε-1‖▽T‖.

(16)

由〈GΦ,Φ〉=0,有

‖▽u‖2+‖▽T‖2≤

C(‖f‖2+‖j‖2),

因此

‖▽u‖≤C(‖f‖Lq+‖j‖Lq),

‖▽T‖≤C(‖f‖Lq+‖j‖Lq),q≥2.(17)

由引理1.2和(17)式有

‖u‖L2q≤C‖▽u‖≤

C(‖f‖Lq+‖j‖Lq).

(18)

由(14)、(15)、(17)和(18)式可以得到

‖u·▽u‖Lq≤C(‖f‖Lq+‖j‖Lq)×

[ε‖Δu‖Lq+Cε-1(‖f‖Lq+‖j‖Lq)]≤

C(ε‖f‖Lq‖Δu‖

2ε-1‖f‖Lq‖j‖Lq+

于是由H?lder不等式,可得

‖u·▽u‖Lq≤C(ε‖f‖Lq‖Δu‖Lq+

ε‖j‖Lq‖Δu‖Lq).

(19)

同理,由(14)式和(16)～(18)式可以得到

‖u·▽T‖Lq≤C(‖f‖Lq+‖j‖Lq)×

[ε‖ΔT‖Lq+Cε-1(‖f‖Lq+‖j‖Lq)]≤

C(ε‖f‖Lq‖ΔT‖

2ε-1‖f‖Lq‖j‖Lq+

于是由H?lder不等式,可得

‖u·▽T‖Lq≤C(ε‖f‖Lq‖ΔT‖Lq+

ε‖j‖Lq‖ΔT‖Lq).

(20)

由引理1.2和(17)式有

‖T‖Lq≤C‖▽T‖≤

C(‖f‖Lq+‖j‖Lq).

(21)

由(13)、(19)～(21)式有

‖u‖W2,q+‖T‖W2,q+‖p‖W1,q≤

C[‖f‖Lq+‖j‖

ε(‖f‖Lq+‖j‖Lq)‖u‖W2,q+

ε(‖f‖Lq+‖j‖Lq)‖T‖W2,q].

(22)

‖u‖W2,q+‖T‖W2,q+‖p‖W1,q≤

C(‖f‖Lq+‖j‖

‖j‖

(23)

最后利用Young不等式和(23)式證得(12)式成立.

最后,引入如下映射

G=L+H:

W2,q(Ω,R3)×W1,q(Ω)→Lq(Ω,R3),

其中

令

于是方程(1)可以寫成

1) 對(duì)任意(f,j)∈Cα(Ω,R3),問(wèn)題(1)存在一個(gè)古典解(u,T,p)∈C2,α(Ω,R2)×C2,α(Ω)×C1,α(Ω);

3) 如果(f,j)∈C∞(Ω,R3),那么問(wèn)題(1)的解(u,T,p)∈C∞(Ω,R4).

證明1) 因?yàn)镃α(Ω)?Lq(Ω),?q≤∞,所以對(duì)任意(f,j)∈Cα(Ω,R3),由定理2.2可知問(wèn)題(1)存在一個(gè)強(qiáng)解

(u,T,p)∈W2,q(Ω,R2)×W2,q(Ω)×

W1,q(Ω), ?q≥2.

由引理1.2知

W2,q(Ω,R2)

于是

(u,T)∈C1,α(Ω,R3),

因此

(u·▽)u∈Cα(Ω,R2), (u·▽)T∈Cα(Ω).

由(9)式中的g∈Cα和(10)式中的h∈Cα,以及引理1.4的結(jié)論1)可知

(u,p)∈C2,α(Ω,R2)×C1,α(Ω).

由引理1.5的結(jié)論1)可知

T∈C2,α(Ω),

即

(u,T,p)∈C2,α(Ω,R2)×C2,α(Ω)×C1,α(Ω).

結(jié)論1)成立.

2) 由于Cα(Ω)?Lq(Ω),?q≤∞,再根據(jù)定理2.3可得,問(wèn)題(1)的解是有限的.

3) 由于(f,j)∈C∞(Ω),有(f,j)∈Wk,q(Ω)(k是任意整數(shù)).根據(jù)定理2.2,可以得到(u,T,p)∈W2,q(Ω,R3)×W1,q(Ω).根據(jù)引理1.4以及引理1.5可得(u,T,p)∈Wk+2,q(Ω,R3)×Wk+1,q(Ω)(k是任意整數(shù)).利用引理1.2可得(u,T,p)∈Ck+1,α(Ω,R3)×Ck,α(Ω)(k是任意整數(shù)).于是(u,T,p)∈C∞(Ω,R4).

3 總結(jié)與展望

本文給出了Boussinesq方程穩(wěn)態(tài)解存在的條件,得到強(qiáng)解的存在性和強(qiáng)解的范數(shù)估計(jì),實(shí)現(xiàn)了解的正則性提升,并得到方程古典解的存在性.

本文考慮的Boussinesq方程形式相對(duì)簡(jiǎn)單,鑒于Boussinesq方程有很多的形式,比如具有水平(垂直)耗散、水平(垂直)熱擴(kuò)散、分?jǐn)?shù)耗散的Boussinesq方程,對(duì)于是否能用弱連續(xù)算子方法考慮其他形式的Boussinesq方程解的存在性與正則性,將是進(jìn)一步所需要考慮的問(wèn)題.此外,本文討論的是穩(wěn)態(tài)Boussinesq方程解的存在性,對(duì)于是否能用弱連續(xù)算子方法研究演化Boussinesq方程解的存在性與正則性,也將是進(jìn)一步所需要考慮的問(wèn)題.