湍流理论学习
1．层流和湍流
粘性流体的运动存在着两种完全不同的流动状态：层流状态和湍流状态。雷诺首先于1883年通过做圆管内流动实验观察到层流与湍流现象。当圆管中流动速度较小时，管中的流线之间层次分明，互不掺混，这样的流动称为层流。当流速增大后，流体作复杂、无规律、
?随机的不定常运动，称为湍流。流动状态与雷诺数Re、下临界雷诺数Rec和上临界雷诺数Rec??时，有关。当Re?Rec时，流动为层流；当Rec?Re?Rec流动为不稳定过渡状态；当Re?Rec时，流动为湍流。
湍流是在连续介质范畴内流体的不规则运动，它有别于物质分子的不规则运动。具体来说，在极不规则的湍流中，流动的最小时间尺度和最小空间尺度都远远大于分子热运动的相应尺度。因此湍流运动产生的质量和能量的输运将远远大于分子热运动产生的宏观输运。 2．湍流的平均化、雷诺粘性应力
经典的湍流理论认为，湍流是一种完全不规则的随机运动，湍流场中的物理量在时间和空间上呈随机分布，不同的瞬时有不同的值，关注某个瞬时的值是没有意义的。因此，雷诺首创用统计平均方法来描述湍流的随机运动，即对各瞬时量进行平均得到有意义的平均值。
从N-S方程出发，利用平均化运算的法则推导平均物理量满足的方程组。只考虑不可压缩流体情形，假设体力可以忽略，此时，N-S方程具有下列形式
?u?u?u1?p??u?u?v?w?????u??t?x?y?z??x??v?v?v1?p??v?u?v?w?????v??t?x?y?z??y? (1) ??w?w?w?w1?p??u?v?w?????w??t?x?y?z??z??u?v?w????0??x?y?z?运用(1)式中的连续性方程，运动方程可改写为
??u?u2?uv?uw1?p????????u??t?x?y?z??x???v?uv?v2?vw1?p????????v??y?z??y??t?x (2) ?2??w??uw??vw??w??1?p???w??t?x?y?z??z??u?v?w????0??x?y?z?对方程组(2)中各式两边进行平均化运算，并利用平均化运算法则得到
____?22??u?u?uv?uw?u?u?v??u?w?1?p????????????u??t?x?y?z?x?y?z??x?____??v?uv?v2?vw?u?v??v?2?v?w?1?p????????????v??t?x?y?z?x?y?z??y (3) ?____?22??w??uw??vw??w??u?w???v?w???w???1?p???w??t?x?y?z?x?y?z??z??u?v?w????0??x?y?z?考虑到方程组(3)的第四式，方程组(3)中的头三个方程可改写成另一种形式，把脉动项
移到右边，得到
____??2?????u???????u?v??????u?w?????u?u?u?u?p?????u?v?w??????u??????t?x?y?z?x?x?y?z???____??2?????v?????v????u?v??????v?w????v?v?v?p????u?v?w??????v???????t?x?y?z?y?x?y?z????____???????w?2????w????u?w??????v?w???w?w?w??p??u?v?w????????????w??x?y?z??z?x?y?z???t??u?v?w????0?x?y?z??(4)
将上式和应力形式的运动方程对比
?其中P是应力张量，则有
dV?divP dtP??pI?2?S?P?
其中I是单位张量。S是平均运动变形速度张量，P?是对称的二阶张量，它的表达式
如下
___________??2??u???u?v???u?w??????xz___________???2????v??????v?w?? (5) yz????uv?___?________???zz????u?w???v?w???w?2???????xy???xx?P???????yxyy??zx????zy由(5)式得到在湍流运动中除了平均运动的粘性应力外，还多了一项由于脉动所引起的
应力，称为雷诺粘性应力。平均流动元所满足的方程(4)称为雷诺方程。上述雷诺方程是不封闭的，方程个数只有4个，而未知函数有10个，即V、p及六个雷诺粘性应力分量。为了使方程组封闭，必须在雷诺应力即平均速度之间建立补充关系式。 3．湍流模型
所谓湍流模型理论，就是根据理论和经验，对雷诺平均运动方程和雷诺应力项建立表达
式或方程，然后对雷诺应力方程的某些项提出尽可能合理的模型和假设，使方程组封闭求解的理论。
3.1 零方程模型
所谓零方程模型，就是在雷诺方程和连续方程以外，不需要另外再加任何方程式来使方程组封闭。即雷诺粘性应力能直接用某些物理量和物理常数表达出来，没有引进高阶微分方程，又称代数模型。
(1) Boussinesq涡粘性系数法
一般的方法利用Boussinesq假设把雷诺压力和平均速度梯度联系起来：
??ui?uj?2??u???ui?u?j??t?????k??ti??ij (6) ???x??xi??j?xi?3?_____(2) Cebeci-Smith模型
该模型通常写成涡粘性系数形式，在内层有
?y??U?m?(?y)2?1?exp(?)?() (7)
lD??y?式中lD是衰减长度系数，当壁面不存在质量引射时为常数。在外层有
2?m?0.0168U?*? (8)
式中?为边界层位移厚度，U是外部自由流速度，?为间隙因子。
3.2 一方程模型
(1) The Spalart-Allmaras模型
Spalart-Allmaras模型的变量中v是湍流动粘滞率除了近壁区域，方程是：
*??1(?v)?(?vui)?Gv??t?xi?v???????xj??v??v2???(???v)??Cb2?()??Yv?Sv (9)
?xj??xj?????这里Gv是湍流粘度生成的，Yv是被湍流粘度消去，发生在近壁区域。
对于解决动力漩涡粘性，S-A模型是相对简单的方程，不必要去计算和剪应力层厚度相
关的长度尺度。
(2) 简化的二方程模型
由于?方程难模化和精度低，可以在保留k方程的基础上舍弃?方程，将?表示为
k3/2?? (9)
l式中，l可视为混合长度，必须根据具体流场情况另外给出。 c 二方程模型
(1) 标准k??模型
标准k??模型是个半经验公式，主要是基于湍流动能和扩散率。k方程是个精确方程，?方程是个由经验公式导出的方程。
湍流动能方程k方程为
???(?k)?(?kui)??t?xi?xj扩散方程?方程为
??t?k??(??)??Gk?Gb????YM?Sk (10)
??x?kj?????t?????2?(??)??C1?(Gk?C3?Gb)?C2???S?(11)
??xkk??j???方程中Gk表示由层流速度梯度而产生的湍流动能，Gb是由浮力产生的湍流动能，YM由???(??)?(??ui)??t?xi?xj于在可压缩湍流中，过渡的扩散产生的波动。 (2) 标准k??模型
标准k??模型是一种经验模型，是基于湍流能量方程和扩散速率方程。
??k???Gk?Yk?Sk (12) ??k?x??j?????????(??)?(??ui)????G??Y??S? (13) ?????t?xi?xj??xj?在方程中，Gk是由层流速度梯度而产生的湍流动能。G?是由?方程产生的，?k和??表明了k和?的扩散率。Yk和Y?由于扩散产生的湍流。
???(?k)?(?kui)??t?xi?xj4．雷诺应力模式
雷诺应力模式简称为RSM，由于雷诺平均运动方程所包含的雷诺应力为一未知项，所以最直接的办法是建立雷诺应力vi?v?j的方程。 (1) 雷诺应力方程和湍流动能方程
得到vi?v?j的方程步骤为：把N-S方程化成
__________D(vivj)Dt的形式，将瞬时速度表示成平均速
度与脉动速度之和，即vi?vi?vi?、vj?vj?v?将其代入j，平均运动的形式
D(vivj)Dt形式方程，减去
D(vivj)DtD(vivj)Dt_____，最终得到雷诺应力vi?v?j方程为
_____D(vi?v?j)Dt_____??vi?v?j?t_____?vk?vi?v?j?xk_____?__________?_____?v?jp?_______?vi?v?j??vi?p??????????ik?vi?v?jvk??jk??xk????xk?????????????????????湍流扩散项分子扩散项??__?__?_____??_____???v?v?v????v?v?vpjjj???i??2?i??vi?vk?v?jvk??i???xk?xk??xk?xk????????xj?xi?????????????????????????耗散项产生项压力变形项_____________________________ (14)
在模式理论中要用到湍动能k?v?/2，在式(14)中令i=j得到
___2iDk?k?k??vjDt?t?xj????xj?1?1????????????vp?vvv?2?vsiijkiij??vivksij?2?sijsij?2???___________________________ (15)
这就是湍动能方程。
(2) 经模化后的雷诺应力方程和湍动能方程
?、vi?p?等，因此，需要在雷诺应力方程(14)和湍动能方程(15)中存在未知项，如vi?v?jvk建立这些未知量的表达式，这一过程称为湍流模式建立的过程。最终经模化后的雷诺应力方
程和湍动能方程为：
______________________??2Dvi?v?j?vi?v?j???k?vi?v?j2?Ck???Pij??ij?Dt?xl???xl?xl?3?? (16) ????_____22????C1?vi?v?j??ijk??C2?Pij??ijPk?k?33???_____Dk??k2?k?k??C???k??Pk?? (17) Dt?xl???xl?xl?其中
?Pij??????v_________?j_____?__??v?i??vi?vk?v?jv,P_____?v?k??vi?vl?i ??xkk?xk????xl式中经验常数Ck?0.09~0.11,C1?1.5~2.2,C2?0.4~0.5。
以上方程中出项的湍流耗散率?仍是未知项，建立?方程，进行模化得到
D?2__???_____Dt?x??Ck???????vi??2???x????C?1vi?vl??C?2 ll?xl?k?xlk式中经验常数C??0.07~0.09,C?1?1.41~1.45,C?2?1.9~1.92。
联立封闭微分方程组，这就是雷诺应力封闭模式。
(18) (19)