Past exam of the advanced fluid dynamics course of the University of Oxford / 2018 / 1

Let $\hat{\mathbf{u}}$ and $\hat{p}$ be the complex plane wave amplitudes and set $k_x=k_{\perp }$, $k_z=k_{\parallel }$. The linearized Euler equations in the rotating frame areHere $p$ is pressure per unit mass density, including the centrifugal potential. The last line is incompressibility, not an additional equation of state.

Taking the curl eliminates pressure. For constant $\mathbf{\Omega }$ and a divergence-free vector field, $\nabla \times (2\mathbf{u}\times \mathbf{\Omega })=2(\mathbf{\Omega }\cdot \nabla )\mathbf{u}$. Thus the inertial-wave vorticity equation givesApply $\mathbf{k}\times$ once more and use $\mathbf{k}\times (\mathbf{k}\times \hat{\mathbf{u}})=-k^2\hat{\mathbf{u}}$. For a nonzero oscillatory amplitude,so the inertial-wave dispersion relation isThe restoring influence is the Coriolis acceleration. In particular $|\omega |\le 2|\Omega |$. For $k_{\parallel }=0$, the projected linear dynamics have zero frequency; these are steady balanced disturbances rather than the nonzero-frequency branch used in the division above. The spatially uniform $k=0$ mode must be treated separately from this dispersion relation.

The group velocity, which gives wave-packet and energy propagation, is the gradient of the dispersion relation:For $\mathbf{k}=(k_{\perp },0,k_{\parallel })$ this becomesThus for $ϵ=|k_{\parallel }|/k_{\perp }\ll 1$, energy travels almost along the rotation axis, with a small horizontal component. By contrast, the normal-to-phase-front phase velocity ${\mathbf{v}}_{\mathrm{ph}}=\omega \mathbf{k}/k^2$ is almost horizontal. Indeed $\mathbf{k}\cdot {\mathbf{v}}_g=0$, so these phase and group directions are perpendicular. Specifying which meaning of propagation is intended matters for an inertial wave.

For the wave polarization, incompressibility and the $y$ equation in part a giveThe horizontal component parallel to ${\mathbf{k}}_{\perp }$ is small, but $|{\hat{u}}_y|\simeq |{\hat{u}}_z|$: the other horizontal and vertical components are in quadrature. More exactly, introduce the orthonormal transverse basisThen $\hat{\mathbf{u}}=A({\mathbf{e}}_1-\sigma i{\mathbf{e}}_2)$. The physical velocity at a fixed point traces a circle in the plane perpendicular to $\mathbf{k}$; its handedness reverses between branches. In the anisotropic limit this plane is nearly vertical. The group velocity and polarization of an inertial wave therefore have quite different geometric roles.

Let $L_{\perp }=k_{\perp }^{-1}$ and let $U$ be the common size of horizontal and vertical velocity perturbations, as found in part b. The linear inertial wave time and nonlinear advection time areThe critical balance in a rapidly rotating fluid is thereforeThe Rossby number is of order $ϵ$. Horizontal geostrophic flow balances the leading pressure gradient against Coriolis acceleration, soThese are dimensionally comparable statements, rather than comparisons of a dimensional frequency with a dimensional velocity. Equivalently, use $L_{\perp }$ and $(2\Omega {)}^{-1}$ as units. ThenThus rotation stays $O(1)$ in these units, while the motions evolve on the slower $O({ϵ}^{-1})$ time. Vertical advection $u_{\parallel }{\partial }_z$ is smaller than horizontal advection by $ϵ$.

Write ${\mathbf{u}}_{\perp }={\mathbf{u}}_{\perp }^{(0)}+{\mathbf{u}}_{\perp }^{(1)}+\cdots$, with ${\mathbf{u}}_{\perp }^{(1)}=O(ϵU)$, and $p=p_0+p_1+\cdots$. At leading horizontal order,Define the stream function $\Phi =p_0/(2\Omega )$, including its horizontally uniform part in this convention. ThenThe next-order incompressibility condition is not zero horizontal divergence:This term has size $ϵk_{\perp }U$ and is essential: its product with the large rotation rate is a leading contribution to the vertical vorticity equation.

To see that contribution explicitly, write $D_{\perp }={\partial }_t+{\mathbf{u}}_{\perp }^{(0)}\cdot {\nabla }_{\perp }$. At the first horizontal inertial order,Taking its vertical curl eliminates $p_1$. The two-dimensional advection identity and the two divergence orders giveEquivalently the ageostrophic flow satisfieswhose divergence is $-D_{\perp }\zeta /(2\Omega )$. This explicitly verifies that the next-order divergence has the required sign.

The vertical momentum equation, at the same slow order, isFinally ${\mathbf{u}}_{\perp }^{(0)}\cdot {\nabla }_{\perp }f=\{\Phi ,f\}$, where the canonical two-dimensional Poisson bracket is $\{f,g\}=f_x{g}_y-f_y{g}_x$. The reduced equations for an anisotropic rotating fluid areThe omitted vertical-advection, ageostrophic-advection and relative-vorticity stretching terms are smaller by $O(ϵ)$ in the equations retained here. Linearizing these equations recovers ${\omega }^2=4{\Omega }^2{k}_{\parallel }^2/k_{\perp }^2$, the anisotropic limit of part a.

Within the reduced equations for an anisotropic rotating fluid, a flow that remains purely horizontal has $u_{\parallel }=0$ and ${\partial }_t{u}_{\parallel }=0$. The second reduced equation therefore requiresThus the retained horizontal motion is two-dimensional. This is a Taylor–Proudman theorem consequence in the slow balanced ordering; setting $u_{\parallel }=0$ at only one instant would not suffice, because vertical motion could subsequently develop.

The word “exactly” applies to the reduced model. It is not a theorem for every finite-Rossby number solution satisfying only $\Omega \gg kU$. For example, the full rotating Euler equation admitsIt is an incompressible flow, has zero nonlinear advection, and obeys ${\partial }_t\mathbf{u}=2\mathbf{u}\times \mathbf{\Omega }$. It is purely horizontal but depends on $z$, and $\Omega \gg kA$ can hold. This fast inertial oscillation lies outside the slow anisotropic critical balance used in part c; it shows why that ordering cannot be dropped from the conclusion.

For a uniform axial magnetic field, the anisotropic Shear Alfvén wave is predominantly horizontally polarized, whereas the pseudo-Alfvén wave is predominantly vertically polarized. Without rotation both have characteristic frequency ${\omega }_A\sim |k_{\parallel }|v_A$. Pressure projection makes the effective Coriolis acceleration couple these two transverse wave polarization states at the inertial-wave rateThus the useful comparison is not $\Omega$ alone, nor anisotropy alone, butIn this limit rotational coupling of shear and pseudo-Alfvén waves is negligible on an Alfvén propagation time, so the two familiar linear polarizations are decoupled to leading order. The factor $k_{\parallel }$ cancels: simply making both disturbances more anisotropic does not by itself remove their relative rotational coupling.

For a nonlinear reduction one must also make the rotational mixing slow compared with the nonlinear evolution being retained. Under Alfvénic critical balance, $|k_{\parallel }|v_A\sim k_{\perp }U$, the same inequality gives ${\omega }_I\ll k_{\perp }U$. Rapid rotation relative to the horizontal nonlinear rate can therefore coexist with leading Alfvénic decoupling in the windowHigh electrical conductivity makes magnetic diffusion slow, and the corresponding inviscid approximation is also assumed. If $2\Omega \gtrsim k_{\perp }{v}_A$, rotational and magnetic restoring effects are comparable and the modes are mixed magneto-Coriolis disturbances; in the rotation-dominated limit they separate into fast inertial-like and slow magnetostrophic branches, rather than into independently evolving shear and pseudo-Alfvén polarizations.

At finite amplitude, independence is generally one-way in the usual anisotropic ideal magnetohydrodynamics reduction: the Alfvénic subsystem has no leading feedback from pseudo-Alfvénic perturbations, but it can still advect them. The small rotational-coupling condition removes an additional two-way coupling, rather than all nonlinear advection.

This is an asymptotic finite-timescale statement. The two nonrotating Alfvén polarizations are degenerate, so even a weak nonzero rotational coupling can rotate their superposition over the much longer time ${\omega }_I^{-1}$. Exact decoupling for arbitrarily long times does not follow from a small ratio alone. The strictly $k_{\parallel }=0$ nonpropagating limit is separate.