固体火箭发动机燃烧室冷态突扩气固多相流动的数值模拟

基于欧拉-欧拉双流体方法,建立固体火箭发动机燃烧室轴对称突扩气固多相流动的统一二阶矩两相湍流模型,采用k-ε双方程模拟气相湍流流动和颗粒相湍流雷诺应力kpg方程描述大尺度颗粒湍流运动、小尺度颗粒各向异性弥散特性以及气固相间作用。数值模拟了突扩流动中的颗粒轴向、径向平均速度和脉动速度以及颗粒-气相脉动速度关联轴向分量的分布特性,计算结果与试验吻合较好。研究结果表明:颗粒轴向速度脉动大于径向速度脉动,轴向平均速度梯度影响脉动速度峰值的发生位置和大小;两相脉动速度关联轴向分量大约是径向分量的2倍。     

铝镁推进剂固冲发动机两相燃烧数值模拟

针对铝镁推进剂中心进气固冲发动机,湍流模型采用Reynolds应力方程模型,气相燃烧采用涡耗散模型,两相流采用颗粒随机轨道模型,铝颗粒燃烧采用Brooks燃烧模型,对二次燃烧和流动进行了三维、两相和化学反应流场数值模拟,对比和分析了冲压空气与一次燃气无旋、同向及反向旋转3种工况下补燃室燃烧流动特性及燃烧效率,并对不同工...     

微型固体姿控发动机微喷管内气粒两相流动规律的CFD-DSMC研究

微型固体姿控发动机在航天领域具有广泛的应用前景。以基于MEMS技术的微喷管为研究对象,首先通过计算微喷管中的克努森数,得到了微喷管中的气相流动状态;然后,采用CFD-DSMC方法,模拟了微喷管中的气粒两相流动,并研究了颗粒相质量分数和粒径对气相流动的影响。结果表明,在所研究的来流条件下,微喷管中的连续介质假设是成立的;气相与颗粒相间的动量和能量交换,导致气相马赫数降低、温度升高,同时也导致颗粒相速度增加、温度降低;颗粒相质量分数和粒径均能显著影响气相的马赫数和温度。     

固体火箭发动机喷管内气粒两相流动的CFD-DSMC模拟

气相采用Jameson中心差分格式和显式Runge-Kutta法求解可压缩的雷诺时均N-S方程,颗粒相采用DSMC方法描述,采用PSIC耦合算法模拟了JPL喷管中的气粒两相流动,得到了气相和颗粒相的参数分布规律.结果表明,流场中颗粒相含量越高,颗粒相对气相的阻碍和加热效应越显著,导致气相马赫数降低、温度升高；在相同颗粒...     

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

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

随机轨道模型在喷管两相流计算中的应用

应用随机轨道模型对喷管内气-固两相流动进行了数值模拟,结合实验测量结果对使用随机轨道模型和确定轨道模型两种方法的计算结果进行了对比分析,研究了不同固相参数(固相质量比和粒子平均直径)下的固体火箭发动机喷管性能。研究发现,随机轨道模型对实际流动现象的模拟优于确定轨道模型,随着固相质量比和粒子平均直径的增加,喷管两相流损失增加,冲量系数下降。     

射流推力矢量喷管中粒子运动轨迹的数值模拟

用颗粒轨道模型对基于激波控制的二维收缩-扩张(2-DCD)喷管中的气固两相流动进行了数值模拟,研究了1～60 μm不同直径颗粒的运动轨迹.研究表明,颗粒尺寸越小,粒子的随流性越好;颗粒尺寸越大,主流中的粒子与收缩段壁面碰撞越剧烈,无粒子区越大,二次流中的粒子对下壁面的碰撞越剧烈.     

纳米铁基金属燃料发动机喷管两相流特性数值模拟

针对纳米铁基金属燃料发动机技术,建立颗粒燃烧模型与喷管两相流动控制模型,基于VC++软件自主编程,对纳米铁基金属燃料发动机喷管内颗粒相速度、颗粒温度、颗粒粒径、质量传递速率、颗粒速度滞后、温度滞后、流体密度、压强、雷诺数、马赫数、气相速度、气相温度等参数进行模拟仿真,重点研究了0.4～1.0 μm颗粒粒径和10～40凝...     

潜入式摆动喷管两相内流场数值模拟

气相采用三维薄层近似N-S方程，固相采用拉格朗日坐标系下颗粒轨道模型，计算方法采用显式ENO差分格式和分区多块网格技术，建立了固体火箭发动机内流场模拟软件。对轴对称JPL喷管进行模拟，计算结果与文献进行比较，符合较好。推广到实际三维装药型面时，发现颗粒数目增加导致计算效率降低。采用软件中的气相计算模块对潜入式摆动喷管内流场进行数值模拟，比较了喷管摆动角在发动机喉道附近引起流动变化，研究了非对称流动对喷管内表面压强的影响。     

双燃速星孔药柱长尾喷管发动机三维两相流场数值模拟

采用欧拉-拉格朗日两相方法模拟了带长尾喷管的固体火箭发动机三维两相内流场。在同位网格基础上用有限体积法离散N-S方程,采用不完全LU分解预处理B iCGStab算法求解线性代数方程组,通过可压缩SIMPLE算法求解气相流场。用PSIC方法进行气粒耦合计算,得到了气粒两相的速度、温度等参数的分布以及不同尺寸的A l2O3粒子运动轨迹。计算结果表明,在发动机工作初期气粒两相流场呈现强三维特征,长尾管中后段、喷管收敛段以及燃烧室后封头等部位是烧蚀最严重的部位。     

Numerical investigation on gas–particle flows in horizontal channel under the reduced gravity environments

A unified second-order-moment gas–particle two-phase turbulent model incorporated with kinetic theory of granular flows (USM-θ) is developed to study the particle dispersion behavior of dense gas–particle flows in horizontal channel with 6.96 μm wall roughness and with earth, lunar and microgravity environments, respectively. Anisotropy of gas and particle two-phase stresses and the interaction between two-phase stresses are fully considered by constructing two-phase Reynolds stress model and the transport equation of two-phase stress correlation. The flow behavior of particles in a horizontal channel of Kussin and Sommerfeld [12] experiments is numerically simulated. Results show that the reduced gravity conditions affect the particle concentration distribution, particle velocity and fluctuation velocity, particle temperature, axial–axial fluctuation velocity correlation of gas and particle and particle collision frequency. Under microgravity conditions, particle temperature and collision frequency are much less than those of earth and lunar gravity. Compared with earth gravity, anisotropic of two-phase flow and sedimentation are weaker.     

高温固体颗粒气动阻力系数数值计算

在发动机喷管的气粒两相流场中,金属颗粒阻力的准确计算对整个流场的模拟有着重要作用.而高温颗粒与周围气相的传热会改变气相的性质,进而影响颗粒的气动阻力.采用k-ε湍流模型和强化壁面处理方法,计算了雷诺数Re在0～10 000范围内,颗粒与气相不同温差情况下球形单颗粒的气动阻力系数.颗粒温度升高,颗粒所受的气动阻力随之增加,但是增加的幅度随Re的增加而减小;根据计算结果,拟合了以颗粒雷诺数为自变量的高温颗粒气动阻力系数的计算公式.通过与低Re下的实验数据的对比,验证了所拟合公式的适用性和可信性.     

不同重力环境下空间机械臂神经自适应鲁棒控制

针对空间机械臂从地面装调到空间应用过程中重力项的变化问题,提出了一种神经网络自适应鲁棒补偿控制策略用于空间机械臂的末端控制,从而实现在地面重力环境下装调好的空间机械臂在空间微重力环境下实现在轨操控任务。通过神经网络在线建模来逼近系统模型中变化的重力项,逼近误差及系统的不确定性通过自适应鲁棒控制器来补偿。该控制策略不依赖于系统的模型,避免了回归矩阵的复杂计算及未知参数的估计,降低了计算量。基于李亚普诺夫理论证明了闭环系统的渐近稳定性。仿真结果表明该控制器对不同重力环境下空间机械臂的末端控制均能达到较高的控制精度,具有重要的理论研究和工程应用价值。     

防空导弹直接力／气动力复合控制系统设计

直接力和气动力复合控制技术是先进防空导弹对付高速机动目标、兼有反导能力以 及发 展成为动能拦截弹最现实有效的技术途径。针对力矩式复合控制方法,为减小由于弹体滚动 所引起的耦合影响,建立了俯仰－偏航通道短周期运动模型;在弹体慢旋条件下给出了参考 信号状态空间模型,利用线性二次型最优跟踪控制方法设计了复合控制姿态稳定系统,并基 于分离原理设计了参考信号状态观测器。仿真结果表明系统对制导指令的响应时间小于0.1 秒,对螺旋机动弹道目标的拦截脱靶量为0.34米,说明该复合控制方法具备直接碰撞目标的 能力。〖JP〗     

燃气喷射角度对含硼固体火箭超燃冲压发动机补燃室燃烧效率的影响

固体火箭燃气超燃冲压发动机具有高比冲、结构简单、流量易调节等优点,然而在超音速空气流的补燃室中,如何让燃料更好地与空气掺混,增加颗粒停留时间,在较短时间内释放出更多的燃烧焓成为目前研究的重点。采用Realiazble k-ε湍流模型,单步涡团耗散模型,在King的硼颗粒点火燃烧模型的基础上考虑了硼颗粒在高速气流当中的气动剥离效应,利用龙格-库塔算法迭代计算硼颗粒点火燃烧过程,对燃气进气方向与轴向夹角从45°~180°的10种进气方式下的补燃室进行了三维两相燃烧流动计算,分析了各种进气角下的燃气燃烧效率、硼颗粒燃烧效率以及总燃烧效率。结果表明:当一次燃气喷射角度与轴向夹角逐渐增加时,燃气与颗粒燃烧效率逐渐增加,并在180°时燃烧效率和比冲为最高。     

Numerical investigation on hydrodynamics of gas-binary particles flows using a second-order-moment turbulent model

A binary-particles Reynolds stress equation is developed on the basis of single second-order-moment turbulence model. In this model, particle–particle two-phase fluctuation velocity correlation transportation is proposed to fully reflect their anisotropic characters, as well as each of particle interaction with gas phase is taken into account. Swirling gas-particle flow experiments data is used to validate this model simulation in term of the reduced codes of the mono-disperse particle gas-particle flows. Numerical simulations are performed for the gas-binary particles mixtures turbulent flows and results shows that the axial Reynolds stress in near wall regions is approximately 3.0 times greater than the tangential direction and the tangential Reynolds stress in central axial region is approximately 5.0 times greater than the axial direction. The Reynolds stress in the axial direction and tangential direction are redistributed due to binary particle–particle collision. The ratio kinetic to particles collision contributions in binary mixture system are varied with the different particle density and the diameter size compositions.     

对欧空局轨道碎片模型的评价

文章评价了颗粒大于1cm的碎片模型,它包含大于10cm的一组粒子和介于1～10cm之间的粒子,它们是由空间暴露模拟得到的。接着讨论了尺寸范围介于0.1～10mm之间的小粒子束流。假定这些粒子主要是由小粒子与卫星碰撞产生的。这种碰撞主要发生在450～500kin的高度范围内(空间站高度),碰撞与否还取决于轨道高度和离心率。     

空间拦截攻击区和威胁区仿真研究

基于冲量变轨理论,建立了天基动能拦截器通过一次变轨机动拦截目标的轨道机动模型,提出了攻击区和威胁区的概念,并给出了确定攻击区和威胁区的数值搜索算法.通过仿真计算,分别对不同的拦截时间和机动能力约束条件下攻击区和威胁区进行了仿真分析,结果表明,利用本算法能够方便地得到动能拦截器的攻击区和目标的威胁区,为拦截的实现提供必要的决策条件.     

滑盘式流量调节燃气发生器动态特性分析

分析了滑盘阀流量调节燃气发生器的工作原理,建立了滑盘式流量调节燃气发生器工作过程的动态模型,得出了燃气发生器的流量传递函数,并分析了系统增益、时间常数的影响因素,并通过实例仿真验证了理论推导结果.结果表明,调节系统为一个非最小相位系统,具有负调特性,为调节系统的进一步设计提供了参考.