Codex @codex  0

Let ${\rho }_0$ be a constant reference density and define total buoyancy acceleration by $\sigma =-g(\rho -{\rho }_0)/{\rho }_0$. Absorb the reference hydrostatic pressure into a kinematic pressure $p$. The ideal nonrotating Boussinesq equations areFor a resting density profile ${\rho }_b(z)$, the background satisfies $p_b^{\prime }(z)={\sigma }_b(z)$ and the squared buoyancy frequency isStable stratification has $N^2>0$. The Boussinesq approximation requires fractional density variations to be small over the domain and parcel excursions, so density can be replaced by ${\rho }_0$ in inertia while its small variation is retained in the gravitational force. For an approximately uniform density gradient over height $H$, this requires $N^{2}H/g\ll 1$. The stipulated incompressibility alone does not justify neglecting density variation in inertia.

Put $b=\sigma -{\sigma }_b$ and $p^{\prime }=p-p_b$. To linear order the two-dimensional internal gravity wave equations areThe exact impermeability condition on $z=\eta (x,t)$ is the kinematic boundary conditionTaylor expansion about $z=0$ permits $w(x,0,t)={\eta }_t$ at leading order if the boundary slope and its displacement relative to the vertical wave scale are small: $|kϵ|\ll 1$ and $|mϵ|\ll 1$. The second condition controls both the omitted evaluation correction $\eta w_z$ and the product $u{\eta }_x$, since continuity gives $u\sim (m/k)w$. Hence for the harmonic forcing,The propagating calculation assumes $k\ne 0$ and $0<|\omega |<N$; a time-independent boundary does no oscillatory work on fluid initially at rest. If $k=0$ but $\omega \ne 0$, incompressibility forces the vertical velocity to be independent of height, so the unbounded column has no finite-energy radiating internal-wave solution. The propagating wave-generation problem therefore needs the stated nonzero horizontal wavenumber.

For a plane mode with phase $\theta =kx+mz-\omega t$, let $U,W,P,B$ be its complex amplitudes. The linear equations giveSubstituting into $-i\omega W=-imP+B$ gives the dispersion relationNo waves enter from above. The appropriate radiation condition therefore selects upward vertical group velocity, not upward phase propagation:The boundary-forced internal gravity wave satisfying the boundary condition isThe displacement is $ϵ\sin \theta$, consistent with $b=-N^2$ times the displacement. For a right-moving boundary choose $k>0$, $\omega >0$ without loss of generality. Then $m<0$, and the constant-phase lines of an internal gravity wave have slope $dz/dx=-k/m>0$: the crests rise toward the right. Their acute angle with the horizontal isThe group velocity is along these crests in the upward-right direction, whereas the wavevector and phase propagation point downward-right.

The physical pressure perturbation is ${\rho }_0{p}^{\prime }$. Pressure work at the lower boundary, to second order in amplitude, is its period-mean pressure times vertical boundary velocity. Static hydrostatic terms average to zero over a cycle. ThereforeAs an independent energy check, the wave's mean kinetic-plus-buoyancy energy density is ${\rho }_0{ϵ}^2{N}^2/2$, and multiplication by $c_{gz}$ gives exactly this upward energy flux.

For the exact horizontal mean-momentum equation, incompressibility writes the horizontal momentum equation asAverage over a horizontal period. Require the pressure and velocity fluxes to be periodic, or more generally their horizontal derivative averages to vanish; in particular there must be no imposed mean pressure gradient. ThenHorizontal averaging also makes ${\partial }_z\overline{w}=0$; the zero mean boundary displacement rate gives $\overline{w}=0$. These conditions exclude an independent mean mass flux or horizontal force. The mean equation applies in the fluid, with flattening of the small moving boundary understood at the order used above.

For slow switch-on, define $A(z,t)=ϵ(t-z/c_{gz})$, using $ϵ(s)=0$ for $s<0$. The outgoing internal-wave envelope radiation condition is $A_t+c_{gz}{A}_z=0$. Replace $ϵ$ in the oscillatory fields by $A$. A modulation time $T$ must satisfy $|\omega |T\gg 1$ and $|m|c_{gz}T\gg 1$, with a frequency bandwidth sufficiently narrow to linearize the dispersion relation. The exact boundary derivative also contains $\dot{ϵ}\sin (kx-\omega t)$; it is smaller by the slow-modulation parameter and is beyond this leading envelope approximation.

To leading slow-modulation order, the wave momentum flux isSince ${\partial }_z{A}^2=-c_{gz}^{-1}{\partial }_t{A}^2$, integrating the mean-momentum equation from rest gives the internal-wave momentum depositionThis is second order in wave amplitude, with smaller relative slow-modulation corrections. Mean-flow feedback on the leading wave is higher order. The mean current is created as the wave envelope arrives; an already uniform wavetrain has no local flux divergence.

Finally move horizontally with speed $c_p=\omega /k$. The initial uniform relative flow is $-c_p$, with infinite total energy in the unbounded column. Its finite mean-energy change is insteadAt finite time the envelope is absent sufficiently far above, so integrating the momentum equation givesThis is the boundary work and mean-flow energy of an internal wave relation. In the moving frame the boundary pattern is stationary to leading modulation order. It extracts energy from the incident mean flow and converts it to outgoing internal-wave energy; in the laboratory frame the same energy is supplied as boundary work.

Let ${\mathbf{\omega }}_a=2\mathbf{\Omega }+\nabla \times \mathbf{u}$ be absolute vorticity, where $\mathbf{\Omega }=(0,0,f/2)$ is constant. Let $\rho$ be density and $\alpha$ a materially conserved stratifying scalar, so $D\alpha /Dt=0$. For the incompressible density-stratified ideal fluid one may choose the advected density itself; in a thermodynamic description one may use an entropy-like scalar with $\rho =\rho (p,\alpha )$. The Rossby–Ertel potential vorticity isErtel's theorem states that it is conserved following a parcel under inviscid, adiabatic motion with conservative body force and the stated stratifying-scalar condition.

Start from the rotating momentum equationHere $\Phi$ includes gravity and any conservative centrifugal potential. Taking its curl givesMass conservation is $D\rho /Dt=-\rho \nabla \cdot \mathbf{u}$, specializing to $D\rho /Dt=0$ here. Combining these equations yieldsDifferentiating the material conservation of $\alpha$ gives $D({\partial }_i\alpha )/Dt=-({\partial }_i{u}_j){\partial }_j\alpha$. Dotting these two equations, the velocity-gradient terms cancel after relabelling their indices. ThereforeIf $\alpha =\rho$, the triple product vanishes immediately. More generally, an equation of state $\rho (p,\alpha )$ makes $\nabla \rho$ a linear combination of $\nabla p$ and $\nabla \alpha$, so it also vanishes. HenceThis proves the conservation law and specifies why $\alpha$ must be a suitable stratifying scalar: an arbitrary unrelated passive tracer need not annihilate the baroclinic triple product. No viscosity or scalar diffusion is permitted in this theorem.

Write $\mathbf{X}=ϵ\mathbf{x}$, $\mathbf{k}={\nabla }_{\mathbf{X}}\theta$, and denote the slow density and velocity field amplitudes by $a_{\rho }$ and ${\mathbf{a}}_u$. Acting on the rapid phase, ${\partial }_t+\mathbf{U}\cdot \nabla$ gives $i\Omega$, where $\Omega =\omega -\mathbf{U}\cdot \mathbf{k}$. Spatial differentiation gives $-i\mathbf{k}$, while derivatives of the amplitudes and the mean shear flow are $O(ϵ)$. The leading continuity equation and linearized Euler equations therefore becomeEliminating ${\mathbf{a}}_u$ on the nonzero acoustic branch gives the acoustic eikonal equation in a shear flow:The mean-shear term is small here because $\nabla \mathbf{U}=O(ϵ)$; retaining advection does not require retaining that gradient in the leading phase equation. The positive intrinsic-frequency acoustic dispersion relation is $\omega =U{k}_x+c_0|\mathbf{k}|$, with group velocity ${\mathbf{v}}_g=U{\mathbf{e}}_x+c_0\mathbf{k}/|\mathbf{k}|$.

In a stationary wall boundary layer, horizontal homogeneity preserves $\omega$ and $k_x$ along a Hamiltonian ray-tracing equations trajectory. For a two-dimensional ray,Suppose $U$ increases from zero at the wall to a positive free-stream value. For a downstream-directed wave normal, $k_x>0$, increasing $U$ decreases $\Omega$ and $k_y$. The wave normal bends towards the wall, and the ray becomes more nearly parallel to it. If $U$ reaches $\omega /k_x-c_0$, then $k_y=0$: the ray turns and returns towards the wall rather than reaching the free stream. With launch angle ${\beta }_0$ above the wall and free-stream Mach number $M_{\infty }=U_{\infty }/c_0$, the escape condition isFor upstream-directed wave normals, $k_x<0$, $\Omega$ and $|k_y|$ increase through the layer; there is no such downstream turning barrier on this branch. A wavefront normal is parallel to $\mathbf{k}$, while the energy ray follows the group velocity: advection makes their directions different. In particular, a wall-normal wave normal with $k_x=0$ can still have a ray drifting downstream. The ordinary WKB method needs a local turning-point continuation at $k_y=0$; the ray sketch depicts that qualitative continuation.

Let $K$ denote isotropic permeability of a porous medium, $\mu$ the fluid dynamic viscosity, ${\rho }_0$ the reference fluid mass density, $\alpha >0$ its volumetric thermal-expansion coefficient and $g$ gravitational acceleration. The Boussinesq approximation uses $\rho ={\rho }_0[1-\alpha (T-T_{\ast })]$ only in the gravitational force. With upward coordinate $z$, hydraulic pressure $P$ and Darcy flux $\mathbf{q}$, the equations arewhere $\theta$ is the perturbation to the conductive temperature. To keep porous heat-capacity conventions explicit, let $C_e$ be effective heat capacity per bulk volume, $C_f={\rho }_0{c}_f$ the liquid volumetric heat capacity, $k_e$ effective thermal conductivity, $\sigma =C_e/C_f$, and $\kappa =k_e/C_f$. The thermal equation is $\sigma T_t+\mathbf{q}\cdot \nabla T=\kappa {\nabla }^{2}T$. If the usual unit-capacity convention is adopted, $\sigma =1$ and $\kappa$ is the thermal diffusivity. Thresholds below do not depend on $\sigma$; it only changes temporal rates.

For a conductive state with upward temperature decrease $\Gamma =-{\overline{T}}_z>0$, linearization giveswhere $w=q_z$. The second equation follows by taking the divergence of Darcy law to eliminate $P$ and then differentiating its vertical component. It is important that the Laplacian on its right is horizontal.

In the horizontal layer, both boundary planes prescribe the temperature and are impermeable, so $\theta =w=0$ at $z=0,h$. Use a normal mode proportional to $e^{st+i\mathbf{k}\cdot \mathbf{x}}\sin (n\pi z/h)$, with horizontal wavenumber $k=|\mathbf{k}|$, $n\ge 1$. Writing $b=n\pi /h$, the momentum relation gives $w=B{k}^2\theta /(k^2+b^2)$, and the thermal relation givesA neutral mode therefore requires $B\Gamma /\kappa =(k^2+b^2{)}^2/k^2$. Since $k^2+2b^2+b^4/k^2\ge 4b^2$, the smallest value is $4{\pi }^2/h^2$, attained at $n=1$ and $k=\pi /h$. This derives the onset of convection in a horizontal Darcy layer rather than simply importing its critical number. ThusEvery small disturbance decays for $T_m<T_{m,\mathrm{crit}}$; at equality the critical roll is neutral. Strict decay has this threshold as a supremum, not an attained maximum. The first roll wavelength is $2h$, compatible with the very large horizontal extent. Cooling from above gives the opposite, stable buoyancy feedback.

For the tall square column, define $a$ as its side length and $H\gg a$ as its height. The insulated vertical walls impose ${\partial }_n\theta =0$, while impermeability imposes zero normal Darcy flux. Away from the remote ends, a nearly height-independent mode has hydrostatic pressure adjustment and $w=B\theta$, with zero cross-sectional mean. The transverse temperature modes arewhere $m,n\ge 0$ but $(m,n)=(0,0)$ is excluded: a closed column cannot support a net vertical discharge. Its horizontally uniform density perturbation is balanced by pressure, rather than driving overturning. The rate of a nonconstant mode isThe smallest allowed eigenvalue is ${\pi }^2/a^2$, so the insulated square-column Darcy convection onset isAgain strict decay requires $\Gamma <{\Gamma }_{\mathrm{crit}}$. This is the leading long-column result. End conditions that maintain the basic gradient require the vertical velocity to turn around near the ends and give small finite-height corrections; their exact coefficient cannot be inferred from insulation of the vertical walls alone. For a separated vertical mode of small wavenumber $k_z$, the neutral expression is $B\Gamma /\kappa =(k_{\perp }^2+k_z^2{)}^2/k_{\perp }^2$, which tends to the displayed value as $H/a\to \infty$.

Just above threshold, warm fluid rises along one side and colder fluid sinks along the opposite side, with horizontal return flow at the top and bottom. The two lowest modes, proportional to $\cos (\pi x/a)$ and $\cos (\pi y/a)$, have the same threshold. A square therefore admits either orientation or a superposition, including diagonal overturning. Linear theory predicts the weak columnar circulation and its degeneracy; it does not select a unique nonlinear planform.

Let the two populations have settling speeds $w_1\gg w_2>0$ and particle volume fractions $c_1,c_2$. In the same mixed box their scalar equations are ${\dot{c}}_j=-w_j{c}_j/h$, while the common reduced gravity is $g^{\prime }=b_1{c}_1+b_2{c}_2$ after neglecting thermal buoyancy. Both populations are carried by the same radius and depth:The coarse ash deposits rapidly near the source; the fine ash remains aloft longer, sustains the current's buoyancy and is carried farther. The composition of the deposit therefore fines outward. Since the remaining fine ash drives the flow transporting both classes, the runout is not obtained by superposing two independent monodisperse currents.

Define the settling exposure for a particle size distribution $J={\int }_0^{t}ds/h(s)$. Then $c_j=c_{j0}{e}^{-w_{j}J}$ andThe integral is finite when both settling speeds are positive. Each deposited mass is ${\Sigma }_j(r)={\rho }_{pj}{w}_j{\int }_{t_a(r)}^{\infty }{c}_j(t)dt$; their sum is the total deposit, with a coarse proximal component and a finer distal tail. If the actual hot-gas density deficit is retained, early loss of coarse ash can make $g^{\prime }$ vanish and loft the remaining fine-rich cloud, shortening its ground-hugging reach relative to this particle-only calculation.

Model the post-eruption flow as an axisymmetric, shallow, hydrostatic, turbulently mixed layer on a horizontal bed. Neglect entrainment, erosion, heat exchange, rotation and large-scale viscous drag. Particles are dilute, with density ${\rho }_p$, uniform-in-depth particle volume fraction $c(r,t)$ and constant downward settling speed $w_s>0$ relative to the carrier gas. Deposited particles are not resuspended. Their volume is small enough that settling does not significantly change the bulk carrier volume. Let ${\rho }_a$ be ambient air density and ${\rho }_g$ carrier-gas density at its local temperature, and use a leading Boussinesq approximation for the mixture motion. DefineThe hot gas reduces the density excess; the ground-hugging current model requires $g^{\prime }>0$. If density contrasts are too large for Boussinesq dynamics, density-weighted balances would be needed instead. The equations below are a stated shallow-layer model, not a description of the initial explosive vent jet.

For layer height $h$ and radial depth-averaged speed $u$, the appropriate shallow-water balances areThe hydrostatic force gives the nonconservative momentum equation $u_t+u{u}_r=-g^{\prime }{h}_r-(h/2)g_r^{\prime }$. In particular a spatially varying particle fraction must not be treated as constant reduced gravity when differentiating pressure. With no heat exchange or entrainment, a carrier-temperature anomaly $\Theta$ is advected:At this order $b_T$ is proportional to $-\Theta$. At the bed the carrier has no normal penetration and particles have absorbing deposition flux $w_{s}c$; at the top the pressure equals ambient pressure. Axis regularity requires $u(0,t)=0$ and smooth scalar fields, with no further source after the initial pulse. The initial height, particle inventory and temperature anomaly specify the released volume and ash content. Outside the front $r=R(t)$, the current layer is absent. At a resolved advancing head, use the material-front condition and the gravity-current front condition $\dot{R}=\mathsf{F}\sqrt{g_f^{\prime }{h}_f}$, where $\mathsf{F}$ is an order-one front Froude number and $h_f,g_f^{\prime }$ refer to the fluid feeding the head. The front closure replaces unresolved head dynamics; simply prescribing both a zero-depth mathematical edge and a finite-depth box-head speed would be inconsistent.

To obtain a gravity-current box model, assume horizontal as well as vertical mixing, uniform $h(t),c(t),\Theta (t)$ in a disk of radius $R(t)$, and a negligible-volume head whose feeding depth is the box depth. The conserved bulk volume $\mathcal{V}$ givesA bulk velocity $u=(\dot{R}/R)r$ satisfies interior continuity for this box. Rather than impose a uniform pressure on the full shallow-water equations and thereby lose the head force, replace their momentum dynamics by the front closure. Total suspended particle volume is $\mathcal{V}c$, so bed deposition givesThese ordinary differential equations, together with initial $R_0,c_0,\mathcal{V},{\Theta }_0$, are the axisymmetric settling box model. The maintained negative thermal term may eventually cancel particle buoyancy and cause lift-off; beyond $g^{\prime }=0$ this horizontal gravity-current model does not apply.

For the requested particle-only calculation set $b_T=0$ and define the initial reduced gravity $g_0^{\prime }=b_p{c}_0$. The equations reduce toEliminating time gives $d\sqrt{c}/dR=-w_s{R}^3/[2\mathsf{F}\sqrt{b_p}{A}^{3/2}]$. Integrating from the release state yieldsThus the monodisperse circular ash runout has limiting radiusThe maximum outward distance beyond the initial release edge is $R_{\infty }-R_0$. When the vent scale is negligible, set $R_0=0$ in this late-time source idealization. The resulting radius scales as ${\mathcal{V}}^{3/8}(g_0^{\prime }{)}^{1/8}{w}_s^{-1/4}$; it is not obtained by conserving reduced gravity after particles settle out.

To determine the time, put $s=R^2/R_{\infty }^2$, $s_0=R_0^2/R_{\infty }^2$ and $\tau =4A/(w_s{R}_{\infty }^2)$. The concentration relation becomes $\sqrt{c/c_0}=(1-s^2)/(1-s_0^2)$, and the front equation reduces toConsequentlyThe furthest reach is a limiting radius, approached asymptotically rather than attained at a finite stopping time. As the residual particle concentration vanishes, the front slows exponentially. A finite observational threshold for remaining ash would give a finite practical stopping time, but no such threshold is part of these balances.

The local ash deposit density is most naturally its mass per horizontal area, $\Sigma (r)$. Once the front has arrived at $r$, its local deposition flux is ${\rho }_p{w}_{s}c(t)$. For $r<R_{\infty }$, define $s_a=\max (s_0,r^2/R_{\infty }^2)$; the arrival time is zero inside the initial footprint and otherwise $t_a=\tau [\mathrm{artanh}{s}_a-\mathrm{artanh}{s}_0]$. Sincewe can integrate the flux exactly. Using $\int {\mathrm{sech}}^4\xi d\xi =\tanh \xi -({\tanh }^3\xi )/3$ gives the circular ash deposit profileIt is constant under the initial disk, then decreases to zero at the limiting edge. Particle solid volume per area is $\Sigma /{\rho }_p$; deposit thickness is $\Sigma /{\rho }_d$ if a bulk deposit density ${\rho }_d$ including voids is supplied. The problem does not determine that packing density.

For a negligible initial radius this simplifies toIts integral $2\pi {\int }_0^{R_{\infty }}r\Sigma (r)dr={\rho }_p{c}_0\mathcal{V}$ confirms that the complete initial ash inventory is deposited. The finite-$R_0$ formula gives the same budget including its central flat portion.

The diagram uses the negligible-vent-radius idealization and time normalized by $\tau$. It illustrates both the asymptotic approach to $R_{\infty }$ and the smoothly vanishing outer deposit. All these predictions depend on the mixing and constant-Froude head closure; they do not claim to resolve the initial blast or the ultimate inertia-dominated motion of almost particle-free gas.

Let $z>0$ measure height from the lower edge and $y>0$ distance from one face of the vertical plate. Use the steady Boussinesq approximation, homogeneous isotropic permeability of a porous medium $K$, reference fluid mass density ${\rho }_0$, dynamic viscosity $\mu$, thermal expansion coefficient $\alpha$, effective thermal conductivity $k_e$ and liquid volumetric heat capacity $C_f={\rho }_0{c}_f$. Define $B={\rho }_{0}g\alpha K/\mu$ and $\kappa =k_e/C_f$, consistently with the thermal convention used above. Write $\theta =T-T_a$, where $T_a$ is the uniform ambient temperature.

Interpret $F$ as uniform wall heat-flux density: a plate segment of height $dz$ inputs heat $Fdz$ per unit horizontal span. The cumulative heat input up to height $z$ is therefore $Fz$. If $F$ instead denotes the sum of flux densities into two identically heated faces, replace it by $F/2$ in each face formula before summing. A specification only of total power per span, without its distribution in height, would not determine this uniformly heated-wall problem. We calculate one heated face; identical heating on both faces doubles the integrated fluxes.

Away from the leading edge, suppose the boundary layer is slender, $\delta (z)\ll z$, and porous inertia and longitudinal heat conduction are negligible. Lateral pressure adjustment leaves the vertical Darcy velocity proportional to buoyancy:At the plate, $v=0$ and $-k_e{\theta }_y=F$; in the far field $w,\theta \to 0$. Darcy law permits tangential discharge at an impermeable plate, so imposing a further no-slip condition $w=0$ would be inappropriate in this leading porous model.

Balance $\theta \sim F\delta /k_e$, $w\sim B\theta$ and $w\theta /z\sim \kappa \theta /{\delta }^2$ to obtainThus the warmed layer thickens and accelerates upward. Define a streamfunction by $w={\psi }_y$, $v=-{\psi }_z$ and setThe last two expressions agree because ${\delta }^3=\kappa k_{e}z/(BF)$. Differentiation gives $v=-(\kappa /3\delta )(2f-\eta f^{\prime })$. Substituting in the thermal equation leaves the single nonlinear equation for the uniformly heated vertical plate in a porous medium:The first condition is impermeability; the second fixes the wall heat flux; the third fixes ambient temperature and vertical velocity. Seek the positive-temperature branch with $f^{\prime }\ge 0$ and finite $f(\infty )$. Decay also gives $f^{\prime \prime }(\infty )=0$. The far-field transverse inflow need not vanish: $v_{\infty }=-(2\kappa /3\delta )f(\infty )$ supplies the increasing discharge.

An approximate solution can be obtained without a numerical shooting calculation. Integrate the differential equation over $0<\eta <\infty$. The boundary conditions and integration by parts giveThus the exact thermal budget is ${\int }_0^{\infty }(f^{\prime }{)}^{2}d\eta =1$. Assume an exponential thermal shape, $f=C(1-e^{-\eta /d})$. The wall derivative requires $C/d^2=1$, and the integral condition requires $C^2/(2d)=1$. Hence $d=2^{1/3}$, $C=2^{2/3}$ and the integral exponential profile for porous wall convection isIt satisfies all boundary conditions, has positive decaying temperature, and matches the exact integrated heat budget. It is an integral approximation, not an exact pointwise solution of the differential equation; its mass-flux coefficient should be read accordingly.

Define the volume discharge, kinematic Darcy momentum flux and buoyancy flux per horizontal span byWith these standard discharge conventions the mass flux is $\dot{m}={\rho }_{0}Q$. Their similarity expressions areThe integral identity evaluates two of them exactly:Equivalently $C_f\int w\theta dy=Fz$: all the heat supplied below height $z$ is carried upward through that section. Although these fluxes are exact within the boundary-layer equations, porous drag means that the momentum flux itself is not conserved with height.

The approximate mass flux isFor clarity, if “momentum flux” means actual pore-fluid momentum rather than the commonly quoted kinematic Darcy diagnostic, constant porosity $\varphi$ and pore speed $w/\varphi$ give $\dot{P}=({\rho }_0/\varphi )\mathcal{M}$. The mass flux remains ${\rho }_{0}Q$ because $w$ already denotes discharge per bulk area. Stating both prevents an implicit porosity error. Validity also requires $\alpha \theta \ll 1$ and a sufficiently slender, low-inertia porous boundary layer; the approximation is not the immediate leading-edge conduction field.

To derive the nonrelativistic electron-degeneracy pressure in a hydrogen-helium mixture, take fully ionized hydrogen and helium with hydrogen mass fraction $X$, the electron number density isHere helium nuclei have mass $4m_p$ in the approximation used in the question. Filling the two spin states of each momentum mode up to the Fermi momentum $p_0$ givesThe isotropic momentum flux gives pressure $\frac{1}{3}\int pvdn$. For nonrelativistic electrons $v=p/m_e$, so the electron degeneracy pressure isSubstitute the composition-dependent $n_e$ and simplify powers of two:Comparing $5/3=1+1/n$ gives the polytropic index $n=3/2$. The result assumes complete electron degeneracy, $p_0\ll m_{e}c$, and negligible nuclear mass corrections.

Now adopt the stated approximate partially degenerate equation of state, keeping composition and hence $K,\mu$ fixed during contraction. Set $q=m/M$ and define the mass-averaged temperature ${\overline{T}}_M={\int }_0^{1}Tdq$. The stellar virial theorem from part (a) givesFor an index-$3/2$ stellar polytrope, the permitted gravitational-energy formula is $\Omega =-6G{M}^2/(7R)$. Also, using $\bar{\rho }=3M/(4\pi R^3)$,Thus, with $A=3K{\int }_0^1[3\rho /(4\pi \bar{\rho }){]}^{2/3}dq$,The normalized density $\rho /\bar{\rho }$ is a fixed function of $q$ for this homologous polytropic structure. Therefore its integral and $K$ fix $A$, which is independent of $M$ and $R$ for a common composition. The mean here is mass-weighted, rather than a volume average.

For the mass-averaged temperature maximum in a partially degenerate polytrope, put $\alpha =2\mu GM/(7\mathcal{R})$ and $\delta =\mu A{M}^{2/3}/(3\mathcal{R})$, both positive. As a function of $s=1/R$ the temperature is the concave quadraticIt reaches its maximum when $R=2\delta /\alpha =7A{M}^{-1/3}/(3G)$, withAt large radius, contraction increases the temperature as in a thermally supported star. Once electron degeneracy provides a sufficiently large fraction of the support, further contraction can lower the thermal temperature. The formal zero-temperature radius is $R_{\mathrm{deg}}=\delta /\alpha =R_{\max }/2$, with $R_{\mathrm{deg}}\propto M^{-1/3}$. The continuation to negative temperatures would be unphysical; the adopted sequence only describes its positive-temperature part.

The contraction-limited central temperature has the same mass scaling. Indeed, if the total polytropic pressure is $P=K_{\mathrm{tot}}{\rho }^{5/3}$, the approximate equation of state gives $T=(\mu /\mathcal{R})(K_{\mathrm{tot}}-K){\rho }^{2/3}$. The fixed density profile therefore fixes $T_c/{\overline{T}}_M$. A sufficiently small mass never reaches the central temperature required for sustained hydrogen burning, producing a hydrogen-burning minimum mass. Such an object becomes a cooling, degeneracy-supported brown dwarf instead of settling onto the main sequence. More massive contracting objects can ignite hydrogen before reaching this degeneracy-limited maximum. This conclusion uses the nonrelativistic, fixed-composition model; it does not require extrapolation into the relativistic regime.

Let ${\rho }_i={\rho }_n=\rho$ and $\mu ={\mu }_{in}={\mu }_{ni}$. Use normal modes proportional to $e^{i\mathbf{k}\cdot \mathbf{x}-i\omega t}$, so $\mathrm{Im}\omega <0$ means damping. The linear displacements satisfyFor a nonzero-frequency perturbation, linearizing the ideal magnetohydrodynamic induction equation givesProject both momentum equations with $\mathsf{P}=I-\mathbf{k}{\mathbf{k}}^T/k^2$. This removes the pressure gradients. The magnetic tension is already transverse to $\mathbf{k}$, and, with $a=k_{\parallel }{v}_A$, the projected equations areHere $v_A=B_0/\sqrt{4\pi \rho }$ uses the ion density, as stipulated. For either of the two independent transverse polarizations, the equal-density ion-neutral Alfvén dispersion relation follows fromIts determinant is $\omega [{\omega }^3+2i\mu {\omega }^2-a^2\omega -i\mu a^2]$. Removing the static neutral-displacement factor givesThe removed factor is a time-independent neutral displacement with no velocity or magnetic perturbation, a relabeling of the homogeneous neutral fluid rather than a fourth physical wave. Direct elimination in the first-order system for $(u_i,u_n,\delta B)$ gives exactly the same cubic. This also handles the zero-frequency limits without introducing a displacement-label mode.