Parallel DD-MIC(0) Preconditioning of Rotated Trilinear FEM Elasticity Systems
نویسنده
چکیده
A new parallel preconditioning algorithm for 3D nonconforming FEM elasticity systems is presented. The preconditioner is constructed in two steps. First displacement decomposition of the stiffness matrix is used. Than MIC(0) factorization is applied to a proper auxiliary M-matrix to get an approximate factorization of the obtained blockdiagonal matrix. The auxiliary matrix has a special block structure – its diagonal blocks are diagonal matrices themselves. This allows the solution of the preconditioning system to be performed efficiently in parallel. Estimates for the parallel times, speedups and efficiencies are derived. The performed parallel tests are in total agreement with them. The robustness of the proposed algorithm is confirmed by the presented experiments solving problems with strong coefficient jumps.
منابع مشابه
Scallability Analysis of Parallel Mic(0) Preconditioning Algorithm for 3d Elliptic Problems
Novel parallel algorithms for the solution of large FEM linear systems arising from second order elliptic partial differential equations in 3D are presented. The problem is discretized by rotated trilinear nonconforming Rannacher–Turek finite elements. The resulting symmetric positive definite system of equations Ax = f is solved by the preconditioned conjugate gradient algorithm. The precondit...
متن کاملComparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems
In this study, the topics of grid generation and FEM applications are studied together following their natural synergy. We consider the following three grid generators: Triangle, NETGEN and Gmsh. The quantitative analysis is based on the number of elements/nodes needed to obtain a triangulation of a given domain, satisfying a certain minimal angle condition. After that, the performance of two d...
متن کاملParallel MIC(0) preconditioning of 3D elliptic problems discretized by Rannacher-Turek finite elements
Novel parallel algorithms for the solution of large FEM linear systems arising from second order elliptic partial differential equations in 3D are presented. The problem is discretized by rotated trilinear nonconforming Rannacher–Turek finite elements. The resulting symmetric positive definite system of equations Ax = f is solved by the preconditioned conjugate gradient algorithm. The precondit...
متن کاملParallel MIC(0) Preconditioning for Numerical Upscaling of Anisotropic Linear Elastic Materials
Abstract. Numerical homogenization is used for upscaling of the linear elasticity tensor of strongly heterogeneous microstructures. The implemented 3D algorithm is described in terms of six auxiliary elastic problems for the reference volume element (RVE). Rotated trilinear RannacherTurek finite elements are used for discretization of the involved subproblems. A parallel PCG method is implement...
متن کامل