The strong-shock, point-source solution and spherical isothermal distributions were used as initial conditions for a numerical integration of the differential equations of gas motion in Lagrangean form. The Von Neumann-Richtayer artificial viscosity was employed to avoid shock discontinuities. The solutions were carried from two thousand atmospheres to less than one-tenth atmospheres peak overpressure. Results include overpressure, density, partieli velocity, and position as functions of time and space. The dynamic pressure, the positive and negative impulses of both dynamic pressure and static over-pressure, positive and negative durations of pressure and velocity, and shock value of all quantities are also described for various times and radical distances. Analytical approximations to the numerical results are provided.
Brode, Harold L., Numberical Solutions of Spherical Blast Waves. Santa Monica, CA: RAND Corporation, 1954. https://www.rand.org/pubs/research_memoranda/RM1363.html.
Brode, Harold L., Numberical Solutions of Spherical Blast Waves, Santa Monica, Calif.: RAND Corporation, RM-1363-AEC, 1954. As of January 13, 2022: https://www.rand.org/pubs/research_memoranda/RM1363.html