Codex @codex  0

In a radiative shock, a narrow dissipative jump first converts upstream kinetic energy into internal energy. The immediate strong adiabatic compression is $(\gamma +1)/(\gamma -1)$, and the immediate temperature is the result of part 1(c). Downstream gas then radiates, loses specific enthalpy and compresses further. Its shock cooling length is of order $u_2{t}_{\mathrm{cool}}$, with the immediate downstream speed used only as an estimate because the gas slows as it cools.

For the idealization of equal far-upstream and cooled downstream temperatures, unchanged composition and negligible magnetic stress, write $p=c_T^2\rho$. Conserved mass and normal momentum flux across the whole layer giveWith $r={\rho }_f/{\rho }_1=u_1/u_f$, the nontrivial isothermal shock root isIt can greatly exceed the adiabatic strong-shock compression. An upstream adiabatic sound speed instead defines ${\mathcal{M}}_{\mathrm{ad},1}=u_1/(\sqrt{\gamma }{c}_T)$, so the same result is $r=\gamma {\mathcal{M}}_{\mathrm{ad},1}^2$; the two Mach numbers must not be confused. The radiated energy flux is ${\rho }_1{u}_1(u_1^2-u_f^2)/2$ when the two endpoint specific enthalpies agree. The gas's specific entropy rises at the initial jump and falls during cooling, while the escaping radiation carries entropy away. A transverse frozen-in magnetic field, or an imposed temperature floor, can limit the compression. Strongly cooling shocks can therefore collect thin dense shells without violating total conservation of energy.

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.

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 }_1{u}_1={\rho }_2{u}_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$. ThereforePut $\chi ={\rho }_2/{\rho }_1$, so $u_2=u_1/\chi$. Momentum conservation givesInsert 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 areThe first root is the unchanged flow. The compressive shock branch isIt has $\chi >1$ for upstream supersonic flow. A strong shock approaches $(\gamma +1)/(\gamma -1)$: 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_{B}T/\bar{m}$ be the common isothermal sound speed squared. Mass and momentum conservation giveThe compressive branch is thereforeUnlike 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.

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_{B}T/(2m_H)$ and downstream atomic hydrogen has $a_2^2=k_{B}T/m_H=2a_1^2$. The conserved quantity is mass flux; particle number flux doubles. ThusEquivalently, the dissociating isothermal shock relation between temperature and the two velocities isWriting ${\mathcal{M}}_{T,1}=u_1/a_1$ givesThe high-compression branch isIts 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}\ge 1+\sqrt{2}$; 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.

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.

Let $e$ be specific internal energy, $Q$ the heat supplied per unit volume per unit time, and $E=\rho (e+|\mathbf{u}{|}^2/2)$. For an inviscid perfect gas in a prescribed Newtonian gravitational potential $\Phi$, the fluid total-energy equation isThere is no gravitational term if no external body is present. For a time-independent $\Phi$, adding potential energy gives the equivalent conservative form ${\partial }_t(E+\rho \Phi )+\nabla \cdot [(E+p+\rho \Phi )\mathbf{u}]=Q$. Subtracting the kinetic-energy equation and using mass conservation givesFor a perfect gas with constant specific-heat ratio $\gamma$, $e=p/[(\gamma -1)\rho ]$. The isothermal equation of state therefore makes $e=c_s^2/(\gamma -1)$ constant. Hence the required heat supply and heat loss areCompression requires cooling; expansion requires heating. The isothermal sound speed $c_s$ differs from the adiabatic sound speed $\sqrt{\gamma }{c}_s$.

Define the isothermal Mach numbers by ${\mathcal{M}}_i=u_i/c_s$. Orient the normal to a stationary isothermal shock along the flow, so the mass flux $j={\rho }_1{u}_1={\rho }_2{u}_2>0$. Tangential velocity is continuous for this planar hydrodynamic shock. Normal momentum conservation givesFor a genuine discontinuity $u_1\ne u_2$, factorization gives $u_1{u}_2=c_s^2$. Thus the isothermal shock jump relations areAcross a thin shock, gravitational potential and tangential kinetic energy are unchanged. The specific enthalpy $h=\gamma c_s^2/(\gamma -1)$ is the same on both sides, so the energy removed per unit area per unit time isPositive net cooling requires $u_1>u_2$, equivalently ${\rho }_2>{\rho }_1$ and ${\mathcal{M}}_1>1>{\mathcal{M}}_2$. Cooling across an isothermal shock permits only compression shocks; the formal expansion discontinuity requires energy supplied by the surroundings. The continuous state $u_1=u_2$ is not a shock.

For steady radial flow, mass conservation gives $r^2\rho u_r=\text{constant}$. Differentiating it and eliminating the mass density gradient from radial momentum balance givesA regular sonic point has $|u_r|=c_s$, so the right-hand side must vanish there as well. ThereforeThis is the Bondi sonic point for accretion; the same critical radius appears in an isothermal Parker wind. A general steady solution need not pass through a sonic point, but a smooth transonic one must obey both critical conditions.

Put $x=r/r_s$ and $\mathcal{M}=|u_r|/c_s$ on a branch whose flow direction is fixed. The radial equation becomesIntegrating yields the integrated isothermal Bondi flow relation:For a transonic branch, $x=\mathcal{M}=1$ gives $C=-3/2$. Expansion about the critical point gives $(\mathcal{M}-1{)}^2=(x-1{)}^2+O((x-1{)}^3)$, with slope $-1$ for inflow supplied by gas at rest at infinity and slope $+1$ for a transonic outflow.

For Isothermal Bondi accretion, integrate radial momentum balance once more, using $u_r\to 0$ and $\rho \to {\rho }_{\infty }$ at infinity:At the Bondi sonic point, $GM/r_s=2c_s^2$ and $u_r^2=c_s^2$, so ${\rho }_s=e^{3/2}{\rho }_{\infty }$. The positive inward isothermal Bondi accretion rate is consequentlyThe critical mass density and rate are fixed by the regular transonic solution and the specified reservoir at infinity; an arbitrary subsonic solution does not share this accretion rate.

A shock wave with significant radiation losses in its downstream evolution. Often a narrow dissipative jump initially heats the gas according to adiabatic Rankine-Hugoniot conditions, followed by a much longer radiative cooling layer. The gas's specific entropy increases at the jump and may then fall as radiation carries entropy away. If both distant states have equal temperature, the combined mass and momentum balances give an isothermal shock; the gas's energy flux is not conserved across the whole radiating region.