二维对流扩散方程非均匀网格上的高精度多重网格方法.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083294 二维对流扩散方程非均匀网格上的高精度 多重网格方法 曹富军 1 葛永斌 1 马廷福 2 (1 宁夏大学 应用数学与力学研究所 750021 ； 2 中国矿业大学 银川学院 750021) (Tel:Email: cfj503@163.com) 摘 要：基于非均匀网格上的四阶紧致差分格式并结合多重网格方法求解定常二维对流扩散方程.为了 有效实施多重网格方法，采用面积率构造插值算子和限制算子．通过求解具有边界层的数值算例，验 证了本文基于非均匀网格多重网格算法的精确性和有效性． 关键词：非均匀网格；高阶紧致格式；多重网格方法；插值算子；限制算子 0 引言 高精度紧致差分格式由于精度高并且具有小的离散子域、稳定性好、边界条件容易 处理等特点越来越受到人们的重视．已经发展了针对泊松方程、对流扩散方程以及涡量 流函数变量 Navier-Stokes 方程组的高精度紧致差分格式．大多数格式都是在均匀网格上 实施的．在实际流场计算中, 经常遇到待求物理量如浓度、温度等急剧变化或空间分布 不均匀、大梯度问题、边界层问题或者局部奇性等问题, 使均匀网格格式的计算精度受 到很大影响．比较合理的做法是在大梯度或者边界层区域内多分布些计算点，而在小梯 度或者物理量变化比较平缓的区域内少分布计算节点，这样既可以兼顾算法的稳定性和 计算结果的精确性，又可节省计算量[1] ．因此发展稳定高效的非均匀网格高精度紧致差 分格式已成为有限差分方法发展的一大趋势． Spotz 和 Carey[2]采用坐标变换推导出了一维和二维对流扩散方程非均匀网格上高阶 紧致格式；Ge 和 Zhang[3]将该方法推广，给出了二维对流扩散方程的四阶紧致格式并结 合多重网格方法求解边界层问题；Wang 等人[4]对坐标变换后的涡量-流函数变量的二维 不可压 N-S 方程采用 Maple 数学软件推导了四阶紧致格式；王烜和杨志峰[5]采用变量代 换的方法并且将基于均匀网格提出的优化差分法和反演差分法推广到非均匀网格中, 提 出了一种有效求解定常及非定常非线性对流扩散问题的高精度差分格式. 由于该方法采 用了坐标变换，通过可逆的坐标变换把物理区域上的非均匀网格变成计算区域上的均匀 网格，因此计算仍然是在均匀网格上实施的．该方法在具体问题求解中获得成功，但是 也有缺点．首先，坐标变换往往会导致控制方程会增加一些项，如交叉导数项，这无疑 _________________________ 资助项目：国家自然科学基金资助项目(No 会增加其计算复杂度．其次，如果变换不能显式给出而是通过求解某些微分方程得到， 就会增加计算量．因此，发展不通过坐标变换而直接在求解区域上采用非均匀网格的计 算方法具有重要意义．Kalita 等人[6]针对定常二维对流扩散方程，不采用坐标变换而直 接在物理区域上建立了基于非均匀网格的四阶紧致格式，数值结果充分显示了该方法在 提高计算精度和减少计算量方面的优势． 另一方面，多重网格方法已经从理论上被证明至少对于求解线性椭圆型问题是一种 最优化的数值计算方法．近几年，Gupta 和 Zhang 等人[7-9]在将高精度紧致差分格式和多 重网格方法结合使用的理论分析和数值分析方面做了很好的工作．葛永斌等人[10-12]将这 方面的工作推广到非椭圆方程以及 N-S 方程组的研究中，做出了一些有益的尝试．由于 差分格式往往是建立在均匀网格上的，所以其多重网格算法也是在均匀网格上实施 的．Ge 和 Zhang [3]提出了在非均匀网格上适合边界层对流扩散问题求解的基于高精度紧 致格式的多重网格算法，由于采用了坐标变换，因此多重网格算法实际上是在均匀网格 上实施的．本文利用文献[6]提出的基于非均匀网格上的四阶紧致差分格式，然后，采用 面积率构造适合非均匀网格的多重网格插值算子和限制算子．从而建立其适用的多重网 格算法来求解具有边界层的二维定常对流扩散方程． 1 紧致差分格式 二维定常对流扩散方程为：