On the preconditioned conjugate gradient solution of a Stokes problem with Robin-type boundary conditions

We discuss the solution by a preconditioned conjugate gradient algorithm of a Stokes problem with Robin-type boundary conditions. This Stokes problem is encountered when applying an appropriate operator-splitting scheme to the time-discretization of a system modeling the interaction of an incompressible viscous fluid with a deformable thin elastic structure. The main contribution of this Note is the identification of a preconditioner operating in the pressure space that reduces substantially the number of iterations when compared to a conjugate gradient algorithm equipped with the canonical scalar product of L2. The results of numerical experiments show the validity of our approach. To cite this article: R. Glowinski, G. Guidoboni, C. R. Acad. Sci. Paris, Ser. I 347 (2009).

RésuméDans cette Note on étudie la résolution d'un problème de Stokes avec conditions aux limites du type Robin par un algorithme de gradient conjugué préconditionné. On rencontre ce type de problèmes lorsque l'on applique certains schémas de discrétisation en temps par décomposition d'opérateurs à la résolution numérique du système d'équations aux derivées partielles modélisant l'interaction d'un fluide visqueux incompressible avec une structure élastique mince. La contribution principale de cette Note est l'identification d'un opérateur de préconditionnement agissant dans l'espace des pressions (en l'occurrence L2(Ω)) ; cet opérateur réduit considérablement le nombre d'itérations nécesssaires à la convergence, lorsque l'on compare à l'algorithme de gradient conjugué opérant dans L2(Ω) muni de son produit scalaire canonique. Les résultats d'essais numériques confirment la validité de notre approche. Pour citer cet article : R. Glowinski, G. Guidoboni, C. R. Acad. Sci. Paris, Ser. I 347 (2009).