改进的SQG贴片方程的有限时间奇点形成

Finite time singularity formation for the modified SQG patch equation

作者: Alexander Kiselev, Lenya Ryzhik, Yao Yao, Andrej Zlatos

类别: math.AP

发布日期: 2015-08-30

arXiv ID: 1508.07613v1

链接: http://arxiv.org/abs/1508.07613v1

详细分析

中文标题: 改进的SQG贴片方程的有限时间奇点形成

1. 研究对象和背景

本文研究的是二维改进的表面准地转（SQG）贴片方程在半平面上的动力学行为，特别是其解的全局正则性与有限时间奇点形成的问题。

背景: 二维不可压缩欧拉方程以其全局正则解而闻名。然而，SQG方程，作为一种常用于大气科学模型的方程，其Biot-Savart定律的正则性比欧拉方程低一个导数，因此其解的全局正则性问题仍是未解决的挑战。

本文考虑了一族活跃标量方程，它们通过一个参数 $\alpha \in [0, 1/2]$ 插值于二维欧拉方程（对应 $\alpha=0$）和SQG方程（对应 $\alpha=1/2$）之间。这族方程的Biot-Savart定律内核的奇异性程度由 $\alpha$ 控制。

方程组可以表述如下。在半平面 $D := \mathbb{R} \times \mathbb{R}^+$ 上，速度场 $u$ 由以下Biot-Savart定律给出：

其中 $\omega$ 是活跃标量，$\nabla^\perp = (-\partial_{x_2}, \partial_{x_1})$ 是旋度算子。上述公式在有界区域内通常通过傅里叶变换或格林函数定义。论文中，对于在半平面上的狄利克雷拉普拉斯算子，作者给出了更明确的积分形式：

其中 $x = (x_1, x_2)$，$y = (y_1, y_2)$，$\bar{y} = (y_1, -y_2)$ 是 $y$ 关于 $x_1$ 轴的镜像。向量场 $u$ 满足散度自由条件，并且切向于边界，即 $u_2(x,t) = 0$ 当 $x_2 = 0$。

活跃标量 $\omega$ 的演化遵循对流方程：

这与二维欧拉方程的涡度形式相同。

文章特别关注涡旋贴片（vortex patches）解。涡旋贴片是指活跃标量 $\omega$ 是若干具有光滑边界的有界开集 $\Omega_k$ 的指示函数的线性组合，即 $ \omega(x,t) = \sum_{k=1}^N \theta_k \chi_{\Omega_k(t)}(x) $ 其中 $\theta_k$ 是常数，$\Omega_k(t)$ 是随时间演化的区域，$\chi_{\Omega_k(t)}$ 是 $\Omega_k(t)$ 的指示函数。涡旋贴片模型模拟了流体中涡度有剧烈变化的流况，这在自然界中很常见。

在研究的区域方面，文章考虑的是半平面 $D := \mathbb{R} \times \mathbb{R}^+$，这增加了处理边界的复杂性。对于二维欧拉方程，虽然泊松方程和涡度方程在整个平面上或周期性区域上已得到充分研究，但在有边界区域（如半平面）上，特别是当贴片初始就接触边界时，全局正则性的证明更加微妙。

函数空间方面，文章同时讨论了 $C^{1,\gamma}$ 类型的贴片解（对应二维欧拉方程）和 $H^3$ 类型的贴片解（对应改进的SQG方程）。$H^3$ 范数比 $C^{1,\gamma}$ 范数要求更高的正则性，并且是局部存在性定理所使用的函数空间。

研究的现象: 文章的核心在于研究有限时间奇点形成的现象。对于二维欧拉方程（$\alpha=0$），已知其贴片解可以保持全局光滑。然而，对于 $\alpha > 0$ 的改进SQG贴片方程，虽然比欧拉方程更“奇异”，但作者发现存在初始数据会导致解在有限时间内失去光滑性，或者贴片边界发生自相交，或者不同贴片互相接触。这表明在 $\alpha$ 穿越0点时，解的行为发生了一个相变。

2. 主要定理或主要结果

文章提出了两个主要结果，通过对比 $\alpha=0$ 和 $\alpha>0$ 的情况，阐述了这种相变：

Theorem 1.2 (全局正则性 for the 2D Euler case $\alpha=0$). 假设 $\alpha = 0$ 且 $\gamma \in (0,1]$。对于任意 $C^{1,\gamma}$ 类型的贴片初始数据 $\omega_0$，存在一个唯一的全局 $C^{1,\gamma}$ 贴片解 $\omega$ 满足方程组 (1.1)-(1.2)。这意味着在二维欧拉方程的情况下，即使贴片初始接触半平面的边界，解也能保持全局光滑。

Theorem 1.3 (有限时间奇点形成 for small $\alpha > 0$). 假设 $\alpha \in (0, 1/4)$。则存在 $H^3$ 类型的贴片初始数据 $\omega_0$，使得方程组 (1.1)-(1.2) 的唯一局部 $H^3$ 贴片解 $\omega$ 在有限时间 $T < \infty$ 内形成奇点。奇点的形成表现为：解的最大存在时间 $T_\omega$ 是有限的。根据Theorem 1.1 (来自伴随论文 [36])，有限最大存在时间的发生意味着奇点形成：要么两个贴片接触，要么一个贴片边界自触，或者贴片边界失去 $H^3$ 正则性。

这两个定理的结合表明，在 $\alpha=0$ 时存在全局正则解，而对于任意小的 $\alpha>0$，都存在导致有限时间奇点形成的初始数据。这有力地支持了将二维欧拉方程视为“临界（critical）”情况的观点。

文章还给出了导致奇点形成的初始数据示意图 (Figure 1)：包含两个符号相反的贴片，关于 $x_1$ 轴对称，并且初始接触 $x_1$ 轴。作者指出，这种几何构型，尤其是涉及边界上的双曲型不动点，与一些其他研究中观察到的奇点形成场景具有相似性。

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

这篇文章的研究方法结合了多种分析工具和技巧，以处理改进的SQG贴片方程在半平面上的动力学。

a) 核心分析学工具:

• 流动映射 (Flow Map): 类似于二维欧拉方程的经典理论（Yudovich理论），文章考虑了由速度场 $u$ 定义的流动映射 $\Phi_t$，并研究了活跃标量 $\omega$ 如何在该映射下被平移（$\omega(x,t) = \omega_0(\Phi_t^{-1}(x), 0)$）。然而，对于 $\alpha > 0$，速度场 $u$ 在贴片边界上仅有Hölder连续性，这导致流动映射可能不再是唯一的，因此文章采用了替代定义（Definition 2.2）。
• Hölder 不等式与 Sobolev 范数: 文章频繁使用 Hölder 范数 $\| \cdot \|_{C^{1,\gamma}}$ 和 Sobolev 范数 $\| \cdot \|_{H^3}$ 来量化解的正则性。在全局正则性证明（Theorem 1.2）中，关键在于控制这些范数随时间的增长。在奇点形成证明（Theorem 1.3）中，则需要展示这些范数会在有限时间内变得无穷大。
• 能量估计 (Energy Estimates): 虽然文中没有直接给出“能量”的显式表达式，但其证明大量依赖于对解的各种范数（特别是与边界曲率、梯度相关的范数）的时间导数进行估计。例如，在证明 Theorem 1.2 时，通过估计 $\nabla w$ 的 $C^\gamma$ 范数（其中 $w = \nabla^\perp \phi$）的时间导数，得到了 $\mathcal{A}_\gamma(t)$（与边界曲率相关）增长速度的界限。
  $$ \mathcal{A}_\gamma'(t) \leq C \mathcal{A}_\gamma(t) (1+\log \mathcal{A}_\gamma(t)) $$
  这种对数型增长表明解是全局正则的。
• 函数空间理论: 对于半平面上的算子，需要用到狄利克雷拉普拉斯算子的性质。此外，文章还借用了 Whitney 延拓定理（[48, Sec 6.2, Theorem 4]）将贴片内的函数延拓到整个平面，以便应用平面上的估计。
• Sobolev嵌入定理: 文章隐式地利用了 Sobolev 嵌入定理，例如 $H^3$ 嵌入到 $C^{1,\gamma}$。
• 几何分析: 对于贴片边界的几何形状（如曲线的曲率），文章使用了如 $\mathcal{A}_\gamma(t)$ 这样的量进行度量。对曲率增长的控制是确保边界保持光滑的关键。
• 复分析技巧: 在对速度场 $u$ 进行估计时，作者使用了复数表示，例如 $\bar{y}$ 表示镜像点，$\perp$ 符号表示旋转。

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

• 不同 $\alpha$ 值下的行为差异: 研究的核心在于对比 $\alpha = 0$（欧拉）和 $\alpha > 0$（改进SQG）。对于 $\alpha = 0$，速度场的对数Lipschitz性质保证了流动映射的唯一性和解的正则性。而对于 $\alpha > 0$，Biot-Savart算子的奇异性增强，导致速度场的正则性下降，从而引发了在有限时间形成奇点的可能性。
• 处理半平面的边界效应: 在半平面上，需要考虑算子关于边界的对称性以及边界对流速的影响。论文中的 Biot-Savart 核 $K_{ij}$ (4.10)-(4.11) 包含了复杂的项，反映了这种边界效应。处理这些项需要仔细的积分估计和几何分析。
• 贴片边界的正则性控制:
  • 全局正则性证明 (Theorem 1.2): 关键在于证明$\mathcal{A}_\gamma(t)$, $\mathcal{A}_\infty(t)$, $\mathcal{A}_{\inf}(t)$ (与曲率、梯度和法线速度相关的量) 的增长是可控的，不会在有限时间内达到无穷大。这主要通过对$w = \nabla^\perp \phi$（其中 $\phi$ 是定义贴片的函数）的 $C^\gamma$ 范数的时间导数进行估计来实现。尤其需要分拆 $(u \cdot \nabla)w$ 并分别处理其项，利用 $(\nabla v)w$ 和 $(\nabla v^\perp)w^\perp$ 的不同性质。
    • 关键技术点在于：
      1. 将 $(u \cdot \nabla)w$ 分解为 $(\nabla v) w$ 和 $(\nabla v^\perp)w^\perp$ 两部分，其中 $v$ 由一个贴片产生，而 $v^\perp$ 由其镜像产生。
      2. 利用 Bocaccio-Constantin [1] 的结果，对 $(\nabla v) w$ 部分得到控制。
      3. 设计一个对称的 $w^\perp$（记为 $w^\perp$)，使得 $(\nabla v^\perp)w^\perp$ 也能得到类似控制。
      4. 通过估计 $(\nabla v^\perp)w^\perp$ 的 $C^\gamma$ 范数，导出 $\mathcal{A}_\gamma'(t)$ 的增长上界，从而证明了全局正则性。
  • 有限时间奇点证明 (Theorem 1.3): 这是文章的主要贡献。作者构造了一个特殊的初始数据，包含两个贴近边界且对称的贴片。核心思想是证明在 $\alpha > 0$ 时，这些贴片有加速向中间（中心线）靠拢的趋势，最终导致在有限时间内接触。
    • 关键技术点在于：
      1. 构造特定的初始数据: chosen to be two patches symmetric with respect to the $x_1$-axis and touching the x-axis (Figure 1).
      2. 估计速度场 $u_1$ (切向速度) 的负值: 在关键区域，证明 $u_1$ 的值是显著负的（Lemma 4.3 and Proposition 4.5）。论文通过将 $u_1$ 分解为“坏”部分（$u_{bad}^1$）和“好”部分（$u_{good}^1$），并分别估计它们，特别是 $u_{bad}^1$ 在特定区域内的负值上界。
      3. 引入时间演化尺度 $X(t)$: 定义了一个时间相关的尺度 $X(t) \propto (T-t)^{1/(1-2\alpha)}$ (Equation 4.3 and 4.4)，这个尺度使得 $X'(t) = -X(t)^{1-2\alpha}$。这个演化规律被用来追踪贴片边界的关键点。
      4. 比较原理: 论文证明了一个比较原理类型的结论：如果初始贴片 $K(0)$ 满足特定包含关系，那么在所有时间 $t \in [0, T)$ 内，$K(t) \subset \tilde{\Omega}(t)$（其中 $K(t)$ 是初始 $K(0)$ 的演化，$\tilde{\Omega}(t)$ 是一个贴片的演化，它们关于 $x_1$ 轴对称）。
      5. 通过 $u_1$ 的负值来证明 $K(t) \subset \tilde{\Omega}(t)$: 证明的核心是利用Proposition 4.5 关于 $u_1$ 的估计。当 $u_1$ 在关键区域是负的时，速度将空间点向左（负 $x_1$ 方向）推移。结合 $X(t)$ 的演化，表明 $K(t)$ 会被“吸入”到 $\tilde{\Omega}(t)$ 内部。
      6. 导致奇点: 当 $K(t) \subset \tilde{\Omega}(t)$ 持续到 $t=T$ 时，意味着 $\Omega(T)$ 和 $\tilde{\Omega}(T)$（贴片及其镜像）在原点接触，这构成了有限时间奇点。
• 对 $\alpha$ 趋于0时的行为分析: 论文的奇点形成论证依赖于 $\alpha>0$。引言中提到，当 $\alpha \to 0$ 时，奇点形成的参数（如 $1/(1-2\alpha)$）会趋于无穷大，这与欧拉方程全局正则性的结果相符。

4. 与之前工作的比较

这篇文章在改进SQG方程和二维欧拉方程的正则性研究领域做出了重要的理论贡献。

a) 文章声称的突破或改进:

• 首个严格的有限时间奇点形成证明: Theorem 1.3 是首个对改进SQG贴片方程（$\alpha>0$）在半平面上，对于任意小的 $\alpha>0$，证明了有限时间奇点形成的严格数学证明。在此之前，虽然有数值模拟（[16]）和一些迹象表明可能存在奇点，但缺乏严格的解析证据。
• 二维欧拉方程在半平面的全局正则性: Theorem 1.2 提供了二维欧拉方程贴片模型在半平面上的全局 $C^{1,\gamma}$ 正则性证明，即使贴片初始接触边界。这填补了一个先前工作（[20]）中可能存在的关于接触边界情况下的时间限制。
• 揭示了 $\alpha=0$ 的临界性: 将这两个结果结合起来，文章明确展示了在 $\alpha=0$ (欧拉方程) 和 $\alpha>0$ (改进SQG方程) 之间存在一个行为上的“相变”。欧拉方程的全局正则性在 $\alpha>0$ 时被打破，出现了有限时间奇点。这为将二维欧拉方程视为“临界”情况提供了坚实的数学基础。

b) 改进的条件和更强的结果:

• 更弱的奇异性（$\alpha > 0$）下的奇点形成: 相较于 SQG 方程 ($\alpha=1/2$ 或更大奇异性，如 $\alpha \in (0, 1]$ 相关研究)，本文关注的是**“轻微奇异”**的改进 SQG 方程 ($\alpha \in (0, 1/4)$)。并且，通过构造特定的、涉及边界接触的初始数据，作者证明了即使在这样的轻微奇异情况下，奇点仍然会形成。这表明奇点形成的条件比先前预期的要宽松。
• 对无散度流（inviscid flow）的理解: 该研究深化了对不可压缩无粘流体（如欧拉方程和SQG方程）中，即使是具有相对简单结构的涡旋贴片在有边界情况下的复杂动力学行为的理解。
• 正则性对比: Theorem 1.2 的全局正则性结果与 Theorem 1.3 的有限时间奇点结果之间形成鲜明对比，突显了 $\alpha$ 参数的微小扰动（从0到正数）对解的长期行为的巨大影响。

与先前工作的对比:

• Yudovich 理论: 对于二维欧拉方程，Yudovich 的理论（[53]）在整个平面上确立了涡旋贴片解的存在性和唯一性。Chemin [7] 进一步证明了其边界保持光滑。此文的 Theorem 1.2 扩展了这些结果到半平面，并且对初始数据接触边界的情况也证明了全局正则性。
• SQG 及改进SQG方程: 对于标准的SQG方程（$\alpha \approx 1/2$），奇点形成问题（[10]）已经被研究，但存在性结果（如[14, 15]）排除了某些特定类型的奇点。作者指出，对SQG和改进SQG方程（$\alpha \in (0, 1]$）光滑初值的问题，即使其比欧拉方程更奇异，但奇点形成的例子是未知的。相反，文章通过贴片模型，找到了一个明确的导致有限时间奇点的构造。
• Caffarelli-Vasseur 和 Kiselev-Nazarov-Volberg 的工作: 这些工作（[3, 34]）证明了临界粘性SQG方程（$\beta=1/2$ 耗散项）的全局正则性。本文研究的是无耗散的情况，并且参数 $\alpha$ 控制的是 Biot-Savart 核的奇异性，而非耗散。
• 伴随工作 [36]: 本文严重依赖于其伴随论文 [36] 中关于局部存在性和唯一性的结果 (Theorem 1.1)，它表明对于 $\alpha \in (0, 1/4)$，存在唯一的局部 $H^3$ 贴片解，并且如果其存在时间有限，则必然形成奇点。本文 Theorem 1.3 将这一局部结果转化为全局有限时间奇点形成的存在性结果。

总的来说，这篇论文通过分析一个参数化的方程族，成功地在一个临界参数值 $\alpha=0$ 处揭示了从全局正则性到有限时间奇点形成的分水岭，为理解流体动力学方程的正则性提供了深刻的见解。