热辐射计内稀薄气体流动的蒙特卡罗模拟

本文应用直接仿真蒙特卡罗(DSMC)方法,研究了Crookes热辐射计内由于叶片两表面的温度不同而引起的稀薄气体的流动,给出了热辐射计叶片高温面温度900K,低温表面温度为300 K,四种不同克努森数(Kn)下流速矢量场、温度场和叶片两面压力差的仿真结果.并与热辐射计在不同Kn数下的转速实验测量结果进行了比较,两者结果基本一致.     

超声速平板绕流DSMC/EPSM混合算法研究

本文将平衡粒模拟方法（EPSM）与蒙特卡罗直接模拟方法（DSMC）结合，构造了可模拟含近连续流区及过渡流区的DSMC／EPSM混合算法。运用混合算法模拟马赫数等于4时超声速平绕流，并与DSMC结果进行比较，证明了DSMC／EPSM混合算法的有效性。将DSMC方法与混合算法的计算效率进行了比较。     

DSMC量子动理学模型在火星再入流动中的应用

为解决化学反应模型高温数据缺乏的难题，探索DSMC方法量子动理学（QK）模型在实际中的应用，本文将该模型进一步应用于火星探测器稀薄气动特性的数值预测。通过计算探路者号在85 km、95 km和110 km高度的稀薄绕流，评估了QK模型的性能和稀薄气体效应的影响规律。结果表明，QK模型不依赖宏观的化学反应速率系数，适用于火星再入流动计算。化学反应及其模型对气动力的影响很小，但对气动热特性的影响不容忽略，考虑化学反应后的驻点热流可以下降约12%~14%。     

合并式电离—复合反应的DSMC模拟

采用ＤＳＭＣ方法研究分子碰撞反应机理,在分析水平上模拟绝热封闭体系事并式电离－复合反应的驰豫过程。     

稀薄气体条件下Rayleigh-Benard系统稳定性的DSMC仿真研究

本文利用DSMC（Direct Simulation Monte Carlo）方法研究了稀薄气体条件下Rayleigh-Benard系统的稳定性问题。考察了在经典的Boussinesq假设不再成立时，Rayleigh-Benard系统稳定性问题的特征。研究指出在气体的稀薄条件下，初场密度梯度的方向和Rayleigh数是影响Rayleigh-Benard系统不稳定性的主要因素。     

基于直接模拟蒙特卡洛方法的真空羽流不确定量化研究

真实流动环境下的真空羽流必然存在着各种不确定性,那么确定性输入条件的数值模拟必然会存在偏差,因此需进一步研究不确定性对羽流流动特征的影响规律。本文采用直接模拟蒙特卡洛(DSMC)方法,对不确定性输入的羽流流场进行模拟;采用稀疏的概率配置点方法对来流、壁面及模型参数等输入不确定量进行描述,对不确定性的传播和输出目标的平均值、方差及不确定度进行计算。研究表明,流场不确定性沿流线传播至流速最大处之后迅速增强,并在声速线前出现骤减的现象;传播至声速线之后,挡板壁面输入不确定性的影响凸显。其中最为显著的是,压力不确定度在挡板驻点位置达到全场最大值,约为输入不确定量(3.54%)的2.1倍。此外,温度跳跃不确定度受到壁面温度不确定性输入的限制而近似保持为一个恒定值,约为输入不确定度的0.8倍。进而,壁面热流不确定度(5.54%)比壁面正应力不确定度(6.25%)略小,切应力不确定度最小(5.07%)。Sobol’全局敏度分析表明,喉道速度和喉道压力的输入不确定性对气动力/热不确定度的贡献是最大的,且远远超过了壁面温度和模拟分子直径不确定性输入的贡献。     

高温高速稀薄流的DSMC算法与流场传热分析

在非结构网格下,研究了热力学碰撞传能的6种方式和化学反应的8种类型,详细给出了上述碰撞传能和化学反应碰撞类型的子程序计算框图.在可变硬球(VHS)分子模型、Borgnakke-Larsen唯象模型、Bird的化学反应几率模型以及壁面CLL(Cercignani-Lampis-Lord)反射模型的基础上,用Fortran语言编制了能够模拟内能松弛、热力学非平衡和化学非平衡的稀薄气体DSMC(direct si mulation Monte-Carlo)源程序,在地球大气层和火星大气层中完成了Ballute减速装置的8个工况计算(其中Knudsen数从0.05变到30.0,飞行Mach数从26.3变到11.2),并与NASA Langley研究中心2007年发表的计算结果作了比较,本文的结果令人满意.另外,这里用|T-T_v|/T去分析热力学非平衡,用Damk hler数去分析化学反应非平衡,用Stan-ton数去分析飞行器壁面的传热效果,所有这些分析对高温部件的热防护设计十分有益.      

脉冲等离子体推力器羽流的一体化粒子模拟

将DSMC(Direct Simulation Monte-Carlo)/PIC(Particle in Cell)流体混合算法与一维MHD双温放电模型相结合,一体化模拟NASA Glenn PPT羽流。对不同出口偏转角的羽流场进行模拟,给出了不同偏转角度下离子、电势、温度的变化情况,并与实验结果进行了比较。结果显示,该模型具有一体化(从工作过程到羽流)预测脉冲等离子体推力器羽流的能力。人工偏转角的加入,增强了羽流的扩散,使得羽流径向扩散、轴向收缩、回流减弱,同时对放电模型带来的不足有一定改善作用。     

空间发动机羽流场的试验研究

为研究空间发动机羽流对卫星的影响,建立羽流污染预报系统,研制了一套羽流试验系统并进行了试验研究.试验台主要包括空间环境模拟系统,电热气体发动机,稳压气源,轴向和径向移动装置,测量系统和温控系统.系统可在满足10-3Pa量级的真空压力和93±5K的背景温度下稳定运行和精确测量.本文完成了表征羽流特征的关键参数-压力场的测量.同时采用CFD和DSMC相结合的方法进行羽流场的数值模拟计算.测量数据与数值计算结果以及相似试验条件下Boyd等人的测量结果进行了对比.     

氧氮混合气体在梯度磁场中流动的数值模拟

用直接蒙特卡洛方法模拟了空气气流在梯度磁场中的氧浓度分布，利用DSMC模型将分子的碰撞和运动解耦，同时将磁场力简化为瞬间作用力，并对模拟结果进行了分析和对比。模拟表明，磁场强度与梯度乘积由100T^2／m增加到800T^2／m时，氧浓度由0．762％增加到3．1％；而当温度由0℃升高到100℃时，氧浓度由0．78％减小到0．25％。另外，随着压力的变化氧浓度存在一个最大值，在模拟条件下，这个最大值出现在0．05MPa(绝对压力)附近。