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 discre...
| Autores: | , , , |
|---|---|
| Tipo de recurso: | artículo |
| Estado: | Versión aceptada para publicación |
| Fecha de publicación: | 2019 |
| País: | España |
| Institución: | Universidad Pública de Navarra |
| Repositorio: | Academica-e. Repositorio Institucional de la Universidad Pública de Navarra |
| OAI Identifier: | oai:academica-e.unavarra.es:2454/35994 |
| Acceso en línea: | https://hdl.handle.net/2454/35994 |
| Access Level: | acceso abierto |
| Palabra clave: | Darcy–Forchheimer Finite volumes Fractured porous media Geometric multigrid |
| Sumario: | 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. Numerical experiments show that the proposed multigrid method is robust with respect to the fracture permeability, the Forchheimer coefficient and the mesh size. The case of several possibly intersecting fractures in a heterogeneous porous medium is also discussed. |
|---|