Second Kelvin–Arnol'd Stability Theorem

• The Second Kelvin–Arnol'd Stability Theorem is a criterion that quantifies nonlinear (Lyapunov) stability of steady ideal fluid flows using energy–Casimir functionals and sharp spectral bounds.
• It extends Arnold's original framework by incorporating functional-analytic bounds and Poincaré-type inequalities to account for complex flow geometries in bounded and multiply-connected domains.
• Applications of the theorem include algebraic stability tests for axisymmetric flows and explicit variational characterizations, providing practical insights into stability boundaries in fluid dynamics.

The Second Kelvin–Arnol’d Stability Theorem provides a quantitative criterion for the nonlinear (Lyapunov) stability of steady flows in two-dimensional ideal fluids under the Euler equations. This result extends Arnold's energy–Casimir stability framework and refines classical criteria by incorporating sharp functional-analytic bounds related to the geometry of the domain and the spectral properties of associated elliptic operators. The theorem is applicable to bounded domains—both simply- and multiply-connected—and yields powerful sufficient conditions that are directly verified by spectral data or Poincaré-type inequalities. Generalizations include axisymmetric flows in annular and pipe geometries, where the theorem admits concise algebraic stability tests.

1. Mathematical Formulation and Main Statement

Let $D \subset \mathbb{R}^2$ be a smooth bounded domain with $N+1$ boundary components $\Gamma_0$ (outer) and $\Gamma_1,\dots,\Gamma_N$ (inner). A steady incompressible ideal flow is characterized by a stream function $\psi^s$ and vorticity $\omega^s$ related through

$-\Delta\psi^s = \omega^s = g(\psi^s)$

with $g \in C^1(\mathbb{R})$, $\psi^s|_{\Gamma_i} = \text{const}$ for each $i$, and $\int_D \psi^s dx = 0$. The key quantity is $g'(\psi^s) = \nabla\omega^s / \nabla\psi^s$ and the central assertion is:

Arnol'd's Second Stability Theorem

Suppose that

$$0 < \min_{D} g'(\psi^s) \leq \max_{D} g'(\psi^s) < C_{ar}$$

with $C_{ar} > 0$. Then the steady flow $v^s = \nabla^\perp\psi^s$ is nonlinearly Lyapunov stable: for any Euler solution with initial data $\psi(0)$ sufficiently close to $\psi^s$ in $W^{2,p}$ norm, $\psi(t)$ remains in an arbitrary $W^{2,p}$ neighborhood of $\psi^s$ for all $t$ (Wang et al., 11 May 2025).

A crucial refinement is that $C_{ar}$ can be chosen as the first eigenvalue $\Lambda_1$ of $-\Delta$ in the space of mean-zero functions with piecewise constant boundary values.

2. Functional-Analytic and Variational Framework

The theorem is cast in a precise functional-analytic setting. The stream function spaces are defined as

$$X = \{u \in H^1(D) : u|_{\Gamma_0} = 0,\, u|_{\Gamma_i} = \text{const},\, i=1\dots N\}$$

$$Y = \{u \in H^2(D) \cap X : \int_D u\,dx = 0,\, \int_{\Gamma_i} \partial_n u\,dS = 0,\, i=1\dots N\}$$

For any $v \in L^p(D)$ and any circulation vector $\gamma \in \mathbb{R}^N$, the elliptic boundary value problem $-\Delta u = v$ admits a unique solution $u \in W^{2,p}(D)$, naturally identifying a bounded, linear operator $P:L^p(D) \rightarrow W^{2,p}(D)$.

The stability condition is reformulated via the energy–Casimir functional: $$E(\omega, \gamma) = \frac{1}{2}\int_D \omega\,P\omega\,dx + \int_D h_\gamma\,\omega\,dx + \frac{1}{2} \gamma^T Q \gamma$$ where $h_\gamma$ encodes boundary circulations and $Q$ is a symmetric matrix on circulation space. The rearrangement class $R_{\omega^s}$ is preserved by the flow, and stability is deduced by showing that $\omega^s$ (or its affine orbit) uniquely maximizes $E$ over $R_{\omega^s}$.

3. Spectral Sharpness and the Role of $\Lambda_1$

The maximal allowable $g'(\psi^s)$ is characterized via eigenvalues of a compact, self-adjoint operator $T$ derived from $P$. Specifically,

$T v = P v - |D|^{-1}\int_D (P v) dx$

acting on $L_0^2(D)$ (mean-zero functions). The Laplacian eigenproblem $-\Delta u = \Lambda u$ on $Y$ corresponds to $v = \Lambda T v$ in $L_0^2(D)$, yielding eigenvalues $0 < \Lambda_1 < \Lambda_2 < \cdots$.

The sharpness is exhibited in the unit disk by constructing a steady flow with $g(s) = \Lambda_1 s$ (where $\Lambda_1 = j_{1,1}^2$, $j_{1,1}$ being the first zero of $J_1$). Instability is seen for non-circular flows attaining $g'(\psi^s) = \Lambda_1$, while a continuum of maximizers persists within the same rearrangement class (Wang et al., 11 May 2025).

4. Generalizations to Multiply-Connected and Axisymmetric Domains

Recent extensions establish the theorem for multiply-connected planar domains $D = D_0 \setminus \cup_{i=1}^N D_i$ with independent circulation invariants. The criteria are relaxed from pointwise bounds on $g'$ to the quadratic form condition

$$\min_{D} g'(\psi) \geq 0 \quad \text{and} \quad \forall u \in Y:\ \int_D(|\nabla u|^2 - g'(\psi) u^2) dx \geq 0$$

which is stable under weak convergence (Wang et al., 2022). Perturbations need not preserve all individual circulations, in contrast to the original theorem.

For axisymmetric flows in annuli or pipes, the theorem is recast as an explicit algebraic criterion. Given base flow $U(r)$, azimuthal wavenumber $n$, and axial wavenumber $k$, the stability condition reads: $0 \leq W_{\beta,N}(r) \leq H(r) \quad \forall r$ where $W_{\beta,N}(r) = r\,Q'(r)/(U(r)-\beta)$ and $H(r)$ is a geometry-dependent constant (Poincaré-type). The KA-II criterion strictly enlarges the stable region compared to the Batchelor–Gill result (KA-I), and is shown to tightly bracket numerically computed neutral surfaces (Deguchi et al., 25 Jan 2026).

5. Energy–Casimir Variational Methodology

The proof exploits the energy–Casimir functional and its maximization on rearrangement classes. For

$$C'(\omega) = -\psi^s \qquad C''(\omega) = -(\psi^s)'(f(\omega))^{-1}$$

the second variation for an admissible perturbation, represented by a stream function $\varphi$ with homogeneous boundary conditions, yields

$$\delta^2 J[\omega^s](\varphi) = \int_D \{|\nabla\varphi|^2 - g'(\psi^s)\varphi^2\} dx$$

Stability follows from coercivity (strict positivity) of $\delta^2 J$ outside neutral symmetry directions (translations, rotations). In the critical case ($\max g'(\psi^s) = \Lambda_1$), the energy gap vanishes on the lowest eigenspace, leading to orbital rather than strict Lyapunov stability.

For multiply-connected and axisymmetric contexts, supporting functionals and auxiliary quadratic forms are introduced to account for harmonic and circulation components (Wang et al., 2022).

6. Applications, Examples, and Stability Boundaries

The theorem yields practical stability maps for model shear flows:

• Sliding Couette flow in annuli: The KA-II criterion extends the stable region beyond classical Rayleigh and Batchelor–Gill bounds, allowing for positive bounded $g'(\psi^s)$ up to $H(r)$ (Deguchi et al., 25 Jan 2026).
• Hagen–Poiseuille pipe flow with heating: The KA-II bound on parameter $\chi$ is sharper, predicting stability for $\chi > 13.18$, whereas KA-I yields $\chi < 4$. Numerical eigenvalue computations confirm KA-II as an effective predictor of stability boundaries.

The Oseen vortex (Gaussian) is shown to satisfy all conditions, both inviscidly and at low viscosity, with explicit weight functions providing uniform coercivity (Gallay et al., 2021).

Domain Type	Stability Condition for $g'(\psi^s)$	Maximizer Type
Simply-connected	$$0 < \min g'(\psi^s) \leq \max g'(\psi^s) < \Lambda_1$$	Unique in rearrangement class
Multiply-connected	$\forall u: \int(|\nabla u|^2 - g' u^2) \geq 0$	Unique up to circulation constraints
Axisymmetric annulus	$0 \leq W_{\beta,N}(r) \leq H(r)$	Energy–Casimir, explicit algebraic

7. Structural Stability and Orbital Dynamics

In cases where the upper bound $\Lambda_1$ is attained, classical uniqueness fails, but structural (orbital) stability persists. The set of maximizers becomes an affine space $\omega^s + E_1$, and solutions may orbit within this set under conserved quantities like the moment of inertia. In disks, imposing $I(\omega) = I(\omega^s)$ restricts the dynamics to rigid rotations within $E_1$ (Wang et al., 11 May 2025).

This suggests that the Second Kelvin–Arnol’d Stability Theorem provides a robust framework for capturing both strict Lyapunov and more general structural (orbital) stability of ideal fluid flows, with spectral bounds and rearrangement constraints forming the core control mechanisms.