台阶后紊流场的大涡模拟

用大涡模拟法对台阶后复杂的紊流场进行了数值研究。用有限体积法求解滤波N－S方程，亚格子紊流模型采用结构模型。数值计算中，对流项，扩散项和时间分别采用三阶SHARP中式，二阶中心差分格式和三阶Runge-Kutta格式，给出了台阶后紊流场的时间统计特性，瞬时特性和紊流强度分布，并与试验结果进行了比较。     

基于蒙特卡罗射线追踪的航天器空间环境因素暴露通量统一计算方法

航天器空间环境效应分析评估的首要条件是获取空间环境因素的暴露通量。文章分析总结了蒙特卡罗方法在原子氧、微流星体、太阳光压、气动力、辐射剂量、材料放气污染等环境因素暴露通量计算中的应用,提出基于蒙特卡罗射线追踪方法实现统一的空间环境因素暴露通量计算,并给出计算流程及初步验证。统一算法的计算结果与理论值(或文献值)符合较好,证明其能够正确计算环境因素暴露通量。     

卫星寿命指标分配方法研究

针对卫星寿命要求如何分解的问题,提出了一种卫星寿命指标分配的方法,在确定卫星耗损寿命要求的基础上,通过开展寿命建模,开发了一种基于蒙特卡罗仿真迭代的寿命指标分配算法,实现了卫星寿命指标的量化分解。该方法可在当前卫星寿命设计工作中推广应用。最后,文章对卫星长寿命设计工作提出了相关建议。     

空投自导鱼雷武器时命中概率计算模型研究

阐述了用蒙特卡罗法计算自导鱼雷武器命中概率的方法，文中避开繁锁的数学公式推导过程，着重说明建模的思路及建模中对关键问题的处理办法，并对建立的模型作了分析和讨论。     

高能(450 kV)X射线参考辐射场附加过滤研究

利用MCNP程序对450kV X光机建立模型，计算各参考辐射质下的原级谱、平均光子能量。对ISO 4037-1中已有范围的管电压与附加过滤片厚度关系曲线进行二次拟合，确定模拟厚度范围。利用MCNP计算各附加过滤组合下的平均光子能量值，并与外推法得出的平均光子能量值进行对比，将相对偏差最小时的附加过滤作为最终确定值。从部分能量点的模拟与测量结果的对比可看出，模拟得出的平均光子能量在2%以内与PTB及中国计量院测量结果保持一致，附加过滤最大相对偏差为18%。     

基于蒙特卡罗法的航空发动机空气系统稳态算法优化

针对目前航空发动机空气系统稳态算法中收敛性依赖初值的问题,将蒙特卡罗方法与流体网络法综合应用到空气系统可压缩流体一维网络计算中,提出了一种新的计算方法Monte Carlo-Fluid Network（MC-FN）。该方法将空气系统简化为由节点和元件组成的网络,借助蒙特卡罗方法获得空气系统内各节点压力分配,再根据空气系统中各元件流阻特性和换热特性计算流量、温度。计算中通过将游动次数比较少的蒙特卡罗方法的计算结果作为流量残差法节点压力、温度的初始值,实现快速求得精确收敛解。与流量残差算法相比,MC-FN方法计算精度不变,收敛速度提升了66.5%;与线性求解法相比,MC-FN方法的计算精度提升了25.2%,收敛速度提升了43.8%。     

热力学非平衡状态下DSMC绝热壁边界的抽样方法

从直接模拟蒙特卡罗仿真DSMC方法的实现途径出发,经过分析推导,给出了一种在考虑分子转动能、振动能等热力学非平衡状态下统计求解绝热壁温的DSMC抽样方法。然后采用锥体外形的算例验证,分析了该方法在实现绝热壁条件时的统计误差。计算结果表明,此种方法能够有效地实现DSMC方法壁面处的绝热壁条件,随着采样步数的增加,其边界条件实现时的统计误差不断降低,该方法收敛后的统计误差可降至0.5%,能够满足工程要求。     

基于蒙特卡罗模拟法的航空发动机装配公差分析

讨论了蒙特卡罗法的基本原理,研究了利用蒙特卡罗模拟法进行装配尺寸链公差分析的实现方法,提出了基于蒙特卡罗模拟法的公差分析流程,并以某型航空发动机高压涡轮的一个径向尺寸链为实例进行了计算。计算结果表明：蒙特卡罗模拟法用于公差分析具有可行性,其精度随着模拟次数的增加而提高,在模拟次数达到200000次以上时,其精度已经较高。蒙特卡罗法作为一种统计分析方法,同时适用于线性和非线性尺寸链。