A Two-Grid Block-Centered Finite Difference Method For Darcy--Forchheimer Flow in Porous Media

Abstract

A two-grid block-centered finite difference method is proposed for solving the two-dimensional Darcy--Forchheimer model describing non-Darcy flow in porous media. To construct the two-grid method we modify the original nonlinear elliptic operator of Darcy--Forchheimer flow to a twice continuously differentiable one by introducing a small and positive parameter $\varepsilon$. By using the two-grid method, solving a nonlinear equation on a fine grid is reduced to solving a nonlinear equation on a coarse grid together with solving a linear equation on a fine grid. Optimal order error estimates for pressure and velocity in discrete $L^2$ norms are obtained. Some numerical examples are given to show the accuracy and efficiency of the presented method.