2022年湍流模型概述 .pdf

《2022年湍流模型概述 .pdf》由会员分享，可在线阅读，更多相关《2022年湍流模型概述 .pdf（17页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

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

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

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

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

5、， 由于保留了 Reynolds 应力的方程， 加上平均运动的方程整个方程组总计15 个方程， 是一个庞大的方程组， 应用这样一个庞大的方程组来解决实际工程问题，计算量很大， 这就极大地限制了二阶矩模式在工程问题中的应用。（2）涡粘性模式名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 17 页 - - - - - - - - - 在工程湍流问题中得到广泛应用的模式是涡粘性模式。这是由Boussinesq 仿照分子粘性的思路提出的，即设Reynolds 应力为，ijijkk

6、ijjiTjikUUUuu32)32(,（）这里jiuuk21是湍动能，T称为涡粘性系数，这是最早提出的基准涡粘性模式，即假设雷诺应力与平均速度应变率成线性关系，当平均速度应变率确定后，六个雷诺应力只需要通过确定一个涡粘性系数T就可完全确定，且涡粘性系数各向同性，可以通过附加的湍流量来模化，比如湍动能k，耗散率，比耗散率以及其它湍流量/k，/2/3kl，kq，根据引入的湍流量的不同，可以得到不同的涡粘性模式，比如常见的k，k-w模式，以及后来不断得到发展的k，q-w，k-l 等模式，涡粘性系数可以分别表示为/2kCT，kCT，kCT，2qCT，. l kCT为了使控制方程封闭，引入多少个附加的

7、湍流量，就要同时求解多少个附加的微分方程，根据求解的附加的微分方程的数目，一般可将涡粘性模式划分为三类：零方程模式， 半方程模型，一方程模式，两方程模式。1） 零方程模式所谓零方程模式是试图直接用平均流动物理量模化T，而不引入任何湍流量（如,k等）。例如， Prandttl的混合长理论就是一种零方程模式：yUlT2（5.7 ）式中l称为混合长。在零方程模式的框架下，得到最为广泛应用的是Baldwin-Lomax 模式22。该模式是对湍流边界层的内层和外层采用不同的混合长假设。这是因为靠近壁面处，湍流脉动受到很大的抑制，含能涡的尺度减小很多，因此长度尺度减小很多；另一方面，在边界层外缘，湍流呈间

8、歇状，质量、动量和能量的输运能力大大下降，即湍流的扩散能力减小。这样，应用混合长理论来确定涡粘性系数在这两个不同的区域应该有不同的形式。Baldwin-Lomax模式的具体数学描述如下。contTcinnTTyyyy)()(（5.8 ）这里cy是ontinn)()(TT的离壁面最小距离y值。对于内层，即cyy，有2T)(linn（5.9 ）名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 17 页 - - - - - - - - - 是涡量，lUjkijk,是长度尺度)/(

9、1 (Ayexykyl（5.10 ）其中 k=0.4 是 Karman常数， A+是模化常数，y是无量纲法向距离：wyUy/而u是摩擦速度，其含义为,wyUu此处下标w表示壁面。对于外层，即cyy，有)()(TyFFklebwakeout（5.11 ）其中)/,min(max2maxmaxmaxFUyCFyFdifwkwakemaxF是下列函数的最大值：)/exp(1()(AyyyF而maxy是)(yF达到最大值的位置。klebF是所谓的Klebanoff间歇函数：16max)(5 .51)(yyCyFklebklebdifU是平均速度分布中最大值和最小值之差。几个模化常数的值如下：.4 .0

10、; ,0.1; 3 .0;02668.0; 0.26KCCCAwkkleb由上述模化关系中可以看出，Reynolds 应力完全地由当时当地的平均流参数用代数关系式所决定。 平均流场的任何变化立刻为当地的湍流所感知，这表明零方程模式是一个平衡态模式， 假定湍流运动永远处于和平均运动的平衡之中。实际上对大多数湍流运动而言，并非如此， 特别是对平均流空间和时间有剧烈变化的情形，再有因为坐标y显式地出现在湍流模式中， 零方程模式不具有张量不变性，当将它应用到复杂几何外形的流动的数值模拟会带来困难。当流动发生分离时，Baldwin-Lomax模式会遇到困难，这是因为在分离点和再附点附近，摩擦速度u为零，

11、此时要引入一些人为的干涉来消除这些困难。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 17 页 - - - - - - - - - 计算实践表明，只要流动是附体的，零方程模式一般都可以较好地确定压强分布，但是摩阻和传热率的估算不够准确，特别是当流动有分离和再附时。这是因为附体流压强分布对湍流应力不敏感。总之，对附体流动，如果只关心压强分布，应用零方程模式通常可以给出满意的结果，而且模式应用起来十分简便。但是对于我们计算摩阻的需求，零方程模式是不能满足要求。对于有分离、再

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

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

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

15、程模式k- 模式是最为人所知和应用最广泛的两方程涡粘性模式，为积分到壁面的不可压缩/可压缩湍流的两方程涡粘性模式，各种不同版本的k- 模式常见于各种文献中，选择Jones-Launder模式作为一般性介绍。k- 模式最初的发展是为了改善混合长(mixing-length)模式和避免复杂流动中湍流长度尺度 (turbulent length scale) 的代数表示 (algebraic prescription)。它求解两个湍流标量 k 和 的输运方程。 k 方程表示湍动能输运方程，方程表示湍动能的耗散率。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - -

16、 - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 17 页 - - - - - - - - - 该模式对较小压力梯度(relatively small pressure gradients)下的自由剪切流(free-shear-layer flows) 具有较好的结果。对于壁面流动(wall bounded flows) ，在零或者小平均压力梯度下，模式结果和实验结果符合得较为一致，但是对大的逆压梯度(adverse pressure gradients)，其结果就不太正确了。另外，在壁面附近，该模式需要壁面衰减函数(wall-damping function

17、s)和较好的网格分布。a. 模式方程雷诺应力的涡粘性模型为32)3(2ijijnnijtjtijkSSuu这 里t为 涡 粘 性 (eddy viscosity)，ijS为 平 均 速 度应变 率 张量 (mean-velocity strain-rate tensor)， 为流体密度，k为湍动能，ij为克罗内克算子(Kronecker delta)。涡粘性定义为湍动能k 和湍流耗散率的函数2kfct基于量纲分析，涡粘性由流体密度，湍流速度尺度 (turbulent velocity scale) 2k和长度尺度 (length-scale) 23k来标度，衰减函数f由湍流雷诺数2Rekt来模

18、化。湍流输运方程可表示成以下形式湍流能量输运方程kijtijjkjjjSxkxkuxtk能量耗散输运方程2221kfcSkcxuxtijtijjjj这里右端项分别表示生成项(production term)耗散项 (dissipation term)和壁面项(wall term)。b. 模式常数和参数模式中各常数的定义为09.0c45.11c92.12c0.1k3.19. 0Prt近壁衰减函数名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 17 页 - - - - - -

19、 - - - )Reexp(3 .01)Re02.01 (4.3exp(222ttff和2Rekt壁面项22ykk和2222yust这里su为平行于壁面的流动速度。c. 边界条件积分到壁面的无滑移边界条件为0k02.1.2 可实现性k模式上述标准k模式，对于高平均切变率流动会出现非物理的结果（例如当7.3/Sk时，其中ijijSSS2）。为了保证模式的可实现性，模式函数C不应该是常数，而应当是平均庆变率的函数。实验表明，对边界层流动和均匀切变流，C的值是非常不同的。为此人们根据可实现性对模式的约束条件，建议采用以下形式的C（Reynolds, 1987, Shih, 1994）kUAACs*0

20、1（5.19 ）式中kijkijijkijkijijijijijijSSU2*而ij是在以角速度k旋转的旋转坐标系中得到的平均旋转速率。.,)6(cos31,cos631ijijkijkijsSSSSSSSWwA（5.20 ）上述关系式中唯一未确定的系数是A。为简单起见，可以设其为常数。对边界层流动。可以取A0 = 4.0 。对其他流动，A0的数值可以调节。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 17 页 - - - - - - - - - 2.1.3 低 Rey

21、nolds 数k上述几种k模式适用于高Reynolds 数情形。但是对近壁区，湍流需诺数很低，对湍流动力学而言粘性效应非常重要，此时湍流Reynolds 数的效应必须加以考虑。我们研究摩阻的计算关注的恰恰是近壁区，因此低Reynolds 数k模式的研究是十分重要的。现将有关结果整理如下：低 Reynolds 数下的涡粘性和k模式方程为)(kkfcT（5.22 ）)23.5()()()(,jijijjkTiitUuukkUkjjTjjTiitSSCkfCSfcU,322211,)()()(式中)32(32)(21,2,ijkkijjiTijjiijjiijijijUUUkuuUUSSSS所有模化

22、常数如下：36exp22.01exp1exp131,2.1,0 .10.1,9 .15,43.0max41225544332211554433221,*321*tijkkijijijijijijksRfRaRaRaRaRafRaRaRaRaRafUSSSkSSUCCCkUAC其中名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 17 页 - - - - - - - - - 22321,)(kRkkRt此处f和21, ff称为阻尼函数，是一个经验公式用来反映近壁区低雷诺数效应

23、，系数iiaa 和列表如下：I 1 2 3 4 5 ia3.310-3-6.010-56.610-7-3.610-98.410-12ia2.5310-3-5.710-56.5510-7-3.610-98.310-122.1.4 常见 k-两方程模式在文献中有许多种k涡粘性模式。为了便于比较，我们将几种常见的k模式作一归类。它们的主要区别一是在和k的方程及其边界条件，另一方面是阻尼函数f的取法。模 式 代 号作者Ch Chien, 1982 LB Lam and Bremhorst, 1981 NT Nagano and Tagawa, 1990 LS Launder and Sharma, 1

24、974 MK Myong and Kasagi, 1988 YS Yang and Shih, 1991 S&L Shih and Lumley, 1993 CMOTT Zhu and Shih, 1995 所有上述八种模式都可以用一个统一的方程组表示：kfCt (5) DxUxkxdtkdjitijikti (6) EfCxUfCxxdtdjitijiti22111 (7) 有关的项，D，E，T列表如下：Model T D E Ch k22yk)5.0exp(22yy名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理

25、 - - - - - - - 第 8 页，共 17 页 - - - - - - - - - LB k0 0 NT k0 0 LS k22yk2222yUTMK k0 0 YS k0 222yUTS&L k0 222yUTCMOTT k0 222yUT阻尼函数21,fff和对不同的模式有不同的表示式。Model f1f2fCh )115.0exp(1y1 )36exp(22.12tRLB )5 .201()1(2165.0tRRek305.1f)exp(12tRNT 4/321 .4126exp1tRy1 22)6exp(1 )5.6(exp(3 .1 yRtLS 2)50/1 (4.3expt

26、R1 )exp(3.12tRMK )70exp(1)(45.31 (yRt1 22)5exp(1 )36exp(921yRt名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 17 页 - - - - - - - - - YS )825004.exp(1432865yeyeyey1 )36exp(22.12tRS&L )825004.exp(1432865kkkkReReReR1 )36exp(22.12tRCMOTT )825004.exp(1432865kkkkReReR

27、eR1 )36exp(22.12tR式中tkRyR和,定义为.,2kRyuyykRtk模式中出现的模化常数分别为Model C1C2CkCh .09 1.35 1.80 1.0 1.3 LB .09 1.44 1.92 1.0 1.3 NT .09 1.45 1.90 1.4 1.3 LS .09 1.44 1.92 1.0 1.3 MK .09 1.40 1.80 1.4 1.3 YS .09 1.44 1.92 1.0 1.3 S&L .09 1.44 1.92 1.0 1.3 CMOTT qE(5.19) 1.44 1.92 1.0 1.3 对不同的模式有不同的处理连界条件的方法：Mod

28、el wkforCB.wforCB.Ch 0 0 LB 0 22yk名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 17 页 - - - - - - - - - NT 0 22ykLS 0 0 MK 0 22ykYS 0 22ykS&L 225.0u4251.0uCMOTT 225.0u4251.0u2.2 其它双方程模式涡粘性系数的量纲为速度长度，当用,k来模化时，它们之间的关系为/2kCT。我们注意到，对标准k模式的方程，在固壁上有奇点问题（壁面上湍动能0k），这是

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

30、ate)的对流输运方程名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 17 页 - - - - - - - - - 已经证明Wilcox k-模式在粘性子层比k- 具有更好的数值稳定性。由于壁面附近，值较大，模式不象k- 模式或者其它两方程模式，它不需要显式的壁面衰减函数。对于比较缓的逆压梯度流动，该模式在对数区域给出的结果和实验数据较为符合。a. 模式方程雷诺应力的涡粘性模型为32)3(2ijijnnijttijkSS这 里t为 涡 粘 性 (eddy viscosi

31、ty)，ijS为 平 均 速 度应变 率 张量 (mean-velocity strain-rate tensor), 为流体密度 , k为湍动能 , ij为克罗内克算子(Kronecker delta)。涡粘性定义为湍动能k 和比耗散率的函数ktk 和的输运方程为k-S)(*ij*tijjtjjxkkuxtk2ij-S)(tijjtjjkxuxtb. 模式常数和参数模式中各常数的定义为954031009*5 .05.0*9. 0Prtc. 边界条件对边界层流动，壁面无滑移边界条件为0k和21)(610y这里 y1为离开壁面第一个点的距离，且y1+1。对称边界条件采用零梯度条件，各种附加的边界

32、条件将在具体流动中讨论。2.2.2 SST 两方程模式 (Menter) k-SST 剪切应力输运(shear-stress-transport)模式在近壁处采用Wilcox k-模式，在边界层边缘boundary layer edges) 和自由剪切层free-shear layers)采用 k-模式k-名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 17 页 - - - - - - - - - 形式 ，其间通过一个混合函数blending function来过渡，属

33、于积分到壁面的不可压缩/可压缩湍流的两方程涡粘性模式。为了有效结合k-和 k- 模式，统一写成k-形式a. 模式方程涡粘性定义为211;maxFakaT这里是涡量的绝对值，31.01a，2F 是混合函数。222500,99.02maxtanhyykFT的形式解决了湍流剪切应力在逆压梯度边界层的输运。k 和由相应的模式输运方程得到。湍动能输运方程kSxkkuxtkijtijjtkjj)(湍流比耗散率方程jjjtjjxxkFPxuxt212)1(2)(上式中最后一项代表交错扩散项(cross-diffusion term)，生成项232ijijnnijSSSPb. 模式常数和参数222214,50

34、0,99.0maxmintanhyCDkyykFk这里20210,2maxjjkxxkCD这里kCD代表 k-模式中的交叉扩散(cross-diffusion)。SST模式常数名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 17 页 - - - - - - - - - 31.01a09. 0*41.0模式参数,k由来表示，用21,分别表示原始k-模式系数和转化的k-模式系数21111F这里,kInner model 系数 : 85.01k5. 01075.01553.0

35、*21*11Outer model 系数 : 0.12k856.020828.02440.0*22*222.2.3 k模式方程为kUuukDtDkjijiiikT, (5.29) 2,1,2CUuukCDtDjijiiiTiiT (5.30) 模化常数为.5 .1,4.0,92.0,122221TTkeCkCC对低雷诺数流动有kkfCT (5.31) )(exp133221RaRaRaf (5.32) 其中1)36/exp(22.0192.1)(22tpRCuukR在k模式的框架下，Speziale(1990)33提出了下列的模式：名师资料总结 - - -精品资料欢迎下载 - - - - -

36、- - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 17 页 - - - - - - - - - kkfCT, (5.33) kUuukDtDkjijiiikT, (5.34) ) 1() 1(2222,1,1,2,2fCUuukCkkDtDjijiiiTiiTiiT (5.35) 模化常数及参数分列如下：Model k12C1C2 CSAA 1.36 1.36 1.36 1.44 )36exp(92183.12tR0.09 2ff9.4exp1y70exp145.31yRt2.2.4 q模式 Coakley(1983)24建议采用如

37、下的q涡粘性模式：kqkqfCT,2 (5.39) 22,qUquuqDtDqjijiiiqT (5.40) 22,1,)(CUUUCCDtDjiijjiiiT (5.41) Model qCC1 C2 fCo 1.0 1.3 .09 045.405.f.92 )0065.exp(1kR名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 17 页 - - - - - - - - - 2.3 双尺度两方程模式选择两个长度尺度，一个是典型大涡的尺度，用下标 p 表示， 一个是小

38、涡尺度，用下标t 表示。大涡的pk和p输运方程为12)()(fcyuykytDkDpTpkTpp (5) 22221)()(fckCpyukCpyytDDppTpppTpp (6) 大涡的输运方程反映了大涡对涡动能产生的贡献，由于它是大尺度的，它和边界条件有关，与平均流的应变率密切相关，直接受到可压缩性的影响。这些在大涡输运方程的模化平均应加以考虑。其中1fc 为湍动能和内能之间的交换，Sarkar 等得到ptTtMyuMfc)(23221，akkMtpt)(221为湍流马赫数，2.0,15.032，2fc 是由于小激波产生的能谱输运的增加。最后有pk和p输运方程为ptTtpkTpMyuMyk

39、ytDkDp)1()()1()(2322 (7) pptTpppTpkMCpCpyukCpyytDDp)()()(223221 (8) 小涡tk和t输运方程为tptkTtykytDkDt)( (9) ttttptTtkCtkCtyytDDt)(221 (10) 小涡的主要贡献是湍动能的耗散。可以注意到，（ 9）式中右端p是大涡的能量耗散，它恰好成为小涡湍动能的来源，反映了湍流能量的级串效应。双尺度模式的涡粘性系数采用了大涡和小涡的平均值，为ptptpTkkkkul)()(2321名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 17 页 - - - - - - - - - 模式参数21)1(CpCp，nnCp12，2111CtCt，1122ptptkkkkCpCt，其中n=1.2， =2.2， =1.05 。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 17 页 - - - - - - - - -