对流扩散方程的变步长摄动有限差分格式 the perturbational finite difference(pfd) scheme of variational step of convection-diffusion equation.pdfVIP

! ! 辑第#$ 卷第% 期! ! ! ! ! ! 水 动 力 学 研 究 与 进 展 ! ! ! ! ! ! ’() ，*+,) #$ ，-+) % ! ! ! #$$. 年. 月! ! ! ! ! ! ! ! /012-3 04 5672076-89:! ! ! ! ! ! ! ! 8; ，#$$. 文章编号：=$$$?@A? （#$$. ）$%$#B%$A 对流扩散方程的变步长摄动 有限差分格式! =，# % = 李桂波 ，! ! 李明军 ，! ! 高智 （=) 中国科学院力学研究所，高温气体动力学开放研究实验室，北京 =$$$@$ ； # ) 中国海洋大学海洋环境学院，青岛#CC$$% ； %) 湘潭大学数学与计算科学学院，长沙? ===$. ） 摘! 要：! 摄动有限差分（D47 ）方法是构造高精度差分格式的一种新方法。变步长摄动有限差分方法是等步长摄动有限差 分方法的发展和推广。对需要局部加密网格的计算问题，变步长D47 格式不需要对自变量进行数学变换，且和等步长D47 格 式一样，具有如下的共同特点：从变步长一阶迎风格式出发，通过把非微商项（对流系数和源项）作变步长摄动展开，展开幂级数 系数通过消去摄动格式修正微分方程的截断误差项求出，由此获得高精度变步长D47 格式。该格式在一、二和三维情况下分别 仅使用三、五和七个基点，且具有迎风性。文中利用变步长D47 格式对对流扩散反应模型方程，变系数方程及EF(G’(H 方程等进 行了数值模拟，并与一阶迎风和二阶中心格式及其问题的精确解作了比较。数值试验表明，与一阶迎风和二阶中心格式相比，变 步长D47 格式具有精度高，稳定性与收敛性好的特点。变步长D47 格式与等步长D47 格式相比，变步长D47 解在薄边界层型 区域的分辨率得到了明显的提高。 关! 键! 词：! 高精度差分格式；对流扩散方程；变步长摄动有限差分方法 中图分类号：! 0%.%I # ! ! ! 文献标识码： !# $#%’%()*+,)- .*,*# /*..#%#,0#（123 ）40#5# +. 6)%*)*+,)- 4#$ +. 0+,6#0*+,7/*..’4*+, #8’)*+, =，# % = 39 JFKL+ ，! ! 39 8KMGNFM ，! ! J0 OPK （=) 9MHQKQFQ’ +R 8’SP;MKSH ，:PKM’H’ S;T’U +R SK’MS’H ，E’KNKMG =$$$@$ ，:PKM; ； # ) :+,,’G’ +R DPHKS;, ;MT VMWK(+MU’MQ;, 0S’;M+G(;XP ，0S’;M 1MKW’(HKQ +R :PKM; ， YKMGT;+ #CC$$% ，:PKM; ； %) SP++, +R 8;QP’U;QKSH ;MT :+UXFQ;QK+M;, SK’MS’ ，ZK;MGQ;M 1MKW’(HKQ ， :P;MGHP; ? ===$. ，:PKM; ） 9(4%)0 ：! [P’ X’(QF(L;QK+M;, RKMKQ’ TKRR’(’MS’ U’QP+T （D47 ）KH ; M’\ U’QP+T +R S+MHQ(FSQKMG PKGP ;SSF(;Q’ TKRR’(’MS’ HSP’U’) [P’ D47 +R W;(K;QK+M;, HQ’X KH