Premium
A two‐level domain decomposition method with accurate interface conditions for the Helmholtz problem
Author(s) -
Astaneh Ali Vaziri,
Guddati Murthy N.
Publication year - 2015
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.5164
Subject(s) - domain decomposition methods , lagrange multiplier , helmholtz equation , helmholtz free energy , mathematics , dirichlet distribution , algorithm , interface (matter) , mathematical optimization , mathematical analysis , computer science , boundary value problem , finite element method , structural engineering , physics , bubble , quantum mechanics , maximum bubble pressure method , parallel computing , engineering
Summary A new and efficient two‐level, non‐overlapping domain decomposition (DD) method is developed for the Helmholtz equation in the two Lagrange multiplier framework. The transmission conditions are designed by utilizing perfectly matched discrete layers (PMDLs), which are a more accurate representation of the exterior Dirichlet‐to‐Neumann map than the polynomial approximations used in the optimized Schwarz method. Another important ingredient affecting the convergence of a DD method, namely, the coarse space augmentation, is also revisited. Specifically, the widely successful approach based on plane waves is modified to that based on interface waves, defined directly on the subdomain boundaries, hence ensuring linear independence and facilitating the estimation of the optimal size for the coarse problem. The effectiveness of both PMDL‐based transmission conditions and interface‐wave‐based coarse space augmentation is illustrated with an array of numerical experiments that include comprehensive scalability studies with respect to frequency, mesh size and the number of subdomains. Copyright © 2015 John Wiley & Sons, Ltd.