高超声速小钝锥边界层非对称转捩研究

在脉动能-比耗散率(ｋ-ω)两方程体系上引入边界层转捩的模型化处理方法,构建高超音速边界层流动的转捩湍流模式。基于新的转捩湍流预测手段研究发现,圆锥小攻角下转捩特性呈现非对称分布———迎风面后移而背风面前移,非对称转捩对于高超音速飞行器的表面热流特性及气动力特性影响显著。     

带减阻杆高超声速飞行器外形气动特性研究

采用高超声速风洞测力试验方法测量钝头飞行器头部减阻杆的高超声速气动特性,研究减阻杆的气动减阻原理,分析了多组不同构型减阻杆的减阻效果。结果表明,减阻杆显著减少了钝头飞行器高超声速的阻力,最大的减阻率达到60％之多;减阻效果与减阻杆构型和迎角状态密切相关;减阻杆会诱发稳定性、“热斑”以及非定常脉动等不利问题。     

飞行器气动/飞行/控制一体化耦合模拟技术

基于非结构重叠网格技术和全局亚迭方法,发展了一种气动、飞行和控制一体化耦合的非定常求解策略,用于飞行器闭环控制的数值虚拟飞行研究。高超声速飞行器按照设计控制律,采用升降舵偏转进行连续变攻角操纵和静不稳定布局稳定控制;数值虚拟飞行逼真显现了上述非定常流动、运动与控制的时序演化过程。数值虚拟飞行与基于气动数据库插值的飞行仿真比较,控制过程和反馈效果十分吻合。     

基于RANS-LES混合方法的翼型大迎角非定常分离流动研究

使用雷诺平均Navier-Stokes方程-大涡模拟（RANS-LES）混合方法中的延迟分离涡模拟（DDES）方法,模拟了NACA 0015翼型在大迎角下的静态绕流和强迫振荡运动并与实验值进行了比较。在大迎角静态翼型大分离流动模拟中,DDES方法捕获了非定常RANS计算未能获得的机翼背风面的涡脱落现象。在所采用的RANS和DDES模型中,基于剪切应力输运（SST）湍流模型的SST-DDES混合方法给出的时均压力系数分布与实验吻合得最好。在大迎角强迫振荡翼型绕流模拟中,DDES方法得到的非定常气动载荷与实验值吻合得很好,正确地反映了最大迎角处阻力和俯仰力矩的阶跃性突变,而非定常RANS计算则给出了完全错误的趋势。     

弯曲和旋转修正的湍流模型在旋转翼身组合弹箭中的应用研究

为了研究高速旋转产生的流线弯曲及壁面强湍流剪切效应,本文采用完全时间相关的非定常N-S方程,对超声速带翼旋转弹箭开展计算,研究了弯曲和旋转修正的湍流模型SARC和SSTRC对弹箭旋转气动特性和流场结构产生的影响。结果表明:对全弹侧向动态特性计算,弯曲和旋转修正的湍流模型与原始模型精度相当,侧向力和力矩旋转导数最大差异6%,4组经验估算公式计算的马格努斯力旋转导数与本文结果误差皆15%。弯曲和旋转修正的湍流模型使物面压力左右两侧同时偏大或偏小,与原始模型相比并没有加剧或削弱不对称效应,这是全弹马格努斯特性变化不大的原因。弯曲和旋转修正湍流模型预测的分离区更大,对分离流动的抑制能力减弱。     

飞行器自激振荡的流动物理与动力学机制

全局亚迭代耦合求解流体动力学方程和刚体动力学方程,研究轴对称飞行器单自由度俯仰失稳运动的非定常特征.高超音速锥-柱-裙飞行器的单自由度俯仰失稳运动发展为极限环形式,周期性运动伴随着波系结构的非定常变化;自激振荡源于静稳定构型的振动回复机制和尾裙激波流动迟滞效应所构成的阻尼机制.基于第二拉格朗日方程和虚功原理,导出能够描述迟滞现象的参数化非线性动力学模型.多尺度近似分析获得参数化运动特征:自激振动过程是拟简谐运动;静稳定构型才能出现自激振荡周期运动;平衡点阻尼是决定运动稳定特性的分叉参数;振幅特性与阻尼非线性相关,频率特性与刚度非线性相关.基于数值结果的参数辨识,非线性模型的理论分析、重构都与数值结果高度一致,从而有效地佐证了自激振荡建模研究的合理性.     

细长三角翼滚转/侧滑耦合运动的数值研究

 采用亚迭代耦合求解流体动力学方程和刚体动力学方程数值研究方法,建立无时间滞后的耦合数值手段来分析研究飞行器非定常的耦合运动特征,研究细长三角翼多自由度非定常的耦合关联运动。研究表明：在流动结构失稳后机体运动逐步形成极限环振荡的自维持运动,滚转力矩系数滞回曲线呈现典型的“双8”稳定形态;在滚转和侧滑两自由度下,三角翼进行“落叶式飘荡”耦合运动,自激滚转振荡更为剧烈。     

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

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

高超声速流动转捩的数值研究

基于线性稳定性的结果，以类似于处理湍流脉动的方式来模化流场脉动，用以研究高超声速的流动转捩，包括转捩起始位置及转捩区长度．研究中采用Faver平均Navier—Stokes方程k-ε湍流模型，分别以Roe二阶流通昔差分分裂格式和Chakravarthy三阶迎风TVD格式进行离散，数值模拟了椭圆锥的高超声速绕流流场。分析比较了不同的湍流转捩机制，并研究了雷诺数的影响规律。结果表明：Bypass转捩机制是传统的高湍流度风洞中模型流动转捩的主控因素。