Stability of the Rankine Vortex and Perimeter Growth in Vortex Patches

Basic Information

• Paper ID: 2511.09772
• Title: Stability of the Rankine Vortex and Perimeter Growth in Vortex Patches
• Author: John Brownfield
• Classification: math.AP (Analysis of PDEs)
• Publication Date: November 2025
• Paper Link: https://arxiv.org/abs/2511.09772

Abstract

This paper establishes that for vorticity functions $\omega: \mathbb{R}^2 \to [0,1]$ with the same total circulation and vorticity center as the Rankine vortex, the $L^1$ deviation from the Rankine patch can be bounded by the pseudo-energy deviation and angular momentum. In the case of $m$-fold symmetry, the dependence on angular momentum can be removed. Using this result, the author confirms previous numerical simulations, demonstrating that simply-connected perturbations of the Rankine vortex exhibit perimeter growth linear in time.

Research Background and Motivation

Problem Background

This paper studies the vorticity formulation of the two-dimensional incompressible Euler equations: $\partial_t\omega + u \cdot \nabla\omega = 0$ where the velocity field $u$ is given by the Biot-Savart law: $u = K * \omega$, $K(x) = \frac{x^\perp}{2\pi|x|^2}$.

Core Problems

1. Loss of Regularity Problem: Although Chemin 1 proved that smooth initial boundaries remain smooth, recent research has focused on constructing solutions that lose initial regularity over infinite time. These solutions are typically stable perturbations of known simple solutions (in weak norm sense).
2. Stability of the Rankine Vortex: The Rankine vortex patch $\omega^* = 1_{B(a,r)}$ is a steady state because the disk is the unique minimizer of angular momentum for fixed total circulation and vorticity center. Previous work (Pulvirenti-Wan 8, Sideris-Vega 6, Tang 7) established quantitative $L^1$ stability estimates, but with limitations.
3. Perimeter Growth: Dritschel's 2 numerical simulations show that vortex patches may exhibit sustained perimeter growth. Drivas-Elgindi-Jeong 3 first constructed multiply-connected vortex patches with linear-time perimeter growth, but the simply-connected case remained unresolved.

Research Motivation

• Improve existing stability estimates: remove initial $L^1$ distance restrictions, relax patch assumptions to $0 \leq \omega \leq 1$
• In the $m$-fold symmetric case, remove dependence on angular momentum and improve exponents
• Construct the first simply-connected vortex patch with linear-time perimeter growth, filling a theoretical gap

Core Contributions

1. Improved Stability Theorems (Theorems 1.2 and 1.3):
   • Under $m$-fold symmetry, prove $\|\omega_t - 1_{B(0,r)}\|_1^2 \leq C(E(1_{B(0,r)}) - E(\omega_0))$, removing angular momentum dependence
   • Under bounded angular momentum, prove $\|\omega_t - 1_{B(0,r)}\|_1^4 \leq C(I,r)(E(1_{B(0,r)}) - E(\omega_0))$
2. Perimeter Growth for Simply-Connected Vortex Patches (Theorem 1.4):
   • First construction of a simply-connected open set $\Omega_0$ whose vortex patch solution satisfies $\text{Perimeter}(\Omega_t) \gtrsim t$
   • Confirms theoretical predictions from Dritschel's numerical simulations
3. Technical Innovations:
   • Utilize rearrangement inequality results from Frank-Lieb 5 (Proposition 1.1)
   • Develop topological argument methods for handling the simply-connected case ("bucket" decomposition technique)

Detailed Methodology

Task Definition

Input: Initial vorticity distribution $\omega_0: \mathbb{R}^2 \to [0,1]$ satisfying:

• Total circulation: $\int \omega_0 = \pi r^2$
• Vorticity center: $\int x\omega_0 = 0$
• Symmetry or angular momentum constraints

Output:

1. $L^1$ stability estimate: upper bound for $\|\omega_t - 1_{B(0,r)}\|_1$
2. Perimeter growth rate: lower bound for $\text{Perimeter}(\Omega_t)$

Stability Analysis Framework

Step One: Rearrangement Inequality (Proposition 1.1)

Using Frank-Lieb's result, for $\omega: \mathbb{R}^n \to [0,1]$: $\int\int_{E^*\times E^*} \ln(|x-y|^{-1})dxdy - \int\int_{\mathbb{R}^n\times\mathbb{R}^n} \omega(x)\ln(|x-y|^{-1})\omega(y)dxdy \geq C_n \inf_{a\in\mathbb{R}^n} \|\omega - 1_{E^*+a}\|_{L^1}^2$

where $E^*$ is the ball centered at the origin with measure $\int \omega$.

Step Two: Geometric Inequalities (Propositions 3.1 and 3.2)

Proposition 3.1 ($m$-fold Symmetry): If $\omega = \omega \circ R_{2\pi/m}$ ($m \geq 2$), then $\inf_{a\in\mathbb{R}^2} \|\omega - 1_{E^*+a}\|_{L^1} \geq \frac{1}{3}\|\omega - 1_{E^*}\|_{L^1}$

Proof Strategy:

• For even $m$, use $\|\omega - 1_{E^*+a}\|_{L^1} = \|\omega - 1_{E^*-a}\|_{L^1}$
• Examine $B(a,r) \cap B(-a,r) \subset B(0,r)$, obtain $\epsilon \geq \delta/2$ through measure calculations
• For odd $m = 2n+1$, construct three balls $B(a_1,r), B(a_2,r), B(a_3,r)$ whose triple intersection is contained in $B(0,r)$
• Apply inclusion-exclusion principle to get $\epsilon \geq \delta/3$

Proposition 3.2 (Bounded Angular Momentum): If $\int |x|^2\omega \leq I$, $\int x\omega = 0$, then $\inf_{a\in\mathbb{R}^2} \|\omega - 1_{B(a,r)}\|_{L^1} \geq C(I,r)\|\omega - 1_{B(0,r)}\|_{L^1}^2$

and the exponent 2 is optimal.

Step Three: Stability Theorem Proof

Combining Propositions 1.1 and 3.1/3.2, using pseudo-energy conservation: $\|\omega_t - 1_{B(0,r)}\|_1^2 \leq C \inf_{a} \|\omega_t - 1_{B(a,r)}\|_1^2 \leq C'(E(1_{B(0,r)}) - E(\omega_t)) = C'(E(1_{B(0,r)}) - E(\omega_0))$

Perimeter Growth Construction

Initial Setup (Section 4.1)

Construct a 3-fold symmetric simply-connected set $\Omega_0$:

• Core: disk of radius approximately $1-\epsilon$, i.e., $B(0, 1-\epsilon)$
• Three slender arms: length $N \gg 1$, total area $\gamma \ll 1$
• Parameter relations: $\gamma N^2 \sim R^3$ (angular momentum), $\gamma \ln(N) \sim \delta$ (pseudo-energy deviation)

Key properties:

1. Total circulation: $\int \omega_0 = \pi$
2. Angular momentum: $\int |x|^2\omega_0 \sim 1 + \gamma N^2$ (large angular momentum ensures $\Omega_t \backslash B(0,R) \neq \emptyset$)
3. Pseudo-energy: $E(1_{B(0,1)}) - E(\omega_0) \lesssim \delta \ll 1$ (ensures $L^1$ stability)

Velocity Field Estimates (Lemma 4.1)

Let $u = u_r e_r + ru_\theta e_\theta$ be the velocity field induced by $\omega_t$, and $u^* = r\mu(r)e_\theta$ be the velocity field induced by $1_{B(0,1)}$, where $\mu(r) = \begin{cases} \frac{1}{2} & r \leq 1 \\ \frac{1}{2r} & r \geq 1 \end{cases}$

Then we have the estimates:

1. $\|u - u^*\|_{L^\infty} \lesssim \tilde{\epsilon}^{1/2}$
2. $\|u_\theta\|_{L^\infty} \lesssim 1$
3. $\|u_\theta - \mu\|_{L^2} \lesssim \ln(R)\tilde{\epsilon}^{1/4}$
4. $\|\frac{1}{r}u_r\|_{L^2} \lesssim \ln(R)\tilde{\epsilon}^{1/4}$

Topological Argument (Section 4.2)

Universal Cover Lifting:

• Lift the flow map $\Phi$ to the universal cover $\mathbb{R}\times\mathbb{R}^+$ of $\mathbb{R}^2\backslash\{0\}$
• The lifted map $\tilde{\Phi} = (\Phi_r, \Phi_\theta)$ satisfies: $\|\Phi_\theta - \theta - t\mu(\Phi_r)\|_{L^2} \lesssim \epsilon \cdot t$

"Bucket" Decomposition Technique:

• Define $M_t^n$ as the first (left to right) radial maximum point of $\partial\tilde{\Omega}_t^n$
• Draw vertical lines from each $M_t^n$, decomposing $\Omega_t^c$ into disjoint "buckets" $B_t^n$
• For a point $z_0 \in \Omega_0^c$, define $N(t)$ as the bucket label containing $Z_t = \tilde{\Phi}(z_0, t)$

Key Lemma 4.2: $|N(T)| \lesssim (1 + \max(1/R, \epsilon)T) \lesssim (\nu + o(1))T$ where $\nu = \max(1/R, \epsilon) \ll 1$.

Proof Strategy:

• Choose time step $\Delta \sim \min(R, \epsilon^{-1})$ so that angular change at radius $\geq R$ is $\leq 1$
• Prove that the evolution of bucket $B_t^n$ in time $\Delta$ is confined to $B_{t+\Delta}^{n-2} \cup \cdots \cup B_{t+\Delta}^{n+2}$
• Partition $[0,T]$ into intervals of length $\Delta$, obtaining $|N(T)| \lesssim 1 + T/\Delta$

Perimeter Growth Proof

1. Find two points $p_0, q_0$ in $\Omega_0^c$ satisfying:
   • $\Phi_\theta(p_0, T) > T/8$, $\Phi_\theta(q_0, T) < T/16$
   • $1 < \Phi_r(p_0, T), \Phi_r(q_0, T) < 20$
2. Apply Lemma 4.2 to find points $p, q$ on $\partial\Omega_0$ satisfying:
   • $\Phi_\theta(q, T) < c_0 T < c_1 T < \Phi_\theta(p, T)$ (angular difference)
   • $\Phi_r(p, T), \Phi_r(q, T) > r_0 > 0$ (away from origin)
3. Use geometric lemma (from 3): the arc of $\partial\Omega_T$ connecting $\Phi(p,T)$ and $\Phi(q,T)$ has length at least $2r_0(c_1-c_0)T - 1 \gtrsim T$

Technical Innovations

1. Geometric-Analytic Synthesis: Transform Frank-Lieb's abstract rearrangement inequality into concrete vorticity stability estimates
2. Symmetry Exploitation: Clever application of $m$-fold symmetry avoids angular momentum dependence, a substantial improvement over Tang 7
3. Topological Methods: The "bucket" decomposition on the universal cover is the key innovation for handling the simply-connected case, solving the problem that the multiply-connected method in 3 could not address
4. Energy-Momentum Balance: Through careful parameter selection ($\gamma, N, R, \delta$), balance large angular momentum (ensuring radial support) with small pseudo-energy deviation (ensuring $L^1$ stability)

Experimental Setup

Theoretical Proof Nature

This is a pure theoretical mathematics paper containing no numerical experiments. All results are obtained through rigorous mathematical proofs.

Construction Parameters

For the construction in Theorem 1.4:

• Symmetry: $m = 3$
• Core radius: $r \approx 1 - \epsilon$
• Arm length: $N \gg 1$
• Arm area: $\gamma \ll 1$
• Parameter relations: $\gamma N^2 \sim R^3$, $\gamma \ln(N) \sim \delta$
• Constraints: $\delta \ll 1$, $R \gg 1$, $\epsilon = \max(\tilde{\epsilon}^{1/2}, \ln(R)\tilde{\epsilon}^{1/4}) \ll 1$

Experimental Results

Main Theoretical Results

Theorem 1.2 ($m$-fold Symmetric Stability):

• There exists a universal constant $C > 0$ such that for $\omega_0: \mathbb{R}^2 \to [0,1]$ satisfying $\omega_0 = \omega_0 \circ R_{2\pi/m}$ ($m \geq 2$) and $\int \omega_0 = \pi r^2$
• The solution satisfies: $\|\omega_t - 1_{B(0,r)}\|_1^2 \leq C(E(1_{B(0,r)}) - E(\omega_0))$ for all time
• Improvement: Compared to Tang 7, removes initial $L^1$ distance restriction, angular momentum dependence, and improves exponent from 4 to 2

Theorem 1.3 (Bounded Angular Momentum Stability):

• For $I, r > 0$, there exists $C = C(I,r) > 0$ such that if $\omega_0$ satisfies $\int |x|^2\omega_0 \leq I$, $\int x\omega_0 = 0$, $\int \omega_0 = \pi r^2$
• The solution satisfies: $\|\omega_t - 1_{B(0,r)}\|_1^4 \leq C(E(1_{B(0,r)}) - E(\omega_0))$
• Improvement: Removes initial $L^1$ distance restriction, relaxes patch assumption

Theorem 1.4 (Linear Perimeter Growth):

• There exists a simply-connected open set $\Omega_0 \subset \mathbb{R}^2$ (smooth boundary) such that the vortex patch solution $\omega_t = 1_{\Omega_t}$ satisfies: $\text{Perimeter}(\Omega_t) \gtrsim t$ for sufficiently large $t$
• Breakthrough: First construction of linear perimeter growth for simply-connected vortex patch

Comparison with Existing Work

Theoretical Findings

1. Power of Symmetry: $m$-fold symmetry ($m \geq 2$) suffices to significantly improve stability estimates, indicating that symmetry plays a fundamental role in vortex dynamics
2. Optimality of Exponent: Proposition 3.2 proves that exponent 2 is optimal in the bounded angular momentum case, revealing the essential difficulty of the problem
3. Topological Constraints: Simply-connectedness restricts the winding speed of points around the origin through the "bucket" decomposition mechanism, which is key to achieving perimeter growth
4. Energy-Geometry Relation: Pseudo-energy deviation directly controls $L^1$ deviation, establishing quantitative connection between energy functional and geometric properties

Related Work

Vortex Patch Stability

1. Pulvirenti-Wan 8 (1985): First used angular momentum and pseudo-energy to establish quantitative $L^1$ stability on bounded domains
2. Sideris-Vega 6 (2009): Extended angular momentum method to the full plane
3. Tang 7 (1987): Used pseudo-energy for full-plane stability, but required initial proximity and bounded angular momentum
4. This paper: Improves exponents, removes restrictions, exploits symmetry

Regularity Theory

1. Yudovich 10 (1963): Proved global unique weak solution existence for $\omega_0 \in L^1 \cap L^\infty$
2. Chemin 1 (1991): Proved smooth initial boundaries remain smooth

Perimeter Growth

1. Dritschel 2 (1988): Numerical simulations show repeated filamentation and perimeter growth of vortex interfaces
2. Drivas-Elgindi-Jeong 3 (2024): First construction of multiply-connected vortex patches with linear perimeter growth, using $L^2$ shear stability estimates
3. This paper: First construction for simply-connected case, develops new topological argument methods

Rearrangement Inequalities

1. Frank-Lieb 5 (2021): Proved spherical clustering theorem based on quantitative rearrangement inequalities
2. Yan-Yao 9 (2022): Pointed out that 5's result implies Proposition 1.1
3. This paper: Applies these abstract results to vortex dynamics

Conclusions and Discussion

Main Conclusions

1. Stability Improvement: Under $m$-fold symmetry, the $L^1$ stability of the Rankine vortex can be controlled solely by pseudo-energy deviation with exponent 2, without angular momentum constraints
2. Simply-Connected Perimeter Growth: First rigorous proof of existence of simply-connected vortex patch with linear-time perimeter growth, confirming Dritschel's numerical simulation predictions
3. Methodological Contributions:
   • Application of Frank-Lieb rearrangement inequality to fluid dynamics
   • Development of "bucket" decomposition topological method on universal covers
   • Establishment of quantitative connections between symmetry, energy, and geometry

Limitations

1. Specificity of Construction: The construction in Theorem 1.4 depends on specific 3-fold symmetric shape (disk with three slender arms); unclear whether more general simply-connected vortex patches exhibit perimeter growth
2. Growth Rate: Only proves linear growth $\text{Perimeter}(\Omega_t) \gtrsim t$, does not provide precise growth rate constants or upper bounds
3. Symmetry Requirement: Theorem 1.2 requires $m \geq 2$ symmetry ($m=1$, i.e., no symmetry case not covered)
4. Angular Momentum Exponent: In Theorem 1.3 with bounded angular momentum, the exponent is 4; although Proposition 3.2 suggests exponent 2 is optimal in some sense, the stability theorem does not achieve this exponent
5. Quantitative Constants: Many estimates use $\lesssim$ notation without explicit constants, limiting practical applicability

Future Directions

1. General Simply-Connected Patches: Can we characterize all simply-connected vortex patches with perimeter growth? Do discrimination criteria exist independent of symmetry?
2. Precise Growth Rates: Determine exact asymptotic behavior of perimeter growth; might there exist superlinear growth $t^\alpha$ ($\alpha > 1$)?
3. Non-Symmetric Case: Establish stability estimates without symmetry assumptions; may require developing new geometric methods
4. Higher Dimensions: Analogous problems for vortex tubes/surfaces in three-dimensional Euler equations
5. Numerical Verification: Numerical simulation of constructions in this paper to observe actual perimeter growth behavior and possible filamentation phenomena

In-Depth Evaluation

Strengths

1. Theoretical Depth:
   • Cleverly combines harmonic analysis (rearrangement inequalities), geometric measure theory (perimeter estimates), and topology (universal covers)
   • Symmetry exploitation demonstrates profound geometric insight
   • Proofs are rigorous with clear logical flow
2. Problem Importance:
   • Resolves long-standing problem of simply-connected perimeter growth in vortex dynamics
   • Stability results significantly improve existing literature (exponent from 4 to 2, removes multiple restrictions)
   • Connects numerical simulation observations with rigorous mathematical theory
3. Method Innovation:
   • "Bucket" decomposition method is innovative technique for handling simply-connected topological constraints
   • Application of Frank-Lieb results opens new research pathways
   • Energy-momentum balance parameter design demonstrates refined analytical technique
4. Writing Quality:
   • Clear structure, progressively develops from motivation to results
   • Figures (Figures 1 and 2) effectively aid understanding
   • Appendix contains complete proofs, enhancing readability

Weaknesses

1. Non-Explicit Construction:
   • Theorem 1.4 construction involves multiple parameters ($N, \gamma, R, \delta, \epsilon$) with relationships given implicitly through $\ll$ and $\gg$
   • Lacks concrete numerical examples or algorithms for parameter selection
2. Technical Dependencies:
   • Critically depends on deep result Frank-Lieb 5 (Proposition 1.1), but paper provides no self-contained proof
   • Lemma 4.1 proof cites Elgindi 4's inequality, increasing reading difficulty
3. Qualitative Nature of Results:
   • Perimeter growth only gives $\gtrsim t$, does not determine constants or provide upper bounds
   • Stability constants $C$ not explicitly computed, limiting practical application
4. Symmetry Restrictions:
   • Theorem 1.2 requires $m \geq 2$, excluding many natural non-symmetric perturbations
   • Theorem 1.4 uses $m=3$, does not explore other symmetries
5. Comparison with Numerical Simulations:
   • While claiming to "confirm" Dritschel's simulations, provides no quantitative comparison
   • Constructed patch (Figure 1) may differ significantly from shapes observed in actual simulations

Impact

1. Theoretical Impact:
   • Provides new tools for stability theory in vortex dynamics (symmetry exploitation, rearrangement inequalities)
   • Resolves existence of simply-connected perimeter growth, potentially inspiring further constructions
   • Methods may generalize to other conservation law systems
2. Practical Value:
   • Provides theoretical foundation for understanding evolution of vortex structures in turbulence
   • Stability estimates applicable to error analysis of numerical methods
   • Perimeter growth mechanism may relate to mixing, dissipation, and other physical phenomena
3. Reproducibility:
   • All results are mathematical theorems, in principle verifiable by any mathematician
   • Construction parameter relationships are clear, enabling numerical implementation
   • Lack of code or concrete numerical examples may hinder rapid reproduction

Applicable Scenarios

1. Fluid Mechanics:
   • Long-time behavior of vortex structures in two-dimensional incompressible flow
   • Theoretical understanding of vortex merging, vortex filamentation phenomena
   • Simplified models for atmospheric and oceanic vortices
2. Mathematical Analysis:
   • Stability theory for nonlinear PDEs
   • Qualitative analysis of conservation law systems
   • Isoperimetric problems in geometric measure theory
3. Numerical Methods:
   • Theoretical foundation for vortex methods
   • Verification benchmark for long-time numerical simulations
   • Design of adaptive mesh strategies (for perimeter growth)
4. Limited Scenarios:
   • Not applicable to three-dimensional or compressible flow
   • Not applicable to bounded domains (methods depend on full space)
   • Not applicable to strong perturbations or far from Rankine vortex

Key References

1 J.-Y. Chemin (1991): Proves smoothness preservation of vortex patch boundaries
2 D. G. Dritschel (1988): Numerical simulations show repeated filamentation of vortex interfaces
3 T. D. Drivas, T. M. Elgindi, I.-J. Jeong (2024): First construction of linear perimeter growth for multiply-connected vortex patches
5 R. L. Frank, E. H. Lieb (2021): Spherical clustering proof via quantitative rearrangement inequalities
7 Y. Tang (1987): Nonlinear stability of vortex patches using pseudo-energy
8 Y. H. Wan, M. Pulvirenti (1985): Nonlinear stability of circular vortex patches
10 V. I. Yudovich (1963): Unsteady flow of ideal incompressible fluid

Overall Assessment: This is a high-quality mathematical analysis paper making important contributions to vortex dynamics. By cleverly combining symmetry, energy methods, and topological arguments, the author not only improves classical stability results but also resolves the long-standing existence problem for simply-connected vortex patch perimeter growth. Despite limitations such as non-explicit construction and undetermined constants, the methodological innovation and theoretical depth make it significant progress in the field. The paper establishes a solid foundation for subsequent research (more general constructions, precise growth rates, numerical verification, etc.).