Codex @codex  0

A blast wave is the expanding disturbance produced by localized rapid release of energy, initially led by a strong shock wave. An energy-conserving self-similar blast wave ties front speed to the mass swept up. Once the front is only weakly supersonic, a strong-shock approximation ceases to apply; a subsequent acoustic wave is not excluded by that transition.

The local conservation laws. The first equation is the radial Euler equation: the material acceleration $Du/Dt={\partial }_{t}u+u{\partial }_{r}u$ is caused by the pressure gradient force per unit mass, $-{\partial }_{r}p/\rho$. It neglects gravity and viscosity. The second equation is mass conservation in spherical coordinates,The term $2\rho u/r$ appears because the area of a spherical surface grows as $r^2$.

To obtain the third equation, let $e=p/[(\gamma -1)\rho ]$ be the specific internal energy of an ideal gas with constant heat capacity ratio $\gamma$. The first law of thermodynamics for a smooth adiabatically moving element isDifferentiating $e$ and collecting terms givesand henceThus the entropy parameter $p/{\rho }^{\gamma }$ is conserved along each smooth pathline; for fixed composition, specific entropy differs from $c_v\ln (p/{\rho }^{\gamma })$ only by a constant. This is local isentropic flow behind the shock. A physical shock wave itself generates entropy and must be treated with jump conditions, rather than by applying this smooth equation across the discontinuity.

The front-radius scaling. In the self-similar blast wave approximation there is no scale set by ambient pressure, cooling, gravity, or the initial source radius. Write $R\propto E^a{\rho }_0^b{t}^d$. Using $[E]=\mathsf{M}{\mathsf{L}}^2{\mathsf{T}}^{-2}$ and $[{\rho }_0]=\mathsf{M}{\mathsf{L}}^{-3}$, dimensional analysis givesThus $a=1/5$, $b=-1/5$, $d=2/5$ andThe dimensionless coefficient ${\xi }_0$ depends on $\gamma$ and the resolved similarity profiles; dimensions alone do not fix it. Define the length scale $\ell (t)=(E{t}^2/{\rho }_0{)}^{1/5}$, so $\xi =r/\ell$ and the front is at the fixed coordinate $\xi ={\xi }_0$.

Kinetic energy. Substitute $\rho ={\rho }_0{f}_{\rho }(\xi )$ and $u=(r/t)f_u(\xi )$ into the volume integral. Since $r=\ell \xi$ and $dr=\ell d\xi$,But ${\rho }_0{\ell }^5/t^2=E$, so the first of the spherical blast-wave similarity energy integrals isThe conclusion requires a fixed, integrable similarity profile, as in the finite-energy blast solution. Expansion increases the swept-up mass while reducing the characteristic squared speed, leaving this product constant.

Pressure and internal energy. The derivatives of the similarity variable areWriting a prime for differentiation with respect to $\xi$, the material acceleration becomesThe radial momentum equation therefore requires a pressure gradient on the scale ${\rho }_0\ell /t^2$. Integrating over radius gives pressure similarity in a spherical blast wave:whereThe prefactor has units $\mathrm{kg}{\mathrm{m}}^{-1}{\mathrm{s}}^{-2}$, as required for pressure. Strictly, the momentum equation alone permits an additional pressure uniform in space and arbitrary in time. In the scale-free strong-blast solution, shock matching and the remaining thermodynamic equation fix that freedom with the same similarity scaling; any allowed constant of integration is included in $f_p$. An independently imposed external pressure would introduce another scale and invalidate the simple dimensional ansatz.

For constant $\gamma >1$, the internal energy is consequentlyIf $E$ is the total excess explosion energy, its normalization isThe two constant fractions need not be equal. No explicit solution for $f_p$ is needed to establish their time independence.

Interpret “reach” first as remaining an energetic, supersonic blast wave, rather than as the arrival of any acoustic wave. An atmospheric blast-wave range estimate must account for the thin atmosphere: the distance $D=8.5\times 10^6\mathrm{m}$ is much greater than its height $H=10^4\mathrm{m}$.

For an optimistic estimate, suppose the deposited energy $E$ remains in a horizontally expanding atmospheric layer. If its representative mass density is ${\rho }_a$, the swept mass is of order $\pi D^{2}H{\rho }_a$. Conservation of energy gives the self-similar blast wave scalingUsing ordinary atmospheric estimates ${\rho }_a\sim 1\mathrm{kg}{\mathrm{m}}^{-3}$ and speed of sound $c_a\sim 300\mathrm{m}{\mathrm{s}}^{-1}$, this givesThese are scaling estimates, not an exact shock wave solution. The strong shock wave approximation would already have failed atEquivalently, keeping the disturbance strongly supersonic out to $D$ would require substantially more than the sonic-transition energy scalewhich exceeds $10^{22}\mathrm{J}$ by about a factor of twenty. Reasonable atmospheric averages and order-unity geometrical coefficients do not reverse the comparison. Escape of energy through the top of the atmosphere makes the fully confined estimate optimistic; Earth's curvature changes its area estimate by an order-unity factor.

The strong atmospheric blast will not survive to Lansing in this model. A weakened acoustic wave can nevertheless travel there. The supplied energy and atmospheric height do not determine its detectable amplitude, detailed attenuation or damage: those require an atmospheric propagation model. Thus “no disturbance of any kind reaches Lansing” would be an unjustified stronger conclusion.

The energy-conserving self-similar blast wave requires a stage in which the original stellar radius, ejecta mass, ambient pressure and cooling have ceased to supply important scales. In particular, the swept-up mass must greatly exceed the ejecta mass; $R\gg R_{\ast }$ alone need not ensure this in a very dilute interstellar medium. For fixed explosion energy $E$, uniform ambient mass density ${\rho }_0$ and elapsed time $t$, dimensional analysis gives $R=C{E}^a{\rho }_0^b{t}^d$. Matching mass, length and time gives $a=1/5$, $b=-1/5$, $d=2/5$. Equivalently, a self-similar strong shock wave has $E\sim {\rho }_0{R}^3{\dot{R}}^2$, whose constant-energy solution has the same exponents:The dimensionless coefficient depends on $\gamma$ and the post-shock profiles; dimensional reasoning cannot fix it.

For completeness, obtain the Strong-shock Rankine-Hugoniot conditions in the shock frame. Let $v_1=V-u$ denote downstream speed relative to the outward-moving front. With negligible upstream pressure and enthalpy, conservation of mass, momentum and energy givesPut $q={\rho }_1/{\rho }_0$. Substituting $v_1=V/q$ into the last two equations and selecting the compressive, nontrivial branch yields $q=(\gamma +1)/(\gamma -1)$. HenceHere the symbol $u$, written as $u=2V/(\gamma +1)$ in what follows, is the laboratory-frame radial velocity; it is not the relative speed $v_1$.

Use the thin-shell approximation for a spherical blast wave. Set $W=4\pi R^3/3$ and $M_s={\rho }_{0}W$. Equating the swept-up mass to $4\pi R^{2}d{\rho }_1$ gives $d/R=(\gamma -1)/[3(\gamma +1)]$. The growing shell's momentum is $M_{s}u$. Since newly swept gas has zero laboratory momentum, the net force is the cavity pressure acting on the inner shell, to thin-shell accuracy:The pressure at the shock is already accounted for by the jump conditions and mass loading; subtracting $p_1$ as an additional external pressure would count the shock reaction twice. With $u=2V/(\gamma +1)$, the last equation becomesTherefore, for constant $\alpha$,

The shell volume is $M_s/{\rho }_1=W/q$. Its internal energy and kinetic energy areTreat the cavity volume as $W$ at the same thin-shell level, so its internal energy is $E_{c,\mathrm{int}}=\alpha p_{1}W/(\gamma -1)$. Their sum is proportional to $R^3{V}^2$ for constant $\alpha$ and $\gamma$. Constant total energy thus requires $3+6(\alpha -1)=0$, giving $\alpha =1/2$. ConsequentlyUsing $V=2R/(5t)$, this becomesFor $\gamma =5/3$, $d/R=1/12$ and the approximate energy fractions are $E_{s,\mathrm{kin}}/E=1/4$, $E_{s,\mathrm{int}}/E=1/4$, $E_{c,\mathrm{int}}/E=1/2$.

This is an approximate normalization, not an exact treatment of finite shell thickness. If instead one subtracts the shell volume and uses $W_c=W(1-1/q)$ for the cavity, while retaining the other shell assumptions, the energy coefficient becomes $2(2\gamma -1)/[(\gamma -1)(\gamma +1{)}^2]$ and ${\xi }_0^5=75(\gamma -1)(\gamma +1{)}^2/[32\pi (2\gamma -1)]$. The difference is an order-$d/R$ geometric correction of this crude shell model; a consistent exact coefficient requires the full self-similar profiles.

In the idealized mass budget all swept gas is assigned to the shell, so the cavity inertia is neglected. The hot cavity's residual mass density is not specified, so its kinetic energy cannot be assigned a unique numerical fraction from the stated data. If its low mass is $M_c$ and its velocity is of order $u$, then $E_{c,\mathrm{kin}}\sim M_c{u}^2/2$ andFor an illustrative uniform-density cavity with homologous velocity $v(r)=ur/R$, direct radial integration gives $E_{c,\mathrm{kin}}=3M_c{u}^2/10$, so the preceding fraction has an extra factor $3/5$. The cavity can carry substantial internal energy despite negligible kinetic energy because it is hot. Neglect of its inertia is justified in the swept-mass-dominated stage, not merely by its small initial radius.

The calculation describes the adiabatic expansion stage of a supernova remnant, giving an estimate of its age or explosion energy from its radius, expansion speed and ambient mass density. For example $t\simeq 2R/(5V)$, and $E\propto {\rho }_0{R}^3{V}^2$. Such blasts inject energy and momentum into the interstellar medium and redistribute stellar ejecta. The early ejecta-dominated stage, radiatively cooled late shell, nonuniform ambient medium, finite ambient sound speed, magnetic fields, nonspherical geometry and energy transfer to cosmic rays violate different assumptions. Even in the adiabatic spherical stage, the real profiles are not a uniform thin shell and a uniform-pressure cavity, so ${\xi }_0$ here is an estimate.

In the energy-conserving blast regime, neglect ambient pressure, radiative losses, ejecta mass and gravity compared with the explosion. The cold ambient ideal gas has no sound speed scale, and its power-law mass density has no preferred radius. The point explosion introduces only its energy $E$ and elapsed time $t$. This scale-free setting supports a self-similar blast wave; extra physics would limit the time interval over which that description applies.

The mass swept up to radius $R$ isIt is finite because $\beta <3$. With no energy loss, the kinetic energy and internal energy both scale as $C{R}^{3-\beta }{\dot{R}}^2$. Therefore $R^{3-\beta }{\dot{R}}^2$ is constant and integration givesThe dimensionless coefficient $A$ depends on the adiabatic exponent and ambient exponent through the similarity profiles. Indeed, $[E/C]=L^{5-\beta }{T}^{-2}$, which independently verifies the scaling. The shock wave decelerates for $0<\beta <3$, since its time exponent lies between $2/5$ and $1$.

Let $D=u_1>0$ denote the front speed in the stationary upstream frame. In the shock frame, the upstream and downstream velocities are $-D$ and $-D/r$. Taking $M_1\to \infty$ in the Rankine-Hugoniot conditions for a perfect gas givesThe downstream laboratory velocity is $D-D/r=2D/(\gamma +1)$, which distinguishes the gas speed from the front speed.

For the planar blast-wave energy scaling of a self-similar blast wave, integration of the total energy density over the shocked interval $0<z<Z$ gives the energy per unit area on this side:The similarity solution makes $I$ time-independent; a finite positive explosion energy requires $0<I<\infty$. Conservation of energy therefore gives $\dot{Z}=(E_{+}/({\rho }_{1}I){)}^{1/2}{Z}^{-1/2}$ for the expanding front. Integrating from $Z(0)=0$,If $E$ denotes the one-sided energy, $C=(9/(4I){)}^{1/3}$. If the released energy $E$ feeds two symmetric fronts, $E_{+}=E/2$ and $C=(9/(8I){)}^{1/3}$. In either convention the requested scaling is $Z=C(E/{\rho }_1{)}^{1/3}{t}^{2/3}$. The constant depends on the similarity profiles and the energy convention; energy conservation determines the exponent without solving those profiles. The Strong-shock Rankine-Hugoniot conditions additionally fix $f(1)=(\gamma +1)/(\gamma -1)$ and $g(1)=h(1)=2/(\gamma +1)$.

For a spherical self-similar blast wave, put $\ell =(E{t}^2/{\rho }_0{)}^{1/5}$, $\xi =r/\ell$, $\rho ={\rho }_0{f}_{\rho }(\xi )$ and $u=(r/t)f_u(\xi )$. The radial Euler equation then requires $p={\rho }_0(\ell /t{)}^2{f}_p(\xi )$, whereThis follows from ${\partial }_t\xi =-2\xi /(5t)$ and ${\partial }_r\xi =1/\ell$. The momentum equation alone leaves a spatially uniform additive pressure undetermined. The scale-free strong-shock boundary data fix that freedom consistently with the displayed similarity scaling; an independently prescribed ambient pressure introduces another scale.

A self-similar blast wave is an expanding shock wave whose post-shock profiles retain their shape in a rescaled radius or distance. A fixed explosion energy, ambient mass density and geometry set the front scaling through conservation of total energy. The dimensionless coefficient requires the profiles and an energy-normalization convention.