混合网格方法研究栅格翼亚跨超声速气动特性

采用基于结构/非结构混合网格的CFD方法对栅格翼的亚跨声速气动特性进行了研究,计算方法经过试验数据的验证,可达到工程精度。对单独栅格翼的研究表明,跨声速壅塞时,流动通过栅格前方亚声速气流减速从栅格外侧溢流而实现流量调节;超声速壅塞时,流动通过在栅格翼上产生网状脱体激波使来流减速至亚声速来进行流量调节,且计算结果与理论估算上下临界马赫数范围一致。对某栅格翼气动布局导弹的计算结果表明,在跨声速壅塞区间内,壅塞有自我调节作用,通过溢流使栅格通道内维持相似流动,从而使栅格气动性能保持基本平稳;同时由于溢流使流量下降,栅格翼升力效率下降,导弹静稳定性减弱。     

非结构网格限制器研究

研究了6种限制器形式及其参数计算方法,其中新的自适应VenkatakrishnanM限制器能够根据当地流场特征自动调节参数从而实现不同区域不同特性.数值试验表明:非结构网格限制器对邻居集合选择方法较敏感,其中共面法比共点法精度高约15%;高马赫数时压缩型限制器比耗散型限制器计算精度高约10%~40%;参变量算法中,方向梯度法得到的流场比增量法平滑;自适应VenkatakrishnanM限制器能够模拟复杂超声速流动.      

超声速内埋武器不同分离方式分析

为给高空高速无人机内埋武器分离方式的选择提供权衡依据,建立了弹舱及弹体模型,采用chimera嵌套网格方法与Menter SST k-ω湍流模式,对舱内重力投放、舱内弹射投放及舱外重力投放三种内埋武器分离方式在高空高速条件下的内埋武器分离过程进行了仿真分析,对比了三种分离方式分离过程的流场特性与弹体运动参数。结果表明,在Ma=3.7,高度为2.5 km的条件下,舱内弹射投放与舱内重力投放能够使内埋武器安全分离,舱外重力投放在高空高速条件下无法完成内埋武器的安全分离,弹体无法快速下落并伴有大幅振荡,威胁无人机的飞行安全。     

基于三维非结构混合网格的离散伴随优化方法

伴随优化方法在优化过程中不受设计变量个数的限制,有巨大的工程应用前景。基于非结构混合网格技术,采用雷诺平均Navier-Stokes方程、离散伴随优化方法以及自由形面变形（FFD）技术,建立了飞行器气动外形优化设计系统。离散伴随方程采用公式推导的方法直接获取,并运用LU-SGS迭代求解。通过与差分结果进行比较的方式,对通量的雅可比矩阵和目标函数敏感导数开展了确认,并对文中黏性通量项、限制器以及湍流黏性系数的近似处理方法做了分析对比,验证了近似处理方法的合理性和可行性。采用建立的优化系统,完成了ONERA M6机翼在跨声速条件下的减阻优化,并比较了有无容积约束下优化结果的差别。模拟验证结果表明,建立的飞行器气动外形优化设计系统具有较高的可靠性和有效性,可以用于三维飞行器外形的减阻优化。     

几何尺寸对高超声速进气道气动性能的影响

为了探索模型缩尺比对高超声速进气道气动性能的影响,对不同缩尺比的二元高超声速进气道开展了数值模拟研究,结果表明：随着缩尺比的增大,进气道流量系数、隔离段出口总压恢复系数和马赫数均逐渐增大,而静压比逐渐减小,且来流马赫数越高,上述参数变化幅度越大。由理论与数值模拟分析可知,上述现象主要是由于不同缩尺比下,进气道当地雷诺数不同,导致进气道附面层相对厚度变化,进而影响进气道气动性能。理论分析了进气道总压恢复系数与缩尺比的定量关系,就进气道而言,进气道进口处附面层相对厚度减小1%,隔离段出口总压恢复系数提高约0.7%。     

级间通道构型对双脉冲发动机燃烧室局部受热的影响

为研究级间通道构型对双脉冲发动机一脉冲燃烧室气流再附着点的位置和对流换热系数的影响,本文采用基于格心的迎风型有限体积法数值求解雷诺平均Navier-Stokes方程组,空间离散采用AUSM-PW矢通量分裂格式,时间推进采用3阶三步TVD型Runge-Kutta显式方法,湍流模型采用适合模拟分离流动的改进SST湍流模型,并通过经验公式计算再附着点处的对流换热系数。结果表明:亚声速后台阶流动和固体火箭发动机燃烧室内流动的最大绝对误差分别为8.9%和5.8%;二脉冲工作的不同时刻,当二脉冲装药内孔直径大于级间通道直径时,级间孔直径每增大9.1%,气流再附着点位置平均减小28.2%,对流换热系数平均下降9.6%;其它条件相同时,级间孔宽度增加对气流再附着点位置和对流换热的影响比较小;级间孔角度变大使得气流再附着点的位置和对流换热系数分别平均下降3.4%和3.1%。     

旋翼气动力数值模拟准确度影响因素计算

针对旋翼气动力CFD计算时的一些常见影响因素的选择问题,进行了数值计算讨论.计算求解器选择自研代码,该代码采用基于共享内存(OMP)的并行化方法,在嵌套网格系统下求解Navier-Stokes方程,空间方向上使用2阶迎风格式,时间方向采用双时间方法.讨论的问题包括背景网格数目对悬停桨尖涡捕捉、对前飞状态激波捕捉和气动力的影响,远场动量理论修正作用对计算结果影响等.对于讨论的不同网格数目情况下,计算的悬停状态旋翼拉力差异在6%左右,扭矩差异2.6%左右,旋翼前飞阻力和扭矩最大差异为2.58%和3.28%.是否采用动量修正旋翼拉力相差3%~5%.通过计算,给出了针对以上因素计算结果差异的范围,可以为旋翼气动力计算工作提供一定的参考.      

基于耦合传热的双脉冲发动机热防护层受热分析

针对双脉冲软隔层通道孔径对双脉冲发动机I脉冲燃烧室热防护层热环境的影响开展了数值仿真研究。对雷诺平均Navier-Strokes方程采用了双时间步LU-SGS迭代方法、AUSMPW迎风格式以及适合模拟分离流动的改进剪切应力输运(SST)湍流模型进行求解,通过保证流固耦合界面上的热流密度连续来实现耦合传热仿真。采用算例验证了算法及程序的精度和可信度。计算结果表明: I脉冲燃烧室热防护层表面的对流换热系数随软隔层通道孔径的增大而减小,最大对流换热系数平均减幅达20.3%,再附着点位置和对流换热系数最大值所在位置均向上游移动,分别平均移动24.2%和23.1%;I脉冲燃烧室热防护层表面的对流换热系数先增大后减小,且回流区内的对流换热系数相对较小,并在再附着点上游达到最大值。     

某型涡扇发动机高压转静子碰摩故障研究

基于某型双转子涡扇发动机高压转、静子在工作中发生的碰摩现象,通过碰摩消除前后的整机振动响应对比,总结了高压转、静子碰摩的典型振动特征,结合碰摩特点及相关振动理论研究,建立了高压压气机转、静子碰摩模型,应用龙格库塔(Rung-Kutta)法求解模型特定转速下碰摩位置振动响应的频谱图。计算与试验结果表明:双转子结构发动机发生转、静子碰摩时靠近碰摩位置的机匣振动响应会出现次谐波、高次谐波和组合谐波成分,且随碰摩接触面积的增加而增加。     

偏心对双级旋流器出口流场影响试验

为了研究双级旋流器在不同偏心量下的冷态流场变化,设计了4种偏心量(分别为0,0.5,1.0,1.7mm偏心)的双级旋流器试验件,并采用粒子图像测速仪(PIV),测量了不同偏心量下流场.试验结果表明:0.5mm偏心下的流场较无偏心下变化较小,但当偏心量大于0.5mm时,双级旋流器出口下游规则的双涡结构变得不对称,并且中心回流区变窄,流场与无偏心下相比差别较大;偏心对双级旋流器出口下游径向速度分布影响较大,在双级旋流器出口X/D=0.5轴向位置附近,当偏心量大于0.5mm时气流径向速度分量变为正向峰值.