二维扩散反应方程的三次样条高精度加权隐式差分格式及其多重网格方法  

High-precision weighted implicit difference format of two-dimensional diffusion-reaction equation with cubic spline and its multigrid method

作　　者：马廷福[1] 金涛[1] 葛永斌[2] 

机构地区：[1]中国矿业大学银川学院,宁夏银川750011 [2]宁夏大学应用数学与力学研究所,宁夏银川750021

出　　处：《兰州理工大学学报》2013年第3期141-146,共6页Journal of Lanzhou University of Technology

基　　金：国家自然科学基金(11061025);霍英东教育基金会高等院校青年教师基金(121105);宁夏高等学校科学技术研究项目的资助

摘　　要：利用二阶微商的三次样条四阶紧致差分逼近公式,推导出两种数值求解二维扩散反应方程的两层9点加权隐式紧致差分格式.当θ=1/2时,该格式在时间和空间方向上分别达到二阶和四阶精度.通过Fourier方法讨论知,当1/2≤θ≤1时,格式是无条件稳定的;当0≤θ<1/2时,格式是条件稳定的.为了克服传统迭代法在求解隐格式方面的困难,差分方程采用多重网格方法进行求解并将本文格式的结果与P-R格式及C-N格式下的结果进行比较.数值实验结果验证本文方法的精确性和可靠性及多重网格方法的效率.Applying second derivative approximation formula with cubic spline and fourth-order compact difference, two kinds of weighted implicit compact difference formats with two levels and nine points were derived for solving two-dimensional diffusion-reaction equation. The two formats could achieve second-or- der accuracy in time and fourth-order accuracy in space respectively when θ= 1/2. By means of the Fourier analysis, it could be concluded that when 1/2≤θ≤1, they would be unconditionally stable. When 0≤θ〈1/2, they would be conditional stable. In order to overcome the ciifficulty of conventional method in sol- ving the implicit format, the multigrid method was adopted to acc1crate convergence rate. The numerical result from two formats presented were compared with that from P-R and C-N formats. The comparison result verified the accuracy and stability of two schemes presented and the efficiency of multigrid method.

关 键 词：扩散反应方程 加权隐式差分格式 高精度 多重网格方法 

分 类 号：O241.82[理学—计算数学]