Stability of the Rankine Vortex and Perimeter Growth in Vortex Patches

基本信息

• 论文ID: 2511.09772
• 标题: Stability of the Rankine Vortex and Perimeter Growth in Vortex Patches
• 作者: John Brownfield
• 分类: math.AP (Analysis of PDEs)
• 发表时间: November 2025
• 论文链接: https://arxiv.org/abs/2511.09772

摘要

本文证明了对于与Rankine涡具有相同总涡量和涡量中心的涡量函数 $\omega: \mathbb{R}^2 \to [0,1]$，其与Rankine patch的 $L^1$ 偏差可由伪能量偏差和角动量的函数界定。在具有 $m$-重对称性的情况下，可以去除对角动量的依赖。利用这一结果，作者证实了先前数值模拟的结果，展示了Rankine涡的单连通扰动具有线性时间的周长增长。

研究背景与动机

问题背景

本文研究二维不可压缩Euler方程的涡量形式： $\partial_t\omega + u \cdot \nabla\omega = 0$ 其中速度场 $u$ 由Biot-Savart定律给出：$u = K * \omega$，$K(x) = \frac{x^\perp}{2\pi|x|^2}$。

核心问题

1. 正则性损失问题：尽管Chemin 1证明了光滑初始边界保持光滑，但近期研究关注构造在无限时间内失去初始正则性的解。这些解通常是已知简单解的稳定扰动（在弱范数意义下）。
2. Rankine涡的稳定性：Rankine涡patch $\omega^* = 1_{B(a,r)}$ 是稳定态，因为圆盘是固定总涡量和涡量中心下角动量的唯一极小值。先前研究（Pulvirenti-Wan 8, Sideris-Vega 6, Tang 7）建立了定量 $L^1$ 稳定性估计，但存在限制。
3. 周长增长：Dritschel 2的数值模拟显示涡patch可能出现持续的周长增长。Drivas-Elgindi-Jeong 3首次构造了具有线性时间周长增长的多连通涡patch，但单连通情况仍未解决。

研究动机

• 改进现有稳定性估计：去除初始 $L^1$ 距离限制，放松patch假设到 $0 \leq \omega \leq 1$
• 在 $m$-重对称情况下，去除对角动量的依赖并改进指数
• 构造首个具有线性时间周长增长的单连通涡patch，填补理论空白

核心贡献

1. 改进的稳定性定理（定理1.2和1.3）：
   • 在 $m$-重对称下，证明 $\|\omega_t - 1_{B(0,r)}\|_1^2 \leq C(E(1_{B(0,r)}) - E(\omega_0))$，去除角动量依赖
   • 在有界角动量下，证明 $\|\omega_t - 1_{B(0,r)}\|_1^4 \leq C(I,r)(E(1_{B(0,r)}) - E(\omega_0))$
2. 单连通涡patch的周长增长（定理1.4）：
   • 首次构造单连通开集 $\Omega_0$，其涡patch解满足 $\text{Perimeter}(\Omega_t) \gtrsim t$
   • 证实了Dritschel数值模拟的理论预测
3. 技术创新：
   • 利用Frank-Lieb 5的重排不等式结果（命题1.1）
   • 发展了处理单连通情况的拓扑论证方法（"桶"分解技术）

方法详解

任务定义

输入：初始涡量分布 $\omega_0: \mathbb{R}^2 \to [0,1]$，满足：

• 总涡量：$\int \omega_0 = \pi r^2$
• 涡量中心：$\int x\omega_0 = 0$
• 对称性或角动量约束

输出：

1. $L^1$ 稳定性估计：$\|\omega_t - 1_{B(0,r)}\|_1$ 的上界
2. 周长增长率：$\text{Perimeter}(\Omega_t)$ 的下界

稳定性分析架构

第一步：重排不等式（命题1.1）

利用Frank-Lieb的结果，对于 $\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$

其中 $E^*$ 是以原点为中心、测度为 $\int \omega$ 的球。

第二步：几何不等式（命题3.1和3.2）

命题3.1（$m$-重对称）：若 $\omega = \omega \circ R_{2\pi/m}$（$m \geq 2$），则 $\inf_{a\in\mathbb{R}^2} \|\omega - 1_{E^*+a}\|_{L^1} \geq \frac{1}{3}\|\omega - 1_{E^*}\|_{L^1}$

证明思路：

• 对于偶数 $m$，利用 $\|\omega - 1_{E^*+a}\|_{L^1} = \|\omega - 1_{E^*-a}\|_{L^1}$
• 考察 $B(a,r) \cap B(-a,r) \subset B(0,r)$，通过测度计算得到 $\epsilon \geq \delta/2$
• 对于奇数 $m = 2n+1$，构造三个球 $B(a_1,r), B(a_2,r), B(a_3,r)$，其三重交集包含于 $B(0,r)$
• 通过包含-排斥原理得到 $\epsilon \geq \delta/3$

命题3.2（有界角动量）：若 $\int |x|^2\omega \leq I$，$\int x\omega = 0$，则 $\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$

且指数2是最优的。

第三步：稳定性定理证明

结合命题1.1和命题3.1/3.2，利用伪能量守恒： $\|\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))$

周长增长构造

初始设置（第4.1节）

构造3-重对称单连通集 $\Omega_0$：

• 核心：半径约为1的圆盘 $B(0, 1-\epsilon)$
• 三个细长臂：长度 $N \gg 1$，总面积 $\gamma \ll 1$
• 参数关系：$\gamma N^2 \sim R^3$（角动量），$\gamma \ln(N) \sim \delta$（伪能量偏差）

关键性质：

1. 总涡量：$\int \omega_0 = \pi$
2. 角动量：$\int |x|^2\omega_0 \sim 1 + \gamma N^2$（大角动量保证 $\Omega_t \backslash B(0,R) \neq \emptyset$）
3. 伪能量：$E(1_{B(0,1)}) - E(\omega_0) \lesssim \delta \ll 1$（保证 $L^1$ 稳定性）

速度场估计（引理4.1）

设 $u = u_r e_r + ru_\theta e_\theta$ 是 $\omega_t$ 诱导的速度场，$u^* = r\mu(r)e_\theta$ 是 $1_{B(0,1)}$ 诱导的速度场，其中 $\mu(r) = \begin{cases} \frac{1}{2} & r \leq 1 \\ \frac{1}{2r} & r \geq 1 \end{cases}$

则有估计：

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}$

拓扑论证（第4.2节）

万有覆盖提升：

• 将流映射 $\Phi$ 提升到 $\mathbb{R}^2\backslash\{0\}$ 的万有覆盖 $\mathbb{R}\times\mathbb{R}^+$
• 提升映射 $\tilde{\Phi} = (\Phi_r, \Phi_\theta)$ 满足： $\|\Phi_\theta - \theta - t\mu(\Phi_r)\|_{L^2} \lesssim \epsilon \cdot t$

"桶"分解技术：

• 定义 $M_t^n$ 为 $\partial\tilde{\Omega}_t^n$ 的首个（从左到右）径向最大值点
• 从每个 $M_t^n$ 向上画线，将 $\Omega_t^c$ 分解为不相交的"桶" $B_t^n$
• 对于点 $z_0 \in \Omega_0^c$，定义 $N(t)$ 为包含 $Z_t = \tilde{\Phi}(z_0, t)$ 的桶标签

关键引理4.2： $|N(T)| \lesssim (1 + \max(1/R, \epsilon)T) \lesssim (\nu + o(1))T$ 其中 $\nu = \max(1/R, \epsilon) \ll 1$。

证明思路：

• 选择时间步长 $\Delta \sim \min(R, \epsilon^{-1})$，使得在半径 $\geq R$ 处角度变化 $\leq 1$
• 证明桶 $B_t^n$ 在时间 $\Delta$ 内的演化被限制在 $B_{t+\Delta}^{n-2} \cup \cdots \cup B_{t+\Delta}^{n+2}$
• 分割 $[0,T]$ 为长度 $\Delta$ 的区间，得到 $|N(T)| \lesssim 1 + T/\Delta$

周长增长证明

1. 在 $\Omega_0^c$ 中找到两点 $p_0, q_0$，满足：
   • $\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. 应用引理4.2，找到 $\partial\Omega_0$ 上的点 $p, q$ 满足：
   • $\Phi_\theta(q, T) < c_0 T < c_1 T < \Phi_\theta(p, T)$（角度差异）
   • $\Phi_r(p, T), \Phi_r(q, T) > r_0 > 0$（远离原点）
3. 利用几何引理（来自3），连接 $\Phi(p,T)$ 和 $\Phi(q,T)$ 的 $\partial\Omega_T$ 段长度至少为 $2r_0(c_1-c_0)T - 1 \gtrsim T$

技术创新点

1. 几何-分析结合：将Frank-Lieb的抽象重排不等式转化为具体的涡量稳定性估计
2. 对称性利用：$m$-重对称的巧妙应用避免了角动量依赖，这是对Tang 7工作的实质性改进
3. 拓扑方法：万有覆盖上的"桶"分解是处理单连通情况的关键创新，解决了3中多连通方法无法推广的难题
4. 能量-动量平衡：通过精心选择参数（$\gamma, N, R, \delta$），平衡了大角动量（保证径向支撑）和小伪能量偏差（保证 $L^1$ 稳定性）

实验设置

理论证明性质

本文为纯理论数学论文，不包含数值实验。所有结果通过严格的数学证明获得。

构造参数

对于定理1.4的构造：

• 对称性：$m = 3$
• 核心半径：$r \approx 1 - \epsilon$
• 臂长度：$N \gg 1$
• 臂面积：$\gamma \ll 1$
• 参数关系：$\gamma N^2 \sim R^3$，$\gamma \ln(N) \sim \delta$
• 约束：$\delta \ll 1$，$R \gg 1$，$\epsilon = \max(\tilde{\epsilon}^{1/2}, \ln(R)\tilde{\epsilon}^{1/4}) \ll 1$

实验结果

主要理论结果

定理1.2（$m$-重对称稳定性）：

• 存在普适常数 $C > 0$，对于满足 $\omega_0 = \omega_0 \circ R_{2\pi/m}$（$m \geq 2$）和 $\int \omega_0 = \pi r^2$ 的 $\omega_0: \mathbb{R}^2 \to [0,1]$
• 解满足：$\|\omega_t - 1_{B(0,r)}\|_1^2 \leq C(E(1_{B(0,r)}) - E(\omega_0))$，对所有时间成立
• 改进：相比Tang 7，去除了初始 $L^1$ 距离限制、角动量依赖，并将指数从4改进到2

定理1.3（有界角动量稳定性）：

• 对于 $I, r > 0$，存在 $C = C(I,r) > 0$，若 $\omega_0$ 满足 $\int |x|^2\omega_0 \leq I$，$\int x\omega_0 = 0$，$\int \omega_0 = \pi r^2$
• 解满足：$\|\omega_t - 1_{B(0,r)}\|_1^4 \leq C(E(1_{B(0,r)}) - E(\omega_0))$
• 改进：去除初始 $L^1$ 距离限制，放松patch假设

定理1.4（线性周长增长）：

• 存在单连通开集 $\Omega_0 \subset \mathbb{R}^2$（光滑边界），使得涡patch解 $\omega_t = 1_{\Omega_t}$ 满足： $\text{Perimeter}(\Omega_t) \gtrsim t$ 对充分大的 $t$ 成立
• 突破：首个单连通涡patch的线性周长增长构造

与现有工作的比较

理论发现

1. 对称性的威力：$m$-重对称（$m \geq 2$）足以显著改进稳定性估计，这表明对称性在涡动力学中起着基础性作用
2. 指数最优性：命题3.2证明了有界角动量情况下指数2是最优的，这揭示了该问题的本质困难
3. 拓扑约束：单连通性通过"桶"分解机制限制了点绕原点的缠绕速度，这是实现周长增长的关键
4. 能量-几何关系：伪能量偏差直接控制 $L^1$ 偏差，建立了能量泛函与几何性质的定量联系

相关工作

涡patch稳定性

1. Pulvirenti-Wan 8 (1985)：首次利用角动量和伪能量在有界域上建立定量 $L^1$ 稳定性
2. Sideris-Vega 6 (2009)：将角动量方法推广到全平面
3. Tang 7 (1987)：使用伪能量在全平面建立稳定性，但需要初始接近性和有界角动量
4. 本文：改进指数、去除限制、利用对称性

正则性理论

1. Yudovich 10 (1963)：证明 $\omega_0 \in L^1 \cap L^\infty$ 下全局唯一弱解存在
2. Chemin 1 (1991)：证明光滑初始边界保持光滑

周长增长

1. Dritschel 2 (1988)：数值模拟显示涡界面的反复丝化和周长增长
2. Drivas-Elgindi-Jeong 3 (2024)：首次构造多连通涡patch的线性周长增长，利用 $L^2$ 剪切稳定性估计
3. 本文：首次构造单连通情况，发展新的拓扑论证方法

重排不等式

1. Frank-Lieb 5 (2021)：证明基于定量重排不等式的球形群集定理
2. Yan-Yao 9 (2022)：指出5的结果蕴含命题1.1
3. 本文：将这些抽象结果应用于涡动力学

结论与讨论

主要结论

1. 稳定性改进：在 $m$-重对称下，Rankine涡的 $L^1$ 稳定性可以仅由伪能量偏差控制，指数为2，无需角动量约束
2. 单连通周长增长：首次严格证明存在单连通涡patch具有线性时间周长增长，证实了Dritschel数值模拟的预测
3. 方法论贡献：
   • 将Frank-Lieb重排不等式应用于流体动力学
   • 发展了万有覆盖上的"桶"分解拓扑方法
   • 建立了对称性、能量和几何之间的定量联系

局限性

1. 构造的特殊性：定理1.4的构造依赖于特定的3-重对称形状（带三个细臂的圆盘），尚不清楚是否存在更一般的单连通涡patch具有周长增长
2. 增长率：仅证明了线性增长 $\text{Perimeter}(\Omega_t) \gtrsim t$，未给出精确的增长率常数或上界
3. 对称性要求：定理1.2需要 $m \geq 2$ 的对称性（$m=1$ 即无对称性的情况未覆盖）
4. 角动量指数：定理1.3中有界角动量情况的指数为4，虽然命题3.2表明指数2在某种意义下最优，但稳定性定理未达到这一指数
5. 定量常数：许多估计使用了 $\lesssim$ 符号，未给出显式常数，这限制了结果的实际应用

未来方向

1. 一般单连通patch：能否刻画所有具有周长增长的单连通涡patch？是否存在不依赖对称性的判别准则？
2. 精确增长率：确定周长增长的精确渐近行为，可能存在 $t^\alpha$（$\alpha > 1$）的超线性增长？
3. 无对称性情况：在无对称性假设下建立稳定性估计，可能需要发展新的几何方法
4. 高维推广：三维Euler方程中涡管/涡面的类似问题
5. 数值验证：对本文构造进行数值模拟，观察周长增长的实际行为和可能的丝化现象

深度评价

优点

1. 理论深度：
   • 巧妙结合了调和分析（重排不等式）、几何测度论（周长估计）和拓扑学（万有覆盖）
   • 对称性的利用展现了深刻的几何洞察
   • 证明技术严谨，逻辑清晰
2. 问题重要性：
   • 解决了涡动力学中长期存在的单连通周长增长问题
   • 稳定性结果显著改进了现有文献（指数从4到2，去除多个限制）
   • 连接了数值模拟观察与严格数学理论
3. 方法创新性：
   • "桶"分解方法是处理单连通拓扑约束的创新技术
   • Frank-Lieb结果的应用开辟了新的研究路径
   • 能量-动量平衡的参数设计体现了精细的分析技巧
4. 写作质量：
   • 结构清晰，从动机到结果逐步展开
   • 图示（图1和图2）有效辅助理解
   • 附录包含完整证明，增强可读性

不足

1. 构造的非显式性：
   • 定理1.4的构造涉及多个参数（$N, \gamma, R, \delta, \epsilon$），它们之间的关系通过 $\ll$ 和 $\gg$ 隐式给出
   • 缺少具体数值示例或参数选择的算法
2. 技术依赖性：
   • 关键依赖Frank-Lieb 5的深刻结果（命题1.1），但本文未提供自包含证明
   • 引理4.1的证明引用了Elgindi 4的不等式，增加了阅读门槛
3. 结果的定性性质：
   • 周长增长仅给出 $\gtrsim t$，未确定常数或给出上界
   • 稳定性常数 $C$ 未显式计算，限制了实际应用
4. 对称性限制：
   • 定理1.2需要 $m \geq 2$，排除了许多自然的非对称扰动
   • 定理1.4使用 $m=3$，未探讨其他对称性的可能性
5. 与数值模拟的比较：
   • 虽然声称"证实"Dritschel的模拟结果，但未提供定量比较
   • 构造的patch（图1）与实际模拟中观察到的形状可能有显著差异

影响力

1. 理论影响：
   • 为涡动力学中的稳定性理论提供了新工具（对称性利用、重排不等式）
   • 解决了单连通周长增长的存在性问题，可能启发更多构造
   • 方法可能推广到其他守恒律系统
2. 实用价值：
   • 为理解湍流中涡结构的演化提供理论基础
   • 稳定性估计可用于数值方法的误差分析
   • 周长增长机制可能与混合、耗散等物理现象相关
3. 可复现性：
   • 所有结果为数学定理，原则上可由任何数学家验证
   • 构造的参数关系清晰，可进行数值实现
   • 缺少代码或具体数值示例可能阻碍快速复现

适用场景

1. 流体力学：
   • 二维不可压缩流动中涡结构的长时间行为
   • 涡合并、涡丝化等现象的理论理解
   • 大气和海洋涡旋的简化模型
2. 数学分析：
   • 非线性PDE的稳定性理论
   • 守恒律系统的定性分析
   • 几何测度论中的等周问题
3. 数值方法：
   • 涡方法（vortex methods）的理论基础
   • 长时间数值模拟的验证基准
   • 自适应网格策略的设计（针对周长增长）
4. 局限场景：
   • 不适用于三维流动或可压缩流动
   • 不适用于有边界的区域（方法依赖全空间）
   • 不适用于强扰动或远离Rankine涡的情况

参考文献（关键文献）

1 J.-Y. Chemin (1991): 证明涡patch边界的光滑性保持
2 D. G. Dritschel (1988): 数值模拟显示涡界面的反复丝化
3 T. D. Drivas, T. M. Elgindi, I.-J. Jeong (2024): 首次构造多连通涡patch的线性周长增长
5 R. L. Frank, E. H. Lieb (2021): 基于定量重排不等式的球形群集证明
7 Y. Tang (1987): 使用伪能量的涡patch非线性稳定性
8 Y. H. Wan, M. Pulvirenti (1985): 圆形涡patch的非线性稳定性
10 V. I. Yudovich (1963): 理想不可压缩流体的非定常流动

总体评价：这是一篇高质量的数学分析论文，在涡动力学领域做出了重要贡献。通过巧妙结合对称性、能量方法和拓扑论证，作者不仅改进了经典的稳定性结果，还解决了单连通涡patch周长增长的存在性问题。尽管存在构造非显式、常数未确定等局限，但其方法论创新和理论深度使其成为该领域的重要进展。论文为后续研究（更一般的构造、精确增长率、数值验证等）奠定了坚实基础。