有限平板绕流Monte—Carlo方法仿真

本文采用直接仿真Monte-Carlo方法求解有限长平板绕流等问题。本方法是通过计算机跟踪仿真分子的运动来实现数值模拟的。仿真分子间的碰撞计算由统计抽样确定。碰撞模型分别选用了硬球分子模型和负幂律分子模型。仿真分子与固壁作用采用由完全扩散反射和镜面反射按比例混合组成的模型。为了检验方法的可靠性,还计算了激波结构和Rayleigh问题等一维流动,二维计算采用同步并行程序。数值结果表明直接仿真Monte-Carlo方法能够较好地模拟稀薄气体力学中的一些问题。对于二维计算,所花费的机时和所需的内存均在国内机器所允许的范围内。     

基于稀薄气-粒两相流的蒙特卡洛颗粒辐射模型研究

在气-粒两相相变模型及液态和固态颗粒碰撞、聚合和分离模型的基础上,发展稀薄条件下考虑颗粒辐射的蒙特卡洛颗粒辐射模型。通过对高超声速稀薄环境中的气-粒两相喷流流场的数值模拟,得到气-粒两相流的流场参数,利用所得流场参数作为颗粒辐射模型的初始参数进行颗粒辐射计算,同时考虑了有无探照发射时的光谱辐射强度。结果表明,在颗粒浓度较大时计算两相稀薄流的流场参数,考虑颗粒辐射是必要的,并且考虑有无探照发射对光谱辐射强度数值的影响。     

气动物理地面模拟试验相似律

依据流动子过程特征量分析方法，导出了高超声速化学反应流动的相似律。将得到的相似律应用于气动物理地面模拟试验，给出了确定地面模拟试验条件的相似准则，分析了地面模拟试验中需要进一步解决的其它基础理论问题。     

过渡区球锥气动力特性的桥函数估算与DSMC仿真研究

采用DSMC(DirectSimulationMonte-Carlo)方法模拟球锥体过渡区流动。论文给出了不同Knudsen数球锥绕流流场的温度与压力分布图,以及在不同Knudsen数、不同马赫数条件下球锥阻力系数和底部阻力系数的分布,并将阻力系数模拟结果与桥函数一体化估算结果进行了比较分析。     

再入三维高超声速粘性流场的数值模拟

基于三维可压Favre质量加权平均N-S方程，采用B-L湍流模型模拟了高超声速条件下，球头双锥再入体周围的流场特征。利用一种改进的时间、空间二阶精度的NND格式，再现了超声速／高超声速条件下再入体周围的复杂流动现象，分析、对比了层流和湍流流场流动现象的差别。     

球头双锥再入体表面脉动压力环境的分析与预测

本文根据超声速和高超声速流动特征分析了球头双锥再入体的表面脉动压力环境特性，并根据超声速和高超声速流动情况下球头双锥的表面压力分布给出一套预测球头双锥再入体表面脉动压力分布的工程估算公式，计算结果表明，预测的均方根脉动压力分布与实验结果基本一致。     

过渡区球头气动力特性DSMC仿真研究

采用DSMC（Direct Simulation Monte-Carlo）方法模拟三维球头过渡区流动。文章给出了球头绕流流场的温度与压力分布图,以及在不同Knudsen数、不同马赫数条件下球头表面压力系数分布曲线,并将阻力系数模拟结果与桥函数公式进行了比较分析。     

气动物理地面模拟试验相似律

依据流动子过程特征量分析方法 ,导出了高超声速化学反应流动的相似律。将得到的相似律应用于气动物理地面模拟试验 ,给出了确定地面模拟试验条件的相似准则 ,分析了地面模拟试验中需要进一步解决的其它基础理论问题。     

高超声速稀薄流的气粒多相流动DSMC算法建模研究

基于直接模拟Monte Carlo(DSMC)方法,构造适用于DSMC算法的固态和液态颗粒碰撞、聚合和分离模型,发展稀薄条件下双向耦合作用的气粒多相流的DSMC算法,在此基础上初步实现高超声速稀薄流环境中的气粒多相喷流流场数值模拟.算例结果表明该方法能为稀薄过渡区气粒多相流动提供一种新的应用研究手段.     

稀薄流热化学非平衡效应的DSMC研究

以钝体绕流驻点线流场为例，应用直接仿真蒙特卡罗法，模拟了五组元混合气体（Ｎ２，Ｏ２，Ｎ，Ｏ，ＮＯ）中的碰撞能松弛，离散反应，交换反应及复合反应，给出了碰撞分子各自由度间的能量交换模式及各类化学反应的选择模式。ＤＳＭＣ仿真结果与宏观粘性激波层方程计算结果，在飞行高度较低时基本相符，表明采用的热化学非平衡非模式是令人满意的，随着飞行高度增高稀薄气体效应增大，二者在激波附近有显著差异，表明高空飞行必须考