Disclaimer: I have no experience with shockwaves, I've just done a fair amount of physics, and the key step of what follows (determining the energy of the shockwave) is just an educated estimate. See future work for details. I will update this based on feedback on assumptions I made.
Outline
First, calculate the kinetic energy produced by the adiabatic expansion. I'll do this by assuming adiabatic expansion, calculating the work done, and setting that equal to the kinetic energy.
Second, we need to determine how much of that kinetic energy goes into producing a shockwave. I think this will require reasonably involved hydrodynamics which I don't know, so I'm going to make estimates in this section (could use improvement).
Third, I assume the re-released air is just going to provide a steady medium for the shockwave to travel through, which means the shockwave travels normally, through an air of constant density and pressure.
Fourth, I'll calculate how the shockwave's pressure decays. I'll assume the thickness of the shockwave is on the order of 2 cm, and doesn't expand as the shockwave propagates.
Finally, I'll find a value for lethal, wounding, and building damage sudden pressure changes and figure out where exactly that threshold is for these explosives.
Kinetic energy calculation
So, to start, we want to calculate energy gained from adiabatic expansion, where PV^1.4 is constant. I'll assume the air in a 10 meter radius is effected, which means we'll have:
PV^1.4 = 101300 Pa * (4/3 * pi * 10m^3 - 4/3 * pi * 5m^3)^1.4 = 1.1268 * 10^10 J
Then, in the final situation, the volume is now just the 10m radius sphere, which means we can calculate the pressure as
P = 1.1269 * 10^10 J / V^1.4 = 1.1269 * 10^10 J / (4/3 * pi * 10m^3)^1.4 = 95702 Pa
Then, assuming we have the same amount of air at the end, and that ideal gas laws hold, this means that
P_i * V_i = nRT_i
and
P_f * V_f = nRT_f
which, we can divide the second equation by the first to get
T_f / T_i = P_f * V_f / (P_i * V_i)
Now, since we know T_i is 293 K, P_i is 101300 Pa, P_f is 95702 Pa, V_i is 3665 m^3, and V_f is 4188 m^3
T_f = (T_i * P_f * V_f) / (P_i * V_i)
T_f = (293K * 95702 Pa * 4188 m^3) / (101300 Pa * 3665 m^3)
T_f = 31 K
This means the change in temperature was T_f - T_i = 262 K
Internal energy is given by
U = cnT
where c is specific heat. Technically, c changes as temperature, pressure, and volume change and depends on the atomic composition for air, but these effects shouldn't change it's order of magnitude, and I don't think they'd vary it outside the range of 10-30 J/(mol*K), so I'll assume c = 20 J/(mol*K). At standard temperature and pressure, a mole of gas takes up 22.4 liters, or 0.0224 m^3. This means that n = 3665 m^3 / 0.0224 m^3/mole = 163616 moles.
thus, change in internal energy is given by:
U = cn*(T_f - T_i) = 20 J/(mol*K) * 163616 mol * (293K - 31 K) = 857 MJ
Which is also the work done by the compressing gas. Thus, the kinetic energy produced via the compressing gas is 857 MJ. First part complete.
Shockwave Energy
I'm sure there's all sorts of complicated hydrodynamics that goes into producing the shockwave. I'd imagine that the shock is produce by the collision of the air from each section expanding inwards to the center, then has some sort of effect as it passes outwards through the gas coming in. I don't know how to calculate that, I don't have the books to let me calculate that, and I don't trust my ability to teach myself it via the internet. What I do know, is that not all the kinetic energy produced by the implosion will contribute to a shockwave. Therefore, I'll assume that when all is said and done, there's a .5m radius sphere at the heart of the vacuum which causes the shockwave, that the kinetic energy is equally spread across the 10m radius sphere (obviously wrong, but I'm not sure of a better approximation), and that all the kinetic energy within the .5m radius sphere is converted into a shockwave. Obviously, I'm flying by the seat of my pants. If someone wants to offer suggestions on improving this section, I'm all ears.
So, the kinetic energy of the shockwave can be calculated as just:
K_shock = K_total * V_shock/V_total
Where the K_total was found earlier to be 857 MJ, the V_total is a 10 m radius sphere (so 4188 m^3), and V_shock is a .5 m radius sphere (so 4/3 * pi * 0.5^3 = 0.524 m^3), so
K_shock = 857 MJ * (0.524 m^3)/(4188 m^3) = 107 kJ
So the energy of the shockwave is 107 kJ
Propogation
So, the energy of 107 kJ will be spread out over the shock-front. The lethal thing is pressure in this case, and we can calculate pressure from kinetic energy as P = K/V. I'll assume that the shockwave expands as a thick surface to a sphere, and estimate the volume as SA * 2 cm. This means our equation for Pressure is
P = E/V
P = 107 kJ/(4 * pi * r^2 * 0.02 m) = 425739/r^2 Pa
Lethality
Based on
physics stack exchange, ~15 psi will cause lung damage, 40 psi will likely be fatal, and ~ 5 psi is likely to cause catastrophic building damage. Taking these as "injurous," "lethal," and "destruction" zones, and converting to Pa, we find:
103421 Pa --> injurous zone
275790 Pa --> lethal zone
34473.8 Pa --> destruction zone
Using our formula for pressure derived earlier, we find an injury sphere of:
r = sqrt(425739Pa * m^2/ 103421 Pa) = 2 meter radius
and a lethal sphere of
r = sqrt(425739Pa * m^2/ 275790 Pa) = 1.25 meter radius
and a building destruction sphere of
r = sqrt (425739Pa * m^2/ 34473.8 Pa) = 3.5 meter radius
Conclusion
So, based on my best calculations, a 5 meter vacuum seal which releases the air back to normal when destroyed will produce a shockwave with a
1.25 meter kill radius and a
2 meter substantial injury radius and a 3.5 meter building destruction radius.
Future Work
There's a lot of not wonderful assumptions I made here. Foremost among them was my calculation for shockwave production. Ultimately, my work says that only 0.0125 % of the energy from the air rushing into the seal is released as a shockwave. If someone knows a better method for calculating the conversion from general kinetic energy to shockwave, I'd be happy to hear it, and the rest of my work should be easy to update.
One reasonable way to calculate the energy would be to estimate a linear distribution of energy, with maximum kinetic energy at the epicenter, and zero kinetic energy at the edges. I think that would increase the kinetic energy to ~0.5%, but I'd need to solve the integral to actually know which.... I'll probably do soon.
For instance, if we assume that 1% of the kinetic energy is converted into a shockwave, suddenly we have a pressure formula resulting in
P = 34098946 Pa * m^2 / r^2
which results in a wounding radius of
r = sqrt(34098946 Pa * m^2 / 103421 Pa) = 18 m
and a kill radius of
r = sqrt(34098946 Pa * m^2 / 275790 Pa) = 11 m
and a destruction radius of
r = sqrt(34098946 Pa * m^2 / 34473.8 Pa) = 31.5 m
If someone knows more than me about shockwave production or other relevant physics, or finds math errors, tag me and I'll update this.