稀薄流利到连续流的气体运动论算法的研究.pdfVIP

第十届全匡高超声速气动力(热)学术交流会 稀薄流到连续流的气体运动论算法研究 李志辉张涵信 中国空气动力研究与发展中心．四川绵阳621000 捅望 本文基于气体运动论模型方程，提出一套直接求解从稀薄流到连续流统一的俺化速度 分布函数方程数值算法．我们从非线性模型Boltzmaml方程出发，引入简化速度分布函薮、使用离 散速度坐标法对速度空间进行离彦(、降维，去掉分布函数对速度分量的连续依麓性：采用对同分裂 法，将简化速度分布蠡数松弛变化方程分解为源项磋擅变化方程、对流运动方程．进抒蔫合计算： 利用^丹刀耗散差分格式求解对流运动方程．发展用于高散速度空间数值积分的高撕一埃尔米特无穷 积分法、等均间隔的牛顿一柯斯复合积分法和高新一勒让营敦值积分法．遗过害观取矩获取钧理空同 各点的流动参数，从而给出一种从稀薄流判连续漉绕任意外形体的气动特性数值算法．为了斜群奉 文的数值方法。并且证实它在求解稀薄气体动力学目题方面的可靠性．我们计算了各十知蝴教 下的一维激渡管问题及二维园柱绕流． 关键词气体分子运动论离教速度坐标法勘，￡埘册模型方程有限差分方法 0引言 航天飞行器再入过程先后经历自由分子流区、过渡漉区、滑流和连续流区。所遢气动环 境十分复杂。为了研究航天飞行器再入过程气动特征．传统的做法是稀薄it用稀薄滚的一套 算法，如基于微观分子气体动力学将分子运动、碰撞耦合模拟，依托羹特一卡罗(Mon据-Co,lo) 方法求解稀薄气体流动闯题【l。21：而连续流有连续流的一套算法，如从宏观气体流动角度出 发，以气体密度、温度、流动速度作为独立变量，数值求解欧拉(Euler)方程、纳维尔—豢 托克斯(Ⅳ一s)方程；但是两者相差很大，且计算结果很难接连。是否能找到一种连接连续 流到稀薄流能有效模拟各个流域气动特性的数值算法，将是本文研究的目的． 子输运现象，它是研究气体动力学的一个主要根据13I。但是，基于分子速度分布函数为变量 从微观角度描述流动特征的玻尔兹曼方程是一个高度非线性积分———徽分方程，精确求肆描 述复杂流动的玻尔兹曼方程未成现实。与其相反，众多学者基于质量、动量、能量守恒定律。 利用分子运动、碰撞趋于平衡态的基本特性，由数学上较简单的统计和碰撞松弛模型代替玻 尔兹曼方程碰撞项，提出了许多类似于玻尔兹曼方程各阶矩的运动论模型方程M．而描述连 数流动，现在研究发现适当修正其碰撞松弛参数，也能定性描述连续流场．因此，作为一个 定性的方法不是求解完全玻尔兹曼方程，而是求解经修正的BGK碰撞模型方程。完全有可 能得到一种计算稀薄流到连续流的统一数值算法。 九十年代，有学者通过BGK模型方程，根据捕捉激波的连续流差分格式思想，基于两 格界面通量守恒，利用气体分子速度分布函数对于Maxwell分布的渐近展开，提出了用于可 压缩流动模拟的一种新型BGK类型格式，并成功地模拟了一维激波结构、：=维激波绕射等 问题[Zsl。对于稀薄气体流动计算，国外有学者使用简化速度分布函数和离散速度坐标法，将 法f9_“1。 本研究中，我们以修正的BGK运动论模型方程为基础，引入简化速度：分布函数，建立 一组从稀薄流到连续流统一的简化速度分布函数方程；将离散速度坐标法应用于分布函数。 以期利用一套在物理空间和时间上连续而在速度空间离散的分布函数来代替原分布函数对速 溢蠢0嘧t÷自逡 度空间的连续依赖性，将气体运动论模型方程化为具有非线性源项的双曲守恒方程；采用时 间分裂法的思想，将分布函数松弛变化过程分解为对流运动和碰撞松弛过程；使用改进的欧 拉方法数值模拟源项碰撞变化方程，利用对于捕捉激波和膨胀波等典型物理现象具有高分辨 积分法，对离散速度空间进行相关的数值积分获取物理空间各点的宏观流动参数，从而给出 一种基于运动论模型方程，能有效解决各流域气动问题的数值算法。为了证实本文数值方法 在求解稀薄气体动力学问题方面的可靠性，我们将各个Knudsen数下的一维激波管问题进行 计算，在连续流区与黎曼精确解进行比较；并将该算法用于园柱绕流的数值模拟。 1 简化速度分布函数方程的建立 考虑气体运动论Boltzmann模型方程【5】