[论文解读] On superconvergence of Runge-Kutta convolution quadrature for the wave equation

Jens Markus Melenk, Alexander Rieder|arXiv (Cornell University)|Apr 1, 2019

Electromagnetic Scattering and Analysis参考文献 45被引用 6

一句话总结

本文通过引入一种时间导数形式，解释了波动方程的Runge-Kutta卷积四次方（RK-CQ）中的超收敛现象，该形式提高了收敛速率。论文证明，使用入射波的时间导数作为输入，可通过一种新型的Dirichlet-to-Impedance映射估计，使频率依赖性降低一个|s|的幂次——多边形区域中带有对数因子——从而实现超越经典界限的改进收敛率。

ABSTRACT

The semidiscretization of a sound soft scattering problem modelled by the wave equation is analyzed. The spatial treatment is done by integral equation methods. Two temporal discretizations based on Runge-Kutta convolution quadrature are compared: one relies on the incoming wave as input data and the other one is based on its temporal derivative. The convergence rate of the latter is shown to be higher than previously established in the literature. Numerical results indicate sharpness of the analysis. The proof hinges on a novel estimate on the Dirichlet-to-Impedance map for certain Helmholtz problems. Namely, the frequency dependence can be lowered by one power of $\abs{s}$(up to a logarithmic term for polygonal domains) compared to the Dirichlet-to-Neumann map.

研究动机与目标

• 解释在RK-CQ处理波动方程时观察到的超收敛现象，即收敛速率超过经典预测。
• 分析两种时间离散化方法：一种使用入射波，另一种使用其时间导数。
• 通过推导Dirichlet-to-Impedance映射的更紧边界，为改进的收敛性建立理论基础。
• 将已知的球形或半空间几何下的边界扩展到一般光滑或多边形区域，且无需凸性假设。
• 通过在非凸（L形）区域上使用高阶Radau IIA格式的数值实验，验证理论发现。

提出的方法

• 使用边界积分方程公式化波动方程的散射问题，并对时间半离散化应用卷积四次方。
• 提出一种改进形式，其中输入数据为入射波的时间导数，导致卷积符号变为s⁻¹K(s)。
• 推导出Dirichlet-to-Impedance映射的新估计，显示频率依赖性降低一个|s|的幂次，多边形区域中带有对数因子。
• 利用改进的边界，证明改进方案的收敛速率相比标准形式提高一个阶次。
• 将理论应用于半离散设置，重点研究时间离散化，同时假设通过边界元法实现空间半离散化。
• 采用Galerkin边界元法进行空间离散化，并通过H⁻¹/²(Γ)范数计算误差，利用算子V(1)实现范数等价性。

实验结果

研究问题

• RQ1 为何在某些形式下，RK-CQ处理波动方程的收敛速率会超过经典预测？
• RQ2 能否通过边界积分算子的改进边界，严格解释超收敛现象？
• RQ3 改进的收敛性是否适用于非凸或多边形几何？几何形状如何影响收敛速率？
• RQ4 与标准形式相比，使用入射波的时间导数作为输入数据如何影响收敛阶？
• RQ5 理论估计是否精确？能否通过高阶时间积分格式的数值实验加以验证？

主要发现

• 使用入射波时间导数的RK-CQ格式的收敛速率相比标准形式提高一个阶次，符号边界中达到q+1−μ+1阶。
• 与Dirichlet-to-Neumann映射相比，Dirichlet-to-Impedance映射的频率依赖性降低一个|s|的幂次，多边形区域中带有对数因子。
• 数值实验验证了预测的收敛速率：三阶段Radau IIA方法达到三阶，五阶段方法达到五阶，且导数形式方法在对数项范围内达到完整的经典阶次。
• 即使在非凸的L形区域中也观察到改进的收敛性，表明理论分析在非凸几何中具有鲁棒性。
• 理论分析是精确的，因为数值结果与预测的收敛阶次一致，验证了Dirichlet-to-Impedance映射改进边界的正确性。
• 超收敛效应的解释在于：s⁻¹加权的Dirichlet-to-Neumann映射被分解为恒等算子与Dirichlet-to-Impedance映射之和，后者具有更优的频率边界。