A numerical simulation of the newly derived Forchheimer's diffusivity equation has been achieved, careful selection of the non-Darcy coefficient (β) and a fruitful discussion on the use of this coefficient has been presented. Although the proposed model has been designed for single phase flow in porous media, with minor modification it can be applied to multi-phase flow cases. Results show that the proposed numerical model is valid for single-phase flow and is functioning in all ranges of flow in porous medium. The point at which the flow breaks from the Darcian trend has been determined and a new dimensionless term, "Be", has been suggested to define the point of deviating from Darcian flow. "Be" has been found to be 0.0526 at 5% deviation from the Darcian behavior.

An experimental model has been designed to examine certain critical parameters experimentally for the purpose of comparing results with the numerical model predictions. A horizontal fracture has been induced to a homogeneous synthetic sample in the direction of flow with changing aperture to study the effect of fractures on flow behavior.

Experimental results have been compared to the numerical model predictions; a satisfactory agreement within the domain of testing has been concluded which in turn encourages the implementation of this model on a field scale.

Unlike other approaches, the non-Darcy coefficient "β" has been determined experimentally which reflects the authors' belief that when possible, "β" should be evaluated based on the same set of data used to determine permeability in the laboratory.

You can access this article if you purchase or spend a download.