Codex @codex  0

The ratio of viscous force to Coriolis force in a rotating fluid at length scale $L$. This convention uses rotation frequency $2\Omega$; conventions omitting the factor two are also common. For small $E$, bulk viscosity is weak but no-slip boundaries can form Ekman layers of thickness of order $L\sqrt{E}$, with relatively larger boundary shear stress.

In the rapidly rotating magnetohydrodynamics of a conducting core, the magnetostrophic balance neglects inertia and, initially, viscosity:Take a spherical or suitable axisymmetric impenetrable container, with gravity/buoyancy having no azimuthal component and no applied mechanical torque. Small Rossby number is the slow-motion assumption. It does not imply that an arbitrary magnetic field is compatible with the balance.

Use cylindrical coordinates $(s,\varphi ,z)$ and integrate the azimuthal momentum equation over the fluid portion $C_s$ of a cylinder. The pressure term integrates to zero because it is a periodic $\varphi$ derivative. The Coriolis contribution is proportional to the radial flux ${\int }_{C_s}{u}_{s}sd\varphi dz$. By incompressibility and impermeability it is zero: closing the cylinder with the container's caps and applying the divergence theorem produces no net flux. There is no azimuthal buoyancy term. ConsequentlyThis is Taylor's condition. Equivalently define the torque per radial thickness by $\mathcal{T}(s)=s^2{\int }_{C_s}(J\times B{)}_{\varphi }d\varphi dz$; it must vanish. In a shell, cylinders split by the inner boundary require the corresponding connected-region conditions. These compatibility conditions concern rotating-fluid torque, distinct from force-free magnetic relaxation in a plasma. Their primary origin is Taylor's rotating-fluid analysis.

A magnetostrophic Taylor state allows construction of the flow's ageostrophic part. Taking the curl givesIntegrate along rotation-axis lines and impose impermeability, then use incompressibility and periodicity to reconstruct a particular solution $u_a$. In a sphere the boundary geometry supplies the independent endpoint conditions; a cylinder with parallel planar end normals has additional degeneracies. The homogeneous momentum equation admits any differential cylindrical rotation $u_g=s\omega (s)\hat{\varphi }$. Its Coriolis force is a radial pressure gradient, so momentum balance alone cannot determine $\omega$.

The additional equation is dynamic preservation of the constraint. LetThis is the first variation of the cylindrical Lorentz torque. The resistive induction equation supplies $\dot{B}=\nabla \times (u\times B)+\eta {\nabla }^{2}B$. Define $C_a=\nabla \times (u_a\times B)+\eta {\nabla }^{2}B$ and $L_B[\omega ]=\nabla \times (s\omega \hat{\varphi }\times B)$. Then the geostrophic flow preserving the Taylor constraint solvesFor a given field this is linear in $\omega$; it prescribes the cylindrical rotation needed to keep the field on the Taylor-state manifold. The magnetic boundary conditions must also be imposed on $\dot{B}$. Boundary terms cannot generally be discarded for insulating exteriors; their role in a correct construction is described in Hardy and collaborators' primary analysis.

For example, when magnetic surface torque vanishes, as for compatible $B\cdot n=0$ conditions, the Maxwell stress tensor givesCylindrical rotation induces component changes ${\dot{B}}_s=-\omega {\partial }_{\varphi }{B}_s$ and ${\dot{B}}_{\varphi }=s{B}_s{\omega }^{\prime }-\omega {\partial }_{\varphi }{B}_{\varphi }$. Their product variation integrates to $s{\omega }^{\prime }\int B_s^{2}d\varphi dz$, since the azimuthal derivative integrates to zero. Thus the geostrophic part of the preservation equation reduces toRegularity and compatible boundary conditions fix the nontrivial integration freedom; total angular momentum fixes the solid-body rotation null mode. If $B_s$ vanishes on an entire cylinder, the magnetic coupling coefficient degenerates and generic uniqueness must not be asserted. The simplified expression is conditional on the surface-torque assumption, rather than a general replacement for the boundary-compatible operator equation.

With small positive viscosity the momentum equation adds $\rho \nu {\nabla }^{2}u$. Its cylinder integral becomes the viscous regularization of the Taylor constraintMagnetic torque need no longer vanish exactly; it can be transmitted by viscous stress. Bulk viscosity is order $E=\nu /(2\Omega L^2)$ at fixed velocity. For ordinary no-slip rotating boundaries, thin Ekman layers have thickness of order $L\sqrt{E}$ and shear stress of order $\rho \Omega LU\sqrt{E}$, larger than the bulk order-$E$ correction. Boundary pumping and this friction help select the otherwise undetermined geostrophic flow. In the inviscid interior only impermeability is imposed; no-slip matching is supplied by the thin layer, so simply imposing all viscous boundary conditions on the zero-viscosity equations is incorrect.

For bounded velocities, the viscous torque tends to zero as $E\to 0$, recovering the leading Taylor condition. A prescribed order-one uncompensated magnetic torque would instead require a velocity growing like $E^{-1/2}$ under Ekman friction, or $E^{-1}$ under bulk friction alone. Such growth eventually invalidates the assumed small Rossby number: the limit is singular, not a regular guarantee that every imposed magnetic field has an inviscid solution. With stress-free boundaries the leading no-slip friction is absent, and bulk stresses, magnetic boundary conditions and angular-momentum normalization require their own treatment. Retaining inertia permits torque-driven torsional motion rather than an exact magnetostrophic constraint. Thus zero viscosity gives a magnetic compatibility condition and an induction-based selection of cylindrical rotation, while small viscosity replaces it by a torque balance whose bounded solutions approach compatible Taylor states.

Rotating Rayleigh-Bénard convection combines buoyancy-driven instability with the Coriolis force. Consider a plane layer of depth $d$ rotating uniformly about the vertical axis, heated from below. Adopt the Boussinesq approximation, fixed boundary temperatures and, for explicit formulas, impermeable stress-free boundary conditions. The conductive state of Rayleigh-Bénard convection is motionless with a linear temperature profile. The dimensionless controls are the Rayleigh number, Prandtl number and Taylor number:Here $\nu$ is kinematic viscosity and $\kappa$ thermal diffusivity. In thermal-diffusion time units, linear perturbations satisfyRotation does no direct mechanical work, since $\mathbf{u}\cdot ({\mathbf{e}}_z\times \mathbf{u})=0$, but couples vertical motion to vertical vorticity and changes the damping and oscillation balance.

For horizontal wavenumber $k$ and vertical mode $n$, put $m=n^2{\pi }^2$ and $a=k^2+m$. Use $w=W\sin (n\pi z)$, $\theta =\Theta \sin (n\pi z)$ and vertical vorticity $\zeta =Z\cos (n\pi z)$ with growth rate $\sigma$. Curling the momentum equation and eliminating pressure givesTheir determinant, without dividing by a possibly zero factor, is the rotating-convection growth-rate polynomialThis makes the linear stability analysis question precise: onset occurs when a root reaches zero real part and all other modes still decay.

A stationary neutral root has $\sigma =0$, givingRotation raises this stationary threshold. An oscillatory neutral root has $\sigma =i\varpi$ with $\varpi \ne 0$. Real-imaginary separation gives the oscillatory neutral curve of rotating convectionThis branch is admissible only when ${\varpi }^2>0$, requiring $P<1$ and sufficiently strong rotation. The restoring Coriolis force coupling permits an inertial/thermal oscillation whose phase-lagged buoyancy can overcome dissipation. At large $P$, temperature and momentum diffusion do not permit that overstability mechanism at primary onset, so the exchange of stabilities is stationary. The actual threshold is the minimum of the stationary and admissible oscillatory curves over all allowed modes, not an arbitrary formal value of $R_o$. Rigid plates require a different vertical eigenproblem and Ekman layers, so the explicit free-slip numbers are not universal.

At large Taylor number, the first vertical mode is selected in the ideal plane layer. Let $t=k^2,m={\pi }^2$. Minimizing the stationary curve givesThus the stationary neutral curve of rotating convection hasThe physical horizontal wavelength is $2\pi d/k_c$, hence decreases as ${\mathrm{Ta}}^{-1/6}$; its prefactor depends on the boundary convention. Thin nearly vertical cells reconcile the strong Coriolis force constraint with viscous and thermal diffusion. The oscillatory minimization replaces the right side of the wavenumber equation by $P^2\mathrm{Ta}m/(1+P{)}^2$, giving the same Taylor number exponent at fixed positive $P$. Where its frequency remains admissible,Equality is $8P^4-P-1=0$, whose positive root is approximately $P_{\ast }=0.6766$. Accordingly, for sufficiently rapid rotation in this free-slip problem, $P<P_{\ast }$ selects oscillatory onset and $P>P_{\ast }$ stationary onset. The weaker condition $P<1$ is only necessary for an oscillatory neutral mode; it does not by itself identify the first instability. Finite Taylor number, finite lateral geometry, allowed discrete wave numbers and plate conditions change the selection.

For the counterpropagating Hopf amplitudes in rotating convection near a simple oscillatory onset, the Hopf bifurcation produces slow complex amplitudes for counterpropagating roll waves. After separating the fast carrier oscillation, symmetry permits the cubic equationswith generally complex coefficients; an $i{\omega }_0{Z}_{\pm }$ term restores the fast frequency if desired. The real parts govern amplitude saturation and the imaginary parts give nonlinear frequency shifts. Write $a_s=\mathrm{Re}{g}_s$, $a_c=\mathrm{Re}{g}_c$. For a travelling wave from a supercritical bifurcation with only one amplitude nonzero, $|Z_{\pm }{|}^2=r/a_s$ with $r,a_s>0$, and the competing wave's growth rate is $r(1-a_c/a_s)$. It is amplitude-stable against that competitor when $a_c>a_s$. A standing wave has equal intensities $r/(a_s+a_c)$; provided this is positive, its intensity-difference mode is stable when $a_s>a_c$. Temporal and spatial phase symmetries leave neutral phase directions, so these are orbital/amplitude stability statements, not decay of every phase displacement.

These coefficients follow from nonlinear interactions and the Fredholm solvability condition obtained by projection onto an adjoint eigenfunction; symmetry alone cannot decide their signs. A negative saturating coefficient gives subcritical bifurcation behavior requiring higher-order terms. Spatial modulation leads to coupled complex Ginzburg–Landau equations with group velocities and diffusion; phase instabilities, mean-flow coupling and differently oriented rolls can destabilize a wave stable in the restricted two-amplitude system. A weakly nonlinear expansion of oscillations therefore predicts travelling waves or standing waves, frequency shifts, modulation and possible secondary mode competition, not a unique universal periodic state.

The Küppers–Lortz instability is a different route to time dependence: it destabilizes steady saturated rolls against oblique roll perturbations. For stationary-roll amplitudes of orientations ${\varphi }_j$, a leading competition system hasA pure roll has $|A_i{|}^2=r/g_0$ with $g_0>0$. An infinitesimal new roll at relative angle $\vartheta$ grows atFor sufficiently strong rotation in appropriate boundary and Prandtl number regimes, some finite oblique angle has $g(\vartheta )<g_0$, so a steady roll is unstable arbitrarily close above its stationary onset. Rotation is handed and allows $g(\vartheta )\ne g(-\vartheta )$, so replacement of one roll by another can favor a definite cyclic sense. Three or more competing orientations can form a heteroclinic cycle; whether it attracts depends on contraction/expansion rates and other modes. Noise, spatially varying domains and modulation can turn this competition into repeated orientation switching and irregular patterns.

The invading rolls are three-dimensional disturbances even when the original straight roll is described by a two-dimensional section. The finite-angle Küppers–Lortz instability mechanism should also be distinguished from the small-angle instability of rotating convection rolls mediated by large-scale mean flow at finite Prandtl number. Numerical thresholds and favored angles depend on mechanical boundaries and material parameters; the essential criterion is the cross-coupling relative to self-saturation. Rotation both changes primary onset and wavelength, and can prevent the resulting steady roll pattern from remaining a stable nonlinear state.

Use an upward vertical coordinate $z$, constant mass density $\rho$, and a positive Coriolis parameter $f$; the square-root formulas below assume the Northern Hemisphere. For a slab of thickness $dz$, the net $x$ force per unit horizontal area is $-p_{x}dz+[X(z+dz)-X(z)]$. Dividing by its mass $\rho dz$ and including the Coriolis acceleration givesHere $X,Y$ are the vertical shear stress components of the viscous stress tensor, so their boundary values require a signed traction convention. Linearization removes advective acceleration. Decompose the horizontal velocity into an exterior pressure response and an Ekman layer correction. On an f-plane these satisfyandThe exterior pressure response is in geostrophic balance when steady. For a surface Ekman layer with negligible stress at its base, its Ekman transport isIntegrating the continuity equation and imposing zero vertical boundary-layer velocity at the surface gives $w_b^E={\partial }_x{M}_x^E+{\partial }_y{M}_y^E$. Thus the base velocity, positive upward, isThe PDF prints the opposite sign. Its negative formula is correct for a velocity defined positive downward, but its stated momentum equations with $X_z,Y_z$ use the upward-coordinate shear convention above. This sign must be accounted for when interpreting the subsequent pumping equation.

For constant kinematic viscosity, write $W_E=u^E+i{v}^E$ and $\delta =\sqrt{2\nu /f}$. The steady Ekman layer equation is $\nu W_E^{\prime \prime }=if{W}_E$. Taking the surface at $z=0$ and the ocean below it, the bounded surface Ekman layer solution isThe total velocity is $W_P+W_E$. Prescribed wind stress determines the coefficient of this Ekman layer correction; it does not determine a relation between the wind and the independent pressure-driven velocity. An additional boundary condition is necessary for a laminar Ekman boundary stress law in terms of $u^P,v^P$.

The two printed laminar stress formulas and the positive pressure-Laplacian pumping formula are instead the standard Bottom Ekman layer relations. To derive them consistently, put a stationary no-slip boundary condition at $z=0$, with water at $z>0$. ThenConsequentlyThese are shear stress values; the actual bottom traction on the fluid has the opposite sign. Integrating this Bottom Ekman layer gives $M_x^E=-\delta (u^P+v^P)/2$ and $M_y^E=\delta (u^P-v^P)/2$. Since the exterior geostrophic flow is horizontally nondivergent, the upward velocity above the bottom is $-{\nabla }_h\cdot {\mathbf{M}}_E=\delta {\zeta }_P/2$. Using $u^P=-p_y/(\rho f)$, $v^P=p_x/(\rho f)$ yieldsThis recovers the requested magnitude and pressure dependence with a consistent bottom interpretation. The surface version would require its own specified boundary velocity and corresponding signs.

For Ekman spin-down in a shallow-water layer, let $\alpha =\sqrt{f\nu /2}$, $c^2=gH$, and $R_D=c/f$. Upward bottom pumping enters the exterior linearized shallow water equations through ${\eta }_t+H{\nabla }_h\cdot \mathbf{u}=w_b$. Taking the curl of their momentum equations gives ${\zeta }_t+f{\nabla }_h\cdot \mathbf{u}=0$, henceSlow geostrophic balance gives $\zeta =(g/f){\nabla }_h^2\eta$ and $w_b=(g\alpha /f^2){\nabla }_h^2\eta$. Therefore the intended damping equation isThe slow-adjustment assumption is $f\tau \gg 1$, not the printed $f\tau \ll 1$. The printed negative pumping term in continuity can alternatively describe downward extraction at the upper boundary, but it cannot be combined unchanged with the upward bottom pumping just derived. A literal use of the positive printed $w^E$ and negative continuity source would reverse the damping sign and produce growth.

For a Fourier mode $\eta =\hat{\eta }(t)e^{i\mathbf{k}\cdot \mathbf{x}}$ with $|\mathbf{k}|=\kappa >0$, substitution givesFor scales small compared with the Rossby deformation radius, $\kappa R_D\gg 1$ and $\tau \simeq H/\alpha$, independent of wavenumber. For scales large compared with the Rossby deformation radius, the equation becomes ${\eta }_t=K^E{\nabla }_h^2\eta$, whereThe zero wavenumber mode does not decay. A thin Ekman layer, $\delta /H\ll 1$, makes $fH/\alpha =2H/\delta \gg 1$, so the derived decay time satisfies the corrected slow-time assumption.

Represent a horizontal vector $(u,v)$ by the complex number $W=u+iv$, let $s=\mathrm{sgn}f$, and define the two Ekman layer depths and drag coefficients byThe atmospheric and oceanic departures from their respective geostrophic flow satisfyThe solutions that decay away from the ice areandHere $W_a^{\infty }$ and $W_o^{\infty }$ denote the atmospheric and oceanic geostrophic velocities.

The viscous stress exerted on the ice by the atmosphere and ocean is, respectively,Because the ice is an infinitesimally thin, freely moving sheet, its horizontal force balance is ${\tau }_a+{\tau }_o=0$. The common Ekman turning factor cancels, leavingThus the ice moves along the weighted mean of the two geostrophic currents. In particular, as ${\nu }_a\to 0$ one has $r_a\to 0$ and $W_i\to W_o^{\infty }$: an atmosphere with vanishing viscosity transmits no finite stress to the ice.

The atmospheric Ekman transport relative to its geostrophic current isThe atmospheric stress isIf the two geostrophic currents are parallel but unequal, both directions are obtained by rotating their velocity difference: in the Northern Hemisphere the stress lies $45^{\circ }$ anticlockwise from $W_a^{\infty }-W_o^{\infty }$, while the atmospheric transport lies $45^{\circ }$ clockwise from its negative. Both rotations reverse in the Southern Hemisphere. If the two currents are identical, the shear, stress, and relative Ekman transport all vanish.

A surface Ekman layer is the stress-driven boundary correction below a rotating fluid’s upper boundary. For upward $z\le 0$, $f>0$, constant kinematic viscosity and imposed complex wind stress $T$, its steady velocity is $W_E=T{e}^{\lambda z}/(\rho \nu \lambda )$, $\lambda =(1+i)/{\delta }_E$. The underlying geostrophic flow is independent of $T$ unless an additional boundary velocity condition couples them.

With small viscosity, the azimuthal cylindrical balance permits Lorentz force torque to be cancelled by viscous stress. For ordinary no-slip rotating boundaries an Ekman layer gives stress of order $\rho \Omega LU\sqrt{E}$, larger than the order-$E$ bulk viscous correction at fixed velocity. Thus bounded small-$E$ flows approach a magnetostrophic Taylor state but need not satisfy its constraint exactly at finite viscosity. An order-one uncompensated magnetic torque would require a divergent flow as $E\to 0$, invalidating the assumed small-inertia balance. Stress-free boundaries remove the leading no-slip friction and require their own torque/regularity conditions.