A Variable-Step IMEX-BDF2 Finite Element Method for the Navier–Stokes Equations: Optimal $L^2$ Error Estimates for Velocity and Pressure

技术摘要：一种用于 Navier–Stokes 方程的变步长 IMEX-BDF2 有限元方法

问题陈述
本文研究了非稳态不可压缩 Navier–Stokes (NS) 方程的数值模拟，这类方程具有内在的多尺度行为，特别是在高雷诺数下。虽然自适应时间步长是平衡精度与效率的公认工具，但针对二阶变步长方法（特别是二阶后向微分公式，即 BDF2）的现有理论框架仍存在显著缺陷。以往关于 NS 方程变步长 BDF2 方案的研究要么仅建立了阶数一的收敛性，要么在非最优范数（如 G-泛函范数）中提供误差估计，或者在 $L^2$ 空间中实现了最优时间阶数但在空间阶数上表现欠佳。此外，在变步长设定下，针对压力场的完整且最优的 $L^2$ 误差估计在文献中一直处于空白。

方法论
作者提出了一种全离散、变步长的隐式-显式 (IMEX) BDF2 有限元方法 (FEM)，采用 Taylor–Hood 元 ($P_2-P_1$)。该方案对对流项使用反对称形式以确保稳定性。关键的方法论组成部分包括：

1. 变时间步长： 该方法利用非均匀时间划分，其步长为 $\tau_k$，相邻步长比为 $r_k = \tau_k/\tau_{k-1}$。施加了一个稳定性条件，即 $0 < r_k \le r_{max} < 4.864$。
2. 初始化： 为了保持全局二阶精度，第一个时间步使用 Crank–Nicolson 格式计算，随后的步骤则采用变步长 BDF2 公式。
3. 离散正交卷积 (DOC) 核： 分析过程利用 DOC 核来处理非均匀时间步长。该框架允许精确量化初始步的影响，并管理误差分析中的时间步变化。
4. 误差分析策略：
   • $H^1$ 估计： 通过数学归纳法建立了速度的最优 $H^1$ 误差估计，并利用 DOC 框架推导出显式常数。
   • $L^2$ 估计（速度）： 为了在不牺牲最优性的前提下将 $H^1$ 估计升级为最优 $L^2$ 估计，作者采用了负范数论证方法。该技术避免了直接对非线性问题应用 Aubin–Nitsche 对偶论证的复杂性。
   • 压力估计： 利用负范数技术和离散 inf-sup 条件，同步导出了压力的最优 $L^2$ 误差估计。

主要贡献
本研究的主要贡献在于理论层面，填补了变步长 BDF2 方法分析中的特定空白：

• 最优 $L^2$ 误差估计： 本文建立了速度在 $L^2$ 范数下在空间和时间上的完全最优先验误差估计。这相对于以往受限于次优空间阶数或非 $L^2$ 范数的变步长分析是一项重大进展。
• 压力收敛性： 本工作提供了变步长 BDF2 设定下压力场首个最优的 $L^2$ 误差估计，这一结果此前在文献中是无法获得的。
• 初始化分析： 作者严谨地证明了使用 Crank–Nicolson 格式进行初始步计算可以保持变步长 BDF2 方法整体的二阶时间精度。
• 显式常数： 分析提供了网格尺寸和时间步约束的显式上界，确保了在特定正则性假设下的误差估计有效性。

结果
理论发现得到了使用 FreeFEM++ 进行的数值实验的支持：

• 精度验证： 在具有精确解的 2D 格点涡问题上的测试证实了预测的收敛速率。该方案在梯度网格和随机时间网格上，对速度和压力均实现了二阶时间收敛。
• 空间收敛性： 在固定时间步长下，该方法在速度的 $L^2$ 范数下表现出三阶收敛，在 $H^1$ 范数下表现出二阶收敛，且压力在 $L^2$ 范数下表现出二阶收敛。
• 鲁棒性： 即使当步长比超过理论界限 4.864 时（测试最高达到 7），该方案仍能保持预期的精度，这表明理论条件对于实际稳定性而言是充分而非必要的。
• 自适应应用： 在 $Re=1000$ 的顶盖驱动流模拟中，展示了自适应时间步长策略的有效性。该方法成功解析了快速流动转变过程中的微小步长，并在流动趋于稳态时增加步长，保持了平滑的步长演化。

意义
本文声称填补了变步长 BDF2 方法在 Navier–Stokes 方程分析中的显著理论空白。通过提供该方案完整的、最优的误差图谱（包括压力场），它巩固了这一广泛使用的变步长方法的理论基础。这项工作验证了变步长 IMEX-BDF2-FEM 可以通过严谨分析实现最优精度，使其成为模拟多尺度流动现象（其中自适应分辨率至关重要）的可靠工具。作者指出，这些结果为未来扩展到长时间稳定性估计及耦合多物理场系统奠定了基础。