湍流模型的选择依据.docx

湍流模型的选择依据.docx

《湍流模型的选择依据.docx》由会员分享，可在线阅读，更多相关《湍流模型的选择依据.docx（14页珍藏版）》请在冰豆网上搜索。

湍流模型的选择依据

解决湍流的模型总计就是那几个方程，Fluent又从工程和数值的角度进行了整理，下面就是这些湍流模型的详细说明。

FLUENT提供了以下湍流模型：

·Spalart-Allmaras模型

·k-e模型

－标准k-e模型

－Renormalization-group（RNG）k-e模型

－带旋流修正k-e模型

·k-ω模型

－标准k-ω模型

－压力修正k-ω模型

雷诺兹压力模型

大漩涡模拟模型

几个湍流模型的比较：

从计算的角度看Spalart-Allmaras模型在FLUENT中是最经济的湍流模型，虽然只有一种方程可以解。

由于要解额外的方程，标准k-e模型比Spalart-Allmaras模型耗费更多的计算机资源。

带旋流修正的k-e模型比标准k-e模型稍微多一点。

由于控制方程中额外的功能和非线性，RNGk-e模型比标准k-e模型多消耗10～15%的CPU时间。

就像k-e模型，k-ω模型也是两个方程的模型，所以计算时间相同。

比较一下k-e模型和k-ω模型，RSM模型因为考虑了雷诺压力而需要更多的CPU时间。

然而高效的程序大大的节约了CPU时间。

RSM模型比k-e模型和k-ω模型要多耗费50～60%的CPU时间，还有15～20%的内存。

除了时间，湍流模型的选择也影响FLUENT的计算。

比如标准k-e模型是专为轻微的扩散设计的，然而RNGk-e模型是为高张力引起的湍流粘度降低而设计的。

这就是RNG模型的缺点。

同样的，RSM模型需要比k-e模型和k-ω模型更多的时间因为它要联合雷诺压力和层流。

概念：

1.雷诺平均：

在雷诺平均中，在瞬态N-S方程中要求的变量已经分解位时均常量和变量。

相似的，像压力和其它的标量

这里

表示一个标量如压力，动能，或粒子浓度。

2.Boussinesq逼近从雷诺压力转化模型：

利用Boussinesq假设把雷诺压力和平均速度梯度联系起来：

Boussinesq假设使用在Spalart-Allmaras模型、k-e模型和k-ω模型中。

这种逼近方法好处是对计算机的要求不高。

在Spalart-Allmaras模型中只有一个额外的方程要解。

k-e模型和k-ω模型中又两个方程要解。

Boussinesq假设的不足之处是假设ut是个等方性标量，这是不严格的。

1.Spalart-Allmaras模型（1equ）：

方程是：

这里Gv是湍流粘度生成的，Yv是被湍流粘度消去，发生在近壁区域。

S~是用户定义的。

注意到湍流动能在Spalart-Allmaras没有被计算，但估计雷诺压力时没有被考虑。

特点：

1）.Spalart-Allmaras模型是设计用于航空领域的，主要是墙壁束缚流动，而且已经显示出和好的效果。

2）。

在原始形式中Spalart-Allmaras模型对于低雷诺数模型是十分有效的，要求边界层中粘性影响的区域被适当的解决。

3）。

不能依靠它去预测均匀衰退，各向同性湍流。

还有要注意的是,单方程的模型经常因为对长度的不敏感而受到批评，例如当流动墙壁束缚变为自由剪切流。

2.标准k-e模型（2equ）：

标准k-e模型的方程

湍流动能方程k，和扩散方程e：

方程中Gk表示由层流速度梯度而产生的湍流动能，计算方法在10.4.4中有介绍。

Gb是由浮力产生的湍流动能，10.4.5中有介绍，YM由于在可压缩湍流中，过渡的扩散产生的波动，10.4.6中有介绍，C1，C2，C3，是常量，σk和σe是k方程和e方程的湍流Prandtl数，Sk和Se是用户定义的。

特点：

标准k-e模型自从被LaunderandSpalding提出之后，就变成工程流场计算中主要的工具了。

适用范围广、经济、合理的精度，这就是为什么它在工业流场和热交换模拟中有如此广泛的应用了。

它是个半经验的公式，是从实验现象中总结出来的。

3.RNGk-e模型（2equ）：

RNGk-e模型的方程

Gk是由层流速度梯度而产生的湍流动能，10.4.4介绍了计算方法，Gb是由浮力而产生的湍流动能，10.4.5介绍了计算方法，YM由于在可压缩湍流中，过渡的扩散产生的波动，10.4.6中有介绍，C1，C2，C3，是常量，ak和ae是k方程和e方程的湍流Prandtl数，Sk和Se是用户定义的。

RNG和标准k-e模型的区别在于：

这里

特点：

RNGk-e模型来源于严格的统计技术。

它和标准k-e模型很相似，但是有以下改进：

·RNG模型在e方程中加了一个条件，有效的改善了精度。

·考虑到了湍流漩涡，提高了在这方面的精度。

·RNG理论为湍流Prandtl数提供了一个解析公式，然而标准k-e模型使用的是用户提供的常数。

·然而标准k-e模型是一种高雷诺数的模型，RNG理论提供了一个考虑低雷诺数流动粘性的解析公式。

这些公式的效用依靠正确的对待近壁区域

这些特点使得RNGk-e模型比标准k-e模型在更广泛的流动中有更高的可信度和精度。

4.带旋流修正的k-e模型（2equ）：

带旋流修正k-e模型的方程

在方程中，Gk是由层流速度梯度而产生的湍流动能，10.4.4介绍了计算方法，Gb是由浮力而产生的湍流动能，10.4.5介绍了计算方法，YM由于在可压缩湍流中，过渡的扩散产生的波动，10.4.6中有介绍，C2，C1e是常量，σk和σe是k方程和e方程的湍流Prandtl数，Sk和Se是用户定义的。

特点：

带旋流修正的k-e模型和RNGk-e模型都显现出比标准k-e模型在强流线弯曲、漩涡和旋转有更好的表现。

由于带旋流修正的k-e模型是新出现的模型，所以现在还没有确凿的证据表明它比RNGk-e模型有更好的表现。

但是最初的研究表明带旋流修正的k-e模型在所有k-e模型中流动分离和复杂二次流有很好的作用。

带旋流修正的k-e模型的一个不足是在主要计算旋转和静态流动区域时不能提供自然的湍流粘度。

这是因为带旋流修正的k-e模型在定义湍流粘度时考虑了平均旋度的影响。

这种额外的旋转影响已经在单一旋转参考系中得到证实，而且表现要好于标准k-e模型。

由于这些修改，把它应用于多重参考系统中需要注意。

5.标准k-ω模型（2equ）：

标准k-ω模型的方程

在方程中，Gk是由层流速度梯度而产生的湍流动能。

Gω是由ω方程产生的。

Tk和Tω表明了k和ω的扩散率。

Yk和Yω由于扩散产生的湍流。

，所有的上面提及的项下面都有介绍。

Sk和Se是用户定义的。

特点：

标准k-ω模型是基于Wilcoxk-ω模型，它是为考虑低雷诺数、可压缩性和剪切流传播而修改的。

Wilcoxk-ω模型预测了自由剪切流传播速率，像尾流、混合流动、平板绕流、圆柱绕流和放射状喷射，因而可以应用于墙壁束缚流动和自由剪切流动。

6.剪切压力传输（SST）k-ω模型（2equ）：

SSTK-

流动方程：

其方程：

和

方程中，

表示湍流的动能，

为

方程，

，

分别代表k与

的有效扩散项

，

分别代表k与

的发散项。

代表正交发散项。

与

用户自定义。

这个公式与标准K-

模型不同，区别在于标准K-

中，

为一常数

而SST模型中，

方程如下：

其中：

特点：

SSTk-ω模型和标准k-ω模型相似，但有以下改进：

·SSTk-ω模型和k-e模型的变形增长于混合功能和双模型加在一起。

混合功能是为近壁区域设计的，这个区域对标准k-ω模型有效，还有自由表面，这对k-e模型的变形有效。

·SSTk-ω模型合并了来源于ω方程中的交叉扩散。

·湍流粘度考虑到了湍流剪应力的传波。

·模型常量不同

这些改进使得SSTk-ω模型比标准k-ω模型在在广泛的流动领域中有更高的精度和可信度。

7.雷诺压力模型（RSM）：

雷诺应力流动方程：

在这些项中，

不需要模型，而

需要建立模型方程使方程组封闭

特点：

由于RSM比单方程和双方程模型更加严格的考虑了流线型弯曲、漩涡、旋转和张力快速变化，它对于复杂流动有更高的精度预测的潜力。

但是这种预测仅仅限于与雷诺压力有关的方程。

压力张力和耗散速率被认为是使RSM模型预测精度降低的主要因素。

RSM模型并不总是因为比简单模型好而花费更多的计算机资源。

但是要考虑雷诺压力的各向异性时，必须用RSM模型。

例如飓风流动、燃烧室高速旋转流、管道中二次流。

8.大涡模拟：

传统的流场计算方法是用N-S方程，即RANS法，在此方法制，所有的湍流流场都可以模拟，其结果可保存。

理论上，LES法处于DNS与RANS之间，大尺寸漩涡用LES法，而小尺寸的漩涡用RANS方程求解，使用LES法的原则如下：

*动量，质量，能量主要由大尺寸漩涡传输

*大涡在流动中期主导作用，它们主要由流动的几何，边界条件来确定。

*小涡不起主导作用（尺寸上），单其解决方法更具有通用性

*当仅有小涡时，更容易建立通用的模型

当解决仅有大涡否则仅有小涡的问题时，所受的限制要比DNS法少的多。

然而在实际工程中，需要很好的网格划分，这需要很大的计算代价，只有计算机硬件性能大幅提高，或者采用并行运算，LES才可能用于实际工程。

如何调用低雷诺数模型

可进行如下操作：

IXy8i8L-LP

  1）、打开Fluent的k-e模型，这是关键。

.iV;[/CHHY!

o

  2）、在Fluent窗口中点击回车键，会出现如下信息：

!

K!

X?

Bh

F&J6n+M`y>

!

m-T+VQr`8G%U7?

adapt/            file/            report/

/`@k*e5uLdS9Wdefine/            grid/            solve/

]9Ju*A|j.oofHdisplay/          parallel/          surface/

#XvO$].o,kaexit              plot/            view/

Q8~3E:

AqtW!

r  3）、复制其中的define并粘贴，然后回车，出现如下信息：

KGt"m!

]U/f;W

P$Y2yT!

G9\/define>`&JR8]Av

boundary-conditions/  materials/          periodic-conditions/

-b~/v.gaKtcustom-field-functions/mixing-planes/      profiles/

0SP}）kVVugrid-interfaces/      models/            units

bFTB3g}zinjections/        operating-conditions/  user-defined/

bu!

d3x,o]$VO/[2T4t  4）、复制其中的models并进行和3）相同的操作，出现如下信息：

/define/models>i7K9`8VN+v

acoustics/          energy?

            steady?

u;~jy-U%C7qq0\crevice-model?

      radiation/          unsteady-1st-order?

cFSkWE0a2iN8vdpm/              solidification-melting?

unsteady-2nd-order?

X,C6@FHeA.O

dynamic-mesh?

        solver/            viscous/

D+T9}k4F6ydS"b  5）复制其中的viscous并进行和3）相同的操作，出现如下信息：

/Ia{'|uW

/define/models/viscous>

E\,F9Mndetached-eddy-simulation?

kw-sst?

            near-wall-treatment/

4Yq_?

~+g-`Minviscid?

          kw-standard?

        reynolds-stress-model?

;S~t&}K

Ske-realizable?

      laminar?

          spalart-allmaras?

%QE（BwZH/Ske-rng?

            large-eddy-simulation?

  turbulence-expert/

B2?

:

eh8`,}2fq-Hke-standard?

        mixing-length?

      user-definedE:

^?

#T）N\9Gy%zWZr

  6）、复制其中的turbulence-expert并进行和3）相同的操作，出现如下信息：

&R.ycNS^9s#KW~

/define/models/viscous/turbulence-expert>5U_}/D&]/y

kato-launder-model?

    low-re-ke?

          turb-non-newtonian?

6\pd#Zwmi（N3]z  7）、此时复制low-re-ke并粘贴，然后回车，出现如下信息：

9G|Bi5s`

/define/models/viscous/turbulence-expert>low-re-ke3BzQvlD（_

Enablethelow-Rek-epsilonturbulencemodel?

[no]

F（e&_"w8g1a/bx  8）此时问你是否要激活低雷诺数k-e模型，输入“Y”然后回车，这样低雷诺数k-e模型就被激活了，在湍流模型中就会出现这一项。

湍流模型选取的准则：

流体是否可压、建立特殊的可行的问题、精度的要求、计算机的能力、时间的限制。

为了选择最好的模型，你需要了解不同条件的适用范围和限制。

       FLUENT软件中提供以下湍流模型：

1Spalart-Allmaras模型；2k-ε模型；3k-ω模型；4雷诺应力模型（RSM）；5大涡模拟模型（LES）。

1.Spalart-Allmaras模型

应用范围：

Spalart-Allmaras模型是设计用于航空领域的，主要是墙壁束缚（wall-bounded）流动，而且已经显示出很好的效果。

在透平机械中的应用也愈加广泛。

在湍流模型中利用Boussinesq逼近，中心问题是怎样计算漩涡粘度。

这个模型被Spalart-Allmaras提出，用来解决因湍流动粘滞率而修改的数量方程。

模型评价：

Spalart-Allmaras模型是相对简单的单方程模型，只需求解湍流粘性的输运方程，不需要求解当地剪切层厚度的长度尺度；由于没有考虑长度尺度的变化，这对一些流动尺度变换比较大的流动问题不太适合；比如平板射流问题，从有壁面影响流动突然变化到自由剪切流，流场尺度变化明显等问题。

Spalart-Allmaras模型中的输运变量在近壁处的梯度要比k-ε中的小，这使得该模型对网格粗糙带来数值误差不太敏感。

Spalart-Allmaras模型不能断定它适用于所有的复杂的工程流体。

例如不能依靠它去预测均匀衰退，各向同性湍流。

2k-ε模型

①标准的k-ε模型：

最简单的完整湍流模型是两个方程的模型，要解两个变量，速度和长度尺度。

在FLUENT中，标准k-ε模型自从被LaunderandSpalding提出之后，就变成工程流场计算中主要的工具了。

适用范围广、经济、合理的精度。

它是个半经验的公式，是从实验现象中总结出来的。

湍动能输运方程是通过精确的方程推导得到，耗散率方程是通过物理推理，数学上模拟相似原型方程得到的。

应用范围：

该模型假设流动为完全湍流，分子粘性的影响可以忽略，此标准κ-ε模型只适合完全湍流的流动过程模拟。

②RNGk-ε模型：

RNGk-ε模型来源于严格的统计技术。

它和标准k-ε模型很相似，但是有以下改进：

a、RNG模型在ε方程中加了一个条件，有效的改善了精度。

b、考虑到了湍流漩涡，提高了在这方面的精度。

c、RNG理论为湍流Prandtl数提供了一个解析公式，然而标准k-ε模型使用的是用户提供的常数。

d、标准k-ε模型是一种高雷诺数的模型，RNG理论提供了一个考虑低雷诺数流动粘性的解析公式。

这些公式的作用取决于正确的对待近壁区域。

这些特点使得RNGk-ε模型比标准k-ε模型在更广泛的流动中有更高的可信度和精度。

③可实现的k-ε模型：

可实现的k-ε模型是近期才出现的，比起标准k-ε模型来有两个主要的不同点：

·可实现的k-ε模型为湍流粘性增加了一个公式。

·为耗散率增加了新的传输方程，这个方程来源于一个为层流速度波动而作的精确方程。

术语“realizable”，意味着模型要确保在雷诺压力中要有数学约束，湍流的连续性。

应用范围：

可实现的k-ε模型直接的好处是对于平板和圆柱射流的发散比率的更精确的预测。

而且它对于旋转流动、强逆压梯度的边界层流动、流动分离和二次流有很好的表现。

可实现的k-ε模型和RNGk-ε模型都显现出比标准k-ε模型在强流线弯曲、漩涡和旋转有更好的表现。

由于带旋流修正的k-ε模型是新出现的模型，所以现在还没有确凿的证据表明它比RNGk-ε模型有更好的表现。

但是最初的研究表明可实现的k-ε模型在所有k-ε模型中流动分离和复杂二次流有很好的作用。

该模型适合的流动类型比较广泛，包括有旋均匀剪切流，自由流（射流和混合层），腔道流动和边界层流动。

对以上流动过程模拟结果都比标准k-ε模型的结果好，特别是可再现k-ε模型对圆口射流和平板射流模拟中，能给出较好的射流扩张。

模型评价：

可实现的k-ε模型的一个不足是在主要计算旋转和静态流动区域时不能提供自然的湍流粘度，这是因为可实现的k-ε模型在定义湍流粘度时考虑了平均旋度的影响。

这种额外的旋转影响已经在单一旋转参考系中得到证实，而且表现要好于标准k-ε模型。

由于这些修改，把它应用于多重参考系统中需要注意。

应用范围：

Wilcoxk-ω模型预测了自由剪切流传播速率，像尾流、混合流动、平板绕流、圆柱绕流和放射状喷射，因而可以应用于墙壁束缚流动和自由剪切流动。

②SSTk-ω模型：

SSTk-ω模型由Menter发展，以便使得在广泛的领域中可以独立于k-ε模型，使得在近壁自由流中k-ω模型有广泛的应用范围和精度。

为了达到此目的，k-ε模型变成了k-ω公式。

SSTk-ω模型和标准的k-ω模型相似，但有以下改进：

·SSTk-ω模型和k-ε模型的变形增长于混合功能和双模型加在一起。

混合功能是为近壁区域设计的，这个区域对标准的k-ω模型有效，还有自由表面，这对k-ε模型的变形有效。

·SSTk-ω模型合并了来源于ω方程中的交叉扩散。

·湍流粘度考虑到了湍流剪应力的传播。

·模型常量不同。

这些改进使得SSTk-ω模型比标准k-ω模型在广泛的流动领域中有更高的精度和可信度。

3k-ω模型

①标准的k-ω模型：

标准的k-ω模型是基于Wilcoxk-ω模型，它是为考虑低雷诺数、可压缩性和剪切流传播而修改的。

标准的k-ε模型的一个变形就是SSTk-ω模型，它在FLUENT中也是可用的

应用范围：

Wilcoxk-ω模型预测了自由剪切流传播速率，像尾流、混合流动、平板绕流、圆柱绕流和放射状喷射，因而可以应用于墙壁束缚流动和自由剪切流动。

②SSTk-ω模型：

SSTk-ω模型由Menter发展，以便使得在广泛的领域中可以独立于k-ε模型，使得在近壁自由流中k-ω模型有广泛的应用范围和精度。

为了达到此目的，k-ε模型变成了k-ω公式。

SSTk-ω模型和标准的k-ω模型相似，但有以下改进：

·SSTk-ω模型和k-ε模型的变形增长于混合功能和双模型加在一起。

混合功能是为近壁区域设计的，这个区域对标准的k-ω模型有效，还有自由表面，这对k-ε模型的变形有效。

·SSTk-ω模型合并了来源于ω方程中的交叉扩散。

·湍流粘度考虑到了湍流剪应力的传播。

·模型常量不同。

这些改进使得SSTk-ω模型比标准k-ω模型在广泛的流动领域中有更高的精度和可信度。

③两个模型的对比

两种模型有相似的形式，有方程k和ω。

SST和标准模型的不同之处是：

·从边界层内部的标准k-ω模型到边界层外部的高雷诺数的k－e模型的逐渐转变。

·考虑到湍流剪应力的影响修改了湍流粘性公式。

一般来说，DES和LES是最为精细的湍流模型，但是它们需要的网格数量大，计算量和内存需求都比较大，计算时间长，目前工程应用较少。

S-A模型适用于翼型计算、壁面边界层流动，不适合射流等自由剪切流问题。

标准K-Epsilon模型有较高的稳定性、经济性和计算精度，应用广泛，适用于高雷诺数湍流，不适合旋流等各相异性等较强的流动。

RNGK-Epsilon模型可以计算低雷诺数湍流，其考虑到旋转效应，对强旋流计算精度有所提供。

RealizableK-Epsilon模型较前两种模型的有点是可以保持雷诺应力与真实湍流一致，可以更加精确的模拟平面和圆形射流的扩散速度，同时在旋流计算、带方向压强梯度的边界层计算和分离流计算等问题中，计算结果更符合真实情况，同时在分离流计算和带二次流的复杂流动计算中也表现出色。

但是此模型在同时存在旋转和静止区的计算中，比如多重参考系、旋转滑移网格计算中，会产生非物理湍流粘性。

因此需要特别注意。

标准K-W模型包含了低雷诺数影响、可压缩性影响和剪切流扩散，适用于尾迹流动、混合层、射流、以及受壁面限制的流动附着边界层湍流和自由剪切流计算。

SSTK-W模型综合了K-W模型在近壁区计算的有点和K-Epsilon模型在远场计算的优点，同时增加了横向耗散导数项，在湍流粘度定义中考虑了湍流剪切应力的输运过程，适用更广，可以用于带逆压梯度的流动计算、翼型计算、跨声速带激波计算等。

雷诺应力模型没有采用涡粘性各向同性假设，在理论上比前面的湍流模型要精确的多，直接求解雷诺应力分量（二维5个，三维7个）输运方程，适用于强旋流动，如龙卷风、旋流燃烧室计算等。