A dual-primal variant of the FETI-H domain decomposition method is designed for the fast, parallel, iterative solution of large-scale systems of complex equations arising from the discretization of acoustic scattering problems formulated in bounded computational domains. The convergence of this iterative solution method, named here FETI-DPH, is shown to scale with the problem size, the number of subdomains, and the wave number. Its solution time is also shown to scale with the problem size. CPU performance results obtained for the acoustic signature analysis in the mid-frequency regime of mockup submarines reveal that the proposed FETI-DPH solver is significantly faster than the previous generation FETI-H solution algorithm.
In the Dual-Primal FETI method, introduced by Farhat et al. [5], the domain is decomposed into non-overlapping subdomains, but the degrees of freedom on crosspoints remain common to all subdomains adjacent to the crosspoint. The continuity of the remaining degrees of freedom on subdomain interfaces is enforced by Lagrange multipliers and all degrees of freedom are eliminated. The resulting dual problem is solved by preconditioned conjugate gradients. We give an algebraic bound on the condition number, assuming only a single inequality in discrete norms, and use the algebraic bound to show that the condition number is bounded by C(1+log 2 (H/h)) for both second and fourth order elliptic selfadjoint problems discretized by conforming finite elements, as well as for a wide class of finite elements for the Reissner-Mindlin plate model.
Mathematics Subject Classification (1991): 65N55
We analyze the convergence of a substructuring iterative method with Lagrange multipliers, proposed recently by Farhat and Roux. The method decomposes finite element discretization of an elliptic boundary value problem into Neumann problems on the subdomains plus a coarse problem for the subdomain nullspace components. For linear conforming elements and preconditioning by the Dirichlet problems on the subdomains, we prove the asymptotic bound on the condition number C (1 + log(H /h)) γ , γ = 2 or 3,where h is the characteristic element size and H subdomain size.
SUMMARYRecently, a discontinuous Galerkin method with plane wave basis functions and Lagrange multiplier degrees of freedom was proposed for the efficient solution of Helmholtz problems in the mid-frequency regime. In this paper, this method is extended to higher-order elements. Performance results obtained for various two-dimensional problems highlight the advantages of these elements over classical higherorder Galerkin elements such as Q 2 and Q 4 for the discretization of interior and exterior Helmholtz problems.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.