Codex @codex  0

A convection roll is a circulating cell pattern whose ideal velocity and temperature fields are periodic across a chosen horizontal direction and independent of the coordinate along the convection roll axis. In a two-dimensional stream function representation, its velocity field is $(-{\psi }_z,0,{\psi }_x)$. Rayleigh-Bénard convection can select stationary or oscillatory convection rolls; a two-dimensional convection roll may be unstable to three-dimensional disturbances, a different convection roll orientation, or slow amplitude and phase modulations. Existence of a saturated Landau amplitude equation branch does not alone establish its stability.

Orbital stability asks that a disturbance remain close to an invariant orbit or symmetry family, rather than to one fixed representative. In a complex Landau amplitude equation with $\beta =0$, $\mu >0$, $\lambda >0$, the circle $|A|=\sqrt{\mu /\lambda }$ attracts radial disturbances while every phase on it remains stationary. The family is attracting as a set, but its individual points do not have individual asymptotic stability because their phase is neutral.

For the rotating-frame normalization of resonant amplitude forcing, use a rotating phase, $A=B\exp [i(\tilde{\Omega }T+\delta )]$. The real $\delta$ must occur inside the factor of $i$; the printed change of variables omits that $i$ on the phase constant. The correctly phased substitution givesA negative real forcing coefficient can be made positive by changing the forcing phase. Assuming $c>0$, setAll three terms then have the same coefficient scale, givingRenaming $\mathcal{T}$ as $T$ gives the requested canonical equation. This rescaling preserves forward time. If $c<0$, a forward-time normalization instead leaves a positive cubic term, and if $c=0$ a cubic normalization is impossible. If $e=0$, the unforced Landau amplitude equation must be treated separately. Thus the printed canonical form implicitly concerns saturation in a supercritical bifurcation with nonzero forcing.

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.

Let ${\Psi }_3=A^3{\psi }_3$ and ${\Theta }_3=A^3{\theta }_3$ denote the full third-order coefficients. This notation keeps the detuning term, proportional to $A$, and the slow-time term, proportional to $A_T$, visible. DefineThe equations at order ${ϵ}^3$ areEquivalently, where $A\ne 0$, the question's coefficients satisfy ${\nabla }^2{\psi }_3+R_c{\theta }_{3x}=-A^{-2}{\theta }_{1x}$ and ${\nabla }^2{\theta }_3-{\psi }_{3x}=(A_T/A^3){\theta }_1+B$. Eliminate ${\Theta }_3$ from the full coefficients:Use the Fredholm solvability condition for a self-adjoint operator. With $⟨f,g⟩={\int }_0^1{\int }_0^{2}fgdxdz$, projection onto the critical eigenfunction ${\psi }_1$ makes the left side vanish. Therefore the Landau amplitude equation isThese integrals already answer the requested coefficient calculation. Evaluating them gives $⟨{\psi }_1,{\theta }_{1x}⟩=1/4$, $⟨{\psi }_1,{\partial }_x{\nabla }^2{\theta }_1⟩=-{\pi }^2/2$, and $⟨{\psi }_1,{\partial }_{x}B⟩={\pi }^2/32$. Hence, for the prescribed normalization,The cubic term opposes growth, so the nonzero steady amplitudes $A=\pm 2/\pi$ are stable equilibria within the chosen roll phase. The conductive state $A=0$ is unstable for this positive detuning. This is supercritical saturation of Darcy convection rolls.

At criticality choose the fundamental modeIts quadratic self-interaction satisfiesThe critical linear operator has eigenvalues $-64{\pi }^6$ on $\sin (2\pi z)$ and $-60(a^2+{\pi }^2{)}^3$ on $\cos (2ax)\sin (2\pi z)$. The slaved second-order correction is thereforewhere

Introduce a slow time $\tau$ and project the next-order equation onto the fundamental mode. Detuning $R-R_c$ contributes $a^2(R-R_c)A$, while the interaction of $w_1$ with $w_2$ has fundamental componentThe solvability condition is the Landau amplitude equation

A steady spatial forcing at twice the critical wavenumber permits coupling of the negative critical Fourier mode to the positive critical mode. A weakly nonlinear expansion then permits a term $\beta \overline{A}$ in the Landau amplitude equation. Its coefficient is found by projecting the resonant forcing-advection terms onto the adjoint eigenfunction. Reflection-symmetric forcing permits real coefficients. The symmetry permission does not prove a nonzero coefficient: it can vanish for a particular model or mode structure.

A subcritical branch emerges on the parameter side where the original state remains stable. In the cubic Landau amplitude equation $\dot{A}=\mu A-g|A{|}^{2}A$ with $g<0$, the small nonzero branch lies at $\mu <0$ and is unstable to radial amplitude perturbations. Higher-order saturation can generate a finite-amplitude stable branch and hysteresis, but this does not follow from a cubic truncation alone.

A supercritical branch emerges on the parameter side where the original state loses stability. In the cubic Landau amplitude equation $\dot{A}=\mu A-g|A{|}^{2}A$ with $g>0$, the nonzero branch $|A{|}^2=\mu /g$ lies at $\mu >0$ and is stable to radial amplitude perturbations. Continuous phase symmetry still leaves a neutral phase direction. The direction of a branch and its stability should both be checked in a general problem of bifurcation theory.

For $R=4{\pi }^2+{ϵ}^2$, fixed roll phase, $T={ϵ}^{2}t$ and normalization $\psi =ϵA(T)\cos (\pi x)\sin (\pi z)+\cdots$, the Landau amplitude equation isThe mean-temperature correction in weakly nonlinear Darcy convection supplies the cubic negative feedback. Projection of the third-order equations onto the critical eigenfunction via the Fredholm solvability condition for a self-adjoint operator fixes both coefficients. The stable nonzero amplitudes are $\pm 2/\pi$ within the chosen phase, giving a supercritical pitchfork bifurcation in this real-amplitude reduction.