在处理反应欧拉方程中源项所引起的刚性问题时,发展了一种IMEX(Implicit-Explicit)型Additive Runge-Kutta半隐式算法,此方法能够节省计算时间和内存;该算法对方程中的非刚性对流项进行显式计算,刚性源项则用半隐式格式计算。为了验证该算法的可行性,通过对典型的马赫反射、二维气相爆轰、爆轰波的传播和爆轰胞格结构进行了数值模拟。模拟结果表明:该算法能够很好地处理由源项引起的刚性问题,准确地反应马赫反射情况,且得到的胞格结构与实验结果十分吻合。
This paper introduces an IMEX (implicit-explicit) Additive Runge-Kutta method,which was applied to resolve the stiff source term into the stiff part and non-stiff part, and the stiff part was solved implicitly while the non-stiff part was handled explicitly. The numerical calculation of the gas- phase detonation show that the method can deal with the stiff problems caused by the source term, ac- curate capture the discontinuity and description of the ;structure of the Mach reflection, meanwhile the features of cellular pattern agree well with the experimental results.