超声速平板绕流DSMC/EPSM混合算法研究

本文将平衡粒模拟方法（EPSM）与蒙特卡罗直接模拟方法（DSMC）结合，构造了可模拟含近连续流区及过渡流区的DSMC／EPSM混合算法。运用混合算法模拟马赫数等于4时超声速平绕流，并与DSMC结果进行比较，证明了DSMC／EPSM混合算法的有效性。将DSMC方法与混合算法的计算效率进行了比较。     

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

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

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

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

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

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

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

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

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

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

模拟真空羽流场的特征线法

介绍了特征线法（MOC法）对真空羽流场进行数值模拟的过程，用MOC法进行了喷管出口流动不均匀情况真空羽流场算例的求解，所得结果与理论分析相符合，与DSMC法计算结果吻合，提高了计算效率，节省了计算时间。结果表明：用MOC法可代替DSMC法对连续介质区进行模拟，并可在用N-S和DSMC耦合数值模拟中确定耦合边界。     

跨流域空气动力学模拟方法与返回舱再入气动研究

针对回收类航天器(返回舱)再入过程所遇跨流域多尺度非平衡绕流问题,综述基于Boltzmann方程碰撞积分物理分析与可计算建模,构造考虑完全气体、转动非平衡、含振动能激发热力学非平衡效应各流域统一Boltzmann模型方程,及由此建立返回舱再入气动力热绕流问题气体动理论统一算法研究进展与算法检验。作为方法间验证结合,进一步简述了融合再入热化学稀薄气体电离非平衡流动DSMC方法、近连续过渡流区N-S/DSMC耦合算法、经滑移边界修正的N-S方程解算器、低密度风洞实验测试等多种空气动力学模拟手段,建立求解Boltzmann模型方程气体动理论统一算法(GKUA)、DSMC、N-S/DSMC、滑移N-S解算器、低密度风洞实验验证补充,适于返回舱再入从外层空间自由分子流到近地面连续流跨流域空气动力学一体化模拟平台。将此平台用于再入H=110～30km各流域球体、高超声速尖前缘中空柱裙、返回式卫星球锥体、飞船返回舱稀薄过渡流以至近连续流区气动力/热与姿态配平绕流问题计算与实验分析比较,证实统一算法在高稀薄流区,与DSMC吻合很好;在连续流区,与(滑移)N-S解算器相一致;在中间过渡带,与N-S/DSMC耦合算法相容;具有全飞行流域很好的计算一致收敛性。简述了跨流域空气动力学几种模拟手段的适应性特点与展望,揭示了返回舱再入跨流域复杂高超声速流动变化规律。     

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

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

解耦N-S/DSMC方法计算推力器真空羽流的边界条件研究

准确模拟羽流流场是合理评估空间推力器真空羽流效应的基础。因羽流流场同时包含连续和稀薄两种流动机理,通常采用解耦方式的Navier-Stokes(N-S)方程和直接模拟蒙特卡罗(DSMC)混合方法对其进行模拟。为保证解耦N-S/DSMC方法应用于羽流计算时的准确性并尽量提高其计算效率,对计算中的DSMC入口和喷管壁面等边界条件设置问题开展了研究。通过与文献中低总压工况的羽流试验数据比较,确定了合理设置DSMC入口边界位置、壁面反射类型、壁温等边界条件的方法。针对将此方法应用于实际推力器工况计算时效率过低的问题,通过多种数值试验,表明从喷管出口处开始沿KnGL为0.05的等值线作为DSMC入口界面可同时保证仿真精度和较高的计算效率;并证明实际工况下壁温设置对羽流场仿真结果影响不大。     

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

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

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

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

一种基于结构化贴体网格的DSMC方法研究

基于DSMC方法,发展了一种结构化的贴体网格的分子运动轨迹跟踪和定位的方法———网格面法向量判断法。在此基础上建立了一套DSMC数值模拟程序,对球双锥和Apollo飞船在过渡区域流动进行了数值模拟,并分析了其过渡区域超音速、高超音速的气动力、气动热等问题。通过将本文的计算结果与参考文献的计算结果和实验结果进行了比较,验证了所发展的方法的可靠性。     

高空羽流场的DSMC计算和实验研究

以小推力发动机的高空羽流场为研究对象，实现了高空小喷管氮气流DSMC方法数值模拟，验证了DSMC方法在高空膨胀流中的可用性，并在现有设备基础上，进行一组高空羽流试验，给出试验与计算的对比结果。结果表明：喷管出口附近能得到满意结果，由于实现设备能力所限，远场试验没有达到预期效果。     

Extension of the low diffusion particle method for near-continuum two-phase flow simulations

The low diffusion (LD) particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo (DSMC) method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.     

姿控发动机内流场及高空羽流流场的DSMC一体化数值模拟

通过引入“平衡抽样方法”有效解决了DSMC求解近连续流区所遇到的碰撞顶计算效率问题;对于过渡流及自由分子流区,构造了适合于DSM仿真的分子碰撞传能模型并配置了一种光滑过渡以减少误差的权函数。在此基础上首次实现了对姿控发动机喷管内流场及包括倒流区在内的高空弱流流场的DSMC方法一体化数值模拟。