求解一维扩散方程的一种高精度紧致差分方法A High-order Compact Difference Method for the One-dimensional Diffusion Equation
杨晓佳,葛永斌
摘要(Abstract):
针对一维扩散方程,空间采用四阶Padé公式,时间采用广义的梯形公式,差分离散得到了一种时间二阶、空间四阶精度的隐式紧致差分格式,其截断误差为O(τ2+h4).通过理论分析证明了此格式是无条件稳定的.最后通过数值实验验证了格式的精确性和可靠性.
关键词(KeyWords): 扩散方程;Padé逼近;紧致格式;广义梯形公式;稳定性
基金项目(Foundation): 国家自然科学基金资助项目(11361045;11161036);; 宁夏高等学校科学技术研究项目(NGY2013019)
作者(Author): 杨晓佳,葛永斌
参考文献(References):
- [1]CHAWLA M M,Al-ZANAIDI M A,EVANS D J.Generalized trapezoidal formulas for parabolic equations[J].International journal of computer mathematics,1999,70(30):429—443.
- [2]SALLAM S,NAIM A M,ABDEL-AZIZ M R.Unconditionally stable C1-cubic spline collocation method for solving parabolic equations[J].International journal of computer mathematics,2004,81(7):813—821.
- [3]周顺兴.解抛物型偏微分方程的高精度差分格式[J].计算数学,1982,4(2):204—213.
- [4]马明书.一维抛物型方程的一个新的高精度显式差分格式[J].数值计算与计算机应用,2001,22(2):156—160.
- [5]葛永斌,田振夫,詹咏,等.求解扩散方程的一种高精度隐式差分方法[J].上海理工大学学报,2005,27(2):107—112.
- [6]詹涌强.求解抛物型方程的一族六点隐式差分格式[J].安徽大学学报(自然科学版),2012,36(4):26—29.
- [7]WEN R H,SHAO H Z.Domain decomposition schemes with high-order accuracy and unconditional stability[J].Applied mathematics and computation,2013,219(11):6170—6181.
- [8]开依沙尔热合曼,努尔买买提黑力力.求解扩散方程的二级四阶隐式Runge-Kutta方法[J].湖北大学学报(自然科学版),2014,36(5):476—480.
- [9]詹涌强,张传林.解抛物型方程的一族高精度隐式差分格式[J].应用数学和力学,2014,35(7):790—797.
- [10]周敏,高学军,董超.解抛物型方程的八点隐式差分格式[J].广东工业大学学报(自然科学版),2014,31(4):71—78.
- [11]LELE S K.Compact finite difference schemes with spectral-like resolution[J].Journal of computational physics,1992,103(1):16—42.
- [12]CHAWLA M M,Al-ZANAIDI M A,EVANS D J.A class of generalized trapezoidal formulas for the numercal integration of y'=f(x,y)[J].International journal of computer mathematics,1996,62(1/2):131—142.