SST湍流模型在强逆压流模拟中应用.docVIP

SST湍流模型在强逆压流模拟中应用 　　摘 要：本文对SST两方程湍流模型提出了一些改进，并将之与Navier-Stokes方程的显式有限元格式相结合，对高雷诺数可压流进行数值模拟。流动方程和湍流方程以一种松散的关联方式求解。时间方向上，流动方程用显式多步Runge-Kutta格式推进，而湍流方程用多步点隐式格式推进以解除源项引起的时间步长限制。数值实践中我们考核了跨音速激波、大攻角等引起的强逆压流动的经典算例，并与实验结果进行了比较，结果是令人满意的。 　　关键词：湍流模型；有限元；N-S方程；强逆压流动 　　中图分类号：O355 文献标志码：A 　　 　　Applications of SST Turbulence Model in Prediction of Flows 　　with Strong Adverse Pressure Gradient 　　TANG Shimin 　　(Changzhou Radio and TV University, Changzhou Jinagsu 213001,China) 　　Abstract：In this paper, we introduce some improvements in the well-known SST turbulence model and combine it with an explicit finite element scheme for Navier-Stokes equations. The flow equations and turbulence equations are solved in loosely coupled manner. The flow equations are advanced using a multi-stage Runge-Kutta time stepping scheme, while the turbulence equations are advanced in a multi-stage point-implicit scheme to alleviate the time-step restriction due to the source term. We test the approach for strong adverse pressure gradient flows induced by transonic shock and large angle of attack. The comparison with experiment data shows that our numerical results are satisfactory. 　　Key words：turbulence model; explicit finite element; Navier-Stokes equations; adverse pressure flows 　　 　　空气动力学中Navier-Stokes方程的一个主要应用就是模拟湍流。对于大规模的分离流、非定常流以及包含多长度尺度的流动，定义代数长度尺度是困难的，B-L、J-K等代数湍流模型的应用变得非常复杂和不确定。显然，数值方法的改进必须跟更一般的湍流模型的发展和应用同步进行。 　　常用的两方程湍流模型是k-ω和模型k-ε。k-ε模型不能准确地推测逆压流动，在粘性底层中也存在数值困难。Wilcox［1］提出的k-ω模型克服了k-ε模型的许多弱点，特别是在近壁面区域中，其性能比k-ε模型有较大的提高；然而，k-ω模型对自由来流的ω值非常敏感。Menter［2］通过混合函数将这两种模型结合起来，既保留了k-ω模型在近壁面区域、k-ε模型在自由剪切层中各自的优良特性，又克服了k-ω模型对来流ω值的敏感性；在此基础上，通过对涡粘系数定义的限制性修正，考虑了逆压边界层中湍流剪切应力的传输效应，进一步设计出SST（Shear Stress Transport）模型。与k-ω和k-ε模型相比，SST模型的性能有了较大提高，尤其对强逆压流中分离开始和分离程度的推测。 　　强逆压边界层流动是湍流模拟的难点之一。本文对SST模型提出了一些改进，并将之与N-S方程的显式有限元解法相结合。湍流方程用多步点隐式格式进行时间推进以解除源项引起的时间步长限制。数值实验中，我们考核了跨音速激波、大攻角等引起的强逆压流动的经典算例。 　　1 SST湍流模型及改进 　　1.1 SST模型 　　Menter首先把k-ε模型转化为k-ω模型的形式。转化形式和原始k-ω模型的差别在于ω方程右端多了一交叉扩散项，且方程的系数不一样。转化过程中还略