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

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

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

本文将平衡粒模拟方法（EPSM）与蒙特卡罗直接模拟方法（DSMC）结合，构造了可模拟含近连续流区及过渡流区的DSMC／EPSM混合算法。运用混合算法模拟马赫数等于4时超声速平绕流，并与DSMC结果进行比较，证明了DSMC／EPSM混合算法的有效性。将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耦合算法相容;具有全飞行流域很好的计算一致收敛性。简述了跨流域空气动力学几种模拟手段的适应性特点与展望,揭示了返回舱再入跨流域复杂高超声速流动变化规律。     

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

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

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

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

基于Boltzmann模型方程不同流区复杂三维绕流HPF并行计算

随着现代航空航天快速发展,是否能建立一套可有效模拟稀薄流到连续流各流域绕流问题统一算法,已成为工程应用部门和学术研究领域关心的问题。通过开展基于Boltzmann简化速度分布函数方程数值求解,发展起从稀薄流到连续流统一算法。借助区域分解并行化方法研究建立气体运动论统一算法并行方案,通过对统一算法进行HPF并行化程序设计及算法考验,拟定不同流域三维球体绕流及类"神舟"返回舱外形体绕流算例,进行HPF并行计算,并将计算结果与有关实验数据、DSMC模拟值进行定量比较、分析。研究表明,所发展的统一算法具有很好的并行独立性,基本达到了线性加速的并行效果,算法负载平衡和并行可扩展性较好,可望建立起新型的能可靠模拟不同流域三维绕流问题的HPF并行算法研究方向。     

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

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

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

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

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

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

返回舱再入跨流域气动及配平特性数值研究

基于直接模拟蒙特卡罗(DSMC)方法,发展流场直角与表面三角形非结构混合网格生成方法和网格自适应技术,构造适于高稀薄流到近连续滑移流多流区共存的变时间步长模拟策略及DSMC区域分解并行计算方案。采用经改进的DSMC方法模拟飞船返回舱再入过程130～70km跨越外层空间自由分子流到稀薄过渡近连续流区热化学非平衡流动,研究分析了跨流域激波过渡带和边界层的强扰动演变规律以及稀薄气体效应、高温真实气体效应对返回舱气动特性和配平特性的影响规律。对比分析了不同壁面反射模型对返回舱配平特性的影响特点,计算的探月试验返回器配平迎角与飞行试验数据一致。计算分析了质心位置偏移对配平迎角的影响机制。计算结果表明:稀薄气体效应和壁面反射模型对90km以上高度的气动力和配平特性影响显著,真实气体效应则对90km以下高度的气动特性影响较大,质心横向偏移对配平迎角影响较大。     

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

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

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

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

DSMC量子动理学化学反应模型的数值模拟

为评估直接模拟蒙特卡罗（DSMC）方法中新近发展的一种量子动理学（QK）化学反应模型，在自研高超声速稀薄气体流动数值模拟软件RariHV中实现该模型，并进行了单个绝热单元内的热泳平衡/非平衡反应测试。结果表明，数值解与理论解吻合良好，QK模型具有很好的预测化学反应流动的能力。为进一步验证QK模型在真实流动中的表现，还计算了高超声速绕圆柱的化学反应流动，QK模型计算的驻点线上温度和组分分布都与总碰撞能（TCE）模型结果吻合良好。相较于TCE模型，QK模型的优势是不再依赖反应速率系数的实验数据，在深空探测等化学反应数据缺乏的领域具有一定的应用前景。     

二维热化学非平衡流动的非结构网格DSMC方法及其应用

研究了二维热化学非平衡流动非结构网格DSMC方法实现的过程。提出了一种新型的高效搜索算法,该算法不仅可以跟踪模拟分子在网格之间的迁移,而且在搜索过程中可以准确判别分子与物面是否相互作用,避免了原有算法中分子表面反射非确定论判据。为加快流场的时间发展历程,设计了适合DSMC方法的动态局部时间步长技术,使其可以应用到定常和非定常流场的计算。利用Fortran90的动态分配内存技术编制了计算程序。最后对过渡流域高超声速圆柱绕流进行了数值试验,计算结果初步验证了该算法的可行性。     

一类二维非结构网格DSMC方法的实现策略及其应用

研究了一类二维非结构网格DSMC方法的实现策略。在数据结构方面，设计了局部化的数据组织方式，节约了内存与计算时间。发展了一种跟踪模拟分子迁移的算法，该算法仅需少量的逻辑运算与代数运算，不仅可以快速跟踪模拟分子在网格之间的迁移，而且可以准确判别分子与物面是否相互作用，搜索过程中的附带信息给出了分子与物面碰撞的精确时间与位置，避免了重新计算。引入碰撞距离的思想，既减少存储又保证正确的模拟结果。在程序编制过程中，我们充分展现了Fortran90高级语言的主要特性，引入动态数组、指针、链表以及派生类型数据，编制了计算程序。最后对过渡流域高超声绕流进行了数值试验。