Viscosity in error upper bound for a consistent splitting scheme of the Navier-Stokes equations

技术摘要：关于一致性分裂法在纳维-斯托克斯方程误差上界中的粘性问题

问题陈述
本文研究了由 Huang 和 Shen [5] 提出的二阶一致性分裂方案在纳维-斯托克斯方程中的鲁棒性。虽然原文献 [5] 中的分析确立了无条件稳定性和误差界，但其依赖于粘性 $\nu$ 固定为单位值的假设。本文探讨的核心问题是：当 $\nu \to 0$（高雷诺数）时，该方案是否仍然保持鲁棒。根据 [2] 中的定义，若速度的误差上界不包含粘性的负幂次，则称该方案是鲁棒的。作者旨在通过在整个推导过程中将 $\nu$ 作为符号参数保留，重新分析误差界，以确定 Huang–Shen 方案是否满足这一标准。

方法论
作者通过系统地追踪所有常数对粘性 $\nu$ 的依赖关系，扩展了 Huang 和 Shen [5] 的证明方法。关键的方法步骤包括：

1. 精化 Stokes 压力估计： 误差分析的一个关键组成部分是 Stokes 压力估计。作者精化了来自 [8] 的一个定理（在 [5] 中被引用为引理 4），以显式量化不等式 $\|\nabla p_s(u)\|^2 \leq (\frac{1}{2} + \varepsilon)\|\Delta u\|^2 + C_S(\varepsilon)\|\nabla u\|^2$ 中的常数 $C_S(\varepsilon)$。他们证明了当 $\varepsilon \to 0^+$ 时，$C_S(\varepsilon)$ 的行为表现为 $O(\varepsilon^{-3})$，这一依赖关系对于后续的误差分析至关重要。
2. 广义稳定性分析： 将随时间变化的 Stokes 方程的稳定性分析扩展到任意 $\nu > 0$。这涉及重新推导方案 (3a)–(3e) 的能量估计，并确定实现无条件稳定性所需的最小整数值稳定参数 $k$。
3. 广义误差分析： 将全纳维-斯托克斯方程的误差分析推广到任意 $\nu$。这需要仔细处理非线性对流项，该项在方案中被显式处理。作者采用了离散 Gronwall 不等式，并通过三次嵌套应用的 Gronwall 引理来追踪依赖于 $\nu$ 的常数的累积。
4. 数值验证： 为了证实分析结果，作者在扰动 Kovasznay 流上进行了数值实验。他们将 GSAV-谱方法与全隐式 FEM-Newton 求解器（作为鲁棒的参考基准）以及该方案对应的随时间变化的 Stokes 对等方案进行了对比。

主要贡献与结果

• 降低参数 $k$ 的要求： 作者证明，参数 $k$ 的稳定性和收敛性要求可以从 [5] 中所述的 $k \geq 5$ 降低到 $k \geq 4$。这种降低改善了方案的截断误差，这一点在示例 1 中得到了数值验证。
• 方案的非鲁棒性： 主要的理论结果是，速度在 $H^1$ 范数下的误差上界依赖于粘性的负幂次。具体而言，误差界 $\|\nabla e^{n+1}\| \leq C_{5.3} \delta t^2$ 中的常数 $C_{5.3}$ 随 $\nu \to 0^+$ 呈现出相对于 $\nu^{-1}$ 的三重指数增长。
  • 常数 $C_1$（与速度稳定性相关）的缩放比例为 $O(\nu^{-3} \exp(\nu^{-4}))$。
  • 常数 $C_2$（与速度误差相关）的缩放比例为 $O(\exp(\exp(\nu^{-4})))$。
  • 最终常数 $C_5$（与 GSAV 一致性相关）的缩放比例为 $O(\exp(\exp(\exp(\nu^{-4})))))$。
    因此，为了保证收敛，所需的时间步长 $\delta t$ 必须至少以 $O(\exp(-\exp(\exp(\nu^{-4})))))$ 的速度缩小，这使得该方案在小粘性下是非鲁棒的。
• 失效定位： 数值实验证实，在高雷诺数（$Re \geq 2 \times 10^4$）下，该方案会发生“爆炸”（中间速度和压力发散）。关键在于，当移除对流项（Stokes 问题）时，该方案在所有测试的雷诺数下均保持稳定。同样，全隐式 Newton 求解器也不会出现爆炸现象。这种对比将非鲁棒性的来源定位在分裂方案中对对流项的显式处理上。

意义与主张
本文声称提供了第一个通过显式追踪粘性依赖关系来对二阶一致性分裂方案进行严格稳定性与误差分析的研究。作者得出结论：虽然 Huang–Shen 方案在弱意义上提供了无条件稳定性，并且在 $k=4$ 时具有更高的精度，但对于高雷诺数下的纳维-斯托克斯方程，它在本质上是非鲁棒的。这种失效源于对对流项的显式处理，这种处理引入了随粘性减小而爆炸的误差界。这项工作强调了如果需要在高雷诺数下实现鲁棒性，则必须采用对流项的隐式处理或替代方案。作者并未提出修复此问题的方案，而是划定了现有一致性分裂方法的局限性。