Balancing domain decomposition for mixed finite elements

Journal Article · Sat Jul 01 04:00:00 EDT 1995 · Mathematics of Computation

DOI:https://doi.org/10.1090/S0025-5718-1995-1297465-9· OSTI ID:135844

Cowsar, L C; Mandel, J; Wheeler, M F

The rate of convergence of the Balancing Domain Decomposition method applied to the mixed finite element discretization of second-order elliptic equations is analyzed. The Balancing Domain Decomposition method, introduced recently by Mandel, is a substructuring method that involves at each iteration the solution of a local problem with Dirichlet data, a local problem with Neumann data, and a {open_quotes}coarse grid{close_quotes} problem to propagate information globally and to insure the consistency of the Neumann problems. It is shown that the condition number grows at worst like the logarithm squared of the ratio of the subdomain size to the element size, in both two and three dimensions and for elements of arbitrary order. The bounds are uniform with respect to coefficient jumps of arbitrary order. The bounds are uniform with respect to coefficient jumps of arbitrary size between subdomains. The key component of this analysis is the demonstration of an equivalence between the norm induced by the bilinear form on the interface and the H{sup 1/2}-norm of an interpolant of the boundary data. Computational results from a message-passing parallel implementation on an INTEL-Delta machine demonstrate the scalability properties of the method and show almost optimal linear observed speed-up for up to 64 processors.

• Sponsoring Organization: USDOE
• DOE Contract Number: FG05-92ER25142
• OSTI ID: 135844
• Journal Information: Mathematics of Computation, Journal Name: Mathematics of Computation Journal Issue: 211 Vol. 64; ISSN 0025-5718; ISSN MCMPAF
• Country of Publication: United States
• Language: English

Related Subjects

99 GENERAL AND MISCELLANEOUS
CALCULATION METHODS
DIFFERENTIAL EQUATIONS
FINITE ELEMENT METHOD
MATHEMATICS