Time quasi-periodic vortex patches of Euler equation in the plane

平面上欧拉方程的时间准周期涡旋斑块

作者: Massimiliano Berti, Zineb Hassainia, Nader Masmoudi

类别: math.AP

发布日期: 2022-02-13

arXiv ID: 2202.06215v2

链接: http://arxiv.org/abs/2202.06215v2

详细分析

中文标题: 平面上欧拉方程的时间准周期涡旋斑块

1. 研究对象和背景

该论文研究的是二维不可压缩欧拉方程（2D Euler equations）在平面上 $\mathbb{R}^2$ 的涡旋斑块（vortex patch）解的存在性问题。欧拉方程描述了无粘性、不可压缩流体的运动，其形式为：

$$ \partial_t \omega + \nabla^\perp \psi \cdot \nabla \omega = 0, \quad \Delta \psi = \omega, $$

其中 $\omega(t,x)$ 是涡度（vorticity），$\psi(t,x)$ 是流函数（stream function），$\nabla^\perp \psi = (-\partial_y \psi, \partial_x \psi)$ 是速度场。涡旋斑块是指涡度 $\omega$ 为某个有界单连通区域 $D_t \subset \mathbb{R}^2$ 的特征函数，即 $\omega(t,\cdot) = \chi_{D_t}(\cdot)$。

论文聚焦于从Kirchhoff 椭圆涡旋（Kirchhoff elliptical vortices）出发，研究其附近的时间准周期解。Kirchhoff 椭圆解是经典精确解，其初始区域为：

$$ D_0 = \left\{ (x,y) \in \mathbb{R}^2 : \frac{x^2}{\gamma} + \gamma y^2 \leq 1 \right\}, \quad \gamma > 1, $$

其中 $\gamma$ 为椭圆的长宽比（aspect ratio），该解以恒定角速度 $\Omega_\gamma = \frac{\gamma - 1}{\gamma + 1}$ 绕原点旋转，即 $D_t = e^{i \Omega_\gamma t} D_0$。

论文的关键在于将涡旋斑块的边界演化问题转化为轮廓动力学方程（contour dynamics equation）：设边界参数化为 $z(t,\theta) = e^{i \Omega t} w(t,\theta)$，则 $w(t,\theta)$ 满足：

$$ \partial_t w(t,\theta) \cdot \partial_\theta w(t,\theta) = i \partial_\theta \psi_\Omega(t, w(t,\theta)), $$

其中 $\psi_\Omega$ 是修正后的流函数，包含旋转项。进一步，作者引入径向扰动变量 $\xi(t,\theta)$，将边界参数化为：

$$ w(t,\theta) = \left(1 + 2\xi(t,\theta)\right)^{1/2} \left( \gamma^{-1/2} \cos \theta + i \gamma^{1/2} \sin \theta \right), $$

从而将原问题转化为关于 $\xi$ 的准线性哈密顿偏微分方程：

$$ \partial_t \xi = \partial_\theta \nabla H_\Omega(\xi), $$

其中哈密顿量 $H_\Omega(\xi) = \frac{1}{2} E(\xi) - \frac{\Omega}{2} J(\xi)$，$E(\xi)$ 为伪能量，$J(\xi)$ 为角动量（prime integral）。该方程定义在函数空间 $H^s(\mathbb{T})$ 上（$\mathbb{T} = \mathbb{R}/2\pi\mathbb{Z}$），且要求解为偶函数（invariant under $\theta \mapsto -\theta$），以保持物理对称性。

值得注意的是，该问题不具有旋转对称性（即不是“动量守恒”系统），且在椭圆情形下，角动量守恒导致线性化系统中存在零特征频率（degenerate mode），这是主要数学障碍。

2. 主要定理或主要结果

论文的核心结果由定理 1.2（Theorem 1.2）给出，其内容如下：

定理 1.2（时间准周期涡旋斑块的存在性）：
对任意整数 $\bar{n} \geq 2$，考虑长宽比区间 $G = [\gamma_1, \gamma_2]$ 如公式 (1.16) 所定义（即当 $\bar{n} = 2$ 时为 $[1,3]$，当 $\bar{n} \geq 3$ 时为 $(\gamma_{\bar{n}}, \gamma_{\bar{n}+1})$，其中 $\gamma_n$ 是使第 $n$ 个法向频率为零的临界值）。令 $S \subset \{\bar{n}+1, \bar{n}+2, \dots\}$ 为有限个整数构成的“切向位点”集合。则存在常数 $\bar{s} > S + \frac{1}{2}$ 和 $\varepsilon_0 \in (0,1)$，使得对任意幅值向量 $\vec{a} = (a_n)_{n \in S}$ 满足 $|a_n| \leq \varepsilon_0$ 且 $a_n \neq 0$，以下结论成立：

  1. 存在一个子集 $G_{\vec{a}} \subset G$，其 Lebesgue 测度渐近完整，即：

    $$ \lim_{|\vec{a}| \to 0} \frac{|G_{\vec{a}}|}{|G|} = 1; $$

  2. 对任意 $\gamma \in G_{\vec{a}}$，存在：

    • 一个常数 $\mu(\gamma, \vec{a}) \in \mathbb{R}$，满足 $\mu \to 0$ 当 $\vec{a} \to 0$；
    • 一个 Diophantine 频率向量 $\omega(\gamma, \vec{a}) = (\omega_n)_{n \in S} \in \mathbb{R}^{|S|}$，满足 $\omega_n \to \Omega_n(\gamma)$ 当 $\vec{a} \to 0$，其中 $\Omega_n(\gamma)$ 是 Kirchhoff 椭圆的线性法向频率：
      $$ \Omega_n(\gamma) = \frac{\gamma - 1}{\gamma + 1} - \frac{1}{2} \left( \frac{1}{\gamma^{1/2}} - \frac{1}{\gamma^{-1/2}} \right) \frac{1}{n^2}, \quad n \geq 1; $$
    • 一个准周期解 $\xi(\omega t, \theta)$，属于 Sobolev 空间 $H^{\bar{s}}(\mathbb{T}^{|S|} \times \mathbb{T})$，且为偶函数，具有如下形式：
      $$ \xi(\omega t, \theta) = \sum_{n \in S} a_n M_n \cos(\omega_n t) \cos(n\theta) + a_n M_{n-1} \sin(\omega_n t) \sin(n\theta) + r(\omega t, \theta), $$
      其中 $M_n$ 由公式 (1.18) 给出，余项 $r$ 满足 $\|r\|_{\bar{s}} \to 0$ 当 $\vec{a} \to 0$。

由此，通过公式 (1.11)，可构造出欧拉方程 (1.1) 的准周期涡旋斑块解 $z(t,\theta)$，其边界围绕椭圆做准周期变形，且不在任何旋转参考系中为稳态。

该定理首次证明了二维欧拉方程中存在任意多个频率的准周期涡旋斑块解，且这些解从 Kirchhoff 椭圆族中“分岔”（bifurcate）而来，其存在性依赖于长宽比 $\gamma$ 属于一个测度渐近完整的集合。

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

a) 核心分析学工具：

  • KAM 理论框架：论文采用经典的 KAM（Kolmogorov-Arnold-Moser）扰动方法，旨在证明在小扰动下，不变环面（即准周期解）能够持久存在。
  • Nash-Moser 迭代：由于方程为准线性（quasi-linear），且存在小除数问题（small divisors），作者采用 Nash-Moser 迭代方案，而非标准的 Banach 不动点定理，以克服正则性损失。
  • 辛几何结构：方程具有哈密顿结构，作者利用辛变换（symplectic change of variables）进行变量约化，这是 KAM 理论的核心工具。
  • 谱分析与特征频率计算：在 $\xi = 0$ 处线性化方程，通过 Fourier 展开得到无穷多个耦合的谐振子系统，并精确计算其频率 $\Omega_n(\gamma)$，识别出零频率模（$n=2$）和椭圆/双曲模的共存。
  • Tame 估计：所有非线性算子（包括积分核）均满足 Tame 估计，这是 Nash-Moser 理论中控制迭代收敛性的关键技术，确保 Sobolev 范数的损失可控。

b) 解决技术难题的主要技巧：

  • 角动量的辛约化（Symplectic Reduction of Angular Momentum）：这是本文最核心的创新。由于角动量 $J(\xi)$ 是守恒量，其线性部分在 $n=2$ 模上诱导出零特征频率，导致 KAM 理论失效。作者引入一个非线性辛变换 $\Phi$，将角动量 $J$ 作为新的辛变量 $\alpha_2$，并构造共轭变量 $\beta_2 = \bar{t}(\xi)$（即流 $\Phi_t^J$ 与 $\beta_2=0$ 截面的首次撞击时间），从而将原系统在 $\xi=0$ 附近“拉直”为：
    $$ \dot{\alpha}_2 = 1, \quad \dot{\beta}_2 = 0, \quad \dot{\alpha}_n = \mu_n \beta_n, \quad \dot{\beta}_n = -\mu_n \alpha_n \quad (n \neq 2), $$
    从而消除退化。这是首次在 PDE 的 KAM 理论中实现此类基于守恒量的辛约化。
  • 非线性流的构造与估计：变换 $\Phi$ 是由角动量生成的非线性输运流 $\Phi_t^J$ 的非线性修正，其逆变换 $\Phi^{-1}$ 和微分 $d\Phi, d\Phi^{-1}$ 的 Tame 估计极其困难，作者通过精细的积分算子分析和傅里叶乘子估计完成。
  • 避免中间余项：由于线性化算子具有显式结构，且在辛变换后仍保持良好形式，作者在一次操作中（Section 10.3）就将线性算子约化为常系数形式，避免了传统 KAM 中多步共轭的复杂余项累积。
  • 处理双曲模与椭圆模共存：对于 $\gamma > 3$，存在有限个双曲模（hyperbolic modes），这在以往涡旋斑块研究中罕见。作者在 Section 11 中通过构造可逆的辛变换，将双曲模纳入不变环面的“稳定/不稳定流形”结构中，成功处理了其对 KAM 迭代的影响。

4. 与之前工作的比较

a) 突破与改进：

  • 首次证明准周期涡旋斑块解的存在性：此前所有关于欧拉方程涡旋斑块的 KAM 结果（如 [40, 48, 42]）均局限于圆形 Rankine 涡旋（$\gamma=1$），此时角动量守恒导致的退化模（$n=2$）不存在，且系统具有旋转对称性。本文首次处理椭圆涡旋（$\gamma > 1$），并克服了其退化模带来的根本性障碍。
  • 突破旋转对称性假设：以往 KAM 结果依赖于方程的旋转不变性（即动量守恒），使得角动量为平凡常数。本文首次在非动量守恒系统中成功实施 KAM，通过将角动量作为动态变量引入，实现辛约化，极大拓展了 KAM 理论的应用范围。
  • 利用内参数 $\gamma$ 满足 Melnikov 条件：在大多数 KAM 对 PDE 的应用中，需依赖外部参数（如表面张力、深度、势能）来避开共振。本文首次利用内禀几何参数 $\gamma$（长宽比）——尽管其对频率的影响是指数小的（$O(e^{-c n})$）——来验证无穷多个非共振条件（Melnikov 条件），这是极强的创新。
  • 首次在无限维 PDE 中实现 Darboux-Carathéodory 辛约化：该思想在有限维力学（如太阳系稳定性）中已有，但本文首次将其推广至无限维哈密顿 PDE，且构造了显式、可估计的非线性辛变换，为后续研究提供了范式。
  • 处理双曲模的共存：以往涡旋斑块 KAM 研究中所有法向模均为椭圆型，而本文处理了椭圆与双曲模并存的情形，为研究混沌行为（如 Jupiter 大红斑）提供了数学基础。

综上，本文不仅证明了新的解的存在性，更重要的是开创了一种全新的 KAM 技术路径，为解决其他具有守恒量退化的哈密顿 PDE（如磁流体、量子流体、三维欧拉）提供了关键工具。