BoltzmannRykov模型的有限体积方法计算.PDFVIP

　 ， ３５ ２ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１４年２月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３５，Ｎｏ．２，Ｆｅｂ．１５，２０１４ 文章编号：１０００⁃０８８７（２０１４）０２⁃０１２１⁃０９ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ Ｂｏｌｔｚｍａｎｎ⁃Ｒｙｋｏｖ模型的有限体积方法计算∗ １ １，２ １ １ 吴俊林 李志辉 彭傲平 蒋新宇 ，　 ，　 ，　 （１． 中国空气动力研究与发展中心 超高速空气动力研究所，四川 绵阳６２１０００； ２． 国家计算流体力学实验室，北京 １００１９１） 摘要：　 构建一种三阶精度的有限体积格式，数值求解考虑转动非平衡影响的Ｂｏｌｔｚｍａｎｎ⁃Ｒｙｋｏｖ模 型方程．针对模型方程的速度空间离散得到各个离散速度坐标点上彼此独立的控制方程组，运用高 阶精度的半离散化有限体积格式在位置空间对离散控制方程进行数值求解，时间项采用三阶 Ｒｕｎｇｅ⁃Ｋｕｔｔａ方法推进，方程右端二体碰撞项采用中心近似技术．该有限体积格式在气体分子对流 运动项上具有三阶精度，同时保证了分布函数的正定性和流通量守恒．计算结果与有限差分方法数 值模拟结果和连续流区非定常激波管问题的Ｒｉｅｍａｎｎ精确解均吻合较好，说明基于有限体积法的 Ｂｏｌｔｚｍａｎｎ⁃Ｒｙｋｏｖ模型方程数值求解过程是正确的． 关　 键　 词：　 Ｒｙｋｏｖ模型方程；　 统一算法；　 有限体积方法 中图分类号：　 Ｏ３５６；Ｖ２１１．２５　 　 　 文献标志码：　 Ａ ｄｏｉ：１０．３８７９／ ｊ．ｉｓｓｎ．１０００⁃０８８７．２０１４．０２．００２ 引　 　 言 稀薄气体动力学的研究主要有确定论方法和概率论方法两种．其中，基于概率统计规律的 Ｍｏｎｔｅ⁃Ｃａｒｌｏ方法将分子运动与碰撞解耦，采用直接模拟仿真分子运动与相互作用规律的方法 描述稀薄气体的流动特征，在稀薄流场数值计算中得到了非常广泛的应用，其代表是由Ｂｉｒｄ ［１］ 发展建立的直接模拟Ｍｏｎｔｅ⁃Ｃａｒｌｏ方法（即ＤＳＭＣ方法） ． 这种方法有很多优点：１） ＤＳＭＣ方法的物理模拟本质使得它能够再现分子运动与碰撞过 程，通过引入内能、化学反应、热辐射、电离复合反应等多种物理作用模型能够较为方便地实现 ［２］ 对复杂物理化学过程的描述；２） 该方法并不需要速度空间的截取边界 ，因此速度点能很好 地独立于物理问题；３） 由于粒子总是集中于分布函数较大的地方，内存空间就不会浪费在实 际分布函数接近于０ 的函数位置，使得基于仿真分子运动碰撞规律的ＤＳＭＣ方法在远离热力 学平衡态的流动问题中具有极大的优势；４） 相对于基于Ｂｏｌｔｚｍａｎｎ方程直接求解的确定论方 法来说，在同等离散自由度数下Ｍｏｎｔｅ⁃Ｃａｒｌｏ方法的计算代价要小得多． 然而，虽然ＤＳＭＣ方法在稀薄气体动力学研究中具有非常广泛的市场，它的一些缺陷也限 ∗ 收稿日期：　 ２０１３⁃０７⁃０１；修订日期：　 ２０１３⁃１２⁃０９ 基金项目：　 国家自然科学基金；国家重点基础研究发展计划（９７３计划）（２０１４ＣＢ７４４１００） ＋ 作者简介：　 吴俊林（１９８５—），男，四川会理人，硕士（Ｔｅｌ： ８６⁃８１６⁃２４６５２６１；Ｅ⁃ｍａｉｌ：ｗｕｊｕｎｌｉｎ１３０＠ ａｌｉｙｕｎ．ｃｏｍ）； 李志辉（１９６８—），男，四川眉山人，研究员，博士，博士生导师（通讯作者． Ｅ⁃ｍａｉｌ：ｚ