二维过渡区含化学反应的喷流干扰流场DSMC/EPSM仿真

构造可仿真含化学反应的稀薄／连续不同流动区域的DSMC／EPSM混合算法，应用该算法，数值模拟过渡区有横向喷流干扰的二维有限平板高超声速化学反应流场，分析流场结构和物面气动系数分布等特性，充分反映稀薄／连续混合区域高超声速流动中的化学反应效应，同时通过与DSMC算法比较，论证DSMC／EPSM混合算法的可行性与计算效率相对提高的优越性。     

高超声速稀薄流的气粒多相流动DSMC算法建模研究

基于直接模拟Monte Carlo(DSMC)方法,构造适用于DSMC算法的固态和液态颗粒碰撞、聚合和分离模型,发展稀薄条件下双向耦合作用的气粒多相流的DSMC算法,在此基础上初步实现高超声速稀薄流环境中的气粒多相喷流流场数值模拟.算例结果表明该方法能为稀薄过渡区气粒多相流动提供一种新的应用研究手段.     

稀薄流过渡区气固两相喷流的建模与数值模拟

从描述粒子运动的细观层次出发,运用DSMC方法,构造了稀薄流过渡区气固相互作用模型和颗粒碰撞模型,数值模拟了有气固两相横向喷流干扰的二维有限平板高超声速流场,分析了流场结构和物面气动系数分布等特性.结果表明,对于不同的来流Knudsen数和颗粒质量载荷等参数,喷流干扰流场存在差异,同时反映了基于DSMC方法的气固相互作用模型是研究稀薄及过渡条件下的气固两相流动的有效方法.     

稀薄过渡流区横向喷流干扰效应数值模拟研究

为了准确预测稀薄过渡流区横向喷流与稀薄大气的干扰流动特征,本文在建立直角与表面非结构网格混合结构的DSMC数值算法以及碰撞网格自适应算法基础上,提出基于MPI的静态随机负载平衡技术,构建了DSMC并行计算代码。计算的不同压力比条件下三维平板模型横向喷流与稀薄大气干扰的分离长度与低密度风洞试验有较好的一致性,验证了本文算法的可靠性。开展了细长钝双锥外形高超声速稀薄来流与超声速/高超声速横向喷流干扰效应的研究,计算分析了不同飞行高度、不同飞行速度、不同飞行攻角、不同喷流推力下复杂流场结构和对气动力特性的影响规律。考察了RCS喷管出口参数不同(均匀/非均匀)对喷口附近分离涡和分离长度的影响。     

侧向喷流与飞行器三维绕流干扰流场特性的DSMC/EPSM数值模拟

采用DSMC/EPSM混合算法,研究了带有侧向喷流的复杂外形高超声速飞行器稀薄区和连续区三维混合流场的干扰特性。采用一种飞行器物面网格与DSMC计算域网格分别标识的方法,通过判断模拟分子与表面碰撞来完成飞行器物面网格与DSMC计算域网格间的信息传递和信息存贮,对于复杂外形飞行器的精确描述的物面网格不需做进一步处理,直接应用于不依赖于飞行器外形的DSMC计算的通用子程序中。引入碰撞数做为混合参数,对流动区域进行划分,并通过网格中碰撞数的统计和处理,有效地将DSMC和EPSM方法结合在一起。仿真计算了三维复杂外形飞行器带喷流流场压力、热流分布量,飞行器表面气动力、气动力矩和气动热参数,对喷流与高超声速空气绕流相互作用以及它们与物面之间的相互扰动进行了分析,证明了采用的方案和技术的有效性。     

稀薄流高超声速钝化前缘逆向射流流场特征研究

为准确预测稀薄过渡区逆向射流对稀薄大气的干扰流动特征，本文采用直接蒙特卡洛算法（DSMC），在自由流马赫数为7的情况下，对稀薄流中逆向射流干扰下的高超声速流动进行了数值模拟。计算中考虑了60～90 km高度的典型大气环境，研究了压比和自由流克努森数对压力、剪切力和热流密度的影响规律。结果表明，逆向射流改变了流场结构，且具有明显的热防护性能，但减阻效果不明显。逆向射流对高超声速稀薄来流的影响随压比的增加明显增强，形成的马赫盘大小也随之增加。高度的增加伴随稀薄效应不断增强，自由流克努森数增大，导致弓形激波厚度随高度增加而增加，壁面处的回流区范围不断减小至消失。这项研究初步揭示了稀薄流中逆向射流/高超声速流动相互作用的机制。     

高超声速连续/稀薄流自适应界面推进重叠网格方法及应用研究

根据高超声速飞行器在过渡流域飞行时绕流流场具有连续流与稀薄流并存的特点,提出了一种适用于混合流场自适应界面推进的多尺度重叠网格方法。用两套独立生成不同尺度的结构网格和非结构网格覆盖整个流场。根据连续/稀薄流界面失效函数的截止值在结构网格内进行自适应"挖洞"。在挖洞后的稀薄流区利用非结构网格DSMC方法进行计算;在连续流区的结构网格内采用基于N-S方程的CFD算法。两种算法进行耦合计算和求解,在连续/稀薄流分界面进行流场信息的插值、交换和处理。根据流场实时耦合计算结果,利用界面失效函数自适应改变两套网格的重叠部分。最后,通过对超声速圆柱和高超声速带扩张角圆管绕流的数值模拟,验证了方法的可行性和有效性。     

高超声速飞行器DSMC／EPSM自适应当地时间步长混合算法

采用优化的DSMC/EPSM混合算法,研究了高超声速飞行器三维复杂流场在近连续区的气动力热特性.发展了一种DSMC/EPSM自适应当地时间步长混合算法,根据网格内仿真分子的碰撞数,自动调整网格中进行仿真分子运动和碰撞计算的时间步长,有效地提高了程序的计算效率和计算精度.用DSMC方法、DSMC/EPSM混合算法和DSMC/EPSM自适应当地时间步长混合算法对比计算了三维复杂飞行器流场压力、热流分布量,飞行器表面气动力、热分布参数,证明了DSMC/EPSM自适应当地时间步长混合算法大幅度缩短了流场稳定所需的CPU时间,并且不降低流场模拟结果的精度.     

超声速和高超声速横向喷流数值模拟

为了解横向喷流的干扰影响和喷口附近的流场结构,采用分块对接网格和“0”型网格技术,数值求解N—S方程来模拟超声速和高超声速流场中横向喷流的干扰流场。对两种尖拱弹身外形的超声速和高超声速喷流干扰流场进行了数值计算,计算结果与风洞实验数据吻合一致。在此基础上,开展了某型导弹多喷构型干扰流场的数值模拟,得出一些横向喷流数值模拟的结论。     

高超声速化学非平衡流动Navier-Stokes/DSMC耦合算法

基于相同的化学反应模型,在已有计算流体力学（CFD）和直接仿真蒙特卡罗（DSMC）方法及程序的基础上,采用Modular Particle-Continuum（MPC）耦合技术,建立了包含化学非平衡Navier-Stokes/DSMC耦合算法。算法结构中DSMC计算区域在CFD计算结果上根据当地克努森数自动选取。发展了适用于流场分区信息交换的亚松弛技术,抑制DSMC方法对CFD计算的影响。把DSMC方法和CFD的应用范围拓展到过渡流区,为复杂飞行器近连续过渡流区高超声速化学非平衡流动数值模拟研究提供了一种工程适用的预测分析手段。通过对二维圆柱高超声速化学非平衡绕流的算例与其他结果的比较研究,表明耦合算法不论在流场结构、流场非平衡现象,还是飞行器表面参数、整体气动力/热特性方面,都能够得到与全DSMC计算吻合的结果,证实了所建立的Navier-Stokes/DSMC耦合计算模型与方法的有效性和可靠性。仿真了某航天器解体碎片在过渡区的化学非平衡流动,得到碎片在过渡区的气动力/热特性,为碎片的陨落计算提供依据。     

超声速平头圆柱绕流DSMC/EPSM混合算法研究

本文将EPSM算法与DSMC方法结合，构造了可模拟含近连续流区及过渡流区的DSMC／EPSM混合算法。运用混合算法模拟马赫数等于5．37时超声速平头圆柱绕流，并与DSMC结果进行比较，证明了DSMC／EPSM混合算法的有效性，并对两种算法的计算效率进行了比较。     

稀薄大气下机械展开式再入飞行器气动特性研究

以机械展开式再入飞行器为研究对象,采用DSMC(直接蒙特卡洛模拟)方法计算分析了其在稀薄大气环境下的气动特性.讨论了飞行器在不同高度下的流场特性和升阻及配平特性,结果表明飞行器流场及受力特性对大气稀薄程度显示出较强的敏感性;进一步分析了连续流过渡至稀薄大气下飞行器的阻力系数,得出了随着来流Knudsen数增加,大气稀薄效应逐渐显现,传统的CFD方法已不再适用于描述流动特征的结论,同时给出了依据来流Knudsen数大小采取CFD和DSMC方法分别进行计算分析的建议.     

高超声速平头圆柱绕流化学反应流场的DSMC/EPSM混合算法研究

本文发展了平衡粒子模拟方法(EPSM)，建立了与高温气体化学反应动力学理论相匹配的EPSM耦合模型，并通过混合参数进行流区的自动识别，将EPSM方法与蒙特卡罗直接模拟方法(DSMC)结合，构造了可模拟化学反应流动的DSMC／EPSM混合算法。应用该算法对汲及化学反应的轴对称情况下高超声速平头圆柱绕流流场进行模拟，将结果与DSMC方法的结果进行比较，验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较，发现混合算法能够大大提高计算效率。     

球双锥外形稀薄过渡流区气动特性研究

采用风洞试验和工程计算相结合的方法,对球双锥类外形的稀薄过渡流区气动特性进行了研究。风洞试验在高超声速低密度风洞中进行,设计了高灵敏度的外式微量天平,解决了温度防护和高精度校准的问题,得到了球双锥外形在高度65km~75km、高超声速流动的气动力数据。在此基础上,结合球双锥细长体外形绕流的特点和DSMC数值模拟结果,从连续流粘性干扰模型、关联参数选取、过渡流搭桥函数等三个方面对现有工程计算方法进行了改进。研究结果表明:稀薄效应对球双锥外形的气动特性,尤其是轴向力系数和升阻比产生非常大的影响;计算与试验结果的对比表明,计算模型的改进效果显著,大幅度提高了预测值的准度。     

电离对高超声速稀薄流飞行器气动热影响

将电离反应模型扩展到（direct simulation Monte Carlo,DSMC）方法中,研究了电离反应效应对高超声速稀薄流飞行器气动热的影响特性.针对稀薄流场中电子出现带来的实际困难,引入“捆绑法”思想处理电子在流场中的运动,并给出了电离反应模型及电离反应处理方法.在以RAM-C Ⅱ飞行器外形为例对增加了电离反应的DSMC代码进行验证的基础上,以“星尘号”探测器外形为研究对象,针对不同飞行高度下5组元混合气体模型（无电离）和11组元混合气体模型（含电离）的化学非平衡流动开展了数值模拟,细致分析和对比了电离反应效应对探测器气动热的影响规律.研究结果表明:采用的电离反应处理方案能够模拟带电离反应的高超声速化学非平衡稀薄流动.在飞行高度为60km时电离反应对探测器气动热的影响最为强烈,使探测器的驻点热流密度降低了5.12%,电离反应对探测器气动热的影响随气体稀薄程度增加而减弱.      

激光光船周围连续/稀薄混合流场的数值模拟

当火箭模式激光推进器在高空稀薄大气中飞行时,需要在激光聚焦区附近喷射稠密气体工质,从而利用激光聚焦击穿气体工质形成等离子体,为光船提供推力。为了分析该光船的推力产生过程,采用连续流体力学计算方法(CFD)和直接数值模拟蒙特卡洛方法(DSMC)耦合的混合算法,数值模拟飞行高度为80 km飞行马赫数为5的激光光船周围高温连续气体和稀薄环境大气耦合的非定常混合流场,并得到推力随时间的演化曲线。结果表明,稀薄效应明显增加了激波厚度,高温效应延迟了推力的产生,并导致冲量耦合系数和比冲的降低。