旋成体射弹倾斜入水运动仿真

为了研究运动参数和弹头外形对弹体斜入水过程的影响规律,采用气液两相流体积分数和水汽空化模型,通过嵌套网格实现刚体三自由度运动学和动力学耦合,模拟了弹体以80～100 m/s速度倾斜入水开空泡阶段的运动过程。经文献实验验证,入水弹体速度与位移的误差为0～6%和-8%～0,转动角度误差为-6%～0。通过对入水速度和入水角度的多工况模拟研究,发现入水速度增大,弹体轴向冲击载荷增大,最大载荷与速度的平方呈线性关系,弹体速度非线性衰减率大;入水角增大,弹体转动角速率减小,运动稳定性强,速度衰减率不受入水角影响。与圆锥头部弹体相比,采用头部阶梯状修型后的弹体的平均速度衰减率、转动角速率和最大轴向冲击载荷分别降低到66.7%、40%和77.2%,显著提高了运动稳定性。     

可压缩流大涡模拟耦合概率密度函数方法中的滤波压力模型

亚格子组分-温度关联项显著影响反应流大涡模拟精度。利用概率密度函数方法中概率等效的特点,发展一种新的滤波压力模型,可以良好封闭亚格子组分-温度关联项。介绍概率密度函数及其耦合求解方法,在已有模型基础上推导建立新的滤波压力模型,并在三维超声速氢气/空气时间发展反应混合层中对不同的滤波压力模型进行数值测试。结果表明,与传统的滤波压力模型相比,新的滤波压力模型可以明显改善反应混合层的模拟准确度。特别地,基于新的滤波压力模型,大涡模拟耦合概率密度函数方法可以较好地模拟链式反应中间微小组分如超氧化氢基等,有望更有效地再现自点火等复杂现象。     

Spalart-Allmaras湍流模型在高超声速气动加热计算中的应用

将AUSM+格式与LU SGS隐式迭代相结合,采用一方程的Spalart Allmaras湍流模型和二阶迎风MUSCL格式,研究了火星"探路者号"实验模型在10马赫来流条件下的气动热问题。采用量热完全和热完全两种气体模型的计算结果与实验数据进行了比较。热完全气体模型模拟的热流分布规律与实验吻合较好,但是对分离区的热流估计过高。在网格合理划分,尤其是满足近壁面的网格分辨率要求的前提下,该方法能够合理反映高超声速流场除分离区以外的气动热特性。     

姿控发动机脉冲工作的管路瞬变特性

采用有限差分格式的特征线方法 ,对某姿控推进系统发动机真实推进剂管路在脉冲工作情况下的瞬变特性进行了数值研究 ,对影响发动机多次脉冲启动的重要参数进行了定量计算 ,对比分析了试验结果与计算结果的一致性。说明该方法可方便地对包括许多管路、阀门和推力室的复杂系统进行动态性能分析 ,所得结果对于液体火箭发动机管路工程设计具有重要参考价值     

激波诱导惰性堆积粉尘运动的建模及数值模拟

基于空间平均的双流体模型 ,引入描述颗粒脉动速度的粒化温度 ,建立颗粒相的本构方程 ,将粒化温度模型推广到解决高速流动问题。采用AUSM+ 有限差分法 ,数值模拟激波在一定厚度的惰性粉尘床中传播及诱导粉尘颗粒运动的过程。结果表明 ,粒化温度模型较好地描述了这一过程 ,证实了此模型在处理稠密颗粒高速流动和流场存在强间断中的可行性     

超声速、高超声速球锥组合体表面压强估算方法及其应用

改进并推广了一种有效的球头单锥组合体在超声速和高超声速流动情况下的表面压力分布算法 ,并利用该算法对若干来流条件下球头双锥组合体表面的压力分布进行了计算。计算结果与数值计算结果以及文献实验数据的比较表明 :本文介绍的方法不失为一种适用范围较宽、准确程度较高的工程算法     

球锥体高超声速绕流的俯仰阻尼导数的数值计算

采用Navier Stokes方程描述球锥体外形非定常振动流场 ,在Etkin理论下给出俯仰阻尼导数的计算公式。对定常流场的计算采用ADI形式的NND格式 ,对非定常流场的计算采用四步Runge Kutta方法 ,并引入变系数残值光顺技术加速收敛。将俯仰阻尼导数计算结果与实验及工程计算进行了比较 ,并数值研究了其随振动中心的变化规律。     

飞船再入舱三维化学非平衡流数值模拟

对飞船再入舱高超声速化学非平衡三维流场进行数值模拟。控制方程为含化学反应源项的全Ｎａｖｉｅｒ－ Ｓｔｏｋｅｓ 方程组, 化学模型为高温空气的七组元模型, 其组元成分为Ｎ２ 、Ｏ２、ＮＯ、Ｎ、Ｏ、ＮＯ＋ 、ｅ－ 。差分格式采用张涵信院士提出的ＮＮＤ激波捕捉格式并用全耦合方法处理流动方程及化学反应方程, 对化学反应源项采用全隐式处理, 运用时间预处理技术以加快收敛。同时, 也对完全气体模型进行计算, 以分析真实气体效应。与实验数据对比, 说明本文结果是可信的。     

高超声速伸缩翼变形飞行器轨迹多目标优化

针对高超声速条件下变形技术的应用模式，对具有伸缩翼的组合式飞行器滑翔弹道进行了多目标优化研究。介绍了伸缩翼的变形模式，给出了不同变形状态下的气动特性；建立了三自由度滑翔轨迹动力学模型和伸缩翼前缘热流计算模型；采用MOEA/D多目标优化算法，以变形条件和飞行攻角为设计变量、以最大射程和最小翼前缘总吸热量为目标函数，进行了多目标优化计算。优化结果表明，MOEA/D计算得到了相对均匀分布的Pareto最优解集，将伸缩翼外形与无变形外形相比，飞行器滑翔段射程得到了显著提高，同时伸缩翼前缘总吸热量有明显的降低。     

反导作战信息流时延模型研究

在现代防空反导作战中,战术弹道导弹的威胁,给防空系统带来了预警时间不足等一些新的问题。为了提高防空系统的反导预警能力,在对反导作战过程和反导预警作战过程进行分析的基础之上,对反导预警信息流程进行了研究。提出了一种边处理-边分发的实时模式,通过对比分析,所提模型缩短了预警信息的传输时间,使作战部队提前预警。