小推力发动机高空羽流场数值模拟

以小推力发动机的高空羽流场为研究对象，完成了氮气流的DSMC方法数值模拟研究，对计算的可靠性进行了实验对比验证，分析了高空羽流场特性及高空稀薄流动的非平衡效应。结果表明，用DSMC方法与加密网格技术结合可有效模拟高空羽流场，且必须计及气体非平衡效应。皮托压力的数值结果与实验符合得很好。     

基于简化Boltzmann方程的超声速流稀薄效应分析

从近连续流到近自由分子流,计算分析超声速流场的稀薄效应。模型Boltzmann方程先后经简化分布函数和离散速度坐标变换后,采用一个隐式通量修正二阶迎风TVD格式差分求解。壁面取漫反射气体分子模型。在Knudsen数为0.001,0.01,0.1,1.0,10等5种情形下,数值模拟二维圆柱氩气体绕流,观察到了不同的波系及尾迹结构。总阻力系数计算值同实验数据基本吻合,碰撞项Shakov动力学模型反映的流场细节略好于修正BGK模型。稀薄效应趋于弱化激波等超声速流动结构,物体影响域增大。     

发动机内剩余推进剂排放对卫星表面污染的研究

应用真空物理的相关理论，分析了在高空超高真空环境中，发动机二次启动前推进剂排放及其在空中迁移的物理过程，给出模型假设和模拟方程，并针对此真空环境及物理过程应用合适的计算方法地方程进行求解。最后计算了推进剂的排放对卫星造成的污染，给出计算域内浓度变化的非定常过程及卫星表面不同时刻的最大污染量，为推进剂排放方案的实施与否提供理论依据。     

返回器稀薄区气动特性工程计算方法的应用研究

中国航天返回器再入过程采用非弹道式再入轨道,其高空稀薄区滞留时间显著增长,稀薄效应对再入飞行的影响显著增强,因此快速而准确地预测返回器稀薄区气动特性变得非常重要。文章分析归纳了多种稀薄区气动特性工程计算方法即桥函数方法的发展和应用,并以 Stardust 返回舱为对象,对比分析了三种不同桥函数的预测精度,给出了在航天返回器气动预测中更为合适的工程方法。结果显示：不同桥函数预测结果差别很大;在不同攻角下,与数值模拟的对比分析表明,局部桥函数方法气动特性预测结果基本优于其它桥函数,尤其是在力矩系数预测上。因此稀薄区气动特性的预测采用局部桥函数较为合适。     

高超声速钝球柱外形表面热流分布研究

文章采用数值求解雷诺平均Navier-Stokes方程方法对高超声速钝球柱地面试验模型进行了计算,比较了不同头部半径和肩部半径条件下表面热流的分布特征。研究表明：钝球柱头部区域表面热流分布与头部半径和肩部半径紧密相关。其中,头部半径对热流分布特征和热流大小均存在显著影响。当头部半径较小时,头部区域热流呈现双峰值分布,热流由驻点峰值逐渐下降,接近肩部时开始上升,至肩部峰值点后又继续下降。随着头部半径的增大,驻点热流减小,肩部热流增加,双峰值分布逐渐演变为单峰值分布,热流由驻点单调上升至肩部峰值点;相对于头部半径,肩部半径主要影响热流大小,驻点热流和肩部热流均随肩部半径的增大而增加,不同肩部半径计算模型的头部热流分布规律基本一致。     

碳/碳材料细观尺度的氧扩散特性与烧蚀分析

将DSMC方法运用在材料细观烧蚀机理分析中,对碳/碳材料在典型烧蚀过程中的简化模型进行了微米量级下的氧扩散特性分析,旨在研究氧分子在材料表面缝隙中的扩散特性及其与壁面的作用规律,计算分析结果验证了纤维烧蚀结构演变特性。结果表明：沿缝隙深度方向,O2与壁面的碰撞频率不断降低;在缝隙入口处烧蚀最快,缝隙深处最慢,碳纤维随着烧蚀过程的进行不断尖化直至达到强度极限产生折断剥离;随着烧蚀过程的不断进行,O2与壁面的碰撞频率增加,材料将烧蚀得更快;在同一烧蚀条件下,O2与壁面发生烧蚀反应的概率越大,将消耗O2越快,致使O2与壁面的碰撞频率降低。     

超低轨航天器气动特性快速预测的试验粒子Monte Carlo方法

超低轨(LEO)卫星气动特性的快速准确计算是对其进行轨道预测和控制的关键输入条件。基于真空技术领域中计算管道分子流率的试验粒子Monte Carlo(TPMC)方法,结合自由分子流理论,发展了一套快速准确预测低轨卫星气动特性的TPMC方法,给出了其模拟步骤及主要关键技术点,并采用该方法模拟了带电池翼超低轨卫星的气动特性和航天器典型构件之间的多次反射效应。结果表明:TPMC方法在计算超低轨航天器气动力、力矩时具有较高的可靠性和对工程复杂外形的适用性;该方法能够准确模拟自由分子流理论无法求解的多次反射问题,给出正确的气动力系数;该方法的计算速度比直接模拟Monte Carlo(DSMC)快3~4个量级,存储量要求也比后者低1~2个量级,是超低轨航天器气动特性快速准确预测的一个理想方法。     

自由分子点源法在非均匀出口羽流场模拟中的应用

为满足高空稀薄流工程估算的精度和效率需求，将Woronovicz自由分子点源模型推广至非均匀出口发动机羽流场的计算模拟。采用分格技术剖分得到多独立点源，以此解决出口平面参数的非均匀性问题。10N双组元发动机羽流的质量和压力计算结果同实验符合良好，验证了自由分子点源模型的有效性。某实验型喷管羽流计算表明，点源模型在远场与DSMC方法基本吻合．两者误差在引入非均匀出口条件修正后得到显著改善。     

火星进入器高空稀薄气动特性

针对火星稀薄大气环境进入器气动特性问题,以类火星科学实验室外形为例,计算分析火星稀薄大气真实气体效应对气动特性的影响,给出火星高空稀薄环境下的气动特性规律。研究发现,随着飞行高度的增加,稀薄度增加,激波脱体距离、激波厚度增大,激波强度减弱,明显的激波结构逐渐消失,流场等值线更趋于圆弧状分布;真实气体效应使得迎风面压缩及背风面膨胀增强,轴向力、法向力及顶点力矩系数等预测结果与完全气体模型预测结果相比绝对值偏大;随着稀薄度增大,轴向力、法向力及顶点力矩系数等绝对值增大,在同样的迎角下,随稀薄度的增加,纵向压心前移,进入器的静稳定性变差。     

地球大气再入返回器辐射加热计算方法研究

针对地球大气再入返回器的辐射加热计算方法开展了研究,采用NASA的高温空气辐射加热预测软件QRAD中用到的四光谱带模型,同时考虑了辐射散热效应及非平衡辐射加热.对于返回舱大底的辐射加热分布特性,采用了基于牛顿压力分布的计算方法.计算结果表明,发展的高超声速再人体辐射加热预测方法结果可靠.