二维波动方程的高精度交替方向隐式方法.pdfVIP

二维波动方程的高精度交替方向隐式方法

马月珍;李小纲;葛永斌

【摘要】基于二阶微商的四阶紧致差商逼近公式及加权平均思想,提出了数值求解

二维波动方程的2种精度分别为O(τ~2+h~4)和O(τ~4+h~4)的交替方向隐式

(ADI)格式,以及与其相匹配的第一个时间层的同阶离散格式,并且通过Fourier方法

分析了格式的稳定性.该方法在沿每个空间方向上只涉及3个网格基架点,因此可以

重复采用TDMA算法,从而大大节省计算时间.数值实验验证了所用方法的精确性和

可靠性.

【期刊名称】《四川师范大学学报（自然科学版）》

【年(卷),期】2010(033)002

【总页数】5页(P179-183)

【关键词】波动方程;高阶紧致格式;交替方向隐式方法;稳定性

【作者】马月珍;李小纲;葛永斌

【作者单位】宁夏大学,应用数学与力学研究所,宁夏,银川,750021;宁夏大学,应用数

学与力学研究所,宁夏,银川,750021;宁夏大学,应用数学与力学研究所,宁夏,银

川,750021

【正文语种】中文

【中图分类】O241.82

有限差分法[1-3]是数值求解偏微分方程的常用方法之一.在航空、气象、海洋、水

利等许多流体力学的问题中,常常遇到双曲型偏微分方程,针对传统的差分离散格式,

普遍有着精度低且受到很强的稳定性条件限制的缺陷,因此发展其高精度且稳定性

好的差分离散方法具有十分重要的意义[4-7].

交替方向隐式(ADI)方法是数值求解该类问题的一种非常有效的方法,它将高维问题

转化为若干一维问题进行求解,而一维问题又可采用高效的TDMA算法,从而可以

大大提高计算效率,节省存储空间.传统的ADI方法如Peaceman-Rachford(P-R)

格式和Beam-Warming[8]格式都是二阶精度的.另一方面,文献[9]提出了求解二

维非定常对流扩散方程的高精度ADI格式,文献[10]提出了求解高维热传导方程的

高精度ADI格式.本文依然根据ADI方法的基本思想,提出两种求解二维波动方程

的高精度紧致ADI差分格式,为此考虑初边值问题：

其中Ω={(x,y):0≤x,y≤1},为ΩΓ的边界,u(x,y,t)为待求未知量,f(x,y,t)为源

项,φ(x,y),ψ(x,y)和g(x,y,t)均为已知函数且具有充分的光滑性.

1高阶紧致ADI格式

用τ表示时间步长,空间取等间距网格,步长用h表示.网格点为

(xi,yj,tn),xi=ih,yj=jh,tn=nτ,i,j=0,1,…,N,h=1/N,n≥0.

1.1ADI(2,4)格式考虑(1)式在n时刻值,对时间导数项采用中心差分,空间导数项

采用Kreiss[11]提出的四阶紧致差分公式:

则有

其中

在空间方向以和的算术平均值代替可得

对(5)式进行整理且略去高阶项可得

为了构造ADI差分格式,采用与文献[9-10]类似的技巧,在(6)式左端加上可得

显然,(7)式与(6)式的截断误差同阶,利用

可将(7)式写为如下形式

引入一个过渡变量则可将(8)式写为

对于过渡变量的边界条件,可以由下式给出

(9)式即为求解二维波动方程的高精度紧致ADI格式,其精度为O(τ2+h4),记为AD

I(2,4).

1.2ADI(4,4)格式对时间和空间导数项均采用四阶紧致差分公式,考虑(1)式在n时

刻值,可得

对上式进行整理且略去高阶项可得

可将(10)式写为如下形式

为了构造ADI差分格式,采用与文献[9-10]类似的技巧,在(11)式左端加上并进行

因式分解,可得

显然,(11)式与(12)式的截断误差同阶,引入一个过渡变量则(12)式可写为

对于过渡变量的边界条件,可以由下式给出

(13)式即为求解二维波动方程的高精度紧致ADI格式,其精度为O(τ4+h4),记为

ADI(4,4).

1.3初始条件的离散因为格式是3层的.即每一次时间推进都需要知道前两个时间

步的值,初始时刻有(2)式精确给出,第一个时间步的值由(3)式给出,因此,须对(3)式

进行离散,下面推出与(9)和(13)式相匹配的第一个时间步的离散格式.

利用Taylor展开式将在处展开可得

利用(1)～(3)式,且略去高阶项,即可得与(9)式相匹配的第一个时间步的离散格式

与(13)式相匹配的第一个