抛物型激光推进光船构型设计与特性分析

以轴对称Navier-Stokes方程组为基础,将球面激光支持的爆轰波模型与流场控制方程组进行耦合求解,数值模拟了抛物型激光推进光船的工作过程。计算中对激光支持的爆轰波阵面进行实时追踪,以获得激光能量吸收源项,同时采用Gupta建立的高温平衡空气模型来计算工质的热力学参数和输运特性。并用所发展的数值模拟程序研究了不同构型设计光船的推进性能。最后分析了流场流动特性和推力生成机制。结果表明,推力面离焦点距离越近,光船获得的冲量耦合系数越大,而峰值推力越小。     

光束扩展、光船倾斜和漂移对光船自由飞行演示的影响

为了研究光船的地面演示自由飞行试验中激光束扩展、光船相对激光束的倾斜和漂移等因素对光船飞行高度的影响,采用了理论分析和数值模拟相结合的方法。分析表明,激光束的扩展将导致光船接收激光功率的减小,光船相对于光束的倾斜会影响激光的聚焦和光船推力方向,但光船的倾斜力矩不会在飞行过程中累积而使光船倾斜过大,而光船的漂移则会改变光船推力以及侧力的大小。光船飞行数值模拟仿真表明,在横向偏移扰动下,光船在没有达到最大可能高度之前将脱离光束从而失去动力。研究结果显示,减小光船的横向偏移和激光束的发散角可以有效提高光船垂直上升的高度。     

来流对吸气式激光推力器冲量耦合系数的影响

针对大气吸气模式激光推进的原理,采用二阶精度的Roe格式,在气体模型为完全气体时,分别对无来流以及来流马赫数为2.5和5的条件下,环聚焦激光推力器内外流场的演化过程进行了数值模拟。结果表明,无来流、来流马赫数为2.5和5时的冲量耦合系数分别为2.42×10-4N.s/J,1.00×10-4N.s/J和1.13×10-4N.s/J,无来流时激光推力器的冲量耦合系数高于同工况下有来流的情况,但来流马赫数为5时由于激波不会传播到推力器前体而产生负推力,因此冲量耦合系数较马赫数为2.5时高。研究结果可为激光推力器的构形设计以及实验模拟高空环境下推力器的飞行情况提供一定的理论指导。     

激光功率和工质对连续激光推力器性能影响的数值模拟

对连续激光维持的等离子体加热推力器内流场建立一种计算模型,考虑的物理机制包括化学反应、高温气体性质、激光吸收、黏性、扩散、热传导以及辐射效应,模型方程形式为可压缩轴对称Navier-Stokes方程,对方程采用SIM-PLEC算法求解.在吸收室压强、聚焦光束形状参数及喷管构型相同的条件下,模拟了氩为工质的推力器比冲随激...     

返回舱再入跨流域气动及配平特性数值研究

基于直接模拟蒙特卡罗(DSMC)方法,发展流场直角与表面三角形非结构混合网格生成方法和网格自适应技术,构造适于高稀薄流到近连续滑移流多流区共存的变时间步长模拟策略及DSMC区域分解并行计算方案。采用经改进的DSMC方法模拟飞船返回舱再入过程130～70km跨越外层空间自由分子流到稀薄过渡近连续流区热化学非平衡流动,研究分析了跨流域激波过渡带和边界层的强扰动演变规律以及稀薄气体效应、高温真实气体效应对返回舱气动特性和配平特性的影响规律。对比分析了不同壁面反射模型对返回舱配平特性的影响特点,计算的探月试验返回器配平迎角与飞行试验数据一致。计算分析了质心位置偏移对配平迎角的影响机制。计算结果表明:稀薄气体效应和壁面反射模型对90km以上高度的气动力和配平特性影响显著,真实气体效应则对90km以下高度的气动特性影响较大,质心横向偏移对配平迎角影响较大。     

氩工质连续激光加热推力器的数值仿真

对氩气工质的连续激光推力器内加热与流动过程进行数值模拟,建立了包含激光吸收、化学反应、高温效应、黏性、扩散、热传导以及辐射效应等物理机制的模型.推力室内二维轴对称流场采用可压缩Navier-Stokes方程描述,用SIMPLEC(semi-implicit method for pressure-linked equa...     

一种求解低速跨流域流固耦合问题的离散统一气体动理学格式

运动边界及流固耦合问题是低速及高速、连续流域及稀薄流域流动中常面临的多场耦合问题。本文在原始离散统一气体动理学格式基础上，采用任意拉格朗日-欧拉方法及动网格技术，并耦合计算结构动力学，实现适用于低速流动运动边界及流固耦合问题的计算方法。在方法验证方面，进行连续流域圆柱绕流强迫振动和自由振动数值模拟，计算结果与参考文献结果吻合良好，准确捕捉圆柱在锁频区和非锁频区的流动特性。在稀薄流域范围，将方法中的连续流D2Q9离散速度模型替换成跨流域离散速度模型，对比低速圆柱绕流强迫振动算例在两种离散速度模型下的数值结果。虽主要关注低速流动问题，但发展的算法仍具备模拟高超声速运动边界问题的能力，对此进行了Ma=5.0稀薄流振动圆柱研究。低速和高速运动边界模拟能力显示了方法在考虑稀薄气体效应流固耦合问题方面的潜在应用价值。     

激光推进运载器大气飞行稳定性影响因素研究

针对环聚焦的Myrabo构型激光推进飞行器在大气中飞行现象,建立了激光推力器冲量模型,给出了光船飞行六自由度动力学模型。设计算法解决了脉冲力作用下,光船动力学方程的解算精度与效率折中问题;仿真分析了激光脉冲频率、脉冲能量和升力线斜率大小对光船运动稳定性影响。结果表明:在垂直方向,光船上升速度呈锯齿状,且存在脉冲频率下限;在侧向,激光脉冲频率越高,回复速度越快,回复曲线也越光滑,升力线斜率的增加能够使光船以更短的时间趋于稳定。     

二维过渡区喷流干扰流场的DSMC/EPS仿真

应用DSMC／EPS混合算法，数值模拟过渡区有RCS横向喷流干扰的二维有限平板高超声速流场，分析了流场结构、分离形态及物面气动系数分布等特性。结果表明，不同：Knudsen数下，喷流干扰流场存在差异，过渡区中的稀薄气体效应得以体现，证明了DSMC／EPS混合算法是研究稀薄，连续混合区域高超声速流动的有效方法。     

基于修正N-S方程本构关系的气体动理论耦合方法

随着稀薄程度的增加,Navier-Stokes方程的线性本构关系难以正确描述稀薄气体输运特性,高阶非线性本构关系往往数学形式极为复杂,对数值求解造成稳定性差等问题。为了发展适宜于近空间飞行器气动特性分析的高超声速稀薄流动模拟方法,本文利用求解Boltzmann模型方程的气体动理论统一算法(Gas Kineitc Unified Algorithm,GKUA)对应力张量、热流等宏观量数值积分求解的优势,提出了一种基于数值修正N-S方程本构关系的气体动理论耦合方法。通过将GKUA获得的应力张量及热流用于修正N-S方程的本构关系,实现了存在局部稀薄效应的流动模拟,并且通过可压缩平板边界层、圆柱绕流问题的数值模拟,验证了方法的有效性。     

基于二维激光支持爆轰波模型的光船工作过程数值模拟

以Raizer建立的一维激光支持的爆轰波模型为基础,根据激光推进中焦点区内入射激光强度的分布情况,建立了二维激光支持的爆轰波模型。在相同的激光入射功率下,分别采用球面激光支持的爆轰波模型、球冠激光支持的爆轰波模型和二维激光支持的爆轰波模型与流场控制方程组进行耦合求解,探讨了焦点区内激光强度分布不同时的光船推进性能。结果表明,在吸收相同激光能量的条件下,三种模型所得推力曲线和冲量耦合系数在误差范围内均吻合,说明能量在焦点区的分布对光船推进性能的影响不大。     

激光参数对光船性能影响分析

采用三温度11组元热化学非平衡空气模型计算了激光能量在等离子体中的沉积过程,并在激光脉冲作用结束8 μs后采用平衡空气模型,完成激光推进光船工作过程的数值模拟,研究了不同入射激光能量和脉冲宽度对光船推进性能的影响。结果表明:当脉冲宽度相同时,入射激光能量越大,所得冲量耦合系数越大;当入射激光脉冲能量相同时,脉冲宽度越小,所得冲量耦合系数越大。将计算所得冲量耦合系数与Schall实验所得结果进行比较,结果非常吻合。     

地基激光器发射微卫星的方案设计与弹道优化计算

1引言随着高功率激光技术的发展,用地基激光器发射近地轨道卫星将成为可能。上世纪90年代初Hum-bleW E和Pierson B L提出用一台地基300 MW的CO2激光器把初始质量为1 900 kg的飞行器发射到近地轨道[1]。飞行器的弹道在初始段有一向后运动的“后摆”段。这种弹道可以充分利用激光     

激光推进辐射击穿温度特性的数值研究

在9种等离子体击穿温度下,数值模拟了二次反射式聚焦系统聚焦情况下激光推力器内流场的演化过程,得到了不同击穿温度对应的能量沉积率、推力峰值、冲量和冲量耦合系数,能量沉积率和推进性能参数在某个等离子体击穿温度值处发生突变。根据空气对激光的逆韧致吸收系数公式,计算了CO2激光辐照下不同等离子体击穿温度对应的空气辐射自由程,发现当等离子体击穿温度为14000K时,辐射自由程为1.4mm,与计算网格的典型尺寸相当,此时入射激光能量在一个网格内以一定效率被吸收,由此确定了基于逆韧致吸收的激光等离子体的击穿温度。     

脉冲爆震火箭发动机和小推力液体火箭发动机的性能对比分析

为了了解脉冲爆震火箭发动机的性能优势,对比了脉冲爆震火箭发动机和小推力液体火箭发动机的推力和比冲,其中脉冲爆震火箭发动机的性能计算采用等容循环计算模型.结果表明:真空状态下,随燃烧室进口温度的升高,比冲增加不大;在推进剂和发动机结构尺寸相同的情况下,脉冲爆震火箭发动机产生的推力比小推力液体火箭发动机的多3.0倍至6.8倍,但比冲相当.     

Thrust measurement method verification and analytical studies on a liquid-fueled pulse detonation engine

In order to test the feasibility of a new thrust stand system based on impulse thrust mea- surement method, a liquid-fueled pulse detonation engine （PDE） is designed and built. Thrust per- formance of the engine is obtained by direct thrust measurement with a force transducer and indirect thrust measurement with an eddy current displacement sensor （ECDS）. These two sets of thrust data are compared with each other to verify the accuracy of the thrust performance. Then thrust data measured by the new thrust stand system are compared with the verified thrust data to test its feasibility. The results indicate that thrust data from the force transducer and ECDS system are consistent with each other within the range of measurement error. Though the thrust data from the impulse thrust measurement system is a litter lower than that from the force transducer due to the axial momentum losses of the detonation jet, the impulse thrust measurement method is valid when applied to measure the averaged thrust of PDE. Analytical models of PDE are also discussed in this paper. The analytical thrust performance is higher than the experimental data due to ignoring the losses during the deflagration to detonation transition process. Effect of equivalence ratio on the engine thrust performance is investigated by utilizing the modified analytical model. Thrust reaches maximum at the equivalence ratio of about 1.1.