三维对流扩散方程的高精度多重网格方法.pdfVIP

第25卷 第 12期 赤 峰 学 院 学报 (自然 科 学 版 ) V0l_25No．12 2009牟 12月 JournalofChifengUniversity(NaturalScienceEdition) Dee．20o9 三维对流扩散方程的高精度多重网格方法 马廷福 ，金 涛 ，葛永斌 (1．中国矿业大学 银川学院，宁夏 银川 750011；2．宁夏大学 应用数学和力学研究所，宁夏 银川 750021) 摘 要：在立方体网格上建立了数值求解三维变系数对流扩散方程的四阶精度 19点紧致差分格式，采用多重网格加速 技术，建立了适用于本文高精度紧致差分格式的多重网格算法，从而大大加快了传统迭代法的收敛速度．数值实验结果表明 本文方法对于不同的网格雷诺数问题，在准确性、稳定性以及减少计算_T-作量方面均明显优于7点中心差分格式． 关键词：三维对流扩散方程；四阶紧致差分格式；多重网格方法 中图分类号：029 文献标识码：A 文章编号：1673—260X(2009)12～0012—03 对流扩散方程的数值计算在计算流体力学领域对流动 知函数且具有充分的光滑性 问题的数值模拟中扮演着一个非常重要的角色．传统的差分 取计算网格为等距网格，步长用h表示．为了方便起见， 离散格式，普遍有着精度低和计算不稳定的缺陷．如传统的 我们分别用数字0到26依次代表网格点 迎风格式，尽管计算是稳定的，但是只有一阶精度 ，从而导 (x,y，z)，(x+h，Y，z)，(x，y+h，z)，(x—h，)，z)，x，Y—h，z)，(x，Y，z+h)，(x，Y，z—h)， 致计算结果受到很大人工粘性的影响而不够精确．而传统的 +h，y+h，z)(x-h，y+h，z)，x—h．，-h，z)，+b，y-h，z)，fx+b，Y，z+h)， 中心差分格式，尽管有两阶精度，但是当求解问题表现为对 ，y+h，z+h)，(x-h，Y，z+h)，，Y—h，z+h)，(x+h，Y，Z—h)，x，y+h，z—h)， 流占优或者网格雷诺数超过某一数值范围时，计算就会发 x(—h，y，z-h)，x，y-h，z—h)，(x+h，y+h，z+h)，(x～h，y+h，z+h)， 散．因此发展对流扩散方程精度高且稳定性好的差分离散方 (x-h，Y—h，z+h)，(x—h，y-h，z+h)，(x+h，Y—h，z+h)，(x+h，y+h，z-h)， 法是 目前研究的一个热点问题． (x-h，y+h，z-h)，(x-h，y-h，z-h)，(x+h，Y—h，z—h)． 对于一维和二维的对流扩散方程此方面的研究已经有 不妨记gij=g，一gJ，g可代表 ，p，q，r,f等．首先 ，将方程(1)改 很多的报道II_，各种新的差分格式在计算精度和稳定性方 写为： 面都有了显著的改善和提高．但实际问题往往是三维的，其 中、、+ +中 g(x，Y，z) (2) 更高的代数复杂性要求使得对三维对流扩散方程的高精度 g(x，Y，z)=f(x，Y，z)一p(x，y，z)中-q(x，Y，z)中一r(x，Y，z)O (3) 差分方法的研究和发展远不如对一维和二维充分．文献!针 fx，Y，z)点处的二阶导数项利用如下四阶紧致差分公式 ： 对一般三维变系数椭圆型方程提出了一种发展其四阶精度 中=[(1十_-8刁6 0+0h(4) 4() 紧致差分格式的思想和方法，但是由于最终给定的公式过 于抽象和复杂而使得人们必须花很多的时间和精力才能由 = (【1+睾_山88中。+0(113 (5)