Self-similar algebraic spiral vortex sheets of 2-D incompressible Euler equations

二维不可压欧拉方程的自相似代数螺旋涡旋片

作者: Feng Shao, Dongyi Wei, Zhifei Zhang

类别: math.AP

发布日期: 2025-05-06

arXiv ID: 2505.03309v1

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

详细分析

**-algebraic spiral vortex sheet solutions to the 2-D incompressible Euler equations are rigorously constructed for the first time.

中文标题: 二维不可压欧拉方程的自相似代数螺旋涡旋片

1. 研究对象和背景

本文研究的是二维不可压欧拉方程（2-D incompressible Euler equations）在涡旋片（vortex sheet）情形下的自相似代数螺旋解。欧拉方程描述无粘性流体的运动，其速度场 $ \mathbf{u} = (u_1, u_2) $ 和压强 $ P $ 满足如下方程组：

$$ \partial_t \mathbf{u} + \mathbf{u} \cdot \nabla \mathbf{u} + \nabla P = 0, \quad \nabla \cdot \mathbf{u} = 0, $$

在初始条件 $ \mathbf{u}|_{t=0} = \mathbf{u}_0 $ 下求解。通过引入涡度 $ \omega = \partial_1 u_2 - \partial_2 u_1 $ 和流函数 $ \psi $（满足 $ \Delta \psi = \omega $，$ \mathbf{u} = \nabla^\perp \psi = (-\partial_2 \psi, \partial_1 \psi) $），该方程可转化为涡度-流函数形式：

$$ \partial_t \omega + \mathbf{u} \cdot \nabla \omega = 0, \quad \mathbf{u} = K * \omega, $$

其中 $ K(x) = \frac{x^\perp}{2\pi |x|^2} $ 是Biot-Savart核，$ x^\perp = (-x_2, x_1) $。

本文关注的是涡旋片（vortex sheet）——一种涡度集中于一条一维曲线上的奇异解，其数学表达为沿曲线的Dirac测度。涡旋片的演化由Birkhoff-Rott方程描述。对于具有 $ m $-重对称性的涡旋片（即涡度分布在 $ m $ 条旋转对称的曲线上），其Birkhoff-Rott方程为：

$$ \partial_t Z^*(t,\Gamma) = \mathrm{p.v.} \frac{1}{2m\pi i} \sum_{k=0}^{m-1} \int_{-\infty}^{0} \frac{d\Gamma'}{Z(t,\Gamma) - \xi_k^m Z(t,\Gamma')}, $$

其中 $ \xi_k^m = e^{2\pi i k / m} $ 是 $ m $ 次单位根，$ Z(t,\Gamma) $ 是曲线的复参数化，$ \Gamma $ 为总环量（circulation）参数，积分在 $ \Gamma' = \Gamma $ 处取主值（principal value）。初始数据为：

$$ Z(0,\Gamma) = - \left( -\Gamma \right)^{\frac{2\mu - 1}{\mu}}, \quad \mu > \frac{1}{2}, \quad \Gamma < 0. $$

该问题的特殊性在于：所研究的曲线是代数螺旋（algebraic spiral），其极坐标形式为 $ r(\theta) \sim \theta^{-\mu} $（当 $ \theta \to \infty $），这种曲线不是弦弧曲线（chord-arc curve），因此其上的Cauchy积分超出了经典奇异积分算子理论的适用范围。此外，解具有自相似结构：$ Z(t,\Gamma) = t^\mu z(t^{1-2\mu} \Gamma) $，意味着解在时间演化中保持几何形状相似，仅缩放。研究区域为整个平面 $ \mathbb{R}^2 \simeq \mathbb{C} $，边界条件隐含于无穷远行为中，函数空间为加权Hölder空间（见后文）。

2. 主要定理或主要结果

本文的核心结果由两个定理构成：

定理 1.1（自相似代数螺旋解的存在性）
设 $ \mu > \frac{1}{2} $，则存在整数 $ m^* \in \mathbb{Z} $，使得对所有 $ m \in \mathbb{Z} \cap (m^*, +\infty) $，Birkhoff-Rott方程：

$$ \partial_t Z^*(t,\Gamma) = \mathrm{p.v.} \frac{1}{2m\pi i} \sum_{k=0}^{m-1} \int_{-\infty}^{0} \frac{d\Gamma'}{Z(t,\Gamma) - \xi_k^m Z(t,\Gamma')}, $$

在初始数据

$$ Z(0,\Gamma) = - \left( -\Gamma \right)^{\frac{2\mu - 1}{\mu}} $$

下，存在一个自相似解 $ Z_m(t,\Gamma) $，满足：

  • 自相似性：$ Z_m(t,\Gamma) = t^\mu z_m(t^{1-2\mu} \Gamma) $，对所有 $ t > 0 $，$ \Gamma < 0 $；
  • 设 $ z_m(\gamma) = r_m(\gamma) e^{i\theta_m(\gamma)} $，则 $ r_m: (-\infty,0) \to \mathbb{R}^+ $ 和 $ \theta_m: (-\infty,0) \to \mathbb{R} $ 均为 $ C^1 $、严格单调双射；
  • 渐近行为：存在 $ \tilde{\beta}_m > 0 $，使得
    $$ \lim_{\theta \to +\infty} \theta^\mu r_m(\theta_m^{-1}(\theta)) = \tilde{\beta}_m, \quad \lim_{\theta \to +\infty} \theta^{\mu+1} (r_m \circ \theta_m^{-1})'(\theta) = -\mu \tilde{\beta}_m. $$

该解在原点附近呈现代数螺旋（algebraic spiral）行为：$ r \sim \theta^{-\mu} $，这是文献中首次严格构造此类解。

定理 1.2（弱解的验证）
在定理1.1的条件下，定义速度场：

$$ u_m(t,x) = \frac{1}{2m\pi i} \sum_{k=0}^{m-1} \int_{-\infty}^{0} \frac{d\Gamma}{(x - \xi_k^m Z_m(t,\Gamma))^*}, \quad x \notin S_m(t) \cup \{0\}, $$

和涡度测度：

$$ \langle \omega_m(t,\cdot), \varphi \rangle = \frac{1}{m} \sum_{k=0}^{m-1} \int_{-\infty}^{0} \varphi(\xi_k^m Z_m(t,\Gamma)) \, d\Gamma, \quad \forall \varphi \in C_c(\mathbb{R}^2), $$

则 $ (u_m, \omega_m) $ 是二维不可压欧拉方程的弱解，初始涡度为：

$$ \omega_m(0,x) = \frac{m}{|x|^{1/\mu}} \sum_{k=0}^{m-1} \delta_{\xi_k^m \Sigma_0}(x), $$

其中 $ \Sigma_0 = \{ x = (x_1, 0) \in \mathbb{R}^2 : x_1 > 0 \} $ 为初始半无限直线。此外，$ u_m \in C_{\mathrm{loc}}([0,\infty); L^2(\mathbb{R}^2)) $，$ \omega_m \in C_{\mathrm{loc}}([0,\infty); H^{-1}(\mathbb{R}^2)) \cap C_{\mathrm{loc}}([0,\infty); \mathcal{M}(\mathbb{R}^2)) $，且具有自相似缩放：$ u_m(t,x) = t^{\mu-1} v_m(t^{-\mu} x) $，其中 $ v_m \in L^2_{\mathrm{loc}}(\mathbb{R}^2) $。

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

a) 核心分析学工具

  • Banach不动点定理：这是证明的核心。作者将方程转化为不动点形式 $ (r,\gamma) = \mathcal{F}_m[r,\gamma] := \mathcal{N}^{-1} \circ \mathcal{I}_m[r,\gamma] $，并在加权Hölder空间中证明映射 $ \mathcal{F}_m $ 是压缩映射。
  • 加权Hölder空间：为处理代数螺旋在无穷远和原点处的奇异性，作者定义了加权空间：
    $$ X := \left\{ (r,\gamma) : \theta^\mu r \in G_1, \, \theta^{\mu+1} r' \in G_0, \, \theta^{2\mu-1} \gamma \in G_1, \, \theta^{2\mu} \gamma' \in G_1 \right\}, $$
    其中 $ G_k $ 是加权Hölder范数空间：
    $$ \|f\|_k := \|\langle \theta \rangle^\alpha f\|_{L^\infty} + [f]_k, \quad [f]_k := \sup_{0 < \theta_2 < \theta_1 < \theta_2 + 1} \frac{|\langle \theta_1 \rangle^k f(\theta_1) - \langle \theta_2 \rangle^k f(\theta_2)|}{|\theta_1 - \theta_2|^\alpha}. $$
    此空间能有效控制 $ r \sim \theta^{-\mu} $ 和 $ \gamma \sim \theta^{1-2\mu} $ 的渐近行为。
  • Frechet导数与线性化：为证明 $ \mathcal{N} $ 的可逆性，作者对非线性算子 $ \mathcal{N} $ 在Kaden解 $ (r_0,\gamma_0) $ 处进行线性化，得到线性算子 $ \mathcal{L} $，并显式构造其逆算子（见引理4.7），从而证明 $ \mathcal{N} $ 局部可逆。
  • Newton-Kantorovich框架：虽然未直接使用完整Newton方法，但其思想（非线性算子的局部可逆性 + 小扰动）贯穿始终。

b) 关键技术与奇异性处理技巧

  • 将Birkhoff-Rott方程分解为局部算子与积分算子：作者将方程重写为 $ \mathcal{N}[r,\gamma] = \mathcal{I}_m[r,\gamma] $，其中 $ \mathcal{N} $ 是局部非线性微分算子（与 $ m $ 无关），而 $ \mathcal{I}_m $ 是包含主值积分的全局算子（依赖 $ m $）。这一分解使得可将 $ \mathcal{I}_m $ 视为对 $ \mathcal{N} $ 的小扰动。
  • Kaden螺旋作为极限解：当 $ m \to \infty $ 时，方程（2.10）形式收敛到Kaden代数螺旋解（2.13）。作者利用这一极限解作为“参考解”，并证明当 $ m $ 足够大时，$ \mathcal{I}_m $ 可被控制为小量，从而通过不动点定理获得解。
  • Cauchy积分的精细估计：这是最核心的技术难点。由于Kaden螺旋 $ r_0(\theta) = \theta^{-\mu} $ 不是弦弧曲线，经典奇异积分理论（如Calderón-Zygmund）失效。作者将积分 $ \mathcal{I}_m $ 分解为振荡积分之和，并采用加权 $ W^{1,\infty} $ 估计逐项控制其范数（见命题2.4），这是对奇异积分在非光滑曲线上分析的创新。
  • 对称性利用：利用 $ m $-重旋转对称性，将 $ m $ 个分支的积分合并为一个积分，大大简化了计算结构。

4. 与之前工作的比较

a) 突破与改进

本文在多个方面实现了对先前工作的实质性突破：

  • 首次构造代数螺旋涡旋片解：此前虽有大量数值模拟（如Krasny [31]、Pullin [39,40]）表明涡旋片在奇点后形成代数螺旋结构，但从未有严格数学证明。本文首次给出了此类解的存在性与定量渐近行为的完整构造。
  • 超越对数螺旋的局限性：Prandtl、Alexander、Kambe等人曾构造过对数螺旋解 $ r \sim e^{-a\theta} $，但Wu [50] 证明对数螺旋是弦弧曲线，无法描述涡旋片在奇点后的快速卷曲行为。本文构造的代数螺旋 $ r \sim \theta^{-\mu} $ 不是弦弧曲线，从而符合物理预期，是更合理的数学模型。
  • 放宽初始数据正则性要求：文献中如[44]虽在初始涡度为 $ \omega(0,x) = r^{-1/\mu} \dot{\omega}(\theta) $ 时证明了弱解存在性，但未给出定量渐近行为。本文在更强的对称性假设（$ m $-重，$ m $ 足够大）下，不仅证明存在性，还精确刻画了 $ r(\theta) \sim \theta^{-\mu} $ 和 $ \gamma(\theta) \sim \theta^{1-2\mu} $ 的渐近常数 $ \tilde{\beta}_m $，这是前所未有的定量结果。
  • 处理非弦弧曲线上的奇异积分：经典理论要求曲线满足Lipschitz或弦弧条件，而本文首次在非弦弧、代数发散的曲线上完成了Cauchy积分的严格估计，为研究更一般奇异涡结构开辟了新途径。
  • 建立在加权Hölder空间中的全局分析框架：相比以往依赖解析函数或Sobolev空间的方法，本文在加权Hölder空间中构建了全新的分析框架，适用于具有幂律奇异性的问题。

综上，本文不仅解决了长期悬而未决的“涡旋片卷曲后形态”的数学验证问题，更在分析工具上实现了从经典奇异积分理论到非光滑奇异曲线分析的范式跃迁。