Codex @codex  0

For a high-Reynolds number thermal, buoyancy force scales as ${\rho }_0{g}^{\prime }{a}^3$ and inertial resistance scales as ${\rho }_0{W}^2{a}^2$. Their balance gives $W^2\propto g^{\prime }a$. It is therefore natural to write a diameter-based Froude number closureThe order-one coefficient $F$ accounts for the thermal's shape, entrainment and flow structure. The main resistance is inertial pressure or form drag and the momentum needed to entrain and accelerate surrounding water, rather than molecular skin friction. Added mass matters when the thermal accelerates, so the Froude closure is a quasi-steady approximation.

In a homogeneous ocean all of the leading buoyancy comes from the bubbles. Parts (a) and (b) giveSubstituting into the speed law yieldsFor $z\gg a_0/\alpha$, this becomesThe square root is essential both physically and dimensionally. This approximation still requires dilute bubbles and retained gas; it is not a valid prediction of infinite speed where the pressure approximation approaches zero.

Let $E(k)$ be the one-dimensional turbulent energy spectrum, normalized so that $\frac{1}{2}⟨|{\mathbf{u}}^{\prime }{|}^2⟩={\int }_0^{\infty }E(k)dk$. Large energy-containing eddies receive energy from the forcing. In the inertial range, nonlinear transfer passes this energy towards smaller scales at a mean rate $ϵ$ per unit mass. Dimensional analysis of Kolmogorov 1941 theory givesThe dissipation range begins near the inverse of the Kolmogorov length scale, where kinematic viscosity removes the transferred kinetic energy. The large-scale spectrum depends on the forcing; the rising low-$k$ branch in this sketch is illustrative, rather than a universal power law.

The tube has neither homogeneous nor isotropic large-scale turbulence: gravity selects a vertical direction, its walls constrain the motion, and its density gradient varies with height. The width $W$ bounds the transverse size of the energy-containing eddies. At sufficiently high Reynolds number, much smaller eddies may become approximately locally isotropic and exhibit an inertial range, even though the full tube flow does not satisfy the homogeneous-isotropic assumptions.

To obtain the Reynolds stress, write $u_i=U_i+u_i^{\prime }$ and $p=P+p^{\prime }$, with $⟨u_i^{\prime }⟩=0$. Averaging the incompressible Navier-Stokes equations uses $⟨u_i{u}_j⟩=U_i{U}_j+⟨u_i^{\prime }{u}_j^{\prime }⟩$ and gives the Reynolds-averaged momentum equationHere the final term is the mean buoyancy force, with the hydrostatic reference removed. Thus averaging the quadratic momentum flux introduces the covariance $⟨u_i^{\prime }{u}_j^{\prime }⟩$; multiplied by ${\rho }_0$, it has the dimensions of a stress. Its divergence transfers momentum between the mean flow and the fluctuations.

An eddy viscosity model represents this Reynolds stress byThe isotropic term can be absorbed into the mean pressure, while the deviatoric term gives a diffusion of mean momentum with turbulent viscosity ${\nu }_T\sim u_{\ast }\ell$, where $u_{\ast }$ and $\ell$ are representative turbulent speed and mixing length. Unlike kinematic viscosity, this coefficient depends on the flow. The closure assumes a local relation between stress and mean strain; anisotropic, buoyancy-driven large eddies may require a more detailed Reynolds stress model.

Write $q={\partial }_z\bar{\rho }>0$ for the unstable, horizontally averaged density gradient. A turbulent displacement over a mixing length of order $W$ encounters a density contrast $\Delta \rho \sim Wq$. Balancing its buoyancy acceleration with the inertial acceleration givesA turbulent diffusivity for density is therefore $D_T=K{u}_{\ast }W=K{W}^2\sqrt{gq/\bar{\rho }}$, absorbing dimensionless factors into $K>0$. The upward density flux is $J_{\rho }=⟨w^{\prime }{\rho }^{\prime }⟩=-D_{T}q$: turbulent exchange carries denser water downwards and lighter water upwards. Defining $F_{\rho }=-J_{\rho }$ as the positive downward flux yields the buoyancy-driven turbulent density diffusion lawThis is a flux per unit horizontal area, with units of mass per area per time. The total density transport across a horizontal section is $W^2{F}_{\rho }$.

The assumptions are established high-Reynolds number turbulence, width-limited eddies, an instantaneous local balance between buoyancy and inertia, approximately constant $K$, negligible molecular diffusion, and zero mean vertical volume flux. The last condition follows from incompressibility and impermeable sidewalls in a closed tube. Density variations and the eddy statistics must also vary slowly enough for a local mixing-length closure to be meaningful.

Horizontal averaging of mass conservation, with upward $z$, gives ${\bar{\rho }}_t=-{\partial }_z{J}_{\rho }={\partial }_z{F}_{\rho }$. Differentiating once more supplies the requested evolution of the density gradient:Under the Boussinesq approximation, replace $\bar{\rho }$ in the inertial prefactor by the constant reference density ${\rho }_0$; retain the changing density gradient inside the flux.

Use the flat-surface Boussinesq approximation, with small thermal expansion, constant properties, negligible surface deformation, and constant surface tension. The top is impermeable and shear-free, while the solid walls are no-slip. Temperature-dependent surface tension would introduce a separate Marangoni stress, which is not included in this buoyancy-driven model. The stated time scale requires $\kappa$ to mean thermal diffusivity; if the material property given is a dimensional thermal conductivity $k$, use $\kappa =k/({\rho }_0{c}_p)$.

Take $z$ positive downwards and define the nondimensional streamfunction by $u={\psi }_z$, $w=-{\psi }_x$. Velocity is scaled by $\kappa /L$, streamfunction by $\kappa$, and temperature is measured from the cold-end surface value. Define the Prandtl number and horizontal-length Rayleigh number bywhere ${\alpha }_T$ is the thermal expansion coefficient. Put $J(\psi ,f)={\psi }_z{f}_x-{\psi }_x{f}_z$, and let ${\nabla }^2={\partial }_x^2+{\partial }_z^2$. Taking the curl of the Boussinesq equations eliminates pressure and gives the full nonlinear vorticity-temperature systemThe buoyancy sign follows from the downward coordinate: warmer fluid has body acceleration $-g{\alpha }_T\Delta TT{\mathbf{e}}_z$.

Choose the same streamfunction constant, zero, on the connected boundary. The mechanical and thermal boundary conditions areThe upper stress-free boundary condition reduces to ${\psi }_{zz}=0$, since $\psi =0$ there already gives ${\psi }_{xx}=0$. An arbitrary additive constant in $T$ changes only the reference hydrostatic pressure.

In steady infinite-Prandtl-number convection, momentum inertia disappears but thermal advection remains:For a shallow lake, with fixed $\mathrm{Ra}$ as $ϵ\to 0$, the leading interior temperature is $T_0=x$. This is an outer solution: it does not satisfy the insulated sidewall derivative right at $x=0,1$. The leading global conductive solution, which resolves those end regions, isIt is harmonic, has $T_c(x,0)=x$, and satisfies the insulated sides and bottom. Away from sidewall layers of horizontal width $O(ϵ)$ it reduces to $x$ with exponentially small end corrections. These layers also turn the flow and enforce the sidewall no-slip boundary condition.

In the shallow interior the dominant momentum derivatives are vertical, so ${\psi }_{0,zzzz}=-\mathrm{Ra}$. Integrating with the surface and bottom conditions gives the shallow horizontal convection with an insulated bottom solution, with $y=z/ϵ$:Surface flow is toward the colder end, with a deeper return flow. The depth-integrated volume flux is zero because ${\psi }_0(ϵ)={\psi }_0(0)=0$. The bulk expressions apply away from the turning regions, not unchanged on the vertical walls. A useful small-advection condition in the end layers is $\mathrm{Ra}{ϵ}^4\ll 1$, automatic for fixed Rayleigh number in the stated thin limit.

The heat flux in these units is $\mathbf{q}=(uT-T_x,wT-T_z)$. For the leading bulk temperature, the horizontal transport through a vertical section per unit out-of-plane width isThus leading advective horizontal transport vanishes: both layers have the same temperature at a given $x$, and their opposing volume fluxes cancel. The leading conductive transport is from the hot end toward the cold end. In dimensional form it is $-{\rho }_0{c}_p\kappa \Delta Tϵ=-k\Delta Tϵ$ per unit width. Heat enters through the warmer portion of the maintained upper surface and leaves through its colder portion, with the conductive end adjustments particularly important near the two insulated sidewalls. It cannot enter through the insulated sides or bottom.

For the advection correction in shallow horizontal convection, set $T=x+\Theta (z)$ in the bulk. Since $w_0=0$ and $T_x=1$, the steady heat equation gives ${\Theta }_{zz}=u_0={\psi }_{0,z}$. The conditions $\Theta (0)=0$ and ${\Theta }_z(ϵ)=0$ imply ${\Theta }_z={\psi }_0$ and henceFor example the bottom bulk temperature is $T(x,ϵ)=x-\mathrm{Ra}{ϵ}^5/320$. The correction makes the returning deeper fluid slightly colder. Its horizontal advective heat transport isThus the transport including this bulk correction is $\boxed{{\mathcal{H}}_x=-ϵ-19{\mathrm{Ra}}^2{ϵ}^9/1451520}$, both terms directed toward the colder end.

Vertical transport must distinguish local flux from net transport. In the leading interior $w_0=T_{0,z}=0$, whereas the temperature correction gives a local downward conductive flux $q_z=-{\Theta }_z=-{\psi }_0>0$. The end circulation gives an upward advective transport: with impermeable sidewalls and $T_0\simeq x$, integration by parts gives ${\int }_0^{1}w{T}_{0}dx=-{\int }_0^{1}x{\psi }_{x}dx={\int }_0^1\psi dx\simeq {\psi }_0(z)<0$. It compensates that downward conduction. More generally, steady heat conservation and insulated sides give $d[{\int }_0^1{q}_{z}dx]/dz=0$, and the insulated impermeable bottom sets the constant to zero. ThereforeThere is no net vertical heat transport through the lake, despite its local recirculating heat flux and its nonzero horizontal heat transport.

Use the steady axisymmetric top-hat plume model, with plume radius $b(z)$, uniform mean upward speed $w(z)>0$ and reduced gravity $g^{\prime }(z)=g[{\rho }_a(z)-{\rho }_p(z)]/{\rho }_0$. The Boussinesq approximation uses constant reference density ${\rho }_0$ in inertia and retains the small density difference only in buoyancy. Assume a quiescent hydrostatic ambient, a slender fully turbulent plume, negligible molecular transport and an entrainment velocity $\alpha w$, with constant positive entrainment coefficient $\alpha$. Define the ambient buoyancy frequency by $N^2=-(g/{\rho }_0)d{\rho }_a/dz$; positive $N^2$ is stable stratification.

The integrated volume, momentum and buoyancy balances of a Boussinesq top-hat plume in a stratified ambient areThe first follows from side entrainment, the second from the integrated buoyancy force, and the third accounts for the changing ambient density. Common factors of $\pi$ are cancelled. Equivalently, the requested radius, speed and reduced-gravity equations areThese equations stop describing an upward light plume once $w$ vanishes, and their light-plume interpretation requires $g^{\prime }>0$.

For nonzero $C$, seek plume similarity in power-law stratification in the form $b=\beta z$, $w=A{z}^n$ and $g^{\prime }=G{z}^m$. The first balance gives $\beta =2\alpha /(n+2)$. Momentum gives $m=2n-1$ and $G=2(n+1)A^2$. The buoyancy balance then requires $3n=2+n+p$, soChoosing $A>0$ gives the exact power-law family, including its normalization:These are the amplitudes and physical branches of a power-law stratified plume. No source flux was specified, so an arbitrary finite-source solution is not determined by $C,p$ alone; the displayed family is an exact similarity solution, rather than the general solution for arbitrary initial fluxes.

With the specific flux convention $Q=b^{2}w$, $M=b^2{w}^2$ and $B=b^{2}w{g}^{\prime }$, their values areThe physical volume flux is $\pi Q$, mass flux is ${\rho }_0\pi Q$, momentum flux is ${\rho }_0\pi M$, and kinematic buoyancy flux is $\pi B$. Thus the “specific mass flux” is the mass flux divided by reference density, with or without the explicitly stated factor $\pi$.

The parameter restrictions follow directly from positivity. Positive radius needs $p>-6$. Positive $g^{\prime }$ and real positive $w$ further require $p>-4$ and $-C/(3p+8)>0$, givingFor stable $C>0$, buoyancy flux decreases with height, while volume and momentum fluxes increase. However, $B\propto z^{4+3p/2}$ diverges at $z=0$ throughout this interval. Also the ambient density gradient and plume speed are singular there. Thus this stable similarity branch can describe a region away from a virtual source, but cannot represent a finite-buoyancy-flux point source extending all the way to $z=0$. For ordinary regular stable stratifications, including positive $p$, a finite-source rising plume instead loses buoyancy and eventually reaches neutral buoyancy and a finite rise limit; the pure-power ansatz does not supply such a solution.

For unstable $C<0$ and $p>-8/3$, the buoyancy flux grows with height and tends to zero at the origin. The plume draws its increasing buoyancy from an unstable ambient rather than a finite positive source buoyancy flux. Such an ambient must be maintained or treated locally, since it is itself convectively unstable. If $p\le -1$, its density profile is singular at the source, while for $p>-1$ its increasing density contrast eventually violates the Boussinesq approximation at sufficiently large height. These are formal local similarity solutions with corresponding physical limits. For $C<0$ and $-6<p<-4$, the algebra also gives a real positive speed but negative $g^{\prime }$: this is a heavy upward flow, not the specified light plume.

At nonzero $C$, the values $p=-6,-4,-8/3$ do not give finite coefficients. Respectively they would require zero entrained volume-flux growth, zero momentum-flux growth despite nonzero buoyancy, or constant buoyancy flux despite nonzero ambient stratification. They are not omitted regular endpoints.

When $C=0$, the ambient is unstratified regardless of the nominal value of $p$. The nontrivial positive-buoyancy point-source solution is the Boussinesq point-source plume:Here $B_0>0$ is the constant specific buoyancy flux with $\pi$ removed. Its $Q\propto z^{5/3}$ and $M\propto z^{4/3}$ vanish at the point source, while $B=B_0$. Its near-source singularity is the usual ideal point-source limit; a finite vent replaces that limit by finite source data or a virtual origin.

Use an inertial gravity-current box model in a horizontal prismatic V-shaped valley. Let its section have area $A(h)=\lambda h^2$ and width $2\lambda h$, where the dimensionless constant $\lambda >0$ is half-width divided by depth. Let $X(t)$ be front position and $h(t)$ the representative uniform current depth. Assume negligible entrainment and volume loss, hydrostatic cross-sections, large Reynolds number, and negligible bed friction and axial bed slope. These assumptions describe the bulk spreading after the short dam-break adjustment; they do not give a universal law for a steep or friction-dominated valley.

For water of fixed volume $W$, conservation gives $W=\lambda h^{2}X$. Close the front speed by the gravity-current front condition $\dot{X}=\mathrm{Fr}\sqrt{gh}$, where $\mathrm{Fr}$ is a constant order-one Froude number for the chosen inertial front model. ThereforeIntegration gives $X^{5/4}=X_0^{5/4}+5K_{w}t/4$. Once the initial reservoir length $X_0$ is negligible, the inertial dam-break current in a triangular valley satisfiesThe time exponent differs from a flat-bottomed channel because a V-shaped section stores volume proportional to $h^2$. A viscous V-shaped channel gravity current obeys a different spreading law, so retaining the inertial assumption is essential.

For the short-lived eruption, interpret $W$ as the bulk mixture volume released during the event. Let $\varphi (t)$ be the well-mixed particle volume fraction, initially ${\varphi }_0$, let ${\rho }_p$ and ${\rho }_a$ be particle and ambient carrier densities, and assume a dilute Boussinesq approximation for the mixture. The particles provide reduced gravityAssume thermal buoyancy, compressibility, erosion and continued injection are negligible, and that the carrier volume remains $W$. Turbulent mixing keeps the particles suspended uniformly while they settle vertically at a constant terminal settling velocity $w_s$. For small spherical particles of radius $a_p$ in carrier viscosity ${\mu }_a$, the Stokes settling velocity is $w_s=2({\rho }_p-{\rho }_a)g{a}_p^2/(9{\mu }_a)$ when the particle Reynolds number is small. If an eruption volume is quoted as solid particle volume rather than mixture volume, replace the $W$ used here by that volume divided by ${\varphi }_0$.

The particle deposition flux per horizontal projected bed area is $w_s\varphi$. The current's projected settling area is $2\lambda hX$, soThe factor two reflects the triangular section's mean depth $h/2$. Put $s=\sqrt{\varphi /{\varphi }_0}$ and $K=\mathrm{Fr}\sqrt{g_0^{\prime }}(W/\lambda {)}^{1/4}$. The front and concentration balances becomeTaking the initial release length negligible compared with the final runout, integration yields the constant-settling runout in a triangular valley:In particular $R\propto W^{3/7}(g_0^{\prime }{)}^{2/7}{w}_s^{-4/7}$ at fixed valley shape and front coefficient. With a finite initial length, replace $R^{7/4}$ by $X_0^{7/4}+7K\sqrt{W/\lambda }/(4w_s)$ in the corresponding integrated relation.

The whole time dependence is determined analytically by the quadratureEarly spreading has $X\simeq (5Kt/4{)}^{4/5}$, the water law with $g$ replaced by $g_0^{\prime }$. At late times the concentration loss slows the flow, and $R-X$ decays exponentially with rate $7K/(4R^{5/4})=w_s/h_R$, where $h_R=\sqrt{W/(\lambda R)}$. The final length $R$ is a limiting runout reached as $t\to \infty$, not a finite-time stopping event in this ideal model.

To calculate the deposit profile of a finite triangular-valley particle current, let $m_d(x)$ be final deposited particle mass per unit downstream length, and $M_0={\rho }_p{\varphi }_{0}W$ the initial particle mass. Deposition at a fixed section starts at arrival time $t_a(x)$, when the front first reaches $x$. With the instantaneous front coordinate denoted $\xi$, use $dt={\xi }^{1/4}d\xi /[K(1-(\xi /R{)}^{7/4})]$ to obtainWriting $\eta =x/R$, this evaluates toand $m_d=0$ beyond $R$. The profile decreases downstream, vanishes quadratically at the runout, and obeys ${\int }_0^R{m}_d(x)dx=M_0$. This integral check ensures that all initially suspended particle mass is accounted for.

For an areal rather than a linear deposit density, the changing valley width must also be specified. At lateral coordinate $y_{\perp }$, the current covers the bed only while $|y_{\perp }|\le \lambda h$. Set $Z=\min (R,\lambda W/y_{\perp }^2)$, with $Z=R$ on the valley axis. The mass per horizontal projected bed area is zero if $Z<x$, and otherwise isOn the axis this reduces to ${\rho }_p{\varphi }_0{h}_R[\frac{7}{5}(1-{\eta }^{5/4})-\frac{7}{12}(1-{\eta }^3)]$. Divide projected areal density by $\sqrt{1+{\lambda }^{-2}}$ to obtain mass per actual sloping bed area. A specified packed deposit density then converts this areal mass to local deposit thickness; the spreading model itself predicts deposited mass, not a packing fraction.

Assume a nonrotating, nonmagnetic equilibrium, inviscid adiabatic perturbations, and no dissipation. Let $g={\varphi }^{\prime }(R)>0$, so hydrostatic equilibrium is $p^{\prime }=-\rho g$. Use $\mathbf{\xi }{e}^{i\omega t}$ for the fluid displacement. To avoid a source ambiguity, write Eulerian fluid perturbations as $p_1,{\rho }_1,{\varphi }_1$, and use $\Delta$ for quantities following the displaced material element, as in Eulerian and Lagrangian fluid perturbations. They are related by $\Delta p=p_1+{\xi }_R{p}^{\prime }$ and similarly for mass density and potential.

The printed identity needs two corrections: its compression denominator must be $\gamma p$, not $\gamma \rho$, and the pressure and potential appearing in this form of the identity must be Eulerian fluid perturbations. The printed denominator has the wrong dimensions for an energy density. The claimed Lagrangian interpretation does not produce the displayed buoyancy-compression split; the following derivation specifies the consistent version.

Linearized mass conservation, the adiabatic relation, and the Poisson equation giveConsequently $p_1=-\gamma p\nabla \cdot \mathbf{\xi }-{\xi }_R{p}^{\prime }$. Linearizing the Euler equations for an inviscid fluid about the static star gives the linear adiabatic stellar oscillation equationsRequire regularity at the center, a vanishing Lagrangian surface pressure, and the usual matching to a decaying exterior gravitational potential. Because $p$ and $\rho$ vanish at the surface, the pressure boundary term vanishes for regular displacements; compactly supported trials also suffice below. The perturbed Newtonian gravitational potential solves Laplace's equation outside the star and must be included there when its field energy is computed.

Multiply the momentum equation by $-{\mathbf{\xi }}^{\ast }$ and integrate. The pressure term is $-{\int }_V{p}_1\nabla \cdot {\mathbf{\xi }}^{\ast }dV$. Let $d=\nabla \cdot \mathbf{\xi }$. The adiabatic relation and hydrostatic equilibrium implyThe pressure and equilibrium-gravity terms then combine pointwise aswhere $N^2=g[(1/\gamma )p^{\prime }/p-{\rho }^{\prime }/\rho ]$ is the stellar buoyancy frequency. The potentially troublesome cross terms cancel exactly.

For the perturbed self-gravity term, integration by parts givesThus the corrected stellar displacement energy identity isThe right side is a real quadratic form. With these conservative boundary conditions, the force operator is a self-adjoint operator. Its Rayleigh quotient characterizes squared mode frequencies, and a negative admissible trial quotient proves an unstable spectral component; it does not require the trial displacement itself to be an exact eigenfunction.

For a generic compressive disturbance of scale $\ell$ in a region with slowly varying equilibrium coefficients, estimate ${\rho }_1\sim \rho \xi /\ell$, $p_1\sim \gamma p\xi /\ell$ and ${\varphi }_1\sim 4\pi G{\rho }_1{\ell }^2$. The self-gravity and compression energy densities have orderswith adiabatic sound speed $c_s^2=\gamma p/\rho$. Under the plausible stellar estimate $c_s^2\sim G\bar{\rho }{a}^2$, global modes can have comparable gravity and pressure contributions, while their ratio for small-scale compressive modes is of order $(\rho /\bar{\rho })(\ell /a{)}^2$. Thus perturbations with $\ell \ll c_s/\sqrt{4\pi G\rho }$ admit the Cowling approximation. This estimate compares two nonzero compressive terms; it cannot alone justify neglecting gravity for a trial deliberately chosen to have $p_1=0$.

For an explicit convective trial, choose $c<d<e<b$ and a nonzero smooth bumpLet $Y$ be a normalized spherical harmonic of degree $\ell \ge 1$, let $L=\ell (\ell +1)$, and let ${\nabla }_{\Omega }$ be the unit-sphere gradient. The pressure-free localized stellar convective trial isIt is supported away from the center and surface, so it satisfies every needed mechanical boundary condition. Since ${\Delta }_{\Omega }Y=-LY$,Therefore $p_1=-\gamma p\nabla \cdot \mathbf{\xi }+\rho g{\xi }_R=0$ exactly. Its buoyancy energy is ${\int }_d^e\rho N^2{R}^2|f{|}^{2}dR<0$, whereas its kinetic norm ${\int }_d^e\rho R^2(|f{|}^2+L|h{|}^2)dR$ is positive. In the Cowling approximation this gives a negative Rayleigh quotient, and hence $N^2<0$ in an interior shell implies convective instability. Retaining perturbed self-gravity only adds a nonpositive term, so the same trial proves instability in the full conservative problem too.

There is also a quantitative small-scale check for this very trial. Its mass density amplitude is ${\rho }_1=(\rho N^2/g)fY$, independent of $\ell$. Write ${\rho }_1=s(R)Y$ and ${\varphi }_1=F(R)Y$. The radial Poisson equation givesBy the Cauchy-Schwarz inequality, $A\le 4\pi G(\int R^4|s{|}^{2}dR{)}^{1/2}(A/L{)}^{1/2}$, and so the angular-degree bound on perturbed stellar self-gravity isIt tends to zero with horizontal scale $R/\ell$, while the strictly negative buoyancy integral is fixed. This justifies the small-scale neglect of gravity even though the trial's compression term vanishes.

The magnetorotational instability extracts energy from differential rotation through magnetic tension. It can destabilize a disk whose specific angular momentum increases outward, even though such a disk satisfies Rayleigh's circulation criterion. The distinction is that an unmagnetized displaced parcel approximately conserves its own angular momentum, whereas a field line transfers angular momentum between parcels. If angular velocity decreases outward, an outward-displaced parcel lags its inner partner. Magnetic tension transfers angular momentum to the outer parcel and removes it from the inner one, making their radial separation increase. A sufficiently stiff magnetic connection suppresses the separation instead. The instability therefore requires a finite range of field-line-bending strengths.

Make this mechanism quantitative in the local shearing sheet of a cylindrical rotation profile $\Omega (r)$. Put $q=-d\log \Omega /d\log r$, and take $0<q<2$ so the unmagnetized flow is stable. Its radial epicyclic frequency is ${\kappa }^2=2(2-q){\Omega }^2$. Let the uniform field be vertical, $B_z{\mathbf{e}}_z$, with constant mass density, and use an axisymmetric horizontal disturbance proportional to $e^{st+ikz}$. Write $\mathbf{h}=\mathbf{b}/\sqrt{{\mu }_0\rho }$ and ${\omega }_A=k{B}_z/\sqrt{{\mu }_0\rho }$. The linearized ideal magnetohydrodynamics equations areThe first pair contains Coriolis acceleration and the epicyclic response; the last equation is the shear-winding of a radial field perturbation into an azimuthal one. Omitting that winding term would remove the driving mechanism.

Eliminating $h_r,h_{\varphi }$ gives $(s^2+{\omega }_A^2)u_r=2\Omega s{u}_{\varphi }$ and $s(s^2+{\omega }_A^2)u_{\varphi }+\Omega [(2-q)s^2-q{\omega }_A^2]u_r=0$. Their determinant is the ideal magnetorotational dispersion relationIts discriminant as a quadratic in $s^2$ is ${\kappa }^4+16{\Omega }^2{\omega }_A^2>0$. Since $2{\omega }_A^2+{\kappa }^2>0$, a positive $s^2$ exists precisely when its constant term is negative:Thus angular velocity, rather than specific angular momentum, must decrease outward. A Keplerian disk has $q=3/2$, ${\kappa }^2={\Omega }^2$ and is hydrodynamically stable, yet magnetic modes with $0<k^2{v}_A^2<3{\Omega }^2$ grow.

The growing root is $s^2=-{\omega }_A^2-{\kappa }^2/2+\frac{1}{2}\sqrt{{\kappa }^4+16{\Omega }^2{\omega }_A^2}$. Differentiation with respect to ${\omega }_A^2$ gives the maximum growth rate of ideal magnetorotational instability:For Keplerian shear these are $3\Omega /4$ and $15{\Omega }^2/16$. The ideal maximum is set by shear; reducing field strength reduces the fastest wavelength instead of reducing this maximum growth rate. This is why the exactly zero-field problem differs from a sequence of weak-field problems with progressively larger wave number. At fixed wave number, growth does tend to zero as the field vanishes. The wavelength must fit the actual disk, and diffusion eventually matters at sufficiently short scales.

The instability has the required angular-momentum-transfer sign. In a growing mode, $u_{\varphi }/u_r=(s^2+{\omega }_A^2)/(2\Omega s)>0$. At the fastest mode this ratio is one, while induction gives $h_{\varphi }=-h_r$ after allowing for their common Fourier phase. Thus the averaged Reynolds stress minus Maxwell stress tensor, $\rho ⟨u_r{u}_{\varphi }-h_r{h}_{\varphi }⟩$, is positive: angular momentum moves outward. Its product with $q\Omega$ is the energy supplied by the background shear. These local mechanisms and coefficient conventions are also checked in the primary MRI stability calculation.

More general fields require separating field-line bending from field gradients and curvature. At leading local order, for a stationary compatible field in an incompressible flow with axisymmetric wavevector $(k_r,0,k_z)$, the bending frequency isPressure elimination projects the radial force by $f=k_z^2/k^2$. The radial momentum equation becomes $s{u}_r-2f\Omega u_{\varphi }=i{\omega }_A{h}_r$; the azimuthal and induction equations retain the form above. Eliminating the amplitudes gives the oblique axisymmetric magnetorotational dispersion relationHere $k_z\ne 0$. The unstable range is $0<{\omega }_A^2<-fd{\Omega }^2/d\log r$. The poloidal magnetic field therefore determines whether an axisymmetric perturbation bends a field line. A locally uniform toroidal magnetic field does not enter $\mathbf{k}\cdot \mathbf{B}$ for axisymmetry, and can contribute magnetic pressure that is projected out of this incompressible leading-order problem. Thus a toroidal component does not suppress the weak-field MRI merely by being present. Conversely a purely toroidal field has no ordinary axisymmetric local bending MRI in this approximation. This does not imply that every globally curved toroidal equilibrium is stable.

A steady general poloidal magnetic field must also be compatible with the equilibrium rotation. The azimuthal ideal magnetohydrodynamic induction equation contains ${\partial }_t{B}_{\varphi }=r{\mathbf{B}}_p\cdot \nabla \Omega$. A stationary axisymmetric equilibrium with no meridional flow consequently requires ${\mathbf{B}}_p\cdot \nabla \Omega =0$, Ferraro's law of isorotation. A radial field crossing cylindrical differential-rotation surfaces would otherwise wind secularly; it cannot simply be substituted into a stationary normal-mode calculation. Boundary conditions and the along-field eigenvalue problem replace a freely chosen local $k$ in nonuniform or bounded equilibria. Strong enough field-line tension can then stabilize all admissible wavelengths.

An explicit axisymmetric effect of nonuniform field strength and curvature is obtained with a purely toroidal field $B_{\varphi }(r){\mathbf{e}}_{\varphi }$, constant mass density and cylindrical rotation. A radial fluid displacement changes the field byThis follows directly from $\mathbf{b}=\nabla \times (\mathbf{\xi }\times \mathbf{B})$ and incompressibility. Its radial curvature force, after absorbing the magnetic-pressure perturbation, is $-2B_{\varphi }{b}_{\varphi }/({\mu }_0\rho r)$. The axisymmetric azimuthal momentum equation is $s{u}_{\varphi }+(2\Omega +r{\Omega }^{\prime })u_r=0$. Combining it with the radial force and then projecting out pressure givesThis is the local Michael criterion for axisymmetric toroidal-field interchange. The new term can either stabilize or destabilize: for $B_{\varphi }\propto r^p$, $\mathcal{D}={\kappa }^2+2(1-p)v_{A\varphi }^2/r^2$. A current-free $B_{\varphi }\propto r^{-1}$ is stabilizing; $B_{\varphi }\propto r$ leaves this axisymmetric interchange coefficient unchanged; a sufficiently steep increase with $p>1$ can destabilize even a centrifugally stable rotation law. That last instability can persist without differential rotation, so its energy source is the magnetic-current equilibrium rather than the ordinary MRI. Nonaxisymmetric current-driven modes are outside the requested essay scope.