Precession of Fluid Ring Dynamics

• Precession of a fluid ring is the dynamic behavior of an annular fluid subjected to simultaneous rotation and precessional forcing.
• It involves interplay among inertial forces, Coriolis effects, viscous dissipation, and boundary-induced mode selection that governs inertial modes.
• The phenomenon spans laminar, resonant, and turbulent regimes and finds applications in geophysical, astrophysical, and quantum hydrodynamic systems.

A precessing fluid ring refers to the distinctive dynamics arising when a fluid in an annular (ring-shaped) geometry is subjected to simultaneous rotation about its own symmetry axis and precession about a secondary axis. This phenomenon encompasses a range of classical and quantum hydrodynamic regimes and is central to studies in geophysical fluid dynamics, astrophysical flows, thin-film rotary systems, and nonequilibrium condensates. The precession of a fluid ring is governed by the interplay between inertial forces, Coriolis effects, viscous dissipation, spin–orbit coupling (where relevant), and boundary-induced mode selection, resulting in a rich spectrum of laminar, resonant, and turbulent flow structures.

1. Governing Equations and Dimensionless Parameters

The precession-driven flow in an annular geometry is typically described by the incompressible Navier–Stokes equations, including both Coriolis and precessional/Poincaré forcing. Considering a rotating cylindrical annulus with angular speed $\Omega_0\hat{z}$ and precessional angular velocity $\bm{\Omega}_p$, the dimensional governing equations are: $$\frac{\partial\bm u}{\partial t} + \bm u\cdot\nabla\bm u + 2\,\bm\Omega_0\times\bm u = -\frac{1}{\rho}\,\nabla p + \nu\nabla^2\bm u + \bm F_{\rm P}, \qquad \nabla\cdot\bm u=0$$ with the Poincaré force,

$$\bm F_{\rm P} = -2(\bm\Omega_p\times\bm\Omega_0)\times\bm r = -2\Omega_0\Omega_p\,s\cos(\phi+\Omega_0 t)\,\hat{z}$$

After nondimensionalization (using height $h$, velocity $U_0 = h\Omega_0$, time $T=1/\Omega_0$), the key control parameters are:

• Ekman number: $E = \nu/(\Omega_0 h^2)$ — ratio of viscous to Coriolis forces
• Poincaré number: $Po = |\Omega_p|/\Omega_0$ — relative strength of precessional forcing
• Outer radius–height ratio: $\Gamma = r_o/h$
• Inner radius–height ratio: $\bm{\Omega}_p$0 These parameters critically determine the mode structure, stability, and transition thresholds of the precessing flow (Liu et al., 2017, Lin et al., 2014).

2. Inertial Modes, Forcing, and Dominant Dynamics

At small $\bm{\Omega}_p$1, the velocity field is naturally decomposed into a complete orthonormal basis of inertial modes. Each mode $\bm{\Omega}_p$2 solves

$\bm{\Omega}_p$3

with boundary conditions and Bessel-function constraints reflecting the annular geometry. The eigenfrequencies $\bm{\Omega}_p$4 are constrained by $\bm{\Omega}_p$5. The lowest-order mode, typically $\bm{\Omega}_p$6, is the most efficiently excited by the Poincaré forcing (azimuthal wavenumber $\bm{\Omega}_p$7 and frequency $\bm{\Omega}_p$8), as its frequency exhibits the smallest detuning from the forcing. Explicitly, for a representative annulus $\bm{\Omega}_p$9: $$\frac{\partial\bm u}{\partial t} + \bm u\cdot\nabla\bm u + 2\,\bm\Omega_0\times\bm u = -\frac{1}{\rho}\,\nabla p + \nu\nabla^2\bm u + \bm F_{\rm P}, \qquad \nabla\cdot\bm u=0$$0 The forced amplitude hierarchy is $$\frac{\partial\bm u}{\partial t} + \bm u\cdot\nabla\bm u + 2\,\bm\Omega_0\times\bm u = -\frac{1}{\rho}\,\nabla p + \nu\nabla^2\bm u + \bm F_{\rm P}, \qquad \nabla\cdot\bm u=0$$1 in $$\frac{\partial\bm u}{\partial t} + \bm u\cdot\nabla\bm u + 2\,\bm\Omega_0\times\bm u = -\frac{1}{\rho}\,\nabla p + \nu\nabla^2\bm u + \bm F_{\rm P}, \qquad \nabla\cdot\bm u=0$$2 (Liu et al., 2017). In wider parameter studies and experiments, this dominant forced mode structure persists until critical thresholds where nonlinear or parametric excitation occurs (Lin et al., 2014).

3. Regimes of Precessional Dynamics and Transition Thresholds

The response of a precessing fluid ring encompasses three principal dynamical regimes:

• Laminar (steady) inertial-mode flow: $$\frac{\partial\bm u}{\partial t} + \bm u\cdot\nabla\bm u + 2\,\bm\Omega_0\times\bm u = -\frac{1}{\rho}\,\nabla p + \nu\nabla^2\bm u + \bm F_{\rm P}, \qquad \nabla\cdot\bm u=0$$3. Dominated by a single forced inertial mode, $$\frac{\partial\bm u}{\partial t} + \bm u\cdot\nabla\bm u + 2\,\bm\Omega_0\times\bm u = -\frac{1}{\rho}\,\nabla p + \nu\nabla^2\bm u + \bm F_{\rm P}, \qquad \nabla\cdot\bm u=0$$4.
• Triadic-resonance–induced instability: In full (non-annular) cylinders or thin annular geometries, for $$\frac{\partial\bm u}{\partial t} + \bm u\cdot\nabla\bm u + 2\,\bm\Omega_0\times\bm u = -\frac{1}{\rho}\,\nabla p + \nu\nabla^2\bm u + \bm F_{\rm P}, \qquad \nabla\cdot\bm u=0$$5, the forced mode can parametrically excite two free inertial modes via triadic resonance, satisfying $$\frac{\partial\bm u}{\partial t} + \bm u\cdot\nabla\bm u + 2\,\bm\Omega_0\times\bm u = -\frac{1}{\rho}\,\nabla p + \nu\nabla^2\bm u + \bm F_{\rm P}, \qquad \nabla\cdot\bm u=0$$6 (Lin et al., 2014). The free modes grow exponentially until rapid collapse transfers energy to small-scale inertial waves.
• Nonlinear, turbulent, and harmonic-rich flow: At $$\frac{\partial\bm u}{\partial t} + \bm u\cdot\nabla\bm u + 2\,\bm\Omega_0\times\bm u = -\frac{1}{\rho}\,\nabla p + \nu\nabla^2\bm u + \bm F_{\rm P}, \qquad \nabla\cdot\bm u=0$$7, triadic resonance ceases; dynamics become dominated by harmonics ($$\frac{\partial\bm u}{\partial t} + \bm u\cdot\nabla\bm u + 2\,\bm\Omega_0\times\bm u = -\frac{1}{\rho}\,\nabla p + \nu\nabla^2\bm u + \bm F_{\rm P}, \qquad \nabla\cdot\bm u=0$$8) and a broad, turbulent state emerges.

In annuli with significant inner boundary, parametric triadic resonance is suppressed — the inner wall enhances viscous damping, restricting the amplitude growth of non-resonant and resonant modes. As a result, the transitional threshold for instability is elevated: laminar flow persists until $$\frac{\partial\bm u}{\partial t} + \bm u\cdot\nabla\bm u + 2\,\bm\Omega_0\times\bm u = -\frac{1}{\rho}\,\nabla p + \nu\nabla^2\bm u + \bm F_{\rm P}, \qquad \nabla\cdot\bm u=0$$9 ($$\bm F_{\rm P} = -2(\bm\Omega_p\times\bm\Omega_0)\times\bm r = -2\Omega_0\Omega_p\,s\cos(\phi+\Omega_0 t)\,\hat{z}$$0 at $$\bm F_{\rm P} = -2(\bm\Omega_p\times\bm\Omega_0)\times\bm r = -2\Omega_0\Omega_p\,s\cos(\phi+\Omega_0 t)\,\hat{z}$$1), and destabilization is driven by boundary-layer nonlinear interactions rather than triadic resonances (Liu et al., 2017).

4. Analytical Theories for Thin Film and Rotating Cylinder Cases

In the regime of a thin fluid film coating a rotating cylinder, the precession effect is captured via amplitude-reduced models analogous to shallow water equations. A two-dimensional “vortex” shallow-water system on the cylinder, after depth-averaging and multiple scales reduction, yields: $$\bm F_{\rm P} = -2(\bm\Omega_p\times\bm\Omega_0)\times\bm r = -2\Omega_0\Omega_p\,s\cos(\phi+\Omega_0 t)\,\hat{z}$$2 where $$\bm F_{\rm P} = -2(\bm\Omega_p\times\bm\Omega_0)\times\bm r = -2\Omega_0\Omega_p\,s\cos(\phi+\Omega_0 t)\,\hat{z}$$3 is the film-thickness perturbation, and $$\bm F_{\rm P} = -2(\bm\Omega_p\times\bm\Omega_0)\times\bm r = -2\Omega_0\Omega_p\,s\cos(\phi+\Omega_0 t)\,\hat{z}$$4 is the precessional drift speed. The precession frequency of the ring is explicitly given by

$$\bm F_{\rm P} = -2(\bm\Omega_p\times\bm\Omega_0)\times\bm r = -2\Omega_0\Omega_p\,s\cos(\phi+\Omega_0 t)\,\hat{z}$$5

with $$\bm F_{\rm P} = -2(\bm\Omega_p\times\bm\Omega_0)\times\bm r = -2\Omega_0\Omega_p\,s\cos(\phi+\Omega_0 t)\,\hat{z}$$6 the mean film radius, $$\bm F_{\rm P} = -2(\bm\Omega_p\times\bm\Omega_0)\times\bm r = -2\Omega_0\Omega_p\,s\cos(\phi+\Omega_0 t)\,\hat{z}$$7 the cylinder radius, and for $$\bm F_{\rm P} = -2(\bm\Omega_p\times\bm\Omega_0)\times\bm r = -2\Omega_0\Omega_p\,s\cos(\phi+\Omega_0 t)\,\hat{z}$$8 (i.e., sufficiently wide ring). Both solitary wave and cnoidal wave solutions are accessible, and the solitary pulses themselves precess azimuthally at $$\bm F_{\rm P} = -2(\bm\Omega_p\times\bm\Omega_0)\times\bm r = -2\Omega_0\Omega_p\,s\cos(\phi+\Omega_0 t)\,\hat{z}$$9, distinct from the rotation rate $h$0 of the cylinder (Zhukov et al., 2013). Painlevé analysis affirms complete integrability in this reduced regime.

5. Precession in Ballistic Ring Condensates with Spin-Orbit Coupling

In quantum hydrodynamics, notably for ballistic polariton condensates in ring traps, precessional phenomena acquire additional structure via spin–orbit coupling. The polariton spinor wavefunction

$h$1

evolves under a Pauli-matrix Hamiltonian with a TE–TM–derived effective Zeeman (spin–orbit interaction) field $h$2. For a ballistic wavepacket with group velocity $h$3,

$h$4

where $h$5 is the TE–TM splitting constant and $h$6 the polariton mass. The experimental signatures include (i) sustained precession of the polarization vector around the ring and (ii) a “zitterbewegung” oscillatory trajectory arising from the non-commutativity of the SOI with the position operator. These effects are quantitatively reproduced by spinor Gross–Pitaevskii models and have been directly observed in microcavity ring lattices (Yao et al., 2022).

6. Boundary Effects and Suppression of Resonant Instabilities

The geometry and boundary conditions of the fluid ring fundamentally alter the eigenmode structure and instability mechanisms. The presence of an inner wall ($h$7) modifies the radial eigenvalue problem for inertial modes, imposing an annular Bessel-matching condition that shifts mode frequencies and amplifies viscous damping. This boundary-induced damping inhibits the growth of all but the least-detuned forced mode, effectively preventing triadic resonance instabilities over a broad parameter range (Liu et al., 2017). Experimental studies confirm that as the relative thickness increases (i.e., as $h$8), this suppression effect is enhanced (Lin et al., 2014).

Parameter	Definition	Physical Role
$h$9	$U_0 = h\Omega_0$0	Viscous/Coriolis ratio
$U_0 = h\Omega_0$1	$U_0 = h\Omega_0$2	Precessional forcing strength
$U_0 = h\Omega_0$3	$U_0 = h\Omega_0$4	Outer radius–height aspect
$U_0 = h\Omega_0$5	$U_0 = h\Omega_0$6	Inner radius–height aspect

7. Summary and Physical Interpretation

The precession of a fluid ring yields a prototypical example of rotationally constrained hydrodynamics modulated by external symmetry-breaking forcing. In the classical annulus, the flow transitions from mode-dominated laminar states to turbulent, harmonic-rich states via thresholds set by $U_0 = h\Omega_0$7 and $U_0 = h\Omega_0$8, with geometry-driven suppression or allowance of triadic parametric instability. In the thin-film limit, precession arises naturally from the dispersion relation of KdV-like equations, with explicit analytical expressions for the precessional frequency. In quantum ring systems, precession re-emerges as spinor precession due to SOI, with both theoretical and experimental confirmation of frequency and trajectory modulation. Collectively, the precession of a fluid ring thus exemplifies the deep coupling between geometry, boundary effects, forcing, and underlying hydrodynamic or quantum structure (Liu et al., 2017, Lin et al., 2014, Zhukov et al., 2013, Yao et al., 2022).