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

示例3：压缩折角DNS 背景： 高超声速飞行器、 激波-边界层干扰 PIV 实验 （易仕和） DNS ： 中截面上的瞬时温度 特点： 全区域计算（转捩+角区） 壁面平均压力分布 角部分离区内的速度剖面 DNS 结果与实验吻合理想，验证了结果的可靠性 瞬时壁面压力时间信号 (虚线： 瞬时信号； 实线：滤波后的信号） 特点： 双时间尺度特征明显 A. 激波振荡的双频特征及低频机制 DMD分解 高频： 剪切振荡 机制 低频： 分离泡 “呼吸”机制 B. 壁温效应对分离泡大小的影响 来流Mach数2.9 不同壁温的压缩折角流动 角部区域的流动 角部区域时均流动的流线 随壁温升高，分离泡变大 分离泡大小随壁温变化规律 理论分析: 结合DNS，给出半经验公式 与DNS吻合较好 （误差<5%） 转捩效应对角部分离区的影响 转捩中期 转捩峰 转捩后期 转捩起始 充分发展 计算域设置图 转捩中期 转捩峰 转捩后期 充分发展 中截面的数值纹影（密度梯度）图 特点： 来流为转捩中期阶段的流动时， 分离泡最小； 来流为充分发展湍流时，分离泡最大 近壁面流向速度分布图 case1 case2 case3 case4 case5 转捩早期： 强三维结构 充分发展湍流：相对均匀 转捩早期， 流动呈大尺度结构，流动抗分离能力强； 转捩后期， 流动破碎成小尺度结构，抗分离能力减弱 示例4： 三维全翼LES 翼面及 z=28.3% 截面上瞬时密度分布 ONEAR-M6 翼， Mach=0.8395 , 特点： 三维全翼 三维全翼 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 平板 改进建议： 外层采用间歇函数减弱湍流粘性系数 内层 外层 * 高精度差分格式及 湍流数值模拟 （三） 李新亮 中国科学院力学研究所 Part 3 可压缩湍流DNS 背景 OpenCFD软件简介 典型可压缩湍流的直接数值模拟示例 湍流模型的评估及改进 湍流燃烧的DNS 6. 小结 湍射流的涡量分布：DNS RANS 瞬时温度场 流动机理、流动控制 评估、改善湍流模型 一、 背景 直接数值模拟 （DNS） 是研究湍流机理、 模型及控制的重要手段 无需使用湍流模型： 准确度高 分辨湍流全部细节： 精度高 无粘流/Euler方程 大型飞机整机流场/气动力计算 1980 势流、面元法 粘性NS方程(RANS)、定常 1990 非定常RANS/DES 2000 2010 LES 2045 DNS 2080 年代 Splart 的 预期 ~m; (~10m) 最小尺度涡： 微米量级； 中小尺度涡： 0.1mm~ 1mm量级 直接数值模拟： 分辨出最小尺度涡； 网格量极为巨大： 计算量、存储量大 对高性能计算需求强烈 DNS对高性能计算需求强烈 二、 高精度CFD软件 OpenCFD OpenCFD: 作者开发的一套高精度、开放的CFD程序 1） OpenCFD-SC : 高精度差分 2） Open