超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

展向振荡对激波/湍流边界层干扰的影响

周期振荡作为一种有效的壁面流动控制手段受到广泛关注,而其对激波/湍流边界层干扰的影响目前鲜有研究。本文采用高精度直接数值模拟（DNS）方法对马赫数2.9、12°激波入射角、强振荡下的激波/湍流边界层干扰进行了系统研究。通过与无振荡工况的定量比较,揭示了展向强振荡对干扰区内复杂流动结构的影响规律及作用机制,如分离泡尺度、物面压力脉动非定常特性、物面剪切的非定常特性及统计特征等。研究发现：在展向强振荡作用下,分离点位置提前,间歇区长度增大;同时由于分离泡内强黏性耗散的影响,展向振荡的穿透高度约为分离泡高度的4%,因而对流动结构不会产生实质影响。但展向强振荡会对壁面附近流动造成显著影响,如强振荡诱导的壁面展向速度远大于流向速度,造成流向剪切与展向剪切之间夹角的概率密度函数峰值从0°偏移到80°~90°之间。物面压力及剪切本征正交分解分析表明,展向振荡会导致模态能量从低阶模态向高阶模态转移,降低低频运动的能量占比,增强再附后Görtler涡等壁面附近旋涡结构的强度。     

钝体稳定射流火焰的火焰面模型数值研究

基于稳态火焰面(SLFM)和交互式非稳态欧拉颗粒火焰面(EPFM)模型对Sydney大学CH4/H2钝体稳定湍流扩散火焰行了数值研究,采用修正的雷诺应力模型(RSM),同时对两种不同规模的甲烷详细化学反应动力学机理进行研究,比较了燃烧模型和反应机理对湍流火焰结构、活性自由基以及氮氧化物预测精度的影响,与实验数据对比结果表明:两种反应机理得到的温度场和主要组分分布基本相同;SLFM模型能对速度场和标量场的分布进行较为准确预测,采用EPFM模型修正后,部分区域OH预测结果更加靠近实验结果;采用EPFM模型对SLFM模型耦合GRI-Mech 211的计算结果修正后,NO量级降低近2倍,预测精度明显改善,与实验结果实现较好的符合,验证了化学反应动力学机理以及非稳态效应对氮氧化物预测精度的影响。     

超声速平板边界层旁路转捩直接数值模拟

采用具有七阶精度的有限差分格式和基于Euler方程组的特征边界条件,通过直接数值模拟研究了由固定相位周期性吹吸扰动诱导的、Ma∞=2.25、Re∞=2.5×107/m(635 000/inch)空间发展平板边界层的旁路转捩过程,并建立了相应的直接数值模拟数据库。数值结果表明:扰动相位的随机性并不是转捩发生的必要条件,固定相位周期性吹吸扰动即可诱发边界层的转捩;当采用固定相位时,随着吹吸频率的增加,转捩位置向下游移动;当扰动频率足够大时,周期性吹吸无法诱导出边界层的转捩。随后,在所得数据库中提取了速度及热力学变量的平均值,与已有的数据库和理论分析吻合良好;速度脉动均方根与公认的不可压缩平板边界层实验基本吻合;热力学变量脉动均方根和雷诺应力与已有数据库的统计结果及实验数据具有相同的分布特征。对比发现,不同转捩的诱导因素对平均转捩位置和热力学变量均方根、雷诺应力等高阶统计量的影响较大。     

无网格法耦合RNG k-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

基于移动最小二乘无网格方法,耦合RNG(Re-Normalisation Group)k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM(Advection Upstream Splitting Method)+-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP(Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNGk-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

γ-Reθ转捩模型在二维低速问题中的应用

介绍了基于当地关联的γ-Reθ转捩模型.通过Spalart建议的环境源项法对自由流湍流度衰减的有效控制,实现了该转捩模型的实用化.利用中国空气动力研究与发展中心的TRIP软件平台上实现的该转捩模型,对S809翼型、NLR7301两段翼型和30P30N三段翼型3个低速翼型绕流问题进行了计算,将不同迎角下的计算结果与试验结...     

可压缩湍流边界层近壁区马蹄涡的演化

采用大涡模拟方法对可压缩槽道湍流进行了数值模拟.通过给定扰动波方法确定进口边界条件,很快得到了完全发展的湍流.利用动态亚格子模型、预处理技术以及高阶格式和方法,得到了精确的数值结果.计算得到的速度剖面和脉动分量与直接数值模拟的解进行了对比,验证了计算结果的可靠性.结果清楚地展示了湍流场中马蹄涡的生成和演化过程,包括对称涡腿、单侧涡腿的马蹄涡结构,以及流场中存在由亚谐波引起的拟序结构的交错现象.分析了湍流边界层近壁区可能出现的复杂的马蹄涡结构,及其涡结构形状也不严格对称的现象.     

用于可压缩壁面湍流LES及DNS入口边界条件的"回收-调节"方法综述

为了减小湍流流动问题的大涡模拟和直接数值模拟的入口湍流脉动边界条件对计算结果的影响,叙述了用于可压缩壁面湍流入口脉动生成的一种方法--"回收-调节"方法(recycling-rescaling method),叙述了这种方法的理论基础及简要推导过程,介绍了几种常用的"回收-调节"方法并指出了实现中的几个关键问题,分析了...     

"回收/调节"方法在混合LES/RANS模拟方法中的应用

采用一种混合大涡模拟/雷诺平均Navier-Stokes(LES/RANS)模拟方法结合三阶加权基本无振荡(WENO)格式对马赫数为2.88的压缩斜坡流动进行了数值模拟,并采用"回收/调节"方法生成入口湍流脉动边界条件.当采用固定入口边界条件时,湍流边界层会缺乏合理的湍流能量,使之抵抗逆压梯度的能力减弱,会严重高估分离...