Inexact and preconditioned Uzawa algorithms for saddle point problems

Variants of the Uzawa algorithm for solving symmetric indefinite linear systems are developed and analyzed. Each step of this algorithm requires the solution of a symmetric positive- definite system of linear equations. It is shown that if this computation is replaced by an approximate solution produced by an arbitrary iterative method, then with relatively modest requirements on the accuracy of the approximate solution, the resulting inexact Uzawa algorithm is convergent, with a convergence rate close to that of the exact algorithm. In addition, it is shown that preconditioning can be used to improve performance. The analysis is illustrated and supplemented using several examples derived from mixed finite element discretization of the Stokes equations.

[1]  J. Gillis,et al.  Matrix Iterative Analysis , 1961 .

[2]  R. P. Kendall,et al.  An Approximate Factorization Procedure for Solving Self-Adjoint Elliptic Difference Equations , 1968 .

[3]  G. Golub,et al.  Linear least squares and quadratic programming , 1969 .

[4]  Louis A. Hageman,et al.  Iterative Solution of Large Linear Systems. , 1971 .

[5]  Isaac Fried.  Bounds on the extremal eigenvalues of the finite element stiffness and mass matrices and their spectral condition number , 1972 .

[6]  I. Babuska.  The finite element method with Lagrangian multipliers , 1973 .

[7]  F. Brezzi.  On the existence, uniqueness and approximation of saddle-point problems arising from lagrangian multipliers , 1974 .

[8]  A. Brandt,et al.  Multigrid Solutions to Elliptic Flow Problems , 1979 .

[9]  Juhani Pitkäranta,et al.  Boundary subspaces for the finite element method with Lagrange multipliers , 1979 .

[10]  P. Raviart,et al.  Finite Element Approximation of the Navier-Stokes Equations , 1979 .

[11]  R. Dembo,et al.  INEXACT NEWTON METHODS , 1982 .

[12]  M. Fortin,et al.  Augmented Lagrangian methods : applications to the numerical solution of boundary-value problems , 1983 .

[13]  M. Fortin,et al.  A stable finite element for the stokes equations , 1984 .

[14]  R. Verfürth.  A combined conjugate gradient - multi-grid algorithm for the numerical solution of the Stokes problem , 1984 .

[15]  F. Musy,et al.  A Fast Solver for the Stokes Equations Using Multigrid with a UZAWA Smoother , 1985 .

[16]  R. Glowinski,et al.  Numerical Methods for Nonlinear Variational Problems , 1985 .

[17]  Jean E. Roberts,et al.  Mixed and hybrid finite element methods , 1987 .

[18]  A. Wathen.  Realistic Eigenvalue Bounds for the Galerkin Mass Matrix , 1987 .

[19]  J. Cahouet,et al.  Some fast 3D finite element solvers for the generalized Stokes problem , 1988 .

[20]  J. Pasciak,et al.  A preconditioning technique for indefinite systems resulting from mixed approximations of elliptic problems , 1988 .

[21]  G. Golub,et al.  The convergence of inexact Chebyshev and Richardson iterative methods for solving linear systems , 1988 .

[22]  R. Bank,et al.  A class of iterative methods for solving saddle point problems , 1989 .

[23]  G. Wittum.  Multi-grid methods for stokes and navier-stokes equations , 1989 .

[24]  W. Queck.  The convergence factor of preconditioned algorithms of the Arrow-Hurwicz type , 1989 .

[25]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[26]  Margaret H. Wright,et al.  Interior methods for constrained optimization , 1992, Acta Numerica.

[27]  Ragnar Winther,et al.  A Preconditioned Iterative Method for Saddlepoint Problems , 1992, SIAM J. Matrix Anal. Appl..

[28]  Anthony T. Patera,et al.  Analysis of Iterative Methods for the Steady and Unsteady Stokes Problem: Application to Spectral Element Discretizations , 1993, SIAM J. Sci. Comput..

[29]  Andrew J. Wathen,et al.  Fast iterative solution of stabilised Stokes systems, part I: using simple diagonal preconditioners , 1993 .

[30]  A. Wathen,et al.  Fast iterative solution of stabilised Stokes systems part II: using general block preconditioners , 1994 .

[31]  David J. Silvester.  Optimal low order finite element methods for incompressible flow , 1994 .