A Cell-centered Lagrangian Scheme Based on Characteristics Theory for Condensed Explosive Detonation

doi:10.19596/j.cnki.1001-246x.7913

Abstract

Abstract: We present a cell-centered Lagrangian scheme for numerical simulation of condensed-explosive detonation. It utilizes finite volume method to discrete detonation equations. Velocity and pressure of grid nodes are obtained with characteristics theory of hyperbolic partial differential equations. They are used to update position of grid nodes and calculate numerical flux of grid cells. Solution of grid nodes obtained by characteristics theory is a "genuinely multi-dimensional" theoretical solution, which is a generalization of one-dimensional Godunov scheme in two-dimensional Riemann problem. Semi-discrete system of detonation equations obtained by finite volume scheme is solved with an implicit-explicit Runge-Kutta method. Convection terms are explicitly treated, and stiff source terms of chemical reactions are implicitly treated. The numerical schemes select ZND model for detonation, JWL equations of state for unreacted explosives and detonation products, and use Ignition-Growth model to simulate evolution process in reaction zone. Typical numerical results show that the cell-centered Lagrangian scheme simulates condensed-explosive detonation problems well.

Key words: condensed-explosive detonation, cell-centered Lagrangian scheme, characteristics theory, chemical reaction rate

CLC Number:

O381

LI Shiyao, YU Ming. A Cell-centered Lagrangian Scheme Based on Characteristics Theory for Condensed Explosive Detonation[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2019, 36(5): 505-516.

References

[1] MADER C L. Numerical modeling of explosives and propellants[M]. 2nd ed. New York:CRC Press, 1998.
[2] BENSON D J. Computational methods in Lagrangian and Eulerian hydrocodes[J]. Computational Methods of Applied Mechanics and Engineering, 1992, 99:235-394.
[3] VON NEUMANN J, RICHTMYER R D. A method for the numerical calculations of hydrodynamical shocks[J]. Journal of Applied Physics, 1950, 21:232-238.
[4] WILKINS M L. Calculation of elastic-plastic flow:Methods in computational physics[M]. New York:Academic Press, 1964.
[5] JONAS Z. Introduction to hydrocodes[M]. Amsterdam:Elsevier Press, 2004.
[6] CARAMANA E J, ROUSCULP C L, BURTON D E. A compatible, energy and symmetry preserving Lagrangian hydrodynamics algorithm in three-dimensional Cartesian geometry[J]. Journal of Computational Physics, 2000, 157:89-119.
[7] GODUNOV S K. Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics[J]. Matematicheskii Sbornik, 1959, 47(3):271-306.
[8] ADDESSIO F L, CARROLL D E, DUKOWICZ J K, et al. CAVEAT:A computer code for fluid dynamics problems with large distortion and internal slip[R]. Los Alamos National Laboratory, UC-32, 1988.
[9] CHENG J, SHU C W. A high order ENO conservative Lagrangian type scheme for the compressible Euler equations[J]. Journal of Computational Physics, 2007, 227:1567-1596.
[10] DESPRÉS B, MAZERAN C. Lagrangian gas dynamics in two dimensions and Lagrangian systems[J]. Archive for Rational Mechanics and Analysis, 2005, 178:327-372.
[11] MAIRE P H. A cell-centered Lagrangian scheme for two-dimensional compressible flow problems[J]. SIAM Journal on Scientific Computing, 2007, 29(4):1781-1824.
[12] MAIRE P H, BREIL J. A second-order cell-centered Lagrangian scheme for two-dimensional compressible flow problems[J]. International Journal for Numerical Methods in Fluids, 2007(56):1417-1423.
[13] MAIRE P H. A high-order cell-centered Lagrangian scheme for two-dimensional compressible fluid flows on unstructured meshes[J]. Journal of Computational Physics, 2009, 228:2391-2425.
[14] MAIRE P H, LOUBERE R. Staggered Lagrangian discretization based on cell-centered Riemann solver and associated hydrodynamics scheme[J]. Communications of Computational Physics, 2011, 10(4):940-978.
[15] ROE P L. Discrete models for the numerical analysis of time-dependent multidimensional gas dynamics[J]. Journal of Computational Physics, 1986,63:458-476.
[16] MORTON K W, ROE P L. Vorticity-preserving Lax-Wendroff-type schemes for the system wave equation[J]. SIAM Journal on Scientific Computing, 2001, 23(1):170-192.
[17] LUNG T B, ROE P L. Toward a reduction of mesh imprinting[J]. International Journal for Numerical Methods in Fluids, 2014, 76:450-470.
[18] PRASAD P. Nonlinear hyperbolic waves in multi-dimensions[M]. New York:Chapman & Hall/CRC, 2001.
[19] LUKÁČOVÁ-MEDVID'OVÁ M, SAIBERTOVÁ J, WARNECKE G. Finite volume evolution Galerkin methods for nonlinear hyperbolic systems[J]. Journal of Computational Physics, 2002, 183:533-562.
[20] SUN Y T, REN Y X. The finite volume local evolution Galerkin method for solving the hyperbolic conservation laws[J]. Journal of Computational Physics, 2009, 228:4945-4960.
[21] SUN Y T, REN Y X, YU M, ZHANG S D. A finite volume method for 2D inviscous Lagrangian hydrodynamics based on characteristics theory[J]. Chinese Journal of Computational Physics, 2011, 28(1):19-26.
[22] SUN Y T, JIA Z P, YU M, REN Y X. A second order Lagrangian scheme based on characteristics theory for two-dimensional compressible flows[J]. Chinese Journal of Computational Physics, 2012, 29(6):791-798.
[23] MENIKOFF, RALPH. JWL equation of state[R]. LA-UR-15-29536, Los Alamos, 2015.
[24] TARVER C M, MCGUIRE E M. Reactive flow modeling of the interaction of TATB detonation waves with inert materials[C]. The 12th Symposium on Detonation, San Diego, 2002.
[25] CHERTOCK A, CUI S, KURGANOV A, WU T. Steady state and sign preserving semi-implicit Runge-Kutta method for ODES with stiff damping term[J]. SIAM Journal on Numerical Analysis, 2015, 53(4):2008-2029.
[26] MADER C L, ZUMBRO J D, FERM E N. Proton radiographic and numerical modeling of colliding, diverging PBX-9502 detonations[C]. The 12th Symposium on Detonation, San Diego, 2002.