Metric-based anisotropic mesh adaptation for 3D acoustic boundary element methods

• Published in: Journal of computational physics Vol. 372; pp. 473 - 499
• Main Authors: Chaillat, Stéphanie, Groth, Samuel P., Loseille, Adrien
• Format: Journal Article
• Language: English
• Published: Cambridge Elsevier Inc 01-11-2018; Elsevier Science Ltd; Elsevier
• Subjects: Acoustic propagation; Acoustic wave scattering; Adaptation; Anisotropic meshes; Anisotropy; Boundary element method; Boundary integral method; Complexity; Computational physics; Convergence; Domains; Errors; Fast BEMs; Galerkin method; Indicators; Integral equations; Mathematical analysis; Mechanics; Mesh adaptation; Physics; Singularities; Solid mechanics; Surface acoustic waves; Wave propagation; Mesh adaptation; Fast BEMs; Acoustic wave scattering; Anisotropic meshes; Boundary element method; Boundary Element Method; fast BEMs
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2018.06.048

•The extension of a metric-based mesh adaptation strategy for the BEM is presented.•The method generates truly anisotropic meshes.•The method is independent of the PDE and discretization technique.•Numerical results confirm the recovery of optimal convergence rates.
This paper details the extension of a metric-based anisotropic mesh adaptation strategy to the boundary element method for problems of 3D acoustic wave propagation. Traditional mesh adaptation strategies for boundary element methods rely on Galerkin discretizations of the boundary integral equations, and the development of appropriate error indicators. They often require the solution of further integral equations. These methods utilize the error indicators to mark elements where the error is above a specified tolerance and then refine these elements. Such an approach cannot lead to anisotropic adaptation regardless of how these elements are refined, since the orientation and shape of current elements cannot be modified.
In contrast, the method proposed here is independent of the discretization technique (e.g., collocation, Galerkin). Furthermore, it completely remeshes at each refinement step, altering the shape, size, and orientation of each element according to an optimal metric based on a numerically recovered Hessian of the boundary solution. The resulting adaptation procedure is truly anisotropic and independent of the complexity of the geometry. We show via a variety of numerical examples that it recovers optimal convergence rates for domains with geometric singularities. In particular, a faster convergence rate is recovered for scattering problems with complex geometries.