DD-MIC(O) preconditioning of rotated trilinear FEM elasticity systems

I. Georgiev; S. Margenov

I. Georgiev Central Laboratory of Parallel Processing, Bulgarian Academy of Sciences
S. Margenov Central Laboratory of Parallel Processing, Bulgarian Academy of Sciences

Abstract

New results about preconditioning of rotated trilinear nonconforming FEM elasticity systems in the case of mesh anisotropy are presented. The solver of the arising linear system is based on the constructed efficient preconditioner of the coupled stiffness matrix. Displacement decomposition of the stiffness matrix is used as a first step of the algorithm. At the second step, modified incomplete factorization MIC(O) with perturbation is applied to a proper auxiliary M-matrix to get an approximate factorization of the obtained block-diagonal matrix. The derived condition number estimates and the presented numerical tests well illustrate the behaviour of the theoretically studied algorithms as well as their robustness for some more realistic benchmark problems.

Keywords

References

[1] O. Axelsson. Iterative Solution Methods. University Press, Cambridge 1996.
[2] O. Axelsson, I. Gustafsson. Iterative methods for the Navier equations of elasticity. Camp. Meth. Appl. Mech. Engin., 15: 241-258, 1978.
[3] R. Beauwens, Y. Notay. Efficient block preconditioning of elasticity problems. Technical Report GANMN 98-01, Universite Libre de Bruxelles, Brussels, Belgium 1998.
[4] R. Blaheta. Displacement decomposition - incomplete factorization preconditioning techniques for linear elasticity problems. Numer. Lin. Alg. Appl., 1: 107-126, 1994.
[5] R.S. Falk. Nonconforming finite element methods for the equations of linear elasticity. Math. Camp., 57: 529-550, 1991.