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

DOI10.1016/J.JCP.2019.02.020zbMATH Open1452.76090arXiv1808.08328OpenAlexW2949074005WikidataQ128302464 ScholiaQ128302464MaRDI QIDQ2218596FDOQ2218596

M. S. Joshaghani, Matthew G. Knepley, J. Chang, K. B. Nakshatrala

Publication date: 15 January 2021

Published in: Journal of Computational Physics (Search for Journal in Brave)

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.

Full work available at URL: https://arxiv.org/abs/1808.08328

Recommendations

• Block-partitioned solvers for coupled poromechanics: a unified framework
• Parallel Multiblock Multigrid Algorithms for Poroelastic Models
• scientific article; zbMATH DE number 1819864
• Characteristic block-centered finite difference method for compressible miscible displacement in porous media
• scientific article; zbMATH DE number 1894298
• Block preconditioners for mixed-dimensional discretization of flow in fractured porous media
• A two-grid block-centered finite difference algorithm for nonlinear compressible Darcy-Forchheimer model in porous media
• Robust block preconditioners for poroelasticity
• A stabilized mixed discontinuous Galerkin formulation for double porosity/permeability model
• Block-centered finite difference methods for non-Fickian flow in porous media

zbMATH Keywords

parallel computingiterative solversflow through porous mediafinite element discretizationsdouble porosity/permeability modelH(div) elements

Mathematics Subject Classification ID

Parallel numerical computation (65Y05) Preconditioners for iterative methods (65F08) Numerical solution of discretized equations for boundary value problems involving PDEs (65N22) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Finite element methods applied to problems in fluid mechanics (76M10) Flows in porous media; filtration; seepage (76S05) Software, source code, etc. for problems pertaining to numerical analysis (65-04) Software, source code, etc. for problems pertaining to fluid mechanics (76-04)

Cites Work

• Title not available (Why is that?)
• Firedrake
• TSFC: A Structure-Preserving Form Compiler
• GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
• Parallel iterative methods for sparse linear systems
• Two families of mixed finite elements for second order elliptic problems
• Mixed finite elements in \(\mathbb{R}^3\)
• Discontinuous Galerkin methods. Theory, computation and applications. 1st international symposium on DGM, Newport, RI, USA, May 24--26, 1999
• Title not available (Why is that?)
• Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata]
• Finite element exterior calculus, homological techniques, and applications
• Finite element exterior calculus: from Hodge theory to numerical stability
• Title not available (Why is that?)
• Mixed Finite Element Methods and Applications
• Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods
• Title not available (Why is that?)
• Computational Methods for Multiphase Flows in Porous Media
• Efficient rectangular mixed finite elements in two and three space variables
• A stabilized mixed discontinuous Galerkin method for Darcy flow
• Variational inequality approach to enforcing the non-negative constraint for advection-diffusion equations
• Large-scale optimization-based non-negative computational framework for diffusion equations: parallel implementation and performance studies
• Modeling flow in porous media with double porosity/permeability: a stabilized mixed formulation, error analysis, and numerical solutions
• Mathematical aspects of multi-porosity continua
• Efficient mesh management in firedrake using PETSc DMPlex
• A scalable variational inequality approach for flow through porous media models with pressure-dependent viscosity
• Stabilized mixed finite elements for deformable porous media with double porosity
• An Algorithm for the Optimization of Finite Element Integration Loops
• Automated generation and symbolic manipulation of tensor product finite elements
• A parallel edge orientation algorithm for quadrilateral meshes
• Comparative Study of Finite Element Methods Using the Time-Accuracy-Size(TAS) Spectrum Analysis

Cited In (6)

• A segregated approach for modeling the electrochemistry in the 3-D microstructure of li-ion batteries and its acceleration using block preconditioners
• Multilevel Schur-complement algorithms for scalable parallel reservoir simulation with temperature variation
• A vertex scheme for two-phase flow in heterogeneous media
• Field-split preconditioned active-set reduced-space algorithm for complex black oil reservoir simulation at large-scale
• Multilevel field-split preconditioners with domain decomposition for steady and unsteady flow problems
• Multilevel Constrained Pressure-Temperature Residual Preconditioners for Large-Scale Non-Isothermal Reservoir Simulation via Restricted Additive Schwarz Algorithms

Uses Software

• hypre
• PETSc
• mpi4py
• Firedrake
• PyOP2
• DMPlex
• Python
• TSFC
• FInAT
• COFFEE
• petsc4py
• FIAT

This page was built for publication: Composable block solvers for the four-field double porosity/permeability model

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2218596)