Next Article in Journal
Assessment of the Ecological Risk from Heavy Metals in the Surface Sediment of River Surma, Bangladesh: Coupled Approach of Monte Carlo Simulation and Multi-Component Statistical Analysis
Previous Article in Journal
Swin Transformer for Complex Coastal Wetland Classification Using the Integration of Sentinel-1 and Sentinel-2 Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Essay

The Relationship between the Darcy and Poiseuille Laws

Dipartimento di Scienze e Tecnologie, University of Sannio, 82100 Benevento, Italy
*
Author to whom correspondence should be addressed.
Water 2022, 14(2), 179; https://doi.org/10.3390/w14020179
Submission received: 15 October 2021 / Revised: 29 December 2021 / Accepted: 4 January 2022 / Published: 10 January 2022
(This article belongs to the Section Hydrogeology)

Abstract

:
The Poiseuille and Darcy laws describe the velocity of groundwater flow under laminar conditions. These laws were deducted empirically in conduit and porous systems, respectively, and are widely used to model the groundwater flow. The analytical relationship between these hydraulic laws has been found by draining a tank-reservoir. Based on equations found, the discharge in a conduit under the Poiseuille law can be transformed in the same amount flowing inside a darcian system, and vice versa. This transformation occurs, for example, in karst aquifers, from the matrix to karst conduits during discharge phases, and from conduits to matrix during recharge phases.

1. Introduction

Under laminar condition of the flow, the most common hydraulic law used to model groundwater flow in saturated media are the Darcy and Poiseuille laws.
The so-called Darcy law comes from the French hydraulic engineer Henry Darcy, who investigated the flow of water through horizontal beds of sand to be used for water infiltration [1]. The general form of the Darcy law expresses the velocity of the groundwater flow, v, in relation to hydraulic head lost, dh, in the length path dl [2]:
v = K d h d l
where K is a hydraulic constant, with dimension of [L × T−1], function of the characteristics of the porous media and of the liquid. The minus sign is needed as the differential term dh is negative; the dimensionless ratio dh/dl is known as hydraulic gradient, i, which controls the flow velocity by a linear law.
Hubbert [3] described as the Equation (1), found empirically, can be derived by the Navier-Stokes equations, providing its physical-based meaning.
The so-called Poiseuille law comes from the French physiologist-physicist Jean Leonard Marie Poiseuille, who investigated the velocity of distilled water flow in small conduits in the 1830s [4]. Initially, he planned to maintain a constant temperature of 10 °C; later, he examined the influence of a wide range of water temperature and extended his studies to a great variety of other liquids [5]. During the 1880s, the German hydraulic engineer Gotthilf Heinrich Ludwig Hagen reached similar results to those of Poiseuille on the flow of water in cylindrical tubes; for this reason, the Poiseuille law is also known as Poiseuille-Hagen law. In their review on the Poiseuille law, Sutera and Skalak [6] stated that the results of Hagen were less extensive and less accurate than Poiseuille’s.
The general form of the Poiseuille law expresses the velocity of a liquid, v, in a tube with constant diameter, d, length L, in relation to the difference of the liquid pressure between the ends of the tube, ΔP:
v = d 2 · Δ P 32 · μ · L
where μ is the dynamic viscosity of the liquid.
Also Equation (2), found empirically, can be derived from physical-based assumptions [7,8].
Both hydraulic laws describe the velocity of groundwater flow in saturated zones, but refer to different physical approaches (conduits and porous systems). To investigate their possible relationship, a simple tank-reservoir is drained by a hole and different physical laws are applied (Figure 1).
The analytical relationship found between these hydraulic laws allows to convert the discharge of a Poiseuille system (conduit medium) in the same amount flowing inside a darcian system (porous medium), and vice versa; scaling processes of the hydraulic sections have been also considered. These changing of the flow conditions occur, for example, in karst aquifers, from the matrix to karst conduits during discharge phases, and from conduits to matrix during recharge phases.

2. Tank-Reservoir Models

In the tank-reservoir of Figure 1, the drainage occurs by a tank-tube, with diameter d and length L, located at the bottom of the tank. The effect of the water viscosity decreases if the diameter of the tube increases; in particular, if the diameter is sufficiently wide, the water viscosity can be ignored and drainage occurs in terms of the Torricelli efflux velocity [9].
If the diameter of the tube is small, the water viscosity effect cannot be ignored, and drainage can be analyzed by the Poiseuille law, which assumes that the flow is laminar. Filling the tank-tube with sand, the drainage can be modelled using the Darcy law, which also assumes a laminar flow condition.
The above different tank-drainage systems have been called Torricelli, Poiseuille and Darcy reservoir, respectively [9].
The discharge-time plot for these different reservoirs has been analytically provided by Fiorillo [9], starting from an initial water height in the reservoir, h0, at time t0 = 0, and applying the principle of mass and volume conservation. Thus, water has been considered as an incompressible liquid; a compressible case [10] would require a different and complex approach.
The water pressure in the tube is directly controlled by the water height, h, in the reservoir; besides, it has been assumed null at the end of tube when the drainage starts at time t0.
For the Poiseuille reservoir, the following equations describe the discharge, Q, from the tank-tube with diameter dP [9]:
Q ( t ) = Q 0 · e α P · t
where:
α P = π · d p 4 · ρ · g 128 · μ · L · A 1
Q 0 = π · d p 4 · ρ · g · h 0 128 · μ · L
The parameter α P has a dimension of [T−1] and expresses the drainage characteristics of the Poiseuille reservoir; the gravity acceleration, g, and the water density, ρ, appear from the substitution of the term ΔP of Equation (2) by the ρ⋅g⋅h.
If the tank-tube, with section area AD (and diameter dD) and length L, is filled by sand with hydraulic conductivity K, the Darcy law can be used to simulate the drainage from the tank-reservoir. The following formula describes the discharge, Q, from the tank [9]:
Q ( t ) = Q 0 · e α D · t
where
Q 0 = K · A D · h 0 L
α D = K · A D L · A 1
The parameter α D has a dimension of [T−1] and expresses the drainage characteristics of the Darcy reservoir.
Equations (3) and (6) are that of Maillet [11], and the constant α is a recession coefficient, which controls the exponential decay of the discharge by the time. The recession coefficient depends on the geometric and hydraulic characteristics of the reservoir (Equations (4) and (8)); it is also the constant linking the discharge, Q(t), to the water volume stored V(t) in a linear reservoir:
Q ( t ) = α · V ( t )

3. Comparison between Darcy and Poiseuille Laws

The fixed geometric and hydraulic conditions of the tank-reservoir of Figure 1 allow for the comparison of the Darcy and Poiseuille laws.
Figure 2 shows an example of the discharge-time plot using the Darcy and the Poiseuille laws with the same recession coefficient ( α D = α P = 0.025 d−1).
Three exponential decays are shown in Figure 2, and correspond to different initial water level heights, h0; the three different curves have the same recession coefficients (Equations (4) and (8)), and would provide straight lines in the semilogarithmic plot with the same slope, which is precisely the recession coefficient [9].
On the right-hand axis of Figure 2 is shown the Reynolds number, Re, computed for the flow in the tube of the Poiseuille reservoir by:
R e = ρ · d P · v μ
For the flow in the tube filled by sand (Darcy reservoir), the Reynolds number has been computed by [12]:
R e = ρ · d 50 · q μ
where d50 is the median grain size of the porous medium as a proxy for the characteristic pore length, and q is the specific discharge (q = QD/AD). In the example of Figure 2, for a value of d50 = 0.02 mm, compatible with the hydraulic conductivity K = 1 × 10−4 m/s of the sand, values of the Reynolds number (Equation (11)) are Re << 1.
Imposing α D = α P of Equations (4) and (8), the relationship between diameter dD and dP can be obtained:
d D 2 = d P 4 · g · ρ 32 · K · µ = d P 4 32 · k i
with ki is the intrinsic permeability of the material, ki = (μ·K)/(ρ·g) with dimension [L2].
The Equation (12) provides the link between the Darcy and Poiseuille laws; it can be written also in terms of section areas:
A D = A P 2 · g · ρ 8 π · K · µ = A P 2 8 π · k i
Of course, the diameter of the Darcy tube, dD, cannot be minor than diameter of Poiseuille tube, dP, as the first is filled by sand which reduces the actual section of the tube. Imposing as the limit condition the same diameter for dD and dP (dDdP), Equations (12) and (13) reduce in:
d D d P = 32 · K · μ g · ρ = 32 · k i
A D A P = 8 π · K · μ g · ρ = 8 π · k i
The Equation (15) would connect the size of the area of the Poiseuille conduit, AP, and the intrinsic permeability of a porous medium, ki.
Figure 3 shows the relation between the dP and dD (Equation (12)) for different hydraulic conductivity, K. The grey field characterizes the zones with dPdD, and does not have physical meaning, and it has been delimited by Equation (14).

4. Discussion

In a smoothed tube, the Poiseuille law can be used for a Reynolds number below a threshold value, Re < 2000 [13]; above this value, the flow progressively changes in the turbulent type.
In the porous medium, the turbulent flow occurs progressively from the laminar condition when the velocity of the flow increases, and it is strictly connected to the structure of the medium. Thus, a threshold cannot be fixed, and values Re < 1 could be required to guaranty the laminar condition of the flow [12].
For the Poiseuille law, since the discharge is proportional to the fourth power of the diameter of the tube, large tubes are very much more conductive than smaller ones. In Figure 2, the single tube with diameter of 2 mm could be replaced by 10,000 tubes with diameter of 0.2 mm to maintain the same drainage of the Poiseuille reservoir; this factor dominates the processes of the drainage and the cave patter construction in karst media [14]. For the Darcy reservoir of Figure 2, the single tube with diameter of 0.2 m could be replaced by 100 tubes with diameter of 0.02 m to maintain the same drainage of the reservoir.
In particular, in the Darcy reservoir there is a linear relationship between the discharge, QD, and the area of the tube, AD, for a fixed hydraulic conductivity, K; that is:
Q D A D
In the Poiseuille reservoir the relationship between the discharge, QP, and the area of the tube, AP, is quadratic, that is:
Q P A P 2
In term of diameters the above relationship they are quadratic and to the fourth, respectively:
Q D d D 2 ;   Q P d P 4
If several small tubes in number n, with same diameter dPn or dDn replace the larger tube dP or dD, for the Poiseuille or Darcy reservoir, respectively, the following relationships are obtained:
d P 4 = n · d P n 4
d D 2 = n · d D n 2
For example, in the Poiseuille reservoir, n = 104 tubes with diameter dPn = 0.0001 m discharge as one tube with diameter dP = 0.001 m. In the Darcy reservoir, n = 102 tubes with diameter dDn = 0.0001 m discharge as one tube with diameter dP = 0.001 m.
The different behavior of the two reservoirs highlighted by the Equations (19) and (20), are only apparent. From Equation (13), considering the relationship:
8 π · k i · A D = A P 2
it appears that there are two contributions areas for the Darcian flow: the section area, AD, and the intrinsic permeability, ki. Thus, the reason of why in the Darcy reservoir the discharge is proportional to the area of the tube, AD (Equation (16)) lies in the fact that the intrinsic permeability, ki, is considered constant. Instead, the Equation (15) would provide the intimate relationship between the size of AP and ki.
The above relationships (Equations (19) and (20)) are shown in Figure 4.
Figure 4 shows how the number, n, of smaller tubes (dPn; dDn) rapidly increase as their diameter decreases. Of course, as a corollary, a single tube can drain the discharge of many smaller tubes. These aspects have an important role in the groundwater flow, both in porous and fractured/conduit systems, with particular emphasis in the karst systems, where a complex conduit network is “immersed” in a low permeability fractured limestone volume [15]. Here the karstification processes lead to a hierarchical conduit network into the aquifer [16], which cause the drainage to be converged to specific points: the karst springs. Therefore, the drainage of a karst spring appears “similar” with that of the venous system that leads the blood flow to the heart from the more peripheral areas. In these systems, a single conduit could require a huge number of smaller conduits to be supplyed.
Table 1 shows the relationship between the conduit diameter, dP, and the conduit number, n, following Equation (19); the total volume occupied by conduits, V, rapidly increases when dP decreases. That is, a large volume of water into the smallest conduits is needed to maintain the discharge of a few wide conduits. For the smallest values of dP, the conduits system could be replaced by a porous medium, where the hydraulic conductivity, K, has been estimated by the Equation (14).
The volume, V, estimated in Table 1, represents the complex of total voids volume in the system, both of conduits net and matrix porosity (which replaces conduit porosity for smallest conduits); its ratio to rock volume defines the porosity of the system. Figure 5 shows the variation of the porosity of the system, ne, in relation to dP, for different rock volumes. The porosity, ne, increases as the conduits diameter decreases, and it decreases when the rock volume increases. To support the discharge of wide conduits by the smallest conduits or matrix, a large rock volume would be required; furthermore, in the smallest conduits or matrix, the water volume is huge compared to that flowing in the wider conduits. The matrix of karst aquifers occupies a zone characterized by the smallest conduits, approximately by dP < 0.0001 mm (or K < 2.3 × 10−3 m/s) and range of porosity of 0.3 < ne < 0.001.

5. Conclusions

The Darcy law has a wide use in groundwater hydrology not only in the porous medium, but also in fractured/conduit systems, where the hydraulic conductivity, K, is considered as a mean value of a wide portion of rock masses. The Poiseuille law describes specifically the flow velocity in the conduit systems, where the diameter of conduits controls the articulation of the flow inside the rock mass.
Based on the tank-reservoir drainage, the link between these two hydraulic laws has been discussed. The comparison of the Darcy law with the Poiseuille law highlights how the discharge in a porous medium depends on two areas: the section area interested by the flow (AD for the Darcy reservoir) and by intrinsic permeability, ki, of the medium. In the Poiseuille law, the flow depends directly on the square of the section area interested by the flow ( A P 2 for the Poiseuille reservoir). Even if these characteristics were already known individually for each hydraulic law, their analytical relationship has been poorly described.
Theoretically, it is possible to switch from a hydraulic law to the other using the Equations (12) and (14), or Equations (13) and (15); the Equations (19) and (20) could help in scaling (or vice versa) the conduit diameter in a net system.
All relationships found are valid under laminar conditions of the flow, which almost always occur in porous aquifers; in the discontinuities of fractured aquifers and in the conduit systems of karst aquifers the Poiseuille law could have a limitation due to the presence of a turbulent flow. In any case, its fourth power law highlights well the natural hierarchization observed in many natural systems, as, for example, from the quasi-porous system of the matrix to the final conduit of karst springs.

Author Contributions

Conceptualization, F.F.; methodology, F.F.; software, F.F., L.E., G.L. and M.P.; writing—original draft preparation, F.F., L.E., G.L. and M.P.; writing—review and editing, F.F., L.E., G.L. and M.P.; supervision, F.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Darcy, H. Les Fontaines Publiques de la Ville de Dijon; Victor Dalmont: Paris, France, 1856. [Google Scholar]
  2. Todd, D.K.; Mays, L.W. Groundwater Hydrology, 3rd ed.; Wiley: Chichester, UK, 2005; 636p. [Google Scholar]
  3. Hubbet, M.K. Darcy’s law and the field equations of the flow of underground fluids. Hydrol. Sci. J. 1957, 2, 23–59. [Google Scholar] [CrossRef] [Green Version]
  4. Poiseuille, J.L.M. Ecoulement des Liquides: Societe Phi/omatique de Paris. In Extraits des Proces-Verbaux des Seances Pendant I’Annee 1838; Rene et Cie: Paris, France, 1838; pp. 1–3, 77–81. [Google Scholar]
  5. Poiseuille, J.L.M. Sur Ie mouvement des liquides de nature differente dans les tubes de tres petits diametres. Ann. Chim. Phys. 1847, 3, 76–110. [Google Scholar]
  6. Sutera, S.; Skalak, R. The hystort of Poiseuille’s law. Annu. Rev. Fluid Mech. 1993, 25, 1–19. [Google Scholar] [CrossRef]
  7. Stokes, G.G. On the theories of the internal friction of fluids in motion, and of the equilibrium and motion of elastic solids. Trans. Camb. Phil. Soc. 1845, 8, 287–341. [Google Scholar]
  8. Vennard, J.K. Elementary Fluid Mechanics; John Wiley and Sons: London, UK, 1940; 372p. [Google Scholar]
  9. Fiorillo, F. Tank-reservoir drainage as a simulation of the recession limb of karst spring hydrographs. Hydrogeol. J. 2011, 19, 1009–1019. [Google Scholar] [CrossRef]
  10. Shepherd, J.; Begeal, D.R. Transient Compressible Flow in Porous Materials; United States Department of Energy: Albuquerque, NM, USA, 1988.
  11. Maillet, E. Essais d’Hydraulique Souterraine et Fluviale [Underground and River Hydrology]; Hermann: Paris, France, 1905; 218p. [Google Scholar]
  12. Bear, J. Dynamics of Fluids in Porous Media; Dover Publications: New York, NY, USA, 1972. [Google Scholar]
  13. Colebrook, C. Turbulent Flow in Pipes, with Particular Reference to the Transition Region between the Smooth and Rough Pipe Laws. J. Inst. Civ. Eng. 1939, 11, 133–156. [Google Scholar] [CrossRef]
  14. Ford, D.; Williams, P. Karst Hydrogeology and Geomorphology; Wiley: Chichester, UK, 2007; 562p. [Google Scholar]
  15. Kiraly, L. Karstification and groundwater flow. In Evolution of Karst: From Prekarst to Cessation; Gabrovsek, F., Ed.; Zalozba ZRC: Postojna-Ljubljana, Slovenia, 2002; pp. 155–190. [Google Scholar]
  16. Kaufmann, G. A model comparison of karst aquifer evolution for different matrix flow formulations. J. Hydrol. 2003, 283, 281–289. [Google Scholar] [CrossRef]
Figure 1. Cylindrical tank-reservoir, with a base A1. Water discharge occurs through a bottom tank-tube, with diameter d and length L, starting from an initial water height h0.
Figure 1. Cylindrical tank-reservoir, with a base A1. Water discharge occurs through a bottom tank-tube, with diameter d and length L, starting from an initial water height h0.
Water 14 00179 g001
Figure 2. Discharge-time plot for different initial water level, h0, using the Darcy reservoir and the Poiseuille reservoir with the same recession coefficient ( α D = α P ). The Reynolds number, Re, refers to the flow into the tube of the Poiseuille reservoir.
Figure 2. Discharge-time plot for different initial water level, h0, using the Darcy reservoir and the Poiseuille reservoir with the same recession coefficient ( α D = α P ). The Reynolds number, Re, refers to the flow into the tube of the Poiseuille reservoir.
Water 14 00179 g002
Figure 3. Relationship between the Poiseuille and Darcy laws; dP, tube diameter of Poiseuille reservoir; dD, tube diameter of Darcy reservoir; K, hydraulic conductivity of sand filling the tube of Darcy reservoir. In the grey field (dPdD) no-relationship exists.
Figure 3. Relationship between the Poiseuille and Darcy laws; dP, tube diameter of Poiseuille reservoir; dD, tube diameter of Darcy reservoir; K, hydraulic conductivity of sand filling the tube of Darcy reservoir. In the grey field (dPdD) no-relationship exists.
Water 14 00179 g003
Figure 4. Number of tubes, n, with the same diameter (dPn; dDn) required to replace a wider tube (dP; dD) for a Poiseuille (a) and Darcy (b) reservoir.
Figure 4. Number of tubes, n, with the same diameter (dPn; dDn) required to replace a wider tube (dP; dD) for a Poiseuille (a) and Darcy (b) reservoir.
Water 14 00179 g004
Figure 5. Relationship between conduit diameter, dP, and the porosity, ne, for different rock volumes (data from Table 1). Hydraulic conductivity, K, has been estimated up to a value of 1 m/s; the zone of karst aquifer matrix is also highlighted.
Figure 5. Relationship between conduit diameter, dP, and the porosity, ne, for different rock volumes (data from Table 1). Hydraulic conductivity, K, has been estimated up to a value of 1 m/s; the zone of karst aquifer matrix is also highlighted.
Water 14 00179 g005
Table 1. Diameter of conduit, dP and relative number of conduits, n, from Equation (19). The volume of all conduits, V, computed by considering constant conduit length, L = 1 m. Hydraulic conductivity, K, and intrinsic permeability, ki, have been estimated by Equation 1 (14).
Table 1. Diameter of conduit, dP and relative number of conduits, n, from Equation (19). The volume of all conduits, V, computed by considering constant conduit length, L = 1 m. Hydraulic conductivity, K, and intrinsic permeability, ki, have been estimated by Equation 1 (14).
dP
m
nV
m3
K
m/s
ki
cm2
117.85 × 10−1--
10−11047.85 × 101--
10−21087.85 × 103--
10−310127.85 × 1052.3 × 10−13.13 × 10−4
10−410167.85 × 1072.3 × 10−33.13 × 10−6
10−510207.85 × 1092.3 × 10−53.13 × 10−8
1 Values of K > 1 m/s (and ki > 10−3 cm2) are outside the natural range of porous medium.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fiorillo, F.; Esposito, L.; Leone, G.; Pagnozzi, M. The Relationship between the Darcy and Poiseuille Laws. Water 2022, 14, 179. https://doi.org/10.3390/w14020179

AMA Style

Fiorillo F, Esposito L, Leone G, Pagnozzi M. The Relationship between the Darcy and Poiseuille Laws. Water. 2022; 14(2):179. https://doi.org/10.3390/w14020179

Chicago/Turabian Style

Fiorillo, Francesco, Libera Esposito, Guido Leone, and Mauro Pagnozzi. 2022. "The Relationship between the Darcy and Poiseuille Laws" Water 14, no. 2: 179. https://doi.org/10.3390/w14020179

APA Style

Fiorillo, F., Esposito, L., Leone, G., & Pagnozzi, M. (2022). The Relationship between the Darcy and Poiseuille Laws. Water, 14(2), 179. https://doi.org/10.3390/w14020179

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop