三维烧蚀内部热响应数值计算研究

对热解型碳化复合材料三维烧蚀内部热响应数值计算关键技术进行了研究。采用碳化层—热解面—原始材料层模型,将热解气体与碳化层之间的对流换热处理为源项,通过有限元法建立移动边界条件下温度场求解方程组,采用Gauss Seidel迭代法计算热解气体质量流量和温度场。同时,研究和分析了三维烧蚀移动边界处理方法以及动网格生成方法。由于每个时间步都需要网格重划,烧蚀热防护数值计算对存储效率和计算效率要求较高,本文研究了内部热响应计算中影响存储效率和计算效率的主要因素,并提出了相应的压缩存储方案和求解方案。计算结果表明,移动边界处理方法准确合理;存储方案的存储效率较高;保持刚度矩阵和形函数矩阵正定对称性可以加快温度场计算的收敛速度。     

转捩对压缩拐角激波/边界层干扰分离泡的影响

为了研究转捩对压缩拐角内分离泡结构的影响,进行了来流马赫数2.9,24°压缩拐角激波/转捩边界层干扰的直接数值模拟（DNS）。通过在拐角上游平板的不同流向位置处添加周期性吹吸扰动激发流动转捩,使得转捩不同阶段进入拐角入口,从而在拐角内产生激波/转捩边界层的相互干扰。计算得到的平均速度剖面、壁面压力分布以及分离泡大小与风洞试验及以往直接数值模拟结果吻合较好,验证了计算结果的可靠性。研究了转捩过程对角部干扰区内分离泡结构的影响规律,分析比较了不同转捩阶段下角部分离区内湍动能的生成、耗散和分配机制。研究结果表明：转捩初期的拟序涡结构对分离泡尺度及形状影响最大,发卡涡包在角部拐点附近发生展向融合,并在角部区域形成湍流斑,此时分离泡尺度最小,形状呈现中间高两边低的山峰型。随着转捩的发展,分离区内湍动能生成和近壁区的耗散逐步降低,此时输运项起到了主要的平衡作用。     

基于边界值的多元混沌发动机性能预测算法

首先提取航空发动机排气温度数据序列的边界值,并且证明了边界值序列具有混沌特征.其次提出了一种基于多元相空间重构的发动机状态混沌预测算法,实现对排气温度的预测.通过检验排气温度预测值是否超过所规定的红线,从而进行发动机的健康状态排查.作为验证实例,使用一组某机型发动机实际飞行数据对该预测算法进行了验证,结果表明:该组合算法降低了预测模型的时间复杂度,并大大提高了疑似异常点的预测精度.该方法可以为这种机型发动机故障预测提供决策依据.      

四阶微分方程两点边值问题解的存在唯一性

讨论了非线性四阶常微分方程y（4）=f（x,y,y’,y”,y”’）在混合两点边值条件y’（a）=0,y”（a）＋y”（6）=0,y（b）=0,y”’（b）=0或y’（a）=0,y”＇（a）＋y＇”（b）=0,y（b）=0,y”（6）=0下,解的存在唯一性。其中f在[a,b]×R4上连续且满足Lipschitz条件。并在推广后的Lipschitz条件与Banach压缩映射原理基础上,得到一些新的存在唯一性结果。     

基于升华法实验研究后掠翼三维边界层的转捩

在西北工业大学的低湍流度风洞,采用升华法研究不同雷诺数下后掠翼上表面的转捩现象。实验发现雷诺数较低时,后掠翼上的转捩由流向不稳定触发,转捩位置在最小压力点之后,转捩分界为一条直线;当雷诺数足够大时,转捩由横流驻波不稳定触发,转捩提前到最小压力点之前,转捩分界呈现尖楔形状。该结果表明升华法不但能够较准确地分辨出后掠翼上的转捩位置,还能够区分不同的转捩机理,判断转捩是由流向不稳定还是横流驻波不稳定触发。此外,实验中还发现在横流驻波不稳定增长较大时,升华法能够提供转捩上游区域边界层内的横流不稳定信息;当横流驻波不稳定进一步增强时,该不稳定受模型表面粗糙度的影响较大,萘的喷涂有可能会影响到升华法的结果。     

激光推进辐射击穿温度特性的数值研究

在9种等离子体击穿温度下,数值模拟了二次反射式聚焦系统聚焦情况下激光推力器内流场的演化过程,得到了不同击穿温度对应的能量沉积率、推力峰值、冲量和冲量耦合系数,能量沉积率和推进性能参数在某个等离子体击穿温度值处发生突变。根据空气对激光的逆韧致吸收系数公式,计算了CO2激光辐照下不同等离子体击穿温度对应的空气辐射自由程,发现当等离子体击穿温度为14000K时,辐射自由程为1.4mm,与计算网格的典型尺寸相当,此时入射激光能量在一个网格内以一定效率被吸收,由此确定了基于逆韧致吸收的激光等离子体的击穿温度。     

弹性整流罩分离的流固耦合仿真方法

发展了一种流固耦合的动力学数值仿真方法,用以模拟弹性结构在其运动诱导的非定常气动载荷作用下的动态响应.针对给定的运动规律,对不可压流场采用侵入式边界方法求解了非定常气动载荷,其中未计及弹性结构变形对流场的影响.对非定常气动载荷作用下的薄壁弹性结构进行了瞬态动力学数值模拟.以地面试验状态的整流罩分离为例,对刚性体及是否考虑表面气体压力作用的柔性体等情形下的整流罩的动力学行为进行了对比分析,验证了发展的流固耦合动力学仿真方法的可行性.     

基于翼面压力的飞行器气动力感知技术与自由飞验证

飞行器在大迎角飞行状态下其复杂绕流流动会导致非指令运动出现,严重影响了飞行器的操纵性与飞行安全。现有以惯性元件为核心的机载设备无法直接提供非定常气动力参数,而如何实时感知飞行器大迎角状态的非定常气动力/力矩,是抑制非指令运动现象的核心所在,将是未来战机设计中亟待解决的空气动力学和飞行控制的关键问题之一。针对上述问题,提出基于翼面压力信息获取特征截面滚转力矩系数C_lsec,估算飞行器全机在大迎角状态下的非定常气动力矩,进而判断飞行器的滚转运动的设想。风洞和飞行试验研究结果表明：对于80°/48°双三角翼,0.8c特征截面滚转力矩系数C_lsec与模型滚转力矩存在相关性;在飞行器进行大迎角平飞动作时,非指令滚转运动下的滚转力矩系数C_lsec大幅增加;C_lsec能够比惯性传感器提前预测模型的滚转运动趋势,可为飞行器大迎角状态的非指令运动抑制提供一定数据依据。     

预处理特征型边界条件及其回流处理技术

通过特征型预处理Euler方程推导出预处理特征型的边界条件处理方法,并发展了预处理的回流处理技术,增强计算的鲁棒性和精度。通过特征型预处理Euler方程的特征关系得到边界条件处理方法;在亚声速速入口和出口边界,根据其特征方程确定回流边界处理方法。通过数值模拟发现:带回流的特征型边界处理方法比Turkel方法收敛快1.5倍以上,而计算精度与SIMPLE方法相当;在出口带有回流的计算中,其压力误差比Fluent降低68.3%。由此可知带有同流处理的特征型边界处理方法具有较高的鲁棒性和精度,在实际的工程计算中有一定的应用价值。