基于附加Runge-Kutta方法的高精度气相爆轰数值模拟

李健; 郝莉; 宁建国

doi:10.11858/gywlxb.2013.02.010

基于附加Runge-Kutta方法的高精度气相爆轰数值模拟

doi: 10.11858/gywlxb.2013.02.010

李健¹,
郝莉^2, ,,
宁建国¹

1.
北京理工大学爆炸科学与技术国家重点实验室，北京 100081；
2.
北京建筑工程学院理学院，北京 100044

详细信息

通讯作者:
郝莉 E-mail:haoli_haoli@126.com

计量
- 文章访问数: 6848
- HTML全文浏览量: 441
- PDF下载量: 418
出版历程
- 收稿日期: 2013-03-20
- 修回日期: 2013-03-20
- 发布日期: 2013-04-15

Additive Runge-Kutta Methods for Accurate H2-O2-Ar Detonation Simulation

LI Jian¹,
HAO Li^{2
, ,},
NING Jian-Guo¹

1.
State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China;
2.
Science School, Beijing University of Civil Engineering and Architecture, Beijing 100044, China

摘要

摘要: 多组分反应欧拉方程是气相爆轰的控制方程，其源项存在刚性。应用显式格式求解，由于格式稳定性条件的限制，时间步长远小于CFL（Courant-Friedrichs-Lewy）条件得出的时间步长，给求解带来很大困难。引入一种IMEX（Implicit-Explicit）型附加Runge-Kutta算法，对非刚性对流项进行显式处理，源项用半隐式格式处理，整体具有高精度和L-稳定性，并进行算例考核。在此基础上，利用基元反应模型对气相爆轰问题进行了数值研究。结果表明，该算法能够很好地处理源项引起的刚性问题，准确地捕捉和描述爆轰波的复杂结构和典型特征，同时也能够很好地模拟爆轰波在楔面上的马赫反射问题。
- 源项 /
- 附加Runge-Kutta方法 /
- 刚性 /
- 气相爆轰
Abstract: 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.
- source term /
- additive Runge-Kutta methods /
- stiffness /
- gaseous detonation

HTML全文

参考文献(22)

Zhong X L. Additive semi-implicit Runge-Kutta methods for computing high-speed no equilibrium reactive flows [J]. J Comput Phys, 1996, 128(1): 19-31.

Yanenko N N. The Method of Fractional Steps [M]. New York/Berlin: Springer-Verlag, 1971.

Oran E S, Boris J P. Numerical Simulation of Reactive Flow [M]. New York: Elsevier Science, 1987.

Ma T B, Wang J, Ning J G. A hybrid VOF and PIC multi-material interface treatment method and its application in the penetration [J]. Science China Physics, Mechanics Astronomy, 2009, 39(9): 1185-1194.

Fei G L, Ma T B, Hao L. Large-scale high performance computation on 3D explosion and shock problems [J]. Appl Math Mech (Eng Ed), 2011, 32(3): 375-382.

Hu X Y, Zhang D L, Jiang Z L. Numerical simulation of gaseous detonation with a detailed chemical reaction model [J]. Acta Aerodynamics Sinica, 2003, 21(1): 59-66. (in Chinese)

胡湘渝, 张德良, 姜宗林. 气相爆轰基元反应模型数值模拟 [J]. 空气动力学学报, 2003, 21(1): 59-66.

Kennedy C A, Carpenter M H. Additive Runge-Kutta schemes for convection-diffusion-reaction equations [J]. Appl Numer Math, 2003, 44(1/2): 139-181.

Xia Y, Xu Y, Shu C W. Efficient time discretization for local discontinuous Galerkin methods [J]. Discrete and Continuous Dynamics Systems-Series B, 2007, 18(3): 677-693.

Westbrook C K. Hydrogen oxidation kinetics in gaseous detonation [J]. Combus Sci Technol, 1982, 29(1/2): 67-81.

Arajo A L, Murua A, Sanz-Serna J M. Symplectic methods based on decompositions [J]. SIAM J Numer Analys, 1997, 34(5): 1926-1947.

Zhang Y J, Liu Y D, Zhang Z H, et al. Computation on deflagration of hydrogen air mixture inside a duct and induced shell response [J]. Science China Physics, Mechanics Astronomy, 2009, 39(9): 1173-1184.

Jiang G S, Shu C W. Efficient implementation of weighted ENO schemes [J]. J Comput Phys, 1996, 126(1): 202-228.

Wang C, Ma T B, Lu J. Influence of obstacle disturbance in a duct on explosion characteristics of Coal Gas [J]. Science China Physics, Mechanics Astronomy, 2009, 39(9): 1248-1257.

Ascher U M, Ruuth S J, Spiteri R J. Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations [J]. Appl Numer Math, 1997, 25(2/3): 151-167.

Calvo M P, de Frutos J, Novo J. Linearly implicit Runge-Kutta methods for advection-reaction-diffusion equations [J]. Appl Numer Math, 2001, 37(4): 535-549.

Fritzen P, Wittekindt J. Numerical solution of viscoplastic constitutive equations with internal state variables, Part Ⅰ: Algorithms and implementation [J]. Math Methods Appl Sci, 1997, 20(16): 1411-1425.

Wang G, Wang J T, Liu K X. New numerical algorithms in SUPER CE/SE and their applications in explosion mechanics [J]. Science China Physics, Mechanics Astronomy, 2009, 39(9): 1214-1220.

Hu X Y. On the structures of gaseous detonation waves [R]. Beijing: The Institute of Mechanics, Chinese Academy of Sciences, 2001. (in Chinese)

胡湘渝. 气相爆轰波结构研究 [R]. 北京: 中国科学院力学研究所, 2001.

Wang G, Zhang D L, Liu K X, et al. An improved CE/SE scheme for numerical simulation of gaseous and two-phase detonations [J]. Comput Fluids, 2010, 39(1): 168-177.

Guo C M, Zhang D L, Xie W. The Mach reflection of a detonation based on soot track measurements [J]. Combustion and Flame, 2001, 127(3): 2051-2058.

施引文献

资源附件(0)

访问统计