Monolithic multigrid for a reduced-quadrature discretization of poroelasticity

by   James H. Adler, et al.

Advanced finite-element discretizations and preconditioners for models of poroelasticity have attracted significant attention in recent years. The equations of poroelasticity offer significant challenges in both areas, due to the potentially strong coupling between unknowns in the system, saddle-point structure, and the need to account for wide ranges of parameter values, including limiting behavior such as incompressible elasticity. This paper was motivated by an attempt to develop monolithic multigrid preconditioners for the discretization developed in [48]; we show here why this is a difficult task and, as a result, we modify the discretization in [48] through the use of a reduced quadrature approximation, yielding a more "solver-friendly" discretization. Local Fourier analysis is used to optimize parameters in the resulting monolithic multigrid method, allowing a fair comparison between the performance and costs of methods based on Vanka and Braess-Sarazin relaxation. Numerical results are presented to validate the LFA predictions and demonstrate efficiency of the algorithms. Finally, a comparison to existing block-factorization preconditioners is also given.



There are no comments yet.


page 6


A new iterative method for solving a class of two-by-two block complex linear systems

We present an iterative method for solving the system arisen from finite...

Low-order preconditioning of the Stokes equations

Low-order finite-element discretizations are well-known to provide effec...

A local Fourier analysis of additive Vanka relaxation for the Stokes equations

Multigrid methods are popular solution algorithms for many discretized P...

Optimization of two-level methods for DG discretizations of reaction-diffusion equations

We analyze and optimize two-level methods applied to a symmetric interio...

A novel multigrid method for elliptic distributed control problems

Large linear systems of saddle-point type have arisen in a wide variety ...

A conservative, physically compatible discretization for turbidity currents

The recently introduced Mass Energy Enstrophy and Vorticity conserving (...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.