可压缩粘性流动笛卡尔网格虚拟单元方法研究.pdfVIP

第 卷 第 期 空 气 动 力 学 学 报 ， ３２ ６ Ｖｏｌ．３２ Ｎｏ．６ 年 月 ACTAAERODYNAMICASINICA ， ２０１４ １２ Ｄｅｃ．２０１４ 文章编号： （ ） ０２５８-１８２５ ２０１４ ０６-０７４８-０７ 可压缩粘性流动笛卡尔网格虚拟单元方法研究 沈志伟 ，赵 宁 ，胡 偶１ １ ２ （ 南京航空航天大学航空宇航学院，江苏南京 ； 中航工业直升机设计研究所，江西景德镇 ） １． ２１００１６ ２． ３３３０００ 摘 要：以二维高雷诺数可压缩粘性流动问题为背景，提出了一种全新的笛卡尔网格虚拟单元方法。基于壁面函 数基本假设，构造了壁面函数 虚拟单元方法（ ），用于定义湍流壁面边界条件。引入参考点的概念计算虚 - WＦ-ＧＣM 拟单元上的基本变量与湍流变量值，定义了“非贴体”笛卡尔网格下的湍流壁面边界条件，并通过壁面函数模型修 正近壁面单元与界面单元。基于自适应笛卡尔网格体系，采用发展的具有二阶精度的格心格式有限体积求解器， 数值模拟了跨音速 ＲＡＥ２８２２ 翼型绕流问题与超音速圆柱绕流问题，计算结果与实验值吻合良好，显示了 WＦ- ＧＣM对高雷诺数可压缩粘性问题是有效的。 关键词：虚拟单元方法；壁面函数；自适应笛卡尔网格；湍流 中图分类号： 文献标识码： ： ／ Ｖ２１１．３ Ａ doi １０．７６３８ ｋ ｄｌｘｘｂ ２０１４．０１０５ ｑ - 可能形成微小的网格单元，造成方程系统的刚性问 ０ 引 言 题，影响流场的收敛特性并在物面边界处产生流场解 网格生成是采用计算流体力学（ 的非物理振荡。为此 ［］２ 提出了融合单元 Ｃｏｍ ｕｔａｔｉｏｎａｌ Ｕｄａ ｋｕｍａｒ ｐ ｙ ， ）进行数值模拟的先行工作，目前 处理技术， 等人［］３ 提出了网格缝合处理 ＦｌｕｉｄＤ ｎａｍｉｃｓ ＣＦＤ Mｕｎｉｋｒｉｓｈｎａ ｙ 较为成熟的网格形式有结构网格、非结构网格和笛卡 技术对切割单元方法进行了改进。另一方面，浸入类 尔网格。笛卡尔网格是一种区别于前两者的非贴体网 边界处理方法也引起越来越多的关注，其与切割单元 格类型。由于早期多采用结构化的直角网格形式，在 方法相比避免了复杂的几何求交运算以及小网格单 处理复杂外形时所需要的网格数目巨大，并未得到广 元出现所带