封闭方腔内自然对流问题的高精度紧致差分格式.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083293 封闭方腔内自然对流问题的 高精度紧致差分格式 1 1 1 2 金 涛 葛永斌 马月珍 田振夫 (1. 宁夏大学应用数学与力学研究所，银川 750021 ； 2. 复旦大学力学与工程科学系，上海 200433) (Tel: Email: yuezhenm@163.com) 摘要：提出了数值求解二维非定常不可压涡量-流函数 Navier-Stokes/Boussinesq 方程组的高精度紧 致差分格式,格式空间为四阶精度,时间为二阶精度,并且是无条件稳定的.为了验证本文高精度紧致差 分格式的精确性和可靠性,对有解析解的二维非定常不可压 Navier-Stokes/Boussinesq 方程组的 Dirichlet 问题和典型的封闭方腔自然对流问题进行了数值模拟. 关键词：不可压Navier-Stokes/Boussinesq 方程组,涡量-流函数方法,高阶紧致差分格式,自然对流 0 前言 二十世纪八九十年代,高精度紧致差分格式由于具有使用节点少、精度高且使求解 问题的边界处理简单等优点倍受人们关注, 已越来越被广泛地应用于不可压Navier- Stokes方程组的数值求解中. Dennis和Hudson[1]提出了一种九点模板的四阶紧致格式求 解了腔体左右壁面有温差的自然对流问题,瑞利数计算到了105 ；Gupta[2]利用四阶紧致差 分格式求解二维稳态不可压涡量-流函数N-S方程组,计算了驱动方腔问题,雷诺数最大 算到2000 ；Chen[3]和Tian[4]分别构造了四阶指数型紧致差分格式求解对流扩散方程,同时 将这种方法推广到求解不可压的N-S方程组；Li等人[5]提出了定常不可压涡量-流函数形 式N-S方程组的四阶紧致格式，并对驱动方腔流问题进行了数值模拟，雷诺数算到7500 ； Choo和Schultz[6]给出了定常不可压涡量-流函数N-S方程组的九点模板的四阶紧致差分 格式，并对驱动方腔和自然对流问题进行了数值模拟,也得到较好的数值结果；Tian和 Ge [7]提出了定常不可压N-S/Boussinesq方程组的四阶紧致差分格式,并对方腔内自然对 流问题进行了数值模拟；葛永斌等人[8]构造了二维非定常不可压涡量-流函数N-S方程组 的一种四阶紧致差分格式,并采用多重网格迭代方法,计算了封闭方腔内自然对流问题. 本文将文献[7]的工作推广到非定常情形,构造了求解二维非定常不可压涡量-流函 数N-S/Boussinesq方程组的一种四阶隐式紧致差分格式.该格式空间为四阶精度,时间为 两阶精度, 并且是无条件稳定的. 1 高精度紧致差分格式 控制流体流动与传热问题的二维非定常不可压 N-S/Boussinesq 方程组，可用涡量- 流函数形式表示为 ψ +ψ =−ω (1) xx yy 资助项目：国家自然科学基金资助项目(No10662006) T +ψ T −ψ T T =+T (2) t y x x y xx yy ω +ψω −ψω Pr(ω =+ω ) +Ra Pr T + f (3) t y x x y xx yy x 其中，ψ 为流函数,ω为涡量,T 为温度,Pr 和Ra