过渡区侧向喷流干扰的并行DSMC数值模拟研究

在高空低密度环境,空间飞行器的控制面效率下降较快,反作用控制系统(RCS)作为改变飞行姿态或轨道的直接力系得到大量应用,为了精确预测高空RCS喷流与稀薄大气的干扰效应,本文建立了直角与表面非结构网格混合结构的DSMC数值算法以及网格自适应算法,保证了碰撞分子是在自适应后的亚网格内选取,提高了计算精度.采用静态随机负载平衡技术构建了并行DSMC代码,计算分析了不同压比条件下的三维平板模型侧向喷流与稀薄大气的干扰流场.计算的复杂流场结构和表面流动特征、分离长度与低密度风洞试验有较好的一致性.在保持来流参数不变的情况下,喷流干扰区的分离长度及喷流透射高度随着喷流压比的增大而增大.随着来流稀薄度的增加,来流对喷流的影响越来越弱,而喷流自身的影响区域却越来越大.     

考虑多种物理效应的钝锥俯仰稳定性参数影响分析

首先采用二阶ROE格式计算了非定常NS方程下高超声速钝锥体的俯仰阻尼导数,通过与文献中内伏牛顿法和实验结果对比验证了本文所采用非定常计算方法的正确性.同时采用二阶ROE格式及双时间步长、刚性旋转动网格技术等,数值分析了NS方程计算条件下高空高马赫数钝锥体模型的强迫振动非定常流场特征,并运用最小二乘法辨识出其俯仰静稳定性导数与俯仰阻尼导数.通过计算结果对比分析了不同高度(含稀薄气体效应)、Maxwell滑移边界条件及五组元化学非平衡模型对钝锥体模型非定常流动特征及静、动导数的影响.结果分析表明,本文所计算的钝锥模型在所假设强迫振动下的俯仰稳定性参数随高度变化剧烈,且考虑滑移边界条件与化学反应会使得俯仰力矩迟滞曲线形态发生较大程度改变.准确计算高空、高马赫数条件下飞行器俯仰阻尼导数应考虑化学非平衡效应和稀薄气体效应的影响,是否考虑这些物理效应会关系到飞行器静、动稳定的程度.     

分子气体稀薄效应的动理学建模

稀薄气体效应,是指气体在特征尺度与其分子平均自由程相当的系统中流动时出现的非平衡效应.相比于单原子气体,分子气体(每个气体分子包含两个及以上原子)流动因同时具有转动、振动等多种自由度的非平衡过程,其稀薄效应更为复杂.分子气体稀薄效应在航空航天、微机电系统和页岩气开采等民生、科技领域广泛存在,而描述该效应的动理学模型与数...     

有局部稀薄气体效应的高超声速流动数值模拟

近空间高超声速飞行器当飞行高度和速度足够高时,其流场计算可能要考虑稀薄气体效应,传统的计算流体力学(CFD)方法预测的阻力和升阻比将不够准确。而现有的模拟稀薄气体流动的计算方法由于其计算量巨大,难以在工程实际中应用。因此需要发展能用于近空间高超声速飞行器流场的可行、可靠的计算方法。陈杰和赵磊在文献[1]中针对边界层中既有强剪切而气体分子自由程又相对较大的情况进行分析,提出了刻画此类局部稀薄效应的无量纲参数Zh,并提出了在传统CFD中通过采用依赖于Zh参数的等效黏性系数考虑局部稀薄效应对阻力计算影响的研究思路。因此,本文尝试将此等效黏性系数纳入CFD模型中,以在70km高空,以马赫数15飞行的小迎角钝平板为例,来检验计算方法是否合理可行。结果表明:和传统的CFD方法所得结果相比,新模型计算的阻力减小,升阻比增加,其改进的方向与现有飞行试验结果定性相符,且所增加的计算时间非常有限,可方便地应用于现有的计算空气动力学中。     

气体稀薄效应对热流计算的影响

临近空间高超声速飞行器的研发越来越受到关注。周恒和张涵信在《中国科学》上的《空气动力学的新问题》一文中分析了在此情况下现有的空气动力学的不足,其中之一就是需要考虑流场中可能出现的局部稀薄气体效应,并且针对临近空间高超声速飞行器边界层进行了详细的分析。陈杰和赵磊在此基础上研究了强剪切下的气体稀薄效应,给出了判别气体稀薄效应的无量纲参数Zh,并提出对传统连续介质模型中的黏性系数基于Zh参数进行修正。本文将通过类似的研究方法,采用DSMC (Direct Simulation Monte Carlo)研究纯导热问题中的气体稀薄效应,通过对粒子速度分布函数的分析提出了相应的刻画气体稀薄效应的参数ZhT,获得了依赖参数ZhT的导热系数修正规律,并进一步将该修正规律纳入CFD (Computational Fluid Dynamics)计算,验证了该算法对圆柱绕流问题表面热流预测的精度。     

激光光船周围连续/稀薄混合流场的数值模拟

当火箭模式激光推进器在高空稀薄大气中飞行时,需要在激光聚焦区附近喷射稠密气体工质,从而利用激光聚焦击穿气体工质形成等离子体,为光船提供推力。为了分析该光船的推力产生过程,采用连续流体力学计算方法(CFD)和直接数值模拟蒙特卡洛方法(DSMC)耦合的混合算法,数值模拟飞行高度为80 km飞行马赫数为5的激光光船周围高温连续气体和稀薄环境大气耦合的非定常混合流场,并得到推力随时间的演化曲线。结果表明,稀薄效应明显增加了激波厚度,高温效应延迟了推力的产生,并导致冲量耦合系数和比冲的降低。     

过渡领域高超声速圆柱绕流直接模拟

本文用直接模拟蒙特卡罗方法模拟了再入速度为７．５ｋｍ／ｓ，高度分别为８０、８５、９０ｋｍ，物壁是完全漫反射有限催化壁的圆柱绕流。在８０ｋｍ考察了壁面催化率的影响。     

基于非线性耦合本构关系的改进边界条件

非线性耦合本构关系（NCCR）模型是在Eu的广义流体动力学方程（GHE）基础上,通过绝热假设、Eu封闭和Myong简化推导出的关于非守恒量（黏性应力与热流）的非线性代数方程,有效拓展了线性的纳维-斯托克斯-傅里叶（NSF）本构模型在非平衡流动中的模拟能力,为快速准确模拟连续与稀薄耦合流动问题提供了强有力的理论工具。针对该模型开展滑移边界条件研究,结合努森层内物理量非线性分布的特点,提出一套在物面处与模型精度相一致的非线性修正滑移边界条件。在有限体积框架下,采用AUSMPW⁺格式和LU-SGS方法以及NCCR的完整耦合求解算法,对不同稀薄程度的高超声速单原子氩气圆柱绕流和平板绕流问题进行数值模拟。研究结果表明,基于NCCR模型的修正边界条件准确刻画出物面努森层内流动的非线性特点,有效提高了固壁滑移边界的精度。采用非线性修正边界的NCCR模型准确预测了连续流、滑移流和过渡流域的物面压力、摩阻与热流系数。     

基于能量守恒的再入轻诱饵鉴别方法

针对稀薄大气中的诱饵干扰,提出了一种基于能量守恒定律的再入轻诱饵鉴别方法.该方法利用目标单位质量机械能的损耗表征质阻比的相对大小,实现对轻诱饵的鉴别,避开了质阻比估计的收敛过程.首先根据不同目标的能量方程,推导出目标单位质量机械能损耗的解析表达式,确定了鉴别量;然后理论推导了鉴别量的方差估计式并设计了鉴别门限;最后通过仿真验证了鉴别方法的有效性.     

稀薄流热化学非平衡效应的DSMC研究

以钝体绕流驻点线流场为例，应用直接仿真蒙特卡罗法，模拟了五组元混合气体（Ｎ２，Ｏ２，Ｎ，Ｏ，ＮＯ）中的碰撞能松弛，离散反应，交换反应及复合反应，给出了碰撞分子各自由度间的能量交换模式及各类化学反应的选择模式。ＤＳＭＣ仿真结果与宏观粘性激波层方程计算结果，在飞行高度较低时基本相符，表明采用的热化学非平衡非模式是令人满意的，随着飞行高度增高稀薄气体效应增大，二者在激波附近有显著差异，表明高空飞行必须考