航天器发动机羽流对敏感器热效应仿真研究

采用差分求解N-S方程与DSMC方法相结合的方法,研究了航天器单台发动机连续工作情况下真空羽流对航天器敏感器的热效应。首先通过求解N-S方程,获得发动机喷管的内流场,然后应用DSMC方法对喷管出口外轴对称羽流场进行计算,最后将轴对称羽流场计算结果作为模拟粒子入口边界条件,在并行计算机平台上进行三维羽流场和热效应计算得到航天器单发动机连续工作情况下羽流场对敏感器的热效应。以两个敏感器为例,对仿真结果进行了分析和比较,并得出了相应结论。     

卫星发动机羽流污染模型的DSMC方法研究

文章通过DSMC（Direct Simulation Monte Carlo,直接模拟蒙特卡罗）方法对卫星发动机羽流进行数值模拟,对仿真方法进行了对比验证,证实了DSMC方法可以比较精确地模拟羽流场。在建立卫星污染模型的过程中,针对污染发动机喷口形状的变化以及出口参数的变化,对DSMC方法模拟结果进行了对比,得到喷口形状和出口参数的变化对模拟结果的影响,并得到最优的卫星敏感表面羽流污染计算模型。     

航天器姿控发动机真空羽流场计算及其扰动分析

航天器姿态控制发动机工作期间,产生的羽流扰动力矩会影响控制精度,对流和辐射产生的热效应会影响热控系统工作,羽流沉积物会影响光学敏感器的精度。文章以某航天器的20 N姿控发动机为例,首先采用工程的MOC方法计算其真空羽流场;然后利用建立的航天器羽流三维冲击模型分析了羽流对太阳电池阵的冲击载荷;最后对发动机羽流脉冲激励下太阳电池阵的响应和姿态控制系统受到的扰动力矩进行了仿真分析。结果表明姿控发动机羽流脉冲激励将对控制精度和稳定度产生不可忽视的影响。     

卫星变轨发动机羽流污染的研究

针对卫星变轨发动机工作时高空羽流污染问题 ,从分子运动论出发 ,采用直接模拟MonteCarlo(DSMC)方法对轴对称羽流场进行模拟 ,发动机喷管出口参数采用内流计算的给定值。首先对计算方法进行了验证 ,然后通过计算给出羽流场压力、密度等参数 ,给出计算位于回流区内的卫星表面敏感部件的羽流污染量的方法 ,并进行了计算。     

空间发动机羽流研究技术发展综述

首先介绍了目前进行空间发动机羽流研究的必要性,同时说明地面试验和数值模拟方法都是研究空间发动机羽流特性的有效手段,两者缺一不可。在此基础上,总结了国内外羽流地面试验关键技术和发展状况。然后,分别总结了国内外最具代表性的空间发动机羽流试验台的组成、真空抽吸方式、主要技术指标和特点,包括美国的J2-A试验舱和CHAFF-IV试验舱,欧洲的CCG羽流污染试验舱和STG低温氦冷羽流试验舱,中国的KM系列空间环境模拟器和PES地面羽流试验台。最后,介绍了与羽流地面试验相关的数值模拟技术的发展,总结了进行羽流数值模拟的模型,重点介绍了常用的DSMC方法的典型应用和基于此方法所开发软件的情况,并针对大密度羽流场和电推进发动机羽流场的特点分别总结了其进行羽流场计算的方法。     

空间CO2羽流压力场的数值计算和实验研究

要解决空间羽流污染问题,首先必须精确掌握羽流场的特征。为此,在真空深冷背景环境下进行了CO2气体模拟发动机羽流压力场的数值计算和实验研究。采用经典的计算流体动力学(CFD)和直接模拟Monte Carlo(DSMC)相结合的方法进行羽流场的数值模拟。试验模拟空间高度在100公里以上,在整个测试过程真空舱内压力维持在10-3Pａ量级,冷阱温度为93±5Ｋ。计算结果与测试数据吻合较好,表明文中计算方法和物理模型可以准确描述真空深冷空间环境下的羽流场特征。     

月球探测器着陆过程羽流热效应数值模拟研究

考虑到月球探测器着陆过程变推力发动机工作可能引发羽流热效应,通过差分求解N S方程与直接模拟蒙特卡洛(DSMC)耦合的方法,针对探测器变推力发动机羽流的连续流区与稀薄流区进行了数值模拟,获得了着陆过程不同状态探测器表面的羽流气动热流密度分布,分析了探测器在不同高度及坡度受到的羽流热效应,研究了探测器因落月惯性与发动机延时拖尾等特殊工况羽流气动热的影响程度。数值模拟结果表明,探测器受到的羽流气动热影响总体随着发动机出口与月面间距离减小而急剧增强,随着距离减小至0.434 m,热流密度最大值为429 kW/m 2 ,而且月面坡度对羽流热效应有减弱作用。研究结果可为探测器关机策略的制定提供支撑,为探测器关键部位的热防护提供输入,为探测器整体的设计优化提供服务。     

卫星回收舱再入过渡流区气动热数值计算

返回式卫星进入地球大气层95km高度左右时,速度达到Ma20,回收舱被高超声速稀薄来流形成的弓形激波环绕,气动加热问题非常明显,准确预测过渡流区的气动热成为一个十分突出的问题。过渡区由于气体空气分子仍然比较密集,基于分子动力学的直接蒙特卡洛模拟方法(DSMC)极其耗费计算资源,而求解Navier-Stokes(N-S)方程的方法误差较大。采用添加二阶滑移条件的N-S方程求解过渡流区气动热,并与开源DSMC2V程序计算结果对比,研究了舱体母线热流、压强系数变化;通过对比两种方法下舱体前缘弓形激波及流场参数变化,分析了滑移条件影响壁面和流场参数的机理。结果表明:回收舱再入过渡区时,钝头驻点区为高压、高热流区,锥身区气动热和壁面压强保持在较低水平。带滑移条件的计算机流体力学方法(ComputationalFluidDynamics,CFD)计算得到的热流、压强系数与DSMC结果吻合良好,具有计算效率高、精度较高的优势。对流场压力、温度、速度等参数分析显示,滑移条件中壁面速度滑移和温度跳跃的加入,改变了壁面流动参数,进而改善了壁面热流和压强的准确模拟能力。滑移条件对外流场参数的影响极小,没有改善稀薄流区流场的模拟能力和激波捕获能力,模拟得到的温度、速度等流场参数及激波位置、激波层厚度等与DSMC结果仍有差别。可以认为,文章采用方法能够满足工程上快速高效预测卫星回收舱再入过渡流区气动热需要。     

微波等离子推力器真空羽流模拟计算

粒子节点－直接模拟蒙特卡罗法（DSMC）是求解稀薄气体流动的数值模拟方法，采用轴对称硬球模型、随机取样频率法对微波等离子推力器（MPT）的真空羽流进行了数值模拟计算，为了提高计算效率和节省计算机内存，采用了随机取样频率法（RSF）和加权技术。对MPT真空羽流的浓度、速度和温度分布进行了计算，并对羽流污染进行了分析。计算结果对以后MPT真空羽流场测试和航天器一体化设计提供了参考。     

脉冲等离子体推力器羽流回流影响分析

将三维MHD双温入口模型的计算结果作为入口条件,运用DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法,模拟实验室PPT样机羽流。验证计算显示,该模型具有模拟脉冲等离子体推力器羽流的能力。对不同初始电压和电容下的羽流场进行了模拟,给出了出口平面返流质量流率的变化情况和20μs羽流中离子和CEX离子的分布情况。计算结果显示,高能量状态对应高质量流率,高质量流率对应高动量,高动量离子和中性粒子对航天器撞击会造成更强影响。     

羽流污染对某型号卫星热控涂层的影响分析

文章以某卫星热控涂层的污染研究为目的,针对20N轨控推力器,在推力器羽流流场模拟的基础上,结合污染沉积过程模拟,计算获得了卫星在轨运行期间热控涂层表面的污染沉积膜的厚度,分析了热控涂层太阳吸收率的退化情况。     

一种卫星双组元推力器羽流热影响分析方法

给出一种分析星上设备受双组元推力器羽流热影响的方法。首先,提取待分析设备三维模型的STL(stereo lithographic)数据;然后,根据STL数据中三角形网格的法向矢量选取有效的三角形网格;最后,利用STL数据和推力器羽流热流场分布函数建立羽流热影响分析模型。不管设备表面能否用解析公式表达,文章提出的方法均能分析推力器羽流对其产生的热影响。     

飞行环境对火箭喷焰微波衰减特性的影响

本文采用一种近场和远场分区耦合、时间相关法和空间推进法相结合的方法计算分析了飞行环境对火箭喷焰微波衰减特性的影响。计算结果表明高空飞行时的喷焰温度低于地面环境下的喷焰温度。但是,由于市容工环境压力的降低使得喷焰的复合效应减弱导致自由电子数密度沿喷焰轴向衰减减慢,使得高空飞行时喷焰对微波的有效作用范围增大。由此可以推断实际飞行时喷焰的微波衰减比在地面环境时更为严重。     

PWS软件应用于探月着陆器羽流效应 数值模拟研究

文章使用自行开发的PWS计算软件就探月着陆器主发动机羽流对着陆缓冲机构产生的气动力和对流换热热效应问题进行了数值模拟研究。PWS软件是一个基于直接模拟蒙特卡罗方法的羽流计算通用软件,其松散的软件架构保证了各个计算模块的相对独立性,使用“与”、“或”逻辑法则的二级构造法来实现边界条件的通用性。软件的计算值与国外CUBRC实验测量值符合得很好。使用PWS软件对着陆器主发动机羽流进行数值模拟计算表明着陆缓冲机构所受的羽流气动压强最大达到2.4 Pa;缓冲机构底盘侧面的对流换热热流密度最大达到2 000 W/m2;当缓冲机构壁面温度从300 K增加到500 K时,其表面对流换热流密度值有所下降,其最大下降值不到5%。     

液体火箭燃烧尾焰冲击干扰效应数值研究

火箭发射时其燃烧尾焰的冲击干扰效应对发射稳定性和发射架、导流槽等地面设施有重要影响。文章采用压力隐式算子分裂算法,通过求解Navier-Stokes方程,对氢氧液体火箭发动机燃烧室内燃烧过程与尾焰流场进行了一体化数值计算,得到了火箭发射后尾焰与地面撞击产生的冲击流场。结果表明:尾焰流场计算模型、方法与结果合理;尾焰冲击干扰效应会大幅度提高地面附近的压力和温度;火箭尾焰撞击地面后,高温区出现在离地面一定距离的高温层内,此时地面附近为低速区;尾焰对其正下面的地面区域产生冲击最大,主要干扰区域集中于半径为15 m的圆形区域。     

万有引力场中带挠性太阳帆板航天器的姿态稳定性

本文讨论带双侧挠性太阳帆板航天器在万有引力场中的姿态运动，建立带挠性帆板航天器的欧拉方程和帆板强迫振动方程。利用Ｇａｌｅｒｋｉｎ方法对动力学方程离散化，利用Ｋｅｌｖｉｎ－Ｔａｉｔ－Ｃｈｅｔａｙｅｖ定量判断航天器在轨道坐标系内相对平衡的稳定性。导出适用于任意阶模态的解析形式稳定性充分条件。     

燃烧室内燃烧模型对尾焰流场及其辐射的影响

基于压力隐式算子分裂(PISO)算法,通过求解Navier-Stokes方程,对燃烧室内一步反应和两步反应模型、无燃烧室三种情形下尾焰流场进行了数值仿真;采用结构和非结构网格并分别用TTM方法和Delaunay三角形方法来生成;利用高温气体高分辨率光谱参数数据库HITEMP对辐射传输方程进行求解,得到三种情形下尾焰中CO2和H2O的光谱辐射亮度分布。仿真结果表明:不同的燃烧模型影响尾焰流场及其辐射,无燃烧室时尾焰辐射较弱。