二维Boltzmann-Rykov模型方程数值算法.pdf

第 3O卷 第 3期 计 算 物 理 Vo1．30．NO．3 2013年 5月 CHINESE JOURNAL OF COMPUTATIONALPHYSICS May，2013 文章编号：1001—246X(2013)03—0326—11 考虑转动能的一维／二维 Boltzmann—Rykov 模型方程数值算法 吴俊林 ， 李志辉 ， 蒋新宇 (1．中国空气动力研究与发展中心超高速研究所，绵阳 621000 2．国家计算流体力学实验室，北京 100191) 摘 要 ：研究考虑转动能的Boltzmann—Rykov模型方程，基于转动 自由度对气体分子速度分布函数矩积分，引入约 化速度分布函数 ，应用离散速度坐标法与数值积分技术 ，将气体运动论模型方程化为在离散速度坐标点处关于三 个约化速度分布函数的联立方程组．应用拓展计算流体力学有限差分方法，数值计算考虑转动 自由度的双原子气 体一维、二维Bohzmann模型方程，得到高、低 Knudsen数一维激波管内流动和二维竖直平板绕流问题的流场，分析 验证考虑转动能的Bohzmann—Rykov模型方程全流域统一算法求解一维 ／二维气体流动问题的可靠性．结果表明，气 体稀薄程度与分子 内自由度对流场具有较大影响，且 Knudsen数较高的稀薄气体流动呈现严重的非平衡流动特 点． 关键词 ：转动 自由度 ；Boltzmann—Rykov模型方程 ；激波管；竖直平板绕流 ；统一算法 中图分类号：0246；0355 文献标识码 ：A 0 引言 描述气体分子运动论的基本方程是 Boltzmann方程 ，它给出速度分布函数对位置空间和时间的变化 率关系．该方程可描述稀薄流到连续流不同流域的气体分子输运现象．因此 ，它可以作为求解各流域气体流 动问题的基本方程．通过求解该方程 ，可以提供远离局部热力学平衡态的气体流动详细特征 I2，可用于研 究往返大气层不同高度下飞行器飞行绕流状态 ．同时，已有研究表 明，Bohzmann方程是研究微尺度流 动 及多种复杂物理过程 的主要依据之一．在 Bohzmann方程于 1872年提 出以后的很长一段时间里 ， 求解这个方程却没有取得很好的进展，主要原因在于Boltzmann方程中的几率密度分布函数，依赖于时间t、 空间坐标 ，和速度 ，属于七维变量．同时，方程右端碰撞项具有非线性 、高维积分属性 ，属于复杂多尺度的 刚性问题 ¨I3．因此，精确求解描述各流域气体流动特征的Bohzmann方程难以实现． 众多学者从微观分子运动论原理出发 ，基于质量 、动量、能量守恒定律 ，利用微观分子运动、碰撞趋于平 衡态的基本特性，用数学上较简单的统计和碰撞松弛模型代替Bohzmann方程碰撞项，提出了许多气体分子 运动论模型方程．1954年由Bhatnagar、Gross和 Krook提出的BGK模型方程 为稀薄气体流动 问题提供了 一 个有效的处理方法．随后 ，人们发展 了多种形式的模型方程 ，如 Holway的椭球统计模型 ¨ 、Shakhov基于 BGK方程修正得到的高阶推广模型方程 、Segal和 Ferziger的多项式模型 以及 由Abe和 Oguehi提 出的 类似于 Shakhov模型的谱系运动论模型方程 等．此外 ，考虑气体分子转动能内自由度影响的模型方程也 被提出来 ，如ES．BGK模型 10]、Morse模型 和 Rykov模型 ．Holway通过将 Bohzmann方程碰撞项分解为 弹性与非弹性碰撞项两部分 ，其 中弹性碰撞时分子 内能保持不变 ，而非弹性碰撞过程则需要考虑内部能级输 运，以此构建了ES—BGK模型方程 ．而Morse从Wang-ChangUhlenbeek(WCU)方程 出发，采用修正的小 扰动技术描述流动的时间尺度 (时间上的Chapman-Enskog渐近展开)，非弹性松弛时间需要按经验性参数处 理，得到计算热传导矢量的Morse模型方程 ．ES—BGK模型和Morse模型在 内部能级处理有相似之处，但 在弹性碰撞项的处理上 Mors