隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

SST-DES在小分离流动数值模拟中的改进

在基于剪切应力运输（SST,Shear Stress Transport）分离涡模拟（DES,Detached Eddy Simulation）应对小分离流动的数值模拟中,上游雷诺平均方法（RANS,Reynolds Averaged Navier Stokes）区域内较大的涡黏性通过对流扩散方程传递到下游大涡模拟方法（LES,Large Eddy Simulation）区域,导致自由剪切层解析湍流解迟迟得不到发展,剪切层失稳延后,出现强烈的灰区效应。为此将SST-DES中体现LES求解的类亚格子模式重新构造成标准k方程亚格子模式,防止多变量在RANS和LES交界面处对湍流黏性系数的混合影响,最终得到了一种改进的SST-DES方法。采用该方法对AS239翼型小分离数值模拟并与DES和延迟分离涡模拟（DDES,Delayed Detached Eddy Simulation）结果进行对比,可以看出该方法在壁面边界层保留了SST-DDES特性的同时,在尾迹区大大提高了LES的可解能力,显著降低了灰区的影响。     

湍流参数对风扇宽频噪声预测影响

对于现代先进涡扇发动机,占比越来越大的风扇/压气机宽频噪声成为亟需解决的问题。基于三维风扇声源模型,开展了湍流输入参数对宽频噪声预测影响的研究。结果表明:后传声功率级比前传声功率级普遍高5 dB左右;当湍流参数增大时,也会较大幅度改变前传和后传的声功率级,但对频谱形状影响较小;声功率级随湍流积分尺度增大而减小,随背景湍流速度、中心线湍流速度、湍流宽度增大而增大,其中背景湍流速度对宽频噪声影响最大。发现湍流参数通过作用于上洗速度谱来影响宽频噪声预测结果。      

基于代理模型的航空燃油泵空化性能分析

基于代理模型的整体敏感度分析方法,对影响燃油泵空化性能的结构参数（即诱导轮出口安放角和叶轮进口安放角）进行了分析,得到影响泵空化性能的主要因素,并对燃油泵的结构参数进行了优化。研究中,采用基于旋转修正的k-ε湍流模型及Zwart空化模型对优化前后的燃油泵空化流场进行了计算。结果表明：在一定变化范围内,叶片安放角对燃油泵的外特性影响较小,对燃油泵空化特性影响较大;其中,叶轮进口安放角对燃油泵空化性能有较大影响,随叶轮进口安放角的增大,燃油泵空化性能呈先下降后上升的变化趋势;优化后诱导轮出口安放角和叶轮进口安放角分别增大了4.4°和3.2°,满足Pareto最优解,燃油泵的空化性能较优化前提高了18%左右。     

UDMH/NTO火箭发动机尾焰流场特性数值仿真

采用AUSM(advection upstream splitting method)空间离散格式,探究了常用的3种k-ε湍流模型和2阶、3阶迎风格式对激波捕捉效果的影响。采用标准k-ε双方程模型和2阶迎风格式对不同高度和不同来流马赫数下的偏二甲肼/四氧化二氮(UDMH/NTO)火箭发动机尾焰流场进行喷管-尾焰流场一体化仿真。复燃反应采用12组分18步化学反应模型,喷管入口参数由热力计算给出。结果表明:随着来流马赫数的增加,波节数逐渐降低;随着飞行高度的上升,尾焰影响区域逐渐扩大;复燃反应造成混合区中的O2、N2质量分数大幅下降,而混合区中O、OH、NO质量分数则有所上升。      

不同湍流度下吹风比对涡轮导向叶片吸力面气膜冷却的影响

采用基于窄带热色液晶测量的瞬态全表面传热测量技术,研究了不同主流湍流度下的吹风比对涡轮导向叶片气膜冷却的影响,获得了叶片吸力面侧圆柱形孔排气膜冷却效率和表面传热系数比的全表面分布数据。结果表明:由于气膜射流与主流掺混的相互作用会随着主流湍流度的变化而变化,因此在主流湍流度不同时,吹风比对气膜冷却效率和表面传热系数比的影响规律会有所不同;主流湍流度较小时,吹风比的增大会显著减弱气膜覆盖效果与气膜冷却效率,但是在大湍流度下,吹风比的影响较弱,尤其是在远下游区域;相同的主流湍流度条件下,吹风比的增大会使得表面传热系数提高,但是在大湍流度下,换热增强效果较弱;相同吹风比下,高湍流度下的表面传热系数比相对较小。      

超声速燃烧中的特征尺度及影响因素

针对超燃冲压发动机典型的飞行条件,分别以氢气、乙烯和煤油(由质量分数为0.8的正十烷和0.2的三甲基苯化学替代)作为燃料,分析了超声速燃烧中的特征尺度及其影响因素,给出了细致的超燃冲压发动机工作范围,并探讨了火焰面模型在超声速燃烧数值模拟中的适用性.结果表明:从氢气、乙烯到煤油,超燃冲压发动机工作范围依次减少,丹姆克尔数(Da)呈现量级的变化,火焰模式以旋涡小火焰为主,其中Taylor尺度起着关键性作用.同时也发现:相对于亚声速燃烧,在超声速燃烧的数值计算中,对能否采用火焰面模型还需要更加仔细的考虑.      

GAO-YONG湍流模型对湍流传热的数值模拟

将GAO-YONG湍流模型应用于湍流传热的研究,分别计算了平板剪切湍流和二维平面冲击射流的湍流传热问题.边界层剪切湍流流动与换热的计算表明:与传统的湍流模型不同,GAO-YONG湍流模型不需要对近壁区域做任何特殊处理(比如壁面函数、低Reynolds数修正等)即可模拟出从壁面到主流区的全部流动与传热情况;另外,对于冲击射流Nusselt数的模拟也得到了与实验符合较好的计算结果,准确地捕捉到了2种冲击高度下流场换热的不同特征,表明了GAO-YONG湍流模型能够较高精度地计算湍流换热.      

用温度修正的k-ε模型研究排气系统红外辐射特性

为研究温度修正的k-ε湍流模型对排气系统红外辐射特性计算精度的改善程度,计算了某轴对称收敛喷管在3~5μm波段的空腔-喷流组合红外辐射特性,并与试验结果对比.喷管流场及温度场采用有限体积法(FVM)求解Navier-Stokes(N-S)方程得到,湍流模型分别采用标准k-ε模型和温度修正的k-ε模型.红外辐射特性的计算采用有限体积法求解吸收-发射性介质条件下的三维辐射传输方程,并考虑大气衰减作用.计算结果表明:温度修正的k-ε模型可有效预测排气系统热喷流的核心区长度,从而获得与试验结果一致的红外辐射强度.与标准k-ε模型相比,采用温度修正的k-ε模型计算得到的红外辐射强度的最大相对误差由64%降至了4.3%,说明温度修正的k-ε模型适用于精确模拟排气系统的红外辐射特性.      

旋转封严篦齿风阻温升的试验研究与数值分析

选取具有典型封严结构的真实尺寸台阶型封严篦齿,通过高速旋转条件下的风阻温升特性试验,进行了进出口压比(1.05~2.8)、封严间隙与齿尖厚度比(相对封严间隙)(2.4~4.0)、雷诺数(1900~28000)及转速(0~12000r/min)等参数对篦齿泄漏流风阻温升的影响研究.并基于试验条件,采用renormalization group(RNG)k-ε湍流模型,对同一篦齿的风阻温升特性进行了数值计算.结果表明:篦齿泄漏流温升随进出口压比、雷诺数、封严间隙与齿尖厚度比的增加而减小,但随转速的增加而增大.数值计算结果与试验数据吻合 良好,台阶型封严篦齿第1道和第2道封严齿内的风阻温升最为显著,前两道封严齿内的温升占总温升的60%. 根据试验结果,提出了计算台阶型篦齿泄漏流温升的经验公式,其计算结果与试验数据具有较好的符合性.