The mathematical study of flow in porous media is typically based on the 1856 empirical result of Henri Darcy. This result, known as Darcy’s law, states that the velocity of a single-phase flow through a porous medium is proportional to the hydraulic gradient. The publication of Darcy’s work has been referred to as “the birth of groundwater hydrology as a quantitative science” (Freeze and Cherry, 1979). Although Darcy’s original equation was found to be valid for slow, steady, one-dimensional, single-phase flow through a homogeneous and isotropic sand, it has been applied in the succeeding 140 years to complex transient flows that involve multiple phases in heterogeneous media. To attain this generality, a modification has been made to the original formula, such that the constant of proportionality between flow and hydraulic gradient is allowed to be a spatially varying function of the system properties. The extended version of Darcy’s law is expressed in the following form: qα=-Kα . Jα (2.1) where qα is the volumetric flow rate per unit area vector of the α-phase fluid, Kα is the hydraulic conductivity tensor of the α-phase and is a function of the viscosity and saturation of the α-phase and of the solid matrix, and Jα is the vector hydraulic gradient that drives the flow. The quantities Jα and Kα account for pressure and gravitational effects as well as the interactions that occur between adjacent phases. Although this generalization is occasionally criticized for its shortcomings, equation (2.1) is considered today to be a fundamental principle in analysis of porous media flows (e.g., McWhorter and Sunada, 1977). If, indeed, Darcy’s experimental result is the birth of quantitative hydrology, a need still remains to build quantitative analysis of porous media flow on a strong theoretical foundation. The problem of unsaturated flow of water has been attacked using experimental and theoretical tools since the early part of this century. Sposito (1986) attributes the beginnings of the study of soil water flow as a subdiscipline of physics to the fundamental work of Buckingham (1907), which uses a saturation-dependent hydraulic conductivity and a capillary potential for the hydraulic gradient.
A fully implicit finite volume scheme for single phase flow with reactive transport in porous media
