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

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

大型航天器再入陨落时太阳翼气动力/热模拟分析

采用直接模拟蒙特卡洛（DSMC）方法对大型航天器离轨再入陨落过程中,其太阳翼帆板在稀薄过渡流域的气动力、气动热特性进行数值模拟,计算中采用流场直角与表面三角形非结构混合网格以及网格自适应技术处理这类复杂外形的流动模拟,考虑内能激发和化学反应来准确模拟气动加热,并基于MPI环境的并行算法解决计算量庞大的难题。通过计算分析太阳翼水平和垂直放置时在不同高度、不同攻角下的复杂流动特征,表明在90km以上高空,太阳翼垂直放置时,飞行器头部脱体激波与帆板脱体激波会产生更强烈、更复杂的激波/激波和激波/边界层的干扰,在气动力和气动热的双重作用下要比水平放置时的太阳翼更快地被撕裂并脱离目标航天器。     

基于修正N-S方程本构关系的气体动理论耦合方法

随着稀薄程度的增加,Navier-Stokes方程的线性本构关系难以正确描述稀薄气体输运特性,高阶非线性本构关系往往数学形式极为复杂,对数值求解造成稳定性差等问题。为了发展适宜于近空间飞行器气动特性分析的高超声速稀薄流动模拟方法,本文利用求解Boltzmann模型方程的气体动理论统一算法(Gas Kineitc Unified Algorithm,GKUA)对应力张量、热流等宏观量数值积分求解的优势,提出了一种基于数值修正N-S方程本构关系的气体动理论耦合方法。通过将GKUA获得的应力张量及热流用于修正N-S方程的本构关系,实现了存在局部稀薄效应的流动模拟,并且通过可压缩平板边界层、圆柱绕流问题的数值模拟,验证了方法的有效性。     

轴对称Boltzmann模型方程统一算法与空天飞行环境喷管流动模拟

拦截机动飞行器周围大范围区域存在主喷/侧喷流/羽流影响,而传统的纳维-斯托克斯（Navier-Stokes, N-S）方程不能很好模拟发动机喷管扩张段出口附近流动情况,需要一种新的方法来处理这种全流域流动问题。为解决该问题,针对特定轴对称喷管内流动,本文通过数学推导确立描述不同克努森数稀薄环境条件下的轴对称喷管内流动Boltzmann模型方程,初步建立适于该模型方程的数值格式与气体动理论统一算法。通过开展同轴圆筒间的定常/非定常旋转流动以及轴对称喷管内流动数值计算研究,发现统一算法计算流场与其他途径得到的结果吻合较好,验证了统一算法在全局克努森数喷管流动模拟的适应性和可靠性。通过与低密度风洞实验对比,喷管出口核心区羽流结构一致,羽流轴线压力分布一致,表明统一算法可以有效解决喷管入口压缩段到扩张段多流域混合,尤其是出口附近稀薄气体真空低压环境流动问题。     

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

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

地基光学图像中本体-帆板结构航天器姿态估计方法

提出一种利用地基光学图像估计空间目标姿态的方法。针对本体-帆板结构航天器外形结构特点，定义了一种不易遮挡的两轴结构特征，设计了地基光学图像中该特征的提取方法；根据航天器成像过程建立了特征点二三维几何投影方程，实现了航天器三维姿态的解算。利用单站序列图像，可获得航天器姿态变化角速度信息。在轨航天器的试验结果表明，本方法能有效地从地基光学图像估计出航天器姿态。     

小尺寸应变天平设计方法研究

Φ0.3m高超声速低密度风洞测力试验模型尺寸小、天平载荷小、流场总温高,天平防热结构设计难,刚度要求和灵敏度要求矛盾突出,测力天平的研制难度大。采用数值仿真和风洞试验相结合的研究方法,对应用于该情况下的内式天平在温度影响和结构设计方面进行了研究。采用有限元方法对天平温度影响和结构应力分布做了详细分析,并对分析结果进行了风洞试验验证。结果表明：在实现天平结构优化设计及提高天平设计的准确性和可靠性方面,有限元分析设计方法是非常有效的技术途径。     

航天器返回地球的气动特性综述

航天器返回地球的飞行过程中,气动特性是实现将宇宙飞行速度减到落地前速度、保证再入飞行得到有效控制以及再入防热安全可靠的关键因素。针对简单旋成体气动外形、半弹道式再入控制、烧蚀防热类返回航天器,综述了返回地球过程中变化的空气流域特性、航天器周围的气体绕流环境、空气与航天器作用产生的动力学与热效应等。系统地给出了该类航天器的再入气动特性参数与飞行性能的共性规律,包括:气动阻力与再入减速、气动升力与再入轨迹控制、配平攻角与飞行稳定性、气动加热与防热,以及再入过程中不同气动特性航天器、气象条件变化等对再入飞行性能的影响规律。为航天器开展返回飞行过程的跨流域气动性能工程研制提供设计参考。     

NPLS流场测量技术在高超声速风洞中纳米粒子跟随性数值仿真

从描述粒子运动的微观层次出发，采用双向耦合技术，建立了一种适用于稀薄条件下两相流动的DSMC数值模拟方法。对相间相互作用进行解耦处理，实现了气固两相间动量和能量相互作用的模拟。采用基于DSMC方法的稀薄两相流双向耦合算法，对NPLS测量技术高超声速流场测量中纳米粒子的跟随性进行了数值研究。通过Φ50nmTiO₂粒子在不同高超声速流场条件下气相-纳米粒子两相流场的仿真，表明在稀薄度很小的流场中，纳米粒子的跟随性很好。而随着流场稀薄度增加，流场中纳米粒子的跟随性降低，纳米粒子在流场中的分布与气相流场分布差异变大，通过NPLS测量得到的激光散射信号不能反映流场结构。     

稀薄流到连续流的气体运动论统一算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数的各流域均适用的气体分子速度分布函数简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而在速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。发展可用于速度空间宏观取矩的离散速度数值积分方法，获取物理空间各点的流动参数，由此发展一套能有效模拟各流域三维绕流问题的气体运动论统一算法。研究气体运动论数值算法所适合的并行方案，基于统一算法的HPF并行实现，建立一套能有效模拟不同流域复杂外形体绕流的HPF并行算法软件。通过对不同Knudsen数的一维、二维、三维气体绕流问题进行数值计算表明，计算结果与有关实验数据及其它途径得到的研究结果吻合较好，证实了本文发展的统一算法在求解稀薄流到连续流不同流域复杂绕流问题方面的可行性。