Numerical simulation of an explosion in the air using a high-order approximation difference scheme
Authors: Menshakov S.S., Taran V.A.
Published in issue: #5(161)/2025
Category: Mechanics | Chapter: Mechanics of Liquid, Gas, and Plasma
The paper considers numerical solution to the problem of an explosion of the TNT charge with the 0.04 m radius. It applies a difference scheme of the higher approximation order to an explosion of the condensed explosive charge in the air. Solution is provided in the flat 2D formulation, where the square computational domain size is equal to 6 m. Preliminary testing of the numerical algorithm was carried out on the specialized test problems reproducing shock jumps, rarefaction waves, contact discontinuities and vortex structures in the two-dimensional domains. Further, they were compared with the known results. The paper shows that the developed algorithm of higher approximation order provides advantages in detailing the shock-wave pattern of the media flow.
EDN LRPRRI
References
[1] Gelfand B.E., Silnikov M.V. Gazovye vzryvy [Gas explosions]. St. Petersburg, Asterion Publ., 2007, 240 p.
[2] Orlenko L.P., ed. Fizika vzryva [Physics of Explosion]. 3rd ed., rev. In 2 vols. Vol. 1. Moscow, Fizmatlit Publ., 2004, 832 p.
[3] Zhdan S.A. Raschet vzryva gazovogo sfericheskogo zaryada v vozdukhe [Calculation of the spherical gas charge explosion in the air]. PMTF — Journal of Applied Mechanics and Technical Physics, 1976, no. 6, pp. 69–74.
[4] Selivanov V.V., Okhitin V.N. Matematicheskoe modelirovanie sluchaynykh gazovykh vzryvov [Mathematical modeling of accidental gas explosions]. FGV — Combustion, Explosion and Shock Waves, 1995, vol. 31, no. 6, pp. 155–161.
[5] Gorev V.A. Sravnenie vozdushnykh voln ot raznykh istochnikov [Comparison of the air explosion waves from different sources]. FGV — Combustion, Explosion and Shock Waves, 1982, vol. 18, no. 1, pp. 94–101.
[6] Sadin D.V., Golikov I.O., Shirokova E.N. Testirovanie gibridnogo metoda krupnykh chastits na dvukhmernykh zadachakh Rimana [Testing the hybrid large-particle method using two-dimensional Riemann problems]. Nauchno-tekhnicheskie vedomosti SPbGPU. Fiziko-matematicheskie nauki — St. Petersburg Polytechnic University Journal: Physics and Mathematics, 2021, vol. 14 (1), pp. 58–71.
[7] Godunov S.K. Raznostnyi metod chislennogo rascheta razryvnykh resheniy gidromekhaniki [Difference method in numerical calculation of discontinuous solutions of the fluid mechanics]. Matematicheskiy sbornik — Sbornik: Mathematics, 1959, vol. 47 (89), no. 3, pp. 271–306.
[8] Liska R., Wendroff B. Comparison of several difference schemes on 1D and 2D test problems for the Euler equations. SIAM Journal on Scientific Computing, 2003, vol. 25, no. 3, pp. 995–1017.
[9] Shi J., Zhang Y.-T., Shu C.-W. Resolution of high order WENO schemes for complicated flow structures. Journal of Computational Physics, 2003, vol. 186, no. 2, pp. 690–696.
[10] Kermani M.J., Gerber A.G., Stockie J.M. Thermodynamically based moisture prediction using Roe’s scheme. In: The 4th Conference of Iranian AeroSpace Society, Amir Kabir University of Technology, Tehran, Iran, 2003, January 27–29.
[11] Titarev V.A., Toro E.F. Finite-volume Weno schemes for three-dimensional conservation laws. Journal of Computational Physics, 2004, vol. 201, pp. 238–260.
[12] Shu C.W. Total-variation-diminishing time discretization’s. SIAM Journal on Scientific and Statistical Computation, 1988, vol. 9, pp. 1073–1084.
[13] Berger M.J., Oliger J. Adaptive mesh refinement for hyperbolic partial differential equations. Journal of Computational Physics, 1984, vol. 53 (3), pp. 484–512.
[14] The HLLC Riemann Solver. Eleuterio TORO. Laboratory of Applied Mathematics. University of Trento, Italy. toro@ing.unitn.it, http://www.ing.unitn.it/toro, August 26, 2012.
[15] Kolpakov V.I., Rubtsov A.A., Ladov S.V. Matematicheskoe modelirovanie kumulyativnykh zaryadov [Mathematical modeling of cumulative charges]. Moscow, BMSTU Publ., 1998, 36 p.