Parallel Computation of Harmonic Waves Using Domain Decomposition: Part 2 — Numerical Implementation and Applications

Abstract

Abstract A numerical implementation of the Domain Decomposition (DD) method for parallel computation of time-harmonic waves is presented. The Finite Element method on unstructured mesh is employed to discretize the boundary-value problem. A nonover-lapping mesh partitioning is used to build an overlapping domain decomposition. The algorithm also defines the new subdomain boundaries, other than the original boundaries of the global computational domain, on which transmission conditions are introduced. The independent subdomain problems can be solved concurrently, one problem per processor, on a parallel computer. The iterative scheme used to update the subdomain transmission conditions is computationally efficient and only swaps the nodal values in the overlapping region between adjacent subdomains. For exterior problems, the communication required to update the conditions on the outer boundary is substantially reduced by employing a modified Dirichlet-to-Neumann (DtN) map with truncated support. Numerical results are presented for a radiation model problem. It is shown that the DtN nonreflecting conditions eliminate spurious reflections and that the iterative scheme leads to a continuous global solution. The DD algorithm also has a very good convergence rate and significantly reduces memory requirement in comparison with a direct solver.