Convergence of a two-parameter hyperbolic relaxation system toward the incompressible Navier-Stokes equations

技术摘要：双参数双曲松弛系统向不可压缩纳维-斯托克斯方程的收敛性

问题陈述
本文研究了在二维（$T^2$）和三维（$T^3$）周期域内，双参数双曲松弛系统 (1.2) 向不可压缩纳维-斯托克斯方程 (1.1) 的严格收敛性。该近似系统 (1.2) 整合了两种经典的范式：人工可压缩性（AC）方法（通过压力的时间演化方程来松弛散度自由约束）以及针对类应力张量 $U$ 的一阶松弛模型。该系统由小参数 $\epsilon$ 和 $\delta$ 控制，分别代表压力和应力张量的松弛时间。

本文解决的主要数学挑战是同时进行极限过程 $(\epsilon, \delta) \to 0$。虽然单参数极限（即 $\epsilon \to 0$ 的 AC 方法或 $\delta \to 0$ 的松弛模型）已有充分文献记载，但同时极限的讨论存在困难，特别是在控制压力误差和处理初始数据扰动方面。具体而言，作者旨在针对两种不同的初始速度扰动情形建立收敛性：

1. 消失扰动（Vanishing perturbations）： 初始差异 $\|u_{\epsilon,\delta}^0 - u^0\|$ 的阶数为 $O((\epsilon+\delta)^\alpha)$。
2. 阶数为 $O(1)$ 的扰动（Order-one perturbations）： 初始差异为 $O(1)$，这一情形对于数值鲁棒性至关重要，但由于松弛系统中的耗散有限，对于基于能量的方法而言传统上难以处理。

方法论
作者根据初始扰动的量级以及空间的维度采用了不同的分析策略。

• 中间仿射系统（针对消失扰动）： 为了在消失扰动情形下推导压力的收敛速率，作者构造了一个中间仿射系统。该系统作为原始松弛系统 (1.2) 与目标纳维-斯托克斯方程 (1.1) 之间的桥梁。通过界定中间系统与极限系统及近似系统之间的差异，作者获得了必要的压力误差估计，从而克服了直接能量法通常只能获得 $O(1)$ 压力误差界的局限性。

• 调制能量与 Bootstrap 参数（针对 $O(1)$ 扰动）： 对于 $O(1)$ 扰动情形，标准的对称双曲能量估计是不够的，因为由张量 $U_{\epsilon,\delta}$ 的阻尼提供的耗散不足以在长时间间隔内抑制非线性效应。相反，作者利用了一种涉及涡度 $\omega_{\epsilon,\delta} = \text{curl } u_{\epsilon,\delta}$ 的调制能量泛函（引入自文献 [4]）。

  • 该方法利用 $\nabla \omega_{\epsilon,\delta}$ 的耗散来控制速度场。
  • 通过 Bootstrap 参数法 来闭合估计。这需要特定的参数缩放关系（特别是 $\epsilon = o(\delta)$）以及对初始数据的约束，即初始速度的散度相对于初始涡度需足够小。
  • 初始层修正（Initial Layer Correction）： 为了处理类应力张量 $U_{\epsilon,\delta}$ 而无需要求其在初始时刻与极限初始数据满足 $o(1)$ 接近，作者在能量估计中引入了初始层修正（指数衰减项）。

• 维度的区别：

  • 在 3D ($T^3$) 中，涡度是一个矢量场。分析识别出了基于初始速度扰动 $L^2$ 范数阶数（$O(\delta^{3/4})$ 与 $O(\delta^{1/4})$）的两种截然不同的“极端情况”，从而导出了定理 2.6 和 2.7。
  • 在 2D ($T^2$) 中，涡度是一个标量。分析显示，对于初始扰动的收敛阶仅存在一种可行的构造方式（定理 2.8）。

主要结果
本文建立了以下严格的收敛结果：

1. 消失扰动情形（定理 2.1, 2.5, 推论 2.3）：

   • 对于足够小的初始扰动，松弛系统的解在与光滑纳维-斯托克斯解相同的时区间内存在（若纳维-斯托克斯解是全局解，则为对数扩展区间）。
   • 证明了速度 $u_{\epsilon,\delta} \to u$ 在 $L^\infty([0,T]; H^1)$ 和 $L^\infty([0,T]; L^2)$ 中收敛。
   • 至关重要的是，在更严格的条件（特别是 $\delta = o(\sqrt{\epsilon})$）下，证明了压力 $p_{\epsilon,\delta}$ 在 $L^\infty([0,T]; H^1)$ 中收敛于 $p$。

2. $O(1)$ 扰动情形（定理 2.6, 2.7, 2.8）：

   • 作者证明了即使在初始速度扰动为 $O(1)$ 时，仍能实现全局时间（或直到极限解的寿命）内的速度收敛。
   • 3D 结果： 两个定理根据初始数据的不同缩放假设，提供了速度在 $L^2$ 和涡度在 $H^1$ 中的收敛界。这些界取决于 $\delta$ 以及比例 $\epsilon/\delta$。
   • 2D 结果： 在特定的初始扰动缩放条件下，也建立了类似的 2D 全局速度收敛性。
   • 在所有 $O(1)$ 情形中，速度场的散度 $\text{div } u_{\epsilon,\delta}$ 被证明保持在较小水平（受 $\epsilon$ 和 $\delta$ 控制），确保了解保持在不可压缩流形附近。

意义与主张
本文声称为双参数双曲松弛模型 (1.2) 在数值模拟中的应用提供了基础性的理论验证。其具体贡献包括：

• 双参数同时极限： 它填补了现有文献中关于同时极限 $(\epsilon, \delta) \to 0$ 的空白，这一场景相比单参数极限研究较少。
• 压力收敛性： 通过引入中间仿射系统，本研究严格推导了消失扰动情形下的压力误差估计，这是以往松弛模型分析中经常缺失的一步。
• 对大扰动的鲁棒性： 将速度收敛扩展到 $O(1)$ 扰动情形，展示了该模型应对显著初始不一致性的鲁棒性，这种不一致在数值启动阶段很常见。这是通过调制能量法和 Bootstrap 参数法克服标准能量估计的局限性来实现的。
• 结构性洞察： 分析强调了 3D 情况相对于 2D 情况在结构上的丰富性，识别出了三维空间中不同的收敛速率区间。

作者指出，虽然在 $O(1)$ 扰动下建立了速度收敛性，但在该特定情形下压力的收敛性仍是未来研究的一个开放性挑战。本研究旨在为该模型在下游任务（如湍流统计）中的应用提供数学上的有效性支持。