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

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

过渡流区N-S/DSMC耦合计算研究

在已有CFD和DSMC方法程序研究基础上,采用MPC耦合处理原理,发展适于流场分区信息交换的亚松弛技术,建立了N-S/DSMC耦合算法.通过对近连续过渡流区不同Kn数条件下钝锥绕流计算与低密度风洞试验结果比较分析,一方面证明了本文所建立的N-S/DSMC耦合计算模型和方法在低Kn数过渡流区的有效性,另一方面剖析了N-S/DSMC耦合计算误差,可望发展一个工程适用的近连续过渡流区高超声速飞行器气动特性预测分析工具.     

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

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

二维过渡区喷流干扰流场的DSMC/EPS仿真

应用DSMC／EPS混合算法，数值模拟过渡区有RCS横向喷流干扰的二维有限平板高超声速流场，分析了流场结构、分离形态及物面气动系数分布等特性。结果表明，不同：Knudsen数下，喷流干扰流场存在差异，过渡区中的稀薄气体效应得以体现，证明了DSMC／EPS混合算法是研究稀薄，连续混合区域高超声速流动的有效方法。     

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

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

高超声速飞行器高温流场数值模拟面临的问题

随着高超声速飞行器目标光辐射和电磁散射特性研究的发展和深入,高温流场特性日益引起人们的关注。由于高温流场特性研究中涉及到非常多的复杂气动现象,如气动加热、烧蚀、辐射、燃烧、化学反应以及湍流等,因此其数值模拟面临着诸多挑战。这里基于连续流计算流体力学(CFD)技术和稀薄气体蒙特卡罗直接仿真(DSMC)方法,从化学物理模型建模、方法稳定性与数值求解效率出发,分析了高超声速飞行器外部绕流、尾迹和发动机喷焰三方面的流场特性数值模拟在不同弹道、热防护手段和飞行流域环境下所面临的问题。在此基础上提出了数值求解技术和化学物理模型建模今后需要发展的方向,为有效提高高超声速高温流场特性数值模拟效率、增加流场特性预测精度提供了指导,从而为研究流场对高超声速飞行器目标光辐射和电磁散射特性影响提供有效的基础数据。     

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

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

高超声速流存在局部稀薄效应的一个判据及相应的流动特性

对于近空间高超声速飞行器的研制,计算流体力学(CFD)起着非常重要的作用。但若流场中存在必须考虑气体稀薄效应的地方,传统的CFD就要做相应的改变,这时首先遇到的问题将是判断是否需要考虑气体的稀薄效应的判据应该是什么?[1]其次就是气体稀薄效应的影响表现在什么地方?如何在CFD中考虑这一效应?本文选取具有代表性的高超声速剪切流为研究对象,采用直接模拟Monte Carlo(DSMC)方法,对剪切强度和稀薄程度不断增强的流动,研究了分子运动速度分布函数以及剪切力的相应变化规律;找到了一个可以判别气体稀薄效应程度的无量纲参数Zh,以及传统连续介质模型下的剪切力和由DSMC所得剪切力随参数Zh的增大而出现的有规律性的差别。     

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

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

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

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

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

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

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

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

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

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

阿波罗指令舱稀薄气体动力学特征的蒙特卡罗数值模拟

本文对阿波罗指令舱稀薄气体动力学特征进行了蒙特卡罗直接模拟（ＤＳＭＣ）。采用中点数值逼近法确定模拟分子与物面碰撞点，利用网格编号将三维计算问题一维化，把流场入口处的分子补充与分子排序有效结合，利用分子横穿流场次数确定流场稳定前的时间步长计算次数。简化了计算过程，大大节约了计算机时和内存。给出了指令舱从自由分子流到过渡区不同Ｋｎｕｄｓｅｎ数、各攻角下的气动特征。模拟结果与有关理论、实验数据相比，证实     

热辐射计内稀薄气体流动的蒙特卡罗模拟

本文应用直接仿真蒙特卡罗(DSMC)方法,研究了Crookes热辐射计内由于叶片两表面的温度不同而引起的稀薄气体的流动,给出了热辐射计叶片高温面温度900K,低温表面温度为300 K,四种不同克努森数(Kn)下流速矢量场、温度场和叶片两面压力差的仿真结果.并与热辐射计在不同Kn数下的转速实验测量结果进行了比较,两者结果基本一致.     

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

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

三维管道超声速气固两相流动的CFD/DSMC仿真

采用NND格式求解气相N-S方程，采用DSMC方法模拟固相颗粒运动，数值模拟了三维管道超声速气固两相流场，分析了固体颗粒对流场结构产生的影响。计算结果表明，将DSMC方法与CFD方法相结合是解决气固两相流动问题的一种有效途径。