湍流模型介绍.doc

《湍流模型介绍.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、湍流模型介绍因为湍流现象是高度复杂的，所以至今还没有一种方法能够全面、准确地对所有流动问题中的湍流现象进行模拟。在涉及湍流的计算中，都要对湍流模型的模拟能力以及计算所需系统资源进行综合考虑后，再选择合适的湍流模型进行模拟。FLUENT中采用的湍流模拟方法包括Spalart-Allmaras模型、standard（标准）k−ε模型、RNG（重整化群）k−ε模型、Realizable（现实）k−ε模型、v2−f模型、RSM（ReynoldsStressModel，雷诺应力模型）模型和LES（LargeEddySimulation，大涡模拟）方法。7.2.1雷诺平均与大涡模

2、拟的对比因为直接求解NS方程非常困难，所以通常用两种办法对湍流进行模拟，即对NS方程进行雷诺平均和滤波处理。这两种方法都会增加新的未知量，因此需要相应增加控制方程的数量，以便保证未知数的数量与方程数量相同，达到封闭方程组的目的。雷诺平均NS方程是流场平均变量的控制方程，其相关的模拟理论被称为湍流模式理论。湍流模式理论假定湍流中的流场变量由一个时均量和一个脉动量组成，以此观点处理NS方程可以得出雷诺平均NS方程（简称RNS方程）。在引入Boussinesq假设，即认为湍流雷诺应力与应变成正比之后，湍流计算就归结为对雷诺应力与应变之间的比例系数（即湍流粘性系数）的计算。根

3、据计算中使用的变量数目和方程数目的不同，湍流模式理论中所包含的湍流模型又被分为二方程模型、一方程模型和零方程模型（代数模型）等大类。FLUENT中使用的三种k−ε模型、Spalart-Allmaras模型、k−ω模型及雷诺应力模型RSM）等都属于湍流模式理论。大涡模拟（LES）方法是通过滤波处理计算湍流的，其主要思想是大涡结构（又称拟序结构）受流场影响较大，小涡则可以认为是各向同性的，因而可以将大涡计算与小涡计算分开处理，并用统一的模型计算小涡。在这个思想下，大涡模拟通过滤波处理，首先将小于某个尺度的旋涡从流场中过滤掉，只计算大涡，然后通过求解附加方程得到小涡的解。过

4、滤尺度一般就取为网格尺度。显然这种方法比直接求解NS方程的DNS方程效率更高，消耗系统资源更少，但却比湍流模式方法更精确。尤其应该注意的是，湍流模式理论无法准确模拟大涡结构，因此在需要模拟大涡结构时，只能采用LES方法1。尽管大涡模拟理论比湍流模式理论更精确，但是因为大涡模拟需要使用高精度的网格，对计算机资源的要求比较高，所以还不能在工程计算中被广泛使用。在绝大多数情况下，湍流计算还要采用湍流模式理论，大涡模拟则可以在计算资源足够丰富的时候尝试使用。7.2.2Spalart-Allmaras模型Spalart-Allmaras模型是一方程模型里面最成功的一个模型，最早

5、被用于有壁面限制情况的流动计算中，特别在存在逆压梯度的流动区域内，对边界层的计算效果较好，因此经常被用于流动分离区附近的计算，后来在涡轮机械的计算中也得到广泛应用。最早的Spalart-Allmaras模型是用于低雷诺数流计算的，特别是在需要准确计算边界层粘性影响的问题中效果较好。FLUENT对Spalart-Allmaras进行了改进，主要改进是可以在网格精度不高时使用壁面函数。在湍流对流场影响不大，同时网格较粗糙时，可以选用这个模型。Spalart-Allmaras模型是一种新出现的湍流模型，在工程应用问题中还没有出现多少成功的算例。如同其他一方程模型一样，Spa

6、lart-Allmaras模型的稳定性也比较差，在计算中1FLUENT中尚未提供DNS计算选项。FLUENT6.1全攻略6采用Spalart-Allmaras模型时需要注意这个特点。7.2.3标准k−ε模型标准k−ε模型由Launder和Spalding提出，模型本身具有的稳定性、经济性和比较高的计算精度使之成为湍流模型中应用范围最广、也最为人熟知的一个模型。标准k−ε模型通过求解湍流动能（k）方程和湍流耗散率（ε）方程，得到k和ε的解，然后再用k和ε的值计算湍流粘度，最终通过Boussinesq假设得到雷诺应力的解。虽然得到了最广泛的使用，但因为标准k−ε模型假定湍

7、流为各向同性的均匀湍流，所以在旋流（swirlflow）等非均匀湍流问题的计算中存在较大误差，因此后来又发展出很多k−ε模型的改进模型，其中包括RNG（重整化群）k−ε模型和Realizable（现实）k−ε模型等衍生模型。7.2.4RNGk−ε模型RNGk−ε模型在形式上类似于标准k−ε模型，但是在计算功能上强于标准k−ε模型，其改进措施主要有：（1）在ε方程中增加了一个附加项，使得在计算速度梯度较大的流场时精度更高。（2）模型中考虑了旋转效应，因此对强旋转流动计算精度也得到提高。（3）模型中包含了计算湍流Prandtl数的解析公式，而不象标准k−