基于半隐式算法的反应欧拉方程数值计算

摘要：

在处理反应欧拉方程中源项所引起的刚性问题时，发展了一种 IM EX （Implicit-Explicit）型Additive Runge-Kutta半隐式算法，此方法能够节省计算时间和内存；该算法对方程中的非刚性对流项进行显式计算，刚性源项则用半隐式格式计算。为了验证该算法的可行性，通过对典型的马赫反射、二维气相爆轰、爆轰波的传播和爆轰胞格结构进行了数值模拟。模拟结果表明：该算法能够很好地处理由源项引起的刚性问题，准确地反应马赫反射情况，且得到的胞格结构与实验结果十分吻合。

关键词：气相爆轰、半隐式算法、刚性、马赫反射、胞格结构

分类号：O381(爆炸力学)

资助基金：国家自然科学基金11221202，11032002，11072025

在线出版日期：2014-11-28（万方平台首次上网日期，不代表论文的发表时间）

页数：共6页

页码：539-544