当前位置：首页 > 湍流的数值模拟
水力输送研究[2] 、飞机设计等的瓶颈之一。对湍流基础研究的进展, 可以直接促进许多实际工程及科学应用的进步。
目前, 湍流数值模拟的方法有: 直接数值模拟( Direct Numerical Simulation, DNS) 、大涡数值模拟( Large Eddy Simulation, LES)和雷诺平均模拟( Reynolds Averaged Navier-Stokes, RANS)。
1、直接模拟(DNS)
直接数值模拟（DNS）特点在湍流尺度下的网格尺寸内不引入任何封闭模型的前提下对Navier-Stokes方程直接求解。这种方法能对湍流流动中最小尺度涡进行求解，要对高度复杂的湍流运动进行直接的数值计算，必须采用很小的时间与空间步长，才能分辨出湍流中详细的空间结构及变化剧烈的时间特性。基于这个原因，DNS目前仅限于相对低的雷诺数中湍流流动模型。另外，利用DNS模型对湍流运动进行直接的数值模拟对计算工具有很高的要求，计算机的内存及计算速度要非常的高，目前DNS模型还无法应用于工程数值计算，还不能解决工程实际问题。
2、大涡模拟(LES)
大涡模拟（LES）是基于网格尺度封闭模型及对大尺度涡进行直接求解N-S
方程，其网格尺度比湍流尺度大，可以模拟湍流发展过程的一些细节，但其计算量仍很大，也仅用于比较简单的剪切流运动及管流。大涡模拟的基础是：湍流的脉动与混合主要是由大尺度的涡造成的，大尺度涡是高度的非各向同性，而且随流动的情形而异。大尺度的涡通过相互作用把能量传递给小尺度的涡，而小尺度的涡旋主要起到耗散能量的作用，几乎是各向同性的。这些对涡旋的认识基础就导致了大涡模拟方法的产生。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数的湍流计算模型，适用于离开壁面一定距离的湍流区域。这里的
本文作者：...共分享92篇相关文档
