气体动理学格式研究进展.pdfVIP

第42 卷 第5 期 力 学 进 展 Vol. 42 No. 5 2012 年9 月25 日 ADVANCES IN MECHANICS Sep. 25, 2012 气体动理学格式研究进展 1 2 李启兵 徐 昆 清华大学航天航空学院, 北京 100084 香港科技大学数学系, 香港九龙清水湾 摘 要 介绍了近年来气体动理学格式 亦简称 格式 的主要研究进展 重 点是高阶精度动理学格式及适合从连续流到稀薄流全流域的统一动理学格式 通过对速度分布函数的高阶 展开和对初值的高阶重构 构造了时间和空间均为三阶精度的气体动理学格式 研究表明 相比于传统的基 于 解的高阶格式 新格式不仅考虑了网格单元界面上物理量的高阶重构 而且在初始场的演化阶段 耦合了流体的对流和黏性扩散 也能够保证解的高阶精度 该研究为高精度计算流体力学 格式的建立提供了一条新的途径 通过分子离散速度空间直接求解 模型方程 在每个时间步长内将宏观量的更新和微观气体分布函数的更新紧密地耦合在一起 建立了适合任意 数的统一格式 相比于已有的直接离散格式具有更高的求解效率 最后 本文还讨论了合理的物理模型对 数值方法的重要性 气体动理学方法的良好性能来自于 模型方程对计算网格单元界面上初始间断 的时间演化的准确描述 气体自由运动与碰撞过程的耦合是十分必要的 通过分析数值激波层内的耗散机制 我们认识到采用 方程的精确 解作为现代可压缩 方法的基础具有根本的缺陷 高马赫数 下的激波失稳现象不可避免 气体动理学格式为构造数值激波结构提供了一个重要的可供参考的物理机制 关键词 气体动理论 高精度动理学格式 统一动理学格式 方程 引 言 对传统直接基于宏观守恒量的格式而言, Eu- ler 方程的 Riemann 解 (及其近似解) 是演化步的 现代航空航天技术的飞速发展, 要求对高速 核心. 以此为基础, 国内外研究学者发展出了众多 飞行器面临的复杂流场进行越来越细致的时空刻 的计算流体力学方法, 在越来越多的领域取得了 画, 在这方面高分辨率、高精度格式日益显示出其 成功应用, 有力地推动了相关学科和工业应用的 独有的优势 - . 高精度有限体积格式由于能保 发展. 应用的发展反过来也对数值方法提出了越 证流场中物理量的守恒性, 且能适用于飞行器面 来越高的要求, 迫切需要开展更为深入的基础研 临的复杂网格系统, 近年来逐渐成为计算流体力 究. 其中, 一个经常被人忽视的问题是, 由于宏观 学的研究热点. 有限体积方法包括单元内物理量 方程的强非线性, 用其描述界面上初始间断的演 的重构、单元界面上重构出的初始物理量随时间 化在理论上存在先天的不足. 由于Riemann 解在 的演化以获得数值通量、通过界面通量计算单元 时间方向上为一阶精度(generalized Riemann prob- 内空间平均值的投影3 个步骤. 大量研究工作集 lem, GRP 求解器 可达到二阶), 需要结合多 中在高阶重构上, 而对初值演化模型的研究非常 步 Runge–Kutta 方法提高时间精度. 通过 Runge– 少. 事实上, 高阶格式不仅需要高阶精度的重构, Kutt