碳纳米管阵列推力器放电特性的数值模拟研究

基于直接蒙特卡洛(DSMC)算法、浸入式有限元算法(IFE)、粒子云及蒙特卡洛碰撞(PIC-MCC)算法以及Ammosov-Delone-Krainov(ADK)隧穿电离模型建立了碳纳米管阵列推力器工质气体和CNT场电离过程的三维仿真模型。利用文献中的实验数据对本文仿真模型的有效性进行了验证,结果表明本文仿真模型与实验结果基本一致。针对碳纳米管阵列推力器的放电特性进行了研究,并重点分析了离子与原子碰撞过程对碳纳米管阵列推力器放电特性的影响。结果表明CNT可以极大地增强局部电场,稳态时场增强因子约为1308;与原子的碰撞使得离子在CNT附近更加集中,降低了局部电场,进而降低了其场增强能力。     

不同直径圆球诱导燃烧的振荡机制与频率特性

为了分析不同直径圆球诱导振荡燃烧的规律，并揭示圆球大小在振荡燃烧现象中所发挥的深层次作用，本文采用二维轴对称Euler方程和基元反应模型，对不同直径的圆球在H2/air预混气体中诱导振荡燃烧的现象开展数值模拟研究。研究发现，振荡频率并不是简单地随直径增大而逐渐从高频向低频连续过渡，而是存在两次突变，形成了超高频、高频以及低频三种振荡燃烧模态。在两种模态间过渡时，振荡达到稳定状态前，会存在一段双频耦合的振荡阶段。三种不同振荡燃烧模态的产生是受到了不同振荡机制的作用，而两种模态间过渡时的双频耦合现象则是两种机制相互竞争的结果。     

火星表面热环境对航天器热控影响分析

火星大气对太阳辐射产生吸收和散射作用，同时还将与火星表面航天器发生对流换热。热设计时难以直接评估对流、辐射和导热三种换热对航天器的影响，从而确定主要的控温途径。在调研火星表面辐射、大气等热环境的基础上，从线性化传热系数和对流辐射比的角度对比分析了辐射、对流和导热对航天器的影响。器表辐射传热系数随光学属性和温度的变化范围为0.3～1.4W/(m2·℃)，对流传热系数随风速变化为0.2～1.5W/(m2·℃)，器内导热传热系数可控制在0.25W/(m2·℃)以下。结果表明，太阳辐射较火星表面和天空辐射而言是主要外热源，航天器表面的辐射和对流换热为两条并联换热途径，两者均可成为主要换热途径，器内导热传热是控制航天器内外隔热的主要可控因素。     

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

跨介质航行器波浪环境入水流场演变和运动特性研究

跨介质航行器是一种既可以在空中飞行又可以在水下潜航的新概念航行器,基于仿生学原理,提出一种通过改变外形实现水空介质跨越的航行器模型,通过入水试验装置和计算流体动力学方法,对航行器带攻角从空气到水的介质跨越过程进行了试验和数值仿真研究,得到了跨介质入水过程航行器的运动姿态和入水空泡形态,并通过数值仿真得到了航行器的升力、阻力、速度和加速度演化规律。同时基于数值模拟方法对有波浪情况和静水情况下航行器入水过程空泡演变以及运动特性进行对比。结果表明:提出的航行器构型在水中具有较好的姿态调整能力,波浪的有无和波高的不同都会对航行体入水运动特性造成影响。     

Initiation characteristics of oblique detonation waves from a finite wedge under argon dilution

In this paper, the flow field characteristics of Oblique Detonation Waves (ODWs) induced by a finite wedge under argon dilution are studied by solving the Euler equations with a detailed chemical model of hydrogen and air. First, the effects of the expansion waves, argon concentration, geometric parameters, and Mach number on the ODW are discussed. The results show that the changes of these parameters may make the oblique detonation not be initiated. Then, the ODW initiation criterion of the finite wedge is summarized, as the characteristic length of the induction zone L_C and the characteristic length of the oblique wedge L_W meet the condition L_C/L_W < 1, the initiation of the ODW occurs; otherwise, it does not occur. What’s more, the Constant Volume Combustion (CVC) theory is applied to study the characteristic length of induction zone. It is found that CVC theory is more suitable for the “smooth transition” type of ODW flow field, the theoretical and numerical characteristic length in induction regions are in good agreement. This work is of great significance for the design of oblique detonation engines for hypersonic vehicles.     

散射对高温气体-碳黑颗粒混合物辐射传输的影响

基于全光谱k分布(Full spectrum k distribution,FSK)模型、MIE理论和有限体积法(Finite volume method,FVM),构建了均温、均质辐射参与性气体-碳黑颗粒混合物介质热辐射传输模型,并分析了碳黑不同尺寸、不同体积浓度以及介质不同路径长度和不同温度条件下,因忽略碳黑颗粒散射所导致的介质热辐射传输特性(如辐射热流、辐射源项)的计算误差。研究结果表明:体积分数不变,增大粒径,计算误差呈现出先增大后减小的趋势;数密度不变,增大粒径,或者粒径不变,增大体积分数,均将使得计算误差相应增大;粒径、体积分数不变,增大路径长度,或者升高介质温度,均将增大计算误差。通常对于含有大颗粒、高碳黑浓度的辐射参与性气体-碳黑颗粒混合物介质,碳黑颗粒散射不能忽略。     

基于有限Fourier级数的航天器转移轨迹设计

针对航天器小推力转移轨迹的初始设计问题，利用基于三阶Fourier级数的设计方法实现了航天器小推力的多圈转移。同时，基于有限Fourier级数的形状法，对具有多个约束条件的小推力多圈转移轨迹进行了优化设计。选取了共面同轴同偏心率的初始和末端轨道位置，对所提出的方法进行了仿真验证。结果表明：与改进逆五阶多项式形状法相比，所提出的方法虽然增加了转移时间，但当转移圈数为5圈、有限Fourier级数的项数为10时，可减少将近75%的转移速度增量，同时大大减小了所需的最大推力加速度的值。     

直升机新型旋翼翼型气动特性与布局分析

利用结构网格计算流体力学(Computational fluid dynamics ,CFD)的翼型气动特性分析方法开展旋翼翼型气动特性计算。通过RAE282，NACA0012，OA212，OA207等翼型压力分布、升力和阻力等特性计算结果与试验结果的对比分析，验证了计算方法的准确性，并进一步完成了HF系列旋翼翼型的气动特性计算。基于翼型的气动特性，采用时间步进自由尾迹的旋翼气动性能分析方法开展旋翼桨叶翼型的气动布局优化设计，对悬停和前飞条件下的旋翼开展计算分析，得到两种条件下的旋翼气动特性。而后通过本文建立的优化方法开展旋翼翼型布局优化设计。     

运载火箭推进剂复杂流动传热问题数值模拟中的模型简化方法

在运载火箭改进优化或新型运载火箭设计过程中,会遇到各种复杂的流动与传热问题。按照真实产品和物理过程直接建立的数值计算模型过于复杂,不同尺度的结构、流动与传热、长时间的飞行过程相互耦合,出现计算时间太长(数年)、收敛困难、程序调试困难等问题。针对初始条件复杂、飞行过程复杂、边界条件复杂3类具体问题,总结了3种相应的模型简化方法,分别是工程经验法、极限参数法、低维近似法,为推进剂复杂流动与传热的数值模拟提供参考。