Past exam of the mathematics course of the University of Cambridge 2014 iii Paper 76 3 Solution Created 2026-10-03 Updated 2026-10-06
Rotating Rayleigh-Bénard convection combines buoyancy-driven instability with the Coriolis force. Consider a plane layer of depth 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 is kinematic viscosity and thermal diffusivity. In thermal-diffusion time units, linear perturbations satisfyRotation does no direct mechanical work, since , but couples vertical motion to vertical vorticity and changes the damping and oscillation balance.
For horizontal wavenumber and vertical mode , put and . Use , and vertical vorticity with growth rate . 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 , givingRotation raises this stationary threshold. An oscillatory neutral root has with . Real-imaginary separation gives the oscillatory neutral curve of rotating convectionThis branch is admissible only when , requiring and sufficiently strong rotation. The restoring Coriolis force coupling permits an inertial/thermal oscillation whose phase-lagged buoyancy can overcome dissipation. At large , 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 . 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 . Minimizing the stationary curve givesThus the stationary neutral curve of rotating convection hasThe physical horizontal wavelength is , hence decreases as ; 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 , giving the same Taylor number exponent at fixed positive . Where its frequency remains admissible,Equality is , whose positive root is approximately . Accordingly, for sufficiently rapid rotation in this free-slip problem, selects oscillatory onset and stationary onset. The weaker condition 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 term restores the fast frequency if desired. The real parts govern amplitude saturation and the imaginary parts give nonlinear frequency shifts. Write , . For a travelling wave from a supercritical bifurcation with only one amplitude nonzero, with , and the competing wave's growth rate is . It is amplitude-stable against that competitor when . A standing wave has equal intensities ; provided this is positive, its intensity-difference mode is stable when . 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 , a leading competition system hasA pure roll has with . An infinitesimal new roll at relative angle grows atFor sufficiently strong rotation in appropriate boundary and Prandtl number regimes, some finite oblique angle has , so a steady roll is unstable arbitrarily close above its stationary onset. Rotation is handed and allows , 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.
Past exam of the mathematics course of the University of Cambridge 2016 iii Paper 331 1 a Solution Created 2026-10-03 Updated 2026-10-06
Introduce the thermal expansion coefficient , so the linear equation of state is . With gravity , 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 , where .
The Boussinesq approximation retains temperature-dependent mass density in buoyancy while replacing it by in inertial coefficients; it requires . 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 .
Use the thermal-diffusion scaling of a convection layerThe Rayleigh number and Prandtl number areThus the nondimensional perturbation equations areHere is kinematic viscosity, is thermal diffusivity, and positive denotes destabilizing heating from below for .