M. S. Joshaghani

M. S. Joshaghani, J. Chang, K. B. Nakshatrala 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.

K. B. Nakshatrala, M. S. Joshaghani

Many processes in nature (e.g., physical and biogeochemical processes in hyporheic zones, and arterial mass transport) occur near the interface of free-porous media. A firm understanding of these processes needs an accurate prescription of flow dynamics near the interface which (in turn) hinges on an appropriate description of interface conditions along the interface of free-porous media. Although the conditions for the flow dynamics at the interface of free-porous media have received considerable attention, many of these studies were empirical and lacked a firm theoretical underpinning. In this paper, we derive a complete and self-consistent set of conditions for flow dynamics at the interface of free-porous media. We first propose a principle of virtual power by incorporating the virtual power expended at the interface of free-porous media. Then by appealing to the calculus of variations, we obtain a complete set of interface conditions for flows in coupled free-porous media. A noteworthy feature of our approach is that the derived interface conditions apply to a wide variety of porous media models. We also show that the two most popular interface conditions -- the Beavers-Joseph condition and the Beavers-Joseph-Saffman condition -- are special cases of the approach presented in this paper. The proposed principle of virtual power also provides a minimum power theorem for a class of flows in coupled free-porous media, which has a similar mathematical structure as the ones enjoyed by flows in uncoupled free and porous media.

M. S. Joshaghani, S. H. S. Joodat, K. B. Nakshatrala

Modeling flow through porous media with multiple pore-networks has now become an active area of research due to recent technological endeavors like geological carbon sequestration and recovery of hydrocarbons from tight rock formations. Herein, we consider the double porosity/permeability (DPP) model, which describes the flow of a single-phase incompressible fluid through a porous medium exhibiting two dominant pore-networks with a possibility of mass transfer across them. We present a stable mixed discontinuous Galerkin (DG) formulation for the DPP model. The formulation enjoys several attractive features. These include: (i) Equal-order interpolation for all the field variables (which is computationally the most convenient) is stable under the proposed formulation. (ii) The stabilization terms are residual-based, and the stabilization parameters do not contain any mesh-dependent parameters. (iii) The formulation is theoretically shown to be consistent, stable, and hence convergent. (iv) The formulation supports non-conforming discretizations and distorted meshes. (v) The DG formulation has improved element-wise (local) mass balance compared to the corresponding continuous formulation. (vi) The proposed formulation can capture physical instabilities in coupled flow and transport problems under the DPP model.