Past exam of the mathematics course of the University of Cambridge / 2018 / iii / Paper 321 / 3

Let $\mathbf{\zeta }=\nabla \times \mathbf{u}$ and $\mathbf{\omega }=\mathbf{\zeta }+2\Omega {\mathbf{e}}_z$ be the relative vorticity and absolute vorticity. The nonlinear acceleration identity rewrites the momentum equation asTaking the curl removes the pressure, kinetic energy gradient and shearing-sheet tidal potential. ThusBoth $\nabla \cdot \mathbf{u}$ and $\nabla \cdot \mathbf{\omega }$ vanish. The vorticity equation is consequentlyFor the absolute-vorticity flux tensor $T_{ij}=u_i{\omega }_j-{\omega }_i{u}_j$, use the flux convention $(\nabla \cdot \mathbf{T}{)}_j={\partial }_i{T}_{ij}$, with the first index denoting transport direction. This givesThe scalar $\theta$ is a buoyancy displacement variable, so dimensional consistency gives it units of length; the printed potential-temperature terminology presupposes this normalization. Its curl source directly forces only the $y,z$ components of absolute vorticity. Hence radial vorticity conservation in a shearing sheet takes the local flux formIts volume integral is constant when the net boundary flux vanishes. However $D_t{\omega }_x=\mathbf{\omega }\cdot \nabla u_x$ still permits vortex stretching; the radial component is not generally a materially conserved scalar in three dimensions. Using divergence on the second index of the stated tensor would reverse the transport sign, so the tensor convention matters.

For the Keplerian shearing sheet background $\mathbf{U}=-3\Omega x{\mathbf{e}}_y/2$, the only spatial derivative is ${\partial }_x{U}_y=-3\Omega /2$. Since $U_x=0$ and all fields are independent of $y$, $\mathbf{U}\cdot \nabla \mathbf{U}=0$ and $\nabla \cdot \mathbf{U}=0$.

The Coriolis acceleration exactly balances the tidal term:With constant pressure and $\theta =0$, the buoyancy displacement variable equation also has zero advection, zero source $U_x$ and zero diffusion. Thereforeis a steady solution. Its absolute vorticity is ${\mathbf{\omega }}_0=(\Omega /2){\mathbf{e}}_z$, since the background relative vorticity is $-3\Omega {\mathbf{e}}_z/2$.

For an axisymmetric vertical mode of a shearing sheet in this Boussinesq approximation, the perturbations have no horizontal gradients, so background-shear advection vanishes. The linear Coriolis acceleration and perturbation advection of the background shear givewhile thermal conduction and incompressibility giveHere $\beta$ is the thermal diffusion rate of a Fourier mode, not the cooling exponent in Question 2. Since $k\ne 0$, incompressibility forces $u_z^{\prime }=0$, and the vertical momentum equation then forces $P^{\prime }=0$ for this nonzero-wavenumber amplitude.

The three remaining amplitudes satisfyA nontrivial solution requires the determinant to vanish, giving the convective overstability dispersion relationNo division by $s$ or $s+\beta$ was needed, so the thermal energy mode of a shearing sheet is retained. The notation $N^2$ denotes a signed buoyancy frequency squared and can be negative for an adverse stratification; it must not be restricted to the square of a real positive frequency.

Treat $n^2=N^2/{\Omega }^2$ as a signed small parameter, with $|n^2|\ll 1$, as printed in the original PDF. At $n^2=0$, the dispersion relation factors asThus the two epicyclic motion modes have $s_0=\pm i\Omega$, while the thermal energy mode of a shearing sheet hasIt decays on the thermal diffusion time ${\beta }^{-1}$. For completeness, its next correction is $s_E=-\beta +\beta N^2/({\beta }^2+{\Omega }^2)+O(\Omega n^4)$ at fixed $\beta /\Omega$.

For either epicyclic motion root, the implicit function theorem applied to the dispersion relation givesHenceSince $\beta >0$, the oscillations undergo overstability precisely whenThe convective overstability growth rate tends to zero both for very slow diffusion and for very fast diffusion. Differentiating $\beta /({\beta }^2+{\Omega }^2)$ shows that its maximum occurs at $\beta =\Omega$, corresponding to $k^2=\Omega /\xi$. ThereforeThis is the growth rate of the perturbation amplitude; a quadratic perturbation energy grows at twice that rate. The convective overstability arises when an adverse radial specific entropy gradient couples to epicyclic motion with a finite thermal conduction lag.