Non-iterative domain decomposition for the Helmholtz equation using the method of difference potentials (original) (raw)
2021, ArXiv
We use the Method of Difference Potentials (MDP) to solve a non-overlapping domain decomposition formulation of the Helmholtz equation. The MDP reduces the Helmholtz equation on each subdomain to a Calderon’s boundary equation with projection on its boundary. The unknowns for the Calderon’s equation are the Dirichlet and Neumann data. Coupling between neighboring subdomains is rendered by applying their respective Calderon’s equations to the same data at the common interface. Solutions on individual subdomains are computed concurrently using a straightforward direct solver. We provide numerical examples demonstrating that our method is insensitive to interior cross-points and mixed boundary conditions, as well as large jumps in the wavenumber for transmission problems, which are known to be problematic for many other Domain Decomposition Methods.
Related papers
A Non-Overlapping Domain Decomposition Method for the Exterior Helmholtz Problem
In this paper, we first show that the domain decomposition methods that are usually efficient for solving elliptic problems typically fail when applied to acoustics problems. Next, we present an alternative domain decomposition algorithm that is better suited for the exterior Helmholtz problem. We describe it in a formalism that can use either one or two Lagrange multiplier fields for solving the corresponding interface problem by a Krylov method. In order to improve convergence and ensure scalability with respect the number of subdomains, we propose two complementary preconditioning techniques. The first preconditioner is based on a spectral analysis of the resulting interface operator and targets the high frequency components of the error. The second preconditioner is based on a coarsening technique, employs plane waves, and addresses the low frequency components of the error. Finally, we show numerically that, using both preconditioners, the convergence rate of the proposed domain decomposition method is quasi independent of the number of elements in the mesh, the number of subdomains, and depends only weakly on the wave number, which makes this method uniquely suitable for solving large scale high frequency exterior acoustics problems.
A new iterative non-overlapping domain decomposition method is proposed for solving the one-and two-dimensional Helmholtz equation on parallel computers. The spectral collocation method is applied to solve the Helmholtz equation in each subdomain based on the Chebyshev approximation, while the patching conditions are imposed at the interfaces between subdomains through a correction, being a linear function of the space coordinates. Convergence analysis is performed for two applications of the proposed method (DDLC and DDNNLC algorithms – the meaning of these abbreviations is explained below) based on the works of Zanolli and Funaro et al. Numerical tests have been performed and results obtained using the proposed method and other iterative algorithms have been compared. Parallel performance of the multi-domain algorithms has been analyzed by decomposing the two-dimensional domain into a number of subdomains in one spatial direction. For the one–dimensional problem, convergence of th...
Loading Preview
Sorry, preview is currently unavailable. You can download the paper by clicking the button above.