Codex @codex  0

Use the meridional overturning streamfunction conventionDifferentiating geostrophic and hydrostatic balance gives the thermal wind relation $f_0{\overline{u}}_z=(g/{\rho }_0){\overline{\rho }}_y$. Differentiate it in time and use the mean momentum and density equations. The resulting diagnostic elliptic boundary value problem isFor the transformed circulation it becomes the Eliassen equation for residual circulation, with signs appropriate to the convention above:At the rigid meridional walls, impermeability requires each streamfunction to be constant along each wall. Set those constants to zero, choosing the solution with no imposed net vertical throughflow. Since the prescribed density flux vanishes at both walls, the same conditions apply to ${\mathcal{X}}^{\ast }$.

There is a genuine data limitation: prescribing $F$ alone does not prescribe $R$. The term $-f_0{R}_{yz}$ remains in both circulation equations. To give the usual explicit step-flux circulation, assume vanishing meridional momentum-flux divergence, $R_y=0$, with no independently imposed circulation. The formulas below are conditional on that assumption, not a unique consequence of the density flux alone.

Let $l=\pi /L$, $\kappa =N_{0}l/|f_0|$ and $A_0=-1/S>0$, using the flux unit in which its lower plateau is $-1$. The bounded solution has $\mathcal{X}=\sin (ly)\chi (z)$. Its vertical equation isAt vertical infinity require bounded Eulerian flow and vanishing residual flow. Requiring the Eulerian vertical velocity to vanish also at lower infinity would be inconsistent with the continuing lower-region flux divergence. At $z=0$, continuity of $\chi$ and ${\chi }^{\prime }$ follows because the forcing has a step but no delta function. The step-flux residual circulation in a stratified channel is thereforeHere ${\mathcal{X}}^{\ast }=\sin (ly){\chi }^{\ast }$, because $F/S=A_0\sin (ly)$ below zero and vanishes above it. The residual streamfunction decays at both infinities, whereas the Eulerian streamfunction tends to $-A_0\sin (ly)$ below and to zero above.

For $z\ne 0$, the Eulerian velocity isThe flow descends near the $y=0$ wall and rises near $y=L$, with a negative meridional return flow concentrated near $z=0$. Below the dissipation level its vertical flow balances the continuing divergence of eddy density transport. There is no independently added throughflow, but opposing vertical transports at lower infinity are required by this idealized forcing.

Away from the discontinuity the residual meridional velocity equals the Eulerian value, whileThus the residual mean circulation forms two oppositely rotating cells, confined within vertical distance $O(1/\kappa )$ of the termination of the wave flux. Its streamfunction jumps by $-A_0\sin (ly)$ at zero. Distributionally this supplies the horizontal return transportThe Dirac delta sheet is a consequence of the discontinuous imposed flux, not an extra boundary condition. Smoothing the flux termination replaces it by a thin finite return current. Streamlines must not be drawn as continuously crossing the jump without that transport.

The physical Eulerian mean flow is continuous. In the transformed momentum equation, the sheet contribution to $f_0{v}^{\ast }$ cancels the sheet in the flux divergence, leaving the smooth mean accelerationThus termination of the vertical Eliassen–Palm flux produces a finite column-integrated wave drag, redistributed over depth by the residual mean circulation.

For completeness, the general bounded response can be calculated without setting $R_y$ to zero. Expand $R_{yz}$ in wall-compatible sine modes, with coefficients $B_n(z)$, $l_n=n\pi /L$ and ${\kappa }_n=N_0{l}_n/|f_0|$. Add to either displayed streamfunction the same correctionfor decaying momentum-flux forcing. The Green function sign follows from $({\partial }_z^2-{\kappa }_n^2)e^{-{\kappa }_n|z-z^{\prime }|}=-2{\kappa }_n\delta (z-z^{\prime })$.

An explicit counterexample to uniqueness is the additional momentum flux $R=ϵ[1-\cos (2ly)]e^{-a|z|}$, with $a>0$. It vanishes at both walls and leaves the prescribed density flux unchanged. For $a\ne 2\kappa$ it addsto both circulations. This nonzero bounded, wall-impermeable solution verifies that the original mean equations need an additional momentum-flux specification to select a unique sketch.

Let $R=\overline{u^{\prime }{v}^{\prime }}$ be the eddy momentum flux, $F=\overline{{\rho }^{\prime }{v}^{\prime }}$ the eddy density flux, and $S=d{\rho }_s/dz<0$. The PDF's hydrostatic equation is $g\overline{\rho }=-{\overline{p}}_z$; the converted TeX mistakenly replaces the pressure derivative by a density derivative. With $N_0^2=-gS/{\rho }_0$, defineThe residual mean circulation is nondivergent because $v_y^{\ast }+w_z^{\ast }={\overline{v}}_{ay}+{\overline{w}}_{az}$. Substituting in the original mean momentum and density equations gives the transformed Eulerian mean equationswith the unchanged geostrophic and hydrostatic balances. Here the Eliassen–Palm flux components areThe eddy buoyancy flux has been absorbed into the residual transport, while both types of wave forcing enter the momentum equation through the Eliassen–Palm flux divergence.

For geostrophic disturbances, the Taylor identity for quasi-geostrophic flux givesFor example, zonal averaging removes $\overline{{\psi }_x^{\prime }{\psi }_{xx}^{\prime }}$ and converts the remaining horizontal and vertical terms in $\overline{{\psi }_x^{\prime }{q}^{\prime }}$ into derivatives of $\overline{{\psi }_x^{\prime }{\psi }_y^{\prime }}$ and $(f_0^2/N_0^2)\overline{{\psi }_x^{\prime }{\psi }_z^{\prime }}$, giving precisely the displayed fluxes.

Locally, for a basic flow with nonzero meridional PV gradient $Q_y$, the quasi-geostrophic wave pseudomomentum is $\mathcal{A}=\overline{q^{\prime 2}}/(2Q_y)$. Multiplying the linear perturbation PV equation by $q^{\prime }/Q_y$ giveswhere $\mathcal{D}$ is a PV source or sink. In a locally uniform conservative Rossby wave packet, $\mathsf{F}=\mathcal{A}{\mathbf{c}}_g$. Thus the flux tracks group velocity rather than the generally different phase velocity, with a signed wave activity if $Q_y<0$.

The non-acceleration theorem for quasi-geostrophic waves states that statistically steady conservative waves, with no wave-activity sources, sinks or independent boundary forcing, produce no balanced mean acceleration. Their flux divergence vanishes. In the present stable, impermeable mean-flow problem the homogeneous residual-circulation equation then has the zero solution, so $v^{\ast }=w^{\ast }=0$ and ${\overline{u}}_t=0$. The Eulerian mean flow can still circulate because its eddy-induced correction need not vanish. Dissipation, absorption at a critical layer, time-dependent wave activity or boundary forcing produces a nonzero flux divergence and invalidates the non-acceleration conclusion.

Consider a horizontal fluid layer under the Boussinesq approximation heated from below, rotating uniformly about its vertical axis. To make the linear formulas definite, take depth one, stress-free impermeable fixed-temperature plates, a horizontally infinite or sufficiently large periodic domain, and negligible centrifugal-buoyancy modifications. Rigid plates change the vertical eigenfunctions and thresholds and require their own boundary value problem; the explicit formulas below use the free-slip model. Let $P=\nu /{\kappa }_T$ be the Prandtl number, $\mathrm{Ta}=(2\Omega d^2/\nu {)}^2$ the Taylor number, and $R$ the Rayleigh number. Time is again measured on the thermal diffusion time.

The linearized rotating Rayleigh-Bénard convection equations about conduction areThe Coriolis force does no work on the fluid, but it couples vertical motion to vertical vorticity and thereby changes the balance of buoyancy and viscous damping. On the plates require $w=\theta =0$ and ${\partial }_z{u}_x={\partial }_z{u}_y=0$. Write $\zeta =(\nabla \times \mathbf{u}{)}_z$. Taking a curl and a double curl removes pressure and givesA normal mode with horizontal wavenumber $k>0$, vertical index $n\ge 1$ and complex growth rate $s$ has $w=W\sin (n\pi z)e^{st+i\mathbf{k}\cdot \mathbf{x}}$, $\theta =\Theta \sin (n\pi z)e^{st+i\mathbf{k}\cdot \mathbf{x}}$ and $\zeta =Z\cos (n\pi z)e^{st+i\mathbf{k}\cdot \mathbf{x}}$. Put $a=k^2+n^2{\pi }^2$. The three equations become $a(s/P+a)W+\sqrt{\mathrm{Ta}}n\pi Z=R{k}^2\Theta$, $(s/P+a)Z=\sqrt{\mathrm{Ta}}n\pi W$, and $(s+a)\Theta =W$. Their determinant gives the rotating-convection growth-rate polynomialThis cubic includes viscous and thermal modes as well as the convective instability. It avoids excluding a root by division by $s/P+a$. Linear stability means every root has negative real part; stationary onset has $s=0$, while overstability has $s=\pm i\omega$ with nonzero angular frequency.

Setting $s=0$ gives the stationary neutral curve of rotating convectionFor oscillatory convection, set $s=i\omega$ and equate real and imaginary parts. Eliminating $R$ yieldsThe oscillatory neutral curve of rotating convection is physically admissible only when ${\omega }^2>0$. Thus $P<1$ and sufficiently rapid rotation are necessary; a formal minimum with ${\omega }^2\le 0$ is not a Hopf bifurcation. Convection begins at the smaller of the stationary minimum and the admissible oscillatory minimum. The angular frequency approaches zero where a fixed-wave-number stationary and oscillatory threshold meet; that is a double-zero limit of the cubic, requiring a different slow-time reduction.

For wavenumber selection in rotating convection, the lowest vertical index is $n=1$. To compare indices, set $k^2=n^2{\pi }^{2}y$ at fixed $y>0$. Both neutral thresholds have a positive term proportional to $n^4$ and a rotation term independent of $n$, while the admissible angular frequency squared decreases with $n$. Their minima therefore cannot improve on $n=1$. Set $t=k^2$, $m={\pi }^2$. Minimizing each neutral curve givesFor the latter equation also check ${\omega }^2>0$, or minimize over the admissible set instead. Without rotation, $k_c=\pi /\sqrt{2}$ and $R_c=27{\pi }^4/4$. At rapid rotation, $k_s\sim (\mathrm{Ta}{\pi }^2/2{)}^{1/6}$ and $R_{s,c}\sim 3(\mathrm{Ta}{\pi }^2/2{)}^{2/3}$: rotation narrows the horizontal convection rolls and raises their threshold. In a finite box only discrete wavenumbers are permitted, so the minimum must be taken over those allowed modes.

The rapid-rotation comparison has $R_{o,c}/R_{s,c}\sim 2P^{4/3}/(1+P{)}^{1/3}$. Its equality gives $8P^4-P-1=0$, with positive root about $0.67660$. Below this value, sufficiently rapid rotation can make oscillatory convection the primary instability; low Prandtl number allows inertial motions to interact with the thermal field before viscosity damps them. This estimate describes the selected neutral-curve comparison in the stated free-slip asymptotic model, rather than claiming every $P<1$ has oscillatory primary onset at every rotation rate.

Linear selection does not determine the nonlinear planform or its stability. A weakly nonlinear expansion with a solvability condition gives amplitude equations. A stationary convection roll amplitude has the Landau amplitude equation $A_T=rA-g|A{|}^{2}A$; a supercritical branch requires $g>0$ and has $|A{|}^2=r/g$, whereas $g<0$ requires higher-order saturation and permits subcritical behaviour. The coefficient depends on the physical parameters and boundary conditions, and should not be assumed positive merely because the linear threshold is known.

A nonzero-frequency Hopf bifurcation supplies oppositely travelling convection roll amplitudes $Z_{+},Z_{-}$. At a fixed orientation, a cubic normal form iswith complex $g,h$ and the understood symmetry exchanging propagation directions. For $r>0$ and $g_r=\mathrm{Re}g>0$, a travelling-roll branch has one nonzero amplitude, of squared modulus $r/g_r$; it is stable to the opposite travelling amplitude when $h_r=\mathrm{Re}h>g_r$. A standing-roll branch has equal squared moduli $r/(g_r+h_r)$, and is stable within this pair when $g_r>h_r$ and $g_r+h_r>0$. The imaginary parts shift angular frequencies. These conditions concern amplitude perturbations in that reduced subspace; orientation and modulation modes still have to be tested. This is the competition of travelling and standing convection rolls.

Near simultaneous stationary and oscillatory thresholds, a codimension-two bifurcation requires retaining both types of mode. After selecting an orientation and fixing a steady spatial phase, a schematic nonresonant steady–Hopf mode interaction isHere $a$ is real, $z$ is the chosen complex oscillatory amplitude, $g_s,h_s$ are real, and $g_o,h_o$ may be complex. Assume $r_s>0$, $g_s>0$ for a radially stable pure steady branch, and $r_o>0$, $g_{o,r}=\mathrm{Re}(g_o)>0$ for a radially stable pure oscillatory branch. Steady convection rolls suppress the oscillatory mode if $r_o-h_{o,r}{r}_s/g_s<0$, where $h_{o,r}=\mathrm{Re}(h_o)$; oscillatory convection rolls suppress the steady mode if $r_s-h_s{r}_o/g_{o,r}<0$. These are transverse tests for branches whose existence and radial stability have already been checked. Mixed states have positive intensities $u=a^2$, $v=|z{|}^2$ solving $g_{s}u+h_{s}v=r_s$, $h_{o,r}u+g_{o,r}v=r_o$. If $D=g_s{g}_{o,r}-h_s{h}_{o,r}\ne 0$, then $u=(r_s{g}_{o,r}-h_s{r}_o)/D$, $v=(g_s{r}_o-h_{o,r}{r}_s)/D$. For positive $u,v$ and the positive self-saturation coefficients above, the intensity Jacobian matrix has negative trace and determinant $4uvD$, so the mixed state is attracting in intensities when $D>0$ and is a saddle equilibrium when $D<0$; phase and other-mode perturbations remain separate tests. Depending on the cross-couplings, one obtains coexistence or competition and bistability; degeneracies or resonances require extra terms. A full travelling/standing-wave competition must retain both $Z_{\pm }$, not just this illustrative one-mode $z$. If the Hopf angular frequency tends to zero, averaging over fast oscillations fails and a double-zero normal form must keep the corresponding two slow variables. A crossing of global minima at distinct wavenumbers is instead an interaction of distinct modes, not automatically the same double-zero problem.

Weakly nonlinear convection rolls can also lose stability spatially. The real Ginzburg–Landau equation $A_T=rA+\xi A_{XX}-g|A{|}^{2}A$, with $r,\xi ,g>0$, has detuned convection rolls with $|A{|}^2=(r-\xi Q^2)/g$. Linear phase modulation gives diffusion coefficient $\xi (r-3\xi Q^2)/(r-\xi Q^2)$, so the Eckhaus instability excludes $Q^2>r/(3\xi )$ even though convection rolls exist up to $Q^2<r/\xi$. Transverse bending and mean-flow couplings impose additional restrictions in the actual rotating layer; the scalar equation is a local longitudinal example, not its complete stability theory.

Most distinctively, rotation breaks mirror symmetry between competing convection roll orientations. For two sets of convection rolls at relative angle $\vartheta$, write ${\dot{A}}_2=r{A}_2-g_0|A_2{|}^2{A}_2-g(\vartheta )|A_1{|}^2{A}_2$ and interchange indices with $\vartheta \mapsto -\vartheta$ for the other equation. The established convection roll $|A_1{|}^2=r/g_0$ is unstable to the new orientation when its linear growth rate $r[1-g(\vartheta )/g_0]$ is positive. The inequality need not be symmetric under $\vartheta \mapsto -\vartheta$ because the imposed rotation supplies handedness. This is the Küppers–Lortz instability: sufficiently rapid rotation can destabilize steady convection rolls immediately above their stationary onset to oblique convection rolls of another orientation. Successive replacements can produce time-dependent convection roll switching or heteroclinic cycles, instead of a stable single convection roll pattern. It requires three-dimensional perturbations even when the original convection roll solution is independent of its axial coordinate.

At finite Prandtl number and stress-free plates, nearly parallel convection rolls can couple resonantly to a slowly damped large-scale Eulerian mean flow. This small-angle instability of rotating convection rolls makes a regular two-roll expansion nonuniform as the angle tends to zero; one must retain the mean-flow mode. It should not be identified automatically with the finite-angle Küppers–Lortz instability or assigned a universal threshold from a one-amplitude equation. Consequently the onset type and selected scale follow from the admissible neutral curves, while persistent convection roll patterns require a separate nonlinear and sideband stability calculation.

For finite Prandtl number and stress-free boundaries, nearly parallel convection rolls can interact with a weakly damped large-scale Eulerian mean flow. The ordinary regular two-roll amplitude expansion becomes nonuniform as the angle tends to zero, so the mean-flow mode must be retained. This mechanism is not automatically the finite-angle Küppers–Lortz instability, and thresholds depend on the stated boundary and parameter regime.

A saturated convection roll is not stable merely because its single-amplitude cubic coefficient is positive. Perturbations include phase sidebands, transverse bending, differently oriented convection rolls, oscillatory modes and large-scale Eulerian mean flow. Their growth rates follow from the corresponding enlarged amplitude equations or the full linearization about the convection roll. Küppers–Lortz instability and small-angle instability of rotating convection rolls are distinct orientation-dependent mechanisms.