CFD2020-第13讲-湍流与转捩2.pptVIP

速度散度云图 （近场声） 超声速射流燃烧DNS * 超声速射流燃烧DNS 动画演示： 速度散度 * 超声速射流燃烧DNS 动画： 生成物H2O分布 动画： 燃烧组分OH分布 * 超声速射流燃烧DNS 组分HO2分布 （自点火相关组分） H2+O2= HO2 + H * 流场中拟序涡分布 （Q等值面，温度染色） * * * * 计算流体力学讲义2020 第十三讲 湍流与转捩（2） 李新亮 知识点： * Copyright by Li Xinliang 1. 湍流的直接数值模拟 转捩过程中的涡动力学机制 转捩过程中的拟序涡结构 基于DNS的湍流模型评估及改进 § 13.1 湍流的直接数值模拟（DNS） 1. DNS: 直接求解Navier-Stokes方程， 无需引入任何湍流模型； 用足够的网格分辨率，分辨湍流的全部尺度细节； 湍流尺度： RANS LES DNS 积分尺度 >> 惯性尺度 >> 粘性尺度 粘性尺度 / 耗散尺度 / Kolmogorov 尺度 湍射流的涡量分布：DNS RANS 2. DNS 是研究湍流机理、 模型及控制的重要手段 无需使用湍流模型： 准确度高 分辨湍流全部细节： 精度高 DNS 高分辨率湍流数据库 湍流机理 湍流控制 湍流模型 湍流的产生、 拟序涡演化、 摩阻/热流/雷诺应力 抑制湍流、 利用湍流 评估模型、 改进模型 无粘流/Euler方程 大型飞机整机流场/气动力计算 1980 势流、面元法 粘性NS方程(RANS)、定常 1990 非定常RANS/DES 2000 2010 LES 2045 DNS 2080 年代 Splart 的 预期 ~m; (~10m) 最小尺度涡： 微米量级； 中小尺度涡： 0.1mm~ 1mm量级 直接数值模拟： 分辨出最小尺度涡； 网格量极为巨大： 计算量、存储量大 对高性能计算需求强烈 3. DNS对计算资源需求强烈 4. DNS 对数值方法要求苛刻 DNS 对数值方法要求： 高分辨率、 低耗散、 低色散 可压缩湍流DNS: 兼具低耗散及激波捕捉特性 矛盾： 低耗散 vs. 抑制振荡（需要耗散） 对于湍流DNS: 格式分辨率更为重要 Mach 6 有攻角钝锥边界层湍流DNS Mach 0.7-6 槽道湍流DNS RAE2822 超临界翼型三维绕流LES及DNS 压缩折角激波-边界层干扰DNS Mach 0.7-8 平板边界层湍流DNS § 13.2 典型可压缩湍流直接数值模拟示例 网格3.2亿 最大网格16亿 边界层湍流DNS, Ma=0.7, 2.25,6 扰动: 二维不稳定波 + 一对三维斜波 (自然转捩） 壁面吹吸扰动 (Bypass) 计算模型：平板边界层 示例1： 平板边界层湍流中的拟序结构 动画演示： 平板边界层拟序结构的形成及演化 流场中涡结构的瞬时演化图（Q=10的等值面） 1）拟序结构的形成及演化规律 应用示例： 平板边界层湍流中的拟序结构 发卡涡形成及发展的涡动力学机制 2） 拟序结构与切应力、壁面摩擦阻力 雷诺瞬时切应力集中在拟序结构周围。 新旧拟序结构干扰摩擦阻力剧烈增加。 发卡涡导致局部高剪切—— 诱发非线性Breakdown 瞬时雷诺切应力等值面 Ma=6 平板边界层 ：准流向涡而不是发卡涡 4) 压缩性效应对拟序结构的影响 Ma=0.7 （平板） Ma=0.7 （尖锥） Ma=2.25 （平板） Ma=6 （平板） 1）背景 飞行器的常用头部外形； 头激波 边界层转捩线 横流 头激波、熵层、横流； 高Reynolds数、高Mach数； 无成熟转捩（湍流）理论及模型； 转捩及湍流机理复杂； DNS: 研究该问题的有力手段 尚无他人DNS报道 示例2：有攻角小头钝锥边界层转捩的DNS 数值计算： 自主开发的程序(OpenCFD-SC) 直接数值模拟； 高精度差分法 对流项: 7阶精度WENO; 粘性项： 8阶精度中心; 时间项： 3阶精度Runge-Kutta Ren(头半径度量) Ma? a(攻角) Tw(壁温) T? 半锥角 10000 6 1度 294K