高超声速乘波飞行器三维流场的并行数值模拟

针对高超声速乘波飞行器三维绕流流场,在基于LINUX+MPI系统的分布式并行计算平台上,并行求解了三维雷诺平均的N-S方程。并行数值方法采用的是有限体积方法(FVM)、OC-TVD差分格式、B-L代数湍流模型及流场分区的并行方法。计算结果表明,所采用的并行数值模拟方法能够求解包含强激波的流场,激波穿越区域边界时无断层、错位等通量不守恒的现象。并行计算效率高,8个处理机计算时的并行加速比达到了6 8。     

基于TOPSIS的战区高层反导威胁评估

在分析战区高层反导威胁评估特点的基础上,建立了以来袭TBM发射点、预测落点、射程等为核心的战区高层反导威胁评估模型指标体系并进行了相应的量化,然后利用基于熵值权重确定的TOPSIS理论对模型进行求解。通过实例证明所提出的战区高层反导威胁评估模型和算法的有效性,对研究美军的战区高层反导武器系统具有一定的借鉴意义。     

弹丸爆炸仿真与破片质量统计软件开发

针对利用试验研究弹丸爆炸过程难度大的问题,采用ANSYS／LS-DYNA有限元软件对弹丸爆炸过程进行了仿真,设计开发了破片质量分布统计软件,再现了弹丸爆炸的宏观过程,并且得到了破片质量分布的定量结果．     

翼栅绕流流场数值模拟及动力性能研究

基于2D不可压缩定常流动的控制方程,对5个NACA4412非对称翼型组成的平面直列翼栅绕流流场进行了数值模拟,对比分析了同一翼型在翼栅中工作和单独绕流时的升、阻力情况,验证了翼栅装置工作的可行性,反映了翼栅绕流的客观规律性。结果表明:翼栅翼型头部流速和压力急剧变化,尾部出现流涡尾迹,前部是负压区,后部是正压区,内部流速大,压力小;在1/3弦长附近翼栅翼型表面压力出现拐点,拐点之前翼型上下表面压力差较大,拐点之后翼型上下表面压力差较小;翼栅内部流场具有相似、相近性,外部流场差异性较大;翼型的压力差是产生升力的主要原因。     

部分服务台同步多重休假的GI/M/c稳态队长

在经典GI/M/c排队中引入部分服务台同步多重休假策略,利用拟单生过程和矩阵几何解的方法,求解系统的稳态队长分布及其条件随机分解。     

装备体系结构优化可行方案空间构建方法

针对装备体系结构方案数量规模较大及传统构建方法效率较低的问题,以单一费用约束条件下的装备体系结构优化问题为例,提出了一种可行空间构建的新方法,试图缩小体系成员的遍历范围,减少运算量,提高体系结构可行方案空间的构建效率。最后,利用Matlab软件对新算法进行数值模拟与仿真,并与传统方法进行了比较,结果表明:新方法的计算效率得到了明显提高。     

复杂条件下岩质边坡施工过程数值模拟

为研究基坑开挖时的变形、稳定性及对坡顶建筑物的影响,以实际工程为例,运用数值模拟分析方法,对含有软弱夹层的顺层岩质基坑边坡在直立开挖与支护过程中的位移场、应力场及对坡顶建筑物影响进行了动态模拟分析。研究结果表明:开挖过程中,岩体卸荷引起变形,支护结构可以限制变形的发展,开挖面处岩体失去原有限制发生卸荷回弹,使得边坡顶部和靠近开挖面处出现拉应力集中区,层状岩体的软弱夹层处会产生较大的剪应力,并随着开挖加深而增大;开挖后边坡顶部已建挡墙附近应力消散,坡顶产生不均匀沉降。     

纯引力轨道验证质量辐射计效应与气体阻尼耦合模型的数值分析

辐射计效应和气体阻尼是纯引力轨道验证质量的重要干扰力,是影响纯引力轨道构造水平的重要因素.在纯引力轨道飞行器中,这两种力分别描述了由腔体中温度梯度和验证质量相对运动引起的气体分子作用,两者从不同角度描述了气体分子作用力,均是气体分子作用力的一部分,而两者的耦合模型则可以反映验证质量受到的气体分子作用力总和.针对耦合模型形式复杂的特点,本文以内编队系统为例,利用数值方法分析了耦合模型中的影响因素,这些因素包括内卫星相对运动速度、内卫星半径、外卫星腔体半径、腔体平均温度、腔体温差和腔体平均压力等.对大量计算结果进行了数据拟合,给出了内卫星气体分子作用力与各物理参数关系的拟合公式,和原始计算结果相比,拟合误差在20％以内.     

液氧与气氧对空气加热器燃烧流场的影响分析

针对一种基于液体火箭发动机燃烧室结构的空气加热器,采用数值仿真技术研究了加热器内部喷雾燃烧、燃气掺混以及出口流场分布等参数.分析对比了采用酒精/液氧/空气与酒精/气氧/空气两种不同氧化剂物态三组元同轴直流式喷嘴所得到的燃烧流场的区别,并通过改变燃烧室特征长度,分析了两种计算工况的加热器的性能差异.结果表明,喷入氧化剂的物态对燃烧流场影响较大,采用液氧喷嘴的火焰较长,气氧喷嘴的火焰分布较宽,且相对于液氧喷嘴,气氧喷嘴的燃烧室前端回流区由于掺混较多的燃气,导致喷注面板附近燃气温度较高,面板承热压力较大.设计的加热器均可保证两种喷嘴的出口流场品质较高,在保证流场出口品质的原则上,气氧喷嘴的燃烧室特征长度可至少小于液氧喷嘴的1/4.     

非定常Euler方程数值计算中高精度格式应用

采用五阶精度加权紧致非线性格式(WCNS)和非定常“双时间步”方法求解非定常Euler方程,模拟NACA0012翼型强迫俯仰振动流场,研究了高精度格式应用到非定常计算时“双时间步”方法物理时间步长、子迭代收敛判据、子迭代步数以及物理时间导数离散方法对计算精度和计算效率的影响.