 Spherical blast wave in a power-law ambient density

For an energy-conserving strong spherical shock wave expanding into mass density $C{r}^{-\beta }$, $0\le \beta <3$, dimensional analysis gives the displayed radius law. Neglect ambient pressure, radiation losses and gravity in this blast-dominated approximation. Write $u=\dot{R}U(\xi )$, $\rho =C{R}^{-\beta }D(\xi )$ and $p=C{R}^{-\beta }{\dot{R}}^{2}P(\xi )$, with $\xi =r/R$. The radial continuity equation, momentum equation and adiabatic pressure equation reduce toThe Strong-shock Rankine-Hugoniot conditions give $U(1)=P(1)=2/(\gamma +1)$ and $D(1)=(\gamma +1)/(\gamma -1)$. The explosion energy fixes the radius coefficient through $E=4\pi C{R}^{3-\beta }{\dot{R}}^2{\int }_0^1[D{U}^2/2+P/(\gamma -1)]{\xi }^{2}d\xi$.