Rankine涡旋的稳定性与涡旋斑块边界的线性增长

Stability of the Rankine Vortex and Perimeter Growth in Vortex Patches

作者: John Brownfield

详细分析

中文标题: Rankine涡旋的稳定性与涡旋斑块边界的线性增长
英文标题: Stability of the Rankine Vortex and Perimeter Growth in Vortex Patches
作者: John Brownfield

1. 研究对象和背景

本文研究二维不可压欧拉方程（2D incompressible Euler equations）在**涡旋斑块（vortex patch）**初值下的长时间动力学行为，核心聚焦于两个相互关联的问题：
（i）Rankine涡旋的定量 $L^1$ 稳定性；
（ii）单连通涡旋斑块边界长度随时间的线性增长现象（即 $\mathrm{Perimeter}(\Omega_t) \gtrsim t$）。

二维不可压欧拉方程的涡度形式在整平面 $\mathbb{R}^2$ 上给出为：

$$ \partial_t \omega + u \cdot \nabla \omega = 0, $$

$$ \omega|_{t=0} = \omega_0, $$

其中速度场 $u$ 由Biot–Savart律通过涡度 $\omega$ 卷积生成：

$$ u = K * \omega, \quad K(x) = \frac{x^\perp}{2\pi |x|^2} = \frac{1}{2\pi |x|^2}(-x_2, x_1). $$

该系统在经典意义下不适定，但 Yudovich [10] 证明：若初始涡度 $\omega_0 \in L^1(\mathbb{R}^2) \cap L^\infty(\mathbb{R}^2)$，则存在全局唯一的弱解 $\omega \in L^\infty_{\mathrm{loc}}([0,\infty); L^1 \cap L^\infty)$。特别地，当 $\omega_0 = \mathbf{1}_{\Omega_0}$ 是某有限测度集合 $\Omega_0 \subset \mathbb{R}^2$ 的特征函数时，解始终保持为特征函数形式 $\omega(t) = \mathbf{1}_{\Omega_t}$，称为涡旋斑块解。由于 $\omega$ 满足输运方程，其支撑集边界 $\partial \Omega_t$ 随拉格朗日流 $\Phi_t$ 演化；Chemin [1] 证明：若 $\partial \Omega_0$ 光滑，则 $\partial \Omega_t$ 对所有有限时间保持光滑。

然而，一个深刻现象是：即使初始边界极其光滑，其长时间演化可能诱导边界正则性的渐近退化——例如边长无界增长、曲率爆破或Filamentation（细丝化）。Dritschel [2] 的数值模拟早已揭示涡旋斑块在演化中自发产生细长结构并持续拉伸；Drivas–Elgindi–Jeong [3] 首次严格构造了具有多连通拓扑的涡旋斑块，其周长满足 $\mathrm{Perimeter}(\Omega_t) \gtrsim t$。但其构造依赖于多个分离的斑块组件，无法推广至单连通情形——这正是本文要突破的核心难点。

本文研究的初值空间为：$\omega_0 : \mathbb{R}^2 \to [0,1]$，即允许“部分填充”的弱斑块（不强制为特征函数），但仍保持总涡量、涡心与角动量等守恒量的定义。关键守恒量包括：

  • 总涡量：$\displaystyle \int_{\mathbb{R}^2} \omega(x)\,dx$，
  • 涡心（center of vorticity）：$\displaystyle \int_{\mathbb{R}^2} x\,\omega(x)\,dx$，
  • 角动量（angular momentum）：$\displaystyle \int_{\mathbb{R}^2} |x|^2 \omega(x)\,dx$，
  • 动能（kinetic energy）：$\displaystyle \int_{\mathbb{R}^2} |u(x)|^2\,dx$（在 $\mathbb{R}^2$ 上对非紧支涡旋发散），
  • 伪能量（pseudo-energy）：
    $$ E(\omega) := \frac{1}{2\pi} \iint_{\mathbb{R}^2 \times \mathbb{R}^2} \omega(x)\ln(|x-y|^{-1})\,\omega(y)\,dx\,dy. $$

伪能量 $E(\omega)$ 是动能的积分重构，在 $\omega \in L^1 \cap L^\infty$ 下有限，且在固定总涡量与涡心条件下，圆盘斑块 $\mathbf{1}_{B(0,r)}$ 是 $E(\omega)$ 的极大值点（由重排不等式与球对称极小化原理保证）。因此，$E(\omega)$ 成为刻画斑块接近圆盘程度的天然“距离”泛函。

本文所有分析均在全平面 $\mathbb{R}^2$ 上进行，不加任何截断或周期化；函数空间取 $L^1 \cap L^\infty$，初值满足 $0 \le \omega_0 \le 1$，无需初始边界光滑性假设（仅需可测）。这一设定显著弱于传统涡旋斑块理论（如 Chemin 要求 $C^{1,\gamma}$ 边界），使稳定性与增长结果更具鲁棒性。

2. 主要定理或主要结果

论文建立了三个核心定理，构成逻辑链：稳定性 → 角动量控制 → 周长线性增长。

Theorem 1.2（$m$-重对称下的稳定性）
存在普适常数 $C > 0$，使得若 $\omega_0 : \mathbb{R}^2 \to [0,1]$ 满足：

  • $m$-重旋转对称：$\omega_0 = \omega_0 \circ R_{2\pi/m}$，其中 $R_\theta$ 为绕原点的旋转，
  • 总涡量 $\displaystyle \int_{\mathbb{R}^2} \omega_0(x)\,dx = \pi r^2$，

则其对应的 2D Euler 解 $\omega_t$ 满足对所有 $t \ge 0$：

$$ \|\omega_t - \mathbf{1}_{B(0,r)}\|_{L^1}^2 \le C\big(E(\mathbf{1}_{B(0,r)}) - E(\omega_0)\big). $$

Theorem 1.3（有界角动量下的稳定性）
给定 $I, r > 0$，存在 $C = C(I,r) > 0$，使得若 $\omega_0 : \mathbb{R}^2 \to [0,1]$ 满足：

  • $\displaystyle \int_{\mathbb{R}^2} |x|^2 \omega_0(x)\,dx \le I$，
  • 涡心在原点：$\displaystyle \int_{\mathbb{R}^2} x\,\omega_0(x)\,dx = 0$，
  • 总涡量 $\displaystyle \int_{\mathbb{R}^2} \omega_0(x)\,dx = \pi r^2$，

则其解满足：

$$ \|\omega_t - \mathbf{1}_{B(0,r)}\|_{L^1}^4 \le C\big(E(\mathbf{1}_{B(0,r)}) - E(\omega_0)\big). $$

Theorem 1.4（单连通斑块的线性周长增长）
存在光滑单连通开集 $\Omega_0 \subset \mathbb{R}^2$，使得其对应的涡旋斑块解 $\omega_t = \mathbf{1}_{\Omega_t}$ 满足：

$$ \mathrm{Perimeter}(\Omega_t) \gtrsim t, \quad \text{当 } t \gg 1. $$

该结果的关键在于构造了一个3-重对称、高角动量、伪能量接近 Rankine 斑块的初值：它由单位圆盘 $B(0,1-\varepsilon)$ 加上三条细长臂（总长 $\sim N \gg 1$，总面积 $\gamma \ll 1$）组成，从而同时满足：

  • $E(\mathbf{1}_{B(0,1)}) - E(\omega_0) \lesssim \gamma \ln N =: \delta$（伪能量偏差小），
  • $\displaystyle \int |x|^2 \omega_0 \sim \gamma N^2$（角动量大），
  • 由 Theorem 1.2，$\|\omega_t - \mathbf{1}_{B(0,1)}\|_{L^1} < \tilde{\varepsilon}$ 对所有 $t$ 成立（$L^1$ 稳定），
  • 由角动量守恒，$\Omega_t \setminus B(0,R) \neq \emptyset$ 对所有 $t$ 成立（径向支撑永不坍缩）。

3. 研究方法、关键技术和核心工具

a) 核心分析工具

本文的稳定性分析建立在**定量重排不等式（quantitative rearrangement inequalities）**之上，特别是 Frank–Lieb [5] 关于对数核 Riesz 势的稳定性估计。其核心引理（Proposition 1.1）断言：对任意 $\omega : \mathbb{R}^2 \to [0,1] \in L^1$，令 $E_* = B(0,r)$ 满足 $|E_*| = \|\omega\|_{L^1}$，则存在 $C>0$ 使得

$$ \iint_{E_* \times E_*} \ln(|x-y|^{-1})\,dx\,dy - \iint_{\mathbb{R}^2 \times \mathbb{R}^2} \omega(x)\ln(|x-y|^{-1})\omega(y)\,dx\,dy \ge C \inf_{a \in \mathbb{R}^2} \|\omega - \mathbf{1}_{E_* + a}\|_{L^1}^2. $$

该不等式将伪能量偏差控制为到最优平移圆盘的 $L^1$ 距离平方，是连接能量泛函与几何偏差的桥梁。其证明依赖于：

  • Riesz 重排不等式：指出对称递减重排最大化卷积型能量；
  • Bathtub 原理：在 $[0,1]$-约束下，极值函数必为特征函数；
  • 分层积分技巧：将 $\ln(|x-y|^{-1})$ 表示为 $\int_0^{2L} \frac{dR}{R} \mathbf{1}_{B_R}(x-y)$，再逐层应用 Frank–Lieb 定理于每个半径 $R$ 的球。

另一支柱是对称性诱导的 $L^1$ 平移稳定性（Proposition 3.1）：若 $\omega$ 具 $m$-重旋转对称性，则

$$ \inf_{a \in \mathbb{R}^2} \|\omega - \mathbf{1}_{B(0,r) + a}\|_{L^1} \ge \frac{1}{3} \|\omega - \mathbf{1}_{B(0,r)}\|_{L^1}. $$

该估计通过几何交集论证完成：利用 $m \ge 2$ 个平移圆盘 $B(a_i,r)$ 的两两交集必含于 $B(0,r)$，从而将平移误差 $\epsilon$ 控制为原始误差 $\delta$ 的分数倍。此性质彻底消除了对角动量有界性的依赖，是 Theorem 1.2 指数提升（平方而非四次）的关键。

在周长增长部分，核心工具是万有覆叠空间上的 lifted 动力学。将 $\mathbb{R}^2 \setminus \{0\}$ 的万有覆叠记为 $\mathbb{R}_+ \times \mathbb{R}$（极坐标 $(r,\theta)$ 的提升 $(r,\tilde{\theta})$，$\tilde{\theta} \in \mathbb{R}$），拉格朗日流 $\Phi_t$ 提升为 $\tilde{\Phi}_t = (\Phi_r, \Phi_\theta)$。关键估计为：

$$ \|\Phi_\theta(p,t) - \theta - t \mu(\Phi_r(p,t))\|_{L^2} \lesssim t \ln(R)^{1/2} \tilde{\varepsilon}^{1/4}, $$

其中 $\mu(r) = \begin{cases} r/2 & r \le 1 \\ 1/(2r) & r \ge 1 \end{cases}$ 是 Rankine 涡旋的角速度剖面。该估计结合了：

  • Biot–Savart 核的 $L^\infty$ 控制：$\|u - u_*\|_{L^\infty} \lesssim \tilde{\varepsilon}^{1/2}$，
  • 角向速度的 $L^2$ 逼近：$\|u_\theta - \mu\|_{L^2} \lesssim \ln(R) \tilde{\varepsilon}^{1/4}$，
  • 径向速度的小性：$|u_r(x)| \lesssim \tilde{\varepsilon}$。

b) 技术性难题的处理技巧

首要难题是从体积分（$L^1$ 稳定）过渡到边界几何（周长增长）。由于边界是零测集，经典 $L^p$ 估计失效。本文采用**桶分解（bucket decomposition）**技巧：在万有覆叠 $\mathbb{R}_+ \times \mathbb{R}$ 中，将 $\partial \widetilde{\Omega}_t$ 的所有提升 $\partial \widetilde{\Omega}_t^{(n)} = \partial \widetilde{\Omega}_t + (0,2\pi n)$ 连成一条简单曲线 $\Gamma_t$；以各提升的最右点 $M_n^{(t)}$ 作垂线，将 $\widetilde{\Omega}_t^c$ 分割为互不相交的“桶” $B_n^{(t)}$。定义 $N(t)$ 为点 $Z_t = \widetilde{\Phi}_t(z_0)$ 所在桶的编号。

关键引理（Lemma 4.2）断言：$|N(T)| \lesssim 1 + \max(1/R, \tilde{\varepsilon}) T$。其证明利用了：

  • 慢界面漂移：因 $|u_r| \lesssim \tilde{\varepsilon}$，桶的左右边界在时间 $\Delta \sim \min(R, \tilde{\varepsilon}^{-1})$ 内最多横向移动 1 单位；
  • 桶间隔离性：任一桶 $B_n^{(t)}$ 在 $\Delta$ 时间内只能扩散至相邻至多 2 个桶内，故 $|N(t+\Delta) - N(t)| \le 2$。

由此，对初值 $p_0 \in A_T$, $q_0 \in B_T$（分别位于不同径向区间），其提升点 $P_T$, $Q_T$ 的缠绕数差 $\Phi_\theta(P_T) - \Phi_\theta(Q_T) \gtrsim T$，而桶编号差受控，从而可回溯找到 $\partial \Omega_0$ 上两点 $p,q$，其轨迹在时间 $T$ 后满足 $\Phi_\theta(p,T) - \Phi_\theta(q,T) \gtrsim T$ 且 $\Phi_r(p,T), \Phi_r(q,T) \gtrsim 1$。最后调用 [3] 的几何引理 3.18：若两点在极坐标下角距 $\gtrsim T$、径距 $\gtrsim 1$，则其间边界弧长 $\gtrsim T$，完成线性增长证明。

4. 与之前工作的比较

a) 突破与改进

本文在三方面实现显著推进：

第一，放宽初值正则性与结构假设。
此前工作（如 Pulvirenti–Wan [8], Tang [7], Sideris–Vega [6]）均要求初值为严格特征函数 $\omega_0 = \mathbf{1}_{\Omega_0}$ 且 $\partial \Omega_0$ 光滑或至少 Lipschitz；而本文允许 $\omega_0 \in [0,1]$ 为任意可测函数，即“模糊斑块”。这使稳定性估计适用于更广的物理初值（如带涡量梯度的过渡层），且避免了边界正则性传播的复杂性。

第二，解除角动量有界性依赖。
Tang [7] 的稳定性估计要求 $\int |x|^2 \omega_0 \le I$，且误差控制为 $L^1$ 距离的四次方；本文 Theorem 1.2 在 $m$-对称下完全取消角动量假设，并将指数提升至二次方，达到信息论意义下的最优（因伪能量为二次型）。这是通过 Proposition 3.1 的纯几何对称论证实现的，不涉及任何谱分析或微分算子估计。

第三，首次实现单连通线性周长增长。
Drivas–Elgindi–Jeong [3] 的构造依赖多连通性：通过多个分离斑块间的剪切流拉伸边界。本文通过高角动量臂状结构创造“径向屏障”，结合万有覆叠桶分解，将剪切效应锁定在单连通边界的全局缠绕上。技术上，[3] 使用 $L^2$ 剪切估计但受限于零测集，而本文用桶编号 $N(t)$ 将点态缠绕转化为可测集运动，绕过测度论障碍。此方法仅依赖 $m$-对称、大角动量、伪能量接近三项，暗示 Dritschel [2] 数值中观察到的细丝化本质是对称性破缺驱动的缠绕放大，而非多连通拓扑专属。