[论文解读] Applying GMRES to the Helmholtz equation with strong trapping: how does the number of iterations depend on the frequency?

Pierre Marchand, Jeffrey Galkowski+2 more|arXiv (Cornell University)|Feb 10, 2021

Electromagnetic Scattering and Analysis参考文献 96被引用 7

一句话总结

本文分析了高频Helmholtz边界积分方程在强局域化条件下GMRES的频率依赖收敛性，证明了在谱假设下，GMRES迭代次数的上界随频率呈多项式增长。该研究首次对局域化Helmholtz问题中的迭代次数缩放行为进行了全面的理论与数值研究。

ABSTRACT

We consider GMRES applied to discretisations of the high-frequency Helmholtz equation with strong trapping; recall that in this situation the problem is exponentially ill-conditioned through an increasing sequence of frequencies. Under certain assumptions about the distribution of the eigenvalues, we prove upper bounds on how the number of GMRES iterations grows with the frequency. Our main focus is on boundary-integral-equation formulations of the exterior Dirichlet and Neumann obstacle problems in 2- and 3-d; for these problems, we investigate numerically the sharpness (in terms of dependence on frequency) of both our bounds and various quantities entering our bounds. This paper is therefore the first comprehensive study of the frequency-dependence of the number of GMRES iterations for Helmholtz boundary-integral equations under trapping.

研究动机与目标

• 研究高频Helmholtz问题中强局域化条件下GMRES迭代次数对频率的依赖关系。
• 解决由局域化几何中准模引起的指数型病态性问题。
• 为Dirichlet和Neumann散射问题的边界积分公式建立GMRES迭代次数的理论上限。
• 通过数值实验验证理论边界的紧致性及其对积分算子谱性质的依赖性。
• 将结果推广至具有超代数小质量准模的一般障碍物，而不仅限于特定马蹄形区域。

提出的方法

• 使用Dirichlet和Neumann条件的边界积分方程来表述Helmholtz外部散射问题。
• 应用伽辽金边界元法对积分方程进行离散化，得到大型稀疏线性系统。
• 使用GMRES求解所得线性系统，并分析迭代次数对波数 $k$ 的依赖性。
• 基于边界积分算子特征值分布的假设，推导出GMRES迭代次数的上界。
• 在二维马蹄形区域及其三维类比结构上进行数值实验，以检验理论边界的紧致性。
• 通过椭圆坐标系中的相空间体积计算，估计局域化轨迹的测度，并将其与谱性质关联。

实验结果

研究问题

• RQ1 在强局域化Helmholtz问题中，GMRES迭代次数随频率 $k$ 增加如何变化？
• RQ2 在对积分算子特征值分布做出假设的前提下，能否建立GMRES迭代次数的理论上限？
• RQ3 在实际中，这些理论边界有多紧致，特别是其对 $k$ 的依赖性如何？
• RQ4 进入边界的量（如谱条件数、准模质量）在多大程度上影响迭代次数？
• RQ5 该分析是否可推广至具有超代数小质量准模的一般障碍物，而不仅限于具有指数小质量准模的障碍物？

主要发现

• 在边界积分算子特征值分布假设下，GMRES迭代次数随频率 $k$ 的增长至多为多项式增长。
• 对于具有指数小质量准模的障碍物（如某些二维马蹄形区域），理论迭代次数上限的增长形式为 $\mathcal{O}(k^{C})$，其中 $C$ 为某常数。
• 数值实验表明，迭代次数随 $k$ 增加，且观测到的增长率与理论边界一致。
• 准模质量（以 $\epsilon(k) = C_1 \exp(-C_2 k)$ 表征）直接影响病态性，从而影响迭代次数。
• 边界在紧致性意义上是紧的，即频率依赖性与观测到的数值行为相匹配，验证了理论框架的有效性。
• 该分析适用于具有超代数小质量准模的一类广泛障碍物，不限于数值实验中研究的具体几何形状。