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 simulate 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 well capture the discontinuity, and accurately describe the detonation complex wave configurations such as the triple wave structure and cellular pattern.
| 科 Family | 属数 Number of genus | 种数 Number of species | 占总种数比例 Percentage of total species (%) | 属 Genus | 种数 Number of species | 占总种数比例 Percentage of total species (%) |
|---|---|---|---|---|---|---|
| 鹅膏菌科Amanitaceae | 2 | 11 | 5.26 | 鹅膏菌属 Amanita | 10 | 4.78 |
| 小菇科 Mycenaceae | 2 | 12 | 5.74 | 丝盖伞属 Inocybe | 5 | 2.39 |
| 多孔菌科 Polyporaceae | 8 | 14 | 6.70 | 蜡蘑属 Laccaria | 5 | 2.39 |
| 红菇科 Russulaceae | 3 | 23 | 11.00 | 小皮伞属 Marasmius | 6 | 2.87 |
| 小菇属 Mycena | 11 | 5.26 | ||||
| 光柄菇属 Pluteus | 5 | 2.39 | ||||
| 红菇属 Russula | 17 | 8.13 | ||||
| 栓菌属 Trametes | 5 | 2.39 |