《对流扩散问题非均匀网格上的部分半粗化多重网格方法》.pdfVIP

第38卷 第4期 江西师范大学学报 (自然科学版) Vo1．38No．4 2014年7月 JournalofJiangxiNormalUniversity(NaturalScience) Ju1．2014 文章编号：1000-5862(2014)04-0403-06 对流扩散问题非均匀网格上的部分半粗化多重网格方法 曹富军 ，袁冬芳 ，葛永斌 (1．中国科技大学数学科学学院，安徽合肥 230000；2．内蒙古科技大学数理与生物工程学院，内蒙古包头 014010； 3．北京应用物理与计算数学研究所计算物理实验室，北京 100088；4．宁夏大学数学计算机学院，宁夏 银川 750021) 摘要：结合非均匀网格上的HOC格式与部分半粗化的多重网格方法对具有边界层的2维对流扩散问题 进行了求解 ，并基于面积率构造了部分半粗化多重网格方法的插值算子和限制算子．数值实验表明：对于 只需要在 1个方向进行网格加密的边界层问题，基于部分半粗化的网格分布策略及多重网格算法可以大 大减少无边界层方向的网格数，从而较完全粗化的网格分布策略及多重网格算法具有更高的计算精度和 求解效率． 关键词 ：非均匀网格；高精度紧致格式；部分半粗化；多重网格；边界层对流扩散 中图分类号：0241 文献标志码 ：A 不少有意义的工作．文献[2]提出了非等距网格上2 0 引言 维泊松方程的网格部分半粗化策略(网格粗化过程 仅在主控方向上进行，而在非主控方向上的网格步 考虑具有如下形式的2维对流扩散方程： 长则保持不变)，数值实验表明：对于各 向异性问 一 口M 一bu +c(x，Y)M +d(x，Y) = ，Y)， 题，由于减少了非主控方向的网格数，该策略下的多 其中 ．Z／(，Y)是未知量，C(，Y)和d(，y)分别为 重网格算法比网格完全粗化策略下的多重网格算 和Y方向的对流系数 ，，(，y)为计算域上足够光滑 法，在不损失精度的前提下，大大节约了计算时间， 的函数． 提高了计算效率．文献 [9]将部分半粗化策略推广 求解对流扩散方程的方法较多，如传统的中心 到非等距网格上3维泊松方程的求解；文献 [10]基 差分格式、迎风格式等，然而，它们都不能够满足实 于J．C．Kalita等提出的求解2维对流扩散方程的高 际计算所需要的精度．紧致差分格式由于精度高并 精度紧致格式，提出相应的完全粗化的多重网格算 且具有小的离散子域、对单元的敏感性小和对边界 法，并通过对边界层问题的求解，展现了非均匀网格 单元的处理无特殊困难等特点越来越受到人们的重 上多重网格算法的效率．文献 [11]将该方法推广到 视 ¨j．然而，大多数紧致格式是基于均匀网格上提 3维，建立了3维非均匀网格上完全粗化多重网格算 出的，对于求解大梯度或者边界层问题，均匀网格上 法的插值算子和限制算子，并求解了3维泊松方程． 的紧致格式很难得到理论上的高精度数值解．因此， 该方法的缺点是在所有坐标轴方向上采用相同的网 J．C．Kalita等…直接在物理区域上建立了求解2维 格数，这对于边界层或大梯度仅在 1个方向的问题， 定常对流扩散方程非均匀网格的高精度紧致差分格 无疑会浪费计算量，因为在该方向并不需要网格加 式，由于该方法在物理区域上直接离散偏微分方程， 密．基于此，本文将设计一种非均匀网格上的网格部 没有采用坐标变换，因此减少了坐标变换带来的额 分半粗化策略，并建立其相应的多重网格算法，然后 外计算量，并且数值结果