Abstract:Low order finite element discretizations of the linear elasticity system suffer increasingly from locking effects and ill-conditioning, when the material approaches the incompressible limit, if only the displacement variable are used. Mixed finite elements using both displacement and pressure variables provide a well-known remedy, but they yield larger and indefinite discrete systems for which the design of scalable and efficient iterative solvers is challenging. Two-level overlapping Schwarz preconditioner for the almost incompressible system of linear elasticity, discretized by mixed finite elements with discontinuous pressures, are constructed and analyzed. The preconditioned systems are accelerated either by a GMRES (generalized minimum residual) method applied to the resulting discrete saddle point problem or by a PCG (preconditioned conjugate gradient) method applied to a positive definite, although extremely ill-conditioned, reformulation of the problem obtained by eliminating all pressure variables on the element level. A novel theoretical analysis of the algorithm for the positive definite reformulation is given by extending some earlier results by Dohrmann and Widlund. The main result of the paper is a bound on the condition number of the algorithm which is cubic in the relative overlap and grows logarithmically with the number of elements across individual subdomains but is otherwise independent of the number of subdomains, their diameters and mesh sizes, and the incompressibility of the material and possible discontinuities of the material parameters across the subdomain interfaces. Numerical results in the plane confirm the theory and also indicate that an analogous result should hold for the saddle point formulation, as well as for spectral element discretizations.

Nonuniformly Weighted Schwarz Smoothers for Spectral Element Multigrid

A Pre-splitting Green's Function Based Hybrid Fast Algorithm for Multiscale Problems

Scalable multigrid methods for immersed finite element methods and immersed isogeometric analysis

Two-level Restricted Additive Schwarz preconditioner based on Multiscale Spectral Generalized FEM for Heterogeneous Helmholtz Problems

Mesh Smoothing for the Spectral Element Method

Parallel Smoothers for Matrix-Based Geometric Multigrid Methods on Locally Refined Meshes Using Multicore CPUs and GPUs

Hybrid and Multiplicative Overlapping Schwarz Algorithms with Standard Coarse Spaces for Mixed Linear Elasticity and Stokes Problems

Fast-convergent two-level restricted additive Schwarz methods based on optimal local approximation spaces

Efficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems

An efficient multigrid method with preconditioned smoother for two-dimensional anisotropic space-fractional diffusion equations

Efficient Multigrid Reduction-in-Time for Method-of-Lines Discretizations of Linear Advection

A Space-Time Multigrid Method for Space-Time Finite Element Discretizations of Parabolic and Hyperbolic PDEs

A Mixed Multiscale Spectral Generalized Finite Element Method

Overlapping Schwarz Methods with A Standard Coarse Space for Almost Incompressible Linear Elasticity

Solving the Poisson equation on small aspect ratio domains using unstructured meshes

On Edge Multiscale Space based Hybrid Schwarz Preconditioner for Helmholtz Problems with Large Wavenumbers

Spectral Methods for Multiscale Stochastic Differential Equations

An energy-efficient GMRES–multigrid solver for space-time finite element computation of dynamic poroelasticity

Achieving $h$- and $p$-robust monolithic multigrid solvers for the Stokes equations

Overlapping schwarz algorithms for almost incompressible linear elasticity

A multigrid reduction framework for domains with symmetries