[论文解读] Acoustic scattering: high frequency boundary element methods and unified transform methods

Simon N. Chandler‐Wilde, S. Langdon|arXiv (Cornell University)|Oct 22, 2014

Numerical methods in engineering参考文献 67被引用 6

一句话总结

本文提出了一种用于高频声散射的混合数值-渐近边界元方法，结合边界元公式与高频渐近分析以提高计算效率。同时，建立了该类方法与统一变换法之间的联系，通过严格的误差界和条件数分析，证明了在衍射光栅问题上SS∗变体的收敛性与稳定性。

ABSTRACT

We describe some recent advances in the numerical solution of acoustic scattering problems. A major focus of the paper is the efficient solution of high frequency scattering problems via hybrid numerical-asymptotic boundary element methods. We also make connections to the unified transform method due to A.S. Fokas and co-authors, analysing particular instances of this method, proposed by J.A. DeSanto and co-authors, for problems of acoustic scattering by diffraction gratings.

研究动机与目标

• 通过整合渐近分析与射线理论洞察，解决标准边界元方法（BEM）在高频下计算效率低下的问题。
• 开发并分析混合数值-渐近BEM公式，自适应地引入波状边界函数以减少自由度。
• 建立混合BEM与统一变换法（UTM）在周期性结构散射问题上的理论与数值联系。
• 为应用于软散射衍射光栅的统一变换法SS∗变体提供收敛性保证与条件数界。
• 通过在不同波数与几何形状下进行的数值实验，展示SS∗方法的鲁棒性与准确性。

提出的方法

• 将声散射问题表述为边界Γ上的边界积分方程，采用满足阻抗或狄利克雷边界条件的亥姆霍兹方程。
• 采用伽辽金法，以入射与出射平面波为基础构造试函数与试空间，具体使用ψ_n = γv(⋅,θ_n)作为试空间Q_N = P_N^Γ的基函数。
• 构造系统矩阵A_N，其元素为a_jm = ∫_Γ^L ψ_n_m ψ_n_j¯ ds，该矩阵在SS∗方法下为厄米特且正定。
• 采用SS∗方法，令χ_m = ψ_n_m = γv(⋅,θ_n_m)，通过Céa引理确保在L²(Γ^L)中实现最优逼近。
• 应用变分公式推导线性系统∑_m ∫_Γ^L χ_m ψ_n_j¯ ds c_m = -2ikL δ_0,n_j，以保证稳定性与收敛性。
• 分析系统矩阵A_N的条件数，表明只要基函数中仅包含实数方向的平面波，条件数便保持较低；但随着N增大，条件数会呈指数增长。

实验结果

研究问题

• RQ1 混合数值-渐近边界元方法是否能显著减少高频声散射问题中的自由度？
• RQ2 统一变换法中的SC、SS与SS∗方法在衍射光栅问题中的精度与条件数方面如何比较？
• RQ3 随着基函数数量N的增加，SS∗方法在L²(Γ^L)中的收敛行为如何？
• RQ4 系统矩阵A_N的条件数如何随N演化？这对数值稳定性有何影响？
• RQ5 统一变换法在多大程度上可与混合BEM框架严格关联并被其涵盖，用于散射问题？

主要发现

• SS∗方法在L²(Γ^L)中对法向导数∂_νu实现了最优逼近，其误差受试空间中最佳逼近的界所控制。
• SS∗方法的系统矩阵A_N为厄米特且正定，确保了有限N下的适定性与数值稳定性。
• 条件数κ_N = cond(A_N)在仅使用实数平面波方向的基函数时保持较低，但随N → ∞而呈指数增长。
• 数值实验表明，三种方法（SC、SS、SS∗）在正弦形轮廓下均能以每波长约一个自由度实现高精度。
• 当对所有充分大的N，有ψ_n ∈ P_N^Γ对所有n ∈ ℤ成立时，理论收敛性可得到保证。
• 由于其有利的矩阵结构与最优逼近性质，SS∗方法在可靠性上优于SC与SS方法。