Additive Runge-Kutta methods for H_2/O_2/Ar detonation with a detailed elementary chemical reaction model  

作　　者：LI Jian REN HuiLan NING JianGuo 

机构地区：[1]State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology

出　　处：《Chinese Science Bulletin》2013年第11期1216-1227,共12页

基　　金：supported by the National Basic Research Program of China(2010CB832706);National Natural Science Foundation of China(11032002, 11172045)

摘　　要：We report here the additive Runge-Kutta methods for computing reactive Euler equations with a stiff source term, and in particular, their applications in gaseous detonation simulations. The source term in gaseous detonation is stiff due to the presence of wide range of time scales during thermal-chemical non-equilibrium reactive processes and some of these time scales are much smaller than that of hydrodynamic flow. The high order, L-stable, additive Runge-Kutta methods proposed in this paper resolved the stiff source term into the stiff part and non-stiff part, in which the stiff part was solved implicitly while the non-stiff part was handled explicitly. The proposed method was successfully applied to simulating the gaseous detonation in a stoichiometric H2/O2 /Ar mixture based on a detailed elementary chemical reaction model comprised of 9 species and 19 elementary reactions. The results showed that the stiffly accurate additive Runge-Kutta methods can capture the discontinuity well, and describe the detonation complex wave configurations accurately such as the triple wave structure and cellular pattern.We report here the additive Runge-Kutta methods for computing reactive Euler equations with a stiff source term, and in particu- lar, their applications in gaseous detonation simulations. The source term in gaseous detonation is stiff due to the presence of wide range of time scales during thermal-chemical non-equilibrium reactive processes and some of these time scales are much smaller than that of hydrodynamic flow. The high order, L-stable, additive Runge-Kutta methods proposed in this paper resolved the stiff source term into the stiff part and non-stiff part, in which the stiff part was solved implicitly while the non-stiff part was handled explicitly. The proposed method was successfully applied to simulating the gaseous detonation in a stoichiometric HE/O2/Ar mix- ture based on a detailed elementary chemical reaction model comprised of 9 species and 19 elementary reactions. The results showed that the stiffly accurate additive Runge-Kutta methods can capture the discontinuity well, and describe the detonation complex wave configurations accurately such as the triple wave structure and cellular pattern.

关 键 词：化学反应模型 添加剂 塔 引爆 气相爆轰 热化学非平衡 应用程序 时间尺度 

分 类 号：O621.25[理学—有机化学]