Perturbation analysis of triadic resonance in columnar vortices: selection rules and the roles of external forcing and critical layers

TL;DR

This paper addresses the intrinsic stability of isolated columnar vortices and whether weakly nonlinear triadic resonance can drive intrinsic instability in the inviscid limit.A multi-scale perturbation framework is developed to derive general amplitude equations for triads, proving that conservative interactions among smooth neutral modes obey Manley-Rowe relations and are constrained by a pseudoenergy criterion in a large-$k$ WKBJ setting.The key findings show that explosive instability is forbidden for smooth neutral modes, while instability can arise via two symmetry-breaking pathways: parametric instability driven by external forcing and non-Hermitian active critical layers that extract energy from the background shear.These results provide theoretical guidance for wake-vortex control, suggesting that mitigation requires either tuned external forcing or engineered critical layers, and they set the stage for future work including viscous regularisation and non-conservative resonances.

Abstract

The remarkable robustness of isolated columnar vortices suggests the existence of fundamental constraints that prevent spontaneous disintegration. In this work, we investigate the weakly nonlinear stability of such flows, demonstrating that the triadic resonance of wave modes is governed by a set of hydrodynamic ``selection rules''. By employing a multi-scale perturbation analysis, we prove that resonant interactions between smooth neutral modes, specifically regular Kelvin waves and discrete critical layer modes with passive singularities, are strictly conservative and confined to the Manley--Rowe relations. Using wave pseudoenergy within a large-$k$ WKBJ framework, we show that these rules topologically prohibit intrinsic instability, analogous to the forbidden transitions of quantum mechanics. Consequently, the breakdown of a columnar vortex requires a specific symmetry-breaking mechanism to overcome this barrier. We identify and analyse two distinct pathways for this violation: (1) \textit{Parametric instability}, where external forcing acts as an active energy pump; using a robust tuning method based on non-degenerate perturbation theory, we generalize classical elliptical instability to arbitrary driving frequencies and identify new instability configurations involving dicrete critical layer modes. (2) \textit{Active critical layers}, where an embedded singularity breaks the Hermitian symmetry of the operator, enabling the extraction of mean-flow energy via a wave-mean resonance. These findings provide a theoretical guidance for flow control, suggesting that the mitigation of aircraft wake vortices requires either tuned external forcing or the ``engineering'' of critical layers (e.g., via thermal stratification) to trigger the forbidden transitions.

Figures (8)

• Figure 1: Stability mechanism via mode degeneracy in the inviscid limit. In the absence of degeneracy, discrete neutral modes (represented by the black circles) are topologically constrained to the imaginary axis (neutral stability) and are only permitted to move vertically along it as parameters vary. When a pair of modes becomes degenerate or nearly degenerate (represented by the white concentric circles), a resonant interaction can unfold in two distinct ways: either via detuning, where the modes separate along the imaginary axis, or via instability, where the modes are pulled symmetrically off the axis into the complex plane, resulting into a conjugate pair of growing and decaying modes.
• Figure 2: Example numerical solutions of the normalised conservative three-wave amplitude equations (\ref{['eqn:3-wave-conservative-universal']}) for distinct interaction coefficient sign signatures ($s_0$, $s_1$, $s_2$). (a--c) Temporal evolution of the normalised wave intensities $|\hat{A}_j|^2$. (d--f) Corresponding 3D phase space trajectories plotting $|\hat{A}_{1,2}|$ against $|\hat{A}_0|$ and the relative phase factor $\cos{(\phi_2-\phi_1-\phi_0)}$. White and black markers indicate the initial and final states of the evolution, respectively. Left and Right Columns: Bounded resonance regimes, where the system is strictly confined to closed periodic orbits on a compact manifold. Middle Column: The explosive resonance regime ($s_0=s_1=-s_2$), where the relative phase locks to a constant value while permitting the system to follow an open trajectory with simultaneous amplitude divergence. See also craik_1986.
• Figure 3: Effect of detuning $\Delta\omega$ on parametric instability as given in (\ref{['eqn:sigma-detuned']}). (a) Doppler-shifted eigenfrequencies $\mathrm{Im}(\sigma_{\text{par},j}) - k_j\omega_0/k_0$ for modes $j=1,2$. The solid and dashed lines represent the upper and lower branches of the four coupled eigenvalues, respectively. Inside the instability band (shaded region), the upper and lower branches of the same mode locks to the same wave frequency, and the eigenvalues merge toward exact resonance where the modes synchronize to a common Doppler-shifted frequency. (b) Real growth rates $\mathrm{Re}(\sigma)$. In the region where frequencies are locked, the growth rates bifurcate to form a symmetric instability bubble (shaded region indicating positive growth). The two wave modes share identical growth rates, leading to the appearance of a single curve for the unstable branch. Outside this region, the frequency mismatch becomes significant, and the coupled modes are no longer in near resonance.
• Figure 4: Characteristics of the inviscid linear modes of the Lamb--Oseen vortex obtained numerically by MLEGS MLEGS. (a) Dispersion curves $\omega(k)$ of helical modes ($m=1$). Solid lines denote Kelvin waves, and dashed lines indicate neutral discrete modes possessing critical layers. The shaded frequency band $\omega_\text{CL} = [-m,0)$ is associated with the critical layer singularity and corresponds to the continuous spectrum. (b--d) Radial velocity component $r\tilde{u}_r(r)$ of selected eigenvectors ($m=1$, $k=3.0$): (b) A regular discrete Kelvin wave, exhibiting a purely real, smooth profile characteristic of global neutral modes; (c) A generic singular critical layer mode from the continuous spectrum. The singularity at the critical radius $r_c$ (vertical dashed line) breaks analytic continuity, resulting in spectral broadening and a distinct phase jump in physical space; (d) A "passive" discrete critical layer mode. Unlike generic critical layer modes, this solution maintains high-order spectral convergence, confirming that the critical layer singularity is effectively suppressed and decoupled from the core dynamics. The eigenvalues associated with the mode profiles shown are marked as circles in the dispersion plot.
• Figure 5: Graphical representations of resonant triads in the Lamb--Oseen vortex. (a) Resonance sustained by a stationary pumping wave with zero axial wavenumber ($m_p, k_p=0,\omega_p=0$). Star symbols mark crossing points between distinct dispersion branches ($k_i = k_j$, $\omega_i = \omega_j$) in resonance with $|m_p| = |m_i - m_j|$; filled circles indicate intersections with the zero-frequency axis, corresponding to a pair of conjugate modes ($\pm m$) in resonance with a static strain field with symmetry $m_p=2m$. (b) General triadic resonance satisfying the frequency matching condition $\omega_1 + \omega_2 = \omega_3$ alongside wavenumber matching. The specific examples represent resonance between modes $(m,k,\omega)$: $(1, 5.6, -1.2) + (2, 0.50, -2.0) = (3, 6.0, -3.3)$ (dash-dotted lines with star markers), and $(1, 4.1, -1.2) + (2, 0.98, -2.1) = (3, 5.1, -3.3)$ (dashed lines with circular markers).

Theorems & Definitions (4)

• Corollary 1
• proof
• Corollary 2
• proof