一种强耦合Spalart-Allmaras湍流模型的RANS方程的高效数值计算方法 An Efficient Numerical Method for Coupling the RANS Equations with Spalart-Allmaras Turbulence Model Equation.pdfVIP

航 空 学 报 2013V01．34No．92007-2018 Sep．25 ActaAeronauticaetAstronauticaSinica ISSN1000-6893ON 11—1929／V hkxb．buaa．edu．cnhkxb@buaa．edu．cn http：d 一种强耦合Spalart—Allmaras湍流模型的 RANS方程的高效数值计算方法 杨小权1’2，杨爱明2一，孙刚2 1．上海飞机设计研究院，上海201210 2．复旦大学力学与工程科学系，上海 200433 是解耦的，也称之为松耦合求解。在松耦合求解过程中，RANS方程和湍流模型方程通常采用不同的数值方法异步求 解。这种求解方式很容易产生因两者计算精度不一致而引起的额外数值耗散。为了消除这种耗散，将RANS方程与 合算法，其中对流项离散采用了Roe格式，时间项的离散采用了隐式LU—SGS(Lower-UpperSymmetricGauss—Seidel)格 式，为了提高计算效率，采用了三层V循环多重网格方法。通过翼型／机翼和振荡翼型／机翼等算例验证了本文发展的 强耦合算法不仅具有较好的收敛性，而且计算精度明显优于松耦合算法，特别对于阻力的预测，强耦合算法更加准确。 关键词：Spalart—Allmaras湍流模型；RANS方程；强耦合算法；多重网格方法；双时间方法 中图分类号：V211．3 文献标识码：A 文章编号：1000—6893(2013)09—2007—12 湍流问题至今仍然是困扰整个流体力学界的 程的求解相对独立；第二类是“强耦合”求解方式， 一个难题，湍流的基本机理至今还没有完全弄清， 即在同一次时间推进中RANS方程和湍流模型 这就决定了各种湍流研究方法必然有各自的局限 方程耦合求解。目前第一类求解方式在工程实际 性。湍流的数值模拟大致可分为三类：直接数值 中得到了广泛的应用；而第二类求解方式还比较 Numerical Simulation，DNS)、大涡 模拟(Direct 少见，但国内外学者已进行了一些有价值的研究。 Simulation，LES)和求解雷诺 模拟(LargeEddy Navi— 平均Navier—Stockes(Reynolds—Averaged方程肛￡湍流模型的RANS方程的稳定性，发现 er—Stokes，RANS)方程。鉴于前两种方法需要显式强耦合算法和显式松耦合算法的收敛性大体 极大的计算资源以及方法本身的问题，目前还无 一致。Liu和Zhang[21发展了基于多重网格的显 法适用于高雷诺数工程问题的模拟，因而求解 式强耦合两方程肛叫湍流模型的RANS方程的 RANS方程自然成为工程实际中的最佳方法。 定常可压缩算法，计算结果发现强耦合求解方法 对于一方程和两方程湍流模型而言，其求解 方式可分为两类：第一类是“松耦合”求解方式，即 is[33对比了耦合两方程k-e和是一cU湍流模型的 在同一次时间推进中RANS方程和湍流模型方RANS方程的3种求解方法(即：隐式强耦合解 11：24 收稿日期：2012-10-30；退修日期：2012-11-22；录用日期：2013-01—04；网络出版时间：2013—05—23 htmI 网络出版地址：WWW．cnki．net／kcms／detail／11．1929．V1124．002 基金项目：国家自然科学基金；复旦