Vanishing viscosity limit for axisymmetric vortex rings

轴对称涡环的零粘性极限

作者: Thierry Gallay, Vladimir Sverak

类别: math.AP, physics.flu-dyn, 35Q30, 76D05, 76D17, 35B07, 35C20, 35B35

发布日期: 2023-01-03

arXiv ID: 2301.01092v3

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

详细分析

中文标题: 轴对称涡环的零粘性极限

1. 研究对象和背景

本论文研究的是三维不可压缩 Navier-Stokes (NS) 方程在低粘性 $\nu > 0$ 情况下的 Cauchy 问题。其研究对象是 轴对称涡环 的动力学行为，特别是在粘性趋于零 ($\nu \to 0$) 时的极限行为。

NS方程的涡度形式 为：

$$ \partial_t \omega + u \cdot \nabla \omega - \omega \cdot \nabla u = \nu \Delta \omega \quad \text{in } \mathbb{R}^3 \times (0, \infty), \quad (1.1) $$

初始条件为：

$$ \omega|_{t=0} = \omega_0 \quad \text{in } \mathbb{R}^3. \quad (1.2) $$

其中，$\omega(x,t)$ 是流体的 涡度 场，速度 $u(x,t)$ 通过 Biot-Savart 定律与 $\omega$ 相关联：$u(x,t) = (4\pi)^{-1} \int_{\mathbb{R}^3} \omega(y,t) \wedge (x-y) |x-y|^{-3} dy$。

本文关注的特殊初始涡度 $\omega_0$ 是一个 无限薄的涡丝 (infinitely thin vortex filament)，其强度为 $\Gamma$，支持在一个 圆周 C 上。更精确地，$\omega_0 = \Gamma \delta_C$ 表示一个强度为 $\Gamma$ 的电流，集中在一个有向圆周 $C \subset \mathbb{R}^3$ 上。

所研究的物理现象是：在粘性流体中，这种初始的无限薄圆周涡丝如何演变为一个具有有限厚度的 粘性涡环 (viscous vortex ring)。这个涡环由于 自感应 (self-induction) 效应，将沿着其对称轴平移。这种效应最早由 Helmholtz 于1858年提出，并由 Kelvin 于1867年定量化，给出了涡环的平移速度 $V$ 的近似公式 (1.4)。对于粘性情况，涡环的厚度预计正比于 $\sqrt{\nu t}$。 Saffman 在此基础上确定了适用于高斯分布涡量的常数 C，从而得到了 Kelvin-Saffman 公式 (1.5)，用于描述粘性涡环的速度 $V(t)$：

$$ V(t) = \frac{\Gamma}{4\pi r_0} \left( \log \frac{8r_0}{\sqrt{\nu t}} - \frac{3}{2}\log 2 - \frac{1}{2}(1-\gamma_E) \right). \quad (1.5) $$

这里 $r_0$ 是涡环的初始半径，$\gamma_E \approx 0.5772...$ 是欧拉常数。

本论文的目标是为圆形涡丝在弱粘性流体中的 双法线运动 (binormal motion) 提供 严格的数学证明 (rigorous justification)，即证明真实解的涡度场 $\omega(x,t)$ 在长时间区间内能被热方程解 $\omega_{lin}(x-a(t),t)$ 很好地近似，其中 $a(t)$ 是由 Kelvin-Saffman 公式决定的平移量。

在PDE研究的区域和条件方面：

  • 区域: $\mathbb{R}^3$ 全空间。
  • 初始条件: $\omega_0 = \Gamma \delta_C$，其中 $C$ 是一个半径为 $r_0$ 的有向圆。这是一个 奇异的初始条件，涡度集中在零维流形上。
  • 对称性: 由于初始条件是轴对称的，论文将分析限制在 轴对称无旋 (axisymmetric flows with no swirl) 的速度和涡度场。这使得三维问题可以在柱坐标 $(r,\theta,z)$ 下简化为二维问题，涡度只具有 $\theta$ 分量 $\omega_\theta$，速度场 $u$ 只有 $r$ 和 $z$ 分量 $u_r, u_z$。通过 Stokes 流函数 $\psi$，速度分量可以表示为 $u_r = -(1/r)\partial_z\psi$ 和 $u_z = (1/r)\partial_r\psi$。
  • 函数空间: 定理 2.1 建立了在 $C^0((0,\infty), L^1(\Omega) \cap L^\infty(\Omega))$ 空间中的唯一解的存在性，其中 $\Omega = \{(r,z) \in \mathbb{R}^2 | r>0\}$。这表明解在 $t>0$ 时是光滑的，但粘性涡环的厚度 $\sqrt{\nu t}$ 需要与 $r_0$ 相比很小。

2. 主要定理或主要结果

本论文的主要结果体现在 定理 1.1 中，并由 定理 2.6 提供了更详细和更强的版本（虽然此处只给出了定理 1.1）。

定理 1.1 (主要结果的概述): 存在无量纲常数 $K > 0$, $R_0 > 0$ 和 $\sigma \in (0, 1/3)$，对于所有 $\Gamma > 0$, $r_0 > 0$, 以及所有满足雷诺数 $Re := \Gamma/\nu \ge R_0$ 的 $\nu > 0$，如果初始涡度 $\omega_0 = \Gamma \delta_C$（$C$ 是半径为 $r_0$ 的有向圆），那么在 [33] 中建立的唯一的轴对称解 $\omega$ 在时间区间 $t \in (0, T_{adv} Re^\sigma)$ 内可以表示为：

$$ \omega(x,t) = \omega_{lin}(x - a(t), t) + \omega_{cor}(x,t), \quad \text{with } \|\omega_{cor}(\cdot,t)\|_{L^1(\mathbb{R}^3)} \le K\Gamma \left(\frac{\sqrt{\nu t}}{r_0}\right)^{1-3\sigma}. \quad (1.7) $$

其中，$\omega_{lin}(\cdot,t) = \exp(\nu t \Delta)\omega_0$ 是以 $\omega_0$ 为初始数据的热方程解，$a(t)$ 描述了涡环沿其对称轴的平移，其轨迹由 Kelvin-Saffman 公式 (1.5) 给出。具体地，如果 $C = \{(r_0 \cos\theta, r_0 \sin\theta, 0); \theta \in [0, 2\pi]\}$ 沿正向定向，则 $a(t) = (0,0,a_3(t))$，其中 $a_3(t) = \int_0^t V(s)ds$，而 $V(s)$ 由 (1.5) 给出。

定理 1.1 的要点：

  • 近似解的构成: 真实解 $\omega(x,t)$ 被一个主要项 $\omega_{lin}(x-a(t),t)$ 加上一个校正项 $\omega_{cor}(x,t)$ 近似。
  • 主要项的意义: $\omega_{lin}$ 是纯粹的线性扩散（热方程）解，其初始位置在 $C$ 上。但它被平移了一个量 $a(t)$，这个平移速度 $V(t)$ 由 Kelvin-Saffman 公式给出，捕捉了涡环的 自感应 运动。
  • 误差估计: 校正项 $\omega_{cor}$ 的 $L^1$ 范数被精确地控制在一个很小的量级，即 $K\Gamma (\sqrt{\nu t}/r_0)^{1-3\sigma}$。这表明在指定的时间区间内，近似是非常好的。
  • 时间区间: 结果在“大时间区间” $T_{adv} Re^\sigma$ 内有效。这里 $T_{adv} = r_0^2/\Gamma$ 是对流时间，$Re = \Gamma/\nu$ 是雷诺数。由于 $Re \gg 1$，这个时间区间比对流时间本身长得多。
  • 粘性效应: 结果适用于低粘性 $\nu > 0$ 但不为零的情况，通过雷诺数 $Re$ 的限制来体现。
  • 改进: 论文提到，在定理 2.6（未给出详细内容）中，指数 $1-3\sigma$ 可以改进到 $1$，如果考虑更高阶的 Kelvin-Saffman 公式校正。

定理 2.1 (已有的全局适定性结果，作为本工作基础): [33] 中证明，对于任意 $\Gamma > 0$, $\nu > 0$, 和 $(r_0,z_0) \in \Omega$，轴对称涡度方程 (2.3) 存在唯一的全局弱解 $\omega_\theta \in C^0((0,\infty), L^1(\Omega) \cap L^\infty(\Omega))$，使得 $\sup_{t>0} \|\omega_\theta(t)\|_{L^1(\Omega)} < \infty$，并且 $\omega_\theta(t) drdz \rightharpoonup \Gamma \delta_{(r_0,z_0)}$ 当 $t \to 0$ 时。此外，存在常数 $C > 0$ (仅依赖于 $\Gamma/\nu$) 使得：

$$ \left\| \omega_\theta(r,z,t) - \frac{\Gamma}{4\pi\nu t} e^{-\frac{(r-r_0)^2+(z-z_0)^2}{4\nu t}} \right\|_{L^1(\Omega)} \le C\Gamma \log\left(\frac{r_0}{\sqrt{\nu t}}+1\right), $$

对所有 $t \in (0, T_{dif})$ 成立，其中 $T_{dif} = r_0^2/\nu$。 这个结果是本文的起点，但它在低粘性极限下信息有限，因为它将涡环固定在初始位置，并且常数 $C$ 依赖于雷诺数，无法处理 $\nu \to 0$ 的情况。本文的目标正是通过考虑涡环的平移来改进这个估计。

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

本论文的证明策略主要包含两大部分：精确近似解的构造 和 近似解的稳定性分析。

a) 证明中使用的核心分析学工具

  1. 自相似变量 (Self-similar variables)： 为了捕捉解在 $t=0$ 处的奇异行为，论文引入了自相似变量。通过以下变换将涡度 $\omega_\theta$ 和 Stokes 流函数 $\psi$ 转换为无量纲的“廓线”函数 $\eta$ 和 $\varphi$：

    $$ \omega_\theta(r,z,t) = \frac{\Gamma}{\nu t} \eta\left(\frac{r - \bar{r}(t)}{\sqrt{\nu t}}, \frac{z - \bar{z}(t)}{\sqrt{\nu t}}, t\right), \quad (2.13) $$

    $$ \psi(r,z,t) = \Gamma \bar{r}(t) \varphi\left(\frac{r - \bar{r}(t)}{\sqrt{\nu t}}, \frac{z - \bar{z}(t)}{\sqrt{\nu t}}, t\right). \quad (2.13) $$

    其中 $(\bar{r}(t), \bar{z}(t))$ 是涡环的瞬时位置，需要随时间确定。新的自相似变量为 $R = (r-\bar{r}(t))/\sqrt{\nu t}$ 和 $Z = (z-\bar{z}(t))/\sqrt{\nu t}$。这个转换使得涡度方程 (1.1) 转化为一个包含 Fokker-Planck 算子 $\mathcal{L}$ 的演化方程 (2.15)：

    $$ \frac{\Gamma}{\nu} \left( \partial_t \eta + \frac{\epsilon}{2} \mathcal{L}\eta + \frac{\epsilon}{2} (R\partial_R + Z\partial_Z + 1)\eta \right) + \frac{1}{1+\epsilon R} \{ \varphi, \eta \} = \frac{\Gamma}{\nu} \left( \frac{\bar{r}'(t)}{\sqrt{\nu t}} \partial_R \eta + \frac{\bar{z}'(t)}{\sqrt{\nu t}} \partial_Z \eta \right). \quad (\text{简化版}) $$

    以及对应的 Stokes 流函数 $\varphi$ 通过 Biot-Savart 定律 (2.20) 与 $\eta$ 关联。 Fokker-Planck 算子为：

    $$ \mathcal{L} = \partial_R^2 + \partial_Z^2 + R\partial_R + Z\partial_Z + 1. \quad (2.16) $$

    自相似变量的引入使得在 $t=0$ 处的奇性（无限薄涡丝）通过平滑的廓线函数 $\eta$ 得到有效描述。

  2. 微扰展开 (Perturbative expansion)： 近似解的构造是基于对廓线函数 $\eta$ 进行的微扰展开。论文使用了两个小参数：$\epsilon = \sqrt{\nu t}/\bar{r}(t)$（表征涡环厚度与半径之比）和 $\delta = \nu/\Gamma = 1/Re$（粘性与环流强度的比值）。廓线函数 $\eta$ 被展开成一个关于 $\epsilon$ 的渐近级数，直到四阶：

    $$ \eta = \eta_0 + \epsilon\eta_1 + \epsilon^2\eta_2 + \epsilon^3\eta_3 + \epsilon^4\eta_4 + \eta_{cor}. $$

    其中 $\eta_j$ 是通过反演包含 $\delta$ 的算子得到的。这种 双参数展开 区别于文献中的其他方法，允许覆盖 $\epsilon^2$ 和 $\delta$ 不可比的情况，适用于非常小和非常大的时间。此外，$\eta_j$ 还可以依赖于 $\log \epsilon$，涡环速度的领先项正是通过消除 $\eta_1$ 中含 $\log \epsilon$ 的项来确定的。

  3. V. I. Arnold 的几何稳定性方法 (V. I. Arnold’s geometric stability methods)： 这是证明近似解稳定性的 核心工具。Arnold 的方法最初是为 无粘性 ($\nu=0$) 的 Euler 方程 提出的，基于哈密顿系统的几何洞察：

    • Euler 方程可以看作是一个哈密顿系统，哈密顿函数是动能。
    • 定常解是辛叶（与协同伴随轨道一致）上的能量临界点。
    • 如果临界点是轨道上的局部最大值或最小值，则对应的定常态是稳定的。 本论文的创新之处在于，它将 Arnold 的几何稳定性方法 推广 到 轻微粘性 ($\nu > 0$) 的情况。虽然对于 $\nu > 0$，Euler 方程背后的几何结构不再被精确保持，但论文发现 Arnold 稳定性分析中出现的 相关二次型 (relevant quadratic forms) 在比 Arnold 理论中使用的子空间更广的子空间上表现良好，并且与 粘性项 (viscous terms) 有着 有利的相互作用 (interact favorably)。这种“有利相互作用”是论文能够处理粘性情况的关键。

  4. 全局适定性结果 (Global well-posedness result)： 论文以 [33] 中建立的 Cauchy 问题 (1.1), (1.2) 的唯一全局解的存在性为基础 (定理 2.1)。这为后续的精确近似和稳定性分析提供了坚实的数学框架。

b) 解决技术性难题（如非线性项或奇异性）时采用的主要技巧

  1. 处理初始奇异性 (Initial Singularity)：

    • 初始涡度 $\omega_0 = \Gamma \delta_C$ 是一个集中在圆上的 Dirac 测度，是高度奇异的。
    • 通过 自相似变量 变换，这种奇异性被“吸收到”了尺度因子 $\Gamma/(\nu t)$ 中，而廓线函数 $\eta$ 在 $t>0$ 时是光滑的。这使得可以在一个平滑的廓线函数上进行分析。
    • 在自相似变量下，涡度方程中自然出现了扩散项 $\mathcal{L}\eta$，它有效地描述了初始奇点如何随时间扩散并形成具有厚度 $\sqrt{\nu t}$ 的涡环。

  2. 处理非线性项 (Nonlinear Terms)：

    • NS 方程是非线性的，其涡度形式包含非线性项 $u \cdot \nabla \omega - \omega \cdot \nabla u$。
    • 在轴对称无旋的设定下，非线性项在柱坐标下简化为 $\frac{1}{1+\epsilon R} \{ \varphi, \eta \}$，其中 $\{\cdot, \cdot\}$ 是泊松括号。
    • 通过 微扰展开，非线性项被分解为不同阶次的线性或准线性贡献。高阶项的计算和控制是构造精确近似解的关键。
    • Biot-Savart 律 (2.4) 和 (2.20) 被用于将速度 $u$（通过流函数 $\psi$ 或 $\varphi$）与涡度 $\omega$（或 $\eta$）关联起来，从而处理非线性项。这个积分算子本身也是非局部的，增加了复杂性。

  3. 处理线性化算子中的大平流项 (Large Advection Terms in the Linearized Operator)：

    • 稳定性分析需要在线性化算子下进行。在低粘性极限下，线性化算子中可能出现很大的平流项，导致潜在的不稳定性。
    • 论文通过 适应 Arnold 的几何稳定性方法 来解决这个问题。Arnold 方法的核心是寻找哈密顿系统中的守恒量（如能量）在特定约束（如辛叶）下的极值点。这些极值点对应于稳定的解。
    • 即使对于 $\nu > 0$ 这种情况，哈密顿结构不严格保留，论文发现 Arnold 方法中的 二次型 仍然表现良好。这些二次型被证明在比原始理论更广的子空间上是稳定的，并且其行为与粘性项有利地结合，从而能够抑制由大平流项引起的不稳定性。这表明粘性项在某种程度上具有 稳定作用。

  4. 多尺度分析和匹配 (Multi-scale Analysis and Matching)：

    • 微扰展开是一个多尺度问题，参数 $\epsilon$ 和 $\delta$ 具有不同的物理含义，并且可能以不同的方式影响解的行为。
    • 论文通过在自相似变量中进行展开，并精确地匹配不同阶次的方程，从而系统地构造了高阶近似解。

  5. 轴对称性 (Axisymmetry)：

    • 利用初始数据的轴对称性将问题从 $\mathbb{R}^3$ 简化到二维平面 $\Omega = \{(r,z) \in \mathbb{R}^2 | r>0\}$。这大大简化了分析，将三维向量场问题转化为标量 $\omega_\theta$ 和流函数 $\psi$ 的问题。
    • 然而，轴对称性也引入了一些技术挑战，例如在 $r=0$ 轴上的边界条件和奇点。

4. 与之前工作的比较

本论文在多个方面对现有研究做出了显著的突破和改进。

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

  1. 严格证明了弱粘性流体中圆形涡丝的双法线运动: 这是本文最核心的贡献。它为 Helmholtz (1858) 和 Kelvin (1867) 提出的经典物理现象，以及 Saffman (1970) 的定量公式，提供了 严格的数学验证。以往的工作通常是形式化的近似或在特定简化假设下进行的。

  2. 处理了奇异的初始条件 ($\rho_0=0$): 论文的初始条件是 $\omega_0 = \Gamma\delta_C$，即一个 无限薄的涡丝（其“小半径” $\rho_0 = 0$）。这比许多现有工作更具挑战性。例如，Brunelli 和 Marchioro [9] 考虑了初始涡度集中在具有有限小半径 $\rho_0 > 0$ 的环面中的情况。本论文通过精确的自相似变量方法有效地处理了这种 $t=0$ 处的强奇异性。

  3. 在更长的时间尺度上对解进行了精确控制: 论文的结果在时间区间 $t \in (0, T_{adv} Re^\sigma)$ 内有效，其中 $\sigma \in (0, 1/3)$ 且 $Re = \Gamma/\nu \gg 1$。这意味着时间区间可以比对流时间 $T_{adv}$ 长得多。

    • 相比之下，Brunelli 和 Marchioro [9] 的结果仅在时间区间 $T_0 \log(r_0/\rho_0)^{-1}$ 内成立，当 $\rho_0 \to 0$ 时，该时间区间会 缩短至零。本论文的结果在 $\rho_0=0$ 的极限下仍然有效，并且提供了更长的有效时间范围。
    • 对于无粘性情况下的涡环，例如 Dávila, Del Pino, Musso, and Wei [22] 关于“蛙跳”运动的研究，其控制时间尺度为 $T/|\log\epsilon|$，这通常比本论文的 $T_{adv} Re^\sigma$ 要短。本论文通过考虑粘性效应，可能利用了粘性的稳定作用。

  4. 粘性效应没有被视为微扰项: 论文明确指出，“扩散效应没有被视为微扰项”。这意味着粘性项 $\nu \Delta \omega$ 是解动力学演化的一个内在组成部分，而非仅作为一个小扰动来处理。通过自相似变量，扩散效应自然地导致涡环厚度增长为 $\sqrt{\nu t}$，并内嵌在主要近似项中。这与 [9] 等工作中，粘性效应被假定为足够小以进行微扰处理（例如，$\nu T \le \rho_0^2$）形成对比。

  5. 创新性地将 Arnold 几何稳定性方法推广到粘性流体: 这是本论文在理论方法上的一个重大突破。Arnold 的方法最初是为无粘性 Euler 方程设计的，其稳定性分析依赖于哈密顿系统和几何结构（如辛叶）。对于 $\nu > 0$ 的 Navier-Stokes 方程，这些几何结构不再精确保持。论文通过发现 Arnold 方法中相关的 二次型在更大的子空间上表现良好，并与 粘性项有利地相互作用，从而成功地将该方法应用于轻微粘性流体。这提供了一个处理粘性流体稳定性问题的新思路，其概念性解释超越了简单的微扰理论。论文 [34] 在一个更简单的情况下首次证明了这一概念，而本论文将其成功应用于更复杂的涡环问题。

  6. 构造了更精确的近似解: 论文通过在自相似变量中进行到 四阶 的微扰展开，构造了一个非常精确的近似解。这种高阶近似对于实现长时间区间内的精确控制至关重要。

  7. 提供了局部归纳近似猜想在圆形情况下的证明: 局部归纳近似 (Local Induction Approximation, LIA) 猜想是关于涡丝运动的一个重要问题。本论文可以被视为该猜想在 粘性版本下的圆形特例 的严格证明。这为研究更一般的曲线提供了重要的基础和启示。

b) 改进主要体现在放宽了先前工作的哪些假设条件或结果更强

  1. 初始数据正则性要求: 论文放宽了对初始涡度正则性的要求，允许初始涡度是 Dirac 测度（无限薄涡丝），而不是像某些工作中那样需要光滑或支持在有限尺寸区域内的函数。

  2. 解的稳定性: 论文建立了近似解的 强稳定性 特性，即使在存在“线性化算子中非常大的平流项”的情况下也能排除潜在的不稳定性。这在以前的文献中是一个难题，特别是在 $\nu \to 0$ 的极限下，稳定时间区间可能趋于零。

  3. 对扩散效应的处理: 本文的扩散效应被视为问题的一个内在组成部分，而不是作为一个小的微扰。这使得论文能够处理更广泛的粘性参数范围，并更准确地捕捉涡环厚度随时间 $\sqrt{\nu t}$ 的增长。

  4. 对雷诺数的要求: 结果在 大雷诺数 ($Re \gg 1$) 条件下成立，这对应于弱粘性流体。许多早期的粘性流体分析往往限于小雷诺数（强粘性或微扰粘性）情况。

总结来说，本论文通过构造高精度的近似解，并创新性地将 Arnold 的几何稳定性方法推广到粘性流体，为圆形涡环在弱粘性流体中的行为提供了迄今为止最严格和最精确的数学描述。它不仅验证了经典的物理现象，还为处理更复杂的三维涡动力学问题提供了新的分析工具和思路。