[论文解读] Isogeometric finite element analysis of time-harmonic exterior acoustic scattering problems

Tahsin Khajah, Xavier Antoine+1 more|arXiv (Cornell University)|Oct 6, 2016

Advanced Numerical Analysis Techniques参考文献 50被引用 13

一句话总结

本论文提出使用基于NURBS的高阶基函数进行等几何有限元分析（IGA），以求解时 harmonic 外部声散射问题。通过在虚构边界上应用吸收边界条件，并将截断误差纳入精确解中，研究成功分离并评估了污染误差。主要发现是，当 p ≥ 3 时，IGA 在高波数（k = 500）下表现出可忽略的污染误差，证明了其在高频声学问题中的鲁棒性，原因在于其高阶精度和精确的几何表示。

ABSTRACT

We present an isogeometric analysis of time-harmonic exterior acoustic problems. The infinite space is truncated by a fictitious boundary and (simple) absorbing boundary conditions are applied. The truncation error is included in the exact solution so that the reported error is an indicator of the performance of the isogeometric analysis, in particular of the related pollution error. Numerical results performed with high-order basis functions (third or fourth orders) showed no visible pollution error even for very high frequencies. This property combined with exact geometrical representation makes isogeometric analysis a very promising platform to solve high-frequency acoustic problems.

研究动机与目标

• 评估基于NURBS的高阶基函数在时谐外部声散射问题中的等几何分析（IGA）性能。
• 通过将吸收边界条件引起的截断误差纳入精确解析解中，分离并量化IGA中的污染误差。
• 研究IGA在高波数（k = 40, 200, 500）下，不同离散化密度（nλ）和基函数阶次（p）时的收敛行为与精度。
• 评估IGA作为高频声学问题稳健平台的适用性，利用其精确的几何表示和高阶近似能力。

提出的方法

• 通过在虚构边界上施加简单的吸收边界条件（ABC），对无限外部区域进行截断，引入连续层面的截断误差。
• 将截断误差纳入解析精确解中，以在数值解中分离出污染误差。
• 使用基于NURBS的等几何有限元方法（IGAFEM）和高阶基函数（p = 1 至 5）对Helmholtz方程进行空间离散化。
• 在有界计算域上，使用Galerkin方法结合高阶NURBS形函数求解Helmholtz问题的弱形式。
• 通过计算精确解与数值解之间的相对 L2 误差 ε₂，量化污染误差和收敛行为。
• 采用平面波入射和模态分析方法，验证不同波数和离散化密度下的方法有效性。

实验结果

研究问题

• RQ1 基于高阶NURBS基函数的等几何分析在时谐外部声散射问题中是否表现出更低的污染误差？
• RQ2 IGA中的污染误差如何依赖于波数 k 和每波长的自由度数（nλ）？
• RQ3 在极高波数（k = 500）下，即使离散化密度较低（例如 nλ = 5），高阶IGA是否仍能保持精度？
• RQ4 IGA中基函数的光滑性与精确的几何表示在多大程度上有助于最小化高频声学问题中的污染误差？

主要发现

• 当 p ≥ 3 时，IGA中的污染误差在高波数下可忽略不计，且随着 k 增大，误差几乎保持不变（例如，k 从 40 增至 200 时）。
• 在 k = 40 且 nλ = 5 时，p = 3 的相对 L2 误差 ε₂ 约为 2–3%，但当 p 提高至 4 时，误差降至约 10⁻³。
• 在 k = 200 且 nλ = 5 时，绝对误差 |u_ex - u_h| 约为 10⁻³ 量级，表明即使在低离散化密度下也具有高精度。
• 在 k = 500 且 nλ = 5 时，p = 3 的误差约为 6%，但当 p 提高至 4 时，误差降低至约 10⁻³，证明了更高阶基函数的有效性。
• 对于 p ≥ 3，误差随 k 的增加并不显著上升，表明污染效应极小；而低阶方法（p = 1, 2）则表现出更强的 k 依赖性。
• IGA中精确几何表示与高阶近似的结合，使其在高频声散射问题中极具适用性。