ISSN:
0271-2091
Keywords:
Stokes equations
;
mixed finite elements
;
stabilization
;
conjugate gradient methods
;
preconditioning
;
Engineering
;
Engineering General
Source:
Wiley InterScience Backfile Collection 1832-2000
Topics:
Mechanical Engineering, Materials Science, Production Engineering, Mining and Metallurgy, Traffic Engineering, Precision Mechanics
Notes:
This paper discusses the influence of the stabilization parameter on the convergence factor of various iterative methods for the solution of the Stokes problem discretized by the so-called locally stabilized Q1-P0 finite element. Our objective is to point out optimal parameters which ensure rapid convergence.The first part of the paper is concerned with the dual formulation of the problem. It gives the theoretical precision and practical developments of our stabilized context Uzawa-type algorithm. We assert that the convergence factor of such a method is majored independently of the mesh size by a function of the stabilization parameter. Moreover, we point out that there exists an optimal value of this parameter that minimizes this upper bound. This gives a theoretical justification of pre-existing numerical results. We show that the optimal parameter can be determined a priori. This is a key point when the method has to be implemented. Finally, we base an interpretation of the iterated penalty method numerical behaviour on some theoretical results about the minimum eigenvalue of the stabilized dual operator. This algorithm involves a penalty parameter and a stabilization parameter and we discuss a strategy for choosing optimal parameters.The mixed formulation of the problem is dealt with in the second part of the paper, which proposes several preconditioned conjugate-gradient-type methods. The indefinite character of the problem makes it intrinsically hard. However, if one chooses a suitable preconditioner, this difficulty is overcome, since the preconditioned operator becomes positive definite. We study the eigenvalue spectrum of the preconditioned operator and thereby the convergence factor of the algorithm. In contrast with the two previous formulations, we show that this convergence factor is majored independently of the stabilization parameter. More precisely, we point out convergence factors comparable with those obtained for Poisson-type problems. Finally, we present a variant of the latter method which uses our so-called macroblock-type preconditioner. A comparison with the simple case of diagonal preconditioning is addressed and the improved performance of the macroblock-type preconditioner is evidenced.Various 2D numerical experiments are given to corroborate the theories presented herein.
Additional Material:
7 Ill.
Type of Medium:
Electronic Resource
URL:
http://dx.doi.org/10.1002/fld.1650201104
Permalink