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

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

基于RANS/LES方法的超声速底部流场数值模拟

分别采用基于两方程k-ω剪切应力输运(SST)湍流模型的延迟DES(DDES)、更改的DDES(MDDES)和改进的DDES(IDDES)方法,并引入可压缩修正,结合三阶MUSCL-Roe和五阶WENO-Roe两种空间离散格式,针对超声速底部的复杂流动现象,开展了数值模拟研究。计算结果表明本文方法能够捕捉到超声速底部流动中丰富的湍流结构,通过分析计算结果对超声速底部的流动机理有了进一步的认识,为下一步的超声速底部流动减阻改进和雷诺平均NavierStokes/大涡模拟(RANS/LES)方法在非定常高可压缩性流动中的应用提供了参考。通过对比分析不同空间离散格式的计算结果研究了数值耗散对计算的影响,五阶WENO-Roe格式的计算结果与实验结果吻合良好;对不同RANS/LES混合方法的计算结果进行了对比分析,结果表明IDDES方法在近壁区的表现优于DDES和MDDES方法。     

沟槽面湍流减阻数值评估方法

为解决雷诺平均Navier-Stokes(RANS)方法模拟沟槽面减阻的实际应用问题,受Wilcox对粗糙度模化方法的启发,通过深入分析k-ω剪切应力输运(SST)湍流模型近壁区ω值的作用效果,发现在黏性底层内增加壁面ω值可使对应近壁区的湍动能、湍流黏度、雷诺应力均减小,这种对近壁区流动特性的作用效果与真实沟槽面一致。以经典的对称V型沟槽面(高度=间距)减阻实验数据为基础,通过引入减阻效果影响因子——压力梯度与偏航角,建立复杂流动条件下沟槽面几何尺寸到壁面ω值的模化函数方程,将数值模拟结果与实验结果对比,验证了模化沟槽面可以达到与真实沟槽面相同的减阻效果,并依此给出沟槽面减阻在工程实践中的应用方法,以C919翼身组合体巡航构型为例,完成其沟槽面减阻的总体设计与评估。     

用强耦合RANS方法模拟旋翼悬停流场

在旋转坐标系下,将Spalart-Allmaras(S-A)一方程湍流模型和Reynolds-averaged NavierStokes(RANS)方程耦合成一个新的RANS方程,并发展了基于多块重叠网格的强耦合RANS求解方法,用于直升机旋翼悬停流场的数值模拟.为了提高计算效率,针对多重网格方法在多块重叠网格上实施的困难,提出了一种基于重叠网格的多重网格实施方法.通过对Caradonna-Tung(C-T)和ONERA 7A旋翼悬停算例验证了发展的强耦合RANS方法和基于重叠网格的多重网格实施方法的有效性.研究结果表明:发展的基于重叠网格的多重网格方法有较高的计算效率,3层网格的加速比约为7.7;强耦合RANS法的计算精度明显高于传统的松耦合RANS方法,特别是在与阻力相关性能参数的预测中,强耦合RANS方法的预测结果更加精确.     

隐式紧耦合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°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

湍流模型离散精度对数值模拟影响的计算分析

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS),通过改变物面法向第一层网格间距,开展了剪切应力输运(SST)两方程模型不同离散精度的数值分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR 7301两段翼型和高速RAE2822翼型,研究内容主要包括湍流模型的二阶精度离散和五阶精度离散两种方式对收敛历程、边界层湍流黏性系数分布、边界层速度分布、压力系数分布以及气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于不同的第一层物面法向网格间距,湍流模型离散精度对低速绕流计算结果有比较明显的影响,对于高速小迎角附着流动计算结果影响不明显;相对于湍流模型二阶精度离散,湍流模型高阶精度离散网格敏感性较弱,具有更高的数值模拟精度,但收敛性略差。     

基于CFD降阶模型的阵风减缓主动控制研究

飞行器飞行时会受到大气紊流的影响,降低飞行品质。阵风减缓控制是改善飞行器飞行性能的关键技术。现有的阵风响应分析多以离散阵风为研究对象,对更加真实描述大气紊流的连续型阵风时域分析关注较少。采用成形滤波器方法将频域形式给出的大气紊流信号转换为时域信号。在跨声速区域内,利用系统辨识技术,基于计算流体力学(CFD)方法建立阵风激励下的气动载荷状态空间降阶模型(ROM)。为方便控制器设计,借助平衡模态法进行模型的进一步降阶。使用模型预测控制(MPC)算法通过控制操纵面偏转实现阵风减缓主动控制。以AGARD445.6标模作为仿真算例,验证基于ROM设计的阵风减缓控制律的有效性。仿真结果表明,在跨声速飞行状态下,模型预测控制器能够在满足操纵面偏转范围的约束下,对连续阵风激励下的翼根弯矩输出进行有效抑制。     

发动机进气帽罩防冰热载荷的数值模拟研究

对某发动机进气帽罩防冰热载荷的计算方法进行研究。利用FLUENT软件对帽罩周围的流场进行模拟,同时采用FLUENT中的离散相模型计算帽罩周围的水滴运动轨迹,得到帽罩表面的水滴撞击效率;利用数值计算结果及流动换热经验关系式分析帽罩表面的各项热流,获得该发动机帽罩的防冰热载荷。计算中还比较了不同湍流模型对防冰热载荷的影响。     

湍流度和吹风比对叶片前缘双出口孔射流冷却效率影响

为了探讨湍流度对一种新型气膜孔射流气膜冷却影响,利用商业软件提供的有限体积法求解N-S方程,对湍流度分别为0.4%、10%和20%时的双出口孔射流冷却效率进行数值模拟。吹风比变化范围为0.5到2.0。首先将圆柱孔射流冷却效率计算结果和实验数据进行了对比,二者吻合较好。计算得到了双出口孔射流气膜冷却下的流场、径向平均冷却效率。结果表明,湍流度和吹风比对冷却效率都有较大影响。湍流度为0.4%和10%时,最高冷却效率在吹风比1.0时获得;湍流度30%时,最高冷却效率在吹风比2.0时获得。吹风比0.5和1.0时,冷却效率随着湍流度的增加而降低;吹风比2.0时,冷却效率随着湍流度的增加而增加。     

结构因素对离心通风器性能影响的数值研究

在高转速下,航空发动机离心通风器内腔存在强烈的湍流流动,而稳定工作后通风器内腔中漩涡呈周期性产生和淹灭。小尺寸颗粒受流体微团湍流随机脉动的影响,产生相对于平均流的随机脉动运动。为了研究通风孔偏心距和辐板顶圆半径对通风器性能的影响,采用离散相模型(DPM)模拟通风器内颗粒的运动轨迹。并应用随机游走(DRW)模型模拟连续相湍流瞬时速度脉动对颗粒轨迹的影响,采用随机涡寿命模型确定随机追踪模型的积分时间。结果表明:辐板顶圆半径的增大有助于提高通风器的分离效率,同时也增大了腔内流通阻力;通风孔偏心距对减小通风阻力的作用明显,但降低了离心通风器的分离效率。