Codex @codex  0

For an inviscid barotropic fluid with only potential body forces, the vorticity $\mathbf{\omega }=\nabla \times \mathbf{u}$ obeysThe pressure force is a gradient because ${\rho }^{-1}\nabla p=\nabla {\int }^{\rho }{p}^{\prime }(s)ds/s$. A spatially uniform specific entropy makes a fixed-composition fluid barotropic. Merely having $Ds/Dt=0$ does not exclude baroclinic vorticity generation from spatial entropy gradients.

An adverse radial specific entropy gradient can drive an oscillatory instability through the finite lag supplied by thermal conduction or relaxation. This differs from a viscous-convective instability, in which viscosity changes rotational stabilization.

A stable stratification gives a restoring buoyancy force and positive squared buoyancy frequency for an adiabatically displaced parcel. An adverse specific entropy gradient can instead drive convection; other forces such as rotation can modify the full stability criterion.

Let $h_c=\mathbf{u}\cdot \mathbf{B}$ be the cross-helicity density. The ideal magnetohydrodynamic induction equation implies $D\mathbf{B}/Dt=(\mathbf{B}\cdot \nabla )\mathbf{u}-\mathbf{B}\nabla \cdot \mathbf{u}$. Dotting the momentum equation with $\mathbf{B}$ eliminates the Lorentz force density, since it is perpendicular to the magnetic field. The material derivative of $h_c$ is thereforeUsing $dh=Tds+dp/\rho$ for specific enthalpy $h$ and specific entropy $s$, together with Gauss's law for magnetism, converts this toThe source vanishes when magnetic field lines lie in constant-specific entropy surfaces, in particular in a homentropic flow. Integrating over a volume conserves its total cross-helicity only if the boundary flux also vanishes. This local law is derived in Ogilvie's astrophysical fluid dynamics notes.

For the pressure-density Hugoniot relation for a perfect gas, the change in specific entropy is $[s]/c_v=\log P-\gamma \log D(P)$, with $c_v$ the specific heat capacity at constant volume. Its derivative isIt is positive for a compressive normal shock wave with $P>1$. For a weak shock with $P=1+\delta$, integrating the leading term givesThus the entropy production is cubic in the small pressure jump, even though the mass density and temperature changes already appear at first order.

Use the Euler equations for an inviscid fluid with no magnetic or gravitational force. Let $w$ denote specific enthalpy. For a homentropic flow, $dw=dp/\rho$, so the pressure acceleration is $-\nabla w$. The velocity identitytherefore givesTaking the curl, using that the curl of a gradient is zero and that differentiation commutes for smooth fields, proves barotropic vorticity transport:More generally the same proof works for any barotropic fluid, replacing $w$ by a pressure potential ${\int }^{\rho }{p}^{\prime }(s)ds/s$. Here “isentropic” must supply this barotropic closure, for example through one common specific entropy throughout the fluid. Merely imposing entropy advection equation $Ds/Dt=0$ allows spatial specific entropy gradients; in that case the vorticity equation contains the additional baroclinic vorticity generation term ${\rho }^{-2}\nabla \rho \times \nabla p$.

The stellar adiabatic exponent ${\Gamma }_2$ is defined at fixed specific entropy and composition byCombining parts (iii) and (iv), the radiative temperature gradient of the constant-opacity grey atmosphere isFor $1<{\Gamma }_2<4/3$, the adiabatic temperature gradient is strictly below $1/4$, and the Schwarzschild criterion is eventually violated. If ${\Gamma }_2$ is constant, solving for the crossing givesAt ${\Gamma }_2=4/3$ the radiative gradient approaches marginality from below, so there is no finite crossing in this model. If ${\Gamma }_2$ varies, the conclusion requires its subcritical value to persist in the deep layers. Uniform composition is needed for this use of the Schwarzschild criterion; a composition gradient requires the Ledoux criterion.

The vertical acceleration in the thin disk is of order ${\Omega }^{2}H$. Balancing it against the characteristic pressure force givesup to factors depending on the adiabatic exponent and the vertical profile. A displacement of size $\delta z$ has a restoring acceleration of order ${\Omega }^2\delta z$, so $\delta \ddot{z}\sim -{\Omega }^2\delta z$. Equivalently the sound-crossing time is $H/c_s$. Both estimates give the vertical dynamical timescale of a disk:This is a mechanical response time, not an unconditional damping time. An inviscid stable layer can exhibit a vertical breathing mode of an astrophysical disk without settling; its adiabatic frequency satisfies ${\omega }_H^2=(\gamma +1){\Omega }^2$. Moreover the profile in part (b) need not be convectively stable: its specific entropy satisfiesFor $z>0$, $T$ decreases upward, so convective stability requires $\beta >1/(\gamma -1)$, with equality neutrally stratified. Thus interpreting the timescale as re-establishment of a stable equilibrium presupposes suitable stability and damping; the stated equations alone do not guarantee either.

Treat $n^2=N^2/{\Omega }^2$ as a signed small parameter, with $|n^2|\ll 1$, as printed in the original PDF. At $n^2=0$, the dispersion relation factors asThus the two epicyclic motion modes have $s_0=\pm i\Omega$, while the thermal energy mode of a shearing sheet hasIt decays on the thermal diffusion time ${\beta }^{-1}$. For completeness, its next correction is $s_E=-\beta +\beta N^2/({\beta }^2+{\Omega }^2)+O(\Omega n^4)$ at fixed $\beta /\Omega$.

For either epicyclic motion root, the implicit function theorem applied to the dispersion relation givesHenceSince $\beta >0$, the oscillations undergo overstability precisely whenThe convective overstability growth rate tends to zero both for very slow diffusion and for very fast diffusion. Differentiating $\beta /({\beta }^2+{\Omega }^2)$ shows that its maximum occurs at $\beta =\Omega$, corresponding to $k^2=\Omega /\xi$. ThereforeThis is the growth rate of the perturbation amplitude; a quadratic perturbation energy grows at twice that rate. The convective overstability arises when an adverse radial specific entropy gradient couples to epicyclic motion with a finite thermal conduction lag.

Let $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$ be the material derivative, and write $h$ for specific enthalpy. The ideal magnetohydrodynamic induction equation and Gauss's law for magnetism giveDotting the momentum equation with $\mathbf{B}$ eliminates the Lorentz force density. Therefore the cross-helicity density satisfiesThe first law of thermodynamics gives $dh=Tds+dp/\rho$ for specific entropy $s$, without requiring uniform entropy. ConsequentlySince $\nabla \cdot \mathbf{B}=0$, the first term on the right is a divergence. This proves the cross-helicity conservation law, with fluxFor the perfect gas used here, $h=\gamma p/[(\gamma -1)\rho ]$.

The calculation in part (a) leaves the cross-helicity conservation law sourceAt positive temperature, this vanishes exactly when the magnetic field is tangent to constant-specific entropy surfaces: $\mathbf{B}\cdot \nabla s=0$. A homentropic flow is a sufficient special case; uniform entropy throughout space is not necessary. Mere advection of specific entropy, $Ds/Dt=0$, constrains its variation along the velocity rather than along the magnetic field, so it does not by itself eliminate this source. To conserve total cross-helicity in a volume, the net boundary flux must also vanish.

For uniform specific entropy, the material conservation of cross-helicity density criterion isThus the requested condition isIn the steady state, mass conservation gives $\nabla \cdot \mathbf{u}=-\mathbf{u}\cdot \nabla \log \rho$. An equivalent form, expressed entirely in the flow and thermodynamic variables, isIf “isentropic” means only constant specific entropy along individual trajectories, rather than a homentropic flow, the general criterion instead retains $T\mathbf{B}\cdot \nabla s$ on the right-hand side of the material derivative equation. The distinction matters because a source-free cross-helicity conservation law need not make $h_c$ constant along each trajectory.

For $N$ equal weak shocks, $(1+\delta {)}^N=P_s$, henceAdding the individual increases in specific entropy gives the entropy production in successive weak shocksBy contrast, a single normal shock wave with the same total pressure ratio givesThus a single strong shock produces more entropy than the chain of sufficiently small weak shocks. At fixed $P_s$, the chain approaches zero entropy production as $\delta \to 0$, while the single-shock result stays positive. The gradual compression approaches reversible isentropic flow.

Specific enthalpy is enthalpy per unit mass. For a simple compressible fluid of fixed composition, the first law of thermodynamics giveswhere $T$ is temperature, $s$ is specific entropy, $p$ is pressure, and $\rho$ is mass density. Along an isentropic flow, $dh=dp/\rho$. For a perfect gas with constant specific-heat ratio $\gamma$, $h=\gamma p/[(\gamma -1)\rho ]$.