García-Archilla 等用 BDF-q（q≤5）时间积分 POD-ROM 求解半线性反应扩散方程，证明并数值验证了 q 阶最优时间收敛率

作者证明了对半线性反应扩散方程使用 BDF-q（1≤q≤5）时间积分的全离散 POD-ROM，在时间方向具有 q 阶最优收敛率，并给出点态时间误差界。 此前 POD 误差分析文献大多只处理隐式 Euler 或至多二阶时间格式，本文把高阶 BDF 时间离散的误差分析扩展到降阶模型。 以稳定性定理（Theorem 1）、相容性定理（Theorem 2）和收敛性定理（Theorem 3、4）给出先验误差界，误差界同时包含特征值尾部、时间误差与空间离散误差项；数值实验在 Brusselator 系统上测得斜率约为 -1.17、-1.96、-2.98、-3.98、-4.93，对应 q=1 至 5。

作者证明任意 BDF-q 格式可写成若干一阶差商的线性组合，并据此把差商纳入快照集合，从而控制时间导数在降阶空间上的投影误差。 这一代数恒等式此前仅对 BDF2 使用过，本文将其推广到 q≤5，使高阶时间格式的降阶误差分析成为可能。 在相容性定理证明中给出显式分解 ∂q u_h^n = (1/Δt)Σ α_j D u_h^{n-j}，并利用快照差商与特征值尾部的关系式（18）得到界（53）。

作者采用 H1_0 投影构造 POD 基，使误差界在特征值尾部意义下达到最优。 相比 L2 投影，H1_0 投影可给出以特征值尾部表示的最优界；文中也说明分析对 L2 投影同样适用，但此时不再需要差商线性组合这一关键步骤。 误差界（55）同时含特征值尾部 Σ_{k>r} λ_k、时间项 (Δt)^{2q} 与空间项 h^{k+1}，作者称其为在特征值尾部、时间误差和空间离散误差三方面均最优。

数值实验显示，随着降阶维数 r 增大，误差由时间积分器主导，需要更高阶 BDF 才能在不增加时间步数的情况下达到最优精度。 这为实践中选择时间格式阶数提供了具体依据：r=26 需二阶、r=32 需三阶、r=42 需四阶、r=50 需五阶。 在 Δt=T/1024 固定下比较 q=1 至 5 的最大 L2 误差（Table 3），并报告 Newton 迭代步数在 2 至 4 之间、平均 2.97。

该结果面向使用 POD-ROM 求解半线性反应扩散类抛物问题的数值分析研究者与科学计算实践者，适用场景是快照来自离线阶段足够精确的有限元近似、时间步满足条件（24）的情形。它使高阶 BDF 时间格式在降阶模型中有了可预期的 q 阶收敛保证，从而支持在保持精度前提下减少时间步数；分析框架也可为其他抛物型降阶模型的高阶时间离散提供参考。

分析假设非线性项全局 Lipschitz 且快照无时间误差，实际中快照的离线时间误差会叠加到误差界上，其影响程度需结合具体问题评估。BDF6 的分析被明确留作单独处理，q>5 的情形尚未覆盖。数值验证集中在 Brusselator 系统与特定参数（ν=0.002、二次有限元、51200 自由度），其他方程与参数下的表现仍是开放问题。此外，达到理论收敛率所需的时间步随 q 增大而减小，实际计算中如何权衡阶数与步长值得进一步观察。