A Two-Grid Block-Centered Finite Difference Algorithm for Nonlinear Compressible Darcy–Forchheimer Model in Porous Media

In this paper, a block-centered finite difference method is proposed to discretize the compressible Darcy–Forchheimer model which describes the high speed non-Darcy flow in porous media. The discretized nonlinear problem on the fine grid is solved by a two-grid algorithm in two steps: first solving a small nonlinear system on the coarse grid, and then solving a nonlinear problem on the fine grid. On the coarse grid, the coupled term of pressure and velocity is approximated by using the fewest number of node values to construct a nonlinear block-centered finite difference scheme. On the fine grid, the original nonlinear term is modified with a small parameter $$\varepsilon $$ε to construct a linear block-centered finite difference scheme. Optimal order error estimates for pressure and velocity are obtained in discrete $$l^\infty (L^2)$$l∞(L2) and $$l^2(L^2)$$l2(L2) norms, respectively. The two-grid block-centered finite difference scheme is proved to be unconditionally convergent without any time step restriction. Some numerical examples are given to testify the accuracy of the proposed method. The numbers of iterations are reported to illustrate the efficiency of the two-grid algorithm.

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

[2]  Thinh Kieu,et al.  Analysis of expanded mixed finite element methods for the generalized forchheimer flows of slightly compressible fluids , 2016 .

[3]  A. Weiser,et al.  On convergence of block-centered finite differences for elliptic-problems , 1988 .

[4]  Hongxing Rui,et al.  A block‐centered finite difference method for Darcy–Forchheimer model with variable Forchheimer number , 2015 .

[5]  Eugenio Aulisa,et al.  Analysis of generalized Forchheimer flows of compressible fluids in porous media , 2009 .

[6]  Yunqing Huang,et al.  A two‐grid method for expanded mixed finite‐element solution of semilinear reaction–diffusion equations , 2003 .

[7]  Yinnian He,et al.  Convergence of three iterative methods based on the finite element discretization for the stationary Navier–Stokes equations☆ , 2009 .

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

[9]  Akif Ibragimov,et al.  Qualitative study of generalized Forchheimer flows with the flux boundary condition , 2012, Advances in Differential Equations.

[10]  Mary F. Wheeler,et al.  A multipoint flux mixed finite element method on distorted quadrilaterals and hexahedra , 2011, Numerische Mathematik.

[11]  Akif Ibragimov,et al.  Structural stability of generalized Forchheimer equations for compressible fluids in porous media , 2010 .

[12]  Kuldeep Chaudhary,et al.  The role of eddies inside pores in the transition from Darcy to Forchheimer flows , 2011 .

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

[14]  Mary F. Wheeler,et al.  Numerical discretization of a Darcy–Forchheimer model , 2008, Numerische Mathematik.

[15]  Mary F. Wheeler,et al.  A Two-Grid Finite Difference Scheme for Nonlinear Parabolic Equations , 1998 .

[16]  Victor Ginting,et al.  Two-grid finite volume element method for linear and nonlinear elliptic problems , 2007, Numerische Mathematik.

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

[18]  Hongxing Rui,et al.  Block-centered finite difference methods for parabolic equation with time-dependent coefficient , 2013, Japan Journal of Industrial and Applied Mathematics.

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

[20]  Peter P. Valko,et al.  MATHEMATICAL FRAMEWORK OF THE WELL PRODUCTIVITY INDEX FOR FAST FORCHHEIMER (NON-DARCY) FLOWS IN POROUS MEDIA , 2009 .

[21]  Wei Liu,et al.  A two-grid algorithm based on expanded mixed element discretizations for strongly nonlinear elliptic equations , 2014, Numerical Algorithms.

[22]  Yanping Chen,et al.  Two-Grid Method for Nonlinear Reaction-Diffusion Equations by Mixed Finite Element Methods , 2011, J. Sci. Comput..

[23]  Weiwei Sun,et al.  Unconditional Convergence and Optimal Error Estimates of a Galerkin-Mixed FEM for Incompressible Miscible Flow in Porous Media , 2012, SIAM J. Numer. Anal..

[24]  Hongxing Rui,et al.  A Block-Centered Finite Difference Method for Slightly Compressible Darcy–Forchheimer Flow in Porous Media , 2017, Journal of Scientific Computing.

[25]  Douglas Ruth,et al.  On the derivation of the Forchheimer equation by means of the averaging theorem , 1992 .

[26]  Yinnian He,et al.  A simplified two-level method for the steady Navier–Stokes equations , 2008 .

[27]  X. Hu,et al.  Two-Grid Methods for Maxwell Eigenvalue Problems , 2014, SIAM J. Numer. Anal..

[28]  Eun-Jae Park.  Mixed finite element methods for generalized Forchheimer flow in porous media , 2005 .

[29]  Mary F. Wheeler,et al.  A Multipoint Flux Mixed Finite Element Method , 2006, SIAM J. Numer. Anal..