科研之友
微信
新浪微博
Facebook
分享链接
摘要
在处理反应欧拉方程中源项所引起的刚性问题时,发展了一种IMEX(Implicit-Explicit)型Additive Runge-Kutta半隐式算法,此方法能够节省计算时间和内存;该算法对方程中的非刚性对流项进行显式计算,刚性源项则用半隐式格式计算。为了验证该算法的可行性,通过对典型的马赫反射、二维气相爆轰、爆轰波的传播和爆轰胞格结构进行了数值模拟。模拟结果表明:该算法能够很好地处理由源项引起的刚性问题,准确地反应马赫反射情况,且得到的胞格结构与实验结果十分吻合。
-
单位爆炸科学与技术国家重点实验室; 北京理工大学
