计算流体力学(中科院力学所)_第14讲-湍流及转捩2

1、1计算流体力学讲义计算流体力学讲义 第十四讲第十四讲 湍流与转捩湍流与转捩 （2） ；力学所主楼；力学所主楼219； 82543801 知识点：知识点： 讲义、课件上传至讲义、课件上传至 (流体中文网）流体中文网） - “流体论坛流体论坛” -“ CFD基础理论基础理论 ”讲课录像及讲义上传至网盘讲课录像及讲义上传至网盘 湍流的模式理论（湍流的模式理论（ RANS） 湍流的大涡模拟湍流的大涡模拟(LES)2知识回顾知识回顾1. 流体力学中的不稳定性流体力学中的不稳定性Kelvin-Helmholtz; Tollmien-Schlichting; Mack ; Richtmyer-Meshkov

2、 ; Reyleigh-Taylor ; Barnard 扰动的线性化控制方程扰动的线性化控制方程 不可压平行流不可压平行流 Orr-Sommerfeld 方程方程O-S方程的求解方程的求解 差分法差分法 （Malik的紧致差分法）的紧致差分法）vyuyuivyRe2222222220)(xBA全局法：全局法： 解出全部特征值解出全部特征值局部法局部法0)(BAF3 14.1 湍流的工程模式理论湍流的工程模式理论RANS 1. 为什么用湍流模型为什么用湍流模型N-S方程适用于湍流，但其解过于复杂方程适用于湍流，但其解过于复杂如果网格分辨率不够，数值解误差较大如果网格分辨率不够，数值解误差较大常

3、用方法常用方法 进行平均，求解平均进行平均，求解平均量满足的方程量满足的方程以不可压缩为例研究以不可压缩为例研究 - 推广的可压缩情况推广的可压缩情况压缩折角流动例例1： 压缩折角流动：压缩折角流动： 如果网格分辨率如果网格分辨率不足，且不用湍流模型，则分离区过大不足，且不用湍流模型，则分离区过大例例2： 有攻角机翼流动，如果分辨率不有攻角机翼流动，如果分辨率不足，且不用湍流模型，则造成足，且不用湍流模型，则造成“非物理非物理分离分离”翼型绕流4jjiijjiiiixxuxpxuutuxu210平进行均平进行均2. Reynolds平均的平均的N-S方程方程以不可压缩以不可压缩N-S方程为例方

4、程为例jjiijjiiiixxuxpxuutuxu210时间平均；时间平均； 空间平均；空间平均； 系综平均系综平均uuuii脉动ijjijijijjiijiuuuuuuuuuuuuuu)(引入平均引入平均：0jiuuRANS比比N-S方程多了该项方程多了该项jiRuu称为称为Reynolds应力应力)(Re1102jijjjiijjiiiiuuxxxuxpxuutuxu5)(102jijjjiijjiiiiuuxxxuxpxuutuxu3. Reynolds平均平均N-S方程的求解方程的求解 Rijiju u 未知量，必须用已知量表未知量，必须用已知量表示才能求解示才能求解湍流模型湍流模型方

5、法方法 1） Boussinesq 涡粘假设涡粘假设 （常用）（常用） kxuxuuuijijjitji32jjitijjiiiixxuxpxuutuxu2)(10与原先方程的唯一区别：与原先方程的唯一区别：改变了粘性系数改变了粘性系数程序实现方便程序实现方便 的计算模型：的计算模型： 0方程（代数模型）：方程（代数模型）： B-L 1方程模型方程模型 : S-A 2 方程模型方程模型 ：tkkSST方法方法2） Reynolds应力模型应力模型jiRuu给出给出的控制方程的控制方程可并入压力项中6 14.2 湍流边界层的结构及平均速度剖面湍流边界层的结构及平均速度剖面 内层1020% 外层钝

6、锥边界层的密度分布内层：内层： 主要受壁面影响主要受壁面影响外层：外层： 受边界层外部影响受边界层外部影响 （压力梯度、外部无粘流（压力梯度、外部无粘流） 速度剖面接近尾迹流速度剖面接近尾迹流内层的速度剖面内层的速度剖面221xyuuuuvputxyxyy 定常小量零压力梯度xyyyu220尾迹流的剖面（亏损律）零压力梯度的平板边界层零压力梯度的平板边界层)(yu7xyyyu220粘性底层区粘性底层区 湍应力可忽略湍应力可忽略5ywwxyyuconstyuyu yu/,/wwyyuu壁面律壁面律总应力保持不变总应力保持不变0)(xyxyy粘性应力粘性应力湍应力湍应力湍流核心区湍流核心区粘性应力

7、可忽略粘性应力可忽略wxyPlantdl混合长模型混合长模型yuyulxy2kyl wxy/wyukyCykuln141. 0kKarman常数C=5.5 （平板）（平板） 5.1壁面律壁面律对数律对数律亏亏损损律律过过渡渡区区8 14.3 常用的涡粘性模型常用的涡粘性模型 14.3.1 零方程模型零方程模型 Baldwin-Lomax (BL)模型模型术语术语“N方程模型方程模型”指计算湍流粘性系数指计算湍流粘性系数 时，使用了时，使用了N个偏微分方程个偏微分方程t零方程模型直接写出零方程模型直接写出 的表达式，简便的表达式，简便tBL模型是模型是 Plantdl混合长模型的推广混合长模型的

8、推广 湍流核心区湍流核心区 过渡区过渡区 粘性子层区粘性子层区壁面律壁面律对数律对数律亏亏损损律律过过渡渡区区内内层层外外层层混合长模型混合长模型kylyult,2)/exp(1 Aykyl0t内层统一内层统一表达式表达式尾迹亏损律尾迹亏损律近壁区趋近于近壁区趋近于0， 远壁区趋近于远壁区趋近于ky26A2)(lint涡量涡量9外层模型外层模型外层特点特点 1） 间歇性间歇性层流层流-湍流交替出现湍流交替出现Klebnoff间歇公式（根据实验间歇公式（根据实验得到的经验公式）：得到的经验公式）：16)/(5 . 51),(yyFKleb1为纯湍流，0为纯层流特点特点2） 类似尾迹流动的亏损律类

9、似尾迹流动的亏损律)/,min(max2maxmaxmaxFUyCFyFdifwkwake)/exp(1 )(AyyyFmaxymaxF yydifU为边界层内的最大速度与最小速度之差为边界层内的最大速度与最小速度之差difU内层外层1wkC10外层模型外层模型=间歇因子间歇因子*亏损律亏损律 wakeklebouttFF16)/(5 . 51),(yyFKleb边界层厚度不易计算16max5 . 51)(yyCyFklebKleb3 . 0klebC用用 估算边估算边界层厚度界层厚度maxyklebCy/max)/exp(1 )(AyyyFmaxymaxF yy内层3 . 0内层及外层的设定

10、内层及外层的设定yt2)(lint wakeklebouttFF内层外层1114.3.2 一方程模型一方程模型 k-方程模型方程模型1） 湍动能方程湍动能方程 （不可压缩）（不可压缩）iiuuk21jjiijjjjijijjupuuuxxkxxuuuxkutk121推导方法推导方法获得扰动量获得扰动量 的方程，的方程，两端乘以两端乘以 并平均即可并平均即可iuiu生成生成 耗散耗散 粘性扩散粘性扩散 湍流扩散湍流扩散jijixuxu2） k方程模型方程模型kxuxuuuijijjitji322kctjjx扩散型dsdVxjjn湍流扩散湍流扩散以湍流粘性系以湍流粘性系数进行的扩散数进行的扩散对于

11、一方程模型（k-模型）内部不会产生，也不会消失11. 009. 0clk2/3由量纲分析得出由量纲分析得出jktjjijijjxkxxuuuxkutk/ijijS1k1214.3.3 一方程模型一方程模型 Spalart-Allmaras (S-A)模型模型构造原则：构造原则： 经验经验 + 量纲分析量纲分析jjjjxxPxutDtD湍流场中标量方湍流场中标量方程的一般形式：程的一般形式： 对流对流 生成生成 耗散耗散 扩散扩散 对流对流 扩散扩散 生成生成 耗散耗散四四大大机机制制假定湍流粘性系数假定湍流粘性系数 满足上述方程满足上述方程假设：假设： 1） 生成项与当地涡量成正比生成项与当地

12、涡量成正比剪切越强，湍流越强：剪切越强，湍流越强： 符合直观符合直观vvScPb1ijijS22/ )/(ijjiijuuuu1355. 01bc关键的参数关键的参数这是湍流模型的这是湍流模型的“主要矛盾主要矛盾”混合长模型也是这么混合长模型也是这么假设的假设的生成项最为关键，对湍流粘性系数的生成项最为关键，对湍流粘性系数的影响最大。影响最大。 更简单的模型（零方程），更简单的模型（零方程），仅保留了生成项仅保留了生成项感想：感想： 生成项很多情况下都是最重要的生成项很多情况下都是最重要的 工资是工资是“生成项生成项”， 消费是消费是“耗散项耗散项”，把钱给其他人（家人、亲朋等）是，把钱给其他

13、人（家人、亲朋等）是“扩散项扩散项” 显然显然“生成项生成项”是最重要的。是最重要的。222vfdkvSS132） 耗散项与到壁面的距离有关，越远耗散越小耗散项与到壁面的距离有关，越远耗散越小yn/mmTurbulent kinetic energy budget012-0.00100.001PTDx=-30.0 mm21dvfcww湍流粘性系数越大，耗散越大湍流粘性系数越大，耗散越大离壁面越近，耗散越大离壁面越近，耗散越大直观/ )1 (/2211bbwckcc3/2,41. 0,622. 02kcb3） 扩散项简单模化扩散项简单模化 以以 层流层流+湍流粘性系数湍流粘性系数 为扩散系数为扩

14、散系数kkxvxD1最终涡粘系数的最终涡粘系数的 的控制方程为的控制方程为kkbjkwwbjjxxcxxdvfcvScxvutv)(12211vgcgcfwww6366312262),(dkSvrrrcrgw14近壁修正近壁修正 保证近壁处湍流粘性系数快速衰减到保证近壁处湍流粘性系数快速衰减到0f01020304000.81v1v/v1vtfvv 衰减函数衰减函数1 . 7,/,131331vvvcvvcf衰减函数衰减函数 的图像的图像1vf最终的湍流粘性系数最终的湍流粘性系数15得到得到 方程后，对方程后，对 求导，乘以求导，乘以 并平均，得到并平均，得到 的方程的方程14

15、.3.4 经典的两方程模型经典的两方程模型 模型模型 k2kctjktjjijijjxkxxuuuxkutk/jijixuxujixu2iujxljljlllliljlijijlijilljlijliljijjxxuxpuxxxuxuxuxuxxuxuuxuxuxuxuxuvxutDtD222222生成耗散扩散主项：主项： 小涡拉伸小涡拉伸粘性耗散粘性耗散ltllilixxkcxuuukcDtD)/(22192. 19 . 1,45. 141. 1, 3 . 121cc近壁区仍需衰减处理近壁区仍需衰减处理 “低低Reynolds数数k- 模型模型”固壁边界条件：固壁边界条件：0, 0k与物理情

16、况不符，近壁需要特殊处理与物理情况不符，近壁需要特殊处理1614.3.5 k- 两方程模型两方程模型用涡量用涡量 方程代替湍流耗散率方程代替湍流耗散率 方程方程/kt由量纲分析得到jkjjijijjxkkcxxuuuxkutk2jktjjijijjxkxkxuuuxkutk*/k*09. 0* 的模化方程：的模化方程：jtjjijijjxkxxuuukxutw*2/2*wk生成 耗散 扩散40/3, 9/5固壁边界条件：固壁边界条件：0k21610y1y第第1个点到壁面的距离个点到壁面的距离k- 模型模型近壁准确性优于近壁准确性优于k- 模型；模型； 但外层预测准确性不如后者但外层预测准确性不

17、如后者1y固壁利用近壁点的信利用近壁点的信息确定涡量息确定涡量1714.3.6 k- SST (Shear-Stress-transport)两方程模式两方程模式近壁：近壁：k- 外层：外层： k- 写成统一的写成统一的 k- 形式形式),max(211Fakt222500,99. 02maxtanhyykFjwtjjjjjxkxxxkFPxut/)1 (2212jktjjijijjxkxkxuuuxkutk/*2wP222214,500,99. 02maxmintanhyCDkyykFky较小时趋近于较小时趋近于k- 模式模式y较大时趋近于较大时趋近于k- 模式模式兼具兼具k- 及及k- 模

18、式的优点，是目前应用最广泛的湍流模型之一模式的优点，是目前应用最广泛的湍流模型之一18 14.4 非涡粘模型非涡粘模型涡粘模型的基本假设涡粘模型的基本假设：kxuxuuuijijjitjit32实际使用时，经实际使用时，经常不考虑该项常不考虑该项不符合物理规律不符合物理规律：（涡）粘性是各向同性；（涡）粘性是各向同性； 雷诺应力为湍流脉动影响雷诺应力为湍流脉动影响各向异性各向异性 舍弃涡粘假设舍弃涡粘假设， 直接针对直接针对 构造模型，构造模型，更为合理更为合理jituu合理的不合理的不一定好用一定好用写出脉动量写出脉动量 的方程，乘以的方程，乘以 并平均，得到雷诺并平均，得到雷诺应力应力 的

19、控制方程：的控制方程：iujujiuuijjilijilijlljilljlililjljlijixuxupxuuuupuuuxxuxuxuuuxuuuuuDtD112生成耗散扩散压力-变形 “再分配”对湍能无影响，不同分量之间再次分配分别模化，即可得到分别模化，即可得到Reynolds应力模型（又称应力模型（又称“二阶矩模型二阶矩模型”）出现三阶统计矩出现三阶统计矩19jilijlljiltuupuuuxD1湍流扩散项湍流扩散项jilkljiltltuuxkcxuuxxD2扩散速度与梯度呈正比假设扩散系数为湍流粘性系数湍流耗散项湍流耗散项模型模型1： 假设湍流耗散是各向同性的假设湍流耗散是各

20、向同性的ijljlixuxu322模型模型2： 考虑各向异性考虑各向异性 大小由速度脉动决定大小由速度脉动决定kuuxuxujiljli222该张量幅值为该张量幅值为分量大小由速度脉分量大小由速度脉动均方根决定动均方根决定6个自由变量的张量个自由变量的张量标量 模型模型3： 各向异性各向异性 耗散率与耗散率与Reynolds本身呈正比本身呈正比kuuxuxujiljli2更为合理出现了出现了湍能湍能耗散率耗散率 （标量），需要单独给出方程（标量），需要单独给出方程 用前文给出的用前文给出的 方程即可方程即可20Reynolds应力模型推导采用了更为理性的方法应力模型推导采用了更为理性的方法 更

21、多的（严格）公式推导更多的（严格）公式推导 更为复杂的公式更为复杂的公式ijjixuxup1压力压力-变形项变形项 （“再分配项再分配项”）模化最为困难模化最为困难 压力压力-速度关联速度关联 实验测量实验测量困难，困难， 数据少数据少 DNS可能发挥很大作用推导思路：推导思路： 脉动压力的控制方程脉动压力的控制方程 压力压力Poisson方程方程 lmmlmlmlmlxuxuuuuuxxp222rdVxuxuuuuuxxpVlmmlmlmlml2412VijijijjilmmlijjimlmlijjirdVxuxuxuxuxuxuxxuuxuxup)2()1(*22411kuukCijjii

22、j321)1(kkijijijPPkC322)2(使使Reynolds应力应力趋近于各向同性趋近于各向同性复杂，计算量大，复杂，计算量大，工程应用不广泛工程应用不广泛（直觉）（直觉）“再分配再分配”的特点的特点 趋近各向同性趋近各向同性21显然，涡粘模型显然，涡粘模型是一种最简单（各向同性）是一种最简单（各向同性）的代数模型。的代数模型。ASM模型考虑了各向异性模型考虑了各向异性效应效应kkijijijjilikljililjljlijiPPkCkuukCxuukcxkuuxuuuxuuuuuDtD3232)(212liljljliijxuuuxuuuP最终，二阶矩模型（雷诺应力模型最终，二阶

23、矩模型（雷诺应力模型RSM）为：）为：仅是其中一种模型，还有仅是其中一种模型，还有其他模型其他模型代数应力模型（代数应力模型（ASM）注：注： 需要与需要与k方程、方程、 方程联立求解方程联立求解某些情况下，某些情况下，RSM的对流项与扩散项可忽略的对流项与扩散项可忽略高剪切流动：高剪切流动： 生成项为主，对生成项为主，对流、扩散项很小流、扩散项很小局部平衡流动：局部平衡流动： 对流与扩散基对流与扩散基本抵消本抵消得到代数模型：得到代数模型：03232)1 (212kkijijjiijPckuukCPCRodi部分保留了对流与扩散项，得到新的代数模型：部分保留了对流与扩散项，得到新的代数模型：

24、kkijijijjiijijkkjiPPkCkuukCPPkuu323232)(21注：注： 仍需要与仍需要与k方程、方程、 方程联立求解方程联立求解ijjijixuxukcuu222 14.5 压缩性对湍流模型的影响压缩性对湍流模型的影响压缩性效应：压缩性效应： 与平均量有关的压缩性效应与平均量有关的压缩性效应 外压缩性效应，非本质压缩性效应外压缩性效应，非本质压缩性效应 与脉动量有关的压缩性效应与脉动量有关的压缩性效应 内压缩性效应，本质压缩性效应，声效应内压缩性效应，本质压缩性效应，声效应Morkovin 假设：假设： 当当Mach数不是很高（例如平板边界层数不是很高（例如平板边界层Ma

25、 近壁温度升高近壁温度升高-密度降低密度降低-平均速度剖面改变平均速度剖面改变举例：举例：常用措施：常用措施： 通过修正（例如密度加权平均）进通过修正（例如密度加权平均）进行弥补行弥补 借用不可压缩的理论借用不可压缩的理论uduuwVD0/常用措施：常用措施：Favre 平均平均密度加权平均密度加权平均 相对于相对于Favre平均的脉动平均的脉动23对可压缩对可压缩N-S方程进行（方程进行（Reynolds）平均，令）平均，令TTTpppuuuiii , ,密度和压力用密度和压力用Reynolds平平均，其他量用均，其他量用Favre平均平均可压缩可压缩N-S方程有大量方程有大量 项利用项利用

26、Favre平均可简化方程推导平均可简化方程推导1） Reynolds平均平均可压缩可压缩N-S方程方程0iiuxt0iiuxt平均0iiuxt0iiutxjiijjijijiuuxxpuuxut Reynolds应力jiijijiijiijTjLjjjjuuuxuuuuqqxHuxEt 21TRRTpkuueEii2112iiHhuuk 1122iiiiuukuu 总能总能=内能内能+动能动能+湍能湍能/ PrLjpjTqcxTjjqu h 层流热流湍流热流242） Reynolds应力及常用模型应力及常用模型12233kijRijTijijijkuu uSkx 涡粘模型涡粘模型 涡粘系数涡粘

27、系数 可利用可利用B-L,S-A,k- , SST等模型计算等模型计算t3）能量方程的模化）能量方程的模化jiijijiijiijTjLjjjjuuuxuuuuqqxHuxEt 21PrTpTjTjcTqx Pr0.89 0.90T湍流热流湍流热流：湍流能量扩散：湍流能量扩散：1/2ijijiitkjkuu uux 以以 层流粘性层流粘性+湍流粘性湍流粘性 扩散扩散 直观直观25湍流模式理论（湍流模式理论（RANS）：）： 计算量较小，但普适性差，很难找到通计算量较小，但普适性差，很难找到通用的模型用的模型 14.6 湍流大涡模拟简介湍流大涡模拟简介原因：原因： 湍流脉动的多尺度性湍流脉动的多

28、尺度性 大尺度脉动：大尺度脉动： 受几何条件、外部因素影受几何条件、外部因素影响强烈。响强烈。 复杂、多态、强各向异性复杂、多态、强各向异性思路：思路： 小尺度脉动受平均流影响较小，更容易模化小尺度脉动受平均流影响较小，更容易模化ijjijixuxukcuu2大涡模拟（大涡模拟（LES）：）： 流动流动= 大尺度流动大尺度流动 + 小尺度脉动小尺度脉动直接求解直接求解通过模型，由大尺度通过模型，由大尺度量给出量给出26kEnergy spectrum10010110210-1010-910-810-710-610-510-410-310-2FE2FE1FF1k*(-5/3)大尺度区 惯性区 耗

29、散区可压均匀各向同性湍流的能谱可压均匀各向同性湍流的能谱受几何条件，外受几何条件，外部因素影响强烈，部因素影响强烈，只能直接求解只能直接求解受外部因素影受外部因素影响较弱，容易响较弱，容易模化模化27滤波滤波 a. 盒式滤波盒式滤波 b. 谱截断滤波谱截断滤波c. Gaussian型滤波型滤波 dfxGxf)(),()(1),(dxGotherwisexifxG02/1),(2233222211/ )()()(62/326),(xxxexG2/2/xx14.5.1 不可压缩湍流的大涡模拟简介不可压缩湍流的大涡模拟简介28设设 的滤波尺度为的滤波尺度为2. 滤波的性质滤波的性质gfgfA.若采用

30、若采用Box 滤波及谱截断滤波则滤波及谱截断滤波则：ff 令：fff则:, 0 fB. 若采用一般的滤波器则：若采用一般的滤波器则：ff ,gfgf如采用如采用Gaussian型滤波有如下性质型滤波有如下性质f 相当于相当于 尺度的滤波尺度的滤波2defddeefdefdexfxxx2222222222/)(32/32/)(6/)(632/)(6/)(632)(26)(6)(6)(dexfx212/)(62/316)(f293. 基本方程基本方程 20()1iiijiijiuxuuuputxx 20()1iiijiijiuxuuuputxx jijiijuuuu大尺度量满足的方程大尺度量满足的

31、方程20()1iiijijiijijuxuuuputxxx )(jijijijiuuuuuuuu滤波：滤波：亚格子亚格子Reynolds应力应力jiijjijijijjiijijijiijuuuuuuuuuuuuuuuuuuuu)(性质：性质：ijijijLuuuuff 由于通常情况下由于通常情况下LES亚格子亚格子Reynolds应力与应力与RANS的的Reynolds应力形式有所区别应力形式有所区别jiRijuu RANSLeonard应力应力特点：特点： 无需模型，可直接计算无需模型，可直接计算304. 亚格子亚格子Reynolds应力模型应力模型jijiijuuuu（1） Smagor

32、insky 模型模型)()(21)(2uSuSuSijij其中其中ijkkijjiijuuuS,32特点：特点： 模型简单，鲁棒性好模型简单，鲁棒性好缺点：缺点： 在层流区耗散过大，在近壁区不适用。在层流区耗散过大，在近壁区不适用。 需要衰减函数需要衰减函数 A. 基本模型基本模型隐式滤波隐式滤波) 1 ()()(22uSuSCijsij涡粘模型25,)/(exp(13AAyD常用的衰减函数：常用的衰减函数：算出算出 后，乘以该后，乘以该函数即可函数即可ij只需将原先的粘性系数 换成t22( )tsCS u31（2）相似模型）相似模型jijiijuuuu假设不同尺度对雷诺应力的贡献是相似的假设

33、不同尺度对雷诺应力的贡献是相似的将上式中的将上式中的 换成换成 得得iuiujijiijuuuu即相似模型即相似模型该模型预测雷诺应力的准确度有所提高该模型预测雷诺应力的准确度有所提高但该模型预测的雷诺应力偏低但该模型预测的雷诺应力偏低 小尺度小尺度uuu大尺度大尺度32（3） 梯度模型梯度模型采用采用Taylor分析的方法找出亚格子应力模型分析的方法找出亚格子应力模型若采用若采用BOX滤波滤波312123222222224( )()1()24kkkf xfdffOxx x)()(241)(.)(21)()(1)(1)(4222222 Oxfxfdxfxfxfdxfxf222222422241

34、1()()2424()1()()241()12jiijijijikjkkkkkijijkkkjjkkkkkuuuuuuuuxxxxuuu uOxxuuOxxxx )(2Ouuii33推导过程并不严密，高阶量推导过程并不严密，高阶量 为必是小量为必是小量从相似模型推导，可以得出同样的公式。从相似模型推导，可以得出同样的公式。)(4O缺点：缺点： 稳定性差稳定性差 Liu et al 1994 建议采用限制器：建议采用限制器：otherwiseuifcjiij001,22241()12jjijkkkkkuuOxxxx 34B. 动力学模型动力学模型采用二次滤波的方法建立亚格子应力模型采用二次滤波的

35、方法建立亚格子应力模型小尺度G-level F-levelGermano 恒等式：F-level 滤波滤波 滤波尺度为滤波尺度为 ，G-level滤波滤波 滤波尺度为滤波尺度为 FG-level滤波滤波：ffkfjijiijuuuujijiijuuuuT ()(1)ijijijijijTuuuuL特点：特点： 该量无需模型，可直接计算该量无需模型，可直接计算FG-level滤波FG滤波F滤波+ G滤波35（1） 动力学涡粘模型动力学涡粘模型 F-levelFG-levelijdijijdijijMCuSuSuSuSkCT) )()()()()(2222ijdijijijMCTLijijijij

36、dMMLMC 预测亚格子雷诺应力的准确性有所提高，改进了层流区及近壁过于耗散预测亚格子雷诺应力的准确性有所提高，改进了层流区及近壁过于耗散的情况。的情况。)()(22uSuSCijdij)()()(22uSuSkCTijdij涡粘系数涡粘系数C动态可调动态可调通过两次滤波，确定该系数通过两次滤波，确定该系数FG滤波，相当于用滤波，相当于用 进行滤波进行滤波k可直接计算，无需模型36(2) 动力学混合模型动力学混合模型 基本模型为相似模型与涡粘模型的混合模型基本模型为相似模型与涡粘模型的混合模型)()(2uSuSCuuuuijdjijiijijijdijijijHMCTL)(jijijijiijuuuuuuuuHijijijijijdMMHLMC)(37（3）动力学）动力学Clark模型模型基本模型为梯度模型与涡粘模型的混合模型基本模型为梯度模型与涡粘模型的混合模型)()(1212222uSuSCuuijdjkikkij)(121()(121222222jkikkjkikkijuukuukHijijijijijdMMHLMC)(385. 近壁处理近壁处理jijiijuuuu显然在近壁处亚格子雷诺应力应当趋于显然在近壁处亚格子雷诺应力应当趋于0，但很多模型却不满足该条件但很多模型却不满足该条件因此需要采用特殊处理（采用衰减函数）因此需要采用特殊处理（采用衰减函数