求解波动方程的高精度紧致显式差分格式

发布时间：2017-12-05 19:03

  本文关键词：求解波动方程的高精度紧致显式差分格式

  更多相关文章： 波动方程 Pade逼近 紧致格式 显式差分 稳定性

【摘要】：本文主要针对波动方程初边值问题的有限差分法进行研究。首先,针对一维波动方程,利用泰勒级数展开公式与原方程推导建立了第一时间层上未知函数值的差分格式；之后,利用Pade逼近离散二阶空间导数项在内部节点上的导数值,时间方向采用中心差分公式进行离散,得到了一种时间二阶、空间四阶精度的紧致显式差分格式,其截断误差为O(τ2+厅4);由于上述格式的时间精度与空间精度不相匹配,本文利用截断误差余项修正的方法,对时间离散进行了改进,改进后格式的截断误差为O(τ4+τ2h2+h4),即格式具有整体四阶精度;然后,通过Fourier分析法分析了两种格式的稳定性,前者的稳定性条件为|α|λ≤(2/3),后者的稳定性条件为|α|λ∈[0,1]∪[2,3]。由于本文格式属于显式差分格式,只需进行一次追赶法求解和一次显式递推计算；最后通过数值实验并与文献中的格式进行对比,验证了本文格式的精确性和稳定性。接下来,将上述一维波动方程的两种格式进行推广,得到了两种求解二维波动方程的高精度紧致显式差分格式,一种格式的截断误差仍为O(τ2+h4),另一种格式的截断误差为O(τ4+τ2h2+h4)。利用Fourier分析法分析了两种格式的稳定性,前者的稳定性条件为|α|λ≤1/3,后者的稳定性条件为|α|λ∈[0,√2/2](?)[1,√6/2]。由于本文格式属于显式差分格式,只需进行两次追赶法求解和一次显式递推计算,无需迭代。最后,通过数值实验,将本文格式的计算结果与文献中格式的计算结果进行了对比,验证了本文格式的精确性和稳定性。最后,将二维波动方程中具有四阶精度的格式推广到三维,推广后的格式仍具有整体四阶精度。通过Fourier分析法得出了该格式的稳定性条件为|α|λ∈[0,√3/3](?)[√6/3,1]。由于本文格式属于显式差分格式,只需进行三次追赶法求解和一次显式递推计算,无需迭代。数值实验结果验证了本文格式的精确性和稳定性。
【学位授予单位】：宁夏大学
【学位级别】：硕士
【学位授予年份】：2016
【分类号】：O241.82

【参考文献】

中国期刊全文数据库 前1条

1 田敏;羊丹平;;波动方程的重叠型区域分解并行有限差分算法[J];山东大学学报(理学版);2007年02期

，

本文编号：1255835