Domain decomposition preconditioners for p and hp finite element approximation of Stokes equations

Abstract Domain decomposition preconditioning techniques are developed in the context of hp finite element approximation of the Stokes problem. Two basic types of preconditioner are considered: a block diagonal scheme based on decoupling the velocity and pressure components, and a scheme based on an indefinite system similar to the original Stokes system. For each type of scheme, theoretical estimates are obtained for the location of the eigenvalues of the preconditioned operators in terms of the polynomial degree, the mesh sizes on the coarse and fine grids, and the inf—sup constant for the method. Theoretical estimates show that the growth of the bounds is modest as the mesh is refined and the polynomial order is increased. The preconditioners are shown to be applicable to various iterative schemes for the Stokes systems. The theoretical bounds are compared with actual quantities obtained in practical computations for several representative problems.

[1]  Axel Klawonn,et al.  An Optimal Preconditioner for a Class of Saddle Point Problems with a Penalty Term , 1995, SIAM J. Sci. Comput..

[2]  Jan Mandel,et al.  Iterative solvers by substructuring for the p -version finite element method , 1990 .

[3]  Tarek P. Mathew,et al.  Schwarz alternating and iterative refinement methods for mixed formulations of elliptic problems, part I: Algorithms and numerical results , 1993 .

[4]  Jinchao Xu,et al.  Iterative Methods by Space Decomposition and Subspace Correction , 1992, SIAM Rev..

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

[6]  Apostol T. Vassilev,et al.  Analysis of the Inexact Uzawa Algorithm for Saddle Point Problems , 1997 .

[7]  Ed Anderson,et al.  LAPACK Users' Guide , 1995 .

[8]  R. Stenberg,et al.  Mixed $hp$ finite element methods for problems in elasticity and Stokes flow , 1996 .

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

[10]  Gregory H. Wannier,et al.  A contribution to the hydrodynamics of lubrication , 1950 .

[11]  Mark Ainsworth,et al.  A Preconditioner Based on Domain Decomposition for H-P Finite-Element Approximation on Quasi-Uniform Meshes , 1996 .

[12]  J. Wang,et al.  Analysis of the Schwarz algorithm for mixed finite elements methods , 1992 .

[13]  I. Babuska,et al.  Efficient preconditioning for the p -version finite element method in two dimensions , 1991 .

[14]  J. Pasciak,et al.  A domain decomposition technique for Stokes problems , 1990 .

[15]  Tarek P. Mathew,et al.  Schwarz alternating and iterative refinement methods for mixed formulations of elliptic problems, part II: Convergence theory , 1993 .

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

[17]  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..

[18]  Abani K. Patra,et al.  Non-overlapping domain decomposition methods for adaptive hp approximations of the Stokes problem with discontinuous pressure fields , 1997 .

[19]  Mark Ainsworth.  A Hierarchical Domain Decomposition Preconditioner for h-P Finite Element Approximation on Locally Refined Meshes , 1996, SIAM J. Sci. Comput..

[20]  H. K. Moffatt.  Viscous and resistive eddies near a sharp corner , 1964, Journal of Fluid Mechanics.

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

[22]  Luca F. Pavarino,et al.  Preconditioned conjugate residual methods for mixed spectral discretizations of elasticity and Stokes problems , 1997 .

[23]  Spencer J. Sherwin,et al.  Hierarchical hp finite elements in hybrid domains , 1997 .