湍流模型概述

1、大多数飞行器都是在高Re数下飞行，概况的流态是湍流。为了准确地确定湍流流态下的摩阻、热流，湍流成为一个重要而困难的研究课题。之南宫帮珍创作创作时间：贰零贰壹年柒月贰叁拾日（一）DNS目前处理湍流数值计算问题有三种方法，第一种方法即所谓直接数值模拟方法（DNS方法），直接求解湍流运动的N-S方程，得到湍流的瞬时流场，即各种尺度的随机运动，可以获得湍流的全部信息。随着现代计算机的发展和先进的数值方法的研究，DNS方法已经成为解决湍流的一种实际的方法。但由于计算机条件的约束，目前只能限于一些低Re数的简单流动，不克不及用于工程应用。目前国际上正在做的湍流直接数值模拟还只限于较低的需诺数（Re200）

2、和非常简单的流动外形，如平板鸿沟层、完全发展的梢道流，以及后台阶流动等。用直接数值模拟方法处理工程中的复杂流动问题，即使是当前最先进的计算机也还差三个量级。（二）LES另一种方法称做大涡模拟方法（LES方法）。这是一种折衷的方法，即对湍流脉动部分直接地模拟，将N-S方程在一个小空间域内进行平均（或称之为滤波），以使从流场中去掉小尺度涡，导出大涡所满足的方程。小涡对大涡的影响会出现在大涡方程中，再通过建立模型（亚格子尺度模型）来模拟小涡的影响。由创作时间：贰零贰壹年柒月贰叁拾日于湍流的大涡结构强烈地依赖于流场的鸿沟形状和鸿沟条件，难以找出普遍的湍流模型来描述具有分歧的鸿沟特征的大涡结构，宜做直接

3、模拟。相反地，小尺度涡对鸿沟条件不存在直接依赖关系，而且一般具有各向同性性质。所以亚格子模型具有更大的普适性，比较容易构造，这是它比雷诺平均方法要优越的地方。自从1970年Deardorff第一次给出具有工程意义的LES计算以来,LES方法已经成为计算湍流的最强有力的工具之一，应用的方向也在逐步扩展，但是仍然受计算机条件等的限制，使之成为解决大量工程问题的成熟方法仍有很长的路要走。（三）RANS目前能够用于工程计算的方法就是模式理论。所谓湍流模式理论，就是依据湍流的理论知识、实验数据或直接数值模拟结果，对Reynolds应力做出各种假设，即假设各种经验的和半经验的本构关系，从而使湍流的平均Re

4、ynolds方程封闭。随着计算流体力学的发展，湍流模式理论也有了很大的进步，有了非常丰硕的成果。从对模式处理的出发点分歧，可以将湍流模式理论分类成两大类：一类称为二阶矩封闭模式，另一类称涡粘性封闭模式。（1）雷诺应力模式所谓二阶矩封闭模式，是从Reynolds应力满足的方程出发，将方程右端未知的项（生成项，扩散项，耗散项等）用平均流动的物理量和湍流的特征尺度暗示出来。典型的平均流动的变量是平均速度和平均温度的空间导数。这种模式理论，由于保存了Reynolds应力所满足的方程，如果模拟的好，可以较好地反映Reynolds应力随空间和时间的变更规律，因而可以较好地反映湍流运动规律。因此，二阶矩模式

5、是一种较高级的模式，但是，由于保存了Reynolds应力的方程，加上平均运动的方程整个方程组总计15个方程，是一个庞大的方程组，应用这样一个庞大的方程组来解决实际工程问题，计算量很大，这就极大地限制了二阶矩模式在工程问题中的应用。(2)涡粘性模式在工程湍流问题中得到广泛应用的模式是涡粘性模式。这是由Boussinesq仿照分子粘性的思路提出的，即设Reynolds应力为，2,、2,八UiUjT(Ui,jUj,i-Uk,kij)Tkij()33这里kjUb是湍动能，T称为涡粘性系数，这是最早提出的基准涡粘性模式，即假设雷诺应力与平均速度应变率成线性关系，当平均速度应变率确定后，六个雷诺应力只需要

6、通过确定一个涡粘性系数t就可完全确定，且涡粘性系数各向同性，可以通过附加的湍流量来模化，比方湍动能k,耗散率，比耗散率以及其它湍流量k/,lk3/2/,q4,根据引入的湍流量的分歧，可以得到分歧的涡粘性模式，比方罕见的k,k-w模式，以及后来不竭得到发展的k,q-w,k-l等模式，涡粘性系数可以分别暗示为,2TCk2/,Tc,Tck,Tcq-,Tckl.为了使控制方程封闭，引入多少个附加的湍流量，就要同时求解多少个附加的微分方程，根据求解的附加的微分方程的数目，一般可将涡粘性模式划分为三类：零方程模式，半方程模型，一方程模式，两方程模式。1)零方程模式所谓零方程模式是试图直接用平均流动物理量模

7、化一而不引入任何湍流量(如k,等)。例如，Prandttl的混合长理论就是l2一种零方程模式：(5.7)式中i称为混合长在零方程模式的框架下，得到最为广泛应用的是Baldwin-Lomax模式22。该模式是对湍流鸿沟层的内层和外层采取分歧的混合长假设。这是因为靠近壁面处，湍流脉动受到很大的抑制，含能涡的尺度减小很多，因此长度尺度减小很多；另一方面，在鸿沟层外缘，湍流呈间歇状，质量、动量和能量的输运能力大大下降，即湍流的扩散能力减小。这样，应用混合长理论来确定涡粘性系数在这两个分歧的区域应该有分歧的形式。Baldwin-Lomax模式的具体数学描述如下。(5.8)(T)innyyc(T)onty

8、yc这里yc是(T)inn(T)ont的离壁面最小距离V值。对于内层，即yyc,有(T)innI'()是涡量，|ijkUk,j,I是长度尺度Iky(1exy(y/A)()其中是Karman常数，A+是模化常数，y是无量纲法向距离:而u是摩擦速度，其含义为此处下标w暗示壁面。(5.11)对于外层，即yyc,有(T)out(FwakeFkieb(y)其中Fmax是下列函数的最大值：而ymax是F(y)达到最大值的位置。Fkieb是所谓的Klebanoff间歇函数：Udf是平均速度分布中最大值和最小值之差。几个模化常数的值如下：由上述模化关系中可以看出，Reynolds应力完全地由当时当地的

9、平均流参数用代数关系式所决定。平均流场的任何变更立刻为当地的湍流所感知，这标明零方程模式是一个平衡态模式，假定湍流运动永远处于和平均运动的平衡之中。实际上对大多数湍流运动而言，并不是如此，特别是对平均流空间和时间有剧烈变更的情形，再有因为坐标y显式地出现在湍流模式中，零Baldwin-Lomax模方程模式不具有张量不变性，当将它应用到复杂几何外形的流动的数值模拟会带来困难。当流动发生分离时，式会遇到困难，这是因为在分离点和再附点附近，摩擦速度u为零，此时要引入一些人为的干涉来消除这些困难。计算实践标明，只要流动是附体的，零方程模式一般都可以较好地确定压强分布，但是摩阻和传热率的估算不敷准确，特

10、别是当流动有分离和再附时。这是因为附体流压强分布对湍流应力不敏感。总之，对附体流动，如果只关心压强分布，应用零方程模式通常可以给出满意的结果，而且模式应用起来十分简便。但是对于我们计算摩阻的需求，零方程模式是不克不及满足要求。对于有分离、再附等复杂流动，零方程模式是不适用的。2）半方程模式为了能计算具有较强压强梯度，特别是较强逆压梯度的非平衡湍流鸿沟层，Johnson-King于1985年提出了一个非平衡代数模型，该模型仍采取涡粘性假设，把涡粘性的分布与最大剪切应力联系在一起，内层涡粘性与外层涡粘性分布用一个指数函数作光滑拟合，外层涡粘性系数作为一个自由参数，由描述最大剪切应力沿流向变更的常微

11、分方程来确定，此常微分方程是由湍流动能方程导出的，故此模型又称为半方程模型。JK模型虽然仍采取涡粘性假设，却包含有雷诺应力模型的特点。由于求解常微分方程比一方程，二方程模型中求解偏微分方程要简单，省时的多，故用JK模型的工作量只略高于通常平衡状态的零方程代数模型的工作量JK模型后又经不竭修正，发展了JK1990A,JK1990J以及JK1992等改进型3)一方程模式Baldwin-Barth(BB)模型是在二方程模型中，将某一导出的应变量作为基本物理量而得到的，应用此一方程模型可防止求解两方程时会遇到的某些数值困难。BB一方程模型所选择的导出应变量为“湍流雷诺数"Rt。BB模型对计算

12、网格的要求低，壁面的网格可以与采取BL代数模型的相当，而不象两方程k-e模型那样要求壁面网格很细，这样就防止了在k-e模型中流场求解的刚性问题。Spalart-Allmaras(SA)模型与BB模型分歧，不是直接利用k-e模型两方程模型加于简化而得，而是从经验和量纲分析出发，由针对简单流动在逐渐弥补发展而适用于带有层流流动的固壁湍流流动的一方程模型，模型中选用的应变量是与涡粘性T相关的量二除在粘性次层外，与T是相等的。上述两种一方程模型具有相似的特点，它们不象代数模型那样需要分为内层模型，外层模型或壁面模型，尾流模型，同时亦不需要沿法向网格寻找最大值，因此易于用到非结构网格中去；但由于在每个时

13、间步长内，需要对整个流场求解一组偏微分方程，故比BL和JK模型更费机时4)两方程模式2.1 k-两方程模式2.1.1 尺度k-两方程模式k-模式是最为人所知和应用最广泛的两方程涡粘性模式，为积分到壁面的不成压缩/可压缩湍流的两方程涡粘性模式，各种分歧版本的k-模式罕见于各种文献中，选择Jones-Launder模式作为一般性介绍。k-模式最初的发展是为了改善混合长(mixing-length)模式和防止复杂流动中湍流长度尺度(turbulentlengthscale)的代数暗示(algebraicprescription)。它求解两个湍流标量k和的输运方程。k方程暗示湍动能输运方程，方程暗示湍

14、动能的耗散率。该模式对较小压力梯度(relativelysmallpressuregradients)下的自由剪切流(free-shear-layerflows)具有较好的结果。对于壁面流动(wallboundedflows),在零或者小平均压力梯度下，模式结果和实验结果符合得较为一致，但是对大的逆压梯度(adversepressuregradients),其结果就不太正确了。另外，在壁面附近，该模式需要壁面衰减函数(wall-dampingfunctions)和较好的网格分布。a.模式方程雷诺应力的涡粘性模型为这里t为涡粘性(eddyviscosity),Sj为平均速度应变率张量(mean-

15、velocitystrain-ratetensor),为流体密度，k为湍动能，j为克罗内克算子(Kroneckerdelta)。涡粘性定义为湍动能k和湍流耗散率的函数基于量纲分析，涡粘性由流体密度，湍流速度尺度(turbulentvelocityscale)k2和长度尺度(length-scale)k3'2/来标度，衰减函数f由湍流雷诺数Retk2/来模化。湍流输运方程可暗示成以下形式湍流能量输运方程能量耗散输运方程这里右端项分别暗示生成项(productionterm)耗散项(dissipationterm)和壁面项(wallterm)。b.模式常数和参数模式中各常数的定义为近壁衰减

16、函数Ret上22Us2yfexp(3.4(10.02Ret)2)2、f210.3exp(Ret)壁面项k小k2和小y这里Us为平行于壁面的流动速度。c.鸿沟条件积分到壁面的无滑移鸿沟条件为可实现性k模式上述尺度k模式，对于高平均切变率流动会出现非物理的结果(例如当Sk/3.7时，其中S丫无离)。为了包管模式的可实现性，模式函数C不该该是常数，而应当是平均庆变率的函数。实验标明，对鸿沟层流动和均匀切变流，C的值是非常分歧的。为此人们根据可实现性对模式的约束条件，建议采取以下形式的C(Reynolds,1987,Shih,1994)1C1-(5.19)*kAoAsU式中而三是在以角速度k旋转的旋转

17、坐标系中得到的平均旋转速率。As,6cos,1cos1晨6w)3SijSjkSkiW3,S/SijSij.(5.20)上述关系式中唯一未确定的系数是A。为简单起见，可以设其为常数。对鸿沟层流动。可以取Ao=4.0。对其他流动，人的数值可以调节。低Reynolds数k上述几种k模式适用于高Reynolds数情形。但是对近壁区，湍流需诺数很低，对湍流动力学而言粘性效应非常重要，此时湍流Reynolds数的效应必须加以考虑。我们研究摩阻的计算关注的恰恰是近壁区，因此低Reynolds数k模式的研究是十分重要的。现将有关结果整理如下:低Reynolds数下的涡粘性和k模式方程为式中所有模化常数如下：其

18、中此处f和fl,f2称为阻尼函数，是一个经验公式用来反映近壁区低雷诺数效应，系数4和8列表如下:I12345-3-5-7-9-12aiX10X10X10X10X10-3-5-7-9-12aiX10X10X10X10X102.1.4罕见k-两方程模式在文献中有许多种k涡粘性模式。为了便于比较，我们将几种罕见的k模式作一归类。它们的主要区别一是在k和的方程及其鸿沟条件，另一方面是阻尼函数f的取法。模式代号作者ChChien,1982LBLamandBremhorst,1981NTNaganoandTagawa,1990LSLaunderandSharma,1974MKMyongandKasagi,

19、1988YSYangandShih,1991S&LShihandLumley,1993CMOTTZhuandShih,1995所有上述八种模式都可以用一个统一的方程组暗示:tCfk(5)dkdtxitkkXi/Dtijxj(6)d_dtXitXi入,1Ui八C1f1tijCxj2f2一E(7)启美的项，D,E,T列表如下：ModelTDEChK2k2y22y-exp(0.5y)LBK00NTk00LSK2y222u2T2yMKk_00YSki02u2T2yS&Lk_022UT2yCMOTTK02u2T2y阻尼函数f,f1和£2对分歧的模式有分歧的暗示式。Modelff

20、1f2Ch1exp(0.115y)LB0.165Rk220.5、(1ek)(1-)RtNT1exp2y1&126Rt3/4dRt11.22exp()363.052、11exp(Rt)fRt21.3exp()6.511exp()26LS3.4exp2(1R/50)2、11.3exp(Rt)MK(11exp(幻)2R21-exp(-)9361y21exp(4)5YSexp(.004y5e63842ey8ey)Rt11.22exp(右)S&Lexp(.004Rk5e63842eRk8eRk)5Rk2CMOTTexp(.004Rk5e3O42e6Rk8e8Rk)5Rk2R11.22ex

21、p(广)36Rt211.22exp(-)136式中Rk,y和R定义为模式中出现的模化常数分别为ModelCC1C2kCh.09LB.09NT.09LS.09MK.09YS.09S&LCMOTT.09Eq(5.19)对分歧的模式有分歧的处理连界条件的方法:ModelB.C.forkwB.C.forwCh002kLB02y2kNT02yLS002kMK02y一2YSk02y4S&L2U0.25u0.251一4CMOTT0.25u20.2512.2其它双方程模式涡粘性系数的量纲为速度x长度，当用k,来模化时，它们之间的关系为TCk2/。我们注意到，对尺度k模式的方程，在固壁上有奇点问

22、题（壁面上湍动能k0）,这是因为模式不尽合理带来的非物理的奇点。此外在计算中由于k,在壁面附近变更剧烈，必须在物面附近将网格划分得非常小，才干得到合理的结果。为了克服这些困难，人们试图寻找其它的湍流量来代替k,。可能的选择有/k,k/,lk3/2/,qn,相应地，涡粘性系数可暗示成：现在就来介绍几种典型的模式：2.2.1 k-w两方程模式(Wilcox)k-模式是最为人所知和应用最广泛的两方程涡粘性模式，为积分到壁面的不成压缩/可压缩湍流的两方程涡粘性模式，最主要文献来自Wilcox。求解湍动能k和它的/k,k/,lk3/2/,q瓜(specificdissipationrate)的对流输运方

23、程已经证明Wilcoxk-模式在粘性子层比k-具有更好的数值稳定性。由于壁面附近，值较大，模式不象k-模式或者其它两方程模式，它不需要显式的壁面衰减函数。对于比较缓的逆压梯度流动，该模式在对数区域给出的结果和实验数据较为符合。a.模式方程雷诺应力的涡粘性模型为这里t为涡粘性(eddyviscosity),Sj为平均速度应变率张量(mean-velocitystrain-ratetensor),为流体密度，k为湍动能，j为克罗内克算子(Kroneckerdelta)。涡粘性定义为湍动能k和比耗散率的函数k和的输运方程为b.模式常数和参数模式中各常数的定义为c.鸿沟条件对鸿沟层流动，壁面无滑移鸿沟

24、条件为k0和10一6一2(yi)这里yi为离开壁面第一个点的距离，且y；1。对称鸿沟条件采取零梯度条件，各种附加的鸿沟条件将在具体流动中讨论。2.2.2 SST两方程模式(Menter)k-SST剪切应力输运(shear-stress-transport)模式在近壁处采取Wilcoxk-模式，在鸿沟层边沿boundarylayeredges)和自由剪切层free-shearlayers)采取k-模式k-形式，其间通过一个混合函数blendingfunction来过渡，属于积分到壁面的不成压缩/可压缩湍流的两方程涡粘性模式。为了有效结合k-和k-模式，统一写成k-形式a.模式方程涡粘性定义为这里

25、是涡量的绝对值，ai0.31,F2是混合函数。T的形式解决了湍流剪切应力在逆压梯度鸿沟层的输运。k和由相应的模式输运方程得到。湍动能输运方程湍流比耗散率方程上式中最后一项代表交错扩散项(cross-diffusionterm)生成项b,模式常数和参数这里CDkmax22,1020XjXj这里CDk代表k-模式中的交叉扩散(cross-diffusion)。SST模式常数由来暗不，用2分别暗示原始k-模式系数和转化的k-模式系数Fii1这里Innermodel系数:Outermodel系数:2.2.3 k模式方程为DkDtk,iUiUjUi,j,i(5.29)DDt,i2一,i,i,iC1UiUjUi,jC2k(5.30)模化常数为对低雷诺数流动有(5.31)2exp(a1Ra2Ra3R3(5.32)其中模式的框架下,Speziale(1990)33提出了下列的模式:(5.33)DkDt二k,ikuiujUi,j(5.34)DDt,i,