劉萍 徐恒睿 楊建榮

1) (電子科技大學(xué)中山學(xué)院電子信息學(xué)院, 中山 528402)

2) (電子科技大學(xué)物理學(xué)院, 成都 610054)

3) (上饒師范學(xué)院物理與電子信息學(xué)院, 上饒 334001)

Boussinesq方程是流體力學(xué)等領(lǐng)域一個(gè)非常重要的方程.本文推導(dǎo)了Boussinesq方程的Lax對(duì).借助于截?cái)郟ainlevé展開, 得到了Boussinesq方程的自B?cklund變換, 以及Boussinesq方程和Schwarzian形式的Boussinesq方程之間的B?cklund變換.探討了Boussinesq方程的非局域?qū)ΨQ, 研究了Boussinesq方程的單參數(shù)群變換和單參數(shù)子群不變解.運(yùn)用Riccati展開法研究了Boussinesq方程, 證明Boussinesq方程具有Riccati展開相容性, 得到了Boussinesq方程的孤立波-橢圓余弦波解.

專題：非線性物理

1 引 言

一般來講, Boussinesq方程可寫為

其中, 下角標(biāo)x和t表示偏微分.Boussinesq方程可以用于描繪淺水波、等離子體、非線性晶格等眾多物理現(xiàn)象[1?5].

由于該方程應(yīng)用廣泛, 一些特殊形式的或者修正的Boussinesq方程被推導(dǎo)出來研究.例如, “壞”Boussinesq方程(也叫不適定Boussinesq方程)的形式為

這個(gè)方程是在1872年由Boussinesq[1]提出來用于描繪淺水波問題的.Benny和 Luke[6]發(fā)現(xiàn)這個(gè)Boussinesq方程非線性弱散色現(xiàn)象的一般近似.“好”Boussinesq方程的形式為

這個(gè)方程是作為描繪弦的非線性振動(dòng)模型提出來的, 也可以用于描繪非線性介質(zhì)材料中的電磁波[7].一種修正的Boussinesq方程的形式為

這個(gè)方程也經(jīng)常被稱為“改進(jìn)的”Boussinesq方程[8],它由流體力學(xué)推導(dǎo)而來, 也可以用于描繪波在磁場中的傳播, 并取代“壞”Boussinesq方程.

很多不同形式的Boussinesq方程, 是方程(1)的特殊形式.本文旨在研究Boussinesq方程(1)的可積性、對(duì)稱性和嚴(yán)格解.在下文中, 如果沒有特殊說明, Boussinesq方程指的是方程 (1).論文結(jié)構(gòu)如下: 在第2節(jié)中, 從一個(gè)簡化的Boussinesq方程的Lax對(duì), 推導(dǎo)出Boussinesq方程(1)的一組Lax 對(duì); 在第 3 節(jié), 對(duì) Boussinesq 方程 (1)進(jìn)行截?cái)?的 Painlevé展 開 , 得 到 Boussinesq 方 程 的B?cklund變換; 第4節(jié)研究了Boussinesq方程的單參數(shù)群變換; 第5節(jié)討論了Boussinesq方程的全點(diǎn)李對(duì)稱性相似解; 第 6節(jié)應(yīng)用 CRE(consistent Riccati expansion, CRE)方法證明了Boussinesq方程的CRE相容性.Boussinesq方程孤立波-周期波在第7節(jié)進(jìn)行了討論; 第8節(jié)是本文的結(jié)論和討論.

2 Boussinesq方程的Lax對(duì)

當(dāng)α=0 ,β=1 ,γ=1/3 時(shí), 方程 (1) 退化成

為了將方程(1)和方程(5)的變量進(jìn)行區(qū)分, 我們將方程(1)中的變量 {u,x,t} 對(duì)應(yīng)地寫成方程(5)中 的 {v,χ,τ} .Weiss[9]通 過 研 究 方 程 (5)的painlevé性質(zhì), 推出了方程 (5)的一組 Lax 對(duì), 其形式如下

方程(1)和方程(5)之間存在標(biāo)度變換

結(jié)合方程(5)的Lax對(duì)(6)式以及標(biāo)度變換, 可以得到方程(1)的Lax對(duì).

定理1(Lax對(duì)定理)

Boussinesq方程(1)具有如下形式Lax對(duì):

這里的l代表譜函數(shù),j表示 {x,t} 的任意函數(shù).

3 與截?cái)郟ainlevé展開相關(guān)聯(lián)的B?cklund變換

截?cái)郟ainlevé展開法, 是分析非線性系統(tǒng)最有效的方法之一[10?12].對(duì) Boussinesq 方程 (1), 可將u展開成

這里的u0,u1,u2和f都是 {x,t} 的函數(shù),f是奇異流函數(shù).將(9)式代入到方程(1)中, 所得到的多項(xiàng)式中,f的所有不同階次的系數(shù)都應(yīng)該為零.由f?6的系數(shù)為零, 可得到

由f?5的系數(shù)為零, 可得

由f?4的系數(shù), 容易得到

將 (10)式–(12)式代入到f?3的系數(shù)中, 得

方程 (13)在 M?bious變換下, 保持形式不變, 因此被稱為Schwarzian形式的Boussinesq方程[9].

將(9)式—(13)式代到方程(1)中, 比較所得方程中f0的系數(shù), 可發(fā)現(xiàn)u0也是 Boussinesq 方程的一個(gè)解, 這表示u=u0是Boussinesq方程的一個(gè) 自 B?cklund 變 換 .而 且 , 對(duì) 以 上 截 斷Painlevé展開進(jìn)行總結(jié), 可得到一個(gè)非自B?cklund變換.

定理2(B?cklund變換定理)

如果f是 Schwarzian形式的 Boussinesq方程 (13)的解, 那么

也是Boussinesq方程(1)的解.

定理3(B?cklund變換定理)

如果f是 Schwarzian形式的 Boussinesq方程 (13)的解, 那么

也是Boussinesq方程(1)的解.

4 單參數(shù)群變換

Boussinesq方程的對(duì)稱σu也相應(yīng)地拓展為滿足下式的四分量對(duì)稱 {σu,σf,σg,σh} ,

對(duì)方程(16), 我們也可以研究它的全點(diǎn)李對(duì)稱.基于這個(gè)目的, 四分量對(duì)稱 {σu,σf,σg,σh} 應(yīng)該滿足Boussinesq方程的線性化的非線性系統(tǒng).按照點(diǎn)李對(duì)稱的方法, 經(jīng)過計(jì)算可得總的對(duì)稱矢量為

各個(gè)對(duì)稱矢量為:

由對(duì)稱矢量(19)式, 可得到六個(gè)單參數(shù)不變子群:

從以上六個(gè)單參數(shù)不變子群, 可到到下列B?cklund變換定理.

定理4(單參數(shù)群變換)

如果{u(x,t),f(x,t),g(x,t),h(x,t)}是拓展的Boussinesq系統(tǒng)(16)的一組解, 則下列函數(shù)也是拓展的Boussinesq系統(tǒng)(16)的一組解,

5 全點(diǎn)李對(duì)稱相似解

對(duì)稱性理論是求解偏微分方程的一種有效系統(tǒng)的方法[13?19].從對(duì)稱矢量 (19)式, 不僅可以得到單參數(shù)不變子群和群不變解, 而且可以得到Boussinesq的相似解和約化方程.將約化方程的嚴(yán)格解和相似解相結(jié)合, 則可以得到所研究系統(tǒng)的嚴(yán)格解.可得到下列四組非平庸情況.

情況 1

在種情這況, 群不變量可寫為

相似解的形式為

情況 2

{σu,σf,σg,σh}包含C4, 而C4是與非局域?qū)ΨQ相關(guān)聯(lián)的, 那么如果令C4=0 , 則相似解會(huì)變得更加簡化.這樣, 相似解為:

與情況一相比, 時(shí)間和空間的對(duì)稱性都沒有改變,因此這種情況的群不變量與情況一相同, 仍為

將 (24b)式代入 (16d)式和 (16e)式, 則變量f和g變成:

將(24b)式代到(16c)式, 可以得到用和F表示的u的表達(dá)式, 將(24b)代入到(16b)式, 可以得到F滿足的約束方程.由于這兩個(gè)式子都很長,此處省略不寫.

情況 3

(18)式和(19)式說明空間x和時(shí)間t的對(duì)稱受到C1的影響.當(dāng)C1=0 時(shí), 群不變量x將比情況一和情況二的群不變量簡單.此時(shí), 群不變量變?yōu)?/p>

相似解為:

其中F(ξ) 滿足

將(29b)式代到(16c)式, 可得到關(guān)于Boussinesq方程的下列B?cklund變換.

定理5 (B?cklund變換定理).

如果F滿足 (30)式, 則 Boussinesq方程的解為

情況 4

這種情況下, 拓展系統(tǒng)(16)的相似解為:

這里, 群不變量x為

將(32b)式代入到(16b)式, 可得到F(ξ) 滿足的約束方程.將(32b)式代入到 (16c)式, 則得到下列定理.

定理6(B?cklund變換定理).

如果F(ξ) 滿足 (32b)式, 則 Boussinesq 方程的解可以寫為

6 Boussinesq方程的CRE相容性

本節(jié)將通過 CRE (consistent Riccati expansion, CRE)方法來討論 Boussinesq方程的嚴(yán)格解[20].Riccati方程的形式為

這里的a0,a1和a2是任意常數(shù).Riccati方程的嚴(yán)格解可寫為

其中,

對(duì)于一個(gè)偏微分系統(tǒng)

我們可假設(shè)它可以展開為

這里的R(w) 是 Riccati方程的嚴(yán)格解.將 (39)式代入到 (38) 式, 并令Ri(w) 的系數(shù)為零, 可得:

如果系統(tǒng) (40)是自洽的, 則展開式 (39)式是“CRE”, 且非線性系統(tǒng) (38)是“CRE”相容系統(tǒng)[20].

為了得到孤立波-周期波碰撞解, 可應(yīng)用CRE方法.CRE方法可被用于證明一個(gè)系統(tǒng)是CRE相容系統(tǒng), 并可用于尋求非線性系統(tǒng)的碰撞波解.對(duì)Boussinesq方程,u可展開成截?cái)嗾归_的形式:

這里,u3,u4,u5和w都是x和t的函數(shù),R(w) 是Riccati方程的一個(gè)解.

將(35)式和(41)式代入到方程(1)中, 并令R(w)所有階次的系數(shù)為零, 可得

這里w滿足

通過CRE和CRE相容性的定義, Boussinesq方程顯然是一個(gè)CRE相容系統(tǒng).基于以上討論, 可得到如下定理:

定理7(CRE相容性定理)

Boussinesq方程是一個(gè)CRE相容系統(tǒng).如果w是相容性條件(43)式的一個(gè)解, 則下列形式的u也是Boussinesq方程的一個(gè)解.

這里的R(w) 和q分別滿足(36)式和(37)式.

7 孤立波-周期波碰撞解

從Boussinesq方程的CRE性質(zhì), 可進(jìn)一步研究Boussinesq方程的嚴(yán)格解.將(36)式代入到(44)式中可得

從(45)式可看到, 如果我們想知道u的具體形式, 那么需要先知道w的表達(dá)式.如果w具有如下形式:

這里k1,k2,ω1,ω2,a3,n和m是常數(shù),Eπ是第三類不完全橢圓積分.將(46)式代入到(43)式中, 并令 s n(k2x+ω2t,m) 的所有不同階次的系數(shù)為零,可發(fā)現(xiàn)參數(shù)應(yīng)該滿足:

這里a4=k1+a3k2.

將(46)式代入到(45)式中, 得:上式中的參數(shù)滿足(47)式或(48)式.

圖1和圖展示了滿足約束關(guān)系(47)的解(49)式.圖1 中的自由參數(shù)選為{n= 0.2,m=0.5,a1= 1,a3= 1,k1= 1,k2= 1,ω2=1,α= ?0.8,β=1}, 圖2中的自由參數(shù)選為{n= 0.2,m= 0.9,a1= 1,a3= 1,k1= 1,k2=1,ω2=1,α= ?0.8,β=1} .圖1 和圖2 展示了亮孤子和周期波的碰撞行為.圖3展示了圖1和圖2中u的密度函數(shù), 圖3(a)對(duì)應(yīng)圖1, 圖3(b)對(duì)應(yīng)圖2.兩種情況的周期波和孤立波的方向是一致的, 而碰撞處的形狀則不相同.

圖1 滿足 (47)式的碰撞波解 (49)式.自由參數(shù)為{n = 0.2,m = 0.5, a1 = 1, a3 = 1, k1 = 1, k2 = 1, w2 = 1, a = –0.8,b = 1}Fig.1.The solution (49) with Formula (47).The free parameters are {n = 0.2, m = 0.5, a1 = 1, a3 = 1, k1 = 1, k2 =1, w2 = 1, a = –0.8, b = 1}.

圖4和圖5展示了滿足參數(shù)限制(48)式的碰撞波解(49)式, 里邊的周期波在扭結(jié)孤立波上運(yùn)動(dòng), 而不是在常數(shù)背景上運(yùn)動(dòng).圖4中的自由參數(shù)選為 {n= 0.4,a1= 1,a2= 1,a3= 2.2,k1= 1,k2= –0.22,ω2=1,α= ?400,β=80} , 其 中(48)式中的m選“+”; 圖5中的自由參數(shù)選為{n= 0.6,a1= 2,a2= 1,a3= 4,k1= 1,k2= –0.12,w2= 0.1,α=?14,β=6} , 其中 (48)式中的m選“–”.圖6展示了圖4和圖5中u的密度函數(shù),圖6(a)對(duì)應(yīng)圖4, 圖6(b)對(duì)應(yīng)圖5.圖6清楚地展示了扭結(jié)孤立波和周期波的碰撞.

圖2 滿足 (47)式的碰撞波解 (49)式.自由參數(shù)為 {n =0.2, m = 0.9, a1 = 1, a3 = 1, k1 = 1, k2 = 1, w2 = 1, a =–0.8, b = 1}Fig.2.The solution (49) with Formula (47).The free parameters are {n = 0.2, m = 0.9, a1 = 1, a3 = 1, k1 = 1, k2 =1, w2 = 1, a = –0.8, b = 1}.

圖3 u 的密度函數(shù)圖.圖 (a)的參數(shù)與圖1 相同, 圖 (b)的參數(shù)與圖2 相同F(xiàn)ig.3.The density of u.The parameters of the Fig.(a) are the same as those of Figure 1 and the parameters of the Fig.(b) are the same as those of Figure 2.

圖4 參數(shù)關(guān)系滿足 (48)式的碰撞波解 (49)式的演化圖.自由參數(shù)為 {n = 0.4, a1 = 1, a2 = 1, a3 = 2.2, k1 = 1, k2 =–0.22, w2 = 1, a = –400, b = 80}Fig.4.The interaction solution (49) with parameter satisfying Formula (48).The free parameters are chosen as {n =0.4, a1 = 1, a2 = 1, a3 = 2.2, k1 = 1, k2 = –0.22, w2 = 1,a = –400, b = 80}.

圖5 參數(shù)關(guān)系滿足 (48)式的碰撞波解 (49)式.自由參數(shù)為{n = 0.6, a1 = 2, a2 = 1, a3 = 4, k1 = 1, k2 = –0.12, w2 =0.1, a = –14, b = 6}Fig.5.The interaction solution (49) with parameter satisfying Formula (48).The free parameters are selected as {n =0.6, a1 = 2, a2 = 1, a3 = 4, k1 = 1, k2 = –0.12, w2 = 0.1,a = –14, b = 6}.

圖6 u 的密度函數(shù)圖.圖 (a)對(duì)應(yīng)圖4, 圖 (b)對(duì)應(yīng)圖5Fig.6.The density of u.The Fig.(a) is related to Fig.4 and the Fig.(b) is corresponding to Fig.5.

8 總結(jié)和討論

本文推導(dǎo)了Boussinesq方程的Lax對(duì), 說明Boussinesq方程是 Lax可積模型.運(yùn)用截?cái)郟ainlevé展開法研究了 Boussinesq方程, 得到了Boussinesq方 程 的 自 B?cklund變 換 , 以 及Boussinesq方 程 和 Schwarzian形 式 的Boussinesq方程之間的非自B?cklund變換.研究了Boussinesq方程的全點(diǎn)李對(duì)稱, 得到了單參數(shù)群變換和單參數(shù)子群不變解.運(yùn)用CRE方法研究了Boussinesq方程, 證明了Boussinesq方程是一個(gè)CRE相容模型, 得到了Boussinesq方程的孤立波-橢圓余弦波碰撞解.Boussinesq方程廣泛地應(yīng)用于描繪流體動(dòng)力學(xué)、電磁學(xué)、等離子體、非線性晶格等物理現(xiàn)象.它作為一個(gè)著名的孤立子方程,各種各樣的激發(fā)模式, 以及它在各種物理情景中的應(yīng)用, 值得不斷深入研究.

感謝樓森岳教授和任博博士的寶貴討論.