高速稀薄流域防热减阻概念构型气动性能

为探究高速稀薄流流域凹腔槽道的气动特性,采用直接模拟蒙特卡罗(DSMC)方法,建立了凹腔深宽比为1,槽道高度分别为0、10、20、30、40 mm以及基准高度为20 mm的唇口钝化半径为2～20 mm的凹腔槽道构型,获得了不同算例的气动热系数和气动力系数的变化情况,分析了不同槽道高度以及不同唇口钝化半径的凹腔槽道构型对...     

Markov随机跳变系统的传感器系统误差估计

 针对Markov随机跳变系统的系统误差估计问题,提出一种基于马尔可夫链蒙特卡罗(MCMC)和最大似然估计相结合的在线系统误差估计方法。利用最大似然估计给出系统误差等效后验概率分布函数,采用Metropolis-Hastings抽样方法从该概率分布函数中进行抽样;利用系统误差估计和状态估计互为因果的关系,采用期望极大化(EM)方法迭代估计出最优的系统误差;分别对时变和时不变系统误差场景进行仿真分析,结果表明,在考虑系统误差统计特性的同时,所提方法对解决目标运动模型难以建立情况下的系统误差估计问题具有可行性和有效性。     

不确定性颤振风险定量分析

 针对小风险颤振定量分析,提出一种概率鲁棒颤振初步分析方法。该方法结合鲁棒颤振分析和概率颤振分析,针对不同层次风险和不确定性水平提供颤振边界特性全面信息,给决策制定提供更好参考。在鲁棒颤振分析中,采用振型迭代的方法,通过结构奇异值( μ )分析和特征值求解得到零风险的鲁棒颤振临界速度。在概率颤振分析中,采用标准蒙特卡罗模拟(MCS)方法得到较大风险时的颤振速度分布和概率临界速度。采用概率比较和二分法估计概率鲁棒颤振不确定度半径和概率鲁棒颤振临界速度,同时采用样本重新利用的随机方法估计颤振风险和不确定性定量关系。通过一个考虑集中质量不确定性的大展弦比双梁式机翼的颤振计算分析表明,确定型μ方法只能用于零风险颤振分析,标准的MCS方法只能用于较大风险颤振速度分布问题。概率鲁棒颤振分析结果表明增加1%颤振风险可以使不确定度半径增大62%或颤振速度边界增大5%,从而放宽了设计要求。     

高温高速稀薄流的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方法的边界采样与统计技术

研究了适合在三维直接模拟蒙特卡罗(DSMC,Direct Simulation Monte Carlo)方法计算中根据入口边界条件合理地布置进入计算区域内粒子参数的边界采样技术.提出了一种根据穿越入口边界、出口边界以及撞击壁面边界的粒子统计数据获取边界位置处宏观流场参数的新型边界统计技术,弥补了DSMC计算中传统统计技术仅能够获取流场内部参数,不能够获取界面处流场参数的不足.这种界面统计技术具有普遍的适用性,不仅限于DSMC方法.以无穷远来流圆球绕流为例进行计算,验证了DSMC方法的边界采样与统计技术的有效性.     

用蒙特卡罗方法研究质子在航天器内部充电中的作用

介绍了航天器内部充电的基本物理机制, 重点研究了质子在内部充电中的作用. 用蒙特卡罗方法模拟质子在介质中的输运过程, 计算了简化的平板介质在一定通量的粒子环境中的内部充电情况. 结果表明, 介质内部最大电场与入射质子能量有关, 当质子能量达到14 MeV时候, 内部电场最强; 当质子与电子的入射数目相同, 并且材料参数一样时, 质子产生的最大电场大于电子产生的最大电场. 选取2004年7月26日TC-2卫星姿控分系统故障前的质子和电子通量数据, 分别计算了二者可能引起的内部最大电场. 计算结果表明,质子产生的最大电场比电子产生的最大电场小2~4个数量级, 并且远小于击穿电场; 在某些极端情况下, 例如质子产生的电场和电子产生的电场方向一致的时候, 电场的叠加会使局地电场得到加强.      

地球与大气反射光对空间光学传感器探测能力的影响仿真研究

现有辐射计算模型不考虑地表的各向异性特征,导致计算结果存在较大误差。针对这一问题,提出一种新的地球与大气反射光对传感器照度贡献的计算模型,考虑了地表反射、大气散射和辐射传输三方面的影响因素;在此基础上提出一种启发式蒙特卡罗重点光线采样方法,该方法根据给定的传感器参数和空间几何参数高效解算地球与大气反射光对空间光学传感器探测能力的影响程度。     

基于置信分布的系统可靠度评估蒙特卡罗方法

大型系统、复杂系统的可靠度的求解问题是可靠性工程中一个热点和难点的问题,针对这个问题,提出了一种基于蒙特卡罗模拟仿真的系统可靠度求解方法.通过利用系统的各个单元产品的可靠度置信分布和系统可靠度函数进行蒙特卡罗仿真,评估出系统可靠度.本方法克服了以往的蒙特卡罗仿真方法处理离散分布时抽取时间样本的困难,并且对系统的结构和复杂程度没有要求.最后,给出计算实例对比说明了该方法的正确性和可用性.      

真空羽流场的DSMC并行数值模拟

对DSMC方法进行小喷管真空羽流数值模拟的计算软件进行了并行化改造。计算表明在多处理器(节点)上实现并行计算能够大幅度地提高运算速度,缩短运行时间,增加计算规模,计算结果与串行计算一致。并对目前串行计算不能完成的真空羽流计算的算例成功地进行了DSMC并行计算,计算结果与理论分析一致。      

某异形卷弧翼弹的蒙特卡罗计算机模拟打靶

卷弧翼弹在实际的训练或作战过程中,由于系统干扰或误差的影响,会存在弹道偏差。蒙特卡罗方法(MCM)利用计算机产生的随机数模拟实际系统中的随机变量,通过随机变量参与的系统的大量仿真,可以得到系统的数学期望意义下的解。通过蒙特卡罗计算机模拟打靶,得到了卷弧翼弹的弹着点散布及命中概率,从结果看,弹着点符合正态分布,其命中概率基本满足实际要求。蒙特卡罗计算机模拟打靶方法对于指导训练或作战有一定的指导意义。