Results 1  10
of
28
Balancing Domain Decomposition
 Comm. Numer. Meth. Engrg
, 1993
"... The NeumannNeumann algorithm is known to be an efficient domain decomposition preconditioner with unstructured subdomains for iterative solution of finite element discretizations of difficult problems with strongly discontinuous coefficients [6]. However, this algorithm suffers from the need to sol ..."
Abstract

Cited by 131 (11 self)
 Add to MetaCart
The NeumannNeumann algorithm is known to be an efficient domain decomposition preconditioner with unstructured subdomains for iterative solution of finite element discretizations of difficult problems with strongly discontinuous coefficients [6]. However, this algorithm suffers from the need to solve in each iteration an inconsistent singular problem for every subdomain, and its convergence deteriorates with increasing number of subdomains due to the lack of a coarse problem to propagate the error globally. We show that the equilibrium conditions for the singular problems on subdomains lead to a simple and natural construction of a coarse problem. The construction is purely algebraic and applies also to systems, such as those that arize in elasticity. A convergence bound independent on the number of subdomains is proved and results of computational tests are reported.
Schwarz Analysis Of Iterative Substructuring Algorithms For Elliptic Problems In Three Dimensions
 SIAM J. Numer. Anal
, 1993
"... . Domain decomposition methods provide powerful preconditioners for the iterative solution of the large systems of algebraic equations that arise in finite element or finite difference approximations of partial differential equations. The preconditioners are constructed from exact or approximate sol ..."
Abstract

Cited by 110 (26 self)
 Add to MetaCart
. Domain decomposition methods provide powerful preconditioners for the iterative solution of the large systems of algebraic equations that arise in finite element or finite difference approximations of partial differential equations. The preconditioners are constructed from exact or approximate solvers for the same partial differential equation restricted to a set of subregions into which the given region has been divided. In addition, the preconditioner is often augmented by a coarse, second level approximation, that provides additional, global exchange of information, and which can enhance the rate of convergence considerably. The iterative substructuring methods, based on decompositions of the region into nonoverlapping subregions, form one of the main families of such algorithms. Many domain decomposition algorithms can conveniently be described and analyzed as Schwarz methods. These algorithms are fully defined in terms of a set of subspaces and auxiliary bilinear forms. A gener...
Balancing Domain Decomposition For Mixed Finite Elements
 Math. Comp
"... . 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 by Mandel in [24], is a substructuring method that involves at each iteratio ..."
Abstract

Cited by 71 (20 self)
 Add to MetaCart
. 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 by Mandel in [24], 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 "coarse grid" 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 size between subdomains. The key component of our analysis is the demonstration of an equivalence between the norm induced by the bilinear form on the interface and the H 1=2 norm of an interpolant of the boundary ...
Balancing domain decomposition for problems with large jumps in coefficients
 Math. Comp
, 1996
"... Abstract. The Balancing Domain Decomposition algorithm uses in each iteration solution of local problems on the subdomains coupled with a coarse problem that is used to propagate the error globally and to guarantee that the possibly singular local problems are consistent. The abstract theory introdu ..."
Abstract

Cited by 57 (10 self)
 Add to MetaCart
Abstract. The Balancing Domain Decomposition algorithm uses in each iteration solution of local problems on the subdomains coupled with a coarse problem that is used to propagate the error globally and to guarantee that the possibly singular local problems are consistent. The abstract theory introduced recently by the firstnamed author is used to develop condition number bounds for conforming linear elements in two and three dimensions. The bounds are independent of arbitrary coefficient jumps between subdomains and of the number of subdomains, and grow only as the squared logarithm of the mesh size h. Computational experiments for two and threedimensional problems confirm the theory. 1.
Optimized Schwarz methods
 SIAM Journal on Numerical Analysis
, 2006
"... Abstract. Optimized Schwarz methods are a new class of Schwarz methods with greatly enhanced convergence properties. They converge uniformly faster than classical Schwarz methods and their convergence rates dare asymptotically much better than the convergence rates of classical Schwarz methods if th ..."
Abstract

Cited by 47 (11 self)
 Add to MetaCart
Abstract. Optimized Schwarz methods are a new class of Schwarz methods with greatly enhanced convergence properties. They converge uniformly faster than classical Schwarz methods and their convergence rates dare asymptotically much better than the convergence rates of classical Schwarz methods if the overlap is of the order of the mesh parameter, which is often the case in practical applications. They achieve this performance by using new transmission conditions between subdomains which greatly enhance the information exchange between subdomains and are motivated by the physics of the underlying problem. We analyze in this paper these new methods for symmetric positive definite problems and show their relation to other modern domain decomposition methods like the new Finite Element Tearing and Interconnect (FETI) variants.
A nonoverlapping domain decomposition method for Maxwellâ€™s equations in three dimensions
 SIAM J. Numer. Anal
"... Abstract. We propose a substructuring preconditioner for solving threedimensional elliptic equations with strongly discontinuous coefficients. The new preconditioner can be viewed as a variant of the classical substructuring preconditioner proposed by Bramble, Pasiack and Schatz (1989), but with muc ..."
Abstract

Cited by 35 (10 self)
 Add to MetaCart
Abstract. We propose a substructuring preconditioner for solving threedimensional elliptic equations with strongly discontinuous coefficients. The new preconditioner can be viewed as a variant of the classical substructuring preconditioner proposed by Bramble, Pasiack and Schatz (1989), but with much simpler coarse solvers. Though the condition number of the preconditioned system may not have a good bound, we are able to show that the convergence rate of the PCG method with such substructuring preconditioner is nearly optimal, and also robust with respect to the (possibly large) jumps of the coefficient in the elliptic equation. 1.
Some Nonoverlapping Domain Decomposition Methods
, 1998
"... . The purpose of this paper is to give a unified investigation of a class of nonoverlapping domain decomposition methods for solving secondorder elliptic problems in two and three dimensions. The methods under scrutiny fall into two major categories: the substructuringtype methods and the Neumann ..."
Abstract

Cited by 35 (6 self)
 Add to MetaCart
. The purpose of this paper is to give a unified investigation of a class of nonoverlapping domain decomposition methods for solving secondorder elliptic problems in two and three dimensions. The methods under scrutiny fall into two major categories: the substructuringtype methods and the NeumannNeumanntype methods. The basic framework used for analysis is the parallel subspace correction method or additive Schwarz method, and other technical tools include localglobal and globallocal techniques. The analyses for both two and threedimensional cases are carried out simultaneously. Some internal relationships between various algorithms are observed and several new variants of the algorithms are also derived. Key words. nonoverlapping domain decomposition, Schur complement, localglobal and globallocal techniques, jumps in coe#cients, substructuring, NeumannNeumann, balancing methods AMS subject classifications. 65N30, 65N55, 65F10 PII. S0036144596306800 1. Introduction. T...
A NeumannNeumann Domain Decomposition Algorithm for Solving Plate and Shell Problems
 SIAM J. NUMER. ANAL
, 1997
"... We present a new NeumannNeumann type preconditioner of large scale linear systems arising from plate and shell problems. The advantage of the new method is a smaller coarse space than those of earlier method of the authors; this improves parallel scalability. A new abstract framework for NeumannNe ..."
Abstract

Cited by 33 (8 self)
 Add to MetaCart
We present a new NeumannNeumann type preconditioner of large scale linear systems arising from plate and shell problems. The advantage of the new method is a smaller coarse space than those of earlier method of the authors; this improves parallel scalability. A new abstract framework for NeumannNeumann preconditioners is used to prove almost optimal convergence properties of the method. The convergence estimates are independent of the number of subdomains, coefficient jumps between subdomains, and depend only polylogarithmically on the number of elements per subdomain. We formulate and prove an approximate parametric variational principle for ReissnerMindlin elements as the plate thickness approaches zero, which makes the results applicable to a large class of nonlocking elements in everyday engineering use. The theoretical results are confirmed by computational experiments on model problems as well as examples from real world engineering practice.
Analysis Of Lagrange Multiplier Based Domain Decomposition
, 1998
"... The convergence of a substructuring iterative method with Lagrange multipliers known as Finite Element Tearing and Interconnecting (FETI) method is analyzed in this thesis. This method, originally proposed by Farhat and Roux, decomposes finite element discretization of an elliptic boundary value pro ..."
Abstract

Cited by 17 (5 self)
 Add to MetaCart
The convergence of a substructuring iterative method with Lagrange multipliers known as Finite Element Tearing and Interconnecting (FETI) method is analyzed in this thesis. This method, originally proposed by Farhat and Roux, decomposes finite element discretization of an elliptic boundary value problem into Neumann problems on the subdomains, plus a coarse problem for the subdomain null space components. For linear conforming elements and preconditioning by Dirichlet problems on the subdomains, the asymptotic bound on the condition number C(1 log(H=h)) fl , where fl = 2 or 3, is proved for a second order problem, h denoting the characteristic element size and H the size of subdomains. A similar method proposed by Park is shown to be equivalent to FETI with a special choice of some components and the bound C(1 log(H=h)) 2 on the condition number is established. Next, the original FETI method is generalized to fourth order plate bending problems. The main idea there is to enfor...