On a Two-Level Parallel MIC(0) Preconditioning of Crouzeix-Raviart Non-conforming FEM Systems
نویسندگان
چکیده
In this paper we analyze a two-level preconditioner for finite element systems arising in approximations of second order elliptic boundary value problems by Crouzeix-Raviart non-conforming triangular linear elements. This study is focused on the efficient implementation of the modified incomplete LU factorization MIC(0) as a preconditioner in the PCG iterative method for the linear algebraic system. A special attention is given to the implementation of the method as a scalable parallel algorithm.
منابع مشابه
MIC(0) Preconditioning of Conforming and Non-conforming FEM Problems on Unstructured Grids⋆
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 the M...
متن کاملGeneralized Aggregation-Based Multilevel Preconditioning of Crouzeix-Raviart FEM Elliptic Problems
Preconditioners based on various multilevel extensions of two-level finite element methods (FEM) are well-known to yield iterative methods of optimal order complexity with respect to the size of the system, as was first shown by Axelsson and Vassilevski [4]. The derivation of optimal convergence rate estimates in this context is mainly governed by the constant γ ∈ (0, 1) in the so-called Cauchy...
متن کاملMathematical Modelling of Flow over Periodic Structures
We examine the influence of roofs’ shapes on the boundary layer of a simplified urban canopy by computing non-stationary Navier-Stokes flow over a periodic pattern. The solution is obtained by means of finite element method (FEM). We use non-conforming Crouzeix-Raviart elements for velocity and piecewise constant elements for pressure. The resulting linear system is solved by the multigrid meth...
متن کامل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 s...
متن کاملLocally Optimized MIC(0) Preconditioning of Rannacher-Turek FEM Systems
In this paper Rannacher-Turek non-conforming rotated bilinear finite elements are applied for the numerical solution of second order elliptic boundary value problems. The preconditioned conjugate gradient method is used for the iterative solution of the arising linear algebraic system. A locally optimized construction for an M-matrix approximation of the global stiffness matrix is the first ste...
متن کامل