189 research outputs found
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
Geometric multigrid methods for Darcy-Forchheimer flow in fractured porous media
In this paper, we present a monolithic multigrid method for the efficient
solution of flow problems in fractured porous media. Specifically, we consider
a mixed-dimensional model which couples Darcy flow in the porous matrix with
Forchheimer flow within the fractures. A suitable finite volume discretization
permits to reduce the coupled problem to a system of nonlinear equations with a
saddle point structure. In order to solve this system, we propose a full
approximation scheme (FAS) multigrid solver that appropriately deals with the
mixed-dimensional nature of the problem by using mixed-dimensional smoothing
and inter-grid transfer operators. Remarkably, the nonlinearity is localized in
the fractures, and no coupling between the porous matrix and the fracture
unknowns is needed in the smoothing procedure. Numerical experiments show that
the proposed multigrid method is robust with respect to the fracture
permeability, the Forchheimer coefficient and the mesh size.Comment: arXiv admin note: text overlap with arXiv:1811.0126
Recommended from our members
Non-Darcy flow through porous media
Since its introduction, Darcy's law has been implemented as a mathematical tool that allows simple calculation and prediction of low velocity subsurface flows. However, turbulence, non-isothermal conditions, as well as other factors can create conditions where Darcy's law does not accurately describe the head and velocity distributions within a given porous matrix. Darcy's law has been widely applied to analytical and numerical modeling of fluid flow through porous matrix, regardless of the hydrogeologic setting. This study attempts to quantify the error incurred by these models through simultaneous numerical modeling of the mass continuity equation using Darcy's law as well as Forchheimer's relation. To this end, results from steady-state and transient Darcy-based and Forchheimer-based numerical models are presented in this study.Geological Science
- …