Corners and stable optimized domain decomposition methods for the Helmholtz problem

We construct a new Absorbing Boundary Condition (ABC) adapted to solving the Helmholtz equation in polygonal domains in dimension two. Quasi-continuity relations are obtained at the corners of the polygonal boundary. This ABC is then used in the context of domain decomposition where various stable algorithms are constructed and analysed. Next, the operator of this ABC is adapted to obtain a transmission operator for the Domain Decomposition Method (DDM) that is well suited for broken line interfaces. For each algorithm, we show the decrease of an adapted quadratic pseudo-energy written on the skeleton of the mesh decomposition , which establishes the stability of these methods. Implementation within a finite element solver (GMSH/GetDP) and numerical tests illustrate the theory.

[1]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[2]  Yassine Boubendir,et al.  Non-overlapping Domain Decomposition Method for a Nodal Finite Element Method , 2006, Numerische Mathematik.

[3]  Sébastien Loisel,et al.  Condition Number Estimates for the Nonoverlapping Optimized Schwarz Method and the 2-Lagrange Multiplier Method for General Domains and Cross Points , 2013, SIAM J. Numer. Anal..

[4]  Martin J. Gander,et al.  Best Robin Parameters for Optimized Schwarz Methods at Cross Points , 2012, SIAM J. Sci. Comput..

[5]  Patrick Joly,et al.  Second-order absorbing boundary conditions for the wave equation: a solution for the corner problem , 1990 .

[6]  Olof B. Widlund,et al.  Extending Theory for Domain Decomposition Algorithms to Irregular Subdomains , 2008 .

[7]  Dirichlet to Neumann map for domains with corners and approximate boundary conditions , 2007 .

[8]  Victorita Dolean,et al.  An introduction to domain decomposition methods - algorithms, theory, and parallel implementation , 2015 .

[9]  Eli Turkel,et al.  Compact high order schemes with gradient-direction derivatives for absorbing boundary conditions , 2015, J. Comput. Phys..

[10]  Patrick Joly,et al.  Quasi-local transmission conditions for non-overlapping domain decomposition methods for the Helmholtz equation , 2014 .

[11]  Martin J. Gander,et al.  On the Applicability of Lions' Energy Estimates in the Analysis of Discrete Optimized Schwarz Methods with Cross Points , 2013, Domain Decomposition Methods in Science and Engineering XX.

[12]  Christophe Geuzaine,et al.  GetDDM: An open framework for testing optimized Schwarz methods for time-harmonic wave problems , 2016, Comput. Phys. Commun..

[13]  Christophe Geuzaine,et al.  A quasi-optimal non-overlapping domain decomposition algorithm for the Helmholtz equation , 2012, J. Comput. Phys..

[14]  Yassine Boubendir,et al.  Coupling of a non‐overlapping domain decomposition method for a nodal finite element method with a boundary element method , 2008 .

[15]  Axel Modave,et al.  Corner treatments for high-order local absorbing boundary conditions in high-frequency acoustic scattering , 2020, J. Comput. Phys..

[16]  A. Dhia,et al.  Perfect transmission invisibility for waveguides with sound hard walls , 2016, 1609.07596.

[17]  Xavier Antoine,et al.  Numerical accuracy of a Padé‐type non‐reflecting boundary condition for the finite element solution of acoustic scattering problems at high‐frequency , 2005 .

[18]  Frédéric Nataf,et al.  Improved Interface Conditions for 2D Domain Decomposition with Corners: Numerical Applications , 2009, J. Sci. Comput..

[19]  Bruno Després,et al.  A Domain Decomposition Method for the Helmholtz equation and related Optimal Control Problems , 1996 .

[20]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .