二级火箭喷流对底部热环境影响的数值模拟

为研究二级火箭工作时对底部热环境的影响,在不同喷管固壁热辐射和飞行高度条件下,用Fluent软件自发动机喷管喉部处计算喷流流场和底部区域内给定平面的对流热流密度,并建立了三维反向蒙特卡罗（RMC）法模型,同时根据HITRAN数据库计算气体红外辐射特性,获得喷流对给定平面的辐射热流密度,两者相加即为喷流到达研究平面的总热流密度。结果表明：热辐射和热对流在底部产生的热效应趋势不同。     

计算火箭发动机尾喷管红外特性的封闭腔-净辐射法

为计算火箭发动机尾喷管的红外特性,建立了封闭腔-净辐射法,推导了光谱透射因子和光谱吸收因子,以及封闭腔表面辐射净热流-表面及气体温度的方程。作为验证,计算了一个充满均匀辐射性介质的轴对称柱形封闭腔的辐射换热,与文献的计算结果比较有很好的一致性;还计算了一火箭发动机尾喷管的辐射特性,分析了该喷管在中波红外波带2~6μm和长波红外波带8~14μm及全波长的有效辐射和净辐射热流。结果表明,在气体辐射作用下,火箭发动机喷管进出口平面的辐射有一定的光谱差异性;壁面有效辐射沿轴向的分布取决于壁温变化,平直段和收敛段有效辐射较强,扩张段有效辐射较弱;中波红外波带2~6μm喷管的辐射热流占全波长辐射能量的60%以上。     

低空多喷管发动机喷焰红外特性研究

对多喷管火箭发动机低空尾焰的红外辐射特性进行了研究.基于热流法建立了喷焰红外传输计算模型,数值模拟了四喷管发动机尾流场2～5 μm光谱的红外特性.获得了光谱、波段辐射强度仿真数据及红外仿真热像,并与单喷管喷焰的红外特性进行了比较.考察了不同的喷管间距、探测方向、飞行高度对喷焰表观辐射强度的影响.结果表明,多喷管尾喷焰红外辐射光谱选择性与单喷管相似,但光谱峰值无线性关系;对于低空四喷管发动机尾喷焰,红外特性随喷管间距增加有所增强,随飞行高度增加,喷焰红外辐射强度提高.     

固体火箭发动机喷管传热与壁面烧蚀的耦合计算分析

为研究某型固体发动机在地面工作过程中喷管的受热与烧蚀,对其工作后140 s内复合喷管壁面受到管内高温喷流辐射与对流加热,以及发动机外部环境辐射与对流冷却条件下的壁面受热与材料热解烧蚀建立一维非稳态热分析模型进行计算分析。其中,喷管材料采用金属基体内衬高硅氧-酚醛复合隔热材料构成,高温喷流对喷管的辐射加热采用非灰参与性介质的封闭腔辐射换热模型计算,对喷管的对流加热采用巴兹公式计算,复合喷管壁面材料升温后的热解分为基体材料升温-基体材料热解-热解层炭化-Si O2熔融-炭化层脱落五个阶段进行分析。研究发现,喷管收敛段和喉部主要受到高温喷流的辐射加热,内壁辐射热流约为对流热流的2.5倍,喉部下游因喷流温度下降,速度激增,内壁对流热流超过辐射热流,在扩张段尾部,内壁的辐射热流再次超过对流热流;发动机工作过程中,喷管收敛段和喉部壁面的高硅氧-酚醛复合隔热材料随时间逐渐被烧蚀,烧蚀厚度随时间上升,喉部烧蚀厚度最大,140 s时烧蚀厚度达到8 mm,平均烧蚀速率为0.057 mm/s;喷管扩张段中后段喷流温度大幅下降,壁面内高硅氧-酚醛复合隔热材料未烧蚀;沿喷管壁面厚度自内向外,壁面温度急剧下降,发动机工作后16 s时,喉部截面处内壁温度达到2700 K,而外壁温度仅为340 K。     

用FVM法计算固体火箭羽流的红外特性

建立了有限体积法(FVM)的理论模型模拟固体火箭羽流的红外特性,研究了燃气组分H2O,CO2,CO,HCL,OH,NO的吸收和发射,以及AL2O3粒子云的吸收、发射和散射,还研究了两相流场中气相和粒子云温度不均匀时对辐射传输方程的修正。计算一单喷管固体火箭的二维轴对称羽流在光谱2~5μm的红外特性,给出羽流辐射的空间分布和光谱分布并与文献计算的结果对比。研究表明:有限体积法能较好地模拟固体火箭羽流的红外特性;粒子辐射占固体火箭羽流红外辐射的主要部分;在光谱2~4μm粒子辐射使羽流辐射增强,而在光谱4~5μm粒子辐射使羽流辐射减弱。     

卫星系统热特性分析

考虑空间轨道外热流、卫星表面自身辐射、热载荷等因素影响,建立卫星温度场计算模型,在采用蒙特卡罗(Monte-Carlo)法求解卫星复杂辐射边界条件的基础上,利用有限容积法对卫星在轨飞行阶段的瞬态温度场进行数值模拟,计算得到卫星瞬态温度场,并考虑其表面自身辐射及空间轨道外热流等因素,建立卫星红外辐射通量计算模型,计算得到不同时刻、不同热载荷情况下的卫星红外辐射通量分布,并简要分析了在轨卫星热控涂层衰减所带来的表面太阳吸收比的变化对卫星温度场的影响。     

轨控发动机羽流红外辐射的数值仿真

针对工作在高真空环境下的轨控发动机,数值模拟了其羽流红外特性。首先计算了考虑化学反应的轨控发动机喷管的内流场和外流场,得到了温度、压力、组分浓度等参数的分布。基于HITRAN 2008和HITEMP 2010数据库采用逐线积分法编程计算了气体光谱吸收系数。最后在此基础上用有限体积法建立的求解辐射传输方程的模型,计算得到羽流红外辐射强度在2~10μm范围内随波长变化的曲线。分析了羽流气体组分、波长、探测角度对羽流光谱辐射的影响,与同类文献中的计算结果进行了比较,结果表明：本文的计算模型和方法能较好地模拟轨控发动机羽流的红外辐射特性。     

计算近地轨道航天器空间外热流的RUD方法

（) 摘 要:根据兰贝特定律推导了能束在辐射半球面上均匀分布的外热流计算公式。采用三角形网格统一表示边界面并根据边界面网格进行遮挡判断。根据面元间互视和遮挡关系,对无遮挡表面所受外热流采用积分法进行计算,被遮挡部分的地球反照和地球红外热流采用能束均匀分布法进行计算。计算结果表明,该方法的计算效率和精度都明显高于蒙特卡洛法;互视和遮挡判断方法准确高效。
关键词:中图分类号： 文献标识码：A 文章编号：DOI:收稿日期：2008-05-09; \ 修回日期：2008-06-17     

收敛-扩张喷管中运用次流推力矢量控制技术的计算研究

采用二阶1TVD格式的有限体积法耦合RNGκ-ε湍流模型和非平衡壁面函数求解二维守恒型雷诺平均N-S方程,数值模拟了运用次流椎力矢量控制技术的二维收敛-扩张喷管中的流动现象。计算结果表明,随着喷管压比的增加,次流喷射所产生的激波向扩张段下游移动,且喷管矢量角不断减小到一固定值;次流压比的增加将导致激波向扩张段上游移动,且喷管矢量角不断增加;次流压比对激波和矢量角的影响比喷管压比更强。     

金属热防护系统瞬态传热数值模拟方法研究

针对由金属蜂窝夹心板和多层纤维绝热材料组成的金属热防护系统(Metallic Thermal Protection Sytems,MTPS)建立了气动加热条件下整体传热计算模型,采用蒙特卡罗方法模拟蜂窝腔内的三维辐射传递和每个纤维层内的一维辐射传递,采用有限体积法对MTPS一体化结构内辐射导热耦合换热能量方程进行求解,对某结构的MTPS瞬态传热及热响应进行了数值模拟分析,结果表明在所给定的典型时变气动热流条件下金属蜂窝夹芯板对MTPS瞬态热响应影响很小,增大蜂窝腔内表面发射率可明显降低蜂窝蒙皮间的温度差,再入过程中时变气体压力对MTPS瞬态热响应有明显的影响.     

复杂系统可靠性鉴定方法

为了确定系统可靠性的验前分布，本文采用蒙特卡罗方法，充公利用系统试验之前元件、整机、分系统的可靠性数据获取系统可靠性的验前分布。另外，本文所讨论的系统可以是非常复杂的系统，包括串联或并联系统、串并联系统及不同分布的模型组成的系统。     

超高真空环境下测试产品的放气分析

本文介绍了材料放气的理论和模型,同时运用费克定理推导出圆柱型测试产品的理论放气率,并介绍了测量放气率的两种常用方法:定容法和小孔流导法。文中指出了运用小孔流导法存在的三个问题,并通过双规流导（TGT）法解决了此类问题。     

Research on Residual Stress Measurement of Magnesium Alloy Cabin Castings

Owing to the nonuniform wall thickness and complex internal structure， the measurement of the residual stress on magnesium alloy cabin castings is complex and difficult extremely， and thus seldom research has focused on the residual stress of magnesium alloy castings. In this paper， the blind-hole method， the X-ray diffraction (XRD) method， and the contour method are used to conduct comprehensive and systematic residual stress tests for a magnesium alloy cabin casting. The results show that the residual stress on the surface of the casting obtained by the blind-hole method is between -20.03 MPa and -71.03 MPa， the residual stress obtained by the XRD method is between -26.01 MPa and -87.11 MPa， while the residual stress obtained by the contour method is between -45.89 MPa and 76.87 MPa. The study can lay a basis for the subsequent research of magnesium alloy cabin castings， and provide a reference for the residual stress test of magnesium alloy castings.     

航天器中频力学环境预示研究进展

中频声振力学环境严重影响着航天器的可靠性和安全性,具有响应复杂、预示难的特点,因此航天器结构的中频动力学问题是当前研究的难点和热点。文章介绍了三类主要的中频预示方法,即改进的确定性方法、改进的统计能量法和混合法的基本理论和应用现状。结合航天器工程需求,对比和分析了几类中频预示方法的计算效率、应用前景和理论局限。最后,提出采用有限元与统计能量分析(FE-SEA)混合法作为国内航天器中频力学环境预示方法,并指出了FE-SEA混合法的研究方向。     

基于控制体积方法的降落伞初始充气模型

考虑降落伞初始充气过程中降落伞变形情况与受力特性,建立了基于控制体积方法的初始充气模型。模型从理想气体状态方程出发,通过伞衣内部气体质量求解伞衣内部各个部分所受压力,将压力结果代入基于凯恩方法建立的伞衣结构模型中,即可得到对应的伞形;模型在计算时能够考虑伞衣结构透气的影响。数值计算的结果表明,模型计算得到的初始充气过程伞衣形状变化与试验吻合较好,初始充气时间与经验公式基本一致。该方法还可用于对降落伞主充气阶段的数值模拟。     

基于无线电波坑道透视算法数学模型的研究

本文从地下电磁波远场场强公式出发,采取地质数学常用数学方法,依据最小二乘法原则进行迭代线性拟合,建立数学模型,采取"F-检验法"检验,旨在通过坑透数据处理,准确推断出工作面的地质构造。     

基于JPEG2000的高光谱图像压缩方案

高光谱图像是三维立体图像，具有“谱像合一”的特点，具有很大的应用价值，但因为其数据量庞大，给存储和传输带来不便。文章针对高光谱图像的特点，在进行8：1压缩条件下，提出了差值法和变换法两种基于JPEG2000的压缩方案，仿真结果表明，这两种压缩方案都能较好的满足高光谱图像压缩的要求。     

刀刃法在轨MTF测量性能分析

文章主要研究采用刀刃法测量在轨遥感器的MTF。刀刃法MTF测量准确程度如何，需要有一个定量的客观的评价，文章用仿真的方法对MTF测量精度进行了检验，对于噪声对MTF的影响进行了研究。     

羽流中固体颗粒在真空环境下分布的数值仿真

通过对气?固两相间动量和能量相互作用解耦处理,建立了一种适于模拟真空环境气固两相混合物羽流的DSMC双向耦合算法与固体颗粒空间输运变化特性的TPMC计算技术。仿真了固体火箭发动机两相羽流流场和固体颗粒在远离发动机喷口数十公里空间扩散运动分布特性,通过将计算结果与典型文献结果及理论分析比较确认,证实本文方法的准确可靠性。结果表明,固体颗粒对气相的扩散有一定的阻滞作用;仅在离发动机喷口一定距离,气相对固体颗粒有较大影响,致颗粒温度下降、速度增加;颗粒温度随发动机喷口距离增大而减小,一直要在远离喷口上百公里颗粒温度才随轴向位置趋于平衡,且不同尺寸的颗粒温度差别较大,对指导外层空间高真空环境气固两相羽流传输影响工程研制具有重要意义。     

空间能量粒子探测方法研究

阐述了空间能量粒子探测的原理和方法。介绍了国外在能量粒子探测方面所采用的半导体望远镜测量法、电场加速法、磁偏转法等方法，分析了它们各自的优缺点。就制约能量粒子探测技术发展的关键因素进行了讨论。并以此为基础，结合我国实际情况提出了一种新的探测方法。