高超声速激波/边界层干扰流动数值模拟研究

在基于Roe格式的全Navier-Stokes方程计算流体力学(CFD)代码中,发展了一种局部熵修正方法,克服了传统熵修正方法在边界层流动模拟中耗散过大的缺点,可用于更加准确的模拟激波/边界层干扰的复杂高超声速流动。对典型高超声速双锥边界层分离与激波干扰的复杂流动进行了模拟,研究了网格密度和熵修正方法耗散性对计算结果的影响。研究表明:高超声速双锥边界层分离与激波干扰流动的数值模拟结果对网格具有很强的敏感性,过稀的网格将产生严重的分离流动预测偏差;低耗散性的局部熵修正方法能更加准确地模拟复杂的高超声速激波与边界层分离流动干扰现象。     

螺旋桨非定常滑流的高效数值模拟研究

基于自主研发的大型"亚跨超CFD软件平台"(TRIP3.0),开发了针对螺旋桨飞机的旋转机械非定常求解模块。控制方程采用曲线坐标系下的非定常Euler方程,综合运用动态拼接结构网格技术和双时间步推进方法,重点模拟了双发螺旋桨飞机在螺旋桨转动时的非定常滑流效应。为了保证计算模拟的高效性,引入多重网格技术和大规模并行计算技术。通过模拟某单独螺旋桨验证算例,验证了该方法的可行性。模拟了不同迎角和螺旋桨转速下全机的气动特性,分析了滑流对流场的影响情况,结果表明:螺旋桨滑流对全机流场和气动特性影响明显,并使得升力系数增大。     

CFD空间精度分析方法及4种典型畸形网格中WCNS格式精度测试

计算流体力学中的精度阶数体现的是数值解随着网格加密向精确解收敛的速度,计算格式、网格质量和边界处理等都会影响计算精度。5阶非线性加权紧致格式(WCNS-E-5)是一种具有较高的精度和分辨率的有限差分格式,对强激波和复杂网格的适应性强。但是,格式的模板自适应选择器和网格的复杂性等都会对最终计算精度造成影响。为了分析精度,首先给出了两种通过数值解反推计算格式精度的方法:一种方法适合已知精确解的情况,另一种方法适合未知精确解的情况。然后在斜交网格、拉伸网格、拐折网格和扭曲网格等4种典型畸形网格中考察了WCNS-E-5格式的精度,考察算例涉及线性波动问题、等熵涡传播问题和超声速边界层问题等。计算结果表明WCNS-E-5在这四种网格中都体现出高阶精度。     

基于N-S方程的剪刀式尾桨前飞状态气动力计算研究

建立了一个适用于计算剪刀式尾桨前飞状态气动力的数值方法。在该方法中,采用惯性坐标系下的非定常Navier-Stokes方程作为主控方程,使用低数值耗散的Roe格式进行空间离散,并应用高效的隐式LU-SGS格式进行时间推进;采用双时间法以模拟尾桨前飞时的非定常现象;使用运动嵌套网格方法以计入尾桨桨叶的旋转运动。应用建立的方法,首先对有试验数据可作对比的Lynx直升机尾桨悬停状态气动性能进行了计算,验证了所建立方法对尾桨气动性能计算的有效性;同时,针对Helishape 7A旋翼进行了前飞状态气动载荷计算,验证了该方法对旋翼(尾桨)前飞气动载荷计算的有效性。最后,对剪刀式尾桨前飞状态进行了气动力计算,并分析了构型参数对剪刀式尾桨气动力的影响。计算结果表明,前飞状态时,相同总距下剪刀式尾桨的平均拉力要大于常规尾桨,且剪刀式尾桨在90°-180°-270°方位角内的桨-涡干扰现象较常规尾桨严重。     

摆线桨非定常气动特性研究

以西北工业大学研制的滚翼飞行器为研究对象,对摆线桨二维简化模型进行了非定常数值模拟。非定常时间推进采用双时间法,计算模型的刚体运动和网格变形采用弹簧近似光滑模型和局部重划模型相结合的动网格技术来处理。重点对四桨叶摆线桨转动过程中的瞬时推力、桨叶气动力及其非定常流场进行了研究,结果表明:摆线桨转动一周中瞬时推力大小随转动方位角变化呈近似正弦规律的周期性波动;随着转速变化,摆线桨转动一周的平均推力方向几乎保持恒定;桨叶瞬时法向力和切向力随方位角变化呈明显的"∞"型迟滞环曲线,且桨叶上行状态的气动力明显大于下行状态气动力;摆线桨桨叶尾迹之间具有很强的非定常干扰,对桨叶气动性能会产生较大影响。     

基于Gama-Theta转捩模型的气冷涡轮气热耦合研究

基于非结构化混合网格,采用有限体积法求解了N-S方程,对流项格式为AUSM+,梯度的求解方法为加权最小二乘法,求解器采用LUSGS隐式求解方法,并进行了预处理和网格重排序。与导热方程耦合,搭建了气热耦合平台,采用直接耦合方法,交界面互相传递温度,并且采用能保证通量守恒的面积加权类的插值方式实现数据传递,将γ-Reθt两方程转捩模型应用在流场程序中。通过MARKII叶片4311和5411工况的实验结果进行了验证,结果表明采用转捩模型对压力分布的影响不大,与实验值吻合较好,通过间歇因子的分布可以看出,γ-Reθt模型成功地预测了从层流到湍流的转捩过程,在层流区采用γ-Reθt模型计算的涡粘系数与真实流动情况更加吻合,由于涡粘系数对传热的影响很大,从而有效地提高了边界层的层流区和转捩区的传热计算的精度,得到的温度和换热系数与实验值吻合更好,但是该模型是和SST模型耦合在一起的,由于SST模型的局限性,所以对激波和边界层干扰区域的模拟产生了7%左右的温度误差。     

基于疲劳分析的起落架形状优化方法

<正>目前,对飞机起落架的疲劳破坏特性和疲劳寿命估算有较多的研究,在结构优化方面多是考虑强度和重量之间的优化关系,而基于疲劳分析的结构优化方法在国内外虽有一些研究,但是有起落架背景的研究较少。结构疲劳优化方法需要解决通过建立相对接近实际有限元分析模型以及起落架各部件之间的连接,应用相对准确的疲劳分析方法对起落架寿命进行估算;在此基础上提出合理的形状优化方案,采用合理的优化方法,进行疲劳优化,以期获得更为精确的优化结果。     

多段翼型局部主动变形流动控制的非定常数值模拟

对30P30N三段翼型失速攻角附近的分离流动进行了数值模拟研究。为了抑制大攻角时背风区的流动分离,在主翼段上表面引入了行波壁变形模型和抛物型局部主动振动模型,利用作者以往发展的动态混合网格技术和相应的非定常计算方法,对变形过程中的非定常分离流动进行了数值模拟,分析了各种变形参数对流动分离的影响。计算结果表明,在适当的条件下,局部主动变形能够抑制翼型背风区的分离,由此可以起到增升减阻的作用,改善翼型的气动性能。     

霜状冰结冰数值模拟研究

提出了一种三维霜冰结冰的数值模拟方法。基于欧拉-拉格朗日法计算空气-水滴两相流流场,获得了表面水滴撞击特性以及撞击量,求解质量守恒和能量守恒方程获得网格单元内的结冰量,三维霜冰的结冰计算程序采用用户自定义（UDF）编写,FLUENT提供的动网格功能实现了结冰后的网格重构。计算了4°攻角,液态水含量为1 g/m3,水滴直径20μm,结冰时间6 min,温度分别为-19.4°和-28.3°时NACA0012翼型的结冰,并与相关文献结果进行了对比,计算结果与文献实验结果吻合较好,验证了方法的合理性和可靠性。分析了不同的粒径、液态水含量、来流速度等对撞击特性以及结冰的影响,随着粒径、来流速度的增加,水滴收集系数、撞击区和结冰厚度增大,随着液态水含量的增加,水滴收集系数和撞击区是不变的,但是结冰厚度增加。     

有限元分析前,后置处理系统

本文介绍了一个基于几何造型系统的有限元分析的前、后置处理系统。前置处理系统的主要功能是对几何造型系统的二维域、三维实体及曲面进行有限元网格剖分、网格结点自动编号，后置处理主要功能是利用图形技术显示有限元分析的输出结果。