理论解用于复杂构型飞行器气动阻力方法研究

针对现有的理论解方法尚不适用于复杂构型飞行器气动力求解的问题,从稀薄气体分子的动力学模型出发,制作等效平面以分析各部分来流的宏观特性,并将入射流和反射流作用分开进行考虑,从而计算出几种典型情况下的气动阻力值。后续验证表明,在分子-壁面间多次碰撞不显著时,此方法能够获得较准确的气动力结果。所提出的气动力计算方法对低轨道飞行器的气动力计算及气动构型优化等工作均具有一定的参考价值。     

稀薄流航天器鼻锥迎风凹腔气动力和气动热性能研究

本文采用直接模拟蒙特卡罗(DSMC)方法,对高超声速稀薄流中航天器鼻锥迎风凹腔气动力与气动热性能进行了数值研究。得到了鼻锥外壁面、凹腔侧壁面以及凹腔底面的热流密度分布,分析了不同凹腔深宽比对鼻锥冷却效率以及凹腔腔体内气体参数的影响;以深宽比为1的凹腔为基准,研究了凹腔唇口钝化半径对航天器气动热与气动力的影响。数值结果表明,高超声速稀薄流中迎风凹腔能够降低鼻锥外壁面的热流密度;当凹腔深宽比达到1之后,凹腔侧壁面热流变化趋于一致,热流密度最低点的轴向位置不随深宽比改变,且凹腔底部热流很小;凹腔近底部气体均由稀薄流转化为连续流,腔内气体压力不断振荡;唇口钝化没有明显优势,虽然可以降低鼻锥峰值热流,但是会带来严重的气动力性能下降。     

热力学非平衡状态下DSMC绝热壁边界的抽样方法

从直接模拟蒙特卡罗仿真DSMC方法的实现途径出发,经过分析推导,给出了一种在考虑分子转动能、振动能等热力学非平衡状态下统计求解绝热壁温的DSMC抽样方法。然后采用锥体外形的算例验证,分析了该方法在实现绝热壁条件时的统计误差。计算结果表明,此种方法能够有效地实现DSMC方法壁面处的绝热壁条件,随着采样步数的增加,其边界条件实现时的统计误差不断降低,该方法收敛后的统计误差可降至0.5%,能够满足工程要求。     

基于稀薄气-粒两相流的蒙特卡洛颗粒辐射模型研究

在气-粒两相相变模型及液态和固态颗粒碰撞、聚合和分离模型的基础上,发展稀薄条件下考虑颗粒辐射的蒙特卡洛颗粒辐射模型。通过对高超声速稀薄环境中的气-粒两相喷流流场的数值模拟,得到气-粒两相流的流场参数,利用所得流场参数作为颗粒辐射模型的初始参数进行颗粒辐射计算,同时考虑了有无探照发射时的光谱辐射强度。结果表明,在颗粒浓度较大时计算两相稀薄流的流场参数,考虑颗粒辐射是必要的,并且考虑有无探照发射对光谱辐射强度数值的影响。     

稀薄气体圆柱绕流的并行化DSMC方法研究

本文针对ＤＳＭＣ方法的内在并行特征构造了两种并行算法,并在ＰＶＭ环境下实现了并行化计算。计算结果表明,实现ＤＳＭＣ方法的并行化能够显著缩短计算时间,提高计算效率,不同并行算法能够达到的加速比存在差异。     

高压捕获翼构型的跨流域气动特性

高压捕获翼(HCW)构型是一种满足高速飞行器高容积、高升力、高升阻比的设计需求的新型气动布局.最近研究表明,HCW构型能够提高飞行器在连续流区的升力和升阻比,缓解飞行器设计中高容积率与高升阻比间的矛盾.为探究该气动布局在过渡流域(70～100 km)的气动特性,以一种楔—平板组合的高压捕获翼原理性构型作为模型,采用直接...     

空间CO2羽流压力场的数值计算和实验研究

要解决空间羽流污染问题,首先必须精确掌握羽流场的特征。为此,在真空深冷背景环境下进行了CO2气体模拟发动机羽流压力场的数值计算和实验研究。采用经典的计算流体动力学(CFD)和直接模拟Monte Carlo(DSMC)相结合的方法进行羽流场的数值模拟。试验模拟空间高度在100公里以上,在整个测试过程真空舱内压力维持在10-3Pａ量级,冷阱温度为93±5Ｋ。计算结果与测试数据吻合较好,表明文中计算方法和物理模型可以准确描述真空深冷空间环境下的羽流场特征。     

霍尔推力器等离子体羽流粒子模拟

建立了霍尔推力器羽流仿真模型,用单元粒子-直接模拟蒙特卡罗（PIC-DSMC）混合方法对SPT-7推力器的流场进行数值模拟,分析背压、扩张角和电子温度对流场的影响。结果表明：背压粒子增加了回流区内离子和高速粒子,加重羽流污染。SPT-70推力器羽流出口处扩张角约为30°。实验数据验证了仿真模型的正确和方法的可行,对电推力器及其羽流污染等的研究有一定的参考价值。     

二维Woods-Burnett方程的稳定性分析

Woods-Burnett方程是Boltzmann方程的二阶近似,是Burnett方程的一种修正,能够描述轻微偏离热力学平衡时的稀薄气体流动.但是Woods-Burnett方程在小扰动下不稳定,这是限制Woods-Burnett方程广泛应用的一个重要原因.本文在一维稳定性分析的基础上,通过线性小扰动理论,首次得到了二维Woods-Burnett方程的稳定性特征方程,并把稳定性方程的解表示在复平面上,得到了二维稳定性特征曲线.通过扰动增长系数和扰动波数的关系,得到二维Woods-Burnett方程的临界努森数为0.130,而一维Woods-Burnett方程的临界努森数为0.184.说明在二维情况下Woods- Burnett方程更加不稳定.     

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

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