Codex @codex  0

Take the outward wavenumber $k>0$; the absolute value in the source makes the two propagation directions equivalent. The phase velocity of the gravity-capillary wave obeysDifferentiate with respect to $k$:The second derivative is positive, and $c_p\to \infty$ at both ends of the positive-$k$ range, so this is the global minimum. The minimum phase speed of a gravity-capillary wave isFor $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$ and the supplied surface tension and mass density,The corresponding wavenumber is $k_{\ast }\simeq 369{\mathrm{m}}^{-1}$, giving wavelength ${\lambda }_{\ast }=2\pi /k_{\ast }\simeq 1.70\mathrm{cm}$. At this minimum, $c_g=d\omega /dk=c_p+kd{c}_p/dk=c_p$. Thus the group velocity and phase velocity agree for this particular ripple, although they differ for general gravity-capillary waves.

A highly conducting solar wind carries the Sun's magnetic field outward by magnetic flux freezing. This can open magnetic field lines that would form closed loops in a static dipole model. Beyond an opening radius $r_o$, approximate the field as radial. Constancy of the magnetic flux in a fixed solid-angle flux tube gives the radial solar-wind magnetic flux lawwhich decays more slowly than the vacuum magnetic dipole field, proportional to $r^{-3}$.

For example, if the field becomes radial already near the solar radius $R_{\odot }$, normalizing both estimates to the same field there givesThus a dipole estimate $10^{-11}\mathrm{T}$ becomes of order $2\times 10^{-9}\mathrm{T}$, comparable to the stated hundredfold enhancement. More generally, matching dipolar decay inside $r_o$ to radial decay outside givesThis is a plausible scaling explanation, not a determination of the actual opening radius. The magnetic field need not be outward everywhere: opposite-polarity sectors can preserve zero net magnetic flux through a closed surface, consistent with Gauss's law for magnetism. The argument concerns the changed geometry of the field in a flowing plasma, rather than changing the dipole exponent by assumption.

The first argument mixes a material fluid element with a fixed spatial volume. For a material volume element following the fluid, its volume is generally time dependent:Therefore the material derivative cannot be taken outside the product with $\delta V$ as in the proposed cancellation. Instead, mass conservation statesThe momentum balance for a material fluid element, using Newton's second law, is consequentlyFor a fixed control volume, fluid crosses its boundary and carries momentum; a momentum flux term must be included. Its local conservative form is ${\partial }_t(\rho \mathbf{u})+\nabla \cdot (\rho \mathbf{u}\otimes \mathbf{u})=-\nabla p$, equivalent to the displayed Euler equations after using the continuity equation. The erroneous equation omits the term $\rho \mathbf{u}\nabla \cdot \mathbf{u}$ needed to convert a material derivative of momentum density into this conservative form.

For the second task, take a homogeneous fluid with constant mass density $\rho$, prescribed gravitational acceleration $\mathbf{g}$, incompressible flow and inviscid flow. Its equations areUnder the replacement $\rho \mapsto \tilde{\rho }=C\rho$, the same velocity solves these equations if $p\mapsto \tilde{p}=Cp$, up to a spatially constant pressure offset. This is density scaling invariance of incompressible Euler flow:The comparison assumes the same geometry and compatible velocity or gravitational free surface conditions; an independently specified nonuniform boundary pressure would supply additional information. The gravitational field is taken as prescribed. If the fluid's own self-gravity must be included, its field changes with mass density, so the invariance need not hold. A single freely moving particle's density-independent gravitational acceleration does not justify neglecting that distinction.

Hydrostatic equilibrium with self-gravity. A resting fluid obeys $\nabla p=-\rho \nabla \Psi$. Spherical symmetry and the Poisson equation for Newtonian gravity giveEliminate the Newtonian gravitational potential derivative to obtainEquivalently, the shell theorem gives ${\Psi }^{\prime }=GM(r)/r^2$ and $M^{\prime }(r)=4\pi r^2\rho$, which leads to the same hydrostatic equilibrium equation.

A uniform-density liquid. With constant $\rho$, $M(r)=4\pi \rho r^3/3$. There is no central point mass, so the interior gravitational acceleration magnitude is $4\pi G\rho r/3$. ThusUse $p(R)=0$ at a vacuum boundary, neglecting surface tension. The uniform-density planetary pressure profile isThe central value can also be writtenIt has units pascals and decreases smoothly from its maximum at the centre to zero at the surface.

An isothermal gas. Here $c$ is the constant isothermal sound speed, with $p=c^2\rho$. Try $p=A/r^2$, so $\rho =A/(c^2{r}^2)$. The two sides of the equilibrium equation becomeThus $A=c^4/(2\pi G)$, giving the singular isothermal sphere,This does not describe a finite gaseous sphere in vacuum. Both pressure and mass density diverge at the centre, and $p(R)>0$ for every finite outer radius $R$, so it cannot meet $p=0$ there. Extending it to infinity removes the finite-radius pressure jump but makes the total mass infinite. It is useful as an idealized profile over an intermediate radial range, or as a truncated model with nonzero external pressure. The vacuum boundary obstruction for an isothermal sphere also applies to the regular, nonsingular isothermal solution; removing the central singularity does not produce a finite vacuum surface.

Steady outflow. In the wind calculation, take the central mass $M$ to be fixed and neglect the wind's own self-gravity. Constant mass flow rate and the isothermal equation of state giveThe steady radial Euler equation isSubstituting the density derivative givesThis is the isothermal Parker wind equation. A smooth sonic point requires both sides to vanish when $u=c$, soMaintaining constant temperature during expansion generally requires heat input; it is an isothermal model, not the adiabatic blast-wave model of Question 5.

The escaping branch at large radius. Integrate the wind equation, using dimensionless logarithm arguments:Set $x=r/r_s$ and $w=u^2/c^2$. ThenFor the accelerating transonic branch, $w(1)=1$, which fixes $C_0=-3$. At $x\gg 1$ on its supersonic branch, $w\to \infty$ and $\ln w$ is smaller than $w$. ConsequentlyA first correction makes the ordering explicit:This is the large-radius asymptotics of an isothermal Parker wind. The leading expression also holds for other supersonic branches with finite $C_0$.

The source implicitly selects an escaping wind. Positive steady outward flow alone does not determine that branch: a subsonic isothermal breeze also satisfies the differential equation. Its $w\ll 1$ limit gives $w\sim e^{-C_0}{x}^{-4}$, hence $u\propto r^{-2}$ and $\rho \to {\rho }_{\infty }>0$. Such a breeze needs nonzero pressure at infinity. Selecting a regular outflow from a subsonic base into vanishing ambient pressure instead gives the accelerating transonic wind used above.

The asymptotic speed has the same units as $c$; the logarithm is dimensionless. Its unlimited growth is a formal feature of indefinite isothermal heating and a prescribed point-mass potential. Real heating, surrounding pressure, and eventually the neglected wind mass limit that extrapolation.

The local conservation laws. The first equation is the radial Euler equation: the material acceleration $Du/Dt={\partial }_{t}u+u{\partial }_{r}u$ is caused by the pressure gradient force per unit mass, $-{\partial }_{r}p/\rho$. It neglects gravity and viscosity. The second equation is mass conservation in spherical coordinates,The term $2\rho u/r$ appears because the area of a spherical surface grows as $r^2$.

To obtain the third equation, let $e=p/[(\gamma -1)\rho ]$ be the specific internal energy of an ideal gas with constant heat capacity ratio $\gamma$. The first law of thermodynamics for a smooth adiabatically moving element isDifferentiating $e$ and collecting terms givesand henceThus the entropy parameter $p/{\rho }^{\gamma }$ is conserved along each smooth pathline; for fixed composition, specific entropy differs from $c_v\ln (p/{\rho }^{\gamma })$ only by a constant. This is local isentropic flow behind the shock. A physical shock wave itself generates entropy and must be treated with jump conditions, rather than by applying this smooth equation across the discontinuity.

The front-radius scaling. In the self-similar blast wave approximation there is no scale set by ambient pressure, cooling, gravity, or the initial source radius. Write $R\propto E^a{\rho }_0^b{t}^d$. Using $[E]=\mathsf{M}{\mathsf{L}}^2{\mathsf{T}}^{-2}$ and $[{\rho }_0]=\mathsf{M}{\mathsf{L}}^{-3}$, dimensional analysis givesThus $a=1/5$, $b=-1/5$, $d=2/5$ andThe dimensionless coefficient ${\xi }_0$ depends on $\gamma$ and the resolved similarity profiles; dimensions alone do not fix it. Define the length scale $\ell (t)=(E{t}^2/{\rho }_0{)}^{1/5}$, so $\xi =r/\ell$ and the front is at the fixed coordinate $\xi ={\xi }_0$.

Kinetic energy. Substitute $\rho ={\rho }_0{f}_{\rho }(\xi )$ and $u=(r/t)f_u(\xi )$ into the volume integral. Since $r=\ell \xi$ and $dr=\ell d\xi$,But ${\rho }_0{\ell }^5/t^2=E$, so the first of the spherical blast-wave similarity energy integrals isThe conclusion requires a fixed, integrable similarity profile, as in the finite-energy blast solution. Expansion increases the swept-up mass while reducing the characteristic squared speed, leaving this product constant.

Pressure and internal energy. The derivatives of the similarity variable areWriting a prime for differentiation with respect to $\xi$, the material acceleration becomesThe radial momentum equation therefore requires a pressure gradient on the scale ${\rho }_0\ell /t^2$. Integrating over radius gives pressure similarity in a spherical blast wave:whereThe prefactor has units $\mathrm{kg}{\mathrm{m}}^{-1}{\mathrm{s}}^{-2}$, as required for pressure. Strictly, the momentum equation alone permits an additional pressure uniform in space and arbitrary in time. In the scale-free strong-blast solution, shock matching and the remaining thermodynamic equation fix that freedom with the same similarity scaling; any allowed constant of integration is included in $f_p$. An independently imposed external pressure would introduce another scale and invalidate the simple dimensional ansatz.

For constant $\gamma >1$, the internal energy is consequentlyIf $E$ is the total excess explosion energy, its normalization isThe two constant fractions need not be equal. No explicit solution for $f_p$ is needed to establish their time independence.

Following a fluid element. Let $\mathbf{X}(t)$ be a pathline, satisfying $\dot{\mathbf{X}}=\mathbf{u}(\mathbf{X},t)$. The material derivative, also called the Lagrangian time derivative, is the rate observed while following that element. The multivariable chain rule givesThe partial derivative holds position fixed; the advection term accounts for motion through a spatially varying field. For a vector field, the same formula applies component by component in fixed Cartesian axes.

Momentum balance. For an inviscid flow, pressure is the only surface stress. A small material volume $dV$ has net pressure force $-\nabla pdV$, obtained by the divergence theorem. Its mass is $\rho dV$ and its material acceleration is $D\mathbf{u}/Dt$. Applying Newton's second law gives the Euler equations for an inviscid fluid,This is the local momentum conservation equation when viscosity and physical body forces are absent. A prescribed body force per unit mass would add $\rho \mathbf{f}$ on the right.

Changing to rotating coordinates. In a rotating frame with constant angular velocity $\mathbf{\Omega }$, the lab acceleration is the rotating-frame material acceleration plus $2\mathbf{\Omega }\times \mathbf{u}+\mathbf{\Omega }\times (\mathbf{\Omega }\times \mathbf{r})$. ThereforeFor uniform rotation,With constant ${\rho }_0$, this gives the Euler equation in a uniformly rotating frame,The term subtracted from $p$ is the equilibrium centrifugal contribution. Dividing the bracket by ${\rho }_0$ defines a reduced pressure $\Pi =p/{\rho }_0-|\mathbf{\Omega }\times \mathbf{r}{|}^2/2$, with units ${\mathrm{m}}^2{\mathrm{s}}^{-2}$. The Coriolis acceleration remains a velocity-dependent term.

Linear disturbances. Linearize about rest in the rotating coordinates, with the equilibrium part of $\Pi$ removed. The advective acceleration is quadratic in disturbance amplitude, so the linearized Euler equations and incompressible flow condition areTaking the curl removes the pressure gradient. The identitygives the inertial-wave vorticity equation,

Plane-wave frequencies. Let $\mathbf{k}\ne 0$ and use the given plane wave convention. Incompressible flow requires $\mathbf{k}\cdot {\mathbf{u}}_0=0$, and the vorticity equation givesFor a nonzero frequency, apply $\mathbf{k}\times$ again, using $\mathbf{k}\times (\mathbf{k}\times {\mathbf{u}}_0)=-k^2{\mathbf{u}}_0$. It follows thatThis is the inertial-wave dispersion relation. If $\theta$ is the angle of the wavevector to the rotation axis, the positive angular frequency is $2|\Omega \cos \theta |$. It is largest, $2|\Omega |$, for a parallel wavevector, decreases to zero for a perpendicular one, and is independent of the magnitude of the wavenumber at fixed direction. The two signs describe the two oscillatory polarization branches. When $k_z=0$, the zero-frequency result is obtained directly from the linear equations rather than by dividing by $\omega$; it describes the nonoscillating balanced limit.

For an inverted free surface, gravity promotes a displacement instead of restoring it. The surface tension term still resists curvature. Hence the gravity-capillary wave relation becomes, with $q=|k|$,For $0<q<k_c$, where $k_c=\sqrt{\rho g/\sigma }$, ${\omega }^2<0$. Writing the displacement as $\eta \propto e^{st+ikx}$ gives the Rayleigh-Taylor instability growth rateShorter wavelengths, $q>k_c$, remain oscillatory because surface tension dominates. The marginal scale is determined by the capillary length,A sufficiently small aperture excludes the long unstable modes: this is finite-aperture stabilization of an inverted liquid surface. The estimate requested is therefore a container opening of order one centimetre, with $1.7\mathrm{cm}$ the marginal full wavelength rather than a geometry-independent maximum diameter.

The exact largest opening cannot be fixed from the stated infinite-surface dispersion relation alone. It requires the aperture shape and contact-line boundary conditions. For example, a one-dimensional span of width $L$ whose contact line is pinned at both ends has a half-sine lowest mode with $q_{\min }=\pi /L$. In that simple model,Different geometries change the numerical factor. This is a small-disturbance stability estimate with negligible air inertia and no leakage or contact-line failure; a guarantee that a real inverted container retains its water needs those additional boundary details.

The exponential factor removes the first spatial derivative. Apply the chain rule explicitly:while ${\partial }_t^2\Delta \rho =e^{-z/(2H)}{\Psi }_{tt}$. Substituting into the supplied density equation, the terms proportional to ${\Psi }_z$ cancel, leavingThis is the atmospheric acoustic cutoff transformation. There is a sign error in the target equation in the original PDF, as well as in the transcription: the last term must be $+{\omega }_0^2\Psi$ for real positive ${\omega }_0$. If the displayed minus sign is retained literally, coefficient matching would require ${\omega }_0^2=-c^2/(4H^2)$, so it would not define a real characteristic oscillation period.

For the corrected equation, a plane wave $\Psi \propto e^{i(kz-\omega t)}$ has dispersion relationThus ${\omega }_0$ is the acoustic cutoff frequency: vertically propagating acoustic waves have $|\omega |\ge {\omega }_0$. The scale height and sound speed have units metres and metres per second, respectively, so ${\omega }_0$ has units ${\mathrm{s}}^{-1}$.

For an order-of-magnitude estimate in the lower terrestrial atmosphere, take $H\simeq 8.0\mathrm{km}$ and $c\simeq 340\mathrm{m}{\mathrm{s}}^{-1}$. These are representative atmospheric values, not additional data specified by the paper. They giveThe isothermal atmosphere approximation is local: a real atmospheric temperature profile and the thermodynamic response of the disturbance modify the precise cutoff and period.

Write the initial mass density as $\rho =mn$. The ideal gas law gives pressure $p=n{k}_{B}T$, and the sphere has $M=4\pi \rho R^3/3$. Here $k_B$ is the Boltzmann constant and $G$ is the gravitational constant. First obtain the gravitational potential energy by adding concentric shells. The mass inside radius $r$ is $M(r)=M{r}^3/R^3$, soIn a small contraction $dR<0$ at fixed total mass, the gravitational energy released is $-dW=3G{M}^2|dR|/(5R^2)$. The positive work required against the gas pressure is $-pdV=4\pi R^{2}p|dR|$. Equating these givesUse $R=(3M/(4\pi mn){)}^{1/3}$ to eliminate the radius. The uniform-sphere pressure-work collapse criterion isThe units are kilograms: $(k_{B}T/(Gm){)}^{3/2}$ has units $\mathrm{k}{\mathrm{g}}^{3/2}{\mathrm{m}}^{-3/2}$, and ${\rho }^{-1/2}$ supplies $\mathrm{k}{\mathrm{g}}^{-1/2}{\mathrm{m}}^{3/2}$. At fixed initial $n$ and $T$, $M>M_{\mathrm{crit}}$ makes gravity stronger than the pressure-work requirement in this model, allowing contraction.

This is a Jeans mass estimate for homologous contraction of an initially uniform sphere. It is not an exact stability limit for a hydrostatic gas sphere: uniform mass density and uniform temperature would give uniform pressure, rather than the gradient required by hydrostatic equilibrium. The numerical coefficient therefore differs from conventions based on a Jeans instability wavelength or a pressure-confined isothermal sphere. Subsequent contraction also depends on heat exchange; isothermal pressure support during gravitational collapse weakens relative to gravity as the radius decreases.

The background and the Jeans prescription. For a resting fluid, ${\mathbf{u}}_0=0$ and continuity is automatic. The remaining hydrostatic equilibrium equations areAn infinite uniform static medium with constant $p_0$ and ${\rho }_0>0$ would require $\nabla {\psi }_0=0$ by force balance, but then its Laplacian would vanish, contradicting the Poisson equation. An isolated static uniform background is therefore not an exact equilibrium of this system.

The Jeans swindle suppresses the gravitational field of the uniform background and retains only the field sourced by perturbations. Formally it replaces the source by $\rho -{\rho }_0$ in the perturbation problem. This can be regarded as a specified compensating background or a local idealization; it does not prove that the original infinite isolated medium can be static.

Linearized equations. Put $\delta \rho =\Delta \rho$, $\delta p=\Delta p$, $\delta \psi =\Delta \psi$, and $\mathbf{v}=\Delta \mathbf{u}$. For uniform ${\rho }_0,p_0$, neglect products of perturbations. Continuity givesThe advective acceleration is quadratic, and the background pressure gradient vanishes, so the Euler equation becomesWith the Jeans prescription and constant sound speed, the gravitational and pressure closures areThe density correction to $1/\rho$ multiplying $\nabla \delta p$ is second order and is consistently omitted.

The longitudinal dispersion relation. Differentiate continuity in time and substitute the divergence of the linearized Euler equation:For a plane wave proportional to $e^{i(\mathbf{k}\cdot \mathbf{x}-\omega t)}$, ${\partial }_t^2\mapsto -{\omega }^2$ and ${\nabla }^2\mapsto -k^2$. ThereforeThe physical Jeans length isSince $G{\rho }_0$ has units inverse square seconds, $k_J$ has units ${\mathrm{m}}^{-1}$ and ${\ell }_J$ metres. This is the compressive density branch of the Jeans instability; transverse velocity disturbances have no pressure or gravity restoring force in this static nonrotating model and are neutral at linear order.

Which scales grow. For $0<k<k_J$, or wavelength $\ell >{\ell }_J$, ${\omega }^2<0$ and one temporal solution grows exponentially. The static Jeans growth-rate spectrum isIts derivative is negative for $0<k<k_J$, so longer wavelengths grow faster. The supremum is $s\to \sqrt{4\pi G{\rho }_0}$ as $k\to 0^{+}$: there is no selected finite fastest wavelength in this infinite-medium calculation. The exactly uniform $k=0$ density change is not an ordinary finite-potential Fourier perturbation with the fixed background subtraction. In a finite cloud, boundaries set the longest allowed scale.

For $k>k_J$, the density perturbation oscillates as a stable acoustic wave. At $k=k_J$ the linear equation is ${\partial }_t^2\delta \rho =0$, so the mode is marginal rather than exponentially growing.

Comparing acoustic and gravitational times. A sound-crossing estimate for a region of size $\ell$ is $t_{\mathrm{sound}}\sim \ell /c_s$. Self-gravitational acceleration scales as $G{\rho }_0\ell$, so collapse over a distance $\ell$ takes $t_{\mathrm{grav}}\sim (G{\rho }_0{)}^{-1/2}$. Equating these scales giveswhich agrees with the Jeans length up to its numerical factor $\sqrt{\pi }$. This is the physical competition between pressure communication and gravitational contraction.

The constants can be checked using the free-fall time of a uniform sphere. For a pressureless shell initially at radius $a$, enclosing fixed mass $M_a$, energy conservation gives ${\dot{r}}^2=2G{M}_a(1/r-1/a)$. Thuswhere $M_a=4\pi {\rho }_0{a}^3/3$. Consequently the Jeans length and gravitational time conventions satisfyThey agree as characteristic scales, not as identical exact times. Alternatively, using the length $k_J^{-1}$ gives $1/(c_s{k}_J)=1/\sqrt{4\pi G{\rho }_0}$, exactly the gravitational inverse frequency entering the dispersion relation. The actual growth time of a marginal Jeans mode is infinite because pressure and gravity cancel there; it should not be confused with a pressureless free-fall time.

Following a parcel. Let $\mathbf{X}(t)$ be a fluid pathline, with $\dot{\mathbf{X}}=\mathbf{u}(\mathbf{X},t)$. The multivariable chain rule gives the material derivative,The first term measures change at fixed position; the advection term accounts for the parcel moving through spatial gradients. The derivative is also called the Lagrangian derivative.

The forces on a parcel. A material volume $dV$ has mass $\rho dV$. The pressure forces on its boundary sum to $-\nabla pdV$, while the gravitational acceleration is $-\nabla \psi$. Newton's second law therefore gives the Euler equations for an inviscid fluid,This neglects viscosity and other body forces. The left side is the parcel's acceleration, not just the time derivative at a fixed spatial location.

Bernoulli conservation. The claimed streamline result needs steady flow and a barotropic fluid, $p=p(\rho )$, or an equivalent enthalpy relation along each isentropic streamline. These hypotheses are necessary; the unsteady Euler equation alone does not imply the stated conservation law. Define the barotropic enthalpy functionDot the steady Euler equation with $\mathbf{u}$, using $\mathbf{u}\cdot [(\mathbf{u}\cdot \nabla )\mathbf{u}]=\mathbf{u}\cdot \nabla (u^2/2)$. This givesThus the Bernoulli function isEvery term has units of energy per mass, ${\mathrm{m}}^2{\mathrm{s}}^{-2}$. Different streamlines can have different constants.

The spherical accretion equation. Write the inward speed as $u(r)>0$, so the velocity vector is $-u(r)\hat{\mathbf{r}}$ and the positive mass accretion rate isAssume a prescribed central potential $\psi =-GM/r$ and negligible gas self-gravity. The steady radial Euler equation isFor a barotropic equation of state, $dp/dr=c_s^{2}d\rho /dr$, with $c_s^2=dp/d\rho$. Eliminate the density derivative to obtainA smooth sonic point requires the right side to vanish where $u=c_s$, henceFor a nonisothermal gas this is an implicit relation involving the local sound speed. In the isothermal case $c_s$ is constant and this is the Bondi sonic point.

The isothermal integral. For $p=c_s^2\rho$, choose the enthalpy reference so that $h_b=c_s^2\ln (\rho /{\rho }_{\infty })$. The reservoir is at rest at infinity and $\psi (\infty )=0$, making its Bernoulli constant zero. ThusorThe logarithm has a dimensionless argument. Constant temperature during compression generally requires heat removal, so this is Isothermal Bondi accretion, rather than adiabatic accretion with $\gamma >1$.

The critical accretion rate. Select the smooth reservoir-connected transonic branch. Evaluate the last relation at $u_s=c_s$ and $GM/r_s=2c_s^2$:Mass conservation then gives the isothermal Bondi accretion rate,The units are $\mathrm{kg}{\mathrm{s}}^{-1}$. Equivalently $\dot{M}=4\pi \lambda G^2{M}^2{\rho }_{\infty }/c_s^3$, with $\lambda =e^{3/2}/4\simeq 1.12$.

The inner-boundary selection of isothermal Bondi flow is important: a standard absorbing accretor inside the sonic radius selects this critical solution. The steady equations and ${\rho }_{\infty }$ alone also admit subcritical subsonic branches supported by a different inner boundary; they do not by themselves fix a unique rate. The usual answer presumes transonic capture.

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.