各种湍流模式的比较与分析

来源：百度文库《各种湍流模式的比较与分析》一文，作者，张晓前。该文由kk zyb与2014年6月17日分享。

---

研究湍流的目的最终是为了预测和控制湍流。为了做到这些必须了解湍流的本质。随着电子计算机的发展，利用数值模拟方法研究湍流成了必不可少的方法。

( j5 }; f5 c* Y, N( i0 m" y

如果从流动控制方程出发，进行最精细的数值模拟，称为直接数值模拟 (DNS，Direct Numerical Simulation)。

只需预测湍流的统计量时，从雷诺平均方程出发的模拟，称为雷诺平均数值模拟 (RANS, Reynolds AverageNumerical Simulation)。

介于DNS和RANS之间的是大涡模拟 (LES, Large EddySimulation)。LES的思想是大尺度脉动用数值方法计算，小尺度脉动对大尺度的作用通过模型假设处理。

( j% i, f; G  w! ^

以下我们主要讨论雷诺平均数值模拟 (RANS) 的封闭模式。

湍流模式的基本思想

+ i; M! u, p2 F0 y' B- t

在雷诺方程中建立雷诺应力和平均运动场间的关系式，就可以使雷诺方程封闭，此模式称为低阶矩模式。如果以雷诺方程和雷诺应力输运方程联立作为湍流统计平衡的控制方程，这种模式称为二阶矩模式。此外，依据封闭关系式的数学方程模式，湍流模式又可分为代数模式或者微分模式。

4 k9 ~% J: D% G$ I$ C7 b

湍流的直接数值模拟 (DNS) 不存在封闭性问题，原则上可以求解所有湍流问题。由于计算机资源的限制，迄今为止DNS只能求低雷诺数的湍流简单问题。雷诺统计平均模式 (BANS) 计算速度快且符合设计精度要求，成为目前工程应用的唯一方法。

$ l0 H& g( _8 [$ y0 J7 u& b

RANS的缺点是缺乏普适性，要使这种方法称为一种简单可靠的方法还需要做很深入的工作。大涡数值模拟方法 (LES) 只对小尺度脉动做封闭模型，因此可能有较宽的适应范围。

, O/ ]! g0 B2 z7 W+ b

% j% l+ _8 u" e; X) g/ T

建立湍流统计模式的一般原理

建立足够的雷诺应力方程组（代数、微分或一般泛函）使得平均运动方程可解。

, P' ?6 }0 v7 V* u7 e7 Q8 z' E

雷诺应力是平均速度场和初边值条件的泛函，待封闭项是平均速度场、边界速度脉动和边界湍动能的泛函。而泛函的具体形式很难用解析、实验或者直接数值模拟的方式确定。因而湍流统计矩具有以下性质：

• 泛函表达式不是局部的。
  
  % o  \4 Q5 q" i
• 脉动场的泛函形式除了与平均速度场有关，还隐含流场的几何边界。
  
  
  

# O6 L4 x$ u: ^! I0 s

综上所述，湍流脉动的统计矩中包含丰富的脉动场的性质，但是他们不可能用解析方法表达出来。因而我们不得不做各种假设，各个假设又有自己的适用范围。应用湍流模式既需要理性又需要经验。
5 v% K3 V$ [" @. Q! U  ]

封闭模式方程的约束条件

脉动速度的高阶矩是张量，统计量之间是张量关系式。张量关系的封闭式必须满足客观性原则：
( B: s% \0 c$ Z

• 张量函数的可表性原则；
  
  
• 关于参照系统的不变性原则；
  
  9 B0 x  R/ c" ?% S' p' \
• 真实性原则；
  
  
• 渐近性原则。
  
  
  

+ `/ v6 I+ n' g4 d5 M4 w9 W. }

" h( k: o7 M# s+ c2 f$ }- z+ q

; V! A/ m* y$ p  P! t

常用的一些封闭模式

; H& I6 Y+ m) p0 P+ K5 _

• 涡粘模式：是工程中广泛应用的模式，表达式和分子粘性类似。
  " G" t, ]' s% N7 H
  
  ( N8 ]; ]( q1 x. r: ?
• 2阶矩模式：也称雷诺应力模式，比涡粘模式多6个雷诺应力的偏微分方程。
  
  
  / g) ]* X" |+ e( [) t% D7 A: W

代数涡粘模式（零方程模式）

1 S- Y; u. \4 u/ a2 a  l; }" j, e

雷诺应力表达式<ui′uj′>=2vTSij-1/3<δij><ui′ui′>。即bussinesqe涡团粘度模式。Sij 为平均切变率张量。目前在工程上广泛应用的代数涡粘模式是Baldwin-Lomax模式。

代数涡粘模式的评价：

• 计算量少，只要附加粘性模块就可以利用通常的N-S数值计算程序，是最受工程师们欢迎的方法。没有普适性，适合对特定的流动做各种修正。
  
  0 H% y( ~3 J  F$ T
  2 l$ t3 M# p7 n$ p

0 f6 I, |  t- L% j

结论：二维薄层湍流预测结果满意，三维复杂湍流基本不能满意。

$ U( h+ q: X  a% H6 e: @

& u& u1 U7 Y; d

标准k-ε模式（两方程模型1）

由涡粘模式发展而来，把涡粘系数和湍动能及湍动能耗散率联系在一起。vT=Cμk²/ε，Cμ为无量纲系数。

3 [6 Y+ g9 J9 Q! _8 t8 U0 ?8 B

上式依据是脉动动量输运的物理机制。vT∝ql，l∝k3/2/ε，q∝√k，q 为脉动速度的特征量，l 为脉动长度尺度。

在k-ε 模式中，k 和ε 分别在他们的输运方程中解出。

湍动能消耗的机制十分复杂，对它的方程做逐项模化几乎是不可能的。现在采用的是基于类比方法，基本思想是湍动能耗散的生成、扩散以及消失等项和湍动能方程中的对应项有类似的机制和公式。其封闭的方程如下：

模型方程的常数如下：
% ~5 P' O7 ]5 j

关于标准k-ε模式的评价，主要缺点是：

• 假定雷诺应力和当地的平均切应变成正比，不能反应雷诺应力沿流向的松弛效应；
  
  
• 模式是各向同性的，不能反应雷诺应力的各向异性，尤其是近壁湍流；
  
  " q3 n. T3 Q2 O7 i4 m4 [
• 不能反应平均涡量的影响。
  
  
  ! u  f/ O5 t' N& Y

8 T, s+ ~- K5 f$ b! G

为了克服模式的缺点，发展了近代的非线性模式。

* E# e2 Y# a. h, q, p

非线性k-ε模式

4 X" N) P# Q0 v+ i3 Q

我们把雷诺应力的泛函表达式近似为代数表达式，并把他做Taylor展开，并保留到2阶项，得如下二次式：

此模式不仅仅是代数意义上的二次式，它包括了涡粘系数的各向异性，历史效应以及平均涡量的影响。此模式能够预测到方管湍流中的二次流，而标准模式只能预测到流向均匀的单向平均流动。

非线性k-ε 模式较之标准k-ε 模式有很大改进，但仍具有涡粘模式的固有缺陷，如没有包括雷诺应力的松弛效应等。同时，在平均切变流很大的流场内，k-ε 模式有可能不满足真实性条件。

最后，理论上k-ε 模式是以湍动能生成和耗散相平衡为基础的，在固壁处分子粘性扩散将在湍动能中起重要作用，实际应用中我们使用壁函数而不是数值积分对其处理。

& R- [. Y( t1 l2 X

壁面律的最大缺陷是难以推广到三维和复杂的湍流中，如果网格分辨率足够，我们还可以对近壁k-ε 模式修正，称为低雷诺数模型。

; f1 v5 e0 P1 ~, {) B& o/ ]

. A; }; s" m( t4 X+ P2 f4 O

k-ω模式（两方程模型2）

+ o  c0 X9 O3 c, B0 ?6 R( b

ω 为湍动能比耗散率。此模式主要求解湍动能及它的比耗散率的对流输运方程。

己经证明，k-ω 两方程模式在粘性子层比k-ε 模式有更好的数值稳定性。同时由于ω在壁面处较大，此模式不需要显式的壁面衰减函数。

雷诺应力的涡粘性方程

vt 为涡粘性，ρ 为流体密度

k 和ω 的输运方程为
; C+ k  N" S) t  y

各个参数为

边界条件

y1为离开避免第一个点的距离。

6 G  j# L6 r  q) w' K2 e, ]$ p

单方程涡粘系数输运模式（一方程模式）

标准k-ε 模式是局部平衡模式，不能准确预测平均流动有剧烈变化的湍流，如流线曲率突然变化、分离流动以及有激波的可压缩湍流。既要保持涡粘模式的简单形式，又要考虑雷诺应力的松弛效应，Spalart提出了单方程涡粘系数模式。仍然采用涡粘形式的雷诺应力公式，但是放弃vT=Cμk²/ε，而直接导出涡粘的输运方程。对平均应变率表达式的平方求质点导数并代入雷诺应力输运方程，再加上雷诺应力输运方程的封闭模式，就能得到涡粘系数的封闭模式。常用的有Baldwin-Barth  (BB) 模式和Spalart-Allmaras (SA) 模式。

涡粘模式的评价：

• 涡粘模式在工程中应用比较方便，主要缺点是它的局部性（Spalart模式除外）。如果从雷诺应力输运方程出发，雷诺应力的历史效应就可以模拟，此封闭模式成为二阶矩模式。
  / W# L% C' X6 g/ ]8 N$ k9 R, ?
  
  + a0 L3 @7 u% [
  

% f) l2 o& f6 d  J: y/ }, r
( ?) l% F: {; w

2阶矩模式的封闭式(RSM)

脉动速度2阶矩的目的就是封闭以下雷诺应力输运方程：

待封闭项为Φij，Dij 和 Eij。

2阶矩模式分别封闭扩散、耗散和压强速度梯度相关项。

( P- N) V! r6 F9 l9 u

/ V- s1 r3 |' h

雷诺应力再分配项的模式

b 为无量纲化的雷诺应力偏张量：

模式常数由典型湍流的实验结果确定，如下：
  T1 X, Z' p- U6 z/ _& t3 a

最后一项称为壁面项，只在极近壁处有较大贡献，通常忽略或者合并到其他项。
, j) ^& K* P, G% ?. ~' N8 F! g
7 ?# x  [. Z0 r: o% k

雷诺应力耗散项的模式

在湍流统计方程的封闭模式中，湍流输运的耗散模式是最难建立的。耗散包括大涡耗散，小涡耗散，湍动能耗散的扩散项，湍动能耗散的耗散项和湍动能耗散的分子扩散项。目前常用各向同性的耗散模型。εij=2/3εδij。这是不得已而为之，但是目前还未提出工程中更实用、更好的耗散模式。
  w) x+ t% C, D5 w


$ ?4 H8 {2 ]. {9 J& r# \5 _) G

雷诺应力扩散项的模式

. _; w0 @$ H1 X$ f8 ]  L# v" Y扩散项包括3阶矩和压强一脉动速度关联，一般情况，后者相对前者小得多，所以把其放在3阶项里一起处理。扩散项通常采用梯度形式封闭，简单的模型采用各向同性形式。
% e* g" U/ U% I0 V* T

考虑各向异性还有其他推广形式。τij 是雷诺应力。模式常数Cs=0.11。

1 x2 \4 U% U  h

2阶矩模式的评价：

• 包括了雷诺应力的发展过程，能很好的预测复杂湍流。
  ( X4 I) s& p& ?" [% r/ {
  
  
  

: G& Q% j/ m2 f+ N2 Z) K8 D

• 近壁雷诺应力各向异性很强，各向同性的耗散模式有待改进。
  ( i& U  u! S2 b. H+ \3 {2 y, Z
  $ i7 V- y6 n9 H. w/ x) V
  

• 2阶矩模式的不足之处在近壁湍流，强旋转湍流中尤其突出。
  
  % v. g, ?7 ?- \, e2 v% r- X! r
  4 i+ B. v& v8 g" Y# w' X5 l

- O  n9 T% T4 E+ _

代数形式的2阶矩模式(ASM)

4 V% B* T' @+ {! f* X* u" k( m) t雷诺应力输运局部平衡假定：不计雷诺应力质点导数，只有生成、耗散和再分配3项。Rodi提出代数应力模式 (ASM )，忽略雷诺应力沿平均轨迹的变化和扩散，得到雷诺应力的隐式代数方程。
5 U! x7 b& `0 W& S8 F

代入各项得实用公式：

ASM模式是基于局部平衡假定。但ASM放弃了雷诺应力和平均切变率间各向同性的假定，因此在准平衡的三维定常湍流的预测要好一些，但在二维并没有多大优势。

; P3 Q5 C2 }& W& Y5 ]

关于湍流模式的综合评述

湍流模式是目前预测复杂湍流的唯一工具。即使直接数值模拟复杂湍流成为现实，快速准确地预测湍流统计特性的模型仍受工程师们的欢迎。

现有统计模式的致命缺点：没有一个模式能准确预测所有湍流。

7 |6 x; J7 F3 V; a- a3 Y

零方程模式可以很好地模拟附着边界层湍流，模式对自由剪切湍流、附着边界层湍流和适度分离湍流有较高的计算精度。

; k2 t; z( y' a2 e. r( o

可能的两种改进：分区模式（流动分区计算技术和并行算法）和利用张量对湍流场进行分类。

关联阅读：湍流的产生及最新研究进展
& Q& @( @$ \. q, `6 ]* {! U. a6 k+ P

湍流的完全再层流化！
# k4 e7 R% D0 R+ r# U

万亿网格！恒星湍流！
7 Q6 H6 e5 B) l. P% _湍流与涡流的区别
8 [& L1 E7 Z* P( g+ }4 h* o

---

声明：本微信转载文章出于非商业性的教育和科研目的，并不意味着支持其观点或证实其内容的真实性。版权归原作者所有，如转载稿涉及版权等问题，请立即联系我们，我们会予以更改或删除相关文章，保证您的权利！

/ y/ }( Y% X2 y% N' D# T1 v, X