On parallel solution of linear elasticity problems. Part III: higher order finite elements
Journal article, 2013

This is the third part of a trilogy on parallel solution of the linear elasticity problem. We consider the separate displacement ordering for a plain isotropic problem with full Dirichlet boundary conditions. The parallel solution methods presented in the first two parts of the trilogy are here generalised to higher order by using hierarchical finite elements. We discuss node numberings on regular grids for high degree of parallelism and even processor load as well as the problem of stability of the modified incomplete Cholesky factorisations used. Several preconditioning techniques for the conjugate gradient method are studied and compared for quadratic finite elements. Bounds for the condition numbers of the corresponding preconditioning methods are derived, and computer experiments are performed in order to confirm the theory and give recommendations on the choice of method. The parallel implementation is performed by message passing interface.

Author

Ivar Gustafsson

Chalmers, Mathematical Sciences, Mathematics

University of Gothenburg

Gunhild Lindskog

University of Gothenburg

Chalmers, Mathematical Sciences, Mathematics

Numerical Linear Algebra with Applications

1070-5325 (ISSN) 10991506 (eISSN)

Vol. 20 5 869-887

Subject Categories

Mathematics

DOI

10.1002/nla.1852

More information

Created

10/7/2017