大瑞利数下二维Benard对流演化过程的DSMC仿真研究

利用DSMC（Direct Simulation Monte Carlo）方法研究了大Rayleigh（瑞利）数条件下稀薄气体的二维Benard对流的特性。考察了流动的非定常演化过程及若干流动特征量随时间的变化特征，得到了与实验较为一致的结论。计算结果表明，当流动控制参数不断增大时，二维Berard对流从热传导状态过渡到热对流状态表现出明显的非定常演化特征。最后当瑞利数增大到一定程度后，注动显示出混沌的特性。     

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

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

稀薄气体圆柱绕流的并行化DSMC方法研究

本文针对ＤＳＭＣ方法的内在并行特征构造了两种并行算法,并在ＰＶＭ环境下实现了并行化计算。计算结果表明,实现ＤＳＭＣ方法的并行化能够显著缩短计算时间,提高计算效率,不同并行算法能够达到的加速比存在差异。     

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

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

外激振动条件下气体减压器工作稳定性仿真研究

为研究运载火箭增压系统中减压器在外激振动条件下的工作稳定性,在减压器模型中通过壳体位移坐标引入振动激励源,建立了可模拟减压器在外激振动条件下工作过程的计算模型,对外激随机振动条件下增压系统动态工作过程进行了仿真.针对仿真发现的减压器下游压力振荡的不稳定现象,对减压器下游不同管路尺寸和减压器不同参数下的系统工作过程进行了对比计算.计算结果显示:增大减压器下游管路尺寸,减小减压器阀芯组件质量与减小阀芯组件与壳体间滑动摩擦力,都可以使压力振荡过程减弱,有利于改善减压器在外激振动过程下的工作稳定性.计算结果为减压器及增压系统的设计和优化提供了参考.      

燃气涡轮发动机压缩系统稳定性数值模拟的研究

论发展了一种预测航空燃气涡轮发动机压缩系统稳定工作边界的数学物理模型。其中压缩系统流动过程的动态模型采用基于平行压气机理论的准一维时间相关模型方程；叶片排对流动过程的影响则应用“激盘－滞后－容积”模型；压缩系统稳定性模型则采用了对整个压缩系统统一判别的方式。计算过程与发动机非设计性能计算相关联，使得压缩系统稳定性分析在真实的发动机运行环境（调节规律）下进行。在均匀进气条件下应用这种方法，对一台九级压气机的涡轮喷气发动机和一台三级低压风扇、五级高压压气机的低涵道比、混合排气、带加力的涡轮风扇发动机在均匀进气和畸变进气条件下的稳定工作边界，分别进行了数值模拟和分析。结果表明，发展的数学物理模型和计算机软件系统可以正确的反映发动机压缩系统的工作状况，用它判别发动机不稳定工作点的重复性和灵敏度都比较好。     

基于数值仿真的时域内人-机闭环系统稳定性研究

针对带电传操纵系统人-机闭环系统稳定性问题,讨论了人-机闭环系统稳定性的研究方法。在分别建立飞机本体、操纵系统和驾驶员模型的基础上,建立了综合考虑非线性因素的横航向补偿控制模型。以横向滚转操纵为例,在时域内对系统进行了动态仿真,研究了驾驶员滞后时间对人-机闭环系统稳定性的影响。仿真结果表明,该方法用于研究人-机闭环系统稳定性是可行的和有效的。     

结冰条件下人-机-环系统的飞行风险概率

以结冰条件下的飞行风险量化概率为研究对象,基于蒙特卡罗飞行仿真实验对结冰条件下人-机-环系统的耦合特性进行了分析,并获取了飞行参数极值样本。构建了飞行风险发生的判定条件;对飞行参数极值样本进行了统计特性分析,验证了其厚尾分布特征。一维分布类型辨识结果表明广义极值分布对相对速度和迎角极值的描述精度最高。为描述二维变量对相关性的各自影响程度,提出了一种新的双参数变权重Copula模型;辨识结果表明该Copula模型能以较高的精度通过假设检验。相关性分析的结果表明相对速度和迎角同时出现极大值和极小值的概率较大。基于二维极值样本的Copula分布模型求出了不同结冰程度下的飞行风险概率值,探讨了飞行风险的非线性增长趋势。     

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

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

一类二维非结构网格DSMC方法的实现策略及其应用

研究了一类二维非结构网格DSMC方法的实现策略。在数据结构方面,设计了局部化的数据组织方式,节约了内存与计算时间。发展了一种跟踪模拟分子迁移的算法,该算法仅需少量的逻辑运算与代数运算,不仅可以快速跟踪模拟分子在网格之间的迁移,而且可以准确判别分子与物面是否相互作用,搜索过程中的附带信息给出了分子与物面碰撞的精确时间与位置,避免了重新计算。引入碰撞距离的思想,既减少存储又保证正确的模拟结果。在程序编制过程中,我们充分展现了Fortran90高级语言的主要特性,引入动态数组、指针、链表以及派生类型数据,编制了计算程序。最后对过渡流域高超声绕流进行了数值试验。     

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

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

超声速平头圆柱绕流DSMC/EPSM混合算法研究

本文将EPSM算法与DSMC方法结合，构造了可模拟含近连续流区及过渡流区的DSMC／EPSM混合算法。运用混合算法模拟马赫数等于5．37时超声速平头圆柱绕流，并与DSMC结果进行比较，证明了DSMC／EPSM混合算法的有效性，并对两种算法的计算效率进行了比较。     

一种稀薄两相流动的数值模拟方法

从描述粒子运动的微观层次出发,采用双向耦合技术,建立了一种适用于稀薄条件下两相流动的DSMC数值模拟方法。对相间相互作用进行解耦处理,实现了气固两相间动量和能量相互作用的模拟。相关算例验证了方法的可靠性。结果表明:在大样本统计平均条件下,该模拟方法能够保证相间动量和能量交换时守恒,固体颗粒尺寸、气体稀薄程度等条件对相间的相互作用有很大的影响。     

过渡流区N-S/DSMC耦合计算研究

在已有CFD和DSMC方法程序研究基础上,采用MPC耦合处理原理,发展适于流场分区信息交换的亚松弛技术,建立了N-S/DSMC耦合算法.通过对近连续过渡流区不同Kn数条件下钝锥绕流计算与低密度风洞试验结果比较分析,一方面证明了本文所建立的N-S/DSMC耦合计算模型和方法在低Kn数过渡流区的有效性,另一方面剖析了N-S/DSMC耦合计算误差,可望发展一个工程适用的近连续过渡流区高超声速飞行器气动特性预测分析工具.     

高超声速平头圆柱绕流化学反应流场的DSMC/EPSM混合算法研究

本文发展了平衡粒子模拟方法(EPSM)，建立了与高温气体化学反应动力学理论相匹配的EPSM耦合模型，并通过混合参数进行流区的自动识别，将EPSM方法与蒙特卡罗直接模拟方法(DSMC)结合，构造了可模拟化学反应流动的DSMC／EPSM混合算法。应用该算法对汲及化学反应的轴对称情况下高超声速平头圆柱绕流流场进行模拟，将结果与DSMC方法的结果进行比较，验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较，发现混合算法能够大大提高计算效率。