二维定常不可压涡量速度NS方程组高精度紧致差分方法.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083292 二维定常不可压涡量-速度 N-S 方程组高 精度紧致差分方法 1 1 2 田巧娴 葛永斌 田振夫 (1.宁夏大学应用数学与力学研究所， 银川 750021；2.复旦大学力学与工程科学系，上海 200433 ) (Email:t12002245129@126.com , Tel 摘要：提出了数值求解二维定常不可压涡量-速度 N-S 方程组高精度紧致差分格式，以及边界条 件的同阶离散格式。针对有解析解的二维定常 N-S 方程组的狄利克雷边值问题和典型的驱动方腔流问 题进行数值模拟，验证了方法的精确性和可靠性。 关键词：不可压 N-S 方程组；涡量-速度方法；高精度紧致格式；有限差分法 0. 引言 不可压 N-S 方程组的数值求解方法受到科研工作者的普遍关注。原始变量法是被普 遍使用的数值方法之一，但是采用这种方法求解，方程组的强耦合性、对流项中非线性 的量以及压力项的处理都给求解过程带来很大的困难。为此，Fasel[1]于 1976 年首次提 出了求解 N-S 方程组的涡量-速度方法，成功地克服了这些困难；Dennis 等人[2]利用该 方法求解了驱动方腔流问题，雷诺数算了 100，并将涡量-速度方法扩展到了求解三维的 稳态流问题中；Farouk 和 Fusegi[3]采用此方法研究了在正方形腔体内和水平环形体中的 自然对流和强迫对流问题；Orlandi[4]结合场方程和边界条件，利用隐式时间步法和交错 网格技术及ADI 迭代方法去解决二维驱动方腔流和后台阶流问题；Cuevirementet 等人[5] 采用有限元方法离散涡量-速度方程，对驱动方腔流问题进行了数值模拟，给出了 Re 数 为 100 和 400 时的三维驱动方腔流问题的结果；Liu[6]将这种方法应用在交错网格上，对 三维驱动方腔的突起和振动两种情况进行了计算；Murugesan 和 Young[7]结合ADI 方法， 对速度 Piosson 方程抛物化后研究了三维驱动方腔问题；另外涡量-速度方法也被用来求 解圆柱绕流问题[8] 。但上述方法均是低精度的，其最高精度均不超过二阶。 本文在已有的求解定常对流扩散方程的高阶紧致格式的基础上，提出了定常不可压 涡量-速度变量的 N-S 方程组的高精度紧致格式，以及边界条件的同阶离散格式。为了 验证本文的高精度紧致差分格式的精确性和可靠性，我们针对有解析解的二维定常 N-S 方程组的狄利克雷边值问题和典型的驱动方腔流问题进行数值实验。 1. 四阶紧致差分格式 二维定常不可压涡量-速度形式的 N-S 方程组可表示为 1 uwx +vwx (wxx =+wyy ) (1) Re u +u =−w (2) xx yy y v +v w (3) xx yy x 其中，u,v 分别是x ,y 方向上的速度分量，Re 为雷诺数。 将二维计算区域在直角坐标系中用一族平行于坐标轴的直线族进行等距网格剖分。 记 x ,y 方向的步长均为 h ，为了方便起见，我们用 0,1,2,3,4,5,6,7,8 分别代表各网格点 (x , y ) 以及与其相邻的 8 个网格点 (x , y ) , (x , y ) , (x , y ) , (x