计算流体力学课件：湍流及转捩.pptVIP

* 湍流扩散项 扩散速度与梯度呈正比 假设扩散系数为湍流粘性系数 湍流耗散项 模型1： 假设湍流耗散是各向同性的 模型2： 考虑各向异性 大小由速度脉动决定 该张量幅值为 分量大小由速度脉动均方根决定 6个自由变量的张量 标量 模型3： 各向异性 耗散率与Reynolds本身呈正比 更为合理 出现了湍能耗散率 （标量），需要单独给出方程 用前文给出的 方程即可 * Reynolds应力模型推导采用了更为理性的方法—— 更多的（严格）公式推导 更为复杂的公式 压力-变形项 （“再分配项”） 模化最为困难—— 压力-速度关联 实验测量困难， 数据少 DNS可能发挥很大作用 推导思路： 脉动压力的控制方程—— 压力Poisson方程 使Reynolds应力趋近于各向同性 复杂，计算量大， 工程应用不广泛 （直觉） “再分配”的特点—— 趋近各向同性 显然，涡粘模型 是一种最简单（各向同性）的代数模型。 ASM模型考虑了各向异性效应 * 最终，二阶矩模型（雷诺应力模型RSM）为： 仅是其中一种模型，还有其他模型 代数应力模型（ASM） 注： 需要与k方程、e方程联立求解 某些情况下，RSM的对流项与扩散项可忽略 高剪切流动： 生成项为主，对流、扩散项很小 局部平衡流动： 对流与扩散基本抵消 得到代数模型： Rodi部分保留了对流与扩散项，得到新的代数模型： 注： 仍需要与k方程、e方程联立求解 * § 11.5 压缩性对湍流模型的影响 压缩性效应： 与平均量有关的压缩性效应—— 外压缩性效应，非本质压缩性效应 与脉动量有关的压缩性效应—— 内压缩性效应，本质压缩性效应，声效应 Morkovin 假设： 当Mach数不是很高（例如平板边界层Ma4.5）时， 主要以外压缩性效应为主 压缩性- 近壁温度升高-密度降低-平均速度剖面改变 举例： 常用措施： 通过修正（例如密度加权平均）进行弥补—— 借用不可压缩的理论 常用措施：Favre 平均 密度加权平均 相对于Favre平均的脉动 * 对可压缩N-S方程进行（Reynolds）平均，令 密度和压力用Reynolds平均，其他量用Favre平均 可压缩N-S方程有大量 项利用Favre平均可简化方程推导 1） Reynolds平均可压缩N-S方程 平均 Reynolds应力 总能=内能+动能+湍能 层流热流 湍流热流 湍流与转捩 知识点： * 流动稳定性 转捩的预测方法 3. 湍流的模式理论（ RANS）: 涡粘模型： 0方程，1方程，2方程 4. 大涡模拟简介 * § 11.1 流动稳定性 一、 稳定性基本概念 常识：流体中的不稳定性 K-H不稳定性 A. K-H (Kelvin-Helmholtz）不稳定性—— 自由剪切流的无粘不稳定性 混合层—— K-H不稳定性 K-H不稳定性的关键： 速度剖面有拐点 Lee-Lin: 速度剖面的拐点是无粘不稳定性的必要条件 流体不禁搓，一搓搓出涡 已知某运动状态； 在此基础上施加微小扰动； 如扰动随时间（或空间）衰减，则称系统稳定，否则为不稳定 * 自然界中 K-H不稳定性图片 智利塞尔扣克岛的卡门涡街 澳大利亚Duval山上空的云 Kelvin–Helmholtz instability clouds in San Francisco 佛兰格尔岛周围的卡门涡街 高速流 低速流 自由剪切层受到扰动界面变形后的情况 K-H不稳定性的产生机理 受阻减速，压力升高，产生高压区 高压导致变形加剧 * B. T-S (Tollmien-Schlichting) 不稳定性——不可压 壁面剪切流的粘性不稳定性 Mack 不稳定性—— 超声速壁面剪切流的不稳定性 不可压边界层速度剖面 （Blasius解）—— 无拐点 可压缩情况—— Mach数足够高时会出现广义拐点—— 出现无粘不稳定性 不可压缩无粘不稳定性—— 需存在拐点 可压缩无粘不稳定性——需存在广义拐点 Mach 6 钝锥（1°攻角） 不同子午面 的分布 超音速平板边界层的不稳定波 第1模态（T-S波） 第2模态 （Mack模态） * 激波 密度界面 R-M (Richtmyer-Meshkov)不稳定性 —— 激波与密度界面作用的斜压效应 惯性约束聚变（ICF）示意图 小知识—— 涡的产生机制： 粘性、 斜压、有旋的外力 激波 密度界面 斜压项 * D. R-T (Reyleigh-Taylor)不稳定性 —