商飞北研讲座2018李新亮part3.pptx

商飞北研培训讲座 2018-10-10高精度差分格式及 可压缩湍流数值模拟（三） 李新亮中国科学院力学研究所mailto:lixl@邮箱： mailto:lixl@lixl@/s/1slfC5Yl网盘： /s/1slfC5Yl/s/1slfC5Yl Part 3可压缩湍流DNS 背景 OpenCFD软件简介 典型可压缩湍流的直接数值模拟示例 湍流模型的评估及改进 湍流燃烧的DNS6. 小结一、 背景 直接数值模拟 （DNS） 是研究湍流机理、 模型及控制的重要手段 无需使用湍流模型： 准确度高 分辨湍流全部细节： 精度高流动机理、流动控制评估、改善湍流模型瞬时温度场RANS湍射流的涡量分布：DNS DNS对高性能计算需求强烈年代2080DNSSplart 的 预期 LES2045~m; (~10m)非定常RANS/DES2010最小尺度涡： 微米量级；中小尺度涡： 0.1mm~ 1mm量级 粘性NS方程(RANS)、定接数值模拟： 分辨出最小尺度涡； 网格量极为巨大： 计算量、存储量大 无粘流/Euler方程1980势流、面元法对高性能计算需求强烈大型飞机整机流场/气动力计算二、 高精度CFD软件 OpenCFDOpenCFD:作者开发的一套高精度、开放的CFD程序1） OpenCFD-SC :高精度差分2） OpenCFD-EC:有限体积3） OpenCFD-Comb: 化学反应计算1）高精度有限差分求解器OpenCFD-SCOpen CFD code for Scientific Computing; 面向科学计算开放的CFD代码算法： 高精度差分特点： 精度高、适用于湍流DNS及LES等 外形相对简单， 但流场复杂 OpenCFD-SC 的大规模测试及计算时间计算机计算规模（核）最大规模时并行效率问题描述2007.12国产并行机4096-5824050%Mach 6 钝锥转捩DNS，计算网格点3.5亿2010.12天河1A256-1638490%Mach 8 平板边界层湍流DNS， 网格点4亿；2013.4国产并行机32768-98万 （含从核）37.5%M6 三维翼湍流LES 2014.12天河-28192-24万65.14%M6 三维翼湍流LES2015.12国产并行机32768-380万（含从核）约60%M6 三维翼湍流LES OpenCFD程序 异构并行移植及优化 GPU系统 （天河-1号） MIC系统 （天河-2号）内存 国产众核系统 （神威太湖之光）2）多块结构网格有限体积程序—— OpenCFD-EC面向工程计算开放的CFD代码(Open CFD code for Engineering Computing)算法： 多块结构网格有限体积 + RANS特点： 适用于复杂外形工程计算差分-有限体积混合方法3） 高精度化学反应模拟程序—— OpenCFD-Comb面向化学反应的开放CFD程序： An OpenCFD code for Combustion高精度差分方法， 适合于化学反应湍流高精度模拟 Chemkin-like 化学反应机理接口 ； Chemkin 热力学模型、 物理模型 （粘性、热导、扩散系数） 高精度差分方法：OMP6, GVC, WENO7 高阶滤波，保证计算稳定性； MPI并行， 可扩展到百万CPU规模；三、 典型可压缩湍流直接数值模拟示例最大网格16亿Mach 0.7-8 平板边界层湍流DNSMach 0.7-6 槽道湍流DNS网格3.2亿RAE2822 超临界翼型三维绕流LES及DNSMach 6 有攻角钝锥边界层湍流DNS压缩折角激波-边界层干扰DNS示例1： 平板边界层湍流中的拟序结构计算模型：平板边界层 边界层湍流DNS, Ma=0.7, 2.25,6 扰动: 二维不稳定波 + 一对三维斜波(自然转捩） 壁面吹吸扰动 (Bypass) 应用示例： 平板边界层湍流中的拟序结构1）拟序结构的形成及演化规律动画演示： 平板边界层拟序结构的形成及演化流场中涡结构的瞬时演化图（Q=10的等值面）发卡涡形成及发展的涡动力学机制2） 拟序结构与切应力、壁面摩擦阻力 雷诺瞬时切应力集中在拟序结构周围。新旧拟序结构干扰摩擦阻力剧烈增加。 发卡涡导致局部高剪切—— 诱发非线性Breakdown瞬时雷诺切应力等值面3) 发卡涡头部脱落发卡涡发生断裂，头部成为孤立自由涡。自由涡脱离主体，沿边界层外层向下游翻滚运动。边界层外层的间歇现象与自由涡的脱落有关。 4) 压缩性效应对拟序结构的影响Ma=0.7 （平板）Ma=0.7 （尖锥）Ma=2.25 （平板）Ma=6 （平板） Ma=6 平板边界层 ：准流向涡而不是发卡涡示例