Capillarity - Permeability - Darcy’s Law and the Field Equations of the Flow of Underground Fluids

In 1856 Henry Darcy described in an appendix to his book, Les Fontaines Publiques de la Ville de Dijon, a series of experiments on the downward flow of water through filter sands, whereby it was established that the rate of flow is given by the equation: q= - K(h2- h2)/I1 in which q is the volume of water crossing unit area in unit time, I is the thickness of the sand, h1 and h2 the heights above a reference level of the water in manometers terminated above and below the sand, respectively, and K a factor of proportionality. This relationship, appropriately, soon became known as Darcy's law. Subsequently many separate attempts have been made to give Darcy's empirical expression a more general physical fornzulation, with the result that .so many mutually inconsistent expressions of what is purported to he Darcy's law have appeared in published literature that sight has often been lost of Darcy's own work and of its significance. In the present paper, therefore, it shall be our purpose to reinform ourselves upon what Darcy himself did, and then to determine the meaning of his results when expressed explicitly in terms of the pertinent physical variables involved. This will be done first by the empirical method used by Darcy himself, and then by direct derivation from the Navier-Stokes equation of motion of viscous fluids. We find in this manner that: q = (N4(p/p)[g - (l/p) gradPI = WE, is a physical expression for Darcy's law, which is valid for liquids generally, and for gases at pressures higher than about 20 atmospheres. Here N is a shape factor and d a characteristic length of the pore structure of the solid, p and IL are the density and viscosity of the fluid, a = (Nd2)(pip) is the volume conductivity of the system, and E = [g - (l/p) grad p] is the impelling force per unit mass acting upon the fluid. It is found also that Darcy's law is valid only for flow velocities such that the inertial forces are negligible as compared with those arising from viscosity. In general, three-dimensional space there exist two superposed physical fields: a field of force of charucteristic vector E, and a field of flow of vector q. The force field is more general than the flow field since it has values in all space capable of being occupied by the fluid. So long as the fluid density is constant or is a function of the pressure only, curl E = 0, E = - grad F where F = gz + dp/p The field of flow, independently of the force field, must satisfy the conservation of mom, leading to the equation of continuity div pq = — f ap/at , where f is the porosity, und t is the time. For steady rnotion ap/at = 0. and div pq = 0 .