A two-grid decoupled algorithm for fracture models

Abstract In this paper, we consider a two-grid decoupled algorithm for two categories of fracture models. In the mixed Darcy/Darcy fracture model, fluid flow in fractures as well as in the surrounding medium is governed by Darcy’s law. And in the mixed Darcy–Forchheimer/Darcy fracture model, flow in fractures is governed by Darcy–Forchheimer’s law while that in surrounding matrix is governed by Darcy’s law. In the proposed two-grid method, we use a coarse grid approximation to the interface coupling conditions for decoupling mixed problems in fractures and surrounding matrix. Error estimates show that the two-grid decoupled algorithm retains the same order of approximation accuracy as the coupled algorithm for the pressure p on H 1 semi-norm and the velocity u on ( L 2 ) 2 norm. Numerical experiments are carried out to verify the accuracy and efficiency of the decoupled algorithm, in which the computation times are reduced greatly compared to the coupled method.

[1]  Jinchao Xu,et al.  Error estimates on a new nonlinear Galerkin method based on two-grid finite elements , 1995 .

[2]  Jinchao Xu,et al.  Numerical Solution to a Mixed Navier-Stokes/Darcy Model by the Two-Grid Approach , 2009, SIAM J. Numer. Anal..

[3]  Jinchao Xu,et al.  A Novel Two-Grid Method for Semilinear Elliptic Equations , 1994, SIAM J. Sci. Comput..

[4]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[5]  Jean E. Roberts,et al.  Modeling fractures as interfaces: a model for Forchheimer fractures , 2008 .

[6]  Hongxing Rui,et al.  A Block-Centered Finite Difference Method for the Darcy-Forchheimer Model , 2012, SIAM J. Numer. Anal..

[7]  Zhengjia Sun,et al.  A block-centered finite difference method for reduced fracture model in Karst aquifer system , 2017, Comput. Math. Appl..

[8]  Jinchao Xu.  Two-grid Discretization Techniques for Linear and Nonlinear PDEs , 1996 .

[9]  Alfio Quarteroni,et al.  Domain Decomposition Methods for Partial Differential Equations , 1999 .

[10]  Wei Liu,et al.  A Two-Grid Block-Centered Finite Difference Method For Darcy-Forchheimer Flow in Porous Media , 2015, SIAM J. Numer. Anal..

[11]  Alessio Fumagalli,et al.  Numerical modelling of multiphase subsurface flow in the presence of fractures , 2011 .

[12]  Peter Knabner,et al.  Mathematical analysis of a discrete fracture model coupling Darcy flow in the matrix with Darcy-Forchheimer flow in the fracture , 2014, ArXiv.

[13]  Jean E. Roberts,et al.  Space-time Domain Decomposition and Mixed Formulation for solving reduced fracture models , 2015 .

[14]  Hongxing Rui,et al.  Mixed Element Method for Two-Dimensional Darcy-Forchheimer Model , 2012, J. Sci. Comput..

[15]  Philippe Angot,et al.  ASYMPTOTIC AND NUMERICAL MODELLING OF FLOWS IN FRACTURED POROUS MEDIA , 2009 .

[16]  Alfio Quarteroni,et al.  A multiscale Darcy–Brinkman model for fluid flow in fractured porous media , 2011, Numerische Mathematik.

[17]  Vincent Martin,et al.  Modeling Fractures and Barriers as Interfaces for Flow in Porous Media , 2005, SIAM J. Sci. Comput..

[18]  Jinchao Xu,et al.  Iterative Methods by Space Decomposition and Subspace Correction , 1992, SIAM Rev..

[19]  Fernando A. Morales,et al.  The narrow fracture approximation by channeled flow , 2010 .

[20]  Jinchao Xu,et al.  A Two-Grid Method of a Mixed Stokes-Darcy Model for Coupling Fluid Flow with Porous Media Flow , 2007, SIAM J. Numer. Anal..

[21]  Vincent Martin,et al.  Modeling fractures as interfaces with nonmatching grids , 2012, Computational Geosciences.