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

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

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

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

固体火箭发动机羽流红外辐射特性研究

考虑了羽流组分间的化学反应,应用Roe-FDS格式求解羽流场,再用羽流场参数作为输入条件,采用离散坐标法和谱带近似法对固体火箭发动机羽流在光谱2～5μm内的红外辐射特性进行了计算,羽流气相辐射强度随波长变化规律的计算结果与实验值基本一致。计算结果表明,不同羽流位置点的光谱辐射强度随波长变化规律一致,纯气相辐射强度在2.7μm和4.3μm处出现峰值,随羽流场轴向和径向尺寸的增大而减小,在其他波段无明显变化规律;气相与粒子总辐射强度随波长的增大而减小,其辐射强度远大于纯气相的辐射强度;观测方向与羽流轴线夹角的方位角增大,红外辐射强度减小。     

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

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

稳态等离子体推力器羽流数值模拟

以国内研制的稳态等离子体推力器羽流为研究对象,建立了基于粒子单元法和直接蒙特卡罗模拟法的数值仿真模型,简化了电子动量方程,运用Boltzmann关系式求解了等离子体电势,研究了羽流电子数密度、离子电流密度及流场温度分布。通过与实验数据和已有数值结果对比,以验证模型和编制的程序。研究表明,羽流出口下游0.3~1.0 m处的羽流中心电子数密度约为2.5×1016~5×1015m-3,电流密度约为65~1 A/m2。文中工作对分析等离子体羽流污染具有参考意义。     

利用10N推力器羽流试验数据建立羽流场数学模型

利用MBB10N推力器羽流试验数据建立推力器羽流场的数学模型并确定其模型参数,以便利用此数学模型分析此类推力器对卫星的羽流干扰力、干扰力矩以及对太阳翼、敏感器的热效应等羽流污染。本文建立的羽流场数学模型不同于Simons羽流场数学模型及羽流特征线方法,但本质上仍满足羽流远流场为自由分子流的假设条件;本文建立的羽流场数学模型的非线性表达式更便于利用试验数据实现系统辨识,而且可达到相当满意的辨识准度。     

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

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

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

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

航天器发动机羽流对太阳电池板力及热效应仿真

文章将差分求解N-S方程与DSMC方法相结合,研究了航天器单台发动机连续工作时的真空羽流对航天器太阳电池板的力效应和热效应。通过求解N-S方程,获得发动机喷管的内流场;再以内流场计算结果作为模拟粒子入口边界条件,应用DSMC方法,在并行计算机平台上进行三维羽流场和力及热效应计算,得到航天器单发动机连续工作情况下羽流场对太阳电池板的力效应和热效应。文章以太阳电池板处于斜45°状态为例,对仿真结果进行了分析。     

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

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

基于刚体转动四元数模型的航天器再定向方案

介绍了一种新的刚体转动的四元数模型。这些模型可以提高再定向控制的对边界条件的适应能力，在此基础上还介绍了这些模型的一些实际应用。     

Determination of parameters of attitude motion of a small communication satellite using the data of measurements of current of solar panels

A communication satellite (small spacecraft) injected into a geosynchronous orbit is considered. Flywheel engines are used to control the rotational spacecraft motion. The spacecraft after the emergency situation has passed into a state of uncontrolled rotation. In this case, no direct telemetric information about parameters of its rotational motion was accessible. As a result, the problem arose to determine the rotational satellite motion according to the available indirect information: current taken from the solar panels. Telemetric measurements of solar panel current obtained on the time interval of a few hours were simultaneously processed by the least squares method integrating the equations of rotational satellite motion. We present the results of processing 10 intervals of the measurement data allowing one to determine the real rotational spacecraft motion and to estimate the total angular momentum of flywheel engines.     

On solving the problems of optimization of trajectories of many-revolution orbit transfers of spacecraft

The problem of optimal control over many-revolution spacecraft orbit transfers between circular coplanar orbits of satellites is considered. The spacecraft flight is controlled by a thrust vector of a jet engine with restricted thrust (JERT). The mass expenditure is minimized at a limited time of flight. The optimal control problem is solved based on the maximum principle. The boundary value problem of the maximum principle is solved numerically using the shooting method. A modified computation scheme of the shooting method is suggested (multi-point shooting), as well as a method (correlated with the scheme) of choosing the initial approximation with the use of a solution to the optimization problem in the impulse formulation. The scheme and method allow one to construct many-revolution spacecraft orbit transfers.     

热变形对一种天线指向机构指向精度的影响分析

热变形对天线指向机构的指向精度具有重要影响,但利用传统的设计方法很难实现对指向精度的准确计算,尤其是热变形对指向精度影响的研究工作更是少见。为此建立了机构的三维实体热-结构耦合有限元分析模型,根据有限元分析结果并结合机构的运动链关系,构建了各关键运动副运动误差的计算方法;在此基础上实现了机构指向误差的计算,研究了运动副名义、最大以及最小配合间隙条件下,温度载荷对各关键运动副运动误差和指向精度的影响。研究结果表明,温度载荷对关键运动副的运动误差以及机构的指向精度具有重要影响,随着温度载荷偏离常温数值的增加,指向误差不断增大,名义配合间隙条件下的指向误差最大。     

攻角变化对超音速进气道再起动特性的影响

以二维非定常可压缩流的N-S方程为控制方程,采用SST k-ω湍流模型对攻角变化引起的超音速进气道再起动过程进行了数值模拟,研究了攻角变化对超音速进气道再起动特性的影响.结果表明,当超音速进气道不起动时,可通过合理改变攻角实现进气道的再起动工作;超音速进气道的再起动攻角随攻角变化速率的增大近似呈线性增加;攻角变化速率较...     

多自由度微振动环境时域波形复现的数值仿真

为满足航天器微振动环境模拟的需要,开展了多自由度微振动时域波形复现控制方法研究。首先,介绍了基于时域波形复现的多自由度微振动环境模拟控制理论方法。其次,针对六自由度微振动激励系统,应用MATLAB软件建立了基于实测传递函数矩阵的多输入多输出微振动激励仿真系统,针对微振动时域波形复现闭环控制过程进行了算法编程,并给出了仿真的闭环控制流程图。最后,通过算例对多自由度微振动时域波形复现进行了数值仿真,以给定的白噪声为输入,模拟对实际存在的系统非线性、测量误差等影响因素的控制效果。仿真结果验证了多自由度微振动时域波形复现控制方法的可行性及有效性,所得结论可以为研究多自由度微振动时域波形复现控制系统提供参考。     

多通道功率放大器中幅度和相位一致性影响分析

文章利用传输矩阵方法,建立了MPA分析模型,并通过引入幅度和相位随机性,结合各路放大器幅度和相位组合概率,对通道隔离度的影响进行了计算和分析。     

国际空间大学董事会主席访问航天科技集团公司

10月14日,国际空间大学董事会主席拉姆·哈迪姆访问中国航天科技集团公司,就集团公司承办国际空间大学2007年暑期进修项目的筹备情况与集团公司人力资源部部长陈学钏、副部长白燕强交换了意见。陈学钏部长向哈迪姆通报了集团公司筹备这次活动的工作进展情况。据陈部长介绍,中国     

Determination of an Estimation of Parameters of Motion of a Spacecraft on the Stage of Its Descent on the Basis of Integration of Inertial and Optoelectronic Navigation Systems

The problem of determining an estimation of current spacecraft motion parameters on the stage of its descent is solved on the basis of readings of an optoelectronic navigation system and additional information from the self-contained inertial navigation system at the most general assumptions on the nature of motion of the object and the statistical properties of the noise of the measuring devices.