论文部分内容阅读
多组分反应欧拉方程是气相爆轰的控制方程，其源项存在刚性。应用显式格式求解，由于格式稳定性条件的限制，时间步长远小于CFL（Courant—Friedrichs—Lewy）条件得出的时间步长，给求解带来很大困难。引入一种IMEX（Implicit—Explicit）型附加Runge—Kutta算法，对非刚性对流项进行显式处理，源项用半隐式格式处理，整体具有高精度和L一稳定性，并进行算例考核。在此基础上，利用基元反应模型对气相爆轰问题进行了数值研究。结果表明，该算法能够很好地处理源项引起的刚性问题，准确地捕捉