Past exam of the mathematics course of the University of Cambridge / 2013 / iii / Paper 52 / 2

Let $e=\rho u^2/2+p/(\gamma -1)$ be the total gas energy density, with $\gamma >1$. The one-dimensional continuity and momentum equations areUsing continuity to put the momentum equation into conservative form givesThe internal energy density $e_{\mathrm{int}}=p/(\gamma -1)$ satisfiesby the adiabatic pressure equation. Multiplying the momentum equation by $u$ and using continuity givesAdding these equations combines the pressure work into $-{\partial }_z(pu)$. The conservation law variables and conservation law fluxes are thereforeThe three rows express conservation of mass, momentum and total energy. Their integral conservation laws remain meaningful across a shock wave, where the differential pressure and velocity equations cannot be applied pointwise.

Integrating each conservation law across a vanishingly thin interval around the stationary shock wave leaves equal conservation law fluxes on its two sides. HenceThese are the Rankine-Hugoniot conditions for a perfect gas. Signed velocities may both be negative when the material travels from positive to negative $z$; no sign change is needed in the conservation law fluxes.

Put $r={\rho }_2/{\rho }_1=u_1/u_2$, $P=p_2/p_1$ and $m=M_1^2={\rho }_1{u}_1^2/(\gamma p_1)$. Momentum conservation givesDivide the energy condition by the nonzero mass flux. Equality of kinetic energy plus specific enthalpy givesSubstituting $P$ and multiplying by $2(\gamma -1)r^2/\gamma$ yieldsThe factor $r-1$ is the continuous, no-shock solution. On the nontrivial normal shock wave branch,An admissible compressive gas shock wave has $M_1>1$, so $r>1$ and $P>1$. The algebraic jump equations alone also allow a reversed expansive discontinuity; the entropy production in a perfect-gas shock excludes that branch. At $M_1=1$ the nontrivial branch joins the continuous solution.

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)$.