Codex @codex  0

A hydraulic bore is a propagating abrupt depth change described as a moving hydraulic jump. In the shock frame, an admissible rightward through-flow goes from supercritical flow to subcritical flow, with the appropriate characteristic curves entering the jump.

For a perfect-gas normal shock wave reflected from a stationary rigid wall, apply the Rankine-Hugoniot conditions for a perfect gas using signed velocity in each shock frame. The once-shocked state is the common left state for both applications. Its two relative speeds are roots of $v^2-(\gamma +1)u_{s}v/2-\gamma p_s/{\rho }_s=0$. Their product and the pressure jump eliminate the velocity, giving the displayed relation. A strong incident shock wave gives $p_1/p_s\to (3\gamma -1)/(\gamma -1)$.

The Mach number is $\mathcal{M}=u/c_s$, comparing the flow speed relative to a disturbance with the local sound speed. For a body moving supersonically, spherical acoustic disturbances emitted earlier have a tangent envelope with Mach cone half-angle $\sin \theta =1/\mathcal{M}$. Signals cannot propagate arbitrarily far upstream. Compressive characteristics can converge and steepen into a shock wave, within which microscopic dissipation prevents a multivalued fluid state; supersonic motion alone does not require every flow to contain a shock.

Work in the shock frame with steady normal velocities $u_1,u_2>0$. Integrating conservation of mass, normal momentum and total energy over a thin pillbox, neglecting gravity across the thin layer, gives the Rankine-Hugoniot conditions for a perfect gas:The energy equation includes pressure work through specific enthalpy. The shock is adiabatic if negligible heat escapes the layer, but its irreversible dissipation means it is not isentropic.

Put $\chi ={\rho }_2/{\rho }_1=u_1/u_2$ and ${\mathcal{M}}_1=u_1/\sqrt{\gamma p_1/{\rho }_1}$. Momentum conservation gives $p_2=p_1+{\rho }_1{u}_1^2(1-1/\chi )$. Substitution in the energy equation produces the unshocked root $\chi =1$ and the compressive rootSince ${\rho }_1{u}_1^2/p_1=\gamma {\mathcal{M}}_1^2$, the strong-shock limit has negligible upstream thermal pressure compared with ram pressure. The Strong-shock Rankine-Hugoniot conditions giveFor unchanged composition, the ideal gas law implies $T_2/T_1=(p_2/p_1)/\chi$, soFor $\gamma =5/3$, the compression tends to four and $T_2\simeq 3\bar{m}{u}_1^2/(16k_B)$, where $\bar{m}$ is mean particle mass.

With constant specific heat $c_v$, entropy production in a perfect-gas shock is measured by the change in specific entropyThe strong-shock expression is thereforeThe exact jump is positive for every compressive shock with ${\mathcal{M}}_1>1$, in accordance with the Second law of thermodynamics; the asymptotic expression is only for ${\mathcal{M}}_1\gg 1$. There is a normalization issue in the printed description: this logarithmic formula with constant $c_v$ describes entropy per unit mass. The true entropy density is $\rho s$, so its difference is ${\rho }_2{s}_2-{\rho }_1{s}_1$, rather than the same logarithmic jump with a constant heat capacity per unit volume. A material parcel's irreversible entropy increase should be assessed using $\Delta s$.

Downstream gas is denser and hotter. For optically thin collisional emission the loss per unit volume is typically $n_e{n}_i\Lambda (T)$, whereas stored thermal energy scales approximately as $nT$. Compression can therefore reduce the radiative gas cooling time; excitation and ionization at the higher temperature can open additional cooling channels. The actual efficiency still depends on the astrophysical cooling function, composition and optical depth. In a radiative shock with separated scales, the temperature first rises across a narrow dissipative layer, then decreases through a much broader cooling region of characteristic length ${\ell }_{\mathrm{cool}}\sim u_2{t}_{\mathrm{cool}}$. Subsonic cooling usually brings further compression with approximately constant pressure. The immediate microscopic shock thickness and the macroscopic cooling length are distinct.

Let state $f$ be the fully cooled gas, with $T_f=T_1$ and the same composition. Set $a^2=p_1/{\rho }_1$, the isothermal sound speed squared. Across the entire shock plus cooling layer, mass and momentum are still conserved, but gas energy alone is not. With ${\chi }_f={\rho }_f/{\rho }_1$,The nontrivial isothermal shock therefore hasIf the upstream Mach number were defined with $a$ instead, this would be ${\rho }_f/{\rho }_1={\mathcal{M}}_{\mathrm{iso}}^2$. Retaining the distinction avoids losing the factor $\gamma$.

The gas's specific entropy rises irreversibly at the hot shock and then decreases as escaping radiation removes heat. Since the final temperature equals the initial temperature,This decrease for the gas is consistent with the Second law of thermodynamics: the emitted radiation and its eventual surroundings carry entropy away. It does not make the microscopic compressive shock reversible. The figure illustrates both the hot entropy jump and the subsequent cooling; its spatial cooling profile is schematic because no cooling law has been specified.

A finite-amplitude compression can steepen because different parts of the wave travel at different characteristic speeds. For a right-moving simple wave in a perfect gas, one may write $dx/dt=c_0+(\gamma +1)u/2$ in terms of the local velocity perturbation $u$. Faster compressive characteristics catch slower ones. When the smooth inviscid description would become multivalued, a physical compression shock replaces it. A rarefaction normally spreads rather than forming an entropy-violating expansion shock.^[ref]

Across a locally planar shock, use velocities relative to its moving front. The Rankine-Hugoniot conditions for a perfect gas conserve mass, momentum and total energy flux:where $h=\gamma p/[(\gamma -1)\rho ]$ is the specific enthalpy. For upstream normal Mach number $M_1>1$,Thus a normal compression shock increases pressure, density and temperature, reduces the normal velocity, and changes a supersonic upstream flow to subsonic downstream flow in the shock frame. In the strong-shock limit, compression approaches $(\gamma +1)/(\gamma -1)$, which is four for a monatomic gas. A shock conserves total energy but generates entropy; an adiabatic shock is not an isentropic compression. For a calorically perfect gas, stagnation temperature is conserved while stagnation pressure falls.^[ref]

Real dissipative processes resolve the narrow front; the jump conditions do not require its detailed internal profile. In a continuum-gas experiment a shock tube can produce a moving shock by releasing a pressure discontinuity. pressure probes separated along the tube measure arrival speed and the pressure jump; distinguishing the shock from the following contact surface is essential. A contact surface can have a density/entropy jump with continuous pressure and normal velocity, whereas a shock carries nonzero mass flux through its front.

Supersonic communication. A point source moving at speed $V>c_s$ outruns its small acoustic waves. After a time $t$, the source has moved $Vt$, while a disturbance emitted at its old position has radius $c_{s}t$. The envelope of these spherical wavefronts is the Mach cone, with half-angle $\mu$ satisfyingInformation is restricted to the downstream cone, so upstream gas cannot smoothly anticipate a supersonic obstacle. Finite compressions can also steepen because their characteristics converge. Such disturbances form a shock wave, across which the fluid variables change rapidly. A perfectly uniform supersonic flow need not contain a shock: the cone expresses the communication constraint, while a disturbance or compression supplies the reason for a discontinuity.

Conserved fluxes. Work in the shock frame, take velocities normal to the front, and integrate the conservation laws through a thin steady control volume of unit area. Mass conservation gives constant flux $j=\rho u$, soThe normal momentum flux is the advected part $\rho u^2$ plus the pressure stress $p$. With negligible body-force impulse across the thin front,For an ideal gas with constant adiabatic index $\gamma$, the specific internal energy is $e=p/[(\gamma -1)\rho ]$. The energy flux iswhere $h_{\mathrm{sp}}$ is the specific enthalpy. No radiation or other heat flux escapes in this adiabatic shock model. Conservation of the energy flux, divided by the nonzero mass flux, givesThese are the Rankine-Hugoniot conditions for a perfect gas. The shock is adiabatic but irreversible: its entropy increases, so an isentropic relation cannot be imposed across the discontinuity.

Strong-shock compression. Put $\chi ={\rho }_2/{\rho }_1=u_1/u_2$ and neglect $p_1$ and the upstream enthalpy compared with ${\rho }_1{u}_1^2$ and $u_1^2$. Momentum conservation gives $p_2={\rho }_1{u}_1^2(1-1/\chi )$. Substitute this into energy conservation:ThusThe root $\chi =1$ describes no shock. The compressive branch of the Strong-shock Rankine-Hugoniot conditions isIt also gives $p_2=2{\rho }_1{u}_1^2/(\gamma +1)$. A monatomic gas has $\gamma =5/3$, hence $\chi =4$.

The finite slab in its initial rest frame. Denote the shock's lab speed by $V_s$. The upstream gas is initially stationary, so $u_1=V_s$ and the downstream shock-frame speed is $u_2=V_s/4$. Transforming to the initial slab frame givesThis is planar shock compression of a finite slab. In the uniform post-shock approximation, breakout occurs at $t_b=L/V_s$. The original front face has then moved $(3V_s/4)t_b=3L/4$, while the shock reaches the originally stationary far face at $L$, leaving thickness $L/4$. Equivalently, the fixed mass per area $\rho L$ becomes $4\rho L_2$. The result assumes a maintained approximately constant shock speed and no important boundary rarefaction during the crossing; an unconfined finite cloud need not remain uniformly compressed.

Escape estimate and the sound-speed convention. Let $M_c$ be the slab mass, $W>0$ its binding-energy magnitude in the star's field, and $E_{\mathrm{th}}$ its initial internal energy. Treat the potential as nearly uniform across the slab. For the bulk-ballistic escape estimate, requireThe Mach-number convention in shock unbinding estimates matters here. With the usual upstream adiabatic sound speed, $c_{\mathrm{ad}}^2=\gamma p/\rho =5p/(3\rho )$, define ${\mathcal{M}}_{\mathrm{ad}}=V_s/c_{\mathrm{ad}}$. ThenIf instead the shock speed is normalized by the isothermal sound speed, $c_T^2=p/\rho$, then ${\mathcal{M}}_T^2=(5/3){\mathcal{M}}_{\mathrm{ad}}^2$ andThe original PDF's “around five” is consistent with the latter convention or an order-of-magnitude sound-speed estimate. It is not the exact coefficient obtained with the conventional adiabatic Mach number and the stated bulk kinetic-energy criterion. Here $W=-E_{\mathrm{grav}}$ if the potential is zero at infinity; using negative gravitational energy directly on the right would be incorrect. Subsequent conversion of shock-generated internal energy to expansion, external pressure, or radiative losses would change the physical escape threshold, so this numerical factor is not universal.

A shock wave forms when a compressive flow steepens or when gas must be decelerated or turned faster than upstream acoustic waves can communicate the change. Astrophysical examples include a supernova remnant sweeping up surrounding gas, a fast stellar wind colliding with ambient material, and supersonic accretion onto an obstacle or onto a dense gas layer. A large Mach number by itself is not sufficient: a uniform supersonic stream need not contain a shock wave.

Choose the shock frame, and denote upstream and downstream normal velocities by $u_1,u_2$. Integrating the local conservation laws over a thin control volume spanning a steady normal shock wave gives equal fluxes on both sides. Mass conservation gives constant mass flux per unit area,Pressure and advected momentum together give the normal momentum flux,For a thin non-radiating layer, with negligible gravitational-potential change across it, conservation of energy gives equal fluxes $j(h+u^2/2)$. A perfect gas has specific enthalpy $h=\gamma p/[(\gamma -1)\rho ]$, soThese are the Rankine-Hugoniot conditions for a perfect gas. Dissipation inside the thin layer produces an entropy increase and selects the physical compressive branch. Efficient radiative cooling would change the energy relation and may instead lead to an isothermal shock.

The Rankine-Hugoniot conditions follow by integrating conservation laws across a thin moving layer, then letting its macroscopic thickness tend to zero. For example, normal momentum conservation in one dimension has the conservative formFor a shock wave moving at speed $V$, integrate over a control interval whose endpoints move with the discontinuity. The Reynolds transport theorem giveswhere $[f]=f_R-f_L$ denotes the difference of limiting states. Mass conservation similarly gives $V[\rho ]=[\rho u]$. If $v=u-V$ is the normal velocity in the shock frame, combining the two relations yieldsThe same construction applied to conservation of energy, for a non-radiating normal shock wave with nonzero mass flux, gives $[h+v^2/2]=0$, where $h$ is the specific enthalpy. An equation of state closes these relations, and the physical compression branch has increasing entropy.

Finite dynamic viscosity replaces the ideal discontinuity by a thin smooth transition. A viscous stress contributes a derivative to the momentum flux inside the layer and converts directed kinetic energy into internal energy. If its gradients vanish in the uniform far states, the same limiting Rankine-Hugoniot conditions still hold there. Balancing nonlinear transport with viscous diffusion gives a thickness of order $\nu /|u_L-u_R|$ in the viscous Burgers equation, where $\nu$ is the kinematic viscosity.

A shock wave is a thin, irreversible compression layer across which fluid density, pressure, temperature and velocity change rapidly, while conserved fluxes satisfy the Rankine-Hugoniot conditions. In the shock frame, a compressive hydrodynamic shock converts ordered upstream kinetic energy into downstream thermal energy and increases entropy. The upstream normal flow is supersonic; the downstream normal flow is subsonic for an ordinary gas-dynamic compression shock.

In an accreting binary star, a supersonic stellar wind can be captured by a compact companion and shocked. An accretion stream striking an accretion disk, an accretion-disk boundary layer, and a magnetically guided accretion column provide other sites where infall energy is dissipated and high-energy radiation is generated. Shocks heat and decelerate material; a disk also needs transport of angular momentum before its material can reach the compact star.

In diffusive shock acceleration, charged particles scatter from magnetic irregularities on both sides of a shock and repeatedly cross the converging flows. The average energy gain is first order in the velocity jump divided by the particle speed. Escape competes with acceleration, producing a nonthermal distribution rather than a single downstream temperature. Supernova remnants are important candidate sources of Galactic cosmic rays; shocks associated with active galactic nuclei, their jets and lobes, and gamma-ray bursts are candidate extragalactic accelerators. Ultrarelativistic energies require sufficient size, magnetic field and lifetime, and losses must not outpace acceleration. In particular a relativistic particle's Larmor radius must fit within the accelerator; the existence of a strong shock alone does not guarantee unlimited particle energy.

Use a steady, plane, nonradiative normal shock wave in a perfect gas, and measure the normal speeds $v_{\mathrm{u}},v_{\mathrm{d}}>0$ in the shock frame. The labels denote upstream and downstream. Conservation of mass, normal momentum and total energy givesandThe enthalpy term includes the pressure work. For very large upstream Mach number, neglect $p_{\mathrm{u}}$ compared with ${\rho }_{\mathrm{u}}{v}_{\mathrm{u}}^2$. If $r={\rho }_{\mathrm{d}}/{\rho }_{\mathrm{u}}$ is the shock compression ratio, the first two equations implySubstitution into the energy equation givesDiscard the no-shock solution $r=1$. For a nonrelativistic monatomic gas, $\gamma =5/3$, so the Strong-shock Rankine-Hugoniot conditions becomeThe factor four requires this adiabatic index and negligible radiative, magnetic and escaping-particle contributions to the flux budget. For example, the same strong, nonradiative calculation with $\gamma =4/3$ gives compression seven, not four.

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.

Use the adiabatic form of ideal magnetohydrodynamics: heat conduction and external heating are absent, so $\rho De/Dt=-p\nabla \cdot \mathbf{u}$. Infinite electrical conductivity and zero viscosity alone would not remove the thermal source terms on the cover page. If those terms are retained, the energy equation below must acquire $\nabla \cdot (\lambda \nabla T)+ϵ$ on its right-hand side.

Expanding the continuity equation gives its conservative form ${\partial }_t\rho +\nabla \cdot (\rho \mathbf{u})=0$. Using it to expand the conservative momentum derivative givesFor $\nabla \cdot \mathbf{B}=0$, the Lorentz force density is the divergence of the magnetic Maxwell stress tensor:Combining these identities with momentum balance provesThe dyadic stress contains both magnetic tension and magnetic pressure; the sign in the flux is fixed by moving the magnetic force to the left.

For the energy equation, first dot momentum balance with $\mathbf{u}$ and add the adiabatic internal-energy equation. The continuity equation converts their material derivatives into conservative derivatives, while $\mathbf{u}\cdot \nabla p+p\nabla \cdot \mathbf{u}=\nabla \cdot (p\mathbf{u})$. The result isThe ideal magnetohydrodynamic induction equation gives the magnetic counterpart. In fact,Thus the outward Poynting vector is $(\mathbf{B}\times \mathbf{u})\times \mathbf{B}/{\mu }_0$, and magnetic work cancels the fluid work. Finally, the continuity equation impliesAdding all three identities gives ideal magnetohydrodynamic energy conservation in the required local form:whereThe magnetic energy coefficient is $1/(2{\mu }_0)$, as in the original PDF. For a self-consistent, time-dependent gravitational field this local balance should not be confused with a global self-gravitating energy integral: the latter uses gravitational field energy, or the familiar one-half factor in $\int \rho \Phi$ after the field is eliminated.

Work in the shock frame of a thin planar discontinuity, with normal speed $u>0$. A smooth gravitational potential has no finite jump through a vanishingly thin front; its advected contribution cancels using the common mass flux. Integrating the conservation laws gives constant normal mass, momentum and energy fluxes, together with the normal-field and tangential-electric-field conditions of the ideal magnetohydrodynamic shock conditions.

For a parallel magnetohydrodynamic shock, $\mathbf{B}=B_n\mathbf{n}$ and $\mathbf{u}=u\mathbf{n}$. The magnetic solenoidal constraint makes $B_n$ the same on both sides. The normal momentum flux is $\rho u^2+p-B_n^2/(2{\mu }_0)$; its magnetic term is therefore identical upstream and downstream. The Poynting vector is zero because $\mathbf{u}\times \mathbf{B}=0$. Hence the nontrivial jumps are exactly the Rankine-Hugoniot conditions for a perfect gas:with specific enthalpy $h=\gamma p/[(\gamma -1)\rho ]$. Physically, motion along the straight field neither bends it nor compresses its transverse flux. The field consequently provides no changed stress or energy flux in this purely parallel geometry. This conclusion does not apply to a shock that develops a transverse magnetic field.

For a perpendicular magnetohydrodynamic shock, the magnetic field is tangential and the velocity normal. Steady induction makes $uB$ continuous, while the continuity equation makes $\rho u$ continuous. Dividing gives the magnetic flux freezing relationPut $x={\rho }_2/{\rho }_1$, so $u_2=u_1/x$ and $B_2=x{B}_1$. Constant normal momentum flux now reads $\rho u^2+p+B^2/(2{\mu }_0)$. Substitution givesThe magnetic correction is negative for compression: part of the incoming ram pressure builds magnetic pressure instead of thermal pressure.

Here the Poynting vector has normal component $u{B}^2/{\mu }_0$. Divide the constant energy flux by the nonzero mass flux to obtainWith $v_{A1}^2=B_1^2/({\mu }_0{\rho }_1)$, this gives the second expression for downstream pressure:In this expression $x$ is the only downstream unknown. Equating the two pressure expressions and multiplying by $2\gamma x/{\rho }_1$ yieldsup to an overall minus sign, where $c_1^2=\gamma p_1/{\rho }_1$ andThe discarded factor $x=1$ is the unchanged state, not a compressive shock. This derives the polynomial for the compression ratio of a perpendicular magnetohydrodynamic shock without omitting its trivial branch.

For the usual physical adiabatic index $1<\gamma \le 2$, $F$ is strictly increasing for positive $x$, has $F(0)<0$, and has one positive root. Sincethat root is greater than one precisely when $u_1^2>c_1^2+v_{A1}^2$. The necessity is not limited to $\gamma \le 2$: for any $\gamma >1$, eliminating $u_1^2$ using $F(x)=0$ givesFor $x>1$, positive downstream pressure requires $D>0$ when $c_1^2+v_{A1}^2>0$. Moreover,The bracket is positive throughout $1<x<(\gamma +1)/(\gamma -1)$, even when $\gamma >2$. The exactly cold unmagnetized limit gives the familiar strong-shock upper compression ratio directly. Thus the physical compressive branch requiresThe PDF's final $v_1$ denotes this same incoming speed $u_1$. The right side is the squared speed of a perpendicular fast magnetosonic wave: the upstream flow must be super-fast-magnetosonic so that this compressive signal cannot run upstream against the incoming flow. At equality the shock amplitude vanishes into a characteristic disturbance.

Use signed velocities in the shock frame; the Mach number depends on their squares. For a genuine nonradiative perfect gas shock wave, the Rankine-Hugoniot conditions for a perfect gas give a common mass flux $j={\rho }_L{u}_L={\rho }_R{u}_R$, a common momentum flux $p+\rho u^2$, and a common specific stagnation enthalpy $\gamma p/[(\gamma -1)\rho ]+u^2/2$. Here $e$ in the given equation of state is specific internal energy, despite the source's informal “energy density” wording.

Put $r={\rho }_R/{\rho }_L$, $P=p_R/p_L$ and $M^2={\rho }_L{u}_L^2/(\gamma p_L)$. Mass and momentum give $u_R=u_L/r$ and $P=1+\gamma M^2(1-1/r)$. Substituting these into the energy relation and clearing denominators givesThe $r=1$ factor is the no-jump solution. On the nontrivial shock wave branch,In particular $r-1=2(M^2-1)/[(\gamma -1)M^2+2]$. With $\gamma >1$, compression is therefore equivalent to $M_L^2>1$. The other Mach number isPositive pressure makes its denominator positive, and $M_R^2-1$ has the opposite sign to $M_L^2-1$. Thus a nontrivial physical shock wave has a supersonic side and a subsonic side; the entropy-admissible flow goes from the former to the latter.

The shock compression ratio also yields $u_L{u}_R=(\gamma -1)u_L^2/(\gamma +1)+2\gamma p_L/[(\gamma +1){\rho }_L]$. Rearranging produces the first velocity identity in the question. The pressure relation gives the second identity in the equivalent formThe numerator here contains mass density, as the original PDF shows; the converted TeX incorrectly replaces that mass density by pressure.

For reflection, the once-shocked gas is the common left state in both applications. Before rebound, its velocity in the shock frame is $a=u_s-U_{+}<0$ and the gas ahead has velocity $-U_{+}$; after rebound, it is $b=u_s+U_{-}>0$ and the gas next to the wall has velocity $U_{-}$. In both cases the right velocity minus the left velocity equals $-u_s$. Applying the first identity therefore gives the same quadratic equation,with roots $a$ and $b$. Their opposite signs distinguish the two roots, so $\boxed{(u_s-U_{+})(u_s+U_{-})=-\gamma p_s/{\rho }_s}$.

Write $q_0=p_0/p_s$, $q_1=p_1/p_s$. The pressure identity applied to each shock wave givesSquaring the root-product relation and eliminating the velocities yieldsand hence the normal shock reflection at a rigid wall relationTaking $p_0/p_s\to 0$ gives $\boxed{p_1/p_s\to (3\gamma -1)/(\gamma -1)}$; for a monatomic gas this limiting amplification is six.

Assume $\gamma >1$. The upstream ideal gas is at rest in the laboratory frame and has mass density ${\rho }_1=C{R}^{-\beta }$. In the shock frame its incident speed is $V$. The strong-shock limit of the Rankine-Hugoniot conditions for a perfect gas gives compression $\chi =(\gamma +1)/(\gamma -1)$ and downstream shock-frame speed $V/\chi$. Transforming back to the laboratory frame gives $u_2=V(1-1/\chi )=2V/(\gamma +1)$.

The normal momentum jump yieldswhere the cold upstream pressure has been neglected. Thus the boundary values in the normalization of part (b) areThe velocity jump uses the relative speed $V-u_2$, not the laboratory speed $u_2$. This distinction is essential to all three shock wave values.