飞机动态推出控制策略仿真优化模型

机场离港运行的无序化造成滑行过程中的长时间排队等待及大量燃油浪费。为了减少燃油消耗和废气排放,在已有的动态推出控制策略基础上提出了阶梯函数控制策略(Step function,SPC)和非线性函数动态推出策略(Nonlinear function,NPC)的一般形式,以离港成本为目标,建立了基于停机位等待惩罚的动态推出控制模型,在不延误的前提下提出了一种基于网格参数优化的蒙特卡洛仿真优化算法。通过北京首都机场实际运行数据对推出过程进行仿真计算,并与无控制策略以及传统N-control策略的最优可达解进行仿真对比,结果表明:在不延误的前提下,提出的推出策略可以更加有效地降低平均滑行道滑行时间,NPC策略的离港运行成本和燃油成本可降低45.52%和54.23%,SPC虽然成本节省劣于NPC策略,但是其简单的操作方式可以为离港推出调度方式的改进提供决策支持。     

结合跳点引导的无人机随机搜索避撞决策方法

针对无人机在城市空域环境和密集交通流下的避撞决策问题，提出马尔科夫决策过程（MDP）和蒙特卡洛树搜索（MCTS）算法对该问题进行建模求解。蒙特卡洛树搜索算法在求解过程中为保证实时性而使其搜索深度受限，容易陷入局部最优，导致在含有静态障碍的场景中无法实现避撞的同时保证全局航迹最优。因此结合跳点搜索算法在全局规划上的优势，建立离散路径点引导无人机并改进奖励函数来权衡飞行路线，在进行动态避撞的同时实现对静态障碍的全局避撞。经过多个实验场景仿真，其结果表明改进后的算法均能在不同场景中获得更好的性能表现。特别是在凹形限飞区空域仿真模型中，改进后的算法相对于原始的蒙特卡洛树搜索算法，其冲突概率降低了36%并且飞行时间缩短47.8%。     

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

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

多航天器协同探测效能影响参数不确定性分析

为解决传统基于蒙特卡洛仿真的探测效能不确定性分析评价方法中时间效率低，参数关系映射单一的问题，提出一种基于神经网络的不确定性分析方法，利用人工神经网络在拟合回归分析上的非线性特性，设计了能够替代复杂系统的神经网络结构，能够通过少量仿真计算结果作为训练样本实现模型的收敛并能有效反映被替代系统的原有特性。选择以一个多航天器的天文观测任务作为典型用例，建立了仿真分析模型，并与蒙特卡洛仿真方法结果进行了对比，验证了此方法的准确性和计算效率。     

Ni,Ti及NiTi合金低能电子散射的Monte Carlo模拟

采用MonteCarlo方法,模拟了低能电子束(能量E0≤5keV)作用下Ni,Ti及NiTi合金固体中的电子散射,分析了3种金属/合金中散射电子的能量与空间分布。研究表明:(1)入射电子能量E0越小,各材料的背散射电子(BSEs)深度、吸收电子(AEs)深度及散射电子能量损失的深度均越小,分辨率越高;(2)相同E0下,各材料的BSEs,AEs及能量损失深度均有NiNiTi>Ti;(3)NiTi合金的BSEs深度、AEs深度及能量损失的深度和分辨率分别近似于Ni,Ti固体相应量的平均值。     

电推进羽流粒子束溅射的蒙特卡罗模拟

为了研究电推进羽流对推力器及航天器造成的溅射腐蚀作用,采用蒙特卡罗方法对羽流等离子体对材料的溅射进行模拟,研究了溅射的级联碰撞机理,并开展了Xe离子正入射、斜入射Cu材料的溅射率、溅射产物的能量分布以及空间角分布的研究.结果显示,蒙特卡罗方法对溅射率和能量分布的计算结果与试验结果吻合较好,误差<30％,但是斜入射溅射率...     

稳态等离子体推进器羽流场数值模拟

采用二维轴对称模型,使用粒子网格法(PIC)和直接模拟蒙特卡洛法(DSMC)相结合的方法,对稳态等离子体推进器(SPT)羽流场进行了数值模拟.采用DSMC方法中的随机取样频率法(RSF)求解粒子碰撞过程,并对比了不同的分配电荷方式、电子运动模型及SPT出口条件时的羽流场.将不同条件计算得到的羽流场中距SPT出口0.1?m,0.5?m及1.0?m处的离子电流密度和电荷密度与实验结果进行了对比,得出在采用面积权重法分配电荷、等熵模型描述电子运动和用实验值设定发动机出口参数时对SPT羽流场数值模拟的电流密度和轴向附近的电荷密度结果与实验结果符合程度较好的结论.     

考虑分布参数依赖的风险混合不确定性分析方法

在定量化风险评估中,针对变量分布参数存在依赖时的混合不确定性传播问题,提出一种考虑分布参数完全依赖、部分依赖以及独立传播情形的双层混合不确定性刻画与传播框架,内外层分布参数不确定性分别用概率分布与可能性分布刻画,采用蒙特卡罗模拟法与模糊扩展原则相结合的数值求解方法。针对认知不确定性分布参数依赖性,构建了统一的认知不确定性分布参数依赖性模型,并给出依赖性系数的概念。为实现认知不确定性分布参数独立性采样,设计了基于D-S证据理论与随机集相结合的不确定性传播算法,相比于概率刻画下的双层蒙特卡罗方法,计算代价有效降低。以某型氢氧发动机贮箱共底漏气率为算例,验证本文方法的有效性与可行性。     

基于蒙特卡洛方法的冰晶撞击特性计算

航空发动机摄入冰晶易发生内部结冰、造成发动机推力损失,严重时引发航空事故。冰晶结冰数值模拟是研究冰晶结冰的重要手段,而撞击收集系数计算是冰晶结冰数值模拟的关键环节。基于NNWICE平台开展冰晶结冰前期的阻力系数、运动轨迹和撞击收集系数的数值模拟与计算研究,实现了基于蒙特卡洛方法的冰晶撞击收集系数的准确计算;在此基础上发展了拉格朗日框架下冰晶粒子运动、传热传质耦合数值计算方法;通过完成二维NACA0012翼型的相关算例,对比分析了不同阻力模型和冰晶相变对撞击收集系数的影响。结果表明阻力模型对冰晶撞击收集系数影响较小,冰晶相变过程会对冰晶运动轨迹及撞击收集系数产生显著影响,为进一步发展冰晶结冰数值模拟计算方法奠定了基础。     

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

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