盛秀蘭, 艾 堯, 吳宏偉

(1.東南大學(xué)數(shù)學(xué)系 江蘇南京210096;2.江蘇廣播電視大學(xué)公共課教學(xué)部 江蘇南京210036)

一個(gè)類(lèi)似Burgers方程的數(shù)值解

盛秀蘭1,2, 艾 堯1, 吳宏偉1

(1.東南大學(xué)數(shù)學(xué)系 江蘇南京210096;2.江蘇廣播電視大學(xué)公共課教學(xué)部 江蘇南京210036)

研究了一個(gè)類(lèi)似Burgers方程的初邊值問(wèn)題的有限差分方法.基于Crank-Nicolson方法,建立了一個(gè)兩層線(xiàn)性化隱式差分格式,討論了差分格式的可解性.利用離散能量估計(jì)方法證明了差分解在最大模意義下關(guān)于時(shí)間和空間的二階收斂性,并用數(shù)值算例驗(yàn)證了理論分析結(jié)果.

類(lèi)似Burgers方程;隱式差分格式;收斂性;離散能量估計(jì)

0 引言

考慮初邊值問(wèn)題[1]:

其中,p=1,2,d>0為粘性常數(shù),f(u)是非線(xiàn)性函數(shù),φ(x),α(t),β(t)為已知光滑函數(shù),且滿(mǎn)足相容性條件.

當(dāng)p=1且f(u)≡0時(shí),方程(1)成為著名的Burgers方程,有很多文獻(xiàn)研究了該問(wèn)題的數(shù)值解法.可以通過(guò)Hop f-Cole變換將方程化為一個(gè)標(biāo)準(zhǔn)的熱方程,然后用高階差分格式求解[2-4].但是當(dāng)f(u)≠0時(shí),此方法不再適用.文[5]研究了當(dāng)f(u)≡0時(shí)保持能量守恒的差分格式,但當(dāng)f(u)非零時(shí)能量不再守恒.文[6]給出了一個(gè)線(xiàn)性化的差分格式,該格式的截?cái)嗾`差是O(τ+h2),但文中沒(méi)有給出理論分析結(jié)果.

通過(guò)分析以上文獻(xiàn),正如[3]指出的,當(dāng)p=1且f(u)=0,如果對(duì)方程(1)直接用差分逼近,得到的線(xiàn)性化格式一般關(guān)于時(shí)間的精度是1階的,要得到關(guān)于時(shí)間是2階的格式,如用Crank-Nicolson格式,通常格式是非線(xiàn)性的或者是3層以上格式.其原因主要是非線(xiàn)性項(xiàng)uux的線(xiàn)性化離散問(wèn)題.當(dāng)帶有非線(xiàn)性強(qiáng)迫項(xiàng)后,要構(gòu)造關(guān)于時(shí)間2階收斂的線(xiàn)性化兩層格式顯得更為困難.

本文研究當(dāng)p=2時(shí)初邊值問(wèn)題(1)～(3)的線(xiàn)性化差分解法.構(gòu)造了一個(gè)關(guān)于時(shí)間和空間2階收斂的線(xiàn)性化兩層差分格式.在處理非線(xiàn)性項(xiàng)f(u)的時(shí)候,采用與文獻(xiàn)[11]類(lèi)似的方法,而對(duì)非線(xiàn)性項(xiàng)u2ux的逼近,提出了一種新的離散方法.證明了差分格式解的存在唯一性和收斂性.給出了差分格式的解在離散最大模意義下收斂階數(shù)為O(τ2+h2).最后給出了兩個(gè)數(shù)值算例,計(jì)算結(jié)果顯示本文提出的方法是有效的.

本文假設(shè):

H1問(wèn)題(1)～(3)存在唯一光滑解,且存在常數(shù)C0,使得對(duì)有

H2函數(shù)f二階可導(dǎo),且存在正常數(shù)C1,δ,使得當(dāng),有

1 記號(hào)及差分格式

取正整數(shù)m,n,記h=1/m,τ=T/n,xi=ih,0≤i≤m,tk=τk,0≤k≤n.定義Ωh={xi|0≤i≤={(xi,tk)|0≤i≤m,0≤k≤n}.設(shè){|0≤i≤m,0≤k≤n}為Ωhτ上的網(wǎng)格函數(shù),引進(jìn)記號(hào):

我們構(gòu)造初邊值問(wèn)題(1)～(3)的差分格式為:

差分格式(4)～(6)是一個(gè)兩層線(xiàn)性化的差分格式,對(duì)每一固定的時(shí)間層k,差分方程的系數(shù)矩陣是3對(duì)角矩陣,可以用追趕法求解.

2 差分格式的可解性和收斂性

引理1[12-13]設(shè)v,w∈Vh,則有

下面的引理2是離散非線(xiàn)性項(xiàng)u2ux的關(guān)鍵.

引理2設(shè)函數(shù)u(t),v(t),w(t)∈C2[0,T],記.則

證明記p(t)=(uvw)(t),則由Taylor展開(kāi)及向前差商得

引理2證畢.

利用引理1和引理2以及能量估計(jì)和數(shù)學(xué)歸納法,可以得到下面的收斂性估計(jì).

定理1設(shè)u(xi,tk)是問(wèn)題(1)～(3)的解,滿(mǎn)足差分方程(4)～(6),記,如果φ(x)∈ C[0,1],且存在常數(shù)C0>0,使得,則當(dāng)h,τ充分小且時(shí),有

其中,C>0,K>0,S>0是與h,τ無(wú)關(guān)的常數(shù).

由定理1可知,當(dāng)步長(zhǎng)h,τ充分小時(shí),對(duì)固定的k,差分格式(4)～(6)的系數(shù)矩陣嚴(yán)格對(duì)角占優(yōu),因此有下面的結(jié)論.

定理2當(dāng)步長(zhǎng)h,τ充分小,差分格式(4)～(6)存在唯一解.

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

利用差分格式(4)～(6)計(jì)算2個(gè)實(shí)例,計(jì)算結(jié)果見(jiàn)表1和表2.

例1

該問(wèn)題的精確解為ex+1.

例2

表1 例1的不同步長(zhǎng)下的誤差和收斂階數(shù)Tab.1 Errorsandconvergencerateaboutthefirstexample

表2 例2的不同步長(zhǎng)下的誤差和收斂階數(shù)Tab.2 Errorsandconvergencerateaboutthesecondexample

[1] Bülent S,˙Idris D.A numerical study of the Burgers’equation[J].J Franklin Institute,2008,345(4):328-348.

[2] Brander O,Hedenfalk J.A new formulation of the general solution to Burgers’equation[J].Wave Motion,1998,28 (4):319-332.

[3] Liao Wenyuan.An imp licit fourth-order compact finite difference scheme fo r one-dimensional Burgers’equation[J].App l Math Comput,2008,206(2):755-764.

[4] Yousuf M.On the classof high order time stepping scheme based on Padé’s app roximations for the numerical solution of Burgers’equation[J].App l Math Comput,2008,205(1):442-453.

[5] Anguelov R,Djoko J K,Lubuma JM S.Energy p roperties p reserving schemes fo r Burgers equation[J].Numer M ethods Partial Differential Equations,2008,24(1):41-59.

[6] Kutluay S,Esen A.A linearized numerical scheme for Burgers-like equations[J].App l Math Comput,2004,156(2): 295-305.

[7] Djoko J K.On the long-time stability of a backward Euler scheme fo r Burgers’equation w ith polynomial force[J].Numer Methods Partial Differential Equations,2008,24(6):1371-1387.

[8] Lewandow ski J L V.A markermethod for the solution of the damped Burgers’equation[J].Numer Methods Partial Differential Equations,2006,22(1):48-68.

[9] 孔令江,張超英,譚惠麗,等.用格子Boltzmann方法模擬kdv-Burgers方程的激波解[J].廣西師范大學(xué)學(xué)報(bào):自然科學(xué)版,2003,21(4):1-4.

[10] 洪寶劍,盧殿臣,張大珩.帶強(qiáng)迫項(xiàng)變系數(shù)組合kdv-Burgers方程的顯式精確解[J].廣西師范大學(xué)學(xué)報(bào):自然科學(xué)版, 2007,25(1):17-20.

[11] 吳宏偉.二維半線(xiàn)性反應(yīng)擴(kuò)散方程的交替方向隱格式[J].計(jì)算數(shù)學(xué),2008,30(4):349-360.

[12] 孫志忠.偏微分方程數(shù)值解[M].北京:科學(xué)出版社,2004.

[13] Zhou Yulin.App lication of Discrete Functional Analysis to the Finite Difference Method[M].Beijing:International Academic Publishers,1990.

Numerical Solution of a Burgers-like Equation

SHENG Xiu-lan1,2, A I Yao1, WU Hong-w ei1
(1.Department of M athem atics,Southeast University,N anjing 210096,China; 2.D ivision of Public Course,Jiangsu Radio and TV University,N anjing 210036,China)

A second-order finite difference scheme for a Burgers-like equation is p roposed and analyzed.The scheme constructed is a double-deck linearized scheme that can be solved easily by Thomas’algo rithm.The solvability of difference scheme is obtained.It is show n that the difference scheme is convergent in maxim um no rm.The convergence o rder is second-o rder in bo th space and time.Some numerical experiments are conducted to illustrate the theoretical results of the p resented method.

Burgers-like equation;imp licit difference scheme;convergence;discrete energy estimate

O 241.82

A

1671-6841(2010)03-0023-04

2009-11-25

國(guó)家自然科學(xué)基金資助項(xiàng)目,編號(hào)10871044.

盛秀蘭(1976-),女,講師,主要從事微分方程數(shù)值解研究,E-mail:113525336@qq.com.