预处理方法对低速充场的数值模拟.pdfVIP

第十届全国计算流体力学会议论文集 预处理方法对低速流场的数值模拟 王振亚 叶友达 国家计算流体力学实验室，北京100083 中国空气动力研究与发展中心．四川绵阳621000 摘要奉文对较复杂外形利用非结构阿梏和有限体积法对低速流场进行数值模拟．对 定常流动采甩预处理方法加速收敛，采用取时间步长结台预处理方法计算非定常流动． 关键词 有限体积方法．非结构阿格．预处理方法，双时间步长 0引言 随着计算流体力学在航空航天以及工业中越来越广泛的应用，人们对它提出了越来越高 的要求： (1)为了满足航空航天中垒机以及工业中建筑群等复杂井型周围流场数值模拟的需要， 促使非结构阿格”“不断发展。 (2)由于一般可压缩计算软件计算低速流动时，会出现收敛速度下降、精度降低的问题， 利用预处理方法”“”求解低速问题成为根有应用意义和理论价值的方向。 (3)本文先对算例进行数值实验。然后针对一些较复杂外形(如多圆柱)数值模拟。此 类流动在实际应用中很常见，所以很有必要研究一下该流场的特征(如涡脱落等)。一般流 动的Re约为104量级，由于在较低Re数与Re数较高时尾部流场大尺度涡脱落现象有一定 相似性【1”，所以本文针对Re=200数值模拟。考虑到Re--200时流动是层流特征，计算中没 有加湍流模型。 1计算网格与数值方法 1．1网格生成 本文中为了模拟壁面附近的边界层采用层推进法生成网格，在此基础上用阵面推进法生 成全场的计算网格。壁面附近嗣格尺度为圆柱直径的o．005，入口边界与圆柱距离为15倍直 径，出口边界与圆柱距离为30倍直径，圈1给出了双圆柱的计算网格示意图。 1．2数值方法 二维守恒型可压缩N-S方程： 塑+堡+盟；上陧+盟1 (1) 甜 缸 却 Re、缸 础7 下标j，v分别代表无粘和有粘通量。采用C．V型控制体(median-dual)，有限体积方法离散 即得： K警+∑【琢】=∑眨1 C2) --47—— 其中，∑表示在控制体的边界积分。 自∞‘ 对流项采用Roe[1”的近似Riemarm方法离散； 万=={[(／8+fL)一f矧△Q】}d目 (3) L分别为右、左特征向量，A为特征值对角矩阵。“一”表示Roe平均变量，并用线性插值 达到二阶精度：Q=Q。+矿·(VQ+出)，妒为限制器9’910 对低速流动，针对方程离散时对流项的雅可比矩阵的刚性，采用预处理方法p]。 上00 0 #RT —!L P 0 0 (4) I’ 占RT o 寺o占矗， 。p 1 1 ’ ’ 暑一p”PVs月r P印j‘J (5) 占=后-，·M； ， M?=max(Mfimtt，rain(1．0，M㈣)) 其中．k是为了区别粘性主导区域与其它区域而i殴gJ--个自由参数： t=一[1，蒜崭]，口=器去 ㈤ 其中：俐=血。面舞忑，Re。为网格雷诺数，ReA=旦等堑，出为网格尺度。 从而方程变为：