Numerical modeling of Forchheimer's equation to describe darcy and non-darcy flow in porous media

This paper adds a new dimension to fluid flow in porous media by replacing Darcy's equation with Forchheimer's equation, which in fact takes in consideration both Darcy and non-Darcy flow behaviors. Mathematical derivation of the diffusivity equation based on the Forchheimer equation has been accomplished for linear flow in both one- and two-dimensional cases. Numerical simulation of the new diffusivity equation has been achieved in both cases of Darcy and non-Darcy (Forchheimer) domains. Both numerical models have been tested and verified giving very reasonable accuracy in describing flow characteristics of Darcy and non-Darcy behaviors. Interestingly, a new dimensionless number “Be” relating β (mainly a function of permeability and porosity), velocity, density and viscosity has been introduced to differentiate between Darcy and non-Darcy flow in porous medium for any types of rock and flowing fluid. This new dimensionless number is far from being considered declaration of turbulence flow in porous medium rather the energy loss is contributed to the nature of both flowing fluid and the porous medium. Theoretically, the break point of non-Darcy behavior from Darcy's has been found at (Be = 0), for practical use it has been determined that non-Darcy flow would start at the point when “Be” become greater than zero. A practical range of permeability with porosity changing accordingly, fluid velocity, density, viscosity and the non-Darcy coefficient (β) estimated from different types of available correlations in the literature have been examined using the numerical models obtained. In all cases “Be” remains the same at the point of flow behavior change to non-Darcy.