A numerical scheme for strong blast wave driven by explosion

After the detonation of a solid high explosive, the material has extremely high pressure keeping the solid density and expands rapidly driving strong shock wave. In order to simulate this blast wave, a stable and accurate numerical scheme is required due to large density and pressure changes in time and space. The compressible fluid equations are solved by a fractional step procedure which consists of the advection phase and non-advection phase. The former employs the Rational function CIP scheme in order to preserve monotone signals, and the latter is solved by interpolated differential operator scheme for achieving the accurate calculation. The procedure is categorized into the fractionally stepped semi-Lagrangian. The accuracy of our scheme is confirmed by checking the one-dimensional plane shock tube problem with 103 times initial density and pressure jump in comparison with the analytic solution. The Sedov–Taylor blast wave problem is also examined in the two-dimensional cylindrical coordinate in order to check the spherical symmetry and the convergence rates. Two- and three-dimensional simulations for the blast waves from the explosion in the underground magazine are carried out. It is found that the numerical results show quantitatively good agreement with the experimental data. Copyright © 2006 John Wiley & Sons, Ltd.

[1]  Feng Xiao,et al.  Constructing oscillation preventing scheme for advection equation by rational function , 1996 .

[2]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[3]  Dimitris Drikakis,et al.  Simulation of blast wave propagation over a cylinder , 1997 .

[4]  A. Harten.  ENO schemes with subcell resolution , 1989 .

[5]  Takayuki Aoki,et al.  Interpolated differential operator (IDO) scheme for solving partial differential equations , 1997 .

[6]  Shen-Min Liang,et al.  Numerical study of spherical blast-wave propagation and reflection , 2002 .

[7]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[8]  Jay P. Boris,et al.  Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works , 1973 .

[9]  B. M. Dobratz,et al.  Properties of chemical explosives and explosive simulants , 1972 .

[10]  D. Oron,et al.  Study of the membrane effect on turbulent mixing measurements in shock tubes , 2000 .

[11]  Huanan Yang,et al.  An artificial compression method for ENO schemes - The slope modification method. [essentially nonoscillatory , 1990 .

[12]  Alain Merlen,et al.  Cut cell strategy for 3-D blast waves numerical simulations , 2001 .

[13]  Charles L. Mader,et al.  Shock and Hot Spot Initiation of Homogeneous Explosives , 1963 .

[14]  Feng Xiao,et al.  Constructing a multi-dimensional oscillation preventing scheme for the advection equation by a rational function , 1996 .

[15]  W. E. Baker.  Explosions in air , 1973 .

[16]  H. Brode.  Numerical Solutions of Spherical Blast Waves , 1955 .

[17]  Takashi Yabe,et al.  A universal solver for hyperbolic equations by cubic-polynomial interpolation I. One-dimensional solver , 1991 .

[18]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[19]  Kenji Takizawa,et al.  Higher-order schemes with CIP method and adaptive Soroban grid towards mesh-free scheme , 2004 .

[20]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[21]  L. Sedov.  Similarity and Dimensional Methods in Mechanics , 1960 .

[22]  Mark L. Wilkins,et al.  Use of artificial viscosity in multidimensional fluid dynamic calculations , 1980 .

[23]  T. Aoki.  MULTI-DIMENSIONAL ADVECTION OF CIP (CUBIC-INTERPOLATED PROPAGATION) SCHEME , 1995 .

[24]  Takayuki Aoki,et al.  Improvement of accuracy and stability in numerically solving hyperbolic equations by IDO (Interpolated Differential Operator) scheme with Runge-Kutta time integration , 2004 .

[25]  G. Taylor.  The formation of a blast wave by a very intense explosion I. Theoretical discussion , 1950, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[26]  Sukumar Chakravarthy,et al.  High Resolution Schemes and the Entropy Condition , 1984 .

[27]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[28]  H. Brode.  Blast Wave from a Spherical Charge , 1959 .

[29]  T. Yabe,et al.  The constrained interpolation profile method for multiphase analysis , 2001 .