[论文解读] A Survey of Trefftz Methods for the Helmholtz Equation

Ralf Hiptmair, Andrea Moiola+1 more|CentAUR (University of Reading)|Jun 15, 2015

Advanced Numerical Methods in Computational Mathematics参考文献 84被引用 8

一句话总结

本综述全面概述了用于Helmholtz方程的Trefftz方法，重点介绍最小二乘法、不连续伽辽金法和超弱变分格式等变分格式。它强调使用如广义调和多项式和平面波等振荡基函数，证明这些基函数在高频波问题中相比多项式具有更优的逼近性能。

ABSTRACT

Trefftz methods are finite element-type schemes whose test and trial functions are (locally) solutions of the targeted differential equation. They are particularly popular for time-harmonic wave problems, as their trial spaces contain oscillating basis functions and may achieve better approximation properties than classical piecewise-polynomial spaces. We review the construction and properties of several Trefftz variational formulations developed for the Helmholtz equation, including least squares, discontinuous Galerkin, ultra weak variational formulation, variational theory of complex rays and wave based methods. The most common discrete Trefftz spaces used for this equation employ generalised harmonic polynomials (circular and spherical waves), plane and evanescent waves, fundamental solutions and multipoles as basis functions; we describe theoretical and computational aspects of these spaces, focusing in particular on their approximation properties. One of the most promising, but not yet well developed, features of Trefftz methods is the use of adaptivity in the choice of the propagation directions for the basis functions. The main difficulties encountered in the implementation are the assembly and the ill-conditioning of linear systems, we briefly survey some strategies that have been proposed to cope with these problems.

研究动机与目标

• 系统性地回顾并分类专用于Helmholtz方程的Trefftz方法，强调其理论基础和数值格式。
• 分析基于广义调和多项式、平面波、基本解和偶极子的离散Trefftz空间的逼近性质。
• 解决实现中的关键挑战，包括病态性和系统组装问题，并评估迭代求解器的预条件策略。
• 探索自适应基函数选择（特别是方向自适应性）在提高精度和效率方面的潜力。
• 评估当前理论和数值验证的状态，识别Trefftz方法发展中存在的空白及未来研究方向。

提出的方法

• 将Trefftz变分格式分为三类主要类型：最小二乘法（在L2范数下最小化残差）、不连续伽辽金法（使用局部分部积分，可选拉格朗日乘子）和加权残差法（用测试函数的迹来检验残差）。
• 回顾使用广义调和多项式（圆形和球形波）、平面波和倏逝波、基本解以及偶极子等基函数的离散Trefftz空间，强调其振荡特性及其在波问题中的适用性。
• 分析Trefftz方法中产生的线性系统的条件性，特别是超弱变分格式（UWVF）中的条件性，并提出缩放技术（如使用贝塞尔函数范数进行归一化）以改善条件数。
• 研究预条件策略，包括基于平面波基函数L2正交化的局部预条件器，以及两层Schwarz区域分解方法，特别关注k-鲁棒性和可扩展性。
• 调查p-型超收敛和hp-型超收敛策略，特别是在奇点区域，通过使用角点函数或混合多项式-Trefftz单元来提升收敛性。
• 讨论使用BiCGStab等迭代求解器配合预条件器的应用，并将基于Trefftz的预条件器与非重叠区域分解方法相联系。

实验结果

研究问题

• RQ1 在Helmholtz方程中，不同Trefftz变分格式（最小二乘法、不连续伽辽金法、超弱变分格式）在稳定性、条件性和收敛性方面如何比较？
• RQ2 基于广义调和多项式、平面波和基本解构建的Trefftz空间具有怎样的逼近性质？与标准多项式有限元相比有何差异？
• RQ3 在高频Helmholtz问题中，基函数选择中的方向自适应性在多大程度上能提高精度并减少自由度？
• RQ4 预条件技术（特别是局部预条件器和两层Schwarz方法）如何影响Trefftz格式中迭代求解器的收敛性？它们对波数k是否具有鲁棒性？
• RQ5 在存在奇点或复杂几何形状的情况下，Trefftz方法的理论和实际限制是什么？如何通过hp-型超收敛或专用基函数来缓解这些问题？

主要发现

• 由于使用了与解的波动特性相匹配的振荡基函数，Trefftz方法在高频Helmholtz问题中相比标准多项式FEM具有更优的逼近性能。
• 通过使用贝塞尔函数范数对广义调和多项式进行缩放，可改善UWVF矩阵的条件性，且该缩放方式在规则网格上比基于平面波的方案具有更优的条件性。
• 数值实验表明，使用基本解作为基函数可在一定程度上降低条件性，但会牺牲精度，且当源中心远离单元时性能会提升。
• 基于平面波基函数L2正交化的局部预条件器效果显著，并与非重叠优化Schwarz区域分解方法相联系，增强了求解器的鲁棒性。
• 用于PWDG（等价于UWVF）的两层Schwarz预条件器对k和平面波方向数具有鲁棒性，但随着子区域数量增加，可扩展性仍受限。
• 在大单元上使用Trefftz函数进行p-型超收敛，特别是在结合角点函数或hp-型超收敛时，相较于多项式FEM在理论和实践上均具有显著优势，尤其在奇点区域。