Steady Contiguous Vortex-Patch Dipole Solutions of the 2D Incompressible Euler Equation

二维不可压缩欧拉方程的稳定连续涡斑偶极子解

作者: De Huang, Jiajun Tong

类别: math.AP, physics.flu-dyn, 35Q31, 35C07, 76B47,

发布日期: 2024-06-14

arXiv ID: 2406.09849v3

链接: http://arxiv.org/abs/2406.09849v3

详细分析

中文标题: 二维不可压缩欧拉方程的稳定连续涡斑偶极子解

1. 研究对象和背景

本论文聚焦于二维不可压缩欧拉方程的稳定行波解，特别是具有连续涡斑偶极子（contiguous vortex-patch dipole）结构的解。研究的方程组是二维不可压缩欧拉方程的涡度形式，其数学表述如下：

$$ \begin{aligned} \omega_t + \mathbf{u} \cdot \nabla \omega &= 0 \\ \mathbf{u} &= -\nabla^{\perp}(-\Delta)^{-1}\omega \end{aligned} $$

其中，$\omega(x, t)$ 是涡度（vorticity），$\mathbf{u}(x, t)$ 是速度场，$x = (x_1, x_2)$ 是空间坐标，$t$ 是时间。 $\nabla^{\perp} = (-\partial_{x_2}, \partial_{x_1})$ 是旋转梯度算子，$(-\Delta)^{-1}$ 是负拉普拉斯算子（Laplacian）的逆。

论文关注的是行波解，形式为 $\omega(\mathbf{x}, t) = \omega_0(\mathbf{x} - \mathbf{c}t)$，其中 $\mathbf{c}$ 是常数（向量）速度。通过变换到一个以速度 $\mathbf{c}$ 移动的合动参考系（co-moving frame），行波解 $\omega_0$ 满足定常二维不可压缩欧拉方程：

$$ \begin{aligned} \nabla^{\perp}\psi \cdot \nabla \omega &= 0 \\ -\Delta \psi &= \omega, \quad \lim_{|\mathbf{x}|\to\infty} -\partial_{x_2}\psi = c \end{aligned} $$

其中 $\psi$ 被称为流函数（stream function）。从这个意义上讲，行波解也被视为定常解，因为其形状不随时间改变。

研究背景可以追溯到流体力学中对涡旋结构（vortex structures）及其动力学的长期探索。尽管三维不可压缩欧拉方程的整体适定性（global well-posedness）仍是一个开放性难题，但二维和三维欧拉方程的定常解的研究在过去一个世纪里一直是该领域的研究热点。特别是，具有紧致支撑涡度（compactly supported vorticity）的非平凡定常解的构造和稳定性分析，很大程度上受到了真实物理实验和模拟刚性涡旋运动的数值模拟所揭示的现象的启发。

在三维情况下，定常涡环（steady vortex rings）的研究已有悠久历史，可以追溯到 Helmholtz 的工作。在二维情况下，具有紧致支撑涡度的定常解也引起了越来越大的兴趣。其中一个重要的非平凡显式解是Lamb-Chaplygin偶极子（Lamb-Chaplygin dipole），其涡度支撑在一个圆盘内，并在圆盘边界处连续消失。这是一个圆形的涡旋偶极子，可以看作是三维 Hill 球形涡旋的二维类比。论文中对“涡旋偶极子”的定义是指在一个非对称二分域（asymmetric bipartite domain）中支撑的、反向旋转的涡旋对（counter-rotating vortex pair）。

然而，论文特别关注的是涡斑解（vortex-patch solutions）。涡斑解的特点是涡度在紧致支撑的独立区域（patch-wise）内为常数。这类解因其能够模拟二维流动中的局部强旋转（localized strong rotations）而受到特别关注。尽管已经存在一些涡斑解的构造结果，例如单个旋转涡斑（V-state）或多涡斑系统，但论文指出，到目前为止，尚未理论上找到一个非平凡的（非径向）紧致支撑的定常涡斑解，其中两个或多个具有不同（常数值）涡度的涡斑区域相互接触。

这种“接触的涡斑”情景，特别是连续涡斑偶极子，早先由 Sadovskii 在其数值工作中提出（[Sad71]）。类似的问题也由 Gol’dshtik [Gol62] 和 Childress [Chi66] 早期考虑过。Sadovskii 的发现激发了大量后续的数值研究，旨在研究这类解的边界正则性（boundary regularity）、唯一性（uniqueness）以及其他性质。然而，对其存在的严格证明一直未能得到。有趣的是，Sadovskii 的场景也被发现在两个反对称涡环迎头碰撞（head-on collision）的大时间渐近剖面（large-time asymptotic profile）中是精确的近似。

因此，本论文的研究对象正是亟待解决的Sadovskii 场景——二维不可压缩欧拉方程的连续涡斑偶极子行波解。这种解的特点是，两个符号相反的涡斑具有共同的边界。对于方程的求解区域，文章主要在整个二维平面 $\mathbb{R}^2$ 上讨论，但为了构造解，特别关注了上半平面 $\mathbb{H} = \{(x_1, x_2) : x_2 \ge 0\}$，涉及零狄利克雷边界条件（zero Dirichlet boundary condition）在 $x_1$ 轴上，以及“远场消失”（vanishing condition in the far field）的条件。

2. 主要定理或主要结果

论文的核心贡献在于严格构造了二维不可压缩欧拉方程的首个稳定的连续涡斑偶极子行波解。这一结果通过以下Theorem 1.1 和 Theorem 2.2 来表述。

Theorem 1.1（摘要部分给出的初步陈述）： 存在一个单连通域 $D \subset \mathbb{R}^2$，该域具有正面积和 $C^1$ 光滑的边界，以原点为中心并关于两个坐标轴对称。在此域 $D$ 内，涡度 $\omega_D(\mathbf{x})$ 定义为：

$$ \omega_D(\mathbf{x}) = \text{sgn}(x_2) \cdot \mathbf{1}_D $$

这是一个关于定常二维不可压缩欧拉方程（1.2）的解，其流函数 $\psi$ 满足：

$$ \psi(x_1, 0) = 0, \quad x_1 \in \mathbb{R} $$

并且：

$$ \lim_{|\mathbf{x}|\to\infty} -\partial_{x_2}\psi = c_D $$

对于某个常数 $c_D > 0$。此外，该域 $D$ 的边界在 $x_1$ 轴之外是解析的（analytic away from the $x_1$-axis）。

Theorem 2.2（对Theorem 1.1的详细阐述）： 存在一个单连通域 $D \subset \mathbb{R}^2$，以原点为中心且关于 $x_1$ 轴和 $x_2$ 轴对称。该域 $D$ 具有正面积和 $C^1$ 光滑的边界。令 $\Omega = D \cap \mathbb{H}$（其中 $\mathbb{H}$ 是上半平面），并令 $c = c_\Omega$ 为由（2.6）给出的常数（速度）。则 $\omega_\Omega = \mathbf{1}_\Omega$ 是（2.2）式在 $\mathbb{H}$ 中的一个解（弱解），流函数为 $\psi = \phi - cx_2$。因此，$\omega_D = \text{sgn}(x_2) \cdot \mathbf{1}_D$ 是（1.1）式在 $\mathbb{R}^2$ 中的一个解，速度为 $c$。 令 $\Gamma = \partial D \cap \mathbb{H}$。那么 $\Gamma$ 可以由（2.4）式描述，其边界由一个偶函数 $f \in C[-1,1] \cap C^1(-1,1)$ 定义，该函数满足以下条件：

  • $f(x_1) = f(-x_1) > 0$ 对于 $x_1 \in (-1,1)$，且 $f(\pm 1) = 0$。
  • $f$ 在 $(-1,1)$ 内是局部解析的。
  • $f'(0) = 0$，对于 $x_1 \in (0,1)$，有 $f'(x_1) < 0$，并且 $\lim_{x_1 \to 1^-} f'(x_1) = -\infty$。

特别地，对 $f'$ 的渐近估计如下：

  • 对于 $x_1 \in [0, 1/2]$，有 $-f'(x_1) \simeq x_1$。
  • 对于 $x_1 \in [1/2, 1)$，有 $-f'(x_1) \simeq \log\left(1 + \frac{1}{1-x_1}\right)$。

其中，$a \simeq b$ 表示存在绝对常数 $C_1, C_2 > 0$ 使得 $C_1 b \le a \le C_2 b$。

基于上述关于 $f'$ 的估计，可以得出 $f$ 的性态：

  • 对于 $x_1 \sim 0$，有 $f(x_1) - f(0) \simeq -x_1^2$。
  • 对于 $x_1 \sim 1$，有 $f(x_1) \simeq (1-x_1)\log\left(1 + \frac{1}{1-x_1}\right)$。

关键结论：

  1. 存在性: 严格证明了具有连续涡斑偶极子结构的二维不可压缩欧拉方程的定常行波解的存在性。
  2. 结构特点: 该解的涡度集中在一个区域 $D$ 内，涡度值在 $D$ 的上半部分为 $+1$，下半部分为 $-1$。该区域 $D$ 关于两个坐标轴对称，且形状不是一个简单的圆盘，而是一个“橄榄形”（olive-shaped）的区域，如 Figure 1.1(a) 所示。
  3. 边界性质: 证明了边界 $\Gamma = \partial D \cap \mathbb{H}$ 由一个偶函数 $f$ 描述。该函数在 $(-1,1)$ 内局部解析，并且在端点 $(\pm 1, 0)$ 处，边界垂直于 $x_1$ 轴，即接触角为 90 度。这个 90 度的接触角与 Sadovskii 的早期数值模拟预测相符。
  4. 与 Lamb-Chaplygin 偶极子的区别: 与圆盘状的 Lamb-Chaplygin 偶极子不同，本文构造的连续涡斑偶极子具有一个非圆形的支撑域，并且涡度是常数（$\pm 1$），而不是像 Lamb-Chaplygin 偶极子那样在边界处连续消失。

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

本论文的核心突破在于其采用了一种新颖的固定点方法（novel fixed-point approach）来构造所寻求的解，这与通常用于构造定常解的变分方法或奇异摄动（bifurcation）方法不同。

a) 核心分析学工具

  1. 固定点定理 (Fixed-Point Theorem):

    • Schauder 固定点定理: 这是证明中使用的主要工具。作者通过将寻找边缘函数 $f$ 的问题转化为一个非线性映射 $\mathbf{R}$ 的不动点问题。为了应用 Schauder 定理，他们首先定义了一个与 $\mathbf{R}$ 相关的连续时间动力系统（continuous-time dynamical system）
      $$ \frac{\partial}{\partial t} f = \mathbf{R}(f) - f $$
      其平衡点即为 $\mathbf{R}$ 的不动点。
    • 他们构造了一个合适的函数空间 $\mathbb{D}$，该空间在 $L^\infty$ 拓扑下是凸的、闭的和紧的。然后证明了解 $f(\mathbf{x}, t)$ 在该动力系统中对于所有 $t \ge 0$ 保持在 $\mathbb{D}$ 中，如果初始状态 $f(\mathbf{x}, 0)$ 也在 $\mathbb{D}$ 中。
    • 最后，他们利用 Schauder 固定点定理证明了该动力系统存在任意正周期的时间周期解。通过将周期趋向于零，他们获得了该动力系统的一个稳态（steady state），这也就是原映射 $\mathbf{R}$ 的一个不动点。

  2. 隐函数定理 (Implicit Function Theorem):

    • 为了定义非线性映射 $\mathbf{R}$，作者首先定义了一个重要的表达式 $F(\mathbf{x}; f)$，它依赖于流函数 $\psi$ 的 $x_2$ 分量及其与速度 $c(f)$ 的关系。
    • 然后，他们利用隐函数定理来定义一个隐式迭代过程 $f \mapsto \tilde{f}$，其中 $\tilde{f}$ 是由方程 $F(x_1, \tilde{f}(x_1); f) = 0$ 所确定的。这个过程定义了映射 $\tilde{f} = \mathbf{R}(f)$。
    • 需要证明 $F$ 具有使得隐函数定理成立的良好性质，例如偏导数的连续性和非奇异性。

  3. 椭圆正则性理论 (Elliptic Regularity Theory):

    • 在处理流函数 $\phi$ 及其导数时，作者明确提到了应用椭圆正则性理论来确保 $\phi \in C^1(\mathbb{H})$。这对证明 $F$ 的连续性以及处理边界附近的奇异性至关重要。
    • 该理论保证了定义在 $\mathbb{H}$ 上的由积分定义的流函数 $\phi$ 及其导数在区域内部和边界上的光滑性。

  4. 非线性积分方程的分析:

    • 流函数 $\phi$ 由一个对数核的积分定义，涉及到所求定义域 $\Omega$ 的边界 $\Gamma$。分析 $\phi$ 及其导数在边界附近的行为是关键。
    • 作者利用了积分表达式，并通过奇点处理（singularity handling）和变量代换来估计积分的行为，尤其是在边界点（如 $x_1 \to \pm 1$）附近。

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

  1. 自由边界问题的转化:

    • 论文最初将寻找连续涡斑偶极子区域 $D$ 的问题，重新表述为寻找描述其上半边界 $\Gamma = \partial D \cap \mathbb{H}$ 的轮廓函数 $f(x_1)$ 的问题。
    • 这个边界 $\Gamma$ 满足流函数 $\psi(\mathbf{x}) = \phi(\mathbf{x}) - cx_2$ 在 $\Gamma$ 上为零的条件，即 $\phi(x_1, f(x_1)) - c f(x_1) = 0$。这将其转化为一个自由边界问题。

  2. 不动点映射的构造与分析:

    • 如上所述，作者的核心技术是将自由边界问题转化为一个非线性映射 $\mathbf{R}$ 的不动点问题。具体来说，他们定义了 $F(\mathbf{x}; f) = \phi_2(\mathbf{x}; f) - c(f)$。寻找解 $f$ 相当于找到 $F(x_1, x_2; f) = 0$ 的轮廓线。
    • 他们仔细分析了 $F$ 在 $x_1$ 轴上的行为（Lemma 3.1），证明了 $F(x_1, 0; f) > 0$ 在 $(-1,1)$ 内，这保证了存在的轮廓线 $f(x_1)$ 满足 $f(\pm 1) = 0$ 且 $f(x_1) > 0$ 在 $(-1,1)$ 内。
    • 他们进一步分析了 $F$ 的偏导数 $\partial_{x_1}F$ 和 $\partial_{x_2}F$ （Lemma 3.2），证明了它们的符号（例如，$\partial_{x_1}F < 0$, $\partial_{x_2}F < 0$ 对于 $x_1, x_2 > 0$），这使得隐函数定理可以用来定义映射 $\mathbf{R}: f \mapsto \tilde{f}$，其中 $\tilde{f}$ 由 $F(x_1, \tilde{f}(x_1); f) = 0$ 确定。
    • 通过分析 $F$ 的相关性质（如单调性），他们成功定义了从 $\mathbf{𝕄}_0$ 到 $\mathbf{𝕄}_1$ 的映射 $\mathbf{R}$。

  3. 控制边界行为的估计:

    • 文中给出了关于轮廓函数 $f$ 的导数 $f'$ 的渐近估计，特别是 $\lim_{x_1 \to 1^-} f'(x_1) = -\infty$。这表明边界在 $(\pm 1, 0)$ 处是垂直的，形成了 90 度的接触角。
    • 这些估计是通过分析 $F$ 的导数和利用积分表示来实现的。例如，Lemma 3.1 的证明中，对 $h'(x)$ 的估计是关键，尽管证明中直接使用了 P.V. 积分和一些不等式，但最终得到了 $h'(x)<0$.
    • Lemma 3.2 中对 $\phi_{x_1}$ 和 $\phi_{x_2}$ 的计算和估计，特别是利用了 $f$ 的单调性和一些积分变换，是验证 $F$ 的单调性以及后续应用隐函数定理的关键。

  4. 构造紧凸集:

    • 为了应用 Schauder 定理，需要一个在 $L^\infty$ 范数下是凸、闭、紧的集合 $\mathbb{D}$，并且该集合是映射 $\mathbf{R}$ 的不变集。
    • 作者通过对 $\mathbf{R}$ 的估计（特别是基于 $F$ 的导数估计），成功地构造了这样一个集合 $\mathbb{D}$。这包括为 $f(x_1, t)$ 找到上界和下界（barriers）。这些边界的选取是基于 $F$ 在边界附近的导数行为。

总而言之，该论文采用了一种基于不动点映射的分析方法，通过将几何的自由边界问题转化为代数方程，再转化为函数空间的映射不动点问题，并最终利用 Schauder 定理证明存在性。这种方法克服了传统方法在处理“连续涡斑偶极子”这种特定几何结构时的困难。

4. 与之前工作的比较

论文在引言部分详细阐述了其研究与先前工作的联系和突破。

a) 论文声称做出的突破和改进

  1. 首次严格构造连续涡斑偶极子解:

    • 突破: 论文首次严格证明了存在具有连续边缘（contiguous boundary）的涡斑偶极子解。这解决了 Sadovskii 早期数值工作提出的一个长期存在的数学问题。
    • 与先前工作对比: 此前虽然有大量关于涡斑解的数值研究（[Sad71], [Pie80], [SS80], [ST82], [WOIZ84], [MST88], [Che93], [FS17], [CG18]），但缺乏严格的数学证明。而之前构造的涡斑解（如 V-state 或其他偶极子）通常不涉及两个涡斑区域完全接触的情况，或者接触区域的边界不是连续的。

  2. 明确的接触角 (90 度):

    • 突破: 论文不仅证明了存在的解，而且精确刻画了在 $(\pm 1, 0)$ 点的边界行为，证明了接触角为 90 度。这是 Sadovskii 数值模拟中观察到的精确现象，现在得到了严格的数学确认。
    • 与先前工作对比: 许多数值研究或尝试性的解析工作可能暗示了这种行为，但缺乏严格证明。

  3. 新颖的数学方法:

    • 改进: 作者采用了一种新颖的固定点方法，而不是常用的变分方法（variational approach）或奇异摄动方法（bifurcation approach）。
    • 与先前工作对比:
      • 变分方法: 例如，构造 Hill 球形涡旋或 Lamb-Chaplygin 偶极子常采用变分框架，将解构造为某种能量泛函的极值点（[Nor75], [Tur83], [Bur88], [Bad98], [HMW20], [GSPSY21]）。这种方法有利于证明动力稳定性（dynamical stability），但可能不利于精确描述边界的几何性质，特别是复杂的接触行为。
      • 奇异摄动方法: 例如，通过扰动简单的定常解来构造复杂解（[Bur82], [HM17], [CL23]）。这种方法适用于解是简单解的微小扰动的情况，但连续涡斑偶极子的构造似乎不属于这类范畴，因为其边界条件（在 $x_1$ 轴处接触）可能难以通过扰动现有闭合解来满足。
    • 作者认为，他们提出的固定点方法非常适合确定涡斑区域的边界，并能精确地导出诸如 90 度接触角这样的几何特征。

  4. 与 Choi-Jeong-Sim (CJS) 的独立工作对比:

    • 论文提到，在他们论文发布的同时，Choi, Jeong, Sim (CJS) 也发表了一项关于 Sadovskii 涡斑解存在性的独立工作 [CJS24]。
    • 不同之处: CJS 的工作采用了完全不同的变分框架。他们的构造方法将 Sadovskii 涡斑解自然地关联到动力稳定性。而本文的固定点构造则更精细地描述了涡斑边界，特别是确认了 90 度接触角。
    • 悬而未决: 由于两者均未证明唯一性，目前尚不清楚他们的构造是否一致。但两项工作都独立地解决了这个长期的数学难题。

b) 论文获得的解的比现有结果强在哪里？

  1. 存在的证明: 最显著的贡献是提供了严格的数学存在性证明，而在此之前只有数值迹象。
  2. 几何特征的刻画:
    • 边界正则性: 论文证明了边界是 $C^1$ 的，并且在 $x_1$ 轴接触点处是垂直的。在 $(-1,1)$ 内部边界是局部解析的。
    • 形状: 证明了支撑域 $D$ 不是一个圆盘，而是一个“橄榄形”的区域，其边界可以通过一个偶函数 $f$ 来描述，并给出了 $f$ 的具体渐近行为（如 $f'(x_1)$ 的行为）。
  3. 构造方法的普适性: 提出的固定点方法本身可能具有更广泛的应用前景，可以用于构造其他具有特定边界行为的涡斑解。
  4. 对 Sadovskii 场景的精确确认: 填补了理论上的空白，确认了 Sadovskii 数值模拟中观察到的许多关键特征。

未能证明的方面: 尽管论文进行了数值实验，并且强烈提示该解在给定归一化条件下（$D\cap\{(x_ 1,0): x_ 1 \in \mathbb{R}\} = \{(x_ 1,0): x_ 1 \in [-1,1]\}$）是唯一的，且支撑域 $D$ 是凸的，但作者承认这些是尚未解决的猜想，留待未来工作。这表明，尽管存在性得到了证明，但该解的唯一性（即使在特定归一化下）和形状（凸性）仍是开放问题。