Composable block solvers for the four-field double porosity/permeability model
arXiv:1808.08328 · doi:10.1016/j.jcp.2019.02.020
Abstract
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.
References in corpus (15)
- Firedrake: automating the finite element method by composing abstractions
- Finite element exterior calculus: from Hodge theory to numerical stability
- TSFC: a structure-preserving form compiler
- Automated generation and symbolic manipulation of tensor product finite elements
- Efficient mesh management in Firedrake using PETSc-DMPlex
- COFFEE: an Optimizing Compiler for Finite Element Local Assembly
- Modeling flow in porous media with double porosity/permeability: A stabilized mixed formulation, error analysis, and numerical solutions
- A parallel edge orientation algorithm for quadrilateral meshes
- Variational inequality approach to enforce the non-negative constraint for advection-diffusion equations
- An algorithm for the optimization of finite element integration loops
- A structure-exploiting numbering algorithm for finite elements on extruded meshes, and its performance evaluation in Firedrake
- A scalable variational inequality approach for flow through porous media models with pressure-dependent viscosity
- Exposing and exploiting structure: optimal code generation for high-order finite element methods
- A stabilized mixed discontinuous Galerkin formulation for double porosity/permeability model
- Large-scale Optimization-based Non-negative Computational Framework for Diffusion Equations: Parallel Implementation and Performance Studies