苏明芳， 何 丽， 胡劲松

(西华大学理学院， 成都 610039)

1 引 言

为研究非线性波在传播中的耗散，人们提出了Benjamin-Bona-Mahony (BBM)方程[1]

ut-uxxt+ux-uxx+uux=0，

(x,t)∈(xL,xR)×(0,T]

(1)

文献[2-3]研究了其解的衰减性.文献[4-6]研究了其解的存在唯一性及收敛性.另一方面，其数值解也引起了众多学者的关注[7-13].

本文考虑BBM方程(1)在如下初边值条件

u(x,0)=u0(x)，x∈[xL,xR]

(2)

u(xL,t)=u(xR,t)=0，t∈[0,T]

(3)

下的数值解，其中u0(x)是已知光滑的函数.问题(1)～(3)具有如下守恒量[14]：

(4)

其中Q(0)为仅与初始条件有关的常数.

文献[14]对问题(1)～(3)提出了在空间层具有四阶理论精度的两层非线性差分格式，但数值求解时需要非线性迭代.文献[15]又对问题(1)～(3)提出了在空间层具有四阶理论精度的三层线性差分格式，但三层格式一般都不是自启动的，且在数值求解时需要储存前两层的数据.本文对方程(1)中的非线性项进行线性化离散处理，构造了一个具有二阶理论精度的两层线性化差分格式，合理地模拟了守恒量(4).在不能得到差分解的最大模估计的情况下，本文综合运用数学归纳法和离散泛函分析方法[16]证明了该格式的收敛性和稳定性，并给出数值算例.

2 差分格式及其守恒律

对问题(1)～(3)考虑如下有限差分格式：

j=1,2,…,J-1;n=1,2,…,N-1

(5)

(6)

(7)

定理2.1差分格式(5)～(7)关于以下离散能量是守恒的：

(8)

其中n=1,2,…,N.

证明 将(5)式两端乘以h后对j从1到J-1求和，由边界条件(7)式和分部求和公式[16]有

(9)

又

(10)

由Qn的定义，将(10)式代入(9)式后两端乘以τ，再对n递推可得(8)式.证毕.

3 差分格式的可解性

定理3.2若时间步长τ充分小，则差分格式(5)～ (7)是唯一可解的.

证明 数学归纳法.显然U0是由初值条件(6)式唯一确定的.假设……