中文

时变边界条件下双曲型 PDE 显式 Runge-Kutta 积分的收敛阶数恢复

数值分析 2026-04-13 v1 数值分析 计算物理

摘要

对于具有时变 Dirichlet 边界数据的双曲初值-边值问题,显式 Runge-Kutta (RK) 积分常出现阶数下降:观察到的收敛阶数低于名义阶数,因为阶段结构与靠近边界的非对称空间闭合相互作用。本文发展一种纯空间的补偿方法,既保留时间积分器,又仅重新设计前两个边界相邻导数算子。对于任意 s 阶段 RK 方法作用于线性平流问题,边界相邻节点的一阶截断误差被证明可容纳由 tableau 决定的分解,其取消可显式给出边界权值的代数条件。解的可行性系数 R(b,c,A)R(\mathbf{b},\mathbf{c},A) 决定是否存在空间补偿机制;结果特别适用于 SSP-RK3,推导闭式条件。受限的差分演化随后识别出 5 点闭合方案,这些方案结合 5 阶中心差分内部模板,可从经典 Taylor 闭合退化为二阶行为中恢复到三阶收敛。一种考虑稳定性的变体在优化中引入特征值惩罚,揭示了阶数恢复与 CFL 鲁棒性之间的权衡。验证涵盖线性平流、manufactured-solution Burgers 流动以及二维分离平移的二维平流。分析阐明了为何时序弱阶数的修复无法解决有限差分边界问题,并指出该框架如何推广至非均匀网格。

关键词

引用

@article{arxiv.2604.09183,
  title  = {Restoring Convergence Order in Explicit Runge-Kutta Integration of Hyperbolic PDE with Time-Dependent Boundary Conditions},
  author = {Giorgio Maria Cavallazzi and Miguel Pérez Cuadrado and Alfredo Pinelli},
  journal= {arXiv preprint arXiv:2604.09183},
  year   = {2026}
}

备注

22 pages, 13 figures