高精度差分格式及可压缩湍流数值模拟（三）.pptVIP

三维全翼 Flow ONEAR-M6 三维翼 上表面的涡结构 三维全翼 Flow ONEAR-M6 三维翼 下表面的涡结构 流场中的拟序涡结构 光滑 条纹 Smoothing wall Ribletted wall 示例5： 条纹壁控制转捩 Mach 3, 6 湍流模型的评估与改进 后验评估 给定外形进行计算，与实验结果比对 先验评估 比RANS 更精细、准确的计算结果(DNS/LES) 作为“精确解”，验证模型 DNS给出 需要典型外形高分辨率湍流场 湍流模型是飞行器气动力热计算的关键之一 最新手段： 基于湍流数据库的先验评估 四、 基于DNS的湍流模型评估及改进 k-e模型： 1） 标准k-e模型 2）“可实现的”k-e模型 （realizable k-e model） Shi TH et al. Computer Fluids, 1995 根据实验结果 不是常数 3） 低Reynolds数k-e模型 近壁区进行了修正，增大了湍流粘性系数 评估方法： 使用DNS得到的湍动能k及湍能耗散率e 及计算 的模型公式 K-e模型的评估：湍流粘性系数及各模型的预测值 标准k-e模型预测值偏高 可实现的k-e模型、低Reynolds数k-e模型准确性较好 建议调整标准k-e模型中的系数Cm, 减小该值 k-e模型评估： Mach 2.25 平板 SA模型： 近壁衰减函数 对流 扩散 耗散 生成 模型系数： SA模型的评估: 预测值与DNS值的对比 近壁区吻合较好； 远离壁面区预测偏高 SA模型： Mach 2.25 平板 改进建议： 外层采用间歇函数减弱湍流粘性系数 内层 外层 BL 模型： 边界层= 内层 + 外层 混合长模型+近壁衰减函数 尾迹模型 （亏损律） 内层 外层： （间歇性） 边界层：内、外层示意图 （Mach 6 平板密度分布） 混合长+衰减函数 亏损律+间歇函数 根据最大值点的位置判断内-外层交界处 BL模型的评估： 模型预测值与DNS结果 （Ma=2.25） BL 模型： Mach 2.25 平板 近壁预测值偏低 外区外缘预测值偏高 改进建议： 内-外区划分算法 内层模型 外层模型 五、超声速射流燃烧DNS 1. Xin-liang Li, Yan Leng, Zhi-wei He. Optimized Sixth-order Monotonicity-Preserving Scheme by Nonlinear Spectral Analysis International Journal for Numerical Methods in Fluids, 2013; 73:560–57 2. J.Li, Z.Zhao, A. Kazakov, F.L.Dryer, An update decomprehensive kinetic model of hydrogen combustion, Int.J.Chem.Kinet.36(2004) 566–575. OMP6 格式[1] （6阶优化保单调） 9组分19反应模型[2] O， O2， H2, H2O, OH, H, HO2, H2O2, N2 粘性系数、 热传导系数、扩散系数采用 Chemkin拟合公式 计算网格 670*459*459 （1.4亿） 时间步长10-8 秒，每步5次化学反应推进； * 计算条件 射流速度： 904m/s ; 射流温度： 305K; 射流组分： 85% H2 + 15% N2 伴流速度: 20m/s 伴流温度： 1150K 伴流组分： 空气 21% O2 + 79% N2 超声速射流燃烧DNS 结果验证 x/D=8 及 x/D=12 处的平均速度剖面 1. Lu Shuqiang, Fan Jianren, Luo Kun, High-fidelity resolution of the characteristic structures of a supersonic hydrogen jet flame with heated co-flow air, Int. J. Hydrogen Energy, 37(4): 3528-3539, 2012 * 结果分析 1） 流场可视化 层流扩散火焰 湍流火焰 * 高精度差分格式及 可压缩湍流数值模拟 （三） 李新亮 中国科学院力学研究所 商飞北研培训讲座 2018-10-10 Part 3 可压缩湍流DNS 背景 OpenCFD软件简介 典型可压缩湍流的直接数值模拟示例 湍流模型的