Results 1  10
of
91
Analysis of the Inexact Uzawa Algorithm for Saddle Point Problems
 SIAM J. Numer. Anal
, 1997
"... . In this paper, we consider the socalled "inexact Uzawa" algorithm for iteratively solving block saddle point problems. Such saddle point problems arise, for example, in finite element and finite difference discretizations of Stokes equations, the equations of elasticity and mixed finite ..."
Abstract

Cited by 72 (2 self)
 Add to MetaCart
. In this paper, we consider the socalled "inexact Uzawa" algorithm for iteratively solving block saddle point problems. Such saddle point problems arise, for example, in finite element and finite difference discretizations of Stokes equations, the equations of elasticity and mixed finite element discretization of second order problems. We consider both the linear and nonlinear variants of the inexact Uzawa algorithm. We show that the linear method always converges as long as the preconditioners defining the algorithm are properly scaled. Bounds for the rate of convergence are provided in terms of the rate of convergence for the preconditioned Uzawa algorithm and the reduction factor corresponding to the preconditioner for the upper left hand block. In the nonlinear case, the inexact Uzawa algorithm is shown to converge provided that the nonlinear process approximating the inverse of the upper left hand block is of sufficient accuracy. Bounds for the nonlinear iteration are given in t...
FirstOrder System Least Squares For SecondOrder Partial Differential Equations: Part I
, 1994
"... . This paper develops ellipticity estimates and discretization error bounds for elliptic equations (with lower order terms) that are reformulated as a leastsquares problem for an equivalent firstorder system. The main result is the proof of ellipticity, which is used in a companion paper to esta ..."
Abstract

Cited by 61 (14 self)
 Add to MetaCart
. This paper develops ellipticity estimates and discretization error bounds for elliptic equations (with lower order terms) that are reformulated as a leastsquares problem for an equivalent firstorder system. The main result is the proof of ellipticity, which is used in a companion paper to establish optimal convergence of multiplicative and additive solvers of the discrete systems. Key words. leastsquares discretization, secondorder elliptic problems, RayleighRitz, finite elements 1. Introduction. The purpose of this paper is to analyse the leastsquares finite element method for secondorder convectiondiffusion equations written as a firstorder system. In general, the standard Galerkin finite element methods applied to nonselfadjoint elliptic equations with significant convection terms exhibit a variety of deficiencies, including oscillations or nonmonotonocity of the solution and poor approximation of its derivatives. A variety of stabilization techniques, such as upwin...
New Estimates for Multilevel Algorithms Including the VCycle
 MATH. COMP
, 1993
"... The purpose of this paper is to provide new estimates for certain multilevel algorithms. In particular, we are concerned with the simple additive multilevel algorithm given in [10] and the standard Vcycle algorithm with one smoothing step per grid. We shall prove that these algorithms have a unifo ..."
Abstract

Cited by 52 (5 self)
 Add to MetaCart
The purpose of this paper is to provide new estimates for certain multilevel algorithms. In particular, we are concerned with the simple additive multilevel algorithm given in [10] and the standard Vcycle algorithm with one smoothing step per grid. We shall prove that these algorithms have a uniform reduction per iteration independent of the mesh sizes and number of levels even on nonconvex domains which do not provide full elliptic regularity. For example, the theory applies to the standard multigrid Vcycle on the Lshaped domain or a domain with a crack and yields a uniform convergence rate. We also prove uniform convergence rates for the multigrid Vcycle for problems with nonuniformly refined meshes. Finally, we give a new multigrid approach for problems on domains with curved boundaries and prove a uniform rate of convergence for the corresponding multigrid Vcycle algorithms.
An Algorithm for Coarsening Unstructured Meshes
 Numer. Math
, 1996
"... . We develop and analyze a procedure for creating a hierarchical basis of continuous piecewise linear polynomials on an arbitrary, unstructured, nonuniform triangular mesh. Using these hierarchical basis functions, we are able to define and analyze corresponding iterative methods for solving the lin ..."
Abstract

Cited by 47 (5 self)
 Add to MetaCart
. We develop and analyze a procedure for creating a hierarchical basis of continuous piecewise linear polynomials on an arbitrary, unstructured, nonuniform triangular mesh. Using these hierarchical basis functions, we are able to define and analyze corresponding iterative methods for solving the linear systems arising from finite element discretizations of elliptic partial differential equations. We show that such iterative methods perform as well as those developed for the usual case of structured, locally refined meshes. In particular, we show that the generalized condition numbers for such iterative methods are of order J 2 , where J is the number of hierarchical basis levels. Key words. Finite element, hierarchical basis, multigrid, unstructured mesh. AMS subject classifications. 65F10, 65N20 1. Introduction. Iterative methods using the hierarchical basis decomposition have proved to be among the most robust for solving broad classes of elliptic partial differential equations, ...
Iterative Techniques for Time Dependent Stokes Problems
 Comput. Math. Appl
, 1994
"... . In this paper, we consider solving the coupled systems of discrete equations which arise from implicit time stepping procedures for the time dependent Stokes equations using a mixed finite element spatial discretization. At each time step, a two by two block system corresponding to a perturbed Sto ..."
Abstract

Cited by 40 (3 self)
 Add to MetaCart
. In this paper, we consider solving the coupled systems of discrete equations which arise from implicit time stepping procedures for the time dependent Stokes equations using a mixed finite element spatial discretization. At each time step, a two by two block system corresponding to a perturbed Stokes problem must be solved. Although there are a number of techniques for iteratively solving this type of block system, to be effective, they require a good preconditioner for the resulting pressure operator (Schur complement). In contrast to the time independent Stokes equations where the pressure operator is well conditioned, the pressure operator for the perturbed system becomes more ill conditioned as the time step is reduced (and/or the Reynolds number is increased). In this paper, we shall describe and analyze preconditioners for the resulting pressure systems. These preconditioners give rise to iterative rates of convergence which are independent of both the mesh size h as well as t...
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 36 (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 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 36 (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.
The Auxiliary Space Method And Optimal Multigrid Preconditioning Techniques For Unstructured Grids
 Computing
, 1996
"... . An abstract framework of auxiliary space method is proposed and, as an application, an optimal multigrid technique is developed for general unstructured grids. The auxiliary space method is a (nonnested) two level preconditioning technique based on a simple relaxation scheme (smoother) and an auxi ..."
Abstract

Cited by 33 (3 self)
 Add to MetaCart
. An abstract framework of auxiliary space method is proposed and, as an application, an optimal multigrid technique is developed for general unstructured grids. The auxiliary space method is a (nonnested) two level preconditioning technique based on a simple relaxation scheme (smoother) and an auxiliary space (that may be roughly understood as a nonnested coarser space). An optimal multigrid preconditioner is then obtained for a discretized partial differential operator defined on an unstructured grid by using an auxiliary space defined on a more structured grid in which a further nested multigrid method can be naturally applied. This new technique make it possible to apply multigrid methods to general unstructured grids without too much more programming effort than traditional solution methods. Some simple examples are also given to illustrate the abstract theory and for instance the Morley finite element space is used as an auxiliary space to construct a preconditioner for Argyris ...
Preconditioning in H(div) and Applications
 Math. Comp
, 1998
"... . We consider the solution of the system of linear algebraic equations which arises from the finite element discretization of boundary value problems associated to the differential operator I \Gamma grad div. The natural setting for such problems is in the Hilbert space H(div) and the variational fo ..."
Abstract

Cited by 32 (4 self)
 Add to MetaCart
. We consider the solution of the system of linear algebraic equations which arises from the finite element discretization of boundary value problems associated to the differential operator I \Gamma grad div. The natural setting for such problems is in the Hilbert space H(div) and the variational formulation is based on the inner product in H(div). We show how to construct preconditioners for these equations using both domain decomposition and multigrid techniques. These preconditioners are shown to be spectrally equivalent to the inverse of the operator. As a consequence, they may be used to precondition iterative methods so that any given error reduction may be achieved in a finite number of iterations, with the number independent of the mesh discretization. We describe applications of these results to the efficient solution of mixed and least squares finite element approximations of elliptic boundary value problems. 1. Introduction The Hilbert space H(div) consists of squareintegr...
Preconditioning discrete approximations of the ReissnerMindlin plate model
 Math. Modelling Numer. Anal
, 1997
"... Abstract. We consider iterative methods for the solution of the linear system of equations arising from the mixed finite element discretization of the Reissner–Mindlin plate model. We show how to construct a symmetric positive definite block diagonal preconditioner such that the resulting linear sys ..."
Abstract

Cited by 28 (10 self)
 Add to MetaCart
Abstract. We consider iterative methods for the solution of the linear system of equations arising from the mixed finite element discretization of the Reissner–Mindlin plate model. We show how to construct a symmetric positive definite block diagonal preconditioner such that the resulting linear system has spectral condition number independent of both the mesh size h and the plate thickness t. We further discuss how this preconditioner may be implemented and then apply it to efficiently solve this indefinite linear system. Although the mixed formulation of the Reissner–Mindlin problem has a saddlepoint structure common to other mixed variational problems, the presence of the small parameter t and the fact that the matrix in the upper left corner of the partition is only positive semidefinite introduces new complications.