廣義變系數(shù)KdV方程的保角能量守恒方法

郭峰, 莊清渠

(華僑大學(xué) 數(shù)學(xué)科學(xué)學(xué)院, 福建 泉州 362021)

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

KdV方程是一類重要的非線性模型,用于描述許多物理現(xiàn)象,如非諧晶體中的聲波、熱等離子體中的磁流體動(dòng)力波及離子聲波等.在各種實(shí)際物理背景中,當(dāng)考慮媒介的非均勻性及邊界的非一致性時(shí),變系數(shù)模型往往比相應(yīng)的常系數(shù)模型更能說明問題.考慮帶依時(shí)阻尼和彌散項(xiàng)的廣義KdV(TDKdV)方程的周期初值問題，即

ut+upux+α(t)u+β(t)uxxx=0,x∈Ω,t∈(0,T],

(1)

u(x,t)=u(x+L,t),x∈Ω,t∈(0,T],

(2)

u(x,0)=u0(x),x∈Ω.

(3)

式(1)～(3)中:p是正整數(shù);Ω=[a,b];T是時(shí)間終點(diǎn)；L=b-a;u0(x)是給定的光滑函數(shù);u(x,t)是相關(guān)波模型的振幅;α(t),β(t)分別表示依時(shí)阻尼系數(shù)和彌散系數(shù).

方程(1)有著廣泛的物理應(yīng)用,如海洋動(dòng)力學(xué)中的海岸波、等離子體物理學(xué)、天體物理學(xué)及流體動(dòng)力學(xué)等[1-3].對(duì)于無耗散的KdV方程,已有的數(shù)值研究方法主要有有限差分法[4-5]、有限元法[6]、譜方法[7-8]等.近年來,保結(jié)構(gòu)算法因能長時(shí)間保持哈密爾頓系統(tǒng)的幾何結(jié)構(gòu)而備受學(xué)者們的關(guān)注.Wang等[9]研究了無耗散KdV方程的局部保結(jié)構(gòu)算法.Brugnano等[10]利用哈密爾頓邊界值方法，構(gòu)造了KdV方程的高階能量守恒格式.對(duì)于具耗散的KdV方程,房少梅[11]在假設(shè)數(shù)值解有界的前提下,僅對(duì)廣義KdV方程的一般譜格式和擬譜格式進(jìn)行解的誤差估計(jì).Guo[12]對(duì)阻尼KdV方程提出一個(gè)保角多辛算法.目前,尚無關(guān)于阻尼KdV方程保能量耗散方面的保角保結(jié)構(gòu)算法.

哈密爾頓系統(tǒng)最具特征的結(jié)構(gòu),就是常常被稱為系統(tǒng)能量的哈密爾頓函數(shù)本身.由于耗散的KdV方程非哈密爾頓系統(tǒng),傳統(tǒng)的保結(jié)構(gòu)方法若直接應(yīng)用到此類問題就會(huì)失去其優(yōu)勢.針對(duì)帶線性阻尼項(xiàng)的保角哈密爾頓系統(tǒng),文獻(xiàn)[13-16]提出能保持其保角守恒律的保角方法.本文根據(jù)這一方法,利用算子分裂技巧,對(duì)方程(1)～(3)構(gòu)造一個(gè)空間上高精度的傅里葉擬譜保角能量守恒格式,該格式在周期邊界條件下能保持系統(tǒng)的離散保角能量守恒律及保角質(zhì)量守恒律.

2 TDKdV方程的保角哈密爾頓系統(tǒng)及保角守恒律

TDKdV方程(1)可以寫成一個(gè)保角哈密爾頓系統(tǒng)[13,15]，即

(4)

式(4)中：H為能量泛函，表達(dá)式為

(5)

保角哈密爾頓系統(tǒng)(4)在合適的邊界條件如周期邊界條件(2)下,具有以下2個(gè)保角性質(zhì).

1) 保角質(zhì)量守恒律

(6)

2) 保角能量守恒律

(7)

3 TDKdV方程的保角能量守恒傅里葉擬譜方法

對(duì)上式求偏導(dǎo)，可得算子IJu(x,t)在節(jié)點(diǎn)xj處的s階偏導(dǎo)數(shù)為

(8)

式(8)中：u=(u(x1,t),u(x2,t),…,u(xJ,t))T；Ds∈RJ×J是s階譜微分矩陣.

D1的元素為

D2的元素為

上式中：j,l=1,2,…,J.

引理1矩陣D1的元素滿足(D1)j=0,j=1,2,…,J,即D1每行元素之和為零.

證明:將u=(1,1,…,1)T,s=1代入式(8),即證明.

在空間上用傅里葉擬譜方法離散方程(4),得到一個(gè)半離散系統(tǒng),即

(9)

式(9)中：

(10)

(11)

將空間上采用傅里葉擬譜方法得到的格式稱為FPEP,采用中心差分方法得到的格式稱為CDEP.

首先,將系統(tǒng)(9)分裂成一個(gè)保守的哈密爾頓系統(tǒng)

(12)

和一個(gè)耗散系統(tǒng)

(13)

然后,在時(shí)間方向上,應(yīng)用二階平均向量場方法(AVF)[18]求解系統(tǒng)(12),得

不失一般性地,設(shè)p=1,即

定理1式(14)具有離散保角質(zhì)量守恒律,即

(16)

證明:注意到D1的反對(duì)稱性,由引理1,對(duì)式(14)關(guān)于空間指標(biāo)j求和,有

從而可得保角質(zhì)量守恒律(16).

定理2式(15)滿足離散保角能量守恒律,即

(17)

4 數(shù)值試驗(yàn)

對(duì)FPEP(15),CDEP(11)及下列傳統(tǒng)AVF方法(TAVF)的解的精度、計(jì)算效率及保角守恒律誤差進(jìn)行比較,并考慮不同類型的阻尼對(duì)保角守恒律的影響,即

TDKdV方程(1)在p=1時(shí),有解析解為

算例1設(shè)阻尼系數(shù)α(t)=0.01，當(dāng)h=0.01,T=10時(shí),FPEP和CDEP的時(shí)間收斂階,如表1所示.τ=0.001,T=10時(shí)，F(xiàn)PEP和CDEP的空間誤差，如表2所示.由表1,2可知：在時(shí)間方向上，F(xiàn)PEP和CDEP均為二階精度；在空間方向上,FPEP的誤差非常小,而CDEP的誤差較大,收斂率為二階.

表1 FPEP和CDEP的時(shí)間收斂階Tab.1 Temporal convergence rates of FPEP and CDEP

表2 FPEP和CDEP的空間誤差Tab.2 Spatial errors of FPEP and CDEP

表3 不同空間步長下2種方法的CPU運(yùn)行時(shí)間Tab.3 CPU times of two methods at different spatial steps s

當(dāng)τ=0.01，T=10時(shí),不同空間步長下2種方法的CPU運(yùn)行時(shí)間，如表3所示.由表3可知：使用FFT技巧的FPEP的運(yùn)算時(shí)間遠(yuǎn)小于CDEP.

在h=0.2,τ=0.001下，分別將3種格式運(yùn)行到T=60時(shí)，其數(shù)值解與精確解的比較，如圖1所示.圖1中：x為空間區(qū)域.由圖1可知:TAVF的誤差最大;FPEP,CDEP的數(shù)值解與精確解吻合得很好.

(a) FPEP (b) CDEP (c) TAVF圖1 3種方法的數(shù)值解與精確解的比較Fig.1 Comparisons of exact solution and numerical solutions with three methods

(a) FPEP (b) CDEP (c) TAVF圖2 3種方法的保角能量守恒律誤差Fig.2 Errors on conformal energy conservation law for there methods

(a) FPEP (b) CDEP (c) TAVF圖3 3種方法的保角質(zhì)量守恒律誤差Fig.3 Errors on conformal mass conservation law for there methods

3種方法在時(shí)間區(qū)間t∈[0,60]上的保角守恒律誤差，如圖2,3所示.圖2,3采用的參數(shù)與圖1相同.由圖2,3可知：FPEP與CDEP均滿足保角能量守恒律和保角質(zhì)量守恒律,而TAVF卻不能滿足.這表明保角能量守恒方法在保持保角守恒律方面比傳統(tǒng)AVF方法更具優(yōu)勢.

算例2當(dāng)h=0.2,τ=0.001時(shí),在不同阻尼系數(shù)下,FPEP方法在時(shí)間區(qū)間t∈[0,60]上的保角守恒律誤差，如圖4,5所示.由圖4,5可知:在不同類型的阻尼下，F(xiàn)PEP均能精確保持保角守恒律.

圖4 FPEP方法的保角能量守恒律誤差Fig.4 Errors on conformal energy conservation law for FPEP method

圖5 FPEP方法的保角質(zhì)量守恒律誤差Fig.5 Errors on conformal mass conservation law for FPEP method

算例3阻尼系數(shù)α(t)=0.02時(shí),FPEP方法下雙孤立子的碰撞及保角守恒律誤差，如圖6所示.圖6中：初始條件取u(x,0)=12(0.15·0.92sech2(0.9(x+12))+0.015·22sech2(2(x+8)))，并假設(shè)h=0.2,τ=0.001.由圖6(b)可知：雙孤立子碰撞后,除了振幅減小,相空間結(jié)構(gòu)仍保持良好，表明碰撞是彈性的.

(a) 雙孤立子的演化 (b) FPEP在不同時(shí)刻的數(shù)值圖像

(c) 保角能量守恒律誤差 (d) 保角質(zhì)量守恒律誤差圖6 FPEP方法下雙孤立子的碰撞及保角守恒律誤差Fig.6 Collision of two solitons by FPEP and errors on conformal conservation laws

5 結(jié)束語

時(shí)間采用二階AVF方法,空間采用傅里葉擬譜逼近,利用Strang分裂,構(gòu)造變系數(shù)廣義KdV方程的一個(gè)保角能量守恒格式.計(jì)算中使用FFT技巧,不但比空間采用二階中心差分的保角能量守恒方法誤差小,而且計(jì)算效率也很高.數(shù)值實(shí)驗(yàn)結(jié)果表明：相比傳統(tǒng)的AVF保能量方法，采用保角能量守恒方法更能體現(xiàn)保角哈密爾頓系統(tǒng)的本質(zhì)屬性,并具有良好的長時(shí)間數(shù)值行為.由于KdV方程是非線性方程,構(gòu)造阻尼KdV方程在時(shí)間上更為高階的保角保結(jié)構(gòu)算法將是一個(gè)不小的挑戰(zhàn).