摘要

A block-centered finite difference scheme is introduced to solve the nonlinear Darcy-Forchheimer equation, in which the velocity and pressure can be approximated simultaneously. The second-order error estimates for both pressure and velocity are established on a nonuniform rectangular grid. Numerical experiments using the scheme show the consistency of the convergence rates of our method with the theoretical analysis.