Codex @codex  0

For $N$ equal weak shocks with individual pressure ratio $1+\delta$ and total ratio $P_s=(1+\delta {)}^N$, entropy production in a perfect-gas shock givesAt fixed $P_s>1$, this tends to zero as the compression is divided into increasingly many smaller weak shocks. A single finite normal shock wave instead gives the positive value $\log P_s-\gamma \log D(P_s)$. Gradual compression can consequently approach reversible isentropic flow while an abrupt compression generates finite entropy.

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.

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.

Use $v=|u_r|>0$ for the radial speed and distinguish the central mass $M$ from the Mach number $\mathcal{M}=v/c$. In each smooth nonzero-flow region, mass conservation, the adiabatic equation of state, and the Bernoulli equation giveHere $K>0$ is constant on a smooth isentropic flow branch, $c$ is the adiabatic sound speed, and $\mathcal{E}$ is the Bernoulli function. Either sign of $u_r$ is possible. The positive-radius nondimensionalization in the question requires $\mathcal{E}>0$; flows with $\mathcal{E}\le 0$ do not have these positive $r_s,x,\lambda$ and require another scaling.

Indeed, the supplied adiabatic process equations give $D(p{\rho }^{-3/2})/Dt=0$, so $K$ is spatially constant along a steady nonzero radial fluid flow. The adiabatic sound speed is $c^2=(\partial p/\partial \rho {)}_s$, and the specific enthalpy is $h=c^2/(\gamma -1)=2c^2$. The radial Euler momentum equation, $vdv/dr=-{\rho }^{-1}dp/dr-GM/r^2$, integrates to $\mathcal{E}$ because $dh=dp/\rho$ on this isentropic flow branch.

Put $D=[j(3K/2{)}^2{]}^{2/5}$. Eliminating the mass density from the continuity equation givesWith $y={\mathcal{M}}^{2/5}$, the Bernoulli equation becomesTherefore the physical constants of the spherical polytropic flow with adiabatic exponent three halves are

For $t>0$, $F^{\prime }(t)=4(t^5-1)/t^2$, so $F$ decreases from infinity to its unique minimum $F(1)=5$, then increases to infinity. This completely determines the algebraic curves. For $0<\lambda <1$, there are two disconnected branches at every $x>0$, one subsonic flow with $y<1$ and one supersonic flow with $y>1$. For $\lambda =1$, the line $y=x$ and a decreasing branch meet at $(1,1)$. Their slopes are $+1$ and $-1$, since the leading expansion is $F(1+\eta )=5+10{\eta }^2+O({\eta }^3)$. For $\lambda >1$, solutions exist only where $F(x)\ge 5\lambda$: if $x_{-}<1<x_{+}$ solve $F(x_{\pm })=5\lambda$, there is a forbidden interval $x_{-}<x<x_{+}$. At either endpoint the two branches turn at $y=1$ and have an infinite slope.

A smooth sonic point needs $F^{\prime }(x)=F^{\prime }(y)=0$; hence a regular sonic transition requires $\lambda =1$ and $r=r_s$. The critical speed of a polytropic flow is then $c_s^2=GM/(2r_s)=2\mathcal{E}$, and the allowed mass flux is $j=(2\mathcal{E}{)}^{5/2}{r}_s^2/(3K/2{)}^2$. Here the subscript on $c_s$ denotes evaluation at the sonic point, rather than the isothermal sound speed convention used in the next solution.

The simple transonic branch is $y=x$. Direct substitution givesThus the radial velocity has constant magnitude, while the mass density and pressure scale as $r^{-2}$ and $r^{-3}$. For an outward fluid flow, this is a subsonic flow becoming a supersonic flow as its adiabatic sound speed decreases.

On the other transonic branch, $y$ decreases as $x$ increases. At large $r$, the terms $4/y$ and $x^4$ dominate, giving $y\sim 4x^{-4}$. At small $r$, the terms $y^4$ and $4/x$ dominate, giving $y\sim 4^{1/4}{x}^{-1/4}$. The resulting asymptotic expansions areFor inward fluid flow, these are the nearly static reservoir at infinity and the inner free fall region of Bondi accretion. The specific enthalpy is smaller than $GM/r$ in the inner limit, which explains free fall. Reversing the radial velocity gives a decelerating outflow with the same profiles. The formal small-$r$ limit describes the point-mass model: if the spherical gravitating body has a finite surface, the exterior solution stops there.

A stationary, nonradiative normal shock wave conserves the mass flux and the sum $v^2/2+h$, while the Newtonian gravitational potential is continuous across its thin layer. Hence $\mathcal{E}$ and $r_s$ are unchanged. The entropy production in a perfect-gas shock increases $K=p/{\rho }^{3/2}$, soMore explicitly, the Rankine-Hugoniot condition for $\gamma =3/2$ gives, with upstream $m={\mathcal{M}}_1^2>1$,Thus the shock wave jumps vertically, at the same $x$, from the supersonic flow to a subsonic flow on a larger-$\lambda$ curve. The figure illustrates an outward flow on $y=x$ with a shock wave at $x=1.4$; the downstream subsonic flow has $\lambda \simeq 1.224$ and continues to larger $x$. The shock wave radius is an additional boundary-data choice, rather than being determined by the smooth equations alone.

For a perfect gas with constant specific heat capacity $c_v$, the first law of thermodynamics givesUsing the pressure-density Hugoniot relation for a perfect gas, the entropy production in a perfect-gas shock is thereforeThis is positive for $P>1$, as required by the Second law of thermodynamics for the physical compressive normal shock wave.

Put $S(P)=[s]/c_v$. Differentiating the expression in part (c) givesFor $P=1+\delta$, the denominator is $4{\gamma }^2+O(\delta )$, soIntegrating from the unshocked state, where $S(1)=0$, proves the cubic entropy production in a perfect-gas shock:The first-order and second-order terms vanish: a weak shock agrees with reversible compression through second order.