方程，其网格尺度比湍流尺度大，可以模拟湍流发展过程的一些细节，但其计算量仍很大，也仅用于比较简单的剪切流运动及管流。大涡模拟的基础是：湍流的脉动与混合主要是由大尺度的涡造成的，大尺度涡是高度的非各向同性，而且随流动的情形而异。大尺度的涡通过相互作用把能量传递给小尺度的涡，而小尺度的涡旋主要起到耗散能量的作用，几乎是各向同性的。这些对涡旋的认识基础就导致了大涡模拟方法的产生。Les大涡模拟采用非稳态的N-S方程直接模拟大尺度涡，但不计算小尺度涡，小涡对大涡的影响通过近似的模拟来考虑，这种影响称为亚格子Reynolds应力模型。大多数亚格子Reynolds模型都是将湍流脉动所造成的影响用一个湍流粘性系数，既粘涡性来描述。
LES对计算机的容量和CPU的要求虽然仍然很高，但是远远低于DNS方法对计算机的要求，因而近年来的研究与应用日趋广泛。
3、Reynolds平均法
许多流体力学的研究和数值模拟的结果表明，可用于工程上现实可行的湍流模拟方法仍然是基于求解Reynolds时均方程及关联量输运方程的湍流模拟方法，即湍流的统观模拟方法。
统观模拟方法的基本思想是用低阶关联量和平均流性质来模拟未知的高阶关联项，从而封闭平均方程组或关联项方程组。虽然这种方法在湍流理论中是最简单的，但是对工程应用而言仍然是相当复杂的。即便如此，在处理工程上的问题时，统观模拟方法仍然是最有效、最经济而且合理的方法。在统观模型中，使用时间最长，积累经验最丰富的是混合长度模型和 K-E模型。其中混合长度模型是最早期和最简单的湍流模型。该模型是建立在层流粘性和湍流粘性的类比、平均运动与湍流的脉动的概念上的。该模型的优点是简单直观、无须增加微分方程。缺点是在模型中忽略了湍流的对流与扩散，对于复杂湍流流动混合长度难以确定。
到目前为止，工程中应用最广泛的是K-E模型。另外针对K-E模型的不足之处，许多学者通过对K-E模型的修正和发展，开始采用雷诺应力模型（DSM）和代数应力模型（ASM）。近年来，DSM模型已用来预报燃烧室及炉内的强旋及浮力流动。很多情况下能够给出优于K-E模型的结果。但是该模型也有不足之处，首
先它对工程预报来说太复杂，其次经验系数太多难以确定，此外，对压力应变项的模拟还有争议。更主要的是，尽管这一模型考虑了各种应变效应，但是其总精度并不总是高于其它模型，这些缺点导致了DSM模型没有得到广泛的应用。总之，虽然从本质上讲DSM模型和ASM模型比K-E模型对湍流流场的模拟更加合理，但DSM和ASM中仍然采用精度不高的E方程，模型中常数的通用性还没有得到广泛的验证，边界条件不好给定，计算也比较复杂。正因为如此，目前用计算解决湍流问题时仍然采用比较成熟的K-E模型。 Reynolds平均法的分类
在这类方程中，将非稳态N-S方程对时间作平均，即把湍流运动看成二个流动的叠加，一是时间平均流动，二是瞬时脉动流动。于是在所得的时均的N-S方程中包含了脉动量乘积的时均值等未知量，称为Reynolds应力，它包括了六个未知量。显然方程的个数小于未知量的个数。要让方程封闭，必须作出假设。
按照假设的不同，Reynolds平均法又可以分为Reynolds应力模型及涡粘模型(湍流粘性系数法)。
1）涡粘模型(湍流粘性系数法)
涡粘模型(湍流粘性系数法)基于Boussinesq假设，它将湍流脉动所造成的附加应力(Reynolds应力)同层流运动应力那样与时均的应变率关联起来，这一假设并无物理基础，且采用各向同性的湍流动力粘度来计算湍流应力，难于考虑旋转流动和表面曲率变化的影响，但以此为基础的湍流模型目前在工程计算却应用最为广泛。
2）Reynolds应力模型
Reynolds应力模型不引入Boussinesq假设，而是试图直接求解Reynolds应力。对此再引入偏微分方程，在建立二个脉动值乘积的时均值方程的过程中，又会引入三个脉动值乘积的时均值，为了封闭，又必须建立微分方程，在建立三个脉动值乘积的时均值方程的过程中，又会引入四个脉动值乘积的时均值，这在理论上是一个不封闭的困难。该模型在四个脉动值乘积这一层次上，加了一个涡量脉动平方平均值的方程式，从而使Reynolds应力方程封闭。这是一个17方程模型，因而对计算机的要求较高，但它克服了涡粘模型的一些缺点，可能是目前最有发展前途的湍流模型。
出于计算量的考虑，将Reynolds应力及热力密度用代数方程而不是微分方程来求解，从而大大减轻了计算量，这就是代数应力模型(ASM)。
四、湍流的数值模拟方法总结。
1、湍流粘性系数法
在湍流粘性系数法是基于Boussinesq假设的，它建立了Reynolds应力相对于平均速度梯度的关系。
在引入Boussinesq假设后，关键问题在于计算湍动粘度 。所谓的粘涡模型，就是把湍动粘度与湍流时均参数联系起来的关系式。依据确定湍动粘度的微分方程数目的多少，涡粘就分为零方程、一方程和双方程模型。
零方程模型，是指不需要微分方程而是用代数关系式把湍流粘性系数与时均值联系起来的模型。多数是基于混合长度理论的。。
在混合长度理论中，湍流粘性系数仅与时均速度场有关，而与湍流的特性参数无关，一方程模型改进了这一缺点。它引入了湍流脉动动能的平方根，作为湍流脉动速度的代表。
在一方程模型中，湍流长度标尺是由经验公式给出的，其实它也应是一个变量，通过微分方程计算。为此，再引入一个耗散率的概念，表示各向同性的小尺度涡的机械能转化为热能的速率。为计算的封闭性，再引入耗散率 的控制方程。这就是k-e二方程模型。
标准的k-e二方程模型假定湍流粘性系数是各向同性的，当它用于强旋流、弯曲壁面或弯曲流线流动时，会产生一定的失真。为此，不少学者提出了对标准k-e模型的修正方案，在Fluent中有RNG k-e模型和Realizable k-e模型。
标准的k-e模型能很好的模拟一般的湍流，RNG k-e模型用于处理高应变率及流线弯曲程度较大的流动，Realizable k-e模型在含有射流和混合流的自由流动、管道内流动、边界层流动以及带有分离的流动时具有优势。
在Fluent中，二方程模型还使用了k-w模型。
必须指出的是，以上模型均是针对湍流发展非常充分的湍流流动来建立的，是针对高Re数的湍流计算模型，适用于离开壁面一定距离的湍流区域。这里的
Re数是以湍流脉动动能的平方根作为速度(又称湍流Re数)计算的。而在Re数比较低的区域，湍流发展不充分，湍流的脉动影响可能不如分子粘性大，在贴近壁面的底层内，流动可能处于层流状态。这时，必须采用特殊的处理，一般有二种解决方法，一种是采用壁面函数法，另一种是采用低Re数的k-e模型。
2、近壁区的处理策略
壁面函数法的基本思想是：对于湍流核心区的流动使用k-e模型求解，而在壁面区不进行求解，直接使用半经验公式将壁面上的物理量与湍流核心区内的求解变量联系起来。它需要把第一个节点布置在对数律层，对第一个节点的值由公式确定。这样，不需要对壁面内的流动进行求解，可直接得到与壁面相邻控制体积的节点变量。各种改进的壁面函数法越来越准确的模拟壁面的相关特性。
为了使k-e模型能够计算到壁面，出现了各种低Re数的k-e模型，它实际上是对k-e模型的控制方程的各个项作出相应的修改，以体现壁面附近流动的各种真实特征。比如：控制方程中的扩散项包括湍流扩散和分子扩散二部分,k控制方程中壁面的脉动动能的耗散是各向异性等影响。
3、Reynolds应力模型(RSM和ASM)
上述的二方程模型都假定湍流粘性系数是各向同性的，这些模型难于反映旋转流动及流动方向表面曲率变化的影响，有必要对湍流脉动应力 及湍流热密度 直接建立微分方程求解。在此模型中，对二个脉动值乘积的时均值方程直接求解，而对三个脉动值乘积的时均值，采用模拟方式计算，不少研究者认为这是目前最有发展前途的湍流模型，这就是Reynolds应力方程模型(RSM)。
为了减轻RSM的计算工作量，将Reynolds应力和热力密度用代数方程式而不是用微分方程来求解，用代数方程去近似的模拟微分方程，这就是代数应力方程模型(ASM)。
对于湍流的计算，初场的选择要合适，初始的k-e分布假定如不合适，会形成一个不合理的场，容易引起迭代发散；为避免由于低阶格式的假扩散误差加剧湍流计算结果的不准确性，建议采用二阶及二阶以上的格式。
结语
本文阐述并分析了三种湍流数值模拟方法的基本思想及其应用特点, 从中能够看出:
1、直接数值模拟( DNS)
仅仅能计算中低Re 数且几何边界简单的湍流流动, 只能做一些探索性工作, 如检验与改进湍流模型等。目前, 从基础理论出发有可能探索出新的高精度的湍流计算途径, 形成适用于湍流各种复杂流场的方法。
2、大涡数值模拟( LES)
在模型构造方法的原理上优于湍流时均方法, 但其理论还处于研究和发展阶段, 至今主要应用在气象和环境科学领域。由于计算机资源不足及亚格子应力模型的不完善, 其在工程问题中的应用还较少, 还不能作为工程设计的工具。当前, 湍流大涡数值模拟需要解决的课题是亚格子模式的改进和复杂几何边界近壁模型的建立。
3、Reynolds平均法
湍流模型多种多样, 没有一个模式能够对所有湍流运动给出满意的预测结果。模拟计算时应根据不同的湍流流动类型选用不同的模型。所有模型采用的都是湍流统计理论, 因此无法绕过方程组不封闭的困难。随着湍流理论其他研究方法的出现, 总体上其重要性相对呈下降趋势, 但雷诺应力模型( RSM) 仍然是一种颇有应用前途的模型, 为提高模型精度, 今后必须探索其所使用的方程的改进方法。
参考文献
（1）李家春 现代流体力学发展的回顾与展望 力学进展，1995.10
（2）路明等 湍流数值模拟方法及其特点分析 河北建筑科技学院学报，2006.6
（3）肖志祥等 湍流模型在复杂流场数值模拟中的应用 计算物理，2003.7 （4）梁政等 论水力旋流器流场数值模拟中湍流模型的选择 天然气工业，2007.3