锥柱型装药固体火箭发动机两相内流场中颗粒运动的数值模拟

研究锥柱型装药固体火箭发动机(SRM)燃烧室和喷管的内流场具有重要的工程意义。针对某型发动机各工作时间下的不同情况,采用高雷诺数下的k-ε湍流模型和欧拉-拉格朗日两相流模型,用全速度SIMPLE方法对方程组进行求解,并用PSIC气固耦合计算方法,对其内流场进行了三维两相一体化计算。最后给出了流场内部分颗粒运动的轨迹图像及颗粒与壁面碰撞的一般规律。在不同的时间下,不同尺寸的颗粒会撞击壁面不同的位置,但主要集中在壳体与喷管的相接部分和喷管的收敛段。局部的旋涡也会对颗粒的运动以及颗粒的碰撞沉积、聚集等造成较大影响。     

激波风洞边界层强制转捩试验研究

针对升力体模型设计了涡流发生器,在中国空气动力研究与发展中心(CARDC)Φ2 m激波风洞上开展风洞试验,研究了高超声速边界层强制转捩问题。试验来流名义马赫数分别为10、12,单位雷诺数分别为2.4×10~/m、2.1×10~6/m,模型攻角10°。试验中应用铂薄膜热流传感器技术和温敏热图(TSP)技术测量了模型表面热流,证明涡流发生器实现了模型边界层强制转捩,使Φ2 m激波风洞拥有了模拟高马赫数低雷诺数湍流边界层的能力。试验结果表明,不同形状不同高度涡流发生器对边界层完全转捩成湍流后的热流影响不明显,由此可提出一种新的激波风洞试验方法,即利用涡流发生器开展相同来流条件下不同边界层流态对模型表面热流等边界层参数分布影响的试验研究。     

轴对称后向台阶不稳定流动及压强振荡数值研究

使用翼柱形装药的战术发动机,在翼槽烧完时会出现轴对称的后向台阶结构,易产生转角涡脱落引起的声涡耦合现象。根据战术发动机中后向台阶流动特点设计简化模型,从线性稳定性理论出发,应用声涡耦合模态预测方法进行分析。采用湍流流动的大涡模拟(LES)方法,进行了较宽来流Re数条件下后向台阶不稳定流动的数值模拟,获得了转角涡脱落产生的一般过程,分析了不同来流条件下的旋涡多尺度运动规律。通过分析流场中关键位置的压强振荡,研究了气动声学现象。结果表明,声涡耦合模态预测方法可对不同速度下可能诱发的声场固有频率进行合理预测。流场头部压强波动的振幅最大,频率与声场轴向声模态相耦合。再发展区出现低于轴向基频的宽峰振荡主要受大尺度旋涡运动的影响,该频率与基频接近时,会造成基频的多峰现象。当近壁面小尺寸旋涡形成频率接近声场某阶固有频率时,该频率及更低阶频率就有可能被激发。     

凝胶推进剂模拟液直圆管流动特性初步研究

对某型凝胶推进剂模拟液在直圆管中的流动特性进行了理论分析和试验研究,得到模拟液在直圆管中的速度分布,定义了表观粘性和平均粘性,进行了流阻分析并定义了非牛顿流体雷诺数。通过试验考察了管长、管径对流阻的影响,得到流阻系数与雷诺数和无量纲速度比与雷诺数的关系,并对影响因素进行了分析。     

重力对载人航天器密封舱内空气对流换热影响的数值分析

地面重力环境中进行航天器密封舱内空气通风换热试验时,由于自然对流的存在导致换热量和温度分布与空间微重力环境中的情况存在偏差。文章针对航天器密封舱,建立了舱内空气对流换热的数值模型,利用数值模拟软件对有无重力时典型工况下的对流换热进行了数值模拟及模拟结果的对比分析。分析表明重力对壁面换热量的影响较大,而对空气温度及分布的影响较小;且重力的影响随空气与壁面温差的增大而增大,随通风流量的增大而减小,舱间通风也会减小重力的影响。因此在重力环境中进行试验时需要对壁面换热量进行修正。     

基于焓平衡法的固体燃料冲压发动机燃速预示模型

修正了固液混合火箭的固体燃料燃面退移速率预示模型中对流换热和汽化热的计算,将其用于固体燃料冲压发动机中,考虑了固体燃料冲压发动机中来流空气总温及装药长度的影响。编写了固体燃料冲压发动机燃面退移速率预示程序。计算结果表明,当来流空气总温大于450 K时,计算误差的平均值不超过6.55%,所完善的模型的计算结果与试验值吻合很好,该模型能准确预示固体燃料冲压发动机中燃面退移速率;燃面退移速率与装药长度的-0.23次方成正比。     

积炭辐射对固体燃料燃烧的影响

本文给出一个研究积炭辐射对固体燃料燃速影响的综合模型,该模型由几个子模型组成:用低雷诺数二方程k—ε模型计算流场输运问题,用旋涡—耗散模型计算湍流及化学反应问题,用二热流方法计算辐射传热问题,除对流和扩散外,输运方程还考虑了积炭的生成及氧化问题,以计算其浓度。该模型对PMMA燃料所作的计算与试验数据进行了对比,表明积炭辐射是确定固体燃料的热解速率最重要的因素。     

高超声速飞行器前缘缝隙流动数值模拟研究

通过分析高超声速飞行器前缘防热瓦结构,建立了一种开缝前缘的简化模型。针对这一模型的流场通过求解三维可压缩Navier Stokes方程进行了数值模拟。研究了缝隙诱导形成的三维旋涡的空间分布特征和旋涡运动对物面气动加热的影响规律。模型圆弧段缝隙肩部倒圆区因存在较强的三维效应形成“常规”高热流区,而缝隙内主旋涡再附致使侧壁上存在一个“非常规”高热流区;模型平直段展向流动诱导缝隙上方出现较强的旋涡运动,同时流动在缝隙倒圆区形成分离涡并于缝隙侧壁面再附,受这些旋涡运动的影响,缝隙肩部倒圆区转变为局部热流低值区,缝隙侧壁上存在局部热流高值区。     

高超声速气动加热关联方法的适应性分析

针对高超声速飞行器热环境预测上有着重要应用的关联方法缺乏构造机制和适应性分析的问题，从驻点热流的Kemp-Riddell公式出发，推导了驻点斯坦顿数和雷诺数之间的经典关联关系，从理论上分析了上述关系的适应性。结果显示斯坦顿数和雷诺数关联关系受马赫数的制约，飞行条件和试验条件马赫数变化范围较大的情况下，更适合采用雷诺数和马赫数的组合变量开展关联分析。当试验来流静焓较低的情况下，过低静焓对热环境的关联特性有不利影响，需要提高试验来流马赫数，远离临界马赫数确保地面数据同飞行条件之间的关联性。在国内4座风洞上开展的风洞试验和DSMC/CFD的数值仿真分析工作，校验了上述结论在连续流区域的正确性。     

液氧贮箱增压过程研究

采用数值模拟方法对液氧贮箱增压过程进行研究.贮箱内流场采用流体体积函数（VOF）多相流模型考虑,选择标准双方程k-ε湍流模型分析湍流效应,气液两项之间的热量、质量转移通过自定义程序（UDF）求解.获得了贮箱压力、排液流量、气垫温度、液氧温度对贮箱内流场温度分布的影响.计算结果表明,在稳定增压过程中,贮箱液面无扰动,贮箱内温度分层分布;各参数变化时,对贮箱内温度分布的影响主要是温度梯度的变化,并且各工况下液面附近和扩散器附近温度梯度基本相同.     

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

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

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.