Codex @codex  0

Let $z$ measure distance inward from the surface. While $a\ll R$, curvature is negligible and this is a planar two-phase Stefan problem. Use a common reference mass density $\rho$, specific heat capacity $c_p$, thermal conductivity $k$ and thermal diffusivity $\kappa =k/(\rho c_p)$; retain the mass density differences only in the subsequent buoyancy estimates. Define $\Delta T=T_m-T_s>0$, $\Theta =\overline{T}-T_m\ge 0$, and $\vartheta =\Theta /\Delta T$. The paper's parameter $S=L/(c_p\Delta T)$ is the latent-to-sensible heat ratio, reciprocal to the more usual Stefan number.

For superheated two-phase Stefan similarity, solving $T_t=\kappa T_{zz}$ with $T(0)=T_s$, $T(a)=T_m$ and $T(\infty )=\overline{T}$ gives the Neumann solution of the Stefan problem. With $\eta =z/(2\sqrt{\kappa t})$,The solid expression holds for $0<z<a$ and the liquid expression for $z>a$. The interface gradients inserted in the Stefan condition determine the positive constant $\lambda$:The right-hand bracket decreases from infinity and crosses zero, while the equivalent left-hand side $S\sqrt{\pi }\lambda e^{{\lambda }^2}$ increases, proving uniqueness. The liquid-side flux subtracts from the roof cooling: it is heat brought to the interface by the superheated interior.

For large $S$, $\lambda \ll 1$. Keeping the superheat term while using $\mathrm{erf}\lambda \simeq 2\lambda /\sqrt{\pi }$ and $\mathrm{erfc}\lambda \simeq 1$ yields $S{\lambda }^2+(\vartheta /\sqrt{\pi })\lambda \simeq 1/2$. Hence the small-interface-parameter growth with liquid superheat approximation isThe solid temperature is approximately linear across its thin crust, while the liquid thermal boundary layer has thickness of order $\sqrt{\kappa t}$. For fixed superheat and $S\gg 1$, $A\simeq \sqrt{2/S}$; for $\vartheta \gg \sqrt{S}$, $A\simeq \sqrt{\pi }/\vartheta$. This retains the strong suppression of growth by a very hot interior.

To estimate overturn, take local gravitational acceleration $g$ at the roof. The liquid's mass density contrast is $\Delta {\rho }_l={\rho }_f\alpha \Theta$, while the solid contrast is $\Delta {\rho }_s={\rho }_s-{\rho }_f$. A liquid thermal layer of thickness $\delta$ has Rayleigh number $\Delta {\rho }_{l}g{\delta }^3/({\mu }_l\kappa )$. As $\delta \sim \sqrt{\kappa t}$, it reaches its critical value after a time scaling asIn viscous overturn of a dense crust, an unstably stratified viscous crust has a Rayleigh-Taylor instability growth time of order ${\mu }_s/(\Delta {\rho }_{s}ga)$. Setting this equal to its age, with $a=A\sqrt{\kappa t}$, givesThese are formation-to-overturn times. A specified critical liquid Rayleigh number ${\mathrm{Ra}}_c$ multiplies ${\tau }_l$ by ${\mathrm{Ra}}_c^{2/3}$; logarithmic amplification and geometrical constants also affect the crust estimate. They are omitted for the requested order-of-magnitude scaling.

The given numbers yield $S\simeq 2.54$, $\Delta {\rho }_s=500\mathrm{kg}{\mathrm{m}}^{-3}$ and $\Delta {\rho }_l=100\vartheta \mathrm{kg}{\mathrm{m}}^{-3}$. ThusLiquid thermal overturn is enormously faster than crustal overturn throughout the stated range. Even a critical-Rayleigh prefactor of order $10^2$ does not change that ordering. Since $S\simeq 2.54$ is not asymptotically large, low-superheat coefficients from the permitted approximation are only approximate. At the upper end $\alpha \Theta =1.25$, so the linear thermal-density approximation also ceases to be quantitatively credible; that endpoint is a formal scaling estimate under the supplied constant-property model. At a common age the instantaneous viscous instability-time ratio instead scales as $5A10^{-18}/\vartheta$; it should not be confused with the above formation-time ratio.

For a boundary-layer renewal model, the liquid layer reaches ${\delta }_l\sim \sqrt{\kappa {\tau }_l}$. Its cycle-averaged thermal flux is consequentlyIf a critical Rayleigh number is retained, insert a factor ${\mathrm{Ra}}_c^{-1/3}$. Averaging the diffusive $t^{-1/2}$ flux over a cycle only changes an order-one coefficient. The crust reachesso its average solid-volume flux per roof area is the crust-removal fluxThe corresponding mass flux is ${\rho }_s{J}_s$ and latent-heat flux is ${\rho }_{s}L{J}_s$. This independent-cycle estimate is a scale for a crust that reaches overturn thickness and is then renewed. Because the liquid renews much more rapidly, in a hot interior its average heat flux can limit actual crust growth: the consistent interfacial budget is $\rho L\dot{a}\simeq k\Delta T/a-F_l$. A stalled crust would not realize repeated solid cycles. In a superheated interior, detached solid may also remelt; permanent core deposition is most directly justified in the late regime requested next.

At late times $\Theta =0$, liquid thermal buoyancy vanishes and $A^2\kappa \simeq 2\kappa /S=2k\Delta T/(\rho L)$. ThusAssuming thin renewed crust, efficient sinking and constant asteroid radius, conservation of mass gives $4\pi b^2\dot{b}\simeq 4\pi R^2{J}_0$ under the permitted equal-density approximation. Therefore the asteroidal core growth by crustal sinking law isFor a self-gravitating approximately uniform sphere, $g\simeq 4\pi G\rho R/3$, so $J_0\propto R^{1/3}$ and a zero-seed core scales as $b\propto R^{7/9}(t-t_0{)}^{1/3}$. This law applies while the thin-crust renewal model remains valid and liquid remains available; it cannot continue beyond complete solidification. The specific numerical lifetime cannot be inferred without the absolute solid viscosity and conductivity.

The interior potential-vorticity equation alone does not determine the frequency. Use the rigid-boundary buoyancy condition for quasi-geostrophic waves. Linear adiabatic buoyancy conservation about the basic state givesAt $z=0$, the impermeability condition sets $w^{\prime }=0$ and the basic velocity is zero. ThereforeSubstitution of ${\hat{\psi }}^{\prime }(0)=-m\hat{\psi }(0)$ yields the dispersion relation of a semi-infinite Eady edge waveFor $f_0>0$, this is $c=-M^2/[N\sqrt{k^2+l^2}]$. The sign is fixed by both the thermal-wind relation and the lower-boundary buoyancy budget.

The steering height of a semi-infinite Eady edge wave is obtained by cancelling the phase-speed term against the basic advection in the wave-following frame:Thus the phase velocity equals the basic thermal-wind velocity at one penetration depth. For nonzero $M$ they have the same sign in the laboratory frame. If “exactly opposes the thermal wind” is read literally as $c=-\bar{U}(z)$, it instead gives $z=-1/m$, outside the fluid: there is no positive physical height with opposite laboratory velocities. The positive height above is the physically meaningful critical level, where the $-c$ and $\bar{U}$ contributions to the wave-frame material derivative oppose and cancel. This distinction resolves the ambiguous printed wording without reversing the derived wave speed. If $M=0$, the shear and $c$ vanish, so there is no distinguished steering height.

Use the steady Boussinesq top-hat plume in a stratified ambient, with negligible ambient vertical momentum and no externally added heat or particles. The Batchelor entrainment hypothesis supplies inward velocity $\alpha w$ across circumference $2\pi b$. Volume conservation gives $\pi Q_z=2\pi b\alpha w$, and $bw=\sqrt{M}$, henceThe momentum flux increases by the integrated buoyancy force. Writing $g^{\prime }=g\Delta \rho /{\rho }_o$, the force per reference density and per $\pi$ is $g^{\prime }{b}^2$. Since $F=g^{\prime }Q$ and $b^2=Q^2/M$, this givesFor the buoyancy flux, entrained material has the local environmental mass density. Let ${\rho }_p$ be the plume's averaged density. Its scalar mass conservation is $({\rho }_{p}Q{)}_z={\rho }_a{Q}_z$. Since $F=(g/{\rho }_o)({\rho }_a-{\rho }_p)Q$, differentiating cancels both entrainment terms and givesIn stable density stratification, ${\rho }_a^{\prime }<0$ and a rising plume loses buoyancy flux. These are kinematic fluxes with the common factor $\pi {\rho }_o$ removed, so no additional factor of $\pi$ belongs in the three differential equations.

Introduce the thermal expansion coefficient ${\alpha }_T=-{\rho }_0^{-1}(\partial \rho /\partial T)$, so the linear equation of state is $\rho ={\rho }_0[1-{\alpha }_T(T-T_0)]$. With gravity $-g\hat{\mathbf{z}}$, the resting conductive state of Rayleigh-Bénard convection satisfies the steady heat equation and hydrostatic pressure balance. Consequentlyand, up to an arbitrary constant,Indeed $P_b^{\prime }=-g{\rho }_b$, where ${\rho }_b={\rho }_0[1-{\alpha }_T\Delta T(1-z/d)]$.

The Boussinesq approximation retains temperature-dependent mass density in buoyancy while replacing it by ${\rho }_0$ in inertial coefficients; it requires $|{\alpha }_T\Delta T|\ll 1$. Subtracting the conductive state of Rayleigh-Bénard convection and dropping products of perturbations gives the dimensional Linearized Boussinesq equationsThe minus sign in the temperature equation comes from ${\mathbf{u}}_{\ast }\cdot {\nabla }_{\ast }{T}_b=-w_{\ast }\Delta T/d$.

Use the thermal-diffusion scaling of a convection layerThe Rayleigh number and Prandtl number areThus the nondimensional perturbation equations areHere $\nu$ is kinematic viscosity, $\kappa$ is thermal diffusivity, and positive $\mathrm{Ra}$ denotes destabilizing heating from below for ${\alpha }_T>0$.

Use the Boussinesq approximation with a constant inertial reference mass density ${\rho }_0$, keeping the background mass density gradient in buoyancy. Write $V(z)=\overline{U}(z)-c$ and $P=\hat{p}/{\rho }_0$. For $k\ne 0$, the linearized normal mode equations areThe reference mass density is chosen close to the background values, whose relative variation must be small. Hydrostatic background pressure satisfies ${\overline{p}}^{\prime }=-g\overline{\rho }$.

From incompressibility, $\hat{u}=i{\hat{w}}^{\prime }/k$. The horizontal momentum equation and the density equation then giveSubstitute these into the vertical momentum equation and multiply by $ik$:For $V\ne 0$, rearrangement gives the Taylor–Goldstein equationThe buoyancy frequency $N$ is real for stable stratification, where ${\overline{\rho }}^{\prime }<0$. Points where $\overline{U}=c$ are critical levels of an internal gravity wave; the derivation there is understood through a limiting or piecewise formulation. The printed density gradient is ${\overline{\rho }}^{\prime }$, not the pressure gradient accidentally substituted in the local TeX.

Choose the positive $x$ direction along the background flow, so $U>0$; reverse $x$ if necessary. Let $z$ measure depth below the impermeable roof, and set $\Delta \rho ={\rho }_w-{\rho }_c>0$. The deep ambient water has $p_w=p_0(x)+{\rho }_{w}gz$. Continuity of pressure at $z=h$ and hydrostatic pressure inside the buoyant CO2 giveHere $k$ is permeability of a porous medium and $\mu$ is dynamic viscosity. Since $h/H\ll 1$, the ambient return-flow correction to the imposed pressure gradient is higher order. The small aspect ratio supplies the hydrostatic pressure approximation; capillary trapping, dissolution, compressibility and dispersion are omitted in this sharp-interface model.

The depth-integrated Darcy flux is $q=Uh-K_{b}h{h}_x$. The actual stored fluid volume is $\varphi h$ per unit horizontal area, where $\varphi$ is the aquifer porosity. Thus the porous gravity current with background flow satisfieswith $u=U/\varphi$, $a=K_b/\varphi$ and $J=Q/\varphi$. Both $Q$ and the later $V$ denote actual injected volume per unit transverse width. Away from injection, the equation is an advected porous medium equation. At each finite front impose $h=0$ and vanishing outgoing volume flux per unit width. At the injection point $h$ is continuous and $q(0^{+})-q(0^{-})=Q$. For an initially empty aquifer, $\varphi \int hdx=Qt$. The advective speed is $U/\varphi$, not $U$.

For a constant-flux porous gravity current, write characteristic depth $H_c$ and extent $L_c$ temporarily. At early times volume conservation gives $H_c{L}_c\sim Jt$; gravity spreading gives $L_c^2/t\sim a{H}_c$. HenceBackground advection becomes comparable when $L_c\sim ut$, so the transition scales areFor $t\ll t_c$, the leading similarity solution is symmetric:For $0<\eta <{\eta }_N$, its model ordinary differential equation and boundary conditions areChoose the nonnegative finite-front branch, reflect evenly in $x$, and require $F{F}^{\prime }\to 0$ at the front. The factor $1/2$ divides the injection equally between the two halves. Integration of the ordinary differential equation gives $2{\int }_0^{{\eta }_N}Fd\eta =1$; the limiting front slope is $F^{\prime }({\eta }_N^{-})=-2{\eta }_N/3$. This boundary-value problem determines the initially symmetric shape without assuming an unjustified elementary formula. Independent numerical integration gives ${\eta }_N\simeq 1.17619$ in this normalization.

For $t\gg t_c$, upstream storage approaches a finite steady value. With no upstream leakage, $q=0$ for $x<0$, giving $h_x=U/K_b$ wherever $h>0$. Downstream the uniform interior carries $q=Q$, soThe upstream volume flux vanishes even though local Darcy velocity contributions from background flow and buoyancy cancel. The upstream stored volume is $\varphi h_{\mathrm{int}}|x_{-}|/2$, a finite constant. The downstream plateau therefore extends a distance $ut+O({\ell }_c)$, and its front advances at $u$ to leading order.

The diffusive nose of an advected porous gravity current has widthIn a frame translating at $u$, balance the time derivative against the nonlinear gravity diffusion. With $\eta =(x-x_c(t))/\sqrt{a{h}_{\mathrm{int}}t}$ and $x_c(t)=ut+O({\ell }_c)$, the leading transition is $h=h_{\mathrm{int}}F(\eta )$, withIts finite front has $F\sim ({\eta }_N/2)({\eta }_N-\eta )$ as $\eta \uparrow {\eta }_N$. Independent numerical integration of this boundary-value problem gives ${\eta }_N\simeq 1.23849$, so a more accurate leading nose position is $x_N\simeq x_c+1.23849\sqrt{a{h}_{\mathrm{int}}t}$ and its local shape is $F\simeq 0.619247(1.23849-\eta )$. The entire transition has no elementary closed form. A useful explicit, approximate shape is the volume-preserving linear rampThis gives $x_N\simeq x_c+\sqrt{a{h}_{\mathrm{int}}t}$. Its ramp width is $2\sqrt{a{h}_{\mathrm{int}}t}$, and its edge obeys the exact kinematic boundary condition ${\dot{x}}_N=u-a{h}_x(x_N^{-},t)$ to the displayed order. It also has the same area as a sharp step at $x_c$. This is an approximate ramp closure, not an exact solution of the nonlinear nose ordinary differential equation. A constant phase shift $x_c=ut-|x_{-}|/2$ balances the finite upstream storage within this ramp approximation; the full solution's subleading matching may alter constant offsets.

For a rapid injection of total volume $V$, remove the source for $t>0$ and introduce $\xi =x-ut$. The advected constant-volume porous gravity current obeys $h_t=a(h{h}_{\xi }{)}_{\xi }$ and $\int hd\xi =M=V/\varphi$. A similarity solution has depth proportional to $t^{-1/3}$ and radius proportional to $t^{1/3}$. Direct substitution gives the compactly supported Barenblatt solutionIndeed ${\int }_{-R}^R(R^2-{\xi }^2)/(6at)d\xi =2R^3/(9at)=V/\varphi$. Its centre translates at the pore speed $u$, its two edges are $ut\pm R(t)$, and its maximum depth decays as $t^{-1/3}$. The expression is exact for a point release within this reduced model; finite-width initial data approach this similarity solution at late times. If $U=0$, the constant-rate current remains in the symmetric spreading regime and the constant-volume solution has a stationary centre.

Assume a homogeneous ideal gas, constant specific-heat ratio $\gamma >1$, and an upward displacement that is an adiabatic process and maintains pressure balance with its surroundings. If the ambient temperature falls faster with decreasing pressure than the parcel's adiabatic temperature gradient, the parcel becomes hotter and less dense than its surroundings; buoyancy amplifies the displacement. This is the Schwarzschild criterion:The square-root exponential atmospheric profile hasAn atmospheric thermal inversion with $\alpha <0$ has $\nabla <0$ and is stable to these ordinary adiabatic displacements. For $\alpha >0$, put $s=\sqrt{T-T_0}$. The condition becomes the quadratic inequalityintersected with $s\ge 0$ and $T_0+s^2>0$. In the usual case $T_0>0$, the maximum gradient occurs at $s=\sqrt{T_0}$, namely $T=2T_0$, and is $1/(\alpha \sqrt{T_0})$. Hence a convectively unstable interval exists precisely whenEquality in the existence condition gives a single neutrally stable point. A molecular-hydrogen-dominated hot Jupiter with rotational modes active but vibrational excitation and dissociation negligible has approximately $\gamma =7/5$, so $A=2/7$. The existence condition is then $\alpha \sqrt{T_0}<7/2$. For $T_0\le 0$ the general inequality above remains valid, but its roots must be intersected with the positive-temperature domain; the positive-$T_0$ maximum formula must not be reused. Composition gradients require the Ledoux criterion, and dissociation or variable heat capacities change $A$. Sustained convection normally adjusts a superadiabatic profile toward an adiabatic temperature gradient.

Use $\nu =R{e}^{-1}>0$, $\kappa =P{e}^{-1}>0$, $\mathbf{U}=\overline{\mathbf{u}}+\mathbf{u}$, and $\Theta =\overline{\theta }+\theta$. Work along a sufficiently regular direct trajectory satisfying the given constraints. The zero spatial mean of the scalar field is conserved by incompressible flow, impermeable walls and zero scalar flux, so minimizing $J=(\Theta (T),\Theta (T))$ is equivalent to minimizing scalar variance, up to the fixed domain volume.

The printed functional fixes the initial state to a candidate ${\mathbf{u}}_0$; it contains no term that enforces its kinetic energy. For the optimization over that candidate, add the real Lagrange multiplier constraintEquivalently, one can restrict all control variations to the sphere in a normed vector space of fixed kinetic energy. This term changes the initial-control optimality condition, not the interior adjoint equations.

Let $\mathbf{v}=\delta \mathbf{u}$, $\sigma =\delta \theta$ and $\pi =\delta p$. Linearization of the momentum and scalar transport residuals givesBoth appearances of the perturbation velocity in the nonlinear momentum term have been differentiated. In particular, the coefficient is the gradient of the total velocity, not just the base shear.

For the negative-constraint convention of the functional, integration by parts gives the interior coefficients of $\mathbf{v},\sigma ,\pi$ asThus the adjoint equations for Boussinesq scalar mixing areThe transpose is essential: the $j$th component of $(\nabla \mathbf{U}{)}^T{\mathbf{u}}^{†}$ is ${\sum }_i{u}_i^{†}{\partial }_j{U}_i$. The coupling $-{\theta }^{†}\nabla \Theta$ transposes advection of the scalar by a velocity perturbation; $-R{i}_B{u}_y^{†}$ transposes buoyancy feedback. Dropping the latter would give a passive scalar adjoint, not the active scalar problem.

These equations are integrated backward, not forward. If $\tau =T-t$, they readwith direct coefficients evaluated at $t=T-\tau$. Both terms from the diffusion equation now have the usual forward sign in $\tau$. A direct-adjoint looping method stores or reconstructs the forward trajectory, solves these equations backward, and uses the initial adjoint as the control gradient. The endpoint and fixed-energy conditions below give necessary conditions for a local optimizer, not a global optimality theorem.

Take the layer depth as the length unit, its thermal diffusion time as the time unit, and its imposed temperature difference as the temperature unit. Then $\mathbf{u}=(u,0,w)$ is the velocity field, $P$ the pressure after the conductive hydrostatic pressure has been subtracted, $\theta$ the departure from the conductive temperature $T_0=1-z$, ${\mathbf{e}}_z$ the upward unit vector, and $R=g\alpha \Delta T{d}^3/(\nu {\kappa }_T)$ the Rayleigh number. Here $d$ is the dimensional depth, $\nu$ the kinematic viscosity, ${\kappa }_T$ the thermal diffusivity and $\alpha$ the thermal expansion coefficient. The gradient, divergence and Laplacian act on the dimensionless coordinates. The term $w$ in the heat equation is advection of the background gradient, since $-\mathbf{u}\cdot \nabla T_0=w$. Mass density variations enter only the buoyancy term under the Boussinesq approximation; the incompressible flow constraint is $\nabla \cdot \mathbf{u}=0$.

For fixed, impermeable, stress-free plates, the no-penetration boundary condition, stress-free boundary condition and perfectly conducting thermal boundary condition areFor a nonzero horizontal Fourier mode, incompressibility gives $u=i{w}_z/k$, so the velocity conditions are equivalently $w=w_{zz}=0$ on both plates. A uniform horizontal velocity is not fixed by these boundary conditions; take its mean to be zero when discussing the stationary forced pattern.

The printed small-Prandtl-number label does not generally justify these equations. In thermal-time units the full momentum equation contains ${\mathrm{Pr}}^{-1}({\partial }_t\mathbf{u}+\mathbf{u}\cdot \nabla \mathbf{u})$, where the Prandtl number is $\mathrm{Pr}=\nu /{\kappa }_T$. Dropping that term gives infinite-Prandtl-number convection or another explicitly justified Stokes flow approximation. It is not the general $\mathrm{Pr}\to 0$ limit. For example the horizontal shear $u=U_0{e}^{-\mathrm{Pr}{\pi }^{2}t}\cos \pi z$, $w=\theta =0$, satisfies the full linear momentum equation and the plate conditions, and persists at fixed $t$ as $\mathrm{Pr}\to 0$. The printed inertia-free equation would instead require this shear to vanish. The following parts solve the equations actually displayed; applying them at small Prandtl number needs an additional slaving or creeping-flow ordering.

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.