Composable block solvers for the four-field double porosity/permeability model

by   M. S. Joshaghani, et al.

The objective of this paper is twofold. First, we propose two composable block solver methodologies to solve the discrete systems that arise from finite element discretizations of the double porosity/permeability (DPP) model. The DPP model, which is a four-field mathematical model, describes the flow of a single-phase incompressible fluid in a porous medium with two distinct pore-networks and with a possibility of mass transfer between them. Using the composable solvers feature available in PETSc and the finite element libraries available under the Firedrake Project, we illustrate two different ways by which one can effectively precondition these large systems of equations. Second, we employ the recently developed performance model called the Time-Accuracy-Size (TAS) spectrum to demonstrate that the proposed composable block solvers are scalable in both the parallel and algorithmic sense. Moreover, we utilize this spectrum analysis to compare the performance of three different finite element discretizations (classical mixed formulation with H(div) elements, stabilized continuous Galerkin mixed formulation, and stabilized discontinuous Galerkin mixed formulation) for the DPP model. Our performance spectrum analysis demonstrates that the composable block solvers are fine choices for any of these three finite element discretizations. Sample computer codes are provided to illustrate how one can easily implement the proposed block solver methodologies through PETSc command line options.


page 1

page 2

page 3

page 4


A new mixed finite-element method for the biharmonic problem

Fourth-order differential equations play an important role in many appli...

A stabilized mixed discontinuous Galerkin formulation for double porosity/permeability model

Modeling flow through porous media with multiple pore-networks has now b...

Slate: extending Firedrake's domain-specific abstraction to hybridized solvers for geoscience and beyond

Within the finite element community, discontinuous Galerkin (DG) and mix...

A field-split preconditioning technique for fluid-structure interaction problems with applications in biomechanics

We present a novel preconditioning technique for Krylov subspace algorit...

A novel block non-symmetric preconditioner for mixed-hybrid finite-element-based flow simulations

In this work we propose a novel block preconditioner, labelled Explicit ...

An iterative generalized Golub-Kahan algorithm for problems in structural mechanics

This paper studies the Craig variant of the Golub-Kahan bidiagonalizatio...

Fourier analysis of membrane locking and unlocking

Membrane locking in finite element approximations of thin beams and shel...

Please sign up or login with your details

Forgot password? Click here to reset