Scalability of a parallel monolithic multilevel solver for poroelasticity

This study investigates a solver for the quasi-static Biot model for soil consolidation. The scheme consists of an extrapolation scheme in time, com- plemented by a scalable monolithic multigrid method for solving the linear systems resulting after spatial discretisation. The key ingredient for the later is a fixed-stress inexact Uzawa smoother that has been suggested and analysed using local Fourier analysis before [8]. The work at hand investigates the parallel properties of the resulting multigrid solver. For a 3D benchmark problem with roughly 400 million degrees of freedom, scalability is demonstrated in a preliminary study on HazelHen. The presented solver framework should be seen as a prototype, and can be extended and generalized, e.g., to non-linear problems easily.