Past exam of the mathematics course of the University of Cambridge / 2012 / iii / Paper 77 / 2

The Lorenz model of thermosolutal convection is a Galerkin method truncation retaining a roll velocity and the leading temperature and concentration modes. The variable $a$ measures the roll velocity; $b$ is the temperature mode correlated with the vertical motion, while $c$ is the horizontally averaged temperature-profile distortion generated by that motion. The corresponding concentration modes are $d$ and $e$. The symbol $d$ in this model is a solute amplitude, not the dimensional layer depth. The terms $ab$ and $ad$ change the mean scalar profiles; $a(1-c)$ and $a(1-e)$ express advection of the remaining background gradients.

The positive parameters are the Prandtl number $\sigma =\nu /\kappa$, the ratio $\tau ={\kappa }_s/\kappa$ of solutal diffusivity to thermal diffusivity, and the geometric vertical-mode decay ratio $\varpi$. In the usual fundamental-roll truncation with time normalized to the thermal roll-mode decay rate, $\varpi =4{\pi }^2/(k^2+{\pi }^2)$; treating it as an independently specified positive truncation parameter leaves the arguments below unchanged.

The parameter $r$ is the destabilizing thermal Rayleigh number divided by the free-slip critical-mode value, and $r_s>0$ measures stabilizing solutal buoyancy in the same thermal-diffusive normalization. For example, using a positive bottom-to-top stabilizing solute difference $\Delta S$ and density coefficient ${\beta }_S$,If a solutal Rayleigh number is instead defined with ${\kappa }_s$ in its denominator, the coefficient $r_s$ here is $\tau$ times that normalized number. Thermal buoyancy drives the roll; solutal buoyancy opposes it and diffuses at a different rate. This is the essential mechanism of thermosolutal convection.

Set the time derivatives to zero. The scalar equations giveFor the nonzero convection branch, division of the velocity equation by $a$ givesIts zero-amplitude limit is $r_e=1+r_s/\tau$. Differentiation yieldsFor $0<\tau <1$, the sharp condition for a subcritical bifurcation at the steady onset is $r_s>{\tau }^3/(1-{\tau }^2)$. The stronger bound printed in the question is sufficient, since $(1-\tau {)}^2<1-{\tau }^2$. Under it the branch initially runs toward smaller $r$, and the unique saddle-node bifurcation hasIt is a genuine positive-amplitude minimum: $r^{\prime \prime }(x)=2r_s\tau (1-{\tau }^2)/({\tau }^2+x{)}^3>0$. Put $h=\sqrt{r_s\tau }$ and $j=\sqrt{1-{\tau }^2}$. At the minimum, $1+x_{\ast }=j(j+h)$ and ${\tau }^2+x_{\ast }=hj$, soThese results describe the steady fold of a thermosolutal Lorenz model. They do not by themselves assert stability against every oscillatory disturbance.

Take $\beta \ge 0$ without loss of generality, since its definition fixes only ${\beta }^2$. With $r_s={\beta }^2\tau$, the steady-branch minimum expands as $r_{\min }=1+2\beta \tau +O({\tau }^2)$. For fixed positive $\sigma$, the source's supplied approximation therefore givesThus the requested inequality follows from the supplied approximation, with equality at first order when $\beta =\delta$. If negative $\beta$ is retained, replace it by $|\beta |$ throughout.

There is an independent source-consistency issue with calling that supplied expression the oscillatory threshold of the displayed model. Linearizing about conduction, the $c,e$ modes decay separately and the $a,b,d$ characteristic polynomial isWriting it as ${\lambda }^3+A{\lambda }^2+B\lambda +C$, a nonzero imaginary pair requires $C=AB$ and $B>0$, by the Routh-Hurwitz stability criterion. Direct algebra gives the oscillatory threshold of a thermosolutal Lorenz modelConsequently, with $\delta =\sigma /(1+\sigma )$,The two factors of $\delta$ in the supplied threshold are interchanged. For example $\sigma =1$, $\beta =2$, $\tau =10^{-3}$ gives a genuine marginal imaginary pair at $r_H=1.004003$, whereas the supplied first-order expression gives $1.0085$. The intended comparison remains true after correction:This comparison applies to a physical Hopf bifurcation when ${\omega }_H^2>0$ and to an interior steady fold when $a_{\ast }^2>0$. For fixed nonzero $\beta$ and fixed $\sigma >0$, both conditions hold for sufficiently small positive $\tau$. At a vanishing squared difference, neglected terms must be examined before drawing a higher-order conclusion.