Solution (source code)

= Solution

The <Rankine-Hugoniot conditions for a perfect gas> express conservation of mass flux, normal momentum flux, and total energy flux through a steady planar <shock wave> in its rest frame. With mass flux $j=\rho_1u_1=\rho_2u_2$, momentum flux is $\rho u^2+p$ and specific transported energy is $u^2/2+e+p/\rho$. The last term is pressure work, so the energy equation contains <specific enthalpy> rather than only <specific internal energy>. Neglect gravity across the thin front, magnetic stresses and radiative energy loss for the adiabatic case.

For a calorically perfect <ideal gas>, $e=p/[(\gamma-1)\rho]$ and $c_s^2=\gamma p/\rho$. Therefore
$$
\boxed{\frac{u_1^2}{2}+\frac{c_{s,1}^2}{\gamma-1}
=\frac{u_2^2}{2}+\frac{c_{s,2}^2}{\gamma-1}.}
$$
Put $\chi=\rho_2/\rho_1$, so $u_2=u_1/\chi$. Momentum conservation gives
$$
\frac{p_2}{p_1}=1+\gamma\mathcal M_1^2\left(1-\frac1\chi\right),
\qquad \mathcal M_1=\frac{u_1}{c_{s,1}}.
$$
Insert this and $c_{s,2}^2/c_{s,1}^2=(p_2/p_1)/\chi$ into the energy equation. After multiplying out, its two factors are
$$
(\chi-1)\left[\big((\gamma-1)\mathcal M_1^2+2\big)\chi
-(\gamma+1)\mathcal M_1^2\right]=0.
$$
The first root is the unchanged flow. The compressive <shock> branch is
$$
\boxed{\chi=\frac{(\gamma+1)\mathcal M_1^2}{(\gamma-1)\mathcal M_1^2+2}
=\frac{\gamma+1}{\gamma-1+2/\mathcal M_1^2}.}
$$
It has $\chi>1$ for upstream supersonic flow. A strong shock approaches $(\gamma+1)/(\gamma-1)$: \b[four for a monatomic gas], or six for a fixed rotationally excited diatomic gas with $\gamma=7/5$. The assumed fixed heat capacity ceases to describe the gas when dissociation or ionization becomes important. Adiabatic here means no net heat lost from the shock-plus-flow energy budget; the irreversible front still increases <entropy>.

For an <isothermal shock> of unchanged composition, let $a^2=p/\rho=k_BT/\bar m$ be the common <isothermal sound speed> squared. Mass and momentum conservation give
$$
u_1^2+a^2=\frac{u_1^2}{\chi}+a^2\chi,
\qquad (\chi-1)(a^2\chi-u_1^2)=0.
$$
The compressive branch is therefore
$$
\boxed{\chi=\mathcal M_{T,1}^2,\qquad \mathcal M_{T,1}=u_1/a.}
$$
Unlike the adiabatic ratio, it has no fixed upper bound as its upstream Mach number grows. Heat removal prevents the pressure increase that would otherwise resist compression. This approximation describes an efficiently cooled jump between equal-temperature far states, not an adiabatic shock with its energy equation left unchanged. For the same upstream physical temperature and velocity, $\mathcal M_{T,1}=\sqrt\gamma\mathcal M_1$; the two sound-speed conventions must not be confused.

\Image[/past-exam-of-the-astrophysical-fluid-dynamics-course-of-the-university-of-cambridge/2015/paper-1-shock-compression.png]
{title=Bounded adiabatic compression and increasing isothermal compression}
{description=Original comparison for fixed composition. The horizontal Mach coordinate uses each model's own upstream sound speed, so equal plotted Mach numbers do not mean identical upstream temperatures and velocities.}

Full dissociation changes the equation of state even when the temperature is held constant. With negligible electron and binding-energy mass corrections, upstream molecular <hydrogen> has $a_1^2=k_BT/(2m_H)$ and downstream atomic <hydrogen> has $a_2^2=k_BT/m_H=2a_1^2$. The conserved quantity is mass flux; particle number flux doubles. Thus
$$
\boxed{\chi=\frac{u_1}{u_2},\qquad
u_1^2+a_1^2=u_1u_2+2a_1^2\frac{u_1}{u_2}.}
$$
Equivalently, the <dissociating isothermal shock> relation between temperature and the two velocities is
$$
\boxed{a_1^2=\frac{u_1u_2(u_1-u_2)}{2u_1-u_2}.}
$$
Writing $\mathcal M_{T,1}=u_1/a_1$ gives
$$
2\chi^2-(\mathcal M_{T,1}^2+1)\chi+\mathcal M_{T,1}^2=0.
$$
The high-compression branch is
$$
\boxed{\chi=\frac{\mathcal M_{T,1}^2+1+
\sqrt{(\mathcal M_{T,1}^2+1)^2-8\mathcal M_{T,1}^2}}4.}
$$
Its strong-shock limit is $\chi\sim\mathcal M_{T,1}^2/2$, rather than $\mathcal M_{T,1}^2$. The pressure per unit mass has doubled, reducing the compression available at a fixed upstream temperature and speed. A real compressive branch in this prescribed fully dissociated, equal-temperature model requires $\mathcal M_{T,1}\ge1+\sqrt2$; the high root has downstream speed no greater than the downstream isothermal sound speed. These mass and momentum conditions do not themselves establish that sufficient energy is available to dissociate every molecule. The chemical energy and any required radiative loss must also be included in a physical energy budget.