广义sine-Gordon方程高精度隐式紧致差分方法.pdf

2015年 5月 计 算 数 学 第 37卷第 2期 Mav．2015 MATHEMATICA NUMERICA SINICA Vo1．37，No．2 广义 sine—Gordon方程高精度 隐式紧致差分方法冰 耿晓月 刘小华 (西南石油大学理学院，成都637001) 摘 要 本文研究一类二维非线性的广义 sine—Gordon(简称SG)方程的有限差分格式．首先构造三 层时间的紧致交替方向隐式差分格式，并用能量分析法证明格式具有二阶时间精度和四阶空间精 度．然后应用改进的Richardson外推算法将时间精度提高到四阶．最后，数值算例证实改进后的 算法在空间和时间上均达到四阶精度． 关键词：SG方程；紧致差分格式；交替方向隐格式；外推法；能量分析法 MR (2000)主题分类：65M06，65N06 1．引 言 SG方程：是数学物理中的重要的非线性发展方程，在超导理论，混沌吸引子，正确 曲率面， 晶格位错的传播等方面有广泛应用 [1--4]，研究该方程的数值方法具有非常重要的理论意义和 实用价值．上世纪90年代末，方樾等人 []提出SG方程的蛙跳差分格式，讨论了三对角循 环线性微分方程 的可解性，并用能量不等式法证明该格式具有二阶精度．Bratsos等分 别运用线性方法 [。]和参数化差分法 []研究 SG方程，这两种方法均有二阶精度；2005年， Bratsos提出2+1维 SG方程的三层时间显式差分格式 ISJ，该方法有二阶精度；2008年，他将 这种方法推广到一维 sG方程，空间精度为二阶，而时间精度从 [7]中的二阶提高到三阶 _9j． 随后用预估格式改进 f9]中的方法，空间精度为二阶，而时间精度从三阶提高到四阶 ．许秋 滨等人 [11]提出了SG方程两个均为二阶精度的隐式差分格式，用离散泛函分析的方法证明 这两种格式的稳定性，并证明了追赶迭代法求解格式的收敛性．胡劲松等人 _12_对文献 [11】中 第一种格式进行修正，在保持二阶精度情况下，提出修正的隐式差分格式，减少了迭代计算量， 并在新格式基础上运用 Richardson外推法，时间精度为二阶，而空间精度达到四阶． 紧致差分格式通过较少计算节点得到高精度的数值解，有效减少计算量和存取量，众多学 者开始把紧致差分方法应用于求解 SG方程中．周艳祖等人 _13_建立了一维 SG方程的三种 差分格式，包括二阶时间精度和四阶空间精度的三层显格式和三层隐格式，以及时间和空间精 度均为四阶的三层紧致格式．Mohebbi等 _14_在空间和时间上分别用紧致差分法和对角隐式 法，得到四阶精度的差分格式．Cui[]用四阶紧致差分算子离散二阶空间导数，用交替方向的 线性算法求解非线性二维 SG方程，空间精度为四阶，时间精度为二阶．Sari等 bj在空间方 向用六阶紧致差分格式，时间方向用三阶 Runge—CKutta格式求解一维SG方程，空间和时间 精度分别达到六阶和三阶． 2014年 9月 18日收到 200 计 算 数 学 2015在 上述文献 中一维 SG方程 已具有较高精度，而二维 SG方程的初边值问题研究还较少，文 献 [8]精度仅为二阶，文献 1【5]空间精度为四阶，而时间精度仅为二阶．本文在 1【5]基础上给 出时间和空间均为四阶的差分格式．首先将 SO方程离散化，用四阶紧致差分格式替代二阶空 间导数；其次加入扰动项，并引入过渡层，将方程化为 Douglas—Gunn交替方向隐格式 (ADI)， 此时方程有四阶空间精度和二阶时间精度；再次通过能量分析法研究格式的收敛性．由于差分 格式的第一层截断误差中存在时间步长的奇次幂，基于两层网格的Richardson外推法无法使 差分格式达到四阶精度，为了加快时间的收敛速度，本文对方程的初始条件进行处理，并运用 改进后基于三层格式的Richards