Results 1  10
of
92
Schwarz Methods of NeumannNeumann Type for ThreeDimensional Elliptic Finite Element Problems
 Comm. Pure Appl. Math
, 1995
"... . Several domain decomposition methods of NeumannNeumann type are considered for solving the large linear systems of algebraic equations that arise from discretizations of elliptic problems by finite elements. We will only consider problems in three dimensions. Several new variants of the basic alg ..."
Abstract

Cited by 90 (21 self)
 Add to MetaCart
. Several domain decomposition methods of NeumannNeumann type are considered for solving the large linear systems of algebraic equations that arise from discretizations of elliptic problems by finite elements. We will only consider problems in three dimensions. Several new variants of the basic algorithm are introduced in a Schwarz method framework that provides tools which have already proven very useful in the design and analysis of other domain decomposition and multilevel methods. The NeumannNeumann algorithms have several advantages over other domain decomposition methods. The subregions, which define the subproblems, only share the boundary degrees of freedom with their neighbors. The subregions can also be of quite arbitrary shape and many of the major components of the preconditioner can be constructed from subprograms available in standard finite element program libraries. However, in its original form, the algorithm lacks a mechanism for global transportation of informatio...
FETI and NeumannNeumann iterative substructuring methods: connections and new results
 Comm. Pure Appl. Math
"... ..."
(Show Context)
Dualprimal FETI methods for threedimensional elliptic problems with heterogeneous coefficients
 SIAM J. Numer. Anal
, 2002
"... Abstract. In this paper, certain iterative substructuring methods with Lagrange multipliers are considered for elliptic problems in three dimensions. The algorithms belong to the family of dual{ primal FETI methods which have recently been introduced and analyzed successfully for elliptic problems i ..."
Abstract

Cited by 66 (12 self)
 Add to MetaCart
(Show Context)
Abstract. In this paper, certain iterative substructuring methods with Lagrange multipliers are considered for elliptic problems in three dimensions. The algorithms belong to the family of dual{ primal FETI methods which have recently been introduced and analyzed successfully for elliptic problems in the plane. The family of algorithms for three dimensions is extended and a full analysis is provided for the new algorithms. Particular attention is paid to nding algorithms with a small primal subspace since that subspace represents the only global part of the dual{primal preconditioner. It is shown that the condition numbers of several of the dual{primal FETI methods can be bounded polylogarithmically as a function of the dimension of the individual subregion problems and that the bounds are otherwise independent of the number of subdomains, the mesh size, and jumps in the coeÆcients. These results closely parallel those for other successful iterative substructuring methods of primal as well as dual type.
Overlapping Schwarz Methods On Unstructured Meshes Using NonMatching Coarse Grids
 Numer. Math
, 1996
"... . We consider two level overlapping Schwarz domain decomposition methods for solving the finite element problems that arise from discretizations of elliptic problems on general unstructured meshes in two and three dimensions. Standard finite element interpolation from the coarse to the fine grid may ..."
Abstract

Cited by 52 (17 self)
 Add to MetaCart
(Show Context)
. We consider two level overlapping Schwarz domain decomposition methods for solving the finite element problems that arise from discretizations of elliptic problems on general unstructured meshes in two and three dimensions. Standard finite element interpolation from the coarse to the fine grid may be used. Our theory requires no assumption on the substructures which constitute the whole domain, so each substructure can be of arbitrary shape and of different size. The global coarse mesh is allowed to be nonnested to the fine grid on which the discrete problem is to be solved and both the coarse meshes and the fine meshes need not be quasiuniform. In addition, the domains defined by the fine and coarse grid need not be identical. The one important constraint is that the closure of the coarse grid must cover any portion of the fine grid boundary for which Neumann boundary conditions are given. In this general setting, our algorithms have the same optimal convergence rate of the usual ...
Rate Of Convergence Of Some Space Decomposition Methods For Linear And Nonlinear Problems
 SIAM J. Numer. Anal
, 1998
"... .<F3.784e+05> Convergence of a space decomposition method is proved for a class of convex programming problems. A space decomposition refers to a method that decomposes a space into a sum of subspaces, which could be a domain decomposition or a multilevel method when applied to partial di#eren ..."
Abstract

Cited by 40 (15 self)
 Add to MetaCart
(Show Context)
.<F3.784e+05> Convergence of a space decomposition method is proved for a class of convex programming problems. A space decomposition refers to a method that decomposes a space into a sum of subspaces, which could be a domain decomposition or a multilevel method when applied to partial di#erential equations. Two algorithms are proposed. Both can be used for linear as well as nonlinear elliptic problems, and they reduce to the standard additive and multiplicative Schwarz methods for linear elliptic problems.<F4.005e+05> Key words.<F3.784e+05> parallel, domain decomposition, nonlinear, elliptic equation, space decomposition<F4.005e+05> AMS subject classifications.<F3.784e+05> 65J10, 65M55, 65Y05<F4.005e+05> PII.<F3.784e+05> S0036142996297461<F4.795e+05> 1. Introduction.<F4.397e+05> We use space decomposition methods to solve a convex programming problem. When the minimization space is suitably decomposed into subspaces, two algorithms are proposed to solve the minimization problem. The f...
Domain decomposition for multiscale PDEs
 Numer. Math
"... We consider additive Schwarz domain decomposition preconditioners for piecewise linear finite element approximations of elliptic PDEs with highly variable coefficients. In contrast to standard analyses, we do not assume that the coefficients can be resolved by a coarse mesh. This situation arises of ..."
Abstract

Cited by 39 (16 self)
 Add to MetaCart
(Show Context)
We consider additive Schwarz domain decomposition preconditioners for piecewise linear finite element approximations of elliptic PDEs with highly variable coefficients. In contrast to standard analyses, we do not assume that the coefficients can be resolved by a coarse mesh. This situation arises often in practice, for example in the computation of flows in heterogeneous porous media, in both the deterministic and (MonteCarlo simulated) stochastic cases. We consider preconditioners which combine local solves on general overlapping subdomains together with a global solve on a general coarse space of functions on a coarse grid. We perform a new analysis of the preconditioned matrix, which shows rather explicitly how its condition number depends on the variable coefficient in the PDE as well as on the coarse mesh and overlap parameters. The classical estimates for this preconditioner with linear coarsening guarantee good conditioning only when the coefficient varies mildly inside the coarse grid elements. By contrast, our new results show that, with a good choice of subdomains and coarse space basis functions, the preconditioner can still be robust even for large coefficient variation inside domains,
Nonlinearly preconditioned inexact Newton algorithms
 SIAM J. Sci. Comput
, 2000
"... Abstract. Inexact Newton algorithms are commonlyused for solving large sparse nonlinear system of equations F (u ∗ ) = 0 arising, for example, from the discretization of partial differential equations. Even with global strategies such as linesearch or trust region, the methods often stagnate at loc ..."
Abstract

Cited by 35 (15 self)
 Add to MetaCart
(Show Context)
Abstract. Inexact Newton algorithms are commonlyused for solving large sparse nonlinear system of equations F (u ∗ ) = 0 arising, for example, from the discretization of partial differential equations. Even with global strategies such as linesearch or trust region, the methods often stagnate at local minima of �F �, especiallyfor problems with unbalanced nonlinearities, because the methods do not have builtin machineryto deal with the unbalanced nonlinearities. To find the same solution u ∗ , one maywant to solve instead an equivalent nonlinearlypreconditioned system F(u ∗ ) = 0 whose nonlinearities are more balanced. In this paper, we propose and studya nonlinear additive Schwarzbased parallel nonlinear preconditioner and show numericallythat the new method converges well even for some difficult problems, such as high Reynolds number flows, where a traditional inexact Newton method fails. Key words. nonlinear preconditioning, inexact Newton methods, Krylov subspace methods, nonlinear additive Schwarz, domain decomposition, nonlinear equations, parallel computing, incompressible
A Polylogarithmic Bound For An Iterative Substructuring Method For Spectral Elements In Three Dimensions
 SIAM J. NUMER. ANAL
, 1993
"... Iterative substructuring methods form an important family of domain decomposition algorithms for elliptic finite element problems. A pversion finite element method based on continuous, piecewise Q p functions is considered for second order elliptic problems in three dimensions; this special method ..."
Abstract

Cited by 32 (4 self)
 Add to MetaCart
Iterative substructuring methods form an important family of domain decomposition algorithms for elliptic finite element problems. A pversion finite element method based on continuous, piecewise Q p functions is considered for second order elliptic problems in three dimensions; this special method can also be viewed as a conforming spectral element method. An iterative method is designed for which the condition number of the relevant operator grows only in proportion to (1 + log p)². This bound is independent of jumps in the coefficient of the elliptic problem across the interfaces between the subregions. Numerical results are also reported which support the theory.
Overlapping nonmatching grids mortar element methods for elliptic problems
 SIAM J. Numer. Anal
, 1999
"... Abstract. In the first part of the paper, we introduce an overlapping mortar finite element methods for solving twodimensional elliptic problems discretized on overlapping nonmatching grids. We prove an optimal error bound and estimate the condition numbers of certain overlapping Schwarz precondit ..."
Abstract

Cited by 30 (7 self)
 Add to MetaCart
(Show Context)
Abstract. In the first part of the paper, we introduce an overlapping mortar finite element methods for solving twodimensional elliptic problems discretized on overlapping nonmatching grids. We prove an optimal error bound and estimate the condition numbers of certain overlapping Schwarz preconditioned systems for the twosubdomain case. We show that the error bound is independent of the size of the overlap and the ratio of the mesh parameters. In the second part, we introduce three additive Schwarz preconditioned conjugate gradient algorithms based on the trivial and harmonic extensions. We provide estimates for the spectral bounds on the condition numbers of the preconditioned operators. We show that although the error bound is independent of the size of the overlap, the condition number does depend on it. Numerical examples are presented to support our theory.
Overlapping Schwarz Methods For Maxwell's Equations In Three Dimensions
 Numer. Math
, 1997
"... . Twolevel overlapping Schwarz methods are considered for finite element problems of 3D Maxwell's equations. N'ed'elec elements built on tetrahedra and hexahedra are considered. Once the relative overlap is fixed, the condition number of the additive Schwarz method is bounded, indepe ..."
Abstract

Cited by 28 (4 self)
 Add to MetaCart
(Show Context)
. Twolevel overlapping Schwarz methods are considered for finite element problems of 3D Maxwell's equations. N'ed'elec elements built on tetrahedra and hexahedra are considered. Once the relative overlap is fixed, the condition number of the additive Schwarz method is bounded, independently of the diameter of the triangulation and the number of subregions. A similar result is obtained for a multiplicative method. These bounds are obtained for quasiuniform triangulations. In addition, for the Dirichlet problem, the convexity of the domain has to be assumed. Our work generalizes wellknown results for conforming finite elements for second order elliptic scalar equations. 1. Introduction. When timedependent Maxwell's equations are considered, the electric field u satisfies the following equation curlcurlu + " @ 2 u @t 2 + oe @u @t = \Gamma @J @t ; in \Omega\Gamma (1) where J(x; t) is the current density and ", , oe describe the electromagnetic properties of the medium. For their...