BBMB方程的Crank-Nicolson差分格式的一種迭代算法

馬維元

(西北民族大學數學與計算機科學學院，甘肅蘭州，730030)

Benjamin-Bona-Mahony-Burgers(BBMB)方程主要是為了描述小振幅的長波在非線性耗散媒介中的傳播而建立的數學模型，在物理中有著重要的意義，對其進行了大量的研究［1－8］，其初邊值問題具體如下:

其中 α，β，γ，T＞0。

K.Omrani等［8］提出了BBMB方程的一種Crank-Nicolson差分格式，但是其求解是通過經典的解非線性方程組的Newton迭代法進行的。眾所周知，解非線性方程組的Newton迭代法其求解具有一般性，但是其計算量大，精度也不一定總是很高，因此筆者提出了求解Crank-Nicolson差分格式的一種迭代算法，然后證明了該算法是二階收斂的，最后通過數值例子說明了筆者所提出的算法是有效的。

1 解Crank-Nicolson差分格式的迭代算法

首先回顧一下K.Omrani在文獻［8］提出的Crank-Nicolson差分格式。

通常設J，N為任意正整數，記h=(R－L)/J，τ=T/N分別為空間步長和時間步長。定義空間:

對于v，w∈W，為了方便起見，引入以下記號:

另外，定義φ:W×W→W的雙線性函數，具體為(φ(v，w))i=

定義C為廣義常數與步長h和τ無關，且常數C可在不同的情況下取到不同的值。對于BBMB方程提出的Crank-Nicolson格式如下:

該差分格式的截斷誤差為O(h2+τ2)。對于(4)式，可以改寫為:

筆者在文獻［9］的基礎上，提出如下的迭代算法:

其中

2 迭代算法的收斂性

下面證明所提出的迭代算法收斂到差分格式。

定理1 假設方程(1)-(3)的解u(x，t)是充分光滑的，當k和τ充分小時，迭代算法(8)以‖·‖∞范數收斂到Crank-Nicolson差分格式(4)-(6)。

當n=0時，有

當n≥1時，有

從而

當k和τ充分小時，可得

［8］可得

利用(11)式，可得如下的估計式成立

將(12)式帶入(10)式，可得

從而

即

可得

因此，迭代算法(8)以‖·‖1，h范數收斂到Crank-Nicolson差分格式(4)-(6)，再利用離散 Soblev不等式［10］，定理1 得證。……