湍流模型概述资料

1、湍流模型概述大多数飞行器都是在高Re数下飞行，表面的流态是湍流。为了准确地确定湍流流态下的摩阻、热流，湍流成为一个重要而困难的研究课题。（一）DNS目前处理湍流数值计算问题有三种方法，第一种方法即所谓直接数值模拟方法（DNS方法），直接求解湍流运动的N-S方程，得到湍流的瞬时流场，即各种尺度的随机运动，可以获得湍流的全部信息。随着现代计算机的发展和先进的数值方法的研究，DNS方法已经成为解决湍流的一种实际的方法。但由于计算机条件的约束，目前只能限于一些低Re数的简单流动，不能用于工程应用。目前国际上正在做的湍流直接数值模拟还只限于较低的 需诺数（Re200）和非常简单的流动外形，如平板边界层、

2、完全发展的槽道流，以及后台阶 流动等。用直接数值模拟方法处理工程中的复杂流动问题，即使是当前最先进的计算机也 还差三个量级。（二）LES另一种方法称做大涡模拟方法（LES方法）。这是一种折衷的方法，即对湍流脉动部分直接地模拟，将 N-S方程在一个小空间域内进行平均（或称之为滤波），以使从流场中 去掉小尺度涡，导出大涡所满足的方程。小涡对大涡的影响会出现在大涡方程中，再通过建立模型（亚格子尺度模型）来模拟小涡的影响。由于湍流的大涡结构强烈地依赖于流场的边界形状和边界条件，难以找出普遍的湍流模型来描述具有不同的边界特征的大涡结 构，宜做直接模拟。相反地，小尺度涡对边界条件不存在直接依赖关系，而且一

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

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

5、，由于 保留了 Reynolds应力的方程，加上平均运动的方程整个方程组总计15个方程，是一个庞大的方程组，应用这样一个庞大的方程组来解决实际工程问题，计算量很大，这就极大地 限制了二阶矩模式在工程问题中的应用。（2）涡粘性模式在工程湍流问题中得到广泛应用的模式是涡粘性模式。这是由Bouss in esq仿照分子粘性的思路提出的，即设 Reynolds应力为，2 2uiujT （Ui,j U j,i -Uk,k ij） 3k ij（）1这里k -UiUj是湍动能，T称为涡粘性系数，这是最早提出的基准涡粘性模2式，即假设雷诺应力与平均速度应变率成线性关系，当平均速度应变率确定后，六个雷诺应力只需

6、要通过确定一个涡粘性系数T就可完全确定，且涡粘性系数各向同性，可以通过附加的湍流量来模化，比如湍动能k,耗散率，比耗散率以及其它湍流量k/ ，l k3/2/ ，q , k，根据引入的湍流量的不同，可以得到不同的涡粘性模式，比如常见的k ，k-w模式，以及后来不断得到发展的k ，q-w，k-l等模式，涡粘性系数可以分别表示为C k2/T C ， T C k ，2T C ， T C 押 kl.为了使控制方程封闭，弓I入多少个附加的湍流量，就要同时求解多少个附加 的微分方程，根据求解的附加的微分方程的数目，一般可将涡粘性模式划分为 三类：零方程模式，半方程模型，一方程模式，两方程模式。1）零方程模式

7、所谓零方程模式是试图直接用平均流动物理量模化T，而不引入任何湍流量（如 k,l2等）。例如，Pran dttl的混合长理论就是一种零方程模式:(5.7 )式中I称为混合长。在零方程模式的框架下，得到最为广泛应用的是Baldwin-Lomax模式22。该模式是对湍流边界层的内层和外层采用不同的混合长假设。这是因为靠近壁面处，湍流脉动受到很 大的抑制，含能涡的尺度减小很多，因此长度尺度减小很多；另一方面，在边界层外缘， 湍流呈间歇状，质量、动量和能量的输运能力大大下降，即湍流的扩散能力减小。这样， 应用混合长理论来确定涡粘性系数在这两个不同的区域应该有不同的形式。Baldwi n-Lomax模式的

8、具体数学描述如下。T(T ) innyy c(T ) ontyyc(5.8)这里yc是(T )inn(T)ont的离壁面最小距离y值。对于内层，即y yc，有(T)innl2(5.9)是涡量，ijk U k,j,l是长度尺度lky(1 exy( y / A)(5.10)其中k=0.4是Karman常数，A+是模化常数，y是无量纲法向距离:y u y/ w而u是摩擦速度，其含义为此处下标w表示壁面。对于外层，即y yc，有T )out(Fwake F kleb( y)(5.11)其中wakemin(ymaxFmax ,Cwk y maxU2dif/FmaxFmax是下列函数的最大值:F(y) y

9、 (1 ex)( y /A )而ymax是F(y)达到最大值的位置。 Fkleb是所谓的Klebanoff间歇函数：1Ckleb y、6Fkleb(y)1 5.5()y maxU dif是平均速度分布中最大值和最小值之差。几个模化常数的值如下：A 26.0; C 0.02668;Ckleb0.3; Cwk 1.0,; K 0.4.由上述模化关系中可以看出，Reyno Ids应力完全地由当时当地的平均流参数用代数关系式所决定。平均流场的任何变化立刻为当地的湍流所感知，这表明零方程模式是一 个平衡态模式，假定湍流运动永远处于和平均运动的平衡之中。实际上对大多数湍流运动 而言，并非如此，特别是对平均

10、流空间和时间有剧烈变化的情形，再有因为坐标y显式地出现在湍流模式中，零方程模式不具有张量不变性，当将它应用到复杂几何外形的流动的 数值模拟会带来困难。当流动发生分离时，Baldwin-Lomax模式会遇到困难，这是因为在分离点和再附点附近，摩擦速度u为零，此时要引入一些人为的干涉来消除这些困难。计算实践表明，只要流动是附体的，零方程模式一般都可以较好地确定压强分布，但是摩阻和传热率的估算不够准确，特别是当流动有分离和再附时。这是因为附体流压强分布对湍流应力不敏感。总之，对附体流动，如果只关心压强分布，应用零方程模式通常可以给出满意的结果，而且模式应用起来十分简便。但是对于我们计算摩阻的需求，零

11、方程模式是不能满足要求。对于有分离、再附等复杂流动，零方程模式是不适用的。2）半方程模式为了能计算具有较强压强梯度，特别是较强逆压梯度的非平衡湍流边界层,Joh nson-Ki ng于1985年提出了一个非平衡代数模型，该模型仍采用涡粘 性假设，把涡粘性的分布与最大剪切应力联系在一起，内层涡粘性与外层涡 粘性分布用一个指数函数作光滑拟合，外层涡粘性系数作为一个自由参数， 由描述最大剪切应力沿流向变化的常微分方程来确定，此常微分方程是由湍 流动能方程导出的，故此模型又称为半方程模型。JK模型虽然仍采用涡粘性假设，却包含有雷诺应力模型的特点。由于求解常微分方程比一方程，二 方程模型中求解偏微分方程

12、要简单，省时的多，故用JK模型的工作量只略高于通常平衡状态的零方程代数模型的工作量JK模型后又经不断修正，发展了 JK1990A，JK1990J以及JK1992等改进型3) 一方程模式Baldw in-Barth(BB)模型是在二方程模型中，将某一导出的应变量作为基 本物理量而得到的，应用此一方程模型可避免求解两方程时会遇到的某些数 值困难。BB 一方程模型所选择的导出应变量为“湍流雷诺数” Rt。BB模型对计 算网格的要求低，壁面的网格可以与采用BL代数模型的相当，而不象两方程k-e模型那样要求壁面网格很细，这样就避免了在k-e模型中流场求解的刚性问题。Spalart-Allmaras(SA

13、)模型与BB模型不同，不是直接利用 k-e模型两方程模型加于简化而得，而是从经验和量纲分析出发，由针对简单流动在逐渐补 充发展而适用于带有层流流动的固壁湍流流动的一方程模型，模型中选用的 应变量是与涡粘性t相关的量，除在粘性次层外，与t是相等的。上述两种一方程模型具有相似的特点，它们不象代数模型那样需要分为 内层模型，外层模型或壁面模型，尾流模型，同时亦不需要沿法向网格寻找 最大值，因此易于用到非结构网格中去；但由于在每个时间步长内，需要对 整个流场求解一组偏微分方程，故比 BL和JK模型更费机时4) 两方程模式2.1 k-两方程模式2.1.1 标准k-两方程模式k-模式是最为人所知和应用最广

14、泛的两方程涡粘性模式，为积分到壁面的不可压缩/可压缩湍流的两方程涡粘性模式，各种不同版本的k-模式常见于各种文献中，选择Jones-Launder模式作为一般性介绍。k-模式最初的发展是为了改善混合长(mixing-length)模式和避免复杂流动中湍流长度尺度(turbulent length scale)的代数表示(algebraic prescription)。它求解两个湍流标量k和的输运方程。k方程表示湍动能输运方程，方程表示湍动能的耗散率。该模式对较小压力梯度(relatively small pressure gradie nts)下的自由剪切流(free-shear-layer

15、flows)具有较好的结果。对于壁面流动 (wall bou nded flows) ，在零或者小平均压力梯度下，模式结果和实验结果符合得较为一致，但是对大的逆压梯度(adverse pressure gradie nts)，其结果就不太正确了。另外，在壁面附近，该模式需要壁面衰减函数(wall-damping functions)和较好的网格分布。a.模式方程雷诺应力的涡粘性模型为tiju u j 2 t (SijSin ij 3)这里t为涡粘性(eddyviscosity)，Sij为平均速度应变 率张量(mean-velocitystrai n-rate ten sor)为湍动能，ij为克

16、罗内克算子 (Kronecker delta) 。涡粘性定义为湍动能 k和湍流耗散率的函数k2.基于量纲分析，涡粘性由流体密度，湍流速度尺度(turbule nt velocity scale)2k 和长度尺度(length-scale)k32.来标度，衰减函数f由湍流雷诺数Retk2.来模化。XjXj这里右端项分别表示生成项(product ion term)耗散项(dissipation term)和壁面项湍流输运方程可表示成以下形式湍流能量输运方程tij Sij能量耗散输运方程tij Sij(wall term) 。b.模式常数和参数模式中各常数的定义为c 0.09 c 11.45c 2

17、1.921.01.3Prt0.9近壁衰减函数Ret上和©22Us2y2f exp( 3.4 (10.02 Ret)2f210.3exp( Ret)壁面项这里Us为平行于壁面的流动速度。c.边界条件积分到壁面的无滑移边界条件为k 002.1.2可实现性模式上述标准k模式，对于高平均切变率流动会出现非物理的结果（例如当Sk/ 3.7 时,其中S .2SijSij）。为了保证模式的可实现性，模式函数C不应该是C的值是(5.19)ijij-,/Sij Sijij ijijij而ij是在以角速度2 ijkijkk旋转的旋转坐标系中得到的平均旋转速率。常数，而应当是平均庆变率的函数。实验表明，对

18、边界层流动和均匀切变流, 非常不同的。为此人们根据可实现性对模式的约束条件，建议采用以下形式的(Rey no Ids, 1987, Shih, 1994)CAoAsU式中1 1cos C- 6w)3W Sij S jk SkiS3Sij Sij .上述关系式中唯一未确定的系数是A。为简单起见，可以设其为常数。对边界层流动。可以取 A = 4.0。对其他流动，Ad的数值可以调节。2.1.3 低 Reynolds 数 k上述几种k模式适用于高 Reyn olds数情形。但是对近壁区，湍流需诺数很低，对(k),t ( Uik),iT)kk,jUiUjUi,j(5.23),j(),t ( Ui ),i

19、T)C2 f 2C3TS,jS,j,j式中S2 Sj Sj ,Sij1畀Uj,i)ijT(Ui,jUj,iU k,k ij)湍流动力学而言粘性效应非常重要，此时湍流Reynolds数的效应必须加以考虑。我们研究摩阻的计算关注的恰恰是近壁区，因此低Rey no Ids 数 k模式的研究是十分重要的。现将有关结果整理如下：低Reyn olds数下的涡粘性和k模式方程为TC fk(k .一)（5.22）所有模化常数如下:4 AsU1SijU3k,k ij1exp23a1 R a2 R a3R1exp23a1R a2Ra3R45a 4 R 85 R45a4R4a5R5C1max0.43,5C21.9,

20、C31.0k1.0,1.2Sk*U. SijS-ijij ij，其中k2(k )2Rtk2此处f和fi, f2称为阻尼函数，是一个经验公式用来反映近壁区低雷诺数效应,系数ai和ai列表如下:I12345-3-5-7-9-12ai3.3X 10-6.0X 106.6 X 10-3.6 X 108.4X 10-3-5-7-9-12ai2.53 X 10-5.7X 106.55 X 10-3.6 X 108.3X 102.1.4 常见k-两方程模式在文献中有许多种 k涡粘性模式。为了便于比较，我们将几种常见的 k 模式作一归类。它们的主要区别一是在k和 的方程及其边界条件，另一方面是阻尼函数f的取法

21、。模式代号作者ChChien, 1982LBLam and Bremhorst, 1981NTNagano and Tagawa, 1990LSLaunder and Sharma, 1974MKMyong and Kasagi, 1988YSYang and Shih, 1991S&LShih and Lumley, 1993CMOTTZhu and Shih, 1995所有上述八种模式都可以用一个统一的方程组表示:d kdtxitkkXiUitijXjD(6)ddtxitXiC 1 f1tijUi C 2 f2Xj-E有关的项，D, E, T列表如下：ModelTDECh2 exp

22、( yLBNTk0022U 2kIk2 TLS22yykMK00k 、2 22UYS%0T2yk2 22US&L0T2yk2 22UCMOTT0T2y阻尼函数f,T1W f2对不冋的模式有不冋的表示式。Modelff1f21exo( 0.115y )Ch11 .22 exo(氏)36LB(10.165Rk ) 2 (120.5)Rt.05fexp(Rt2)NTexp2y261.3exp(RtLS3.4MKYSS&Lexp(1 Rt /50)2(11 exp( y0)exX2e y.004y3 8e 8y5e4)1 exp(2e 6Rk.004 Rk3小85e48e 8Rk )4

23、.1RT5y25Rr21exX.3exp(11-exp(9exp(1 .22 exo(1 .22 exo(羔)2)R2)£36戌)36CMOTT51 exp( .004 Rk5e R.6 3 8 42e 6Rk8e 8Rk )1 .22 exp(R* 236式中Rk, y和Rt定义为Rtk2模式中出现的模化常数分别为ModelCC 1C 2kCh.091.351.801.01.3LB.091.441.921.01.3NT.091.451.901.41.3LS.091.441.921.01.3MK.091.401.801.41.3YS.091.441.921.01.3S&L.0

24、91.441.921.01.3CMOTTEq (5.19)1.441.921.01.3对不同的模式有不同的处理连界条件的方法:ModelB.C. for kwB.C. for2k2y2k2y2k2yYSS&L00.25 u20.251u4CMOTT0.25 u20.251u42.2其它双方程模式涡粘性系数的量纲为速度X长度，当用k,来模化时，它们之间的关系为2t C k / 。我们注意到，对标准k模式的方程，在固壁上有奇点问题(壁面上湍动能k 0)，这是因为模式不尽合理带来的非物理的奇点。此外在计算中由于k,在壁面附近变化剧烈，必须在物面附近将网格划分得非常小，才能得到合理的结果。为了

25、克服这些困难，人们试图寻找其它的湍流量来代替k,。可能的选择有/k,k/ ,l k3/2 / , q .、k,相应地，涡粘性系数可表示成t C -,2T C k ,T C ,T C kl.现在就来介绍几种典型的模式:2.2.1 k-w两方程模式 (Wilcox)k-模式是最为人所知和应用最广泛的两方程涡粘性模式，为积分到壁面的不可压缩/可压缩湍流的两方程涡粘性模式，最主要文献来自Wilcox。求解湍动能 k 和它的/ k, k / , l k3/2 / , q . k, (specificdissipation rate)的对流输运方程已经证明 Wilcox k-模式在粘性子层比k-具有更好的

26、数值稳定性。由于壁面附近，值较大，模式不象 k-模式或者其它两方程模式，它不需要显式的壁面衰减函数。对于比 较缓的逆压梯度流动，该模式在对数区域给出的结果和实验数据较为符合。a. 模式方程雷诺应力的涡粘性模型为tij2 t(Sij Snn ij3) 2 k ij 3这里t为涡粘性(eddy viscosity) , Sij为平均速度应变 率张量(mean-velocitystrain-rate tensor), k 为湍动能，j 为克罗内克算子(Kronecker delta)涡粘性定义为湍动能k和比耗散率的函数t kk和的输运方程为XjUjk (tijSj -XjUj (t)-Xjk tij

27、Sij-b. 模式常数和参数模式中各常数的定义为3401000.50.5Prt0.9c. 边界条件对边界层流动，壁面无滑移边界条件为106(77这里yi为离开壁面第一个点的距离，且y1+<1。对称边界条件采用零梯度条件，各种附加的边界条件将在具体流动中讨论。2.2.2 SST 两方程模式(Menter)k-SST剪切应力输运(shear-stress-transport)模式在近壁处采用Wilcox k-boundary layer edges)禾口自由剪切层 free-shear layers) 采用 k- k-ble nding fun ctio n积分到壁面的不可压缩 /可压缩湍流

28、的两方程涡粘性模式。为了有效结合k- k-模式，统一写成k-形式 a.模式方程涡粘性定义为a1kmax a1 ; F2这里是涡量的绝对值，a10.31，F2是混合函数。虑 500F2tanh max 2,2-0.99 y yt的形式解决了湍流剪切应力在逆压梯度边界层的输运。k和由相应的模式输运方程得到。k,/、kUjk (k t )tXjXj湍流比耗散率方程uj(t)tXjXj湍动能输运方程上式中最后一项代表交错扩散项生成项tij Sijk22 kP2(1FJ Xj Xj(cross-diffusi on term)P 2SjSnn 03Sijb.模式常数和参数Vk50042kcF1tanhm

29、in max0.99yy2CDk y2这里CDk2 2maxk,1020Xj Xj这里CDk代表k-模式中的交叉扩散(cross-diffusion)SST模式常数a10.31*0.090.41模式参数,k ,由来表示，用1, 2分别表示原始k-模式系数和转化的k-模式系数F11 1 12这里k)77Inner model系数：k10.8510.510.07511 *12厂0.553Outer model系数：k21.0 20.85620.082822 *2 2 J0.440223 k 模式方程为其中I、(u Up)0.22exp( R2/36) 1Speziale(1990) 33提出了下列的模式:k(5.33)C21.921模式的框架下,DkTkk,iUiUjU i,jDtkT(5.34)DT2TDt,i2,i,i ,i22T k,i,i(C11)UiUjUi,j(C2f21)k1k(5.35)DkDtT k,ik,iUiUjUi,jk(5.29)DTT2ii iC1 UjUjUjjC2(5.30)Dt,ik模化常数为k22C1C2CT1 ,e20.92, k0.4,T 1.5.对低雷诺数流动有T C f k