基于三维对流扩散方程四阶紧致差分格式的预条件迭代法.pdfVIP

第33卷 第1期 宁夏大学学报 (自然科学版) 2012年3月 Vo1．33No．1 JournalofNingxiaUniversity(NaturalScienceEdition) M ar．2012 文章编号 ：0253—2328(2012)01—0001—05 DOI：CNKI：64—10O6／N1731．004 基于三维对流扩散方程四阶紧致差分格式的 预条件迭代法 袁冬芳，庄 昕，葛永斌 (宁夏大学 应用数学与力学研究所 ，宁夏 银川 750021) 摘 要 ：针对三维对流扩散方程 ，采用四阶紧致差分格式和预条件迭代法进行数值 实验 ，利用带填补数 的不完全 LU 分解 (ILUT (r，s))做预处理器 ，FGMRES(20)做迭代加速器对离散所得方程组进行求解．验证 了四阶紧致差 分格式的计算精度 ，通过比较预条件迭代法与高斯一赛德尔迭代法 以及超松弛迭代法的迭代次数和 CPU时间，充 分显示 了预条件迭代法的高速求解特性． 关键词 ：三维对流扩散方程 ；四阶紧致格式 ；迭代法 ；预条件技术；不完全 LU分解 分类号 ：(中图)O24l_82 (200MR)65M06；65M12 文献标志码 ：A 对于三维问题 的数值模拟，由于对计算机存储 的输运规律．当对流项系数较大时，方程 (1)被认为 空间和算法的求解效率都有更高的要求，因此受到 是对流 占优的，否则是扩散 占优的． 了研究者更多的关注。由于传统 的计算格式精度低 ， 对二维的情形，文献[1—2]已提出了对流扩散 要想得到精确的计算结果 ，需要将网格剖分得很细， 方程和Navier—Stokes方程组的一些四阶紧致格式． 从而离散所得线性方程组的阶数也较高 ，以致于在 这些格式有很好 的数值稳定性 ，并且能获得高精度 现代计算机上仍然难以实现 ，甚至是不可能的．解决 的近似解．文献[3～4]利用多重 网格方法求解 四阶 这一问题的方法之一是利用高精度或谱方法 ，它们 格式 ，得到了很好的数值结果．对于扩散 占优问题， 通常可采用较粗 的网格离散得到较高 的精度 ，由于 利用多重网格方法求解 ，四阶紧致格式 比二阶中心 所得线性方程组 的阶数相对较低 ，因此可以大大减 差分格式更高效_5]．文献[6]构造了非均匀网格上求 少计算量，节约计算 时间．为此 ，本文考虑三维对流 解定常对流扩散方程的高精度紧致差分格式 ，数值 扩散方程的狄利克雷边值问题 ： 验证了格式的高精度和高分辨率的优点．为了得到 所需精度的计算结果，四阶紧致格式 比中心差分格 粤+十磬+十馨+十pC,x,Y，’2)五+十 式更省时间，而且只占用较小的内存． q(x，Y，z)Du- +-r(x，，z)Du— f(x 对三维对流扩散方程 (1)，最常用的格式是七点 ，Y，z)， (1) 中心差分格式和迎风格式．格式通常采用传统 的迭 (z，Y，z)∈n， 代法进行求解，如Jacobi和Gauss—Seidel迭代法等． u(x，Y，)一 g(x，Y，)， (z，Y，)∈a0， 然而求解对流 占优问题时，传统迭代法可能不收敛． 其 中a力为n的边界．这里假定 是 R上的单位立 对一阶迎风格式 ，传统迭代法尽管收敛，但得到的解 方体 区域，未知函数 u(x，，)、对流项系数 P(z，， 精度太低．Greif和Varah对三维常系数对流扩散方 )、q(x，Y，)及