馬健軍, 陳愛敏, 郝怡非(. 四川外國語大學(xué) 國別經(jīng)濟(jì)與國際商務(wù)研究中心, 重慶 40003; . 重慶大學(xué) 城市科技學(xué)院, 重慶4067)

?

Camassa-Holm方程的高階局部中心間斷Galerkin有限元法

馬健軍1, 陳愛敏2, 郝怡非1
(1. 四川外國語大學(xué) 國別經(jīng)濟(jì)與國際商務(wù)研究中心, 重慶 400031; 2. 重慶大學(xué) 城市科技學(xué)院, 重慶402167)

發(fā)展一個(gè)求解具有尖波解的Camassa-Holm方程的高階局部中心間斷Galerkin有限元法,該方法首先將Camassa-Holm方程改寫為一個(gè)守恒律方程和一階方程組的耦合系統(tǒng),然后,使用局部中心間斷Galerkin法求解該守恒律和使用有限元法求解一階方程組,數(shù)值算例用來檢驗(yàn)該方法的精度和有效性.

Camassa-Holm方程; 尖波解; 局部中心間斷Galerkin法; 有限元法

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

Camassa-Holm方程(簡稱C-H方程)是一類十分重要而又特別的新型淺水波方程.1981年,C-H方程由Fuchssteniner和Fokas作為具有雙Hamilton結(jié)構(gòu)的例子給出,隨后,R. Camassa等[1]將其作為淺水波方程重新提出,并發(fā)現(xiàn)了具有的一些特殊性質(zhì)一尖峰孤波解和Blow-up解等,由此引發(fā)了人們對(duì)C-H方程的極大興趣.關(guān)于C-H方程的數(shù)學(xué)理論研究( 從方程的角度) 從20世紀(jì)90年代開始,眾多數(shù)學(xué)家在這方面做了很多重要的研究工作,如孤立波的軌道穩(wěn)定性[2-3]、強(qiáng)解的局部適定性[4]、尖峰孤立子形式弱解的結(jié)構(gòu)[5]、C-H方程光滑孤立波是軌道穩(wěn)定的[6-7]、C-H方程的孤立尖波、孤子類解、周期解[8]和行波解[9]等.

本文考慮如下的Camassa-Holm方程

ut-uxxt+3uux=2uxuxx+uuxxx,

(1)

其中u表示流體的自由表面高度,下標(biāo)ut和ux分別表示對(duì)時(shí)間和空間的偏導(dǎo)數(shù).關(guān)于方程(1) 的基本性質(zhì),文獻(xiàn)[1]中有比較全面的討論.該方程是一個(gè)非線性色散的偏微分方程,除了具有光滑解,還允許有尖波解.由于尖波解的極大值點(diǎn)不光滑,在計(jì)算中可能引起高頻的色散誤差,這是數(shù)值求解該方程的一個(gè)難點(diǎn).因此數(shù)值求解C-H方程吸引了很多研究者的注意.例如,有限差分法[10-11]、有限體積法[12]、局部間斷有限元法[13]、中心間斷Galerkin有限元法[14]等都被使用來模擬該方程的尖波解.本文將提出一個(gè)局部中心間斷Galerkin-有限元法來求解該方程.

2 數(shù)值方法

(2)

(3)

其中,Pk(I)是定義在I上的維度不超過k的多項(xiàng)式集合,{tn}n是對(duì)時(shí)間域的剖分,Δtn=tn+1-tn.方程(1)能夠被改寫為一個(gè)一階系統(tǒng)

qt+f(q,u,r)x=0,

(4)

q=u-rx,

(5)

r=ux,

(6)

(7)

(8)

(9)

(10)

(11)

(12)

3 數(shù)值算例

在這一部分,一些數(shù)值算例被使用來檢驗(yàn)該算法的精度和有效性.在計(jì)算中,總是使用均勻的網(wǎng)格剖分,網(wǎng)格大小為△x,時(shí)間步長△tn=0.1△x,參數(shù)θ=1.

例 1 精度測(cè)試.該算例考慮光滑行波u(x,t)=φ(x-ct),其中φ滿足二階常微分方程

(13)

在計(jì)算中,取α=c=3,方程(13)的初始條件為φ(0)=1,φx(0)=0.由此產(chǎn)生了一個(gè)周期行波,周期為6.469 546 036 35.將四階的Runge-Kutta法在100 000個(gè)節(jié)點(diǎn)上計(jì)算出的近似解做為參考解來計(jì)算本文算法的精度和收斂階,如表1所示.計(jì)算結(jié)果顯示了本文提出的算法對(duì)pk(k=1,2)近似具有k+1階精度.

表 1 u在t=0.3時(shí)的L2誤差和收斂階 Table 3 The L2 tolerance and order of convergence while u at t=0.3

例 2 尖波解.該算例考慮了一個(gè)尖波解.初始條件為u(x,0)=φ(x,x0,α,c),其中

φ

其中,α是周期,取α=30,x0=-5和c=1,計(jì)算區(qū)域是[0,α],離散為320個(gè)均勻網(wǎng)格.時(shí)刻t=0,5,7.5,10的解顯示在圖1中.從計(jì)算結(jié)果可以看出,尖波保持得非常好,該計(jì)算結(jié)果與文獻(xiàn)[14]的結(jié)果也吻合得比較好.

例 3 本算例研究2個(gè)尖波的相互作用,初始條件為u(x,0)=φ(x,-5,30,2)+φ(x,5,30,2).計(jì)算區(qū)域[0,30]離散為320個(gè)均勻網(wǎng)格.時(shí)刻t=0,5,12,18的解顯示在圖2中.從計(jì)算結(jié)果可以看出,2個(gè)尖波在相互作用后保持得非常好,除了尖波的位置在后期有一定的偏移,該計(jì)算結(jié)果與文獻(xiàn)[14]的結(jié)果大體上吻合得比較好(圖1～3中的點(diǎn)均為本文數(shù)值解,線均為文獻(xiàn)[14]的結(jié)果).

例 4 本算例研究3個(gè)尖波的相互作用,初始條件為u(x,0)=φ(x,-5,30,2)+φ(x,-3,30,1)+φ(x,-1,30,0.8).計(jì)算區(qū)域[0,30]離散為320個(gè)均勻網(wǎng)格.時(shí)刻t=0,1,2,3,4,6的解顯示在圖3中.從計(jì)算結(jié)果可以看出,3個(gè)尖波在相互作用后保持得非常好,且與文獻(xiàn)[14]的結(jié)果吻合得非常好.

4 結(jié)論

本文提出一個(gè)高階局部中心間斷迦遼金有限元法來求解一維Camassa-Holm方程的,與文獻(xiàn)[14]不同的是,該方法將Camassa-Holm方程改寫為一個(gè)一階方程組.數(shù)值算例檢驗(yàn)了該方法的精度和有效性.

[1] Camassa R, Holm D D. An integrable shallow water equation with peaked solitons[J]. Phys Rev Lett,1993,71:1161-1164.

[2] Constantin A, Strauss W. Stability of peakons[J]. Commun Pure Appl Math,2000,53(5):603-610.

[3] Zhou Y. Stability of solitary waves for a rod equation[J]. Chaos Solitons Fractals,2004,21(4):977-981.

[4] Danchin R. A note on well-posedness for Camassa-Holm equation[J]. J Differ Eqns,2003,192(2) :429-444.

[5] Constantin A, Escher J. Global weak solutions for a shallow water equation[J]. Indiana Univ Math J,1998,47:1527-1545.

[6] Constantin A, Strauss W. Stability of the Camassa-Holm solitons[J]. J Nonlinear Sci,2002,12:415-422.

[7] Constantin A, Molinet L. Orbital stability of solitary waves for a shallow water equation[J]. Phys,2001,D157:75-89.

[8] 李春海,唐生強(qiáng),黃文韜,等. 一類廣義Camassa-Holm方程的孤立尖波、孤子類解和周期解[J]. 四川師范大學(xué)學(xué)報(bào):自然科學(xué)版,2009,32(5):572-575.

[9] 張克磊,唐生強(qiáng),王兆娟. Camassa-Holm-KP方程的行波解分支[J]. 四川師范大學(xué)學(xué)報(bào):自然科學(xué)版,2010,33(1):58-61.

[10] Coclite G M, Karlsen K H, Risebro N H. A convergent finite difference scheme for the Camassa-Holm equation with generalH1initial data[J]. SIAM J Numer Anal,2008,46:1554-1579.

[11] 常紅. Camassa-Holm方程的守恒有限差分格式[J]. 高等學(xué)校計(jì)算數(shù)學(xué)學(xué)報(bào),2012,34(1):78-86.

[12] Artebrant R, Schroll H J. Numerical simulation of Camassa-Holm peakons by adaptive upwinding[J]. Appl Numer Math,2006,56:695-711.

[13] Xu Y, Shu C W. A Local discontinuous Galerkin method for the Camassa-Holm equation[J]. SIAM J Numer Anal,2008,46:1998-2021.

[14] Li M, Chen A. High order central discontinuous Galerkin-finite element methods for the Camassa-Holm equation[J]. Appl Math Comput,2014,227:237-245.

[15] Cockburn B, Shu C W. The Runge-Kutta discontinuous Galerkin method for conservation laws V:multidimensional systems[J]. J Comput Phys,1998,141:199-224.

2010 MSC：65M60; 35J05

(編輯 鄭月蓉)

High Order Local Central Discontinuous Galerkin-Finite Element Methods for the Camassa-Holm Equation

MA Jianjun1, CHEN Aimin2, HAO Yifei1
(1.ResearchCenterforInternationalBusinessandEconomy,SichuanInternationalStudiesUniversity,Chongqing400031;2.CityCollegeofScienceandTechnology,ChongqingUniversity,Chongqing402167)

In this paper, we develop high order local central discontinuous Galerkin-finite element methods for solving the Camassa-Holm equation which supports peakon solutions. In our numerical approach, we first reformulate the Camassa-Holm equation into a conservation law coupled with a system of first order equations. Then we propose a family of high order numerical methods which discretize the conservation law with local central discontinuous Galerkin methods and the system of first order equations with continuous finite element methods. Numerical tests are presented to illustrate the accuracy and validity of the proposed schemes.

Camassa-Holm equation; peakon solutions; local central discontinuous Galerkin methods; finite element methods

2014-08-29

重慶市自然科學(xué)基金一般項(xiàng)目(CSTC2012JJA00005)

馬健軍(1983—),男,講師,主要從事偏微分方程數(shù)值解的研究,E-mail:majianjun425@163.com

O241.82

A

1001-8395(2015)06-0871-04

10.3969/j.issn.1001-8395.2015.06.016