增升装置微型涡流发生器数值模拟方法研究

 针对数值模拟带微型涡流发生器(VG)的增升装置绕流过程中出现的典型问题,研究相应解决途径.引入面搭接网格技术,避免传统结构网格方法造成的网格数量过大问题;在微型涡流发生器网格上采用特殊边界条件处理以分别对应有/无涡流发生器的状态,排除网格变动对计算结果造成的干扰;研发出了面搭接网格生成技术与特殊边界处理方法结合的求解加装微型涡流发生器的增升装置流动数值模拟方法.研究结果表明:该方法在处理增升装置这类复杂外形、复杂流动问题以及加装涡流发生器的流动控制问题时,具有良好的数值计算精度,能够满足研究需要,为开展涡流发生器精细设计和参数优化的数值模拟研究提供了可靠的技术途径.     

趋于临界马赫数的圆柱跨声速绕流特性分析

通过深化认识趋于临界马赫数Ma_cr的圆柱跨声速绕流特性,明确新型飞行器增升减阻设计的空气动力学理论依据。采用大涡模拟方法数值研究了来流马赫数Ma_∞为0.75和0.85、雷诺数Re为2×10⁵的圆柱跨声速绕流。结果表明:当Ma_∞趋于临界马赫数(Ma_cr≈0.9)时,圆柱的阻力下降且升力系数振荡被抑制;通过力的分解,得知圆柱的阻力减小来自旋涡力的影响,而非可压缩性;圆柱的阻力减小与其背压上升有关;剪切层初始阶段的对流马赫数Ma_c随Ma_∞的增加而增大,而增长率相反,这使得剪切层更为稳定、柱体背压更高。此外,由于Ma_∞=0.85时边界层分离点处的激波和尾迹处的激波向下游推移,使得近尾迹处的湍流脉动减弱,进而导致柱体的表面压力振荡和升力系数振荡被抑制。     

进气压力畸变试验中面平均紊流度的计算

为研究航空发动机进气压力畸变试验中数据处理方法对面平均紊流度计算结果的影响,通过分析涡扇发动机进气压力畸变试验数据,得到了不同取样时间及不同滤波截止频率对面平均紊流度计算的影响结果。分析表明,只有合理选择取样时间及滤波截止频率,才能保证所计算面平均紊流度的准确性。本研究为进气压力畸变试验中面平均紊流度的计算提供了一种参考方法。     

改进S-A湍流模型对横向射流的CFD模拟

通过在Spalart-Allmaras(S-A)湍流模型中加入应变率张量,来改善模型对横向射流(JICF)的计算精度。模拟了速度比分别为0.5和1.5工况下的横向射流。经实验值验证,结果表明:改进后的S-A模型预测精度有所提升;并对湍流模型改进部分对横向射流的作用机制做出了分析。将改进后S-A模型的计算结果与k-ε模型,雷诺应力模型(RSM)和大涡模拟(LES)的结果进行比较,发现S-A模型在速度比为0.5时精度最高,而在速度比为1.5工况下没有明显优势;并对这一现象做出了合理分析。      

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

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

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