1、湍流理论和湍流模型湍流理论和湍流模型西北工业大学西北工业大学2012年年3月月许和勇许和勇绕圆柱的理想流动：绕圆柱的理想流动：(a)无升力流动无升力流动(b)有升力流动有升力流动1(c)有升力流动有升力流动2(d)有升力流动有升力流动30Re44Re4040Re1903.5*105Re3*106103Re2*105绕圆柱的真实流动绕圆柱的真实流动(P257)：Re=1.54Re=26Re=140粘性流体运动的两种流态粘性流体运动的两种流态-层流层流和和湍流湍流雷诺实验雷诺实验：1883年圆管内流动实验年圆管内流动实验层流层流：管中水流稳定地沿轴向运动，流线之间层次分明、互不掺混，：管中水流稳定
2、地沿轴向运动，流线之间层次分明、互不掺混，流体质流体质点没有垂直于主流方向的横向运动点没有垂直于主流方向的横向运动；湍流湍流：流体作复杂的、无规则的、随机的非定常运动，也称紊流：流体作复杂的、无规则的、随机的非定常运动，也称紊流;上临界流速：层流变湍流上临界流速：层流变湍流下临界流速：湍流变层流下临界流速：湍流变层流ccVV dVdVRVcccc,dVdVRVcccc,)2320(cReRe)13800(ccReReRe)13800(cReRe流动为层流流动为层流流动为不稳定的过渡状态流动为不稳定的过渡状态流动为湍流流动为湍流工程上，将工程上，将下临界雷诺数下临界雷诺数作为流态的判断依据。作为
3、流态的判断依据。1.1 湍流的不规则性湍流的不规则性湍流速度场是时间、空间坐标、实验次数的不规则函数湍流速度场是时间、空间坐标、实验次数的不规则函数),(txuuii在不规则湍流中，流动的最小时间尺度和最小空间尺度都远远大于分在不规则湍流中，流动的最小时间尺度和最小空间尺度都远远大于分子热运动的相应尺度，因此湍流运动产生质量和能量的输运远大于分子热运动的相应尺度，因此湍流运动产生质量和能量的输运远大于分子热运动产生的宏观输运，所以湍流场中质量和能量的平均扩散远大子热运动产生的宏观输运，所以湍流场中质量和能量的平均扩散远大于层流扩散。于层流扩散。随机变量的概率随机变量的概率随机变量的概率密度随机
4、变量的概率密度TiNNuP/)(uuPup/ )()(1.2 湍流的统计湍流的统计随机变量随机变量：湍流速度变量：湍流速度变量u的实数集合，可表示为的实数集合，可表示为 u()事件集合事件集合：相同边界条件下不同初场演化出的所有流场状态：相同边界条件下不同初场演化出的所有流场状态系综系综：所有可能实现的事件集合：所有可能实现的事件集合举例：在相同边界条件下，举例：在相同边界条件下，N个真实初始条件产生个真实初始条件产生N个实验流场（理论上个实验流场（理论上N可以无穷大）是一个系综，其中某一次实验称为一个事件。可以无穷大）是一个系综，其中某一次实验称为一个事件。概率的定义概率的定义：规定全系综的
5、测度为：规定全系综的测度为1，则随机变量，则随机变量u的概率的概率P(x)定义为一切定义为一切u0.3时，一般要考虑压缩时，一般要考虑压缩性影响。性影响。普朗特普朗特(Prandtl)数数, Pr= cp/k, cp是比定压热容，是比定压热容，k是导热系数。是导热系数。Pr数是数是动量交动量交换和热交换之比换和热交换之比，大多数气体，大多数气体Pr小于但是接近小于但是接近1。韦伯韦伯(Weber)数数, We=v2L/, 是表面张力。是表面张力。We数是数是惯性力与表面张力之比惯性力与表面张力之比，在，在大液面曲率如毛细流动、空化起始等过程中很重要。大液面曲率如毛细流动、空化起始等过程中很重要
6、。湍流运动的基本方程湍流运动的基本方程2.1 Navier-Stokes方程和湍流方程和湍流ijjiijijifxxuxpxuutu210iixu)()0 ,(),()0 ,(xxxxiiiiPpVu00)(),(PptUuiixx层流向湍流过渡现象是层流向湍流过渡现象是N-S方程初边值问题解的性质在变化：方程初边值问题解的性质在变化：(1) 层流是小雷诺数下层流是小雷诺数下N-S方程初边值问题的唯一解；方程初边值问题的唯一解；(2) 过渡流动过渡流动(转捩转捩)是是N-S方程的分岔解；方程的分岔解；(3) 高雷诺数的湍流是高雷诺数的湍流是N-S方程的渐近方程的渐近(t)不规则解。不规则解。i
7、jjiijijifxxuxpxuutu210iixu2.2 雷诺方程和脉动运动方程雷诺方程和脉动运动方程ijjiijijifxxuxpxuutu210iixu0iixujjijjijjijjijijxuuxuuxuuxuuxuujijijiuuuuuujjijijjijxuuxuuxuuijjijjiijijifxuuxxuxpxuutu12雷诺平均方程雷诺平均方程2雷诺平均方程雷诺平均方程1)(12jijijjjiijijjijiuuuuxxxuxpxuuxuutu0iixu脉动运动方程脉动运动方程2.3 雷诺应力和雷诺应力输运方程雷诺应力和雷诺应力输运方程) )(111uuuuuu111u
8、uuuuu湍流运动动量通量的平均值湍流运动动量通量的平均值=平均运动的动量通量平均运动的动量通量+脉动动量通量的平均值脉动动量通量的平均值(不可压流动，密度恒定不可压流动，密度恒定)雷诺应力与粘性应力有着量级上和本质上的区别：雷诺应力与粘性应力有着量级上和本质上的区别：(1) 雷诺应力雷诺应力 粘性应力粘性应力(2) 分子运动平均自由程分子运动平均自由程 湍流脉动最小特征尺度湍流脉动最小特征尺度(3) 产生机制不同产生机制不同(1) 设想有一层厚度为设想有一层厚度为的湍流剪切层，流向脉动速度的湍流剪切层，流向脉动速度u1是平均速度是平均速度U的的10%左右，横向脉动速度左右，横向脉动速度u2较
9、较u1小一个量级，所以典型雷诺应力小一个量级，所以典型雷诺应力221001. 0Uuu平均分子粘性应力的量级可估计为平均分子粘性应力的量级可估计为/URe001. 0001. 0/001. 0/221UUUUuu二者比值二者比值在高雷诺数时，在高雷诺数时，510Re，二者比值达到，二者比值达到102的量级。的量级。(2)(3)计算计算NACA2412翼型绕流的以下参数翼型绕流的以下参数:(a)后缘处的边界层厚度；后缘处的边界层厚度；(b)翼型的表面摩擦阻力系数。翼型的表面摩擦阻力系数。流动条件为：攻角为流动条件为：攻角为0度，翼型弦长为度，翼型弦长为1.5米，基于弦长的雷诺数米，基于弦长的雷诺
10、数Rec=3.1*106。由右图的实验数据曲线中可知，由右图的实验数据曲线中可知，0度攻度攻角时翼型阻力系数为角时翼型阻力系数为0.0068。(4) 边界层厚度及摩擦阻力系数计算边界层厚度及摩擦阻力系数计算mcc00426. 0101 . 35 . 1*0 . 5Re0 . 56000754. 0101 . 3328. 1Re328. 16cfC0015. 0000754. 0*2fCNet小于实验测量值小于实验测量值0.0068，仅为，仅为22%因此，用全层流计算结果不准确。因此，用全层流计算结果不准确。层流计算：层流计算：全湍流计算全湍流计算mcc0279. 0)101 . 3(5 . 1
11、*37. 0Re37. 05/165/100372. 0)101 . 3(074. 0Re074. 05/165/1cfC00744. 000372. 0*2fCNet实验值实验值0.0068为为“摩擦阻力摩擦阻力+压差阻力压差阻力”之和，实际摩擦阻力小于之和，实际摩擦阻力小于0.0068，所以全湍流的计算预测值偏大较多所以全湍流的计算预测值偏大较多考虑转捩的计算（假设转捩临界雷诺数考虑转捩的计算（假设转捩临界雷诺数Rex,cr=5*105）:51,105RexVcrx6101 . 3RecVc1613. 0101 . 3105ReRe65,1ccrxcxturbulentcfturbulen
12、tcfturbulentcfCcqCSqD)()()(,turbulentxfturbulentxfCxqD)()(11,1,turbulentxfturbulentcfturbulentxfturbulentcfturbulentxfCxqCcqDDD)()()()()(112,1,arlaxfarlaxfCxqDmin,1min,)()(11turbulentxfturbulentcfarlaxfturbulentxfarlaxffCxqCcqCxqDDD)()()()()(1121,1,min,1,min,turbulentxfturbulentcfarlaxffCcxCCcxC)()(
13、)(11,1,min,1turbulentxfturbulentcfarlaxffCcxCCcxC)()()(11,1,min,100188. 0105328. 1Re328. 1)(5min,11xarlaxfC00372. 0)(,turbulentcfC00536. 0)105(074. 0Re074. 0)(5/155/1,11xturbulentxfC003158. 000536. 0*1613. 000372. 000188. 0*1613. 0fC0063. 0003158. 0*2fCNet考虑转捩的计算结果更加接近实验测量值考虑转捩的计算结果更加接近实验测量值0.0068。t
14、urbulentxfarlaxffCcxCcxC)()(21,2min,100188. 0105328. 1Re328. 1)(5min,11xarlaxfC003856888. 0)1026(074. 0Re074. 0)(5/155/1,22xturbulentxfC003538. 0003856888. 0*8387. 000188. 0*1613. 0fC007076. 0003538. 0*2fCNet比实验测量值偏大。比实验测量值偏大。注意：如果按照公式直接计算注意：如果按照公式直接计算 x2 段的湍流摩擦系数，有较大偏差。段的湍流摩擦系数，有较大偏差。turbulentxftur
15、bulentcfarlaxffCcxCCcxC)(2)()(211,1,min,1001328. 0101328. 1Re328. 1)(6min,11xarlaxfC00372. 0)(,turbulentcfC004669. 0)101 (074. 0Re074. 0)(5/165/1,11xturbulentxfC00264. 0004669. 0*3226. 000372. 0001328. 0*3226. 0fC00528. 000264. 0*2fCNet转捩雷诺数扩大一倍后，摩擦阻力系数减小转捩雷诺数扩大一倍后，摩擦阻力系数减小16.2%如果转捩临界雷诺数增大为如果转捩临界雷诺数
16、增大为1*106，则有，则有%2 .16%1000063. 000528. 00063. 0对于流线型翼型，摩擦阻力与压差阻力的比值（对于流线型翼型，摩擦阻力与压差阻力的比值（Cd=Cf+Cp）估算：）估算：文献文献CFD计算结果：计算结果：NACA0012翼型，翼型，Re=3*106，加入转捩模型，加入转捩模型， Cd=0.00623，Cf=0.00534 Cf/Cd=85.7% Cp/Cd=14.3%上例近似计算结果：上例近似计算结果：NACA2412翼型，翼型，Re=3.1*106，加入转捩模型，加入转捩模型， Cd=0.0068，Cf=0.0063 Cf/Cd=92.6% Cp/Cd=
17、7.4%Lombardi, G., Salvetti, M.V. and Pinelli, D. Numerical Evaluation of Airfoil Friction Drag, Journal of Aircraft, 2000, 37(2): 354-3562.3.2 雷诺应力输运方程雷诺应力输运方程2.3.3 湍动能输运过程湍动能输运过程2.3.4 雷诺应力输运过程雷诺应力输运过程111111222111)(2EDdxxdUuutuu121212222221)(EDdxxdUuutuu222222220EDtuu假定二维平均流的速度分布为假定二维平均流的速度分布为12)(ii
18、xUu，其雷诺输运方程为，其雷诺输运方程为2.3.5 脉动压强分布和压强变形率相关的解析表达式脉动压强分布和压强变形率相关的解析表达式1) ()(41)(2ijvmlmlijjiijjiijrduuxuxuxuxup2)(241ijvlmmlijjirduuxuxu2.3.6 湍流统计方程的封闭性讨论湍流统计方程的封闭性讨论(1) 雷诺方程中出现了雷诺应力项雷诺方程中出现了雷诺应力项(2阶速度相关阶速度相关)(2) 在雷诺应力输运方程中又出现了更高阶的统计相关量在雷诺应力输运方程中又出现了更高阶的统计相关量(3) 如果进一步通过如果进一步通过N-S方程导出高阶相关量的演化方程，则将出现更高方程
19、导出高阶相关量的演化方程，则将出现更高阶的相关量阶的相关量(4阶以上阶以上)结论：从结论：从N-S方程导出的湍流统计方程是永远不封闭的。方程导出的湍流统计方程是永远不封闭的。湍流统计理论的主要任务：研究统计方程的封闭方法。湍流统计理论的主要任务：研究统计方程的封闭方法。2.4 不可压缩湍流的标量输运方程不可压缩湍流的标量输运方程-传热和传质传热和传质qxxxutjjjj20iixujjjjjjxuxxxut2运动方程与温度或浓度无关，即标量场是由速度场确定的，而没有标量场运动方程与温度或浓度无关，即标量场是由速度场确定的，而没有标量场对速度场的反馈作用。这种标量疏运过程，称为对速度场的反馈作用
20、。这种标量疏运过程，称为被动标量输运被动标量输运。2.6 涡量的输运与湍流涡量的输运与湍流iiiijijijixxsxut2流体质点变形率张量流体质点变形率张量jijijjiiijjijiijiixxxxsxut2/2/2/2拟涡能方程拟涡能方程定义定义2/ii称为称为拟涡能拟涡能jjijijijiijijijixxxussxut2iiijiijijixxxuxut2fuu2jxt)(21ijjiijxuxus)(21ijjiijxuxurjijijjijjijjijjjijixussssxutjjijijijxxxuu) (222ijjiijijijjijiijiisssxutjijijji
21、ijiijjiijxxxxxuxu222脉动涡量方程脉动涡量方程脉动涡量的拟涡能方程脉动涡量的拟涡能方程(1) 湍流生成项的主要来源是涡管的伸长湍流生成项的主要来源是涡管的伸长(2) 小尺度湍流是由湍涡拉伸产生的小尺度湍流是由湍涡拉伸产生的(3) 脉动涡量的拉伸时维持湍流的主要机制脉动涡量的拉伸时维持湍流的主要机制湍流的一般定义和描述湍流的一般定义和描述1）湍流场中的流体仍可视为连续介质；）湍流场中的流体仍可视为连续介质；2）物理量呈连续变化，即在空间上和时间上是可微的；）物理量呈连续变化，即在空间上和时间上是可微的；3）可采用描述一般流体运动的方法来建立湍流场数学模型；）可采用描述一般流体运
22、动的方法来建立湍流场数学模型；4）湍流场满足）湍流场满足N-S方程方程5）19世纪初以来，世纪初以来，“湍流是一种完全湍流是一种完全不规则的随机运动不规则的随机运动” 雷诺雷诺首创用首创用统计平均方法统计平均方法来描述湍流运动来描述湍流运动6）20世纪世纪70年代开始，年代开始，“湍流并不是完全随机运动，存在一种可检测和显湍流并不是完全随机运动，存在一种可检测和显示的拟序结构，亦称示的拟序结构，亦称大涡拟序结构大涡拟序结构”，仍存争议，仍存争议7）大多数人观点：由各种大小和涡量不同的涡旋叠加而成，某些情况下，）大多数人观点：由各种大小和涡量不同的涡旋叠加而成，某些情况下，流场作完全随机运动；另
23、一些情况下，流场随机运动和拟序结构并存。流场作完全随机运动；另一些情况下，流场随机运动和拟序结构并存。湍流的统计平均湍流的统计平均瞬时值记为瞬时值记为 A,平均值记为平均值记为A),(1),(dttzyxATtzyxATtt(1) (1) 时间平均时间平均T为时间平均的周期，即要求比湍流脉动周期大得多，以保证得到稳定的为时间平均的周期，即要求比湍流脉动周期大得多，以保证得到稳定的平均值，又要求比流体做非定常运动时的特征时间小得多，以免取平均后，平均值，又要求比流体做非定常运动时的特征时间小得多，以免取平均后，抹平整体的非定常性。抹平整体的非定常性。),(1),(dzdydxtzyxAtzyxA
24、(2) (2) 空间平均空间平均(3) (3) 条件采样平均条件采样平均 规定一个条件准则，对符合该准则的数据进行平均规定一个条件准则，对符合该准则的数据进行平均例如规定一个检测函数例如规定一个检测函数01)(tD湍流信号湍流信号层流信号层流信号 )()()(lim11NiiNiiiNttDtAtDA )(1 )()(1 lim11NiiNiiiNttDtAtDA.,)(, 0, 0,dsAAdsBABABABABABABBAAABtAtAzAzAyAyAxAxABABABABAAAAAAccAAAA则流场处于湍流时的平均为则流场处于湍流时的平均为则流场处于层流时的平均为则流场处于层流时的平均
25、为对于对于瞬时量瞬时量、平均量平均量、脉动量脉动量的运算法则如下：的运算法则如下：;BBBAAAiiivvv2 iv对于对于湍流场速度而言，湍流场速度而言，而，而表示表示湍流强度湍流强度。不可压缩湍流平均运动的基本方程不可压缩湍流平均运动的基本方程000)(iixvtVV00)(0iiiiiiiiiiiixvxvxvxvxvxvvxv221jiijijixvxpxvvtv（1）连续方程）连续方程（2）动量方程）动量方程-雷诺平均运动方程雷诺平均运动方程22)()(1)()(jiiijjjiiiixvvxppxvvvvtvv0jjxv221jiijjijjiixvxpxvvxvvtv，由连续方程
26、，由连续方程221jiijjiixvxpxvvtv2222)11)(jijiiijjijjijjijjiiixvxvxpxpxvvxvvxvvxvvtvtvjjijiijijixvvxvxpxvvtv)(1122湍流的湍流的雷诺平均运动方程雷诺平均运动方程与对应的层流运动方程相比，多了最后一项，该项中的与对应的层流运动方程相比，多了最后一项，该项中的jivv称为称为雷诺应力雷诺应力,是唯一的脉动量项，所以可以认为脉动量是通过雷诺应力来影响平均运动的。是唯一的脉动量项，所以可以认为脉动量是通过雷诺应力来影响平均运动的。222211jijiiijjijjijjijjiiixvxvxpxpxvvxv
27、vxvvxvvtvtv222211jijiiijjijjijjijjiiixvxvxpxpxvvxvvxvvxvvtvtv001100022jiijjijjiixvxpxvvxvvtv0jjxvjjijiijjijijixvvxvxpxvvxvvtv)(1122由连续方程由连续方程221jiijijixvxpxvvtvzpypzuyuxuxpxuutu1010)(1222222, 0, 0, 0tpxutu圆管中充分发展的层流和湍流圆管中充分发展的层流和湍流N-S方程张量形式方程张量形式（一）圆管中的层流（一）圆管中的层流层流中流体质点只有沿轴线的流动层流中流体质点只有沿轴线的流动u，而无横向
28、运动，所以，而无横向运动，所以v=w=0.假设管道水平放置，直径不大，管中具有一定压力，所以重力可以忽略。假设管道水平放置，直径不大，管中具有一定压力，所以重力可以忽略。流动恒定，流动恒定，u不随不随x和和t而变，只是而变，只是y和和z的函数，即的函数，即)(2222zuyudxdp所以所以，只有左右均等于同一常数才能成立。，只有左右均等于同一常数才能成立。lplppdp12常数21ppplpzuyu2222是长度为是长度为 l 的水平直管上的压降。的水平直管上的压降。所以所以因为管中流动是对称于因为管中流动是对称于x轴的，所以采用圆柱坐标系来分析更为方便，轴的，所以采用圆柱坐标系来分析更为方
29、便，22222222211rururruzuyu由于由于又因为速度又因为速度u的分布是轴对称的，所以的分布是轴对称的，所以0u0122lpdrdurdrud或或022lprdrdudrudr积分两次，可得积分两次，可得2214ClprlnrCu边界条件：边界条件：r=0时，时，u为有限值，得为有限值，得C1=0；r=d/2时时, u=0, 得得lpdC1622, 所以所以)4(422rdlpu圆管层流的圆管层流的速度分布规律速度分布规律，对称于管轴的，对称于管轴的抛物体。抛物体。（二）圆管中的湍流（二）圆管中的湍流湍流场质点间相互混杂、碰撞，导致运动状况极其复杂。对湍流的研究往湍流场质点间相互
30、混杂、碰撞，导致运动状况极其复杂。对湍流的研究往往是在某些特定条件下，对观测到的现象作某些假定，从而建立有局限性往是在某些特定条件下，对观测到的现象作某些假定，从而建立有局限性的半经验理论，再通过大量实验结果进行修正补充，得出湍流的半经验规的半经验理论，再通过大量实验结果进行修正补充，得出湍流的半经验规律。律。（1）脉动与时均流动）脉动与时均流动利用利用热线风速仪热线风速仪或或激光测速仪激光测速仪来测定湍流流速变化规律。来测定湍流流速变化规律。质点的真实流速质点的真实流速是无规律且瞬息万变的是无规律且瞬息万变的，这种现象称为，这种现象称为脉动脉动。每次实验在一个长的时间内平均后的速度值相同，为
31、每次实验在一个长的时间内平均后的速度值相同，为时均时均值。当湍流场中值。当湍流场中任一空间点上的运动参数的时均值不随时间（这里的时间是指湍流流动的任一空间点上的运动参数的时均值不随时间（这里的时间是指湍流流动的某一过程，而不是时均参数定义中所选定的某一很小的时间段某一过程，而不是时均参数定义中所选定的某一很小的时间段T）变化时，）变化时，称为称为定常湍流流动定常湍流流动，或称为，或称为准定常湍流准定常湍流，否则称为，否则称为非定常湍流非定常湍流。时均法只。时均法只能用来描述对时均值而言的定常湍流流动。能用来描述对时均值而言的定常湍流流动。注意：注意：时均化时均化的概念及准定常湍流流动，完全是人
32、为提出的一种模型，而的概念及准定常湍流流动，完全是人为提出的一种模型，而湍流实质是非定常的。因此在研究湍流的物理实质时，如研究湍流切应力湍流实质是非定常的。因此在研究湍流的物理实质时，如研究湍流切应力及湍流速度分布结构时，就必须考虑脉动的影响。及湍流速度分布结构时，就必须考虑脉动的影响。(2) 湍流流动中的附加切应力湍流流动中的附加切应力-雷诺应力雷诺应力粘性产生的切应力粘性产生的切应力 + 因质点混杂而形成的附加切应力因质点混杂而形成的附加切应力dAdtvmdAuvFFmdtFdAdtuvmuImvudAFtvut0vut0vutvudydUdAuvFm脉动流入脉动流入b层的流体质量层的流体
33、质量动量变化量动量变化量 = 切向力的冲量切向力的冲量脉动流体所受的脉动切向力脉动流体所受的脉动切向力脉动流体脉动流体m对对b层流体的脉动切向力层流体的脉动切向力a、b两层流体之间的脉动切应力两层流体之间的脉动切应力雷诺应力时均值雷诺应力时均值从上往下脉动时，雷诺应力大于零从上往下脉动时，雷诺应力大于零从下往上脉动时，雷诺应力大于零从下往上脉动时，雷诺应力大于零湍流运动中的湍流运动中的总的切应力总的切应力混合长度理论示意图混合长度理论示意图（3）普朗特混合长度理论）普朗特混合长度理论把湍流脉动与气体分子运动相比拟。（把湍流脉动与气体分子运动相比拟。（涡粘假设涡粘假设）粘性切应力由粘性切应力由分
34、子动量交换分子动量交换引起；假定脉动引起的附加切应力也为相同形式：引起；假定脉动引起的附加切应力也为相同形式：dydUdyduttdydUldydUlUUlyUyUU)()()(1dydUlUUu)(2121混合长度理论在于建立湍流运动中的附加切应力与时均流速混合长度理论在于建立湍流运动中的附加切应力与时均流速U之间的关系。之间的关系。引入了一个与分子平均自由程相当的长度引入了一个与分子平均自由程相当的长度 l，质点在走了，质点在走了 l 长度后与新位置长度后与新位置的质点掺混，完成动量交换。的质点掺混，完成动量交换。混合长度理论假定：在混合长度理论假定：在y层处，由于流体质点的横向运动所引起
35、的层处，由于流体质点的横向运动所引起的x方向湍方向湍流脉动速度流脉动速度u的大小为：的大小为：dydUlUdydUlUyUlyUU)()()(2也称为也称为涡黏度涡黏度t当流体质点从上层或下层进入所讨论的那一层时，它们以相对速度当流体质点从上层或下层进入所讨论的那一层时，它们以相对速度u相互接相互接近或离开，由流体连续性原理可知，它们空出来的空间位置必将由相邻的流近或离开，由流体连续性原理可知，它们空出来的空间位置必将由相邻的流体质点来补充，于是引起流体的横向脉动体质点来补充，于是引起流体的横向脉动v，两者相互关联，因此，两者相互关联，因此u与与v的大的大小必为同一量级，故：小必为同一量级，故
36、： ， vuvucv vuvu22)(dydUclvu22)(dydUlvudydUdydUldydUlvut222)(dydUlt2湍流粘性系数湍流粘性系数脉动引起的附加切应力脉动引起的附加切应力(雷诺切应力雷诺切应力)一般来说，混合长度不是常数一般来说，混合长度不是常数横向脉动和纵向脉动符号相反横向脉动和纵向脉动符号相反:代入相关式子，得代入相关式子，得将将 c 归并到尚未确定的归并到尚未确定的 l 中去中去可表示成可表示成，简单剪切湍流简单剪切湍流近代湍流研究的重大进展之一：发现剪切湍流中存在近代湍流研究的重大进展之一：发现剪切湍流中存在拟序结构拟序结构。湍流脉动并非完全不规则的随机过程
37、，而是在不规则的脉动中包含可辨认湍流脉动并非完全不规则的随机过程，而是在不规则的脉动中包含可辨认的有序大尺度运动。的有序大尺度运动。剪切湍流：剪切湍流： 壁湍流壁湍流，即壁面附近的剪切湍流，例如槽道、圆管、边界层湍流；，即壁面附近的剪切湍流，例如槽道、圆管、边界层湍流； 自由剪切湍流自由剪切湍流，例如射流、混合层、远场尾流。，例如射流、混合层、远场尾流。简单剪切湍流：平均剪切流动是简单剪切湍流：平均剪切流动是平行流动平行流动或或准平行流动准平行流动。研究雷诺数很高的壁湍流，研究雷诺数很高的壁湍流，/ReHUmH是直槽宽度之半是直槽宽度之半, 或圆管半径或圆管半径, 或边界层的平均名义厚度或边界
38、层的平均名义厚度, Um是平均特征速度。是平均特征速度。设定坐标系：设定坐标系：x为流动方向，为流动方向，y为垂直壁面方向，为垂直壁面方向，z为平均流动的展向为平均流动的展向进一步假定直槽沿展向无限长，流向单位长度上的平均压降是常数，则进一步假定直槽沿展向无限长，流向单位长度上的平均压降是常数，则 平均运动是定常的平均运动是定常的单向平行直线运动单向平行直线运动：=U(y)i1 脉动速度场在脉动速度场在流向和展向都是统计平均流向和展向都是统计平均的的 壁面上的速度等于零壁面上的速度等于零：无论是平均速度还是脉动速度都为零：无论是平均速度还是脉动速度都为零1、平均运动方程、平均运动方程0122y
39、vuyUxpypyv120zp(4.1a)(4.1b)(4.1c)(4.1b)(02xPvpdxdPxp0(4.1a)CydxdPdydxdPvuyU00当当 y=0时时, C= 0(壁面切应力壁面切应力)ydxdPvuyU00是分子粘性应力和雷诺应力之和，称为是分子粘性应力和雷诺应力之和，称为总切应力总切应力，它是，它是y的线性函数。的线性函数。在槽道的对称轴上在槽道的对称轴上( y = H )，由于平均运动的对称性，分子粘性应力和雷，由于平均运动的对称性，分子粘性应力和雷诺应力都等于零，于是有诺应力都等于零，于是有HdxdP00)1 (0HyvuyU在壁湍流中用壁面切应力定义壁湍流的速度尺
40、度，称为在壁湍流中用壁面切应力定义壁湍流的速度尺度，称为壁面摩擦速度壁面摩擦速度/0u/ReHu/ReHUmmUuReRe )10(21215 . 022220OcUucUuUcfmfmmf平板平板 Cf 范围（范围（0.0030.006）(3.95.5) 10-22、等切应力层、等切应力层1/,0HyyvuyU等切应力层可以进一步分为：等切应力层可以进一步分为：线性底层线性底层+对数层对数层(1) 线性底层线性底层(粘性底层粘性底层)0yUyuuU yU(2) 对数层和对数律对数层和对数律（等雷诺应力层）（等雷诺应力层）ReyHuHyyuy下面讨论近壁下面讨论近壁等雷诺应力层等雷诺应力层中的
41、统计特性中的统计特性)(ykvkvpyyUvu 在贴近壁面区，湍动能在贴近壁面区，湍动能耗散耗散和和扩散扩散相平衡；相平衡； 在稍离壁面且远离中心的流动区域中，在稍离壁面且远离中心的流动区域中，扩散项几乎可以忽略扩散项几乎可以忽略， 生成项生成项和和耗散项耗散项相平衡相平衡yUvu 在壁湍流中存在一个湍动能生成和耗散相平衡的区域。在壁湍流中存在一个湍动能生成和耗散相平衡的区域。由于平衡区远离中心区，可以用壁面参数表示速度梯度、雷诺应力、湍动由于平衡区远离中心区，可以用壁面参数表示速度梯度、雷诺应力、湍动能耗散率的无量纲式如下：能耗散率的无量纲式如下：)(2yguvu)(3yhyuyuyU2uv
42、uyu3(4.9a)(4.9b)(4.9c)(4.8a)(4.8b)(4.8c)积分式积分式(4.9a)，得到平均速度的对数分布，得到平均速度的对数分布(湍流的壁面律湍流的壁面律):ClnyuU)(yfyuyUBlnyu1（1）粘性底层内，流体质点没有混杂，故切应力主要为粘性切应力，附加）粘性底层内，流体质点没有混杂，故切应力主要为粘性切应力，附加切应力近似为零。粘性底层内速度梯度可认为是常数。切应力近似为零。粘性底层内速度梯度可认为是常数。yUdydUyw,22)(,dydUlyt*uw，它具有速度的量纲，称为，它具有速度的量纲，称为壁面摩擦速度壁面摩擦速度，则，则yReyuuUuyU*0*
43、22211)(ClnykuUydykudUdydUykw（2）粘性底层外，湍动剧烈，粘性影响可以忽略不计。）粘性底层外，湍动剧烈，粘性影响可以忽略不计。普朗特假设在近壁处混合长度普朗特假设在近壁处混合长度 l 与离壁面的距离与离壁面的距离y成正比，即成正比，即 l = ky。根据尼。根据尼古拉兹的实验证明，这个规律可以扩展到整个湍流区域。此外还假设在整古拉兹的实验证明，这个规律可以扩展到整个湍流区域。此外还假设在整个湍流区内切应力也为常数个湍流区内切应力也为常数 ，则，则w设设湍流边界层结构图湍流边界层结构图(1) 线性底层，线性底层， 该层内粘性应力远大于雷诺应力。该层内粘性应力远大于雷诺应
44、力。yuy*yyuuyuuyuyuyyuuuuww*2*当地雷诺数当地雷诺数, u+是量纲为是量纲为1的速度，则的速度，则说明速度随说明速度随 y 线性变化，所以称线性变化，所以称线性底层线性底层，由实验得出，该层范围为，由实验得出，该层范围为5yw*0ln1uCC2 .12, 2 . 5ln44. 22 .12,yyuyyuClnyu1(2) 对数律层对数律层由实验结果，当由实验结果，当y+40以后，雷诺切应力与壁面切应力以后，雷诺切应力与壁面切应力 大致相等且近似大致相等且近似为常数，可见粘性切应力可以忽略，其速度分布为为常数，可见粘性切应力可以忽略，其速度分布为式中的式中的C为常数，对光
45、滑壁为常数，对光滑壁C约为约为5.0-5.2,为卡门常数，一般取为卡门常数，一般取0.4-0.41。上。上式说明速度随式说明速度随 y 的增长呈对数关系增长，这就是对数率，满足对数率关系的的增长呈对数关系增长，这就是对数率，满足对数率关系的区域也成为区域也成为对数律层对数律层。在粘性底层和对数律层之间，平均速度分布既非线性的，也非对数的，因为在粘性底层和对数律层之间，平均速度分布既非线性的，也非对数的，因为这里分子粘性应力和雷诺应力属同一量级。介于粘性底层和对数律层之间的这里分子粘性应力和雷诺应力属同一量级。介于粘性底层和对数律层之间的流动区域成为流动区域成为过渡层过渡层，过渡层很薄，工程实用
46、上，常常不计过渡层，而用线，过渡层很薄，工程实用上，常常不计过渡层，而用线性分布和对数律分布组合成内层的平均速度分布。对于直槽湍流，可应用如性分布和对数律分布组合成内层的平均速度分布。对于直槽湍流，可应用如下的平均速度分布：下的平均速度分布：，对于高雷诺数绕流的对于高雷诺数绕流的雷诺平均雷诺平均求解方法，近壁湍流边界层很薄，实际数值计求解方法，近壁湍流边界层很薄，实际数值计算时，壁面网格只能达到等应力区外缘。另一方面，从壁面到等应力区的边算时，壁面网格只能达到等应力区外缘。另一方面，从壁面到等应力区的边缘（缘（y+=30）湍流统计量有剧烈的增加，任何数值方法都无法在一个网格中）湍流统计量有剧烈
47、的增加，任何数值方法都无法在一个网格中近似这种急剧变化。近似这种急剧变化。这时只好放弃数值积分到真实壁面，而是在离开壁面的第一层网格上用壁面这时只好放弃数值积分到真实壁面，而是在离开壁面的第一层网格上用壁面函数作为边界条件，或者说，将雷诺方程和近壁等应力层做渐近衔接，这时函数作为边界条件，或者说，将雷诺方程和近壁等应力层做渐近衔接，这时需要用到需要用到壁面函数壁面函数。Clnyu1壁面律壁面律推导近壁平均速度对数分布律的理论依据是存在等切应力层；而且在雷诺数推导近壁平均速度对数分布律的理论依据是存在等切应力层；而且在雷诺数很大时还存在等雷诺应力层。只要壁面切应力为有限值，上式对于任意平行很大时
48、还存在等雷诺应力层。只要壁面切应力为有限值，上式对于任意平行于壁面的湍流运动都适用，因此上式称为于壁面的湍流运动都适用，因此上式称为湍流的壁面律湍流的壁面律。计算中可以用上式的壁面律代替固壁无滑移条件，即将计算域的第一层网格计算中可以用上式的壁面律代替固壁无滑移条件，即将计算域的第一层网格设置在等应力层中，用上式作为边界条件。壁面剪应力的特征量是摩擦速度，设置在等应力层中，用上式作为边界条件。壁面剪应力的特征量是摩擦速度，它隐含于边界条件中，在数值求解中通过迭代求出。必须指出，上述壁面律它隐含于边界条件中，在数值求解中通过迭代求出。必须指出，上述壁面律只适用于附体边界层。只适用于附体边界层。
49、当壁面摩擦速度很小时，既要求当壁面摩擦速度很小时，既要求1y1y00，又要求，又要求的条件不能满足，的条件不能满足，这时壁面律就不再成立。特别是，接近流动的分离点处，壁面切应力这时壁面律就不再成立。特别是，接近流动的分离点处，壁面切应力0*u，即，即，不能应用以上壁面律来计算平均流速分布。，不能应用以上壁面律来计算平均流速分布。湍流数值模拟方法湍流数值模拟方法 直接数值模拟直接数值模拟(DNS, Direct Numerical Simulation) 雷诺平均数值模拟雷诺平均数值模拟(RANS, Reynolds Averaged Navier-Stokes) 大涡数值模拟大涡数值模拟(LE
50、S, Large Eddy Simulation)雷诺应力的封闭模式不可能是封闭的，而小尺度脉动对大尺度运动的雷诺应力的封闭模式不可能是封闭的，而小尺度脉动对大尺度运动的统计作用可能是普适的。统计作用可能是普适的。直接数值模拟直接数值模拟直接数值模拟的意义：直接数值模拟的意义： 可以获得湍流场的全部信息，而实验测量只能提供有限的流场分布可以获得湍流场的全部信息，而实验测量只能提供有限的流场分布 可以获得实时的流动演化过程，是研究湍流控制方法的有效工具可以获得实时的流动演化过程，是研究湍流控制方法的有效工具 可以评价已有湍流模型，研究改进湍流模型的途径可以评价已有湍流模型，研究改进湍流模型的途径
