A multiscale flux basis for mortar mixed discretizations of reduced Darcy-Forchheimer fracture models

In this paper, a multiscale flux basis algorithm is developed to efficiently solve a flow problem in fractured porous media. Here, we take into account a mixed-dimensional setting of the discrete fracture matrix model, where the fracture network is represented as lower-dimensional object. We assume the linear Darcy model in the rock matrix and the non-linear Forchheimer model in the fractures. In our formulation, we are able to reformulate the matrix-fracture problem to only the fracture network problem and, therefore, significantly reduce the computational cost. The resulting problem is then a non-linear interface problem that can be solved using a fixed-point or Newton-Krylov methods, which in each iteration require several solves of Robin problems in the surrounding rock matrices. To achieve this, the flux exchange (a linear Robin-to-Neumann co-dimensional mapping) between the porous medium and the fracture network is done offline by pre-computing a multiscale flux basis that consists of the flux response from each degree of freedom on the fracture network. This delivers a conserve for the basis that handles the solutions in the rock matrices for each degree of freedom in the fractures pressure space. Then, any Robin sub-domain problems are replaced by linear combinations of the multiscale flux basis during the interface iteration. The proposed approach is, thus, agnostic to the physical model in the fracture network. Numerical experiments demonstrate the computational gains of pre-computing the flux exchange between the porous medium and the fracture network against standard non-linear domain decomposition approaches

On pore-scale modeling and simulation of reactive transport in 3D geometries

Pore-scale modeling and simulation of reactive flow in porous media has a range of diverse applications, and poses a number of research challenges. It is known that the morphology of a porous medium has significant influence on the local flow rate, which can have a substantial impact on the rate of chemical reactions. While there are a large number of papers and software tools dedicated to simulating either fluid flow in 3D computerized tomography (CT) images or reactive flow using pore-network models, little attention to date has been focused on the pore-scale simulation of sorptive transport in 3D CT images, which is the specific focus of this paper. Here we first present an algorithm for the simulation of such reactive flows directly on images, which is implemented in a sophisticated software package. We then use this software to present numerical results in two resolved geometries, illustrating the importance of pore-scale simulation and the flexibility of our software package.Comment: 15 pages, 6 figure

A Lagrange multiplier method for a Stokes-Biot fluid-poroelastic structure interaction model

We study a finite element computational model for solving the coupled problem arising in the interaction between a free fluid and a fluid in a poroelastic medium. The free fluid is governed by the Stokes equations, while the flow in the poroelastic medium is modeled using the Biot poroelasticity system. Equilibrium and kinematic conditions are imposed on the interface. A mixed Darcy formulation is employed, resulting in continuity of flux condition of essential type. A Lagrange multiplier method is employed to impose weakly this condition. A stability and error analysis is performed for the semi-discrete continuous-in-time and the fully discrete formulations. A series of numerical experiments is presented to confirm the theoretical convergence rates and to study the applicability of the method to modeling physical phenomena and the sensitivity of the model with respect to its parameters

Coupling Biot and Navier-Stokes equations for modelling fluid-poroelastic media interaction

The interaction between a fluid and a poroelastic structure is a complex problem that couples the Navier–Stokes equations with the Biot system. The finite element approximation of this problem is involved due to the fact that both subproblems are indefinite. In this work, we first design residual-based stabilization techniques for the Biot system, motivated by the variational multiscale approach. Then, we state the monolithic Navier–Stokes/Biot system with the appropriate transmission conditions at the interface. For the solution of the coupled system, we adopt both monolithic solvers and heterogeneous domain decomposition strategies. Different domain decomposition methods are considered and their convergence is analyzed for a simplified problem. We compare the efficiency of all the methods on a test problem that exhibits a large added-mass effect, as it happens in hemodynamics applications