康 靜，鄧大文

（湘潭大學(xué) 數(shù)學(xué)與計(jì)算科學(xué)學(xué)院，湖南 湘潭 411105）

Boussinesq方程組描述了不可壓流體在有熱對(duì)流時(shí)的流動(dòng)情況，模擬大氣、洋流等實(shí)際現(xiàn)象[1-3]。二維無(wú)粘性無(wú)熱傳導(dǎo)的Boussinesq方程組的全局正則性是個(gè)開(kāi)問(wèn)題，即還不清楚局部光滑解是否一定可以延拓為全局光滑解。R2和二維光滑區(qū)域上不可壓Euler方程組具有全局正則性[4-6]，Kiselev和Zlatos[7]在一個(gè)有尖點(diǎn)的區(qū)域上，對(duì)該方程組構(gòu)造了一類在有限時(shí)間爆破的局部光滑解。對(duì)無(wú)黏性無(wú)熱傳導(dǎo)的Boussinesq方程組，在上述3種區(qū)域上的全局正則性都是開(kāi)問(wèn)題，本文用Kiselev和Zlatos的方法，在他們使用的有尖點(diǎn)區(qū)域上，構(gòu)造一類在有限時(shí)間內(nèi)爆破的局部光滑解。

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

我們先回憶一般二維有界光滑區(qū)域Ω上的Boussinesq方程組

其中 u=(u1(x,t),u2(x,t))是速度場(chǎng)，θ=θ(x,t)∈R是溫度場(chǎng),p=p(x,t)∈R是壓力場(chǎng)，x=(x1,x2)∈Ω，t＞0，e2=(0,1)，n是?Ω的單位外法向量，u在?Ω上滿足u·n=0。本文約定，(1)i指式（1）中的第i條方程，i=1,2,3?？梢园咽剑?）寫(xiě)成渦量-流函數(shù)形式，令ω=為渦量，又由div u=0，知存在流函數(shù) ψ：Ω→R ，使 u=?⊥ψ ,其中 ?⊥:=(-?x2,?x1)。由式 (1)4，可知ψ|?Ω為常數(shù)，不妨設(shè)為0。所以式（1）可寫(xiě)成

可參看文獻(xiàn)[4，8]等對(duì)Euler方程渦量-流函數(shù)形式的描述。由式(2)3和(2)4

其中GΩ(·,·)為區(qū)域 Ω 上的格林函數(shù)的負(fù)，即

其中 KΩ(x,y)=GΩ(x,y)。注意在 Ω×Ω 上，GΩ(·,·)＜ 0 ，GΩ(x,y)|x∈?Ω=0 ，即 ?Ω 是 GΩ(·,y)的0-等值線。所以?xGΩ(x,y)與?Ω垂直，指向Ω外，GΩ(x,y)與?Ω 相切，指向反時(shí)針?lè)较颉Ｔ赗2上式（1）的局部光滑解存在性定理，可參看文獻(xiàn)[10]。我們將用到二維有界光滑區(qū)域上Boussinesq方程組的局部存在性定理：若ω和θ的初值ω0和θ0都在 C1(Ω)，則存在T＞0，使式（2）在 Ω×[0,T]上有C1解?？捎妙愃朴诙SEuler方程的Yudovich定理的證明方法得到，可參看文獻(xiàn)[5]和本文末段的概述。

以下描述本文中的尖點(diǎn)區(qū)域D'。我們沿用文獻(xiàn)[7]中的區(qū)域（圖1）。令

為兩個(gè)分別以(0,1)和(0,-1)為中心的球場(chǎng)形區(qū)域，每個(gè)的高度是 2r＞0（圖1(a)）。注意只是C1,1區(qū)域。通過(guò)在圓弧和直線段的接口處適當(dāng)?shù)馗膭?dòng)圓弧，可得無(wú)窮光滑的區(qū)域。這個(gè)改動(dòng)可以任意小，使和的面積任意接近。最后，令L ，其中 L=(-1,1)×{0}為開(kāi)線段，用D:=表示D'的上半部，如圖1(b)所示。

在D'上考慮式(2)，設(shè)其上的光滑初值ω0和θ0對(duì)x2反對(duì)稱，且在L的一個(gè)鄰域內(nèi)為0，則從光滑區(qū)域上的局部存在定理得到它們有互為正負(fù)的局部光滑解，由Dθ/Dt=0和Dω/Dt=θx1，容易知道對(duì)小的時(shí)間 t，ω(·,t)和 θ(·,t)在 L 的一個(gè)鄰域內(nèi)仍為0，所以兩個(gè)解光滑地粘合成為D'上的局部光滑解。注：這樣的在D'上的解的速度場(chǎng)u在x∈D處的值可由u(x,t)=(x,y)ω(y,t)dy給出。

圖1 尖點(diǎn)區(qū)域

本文用Kiselev和Zlatos[7]的方法，在 D'上對(duì)無(wú)黏性無(wú)熱傳導(dǎo)Boussinesq方程組構(gòu)造在有限時(shí)間內(nèi)爆破的局部光滑解。與文獻(xiàn)[7]中對(duì)二維Euler方程的討論不同，從式(2)1可見(jiàn)，當(dāng)θx1非零時(shí)（若要構(gòu)造式(2)非平凡的解，θx1就不能恒為0），渦量不是簡(jiǎn)單地被搬運(yùn)，而是隨著流體質(zhì)點(diǎn)的流動(dòng)可以變動(dòng)，以致變號(hào)，為估計(jì)邊界上質(zhì)點(diǎn)的移動(dòng)速度增添了難度。本文通過(guò)選取適當(dāng)?shù)某跏紲囟圈?，限制了在任何時(shí)候θx1≠0的區(qū)域的面積，從而解決了問(wèn)題。

2 定理及其證明

在陳述主要定理前，先證一個(gè)引理。

又因?yàn)?D在每個(gè)x∈?D處都滿足內(nèi)球條件，所以由Hopf引理，對(duì)任意的x∈?D，

(?v·n)(x)=(?v/?n)(x)＞0

從而 ?v(x)≠0→，且因 ?D 是緊的，所以

ε：=inf{ ||?v(x)，x∈?D}＞0斷言得證。

以下我們完成引理1的證明。因?yàn)?D是GD(·,y) 和 v(·) 的 0-等 值 線 ，所 以 對(duì) x∈ ?D ，?xGD(x,y)和?v(x)都與?D垂直，都指向外。因此，在?D上，對(duì)y∈

定理1帶有尖點(diǎn)的區(qū)域D'上二維Boussinesq方程組（2）存在有限時(shí)間內(nèi)爆破的局部光滑解。詳細(xì)地，若初始溫度和渦量 θ0, ω0∈C1(Dˉ')對(duì) x2反對(duì)稱，在L的一個(gè)鄰域內(nèi)為0，在D內(nèi)它們非負(fù)，θ0只在足夠小的區(qū)域上非0，ω0在θ0的支撐內(nèi)不恒為零，且ω0在?DL上與尖點(diǎn)P（圖1(b)）的距離足夠小的一點(diǎn)處非0，則ω在有限時(shí)間內(nèi)變得不連續(xù)。

證明：選取對(duì) x2反對(duì)稱的ω0,θ0，從而由Biot-Savart公式知速度場(chǎng)u相對(duì)于x1-軸對(duì)稱，使在某時(shí)間t0＞0前在?D上的質(zhì)點(diǎn)最低以某速度沿反時(shí)針?lè)较蛄鲃?dòng)，且在足夠接近尖點(diǎn)P的某點(diǎn)x'∈?D處ω0(x')＞0，且保證從x'出發(fā)的質(zhì)點(diǎn)流動(dòng)時(shí)渦量不變。則從ω0對(duì) x2的反對(duì)稱性，當(dāng) x'與它相對(duì)于x1-軸的鏡像x''同時(shí)流到L時(shí)，ω就在L處不連續(xù)。以下是詳細(xì)證明。

第一步：選初始渦量和溫度ω0,θ0：D'→R滿足以下條件

(i) ω0,θ0∈C1(Dˉ')，對(duì) x2反對(duì)稱；

(ii)在 D 內(nèi)，ω0,θ0≥0，不恒等于0；

(iii)ω0,θ0各在 L的某開(kāi)鄰域 Nω0,Nθ0內(nèi)為0；

(v)∫W0ω0(y)dy＞0 令，即 ω0在W0內(nèi)不恒為0；

(vi)ω0在?DL離P足夠近的一點(diǎn)處非零（第三步中詳細(xì)描述）。

設(shè)X(x,t)為從x出發(fā)的質(zhì)點(diǎn)軌跡，即滿足的唯一解。令

所以在 DWt上 ω(·,t)≥0 。事實(shí)上，若 ?W0，則存在的一個(gè)鄰域 N，N?W0=φ ，其上 θ≡0 ，從而θx1

≡0，Dω/Dt=0。因?yàn)?X(·,t)是 D 上的同胚，所以

第二步：由光滑區(qū)域上局部光滑解的存在性，可知存在T＞0，使在 D×[0,T]上式（2）有以 ω0,θ0為初值的光滑解，若它在有限時(shí)間內(nèi)爆破，則定理1結(jié)論已成立。以下設(shè)它可以延拓為D上的全局光滑解。

我們斷言：若x∈?D，u(x,t)=∫DKD(x,y)ω(y,t)dy 與 ?D 相切，指向反時(shí)針?lè)较颍掖嬖诔?shù)C＞0，t0＞0，使當(dāng)t∈[0,t0]時(shí)，有|u(x,t)|≥C。

為證明此斷言，取δ＞0，使

以上所涉及的區(qū)域可參看圖2。把速度u分解為

圖2 區(qū)域D，和Wt

以下分別估計(jì)I,II和III。從式（4）后的討論可知，對(duì) x∈?D,y∈D,KD(x,y)與?D相切，指向反時(shí)針?lè)较?。因?yàn)樵?D)Wt上 ω(y ,t)≥0，所以I指向反時(shí)針?lè)较?。?duì)于II,同樣地在積分區(qū)域里ω(y,t)≥0，令

由式（8）可得，|S|≥δ,ω≥δ 。令τ為 ?D 上指向反時(shí)針?lè)较虻膯挝磺邢蛄浚瑒t

由引理1，II是長(zhǎng)度不小于εδ2，跟?D相切，指向反時(shí)針?lè)较虻南蛄?。以下估?jì)III。由積分中值定理，

分以下兩種情況討論：若對(duì)所有 t≥0，ω(y0,t)≥0，則在任何時(shí)間III都指向反時(shí)針?lè)较颍胾(x,t)滿足 |u(x,t)|≥εδ2，指向反時(shí)針?lè)较虻??D 的切向量(此時(shí)取t0為任意正數(shù))。若存在時(shí)間t，使ω(y0(t),t)＜0，則令t0為使

第三步：從第二步可知直到時(shí)間t0，?D上的質(zhì)點(diǎn)沿反時(shí)針?lè)较蛑辽僖苿?dòng)εδ2t0?，F(xiàn)在選x'∈?D,x'在令θ0消失的L的開(kāi)鄰域 Nθ0內(nèi)（從而 x'?W0，ω(X(x',t),t)=ω0(x')），并從x'沿反時(shí)針?lè)较虻郊恻c(diǎn)P的距離小于εδ2t0，那么在t0之前的某時(shí)間t1，從x'出發(fā)的質(zhì)點(diǎn)將流到L中，又由ω對(duì)x2的反對(duì)稱性（從而由Biot-Savart公式u1和u2分別對(duì) x2對(duì)稱和反對(duì)稱），在上x(chóng)'相對(duì)于 x1-軸的鏡像 x''也同時(shí)流到 L。有兩種情況：若 ω0(x')＞0，則ω0(x'')＜0，ω(·,t1)在 L 處將不連續(xù)；若 ω0(x')=0，則對(duì)ω0在x'附近作小改動(dòng)使得ω0(x')＞0，這樣增加ω0后，I仍然指向反時(shí)針?lè)较?。?duì)已選的δ，式（6）依然成立。所以在t0之前的某t1，仍有ω(·,t1)在L處不連續(xù)。定理1證畢。

以下粗略描述一般二維光滑區(qū)域Ω上Boussinesq方程組經(jīng)典解的局部存在性。把式（2）補(bǔ)上u和θ的初值條件(10)5、 (10)6，把式(2)3、 (2)4換成等價(jià)的式(10)3、 (10)4后得到

存在性的證明平行于文獻(xiàn)[5]中對(duì)二維Euler方程的討論。先假定 ω0∈L∞(Ω)，θ0∈C1(Ω)，證明式（10）存在局部弱解，即存在T＞0，及

n=1,2，…。可證式（12）有解，并在區(qū)域Ω×[0,T]內(nèi)收斂到式（11）的解 ω,θ∈L∞(Ω×[0,T])，然后證明它們是式（10）的弱解。進(jìn)一步，當(dāng)ω0∈C1(Ω)時(shí)，可用類似于文獻(xiàn)[5]的方法證明

ω,θ∈C1(Ω×[0,T])。

3 結(jié)論

本文考慮無(wú)黏性無(wú)熱傳導(dǎo)的Boussinesq方程組的全局正則性。通過(guò)粘合兩個(gè)光滑區(qū)域上的光滑解，我們得到一個(gè)在帶尖點(diǎn)區(qū)域上有限時(shí)間失去光滑性的局部光滑解。所以在這類區(qū)域上，這組方程組不具有全局正則性。