侧向喷流干扰流场建立与消退过程数值模拟

针对高超声速轨控侧向喷流的非定常干扰效应问题,应用非定常数值模拟方法和多重网格加速收敛技术,研究了锥-柱-裙外形轨控侧向喷流干扰流场的建立和消退过程,获得了详细的喷流瞬时干扰流场结构特性,分析了法向力放大系数、干扰力矩系数、法向力系数及俯仰力矩系数随时间的变化特性。研究结果表明：在侧向喷流干扰流场建立和消退过程中气动力变化较大,存在峰谷值;法向力放大系数及干扰力矩系数的定常值和非定常时均值之间存在明显差异。     

大气层内导弹侧向喷流干扰流场的数值模拟

采用3阶的WENO格式计算了大气层内带有侧向喷流的超声速导弹绕流流场。计算显示了喷流和绕流干扰的流动结构。喷流上游的高压区域包络弹体,与喷流下游的低压区域一起削弱了喷流产生的侧向力。在文中计算条件下,喷流放大因子小于1,并形成了较小的俯仰力矩。     

过渡区球头气动力特性DSMC仿真研究

采用DSMC（Direct Simulation Monte-Carlo）方法模拟三维球头过渡区流动。文章给出了球头绕流流场的温度与压力分布图,以及在不同Knudsen数、不同马赫数条件下球头表面压力系数分布曲线,并将阻力系数模拟结果与桥函数公式进行了比较分析。     

数值模拟侧向超声速单喷流干扰流场特性

采用数值方法研究了平板上超/高超声速来流与超声速横向喷流相撞引起的复杂干扰流场特性。所建立的单介质冷喷流数值模拟方法,经过了表面多方位压力分布测量结果、纹影显示的激波结构以及表面油流图谱表现的表面分离范围的实验验证。根据数值模拟与实验对比的结果,合理地描述了喷流干扰流场压力分布以及表面、空间结构特性,并分析了压力比对流场结构和特性的影响。     

大攻角侧向喷流流场特性及喷流干扰效应

采用HLLC近似黎曼解格式求解N-S方程,对大攻角条件下超音速尖拱细长旋成体导弹的单侧向喷流干扰流场进行了数值模拟.并对喷流流场压力分布、喷口附近流动及流场旋涡结构进行了研究,分析了大攻角对细长旋成体侧向喷流干扰压力分布及干扰因子的影响.研究结果表明,HLLC格式计算复杂流场具有较好的数值稳定性;侧向喷流与来流的干扰总...     

轴对称横向喷流强干扰流场的巨型机向量并行计算

本文在国产巨型计算机ＹＨ－１和ＹＨ－２上对ＭａｃＣｏｒｍａｃｋ显格式解轴对称横向喷流强干扰流场的计算程序进行了重构和向量优化使其能并行执行。得到该程序在ＹＨ－１上的向量加速比约为６．０，ＹＨ－２单机向量运行速度约为ＹＨ－１上的向量运行速度的３倍。与传统的串行计算相比较，并行计算大大缩短了程序运行的ＣＰＵ时间。     

火星探测器再入RCS喷流干扰效应数值模拟研究

针对火星大气环境下的反作用力控制系统（RCS）喷流干扰问题,采用自研CFD软件数值模拟“火星科学实验室”（MSL）外形在火星大气环境中的喷流干扰效应;并与文献中的高超声速喷流干扰风洞试验和计算结果进行对比验证;继而完成了火星大气环境高超声速（M_∞=5～10）条件下的偏航方向喷流干扰效应数值模拟,获得了不同攻角（α=0°～-30°）、不同马赫数条件下的喷流干扰气动规律：负攻角增大将导致喷流干扰效应增强以及附加偏航力矩增加;来流马赫数对附加干扰力矩的影响较小。研究结果可为火星探测器喷流控制设计提供参考。     

过渡区DSMC数值模拟的一种自适应当地时间步长新方法

采用优化的DSMC方法,研究了高超声速三维复杂外形飞行器在过渡领域飞行的气动力、热特性。对建立在流场网格和物面网格耦合基础上的三维复杂流场DSMC通用程序进行了变时间步长的优化设计,发展了一种自适应当地时间步长新方法,根据网格内仿真分子的碰撞数,自动调整网格中进行仿真分子运动和碰撞计算的时间步长,有效地提高了程序的计算效率。对比计算了三维复杂飞行器流场压力、热流分布量,飞行器表面气动力、热分布参数,证明了自适应当地时间步长新方法大幅度缩短了流场稳定所需的CPU时间,并且不降低流场模拟结果的精度。     

尾部喷流对飞行器阻力影响的数值模拟分析

底阻在弹类飞行器阻力中占比较大,准确预示底阻对于弹类飞行器飞行性能评估至关重要,而发动机尾部喷流对底阻影响明显。采用基于雷诺平均Navier-Stokes方程的流场仿真方法研究了飞行器底部发动机喷流和外流干扰流场特性,主要分析了喷流对飞行器阻力的影响。飞行器安装了两台推力可调的液体火箭发动机,发动机在不同飞行工况下采用不同推力工作。分别研究了亚音速、跨音速和超音速典型飞行工况下弹体底部无喷流状态、单喷管喷流状态和双喷管喷流状态时,飞行器阻力变化情况。结果表明:不同马赫数下,发动机喷流对底部阻力影响情况基本一致,与无喷流情况相比,当发动机工作时,无论单喷管喷流还是双喷管喷流状态,底部发动机喷流引射效应明显,弹体阻力系数明显增加。     

高超声速飞行器复杂流场过渡区DSMC数值模拟的一种新方案

采用DSMC方法,研究了高超声速三维复杂外形飞行器在过渡领域飞行的气动力、热特性.提出了一种对飞行器物面网格与DSMC计算域网格分别标识,进行高超声速复杂流场过渡区DSMC数值模拟的一种新方案,通过判断模拟分子与表面碰撞来完成飞行器物面网格与DSMC计算域网格间的信息传递和信息存贮,对于复杂外形飞行器的精确描述的物面网格不需做进一步处理,直接应用于不依赖于飞行器外形的DSMC计算的通用子程序中.对物面网格的标识、分子的运动及与表面碰撞的判断、流场量的采集的实施方法和细节进行了分析和讨论.仿真计算了三维复杂飞行器流场压力、热流分布量,飞行器表面气动力、气动力矩和气动热参数,证明了采用的方案和技术的有效性.     

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

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

基于MPI的弹道仿真的并行计算研究

介绍了MPI的基本概念,列举了提高并行计算效率的一些措施,并实现了基于MPI消息传递机制的弹道仿真计算模型的并行编程,最后对并行计算和串行计算的结果进行了比较。研究结果表明,将并行计算技术应用于弹道仿真中是可行的,并获得了较高的加速比和计算效率;与单机运行相比,计算机机群可以大大缩短计算时间,提高计算效率,而且可以扩大计算规模,系统适应性和延展性较好。     

用动网格方法模拟导弹发射过程中的燃气射流流场

采用局部非结构网格弹性变形方法和网格再生方法相结合的动网格技术，计算导弹与发射简具有相对运动的燃气射流非定常流场。计算中将安置在发射简内的导弹作为运动实体，随着导弹运动，相应流场计算域边界发生变化。计算结果表明，采用动网格技术模拟导弹发射过程中燃气射流的非定常流场，具有较高的精度。     

导弹轨控喷流气动干扰全空域插值方法研究

导弹侧向喷流气动干扰受海拔高度影响显著,主要原因在于雷诺数和压力比会随海拔高度的变化而变化。受大气参数的影响,雷诺数和压力比的变化规律并不一致。为能在现有气动模型基础上既保证插值精度又减少节点数,通过数值模拟方法研究了雷诺数和压力比对轨控喷流气动干扰的影响。结果表明:海拔高度对喷流气动干扰的影响主要由压力比变化造成,雷诺数变化产生的影响相对较小,且喷流干扰与压力比间存在一定的线性关系。因此,提出了一种基于压力比的喷流气动干扰插值方法。相比于传统的基于高度的插值方法,该方法可显著提高插值精度,并减少节点数,适用于全空域轨控喷流气动建模。     

超声速气流中液体横向射流的气液相互作用过程数值研究

采用Euler-Lagrange方法对来流马赫数为1.94的超声速气流中液体横向射流的气液相互作用过程进行数值研究。计算给出的射流穿透深度、液滴Sauter平均直径(SMD)及液滴速度分布均与实验吻合较好。仿真结果较详细地揭示液体射流喷雾与气流之间的强烈相互作用过程。受液雾影响,射流前形成较强激波,气流依次经过激波及液雾区域,气流速度存在两次下降过程。计算结果揭示,超声速来流可以与射流的液滴轨迹相交,气流经液雾前沿进入液雾区域后,流向往壁面偏折。本文首次发现并提出,由于气液相互作用诱导形成两组反向反转漩涡对,这对于理解两相混合过程具有重要意义。气液相对滑移速度的分析表明,液滴在穿透自由来流并开始转向时受气流作用最为显著,完成转向后气液相互作用逐渐减弱。     

单喷管火箭自由喷流噪声数值模拟

利用非定常燃气流场数值模拟方法,研究超声速喷流核心区及其外围亚声速区燃气流动态分布和流动特性,在此基础上利用欧拉方程条件的伽辽金有限元方法以及FWH方法,实现并完成喷流噪声传播特性、辐射特性数值模拟。喷流噪声数值模拟结果显示：燃气流推进初期,强喷流噪声区域紧随燃气流前锋;燃气流场相对稳定后,强喷流噪声区域主要位于燃气流等能区末稍,一些小尺度试验的燃气流激波系附近也存在较强喷流噪声。这些强喷流噪声主要由燃气流前锋带动的大涡动态卷吸、燃气流强湍流脉动以及激波扰动引起。受数值模拟网格分辨率影响,当前仅能保持中低频段声压级数值模拟结果与实测结果总体接近。     

SPMD并行计算机系统中互连结构的模拟研究与实现

并行计算机系统中网络对整个系统性能影响起到了关键作用，而SPMD并行计算机系统在科学计算，特别是航空航天中比较普及，因此开发了通用SPMD并行系统结构性能评价模拟系统。该文描述的面向对象网络模拟系统是该模拟环境的一部分，采用伪汇编代码的算法程序作为输入，通过配置结构参数、与系统提供的描述不同结构的伪虻编指令函数库编译链接，达到对不同层次细节、不同侧重部件模拟研究的目的。与同类研究相比，本文的研究具有更通用的特点，能够方便地构造不同的动态和静态网络，特别是能够对可重构的多级互连网络进行性能评价研究。     

羽流中固体颗粒在真空环境下分布的数值仿真

通过对气?固两相间动量和能量相互作用解耦处理,建立了一种适于模拟真空环境气固两相混合物羽流的DSMC双向耦合算法与固体颗粒空间输运变化特性的TPMC计算技术。仿真了固体火箭发动机两相羽流流场和固体颗粒在远离发动机喷口数十公里空间扩散运动分布特性,通过将计算结果与典型文献结果及理论分析比较确认,证实本文方法的准确可靠性。结果表明,固体颗粒对气相的扩散有一定的阻滞作用;仅在离发动机喷口一定距离,气相对固体颗粒有较大影响,致颗粒温度下降、速度增加;颗粒温度随发动机喷口距离增大而减小,一直要在远离喷口上百公里颗粒温度才随轴向位置趋于平衡,且不同尺寸的颗粒温度差别较大,对指导外层空间高真空环境气固两相羽流传输影响工程研制具有重要意义。