Darcy's law 849543 220770964 2008-06-21T14:22:49Z Mwtoews 711150 rv In [[fluid dynamics]], '''Darcy's law''' is a [[Phenomenology (science)|phenomologically]] derived [[constitutive equation]] that describes the flow of a [[fluid]] through a [[porous medium]]. The law was formulated by [[Henry Darcy]] based on the results of experiments<ref>H. Darcy, Les Fontaines Publiques de la Ville de Dijon, Dalmont, Paris (1856).</ref> on the flow of [[water]] through beds of [[sand]]. It also forms the scientific basis of [[Permeability (fluid)|fluid permeability]] used in the [[earth science]]s. == Background == Although Darcy's law (an expression of [[conservation of momentum]]) was determined experimentally by Darcy, it has since been derived from the [[Navier-Stokes equations]] via [[Homogenization (mathematics)|homogenization]]. It is analogous to [[Fourier's law]] in the field of [[heat conduction]], [[Ohm's law]] in the field of [[electrical networks]], or [[Fick's law]] in [[diffusion]] theory. One application of Darcy's law is to water flow through an [[aquifer]]. Darcy's law along with the equation of [[conservation of mass]] are equivalent to the [[groundwater flow equation]], one of the basic relationships of [[hydrogeology]]. Darcy's law is also used to describe oil, water, and gas flows through petroleum reservoirs. == Description == [[Image:Darcy's Law.png|thumb|right|300px|Diagram showing definitions and directions for Darcy's law.]] Darcy's law is a simple proportional relationship between the instantaneous discharge rate through a porous medium, the viscosity of the fluid and the pressure drop over a given distance. :<math>Q=\frac{-\kappa A}{\mu}\frac{(P_b-P_a)}{L}</math> The total discharge, <math>Q</math> (units of volume per time, e.g., m³/s) is equal to the product of the [[Permeability (fluid)|permeability]] (<math>\kappa</math> units of area, e.g. m²) of the medium, the cross-sectional area (<math>A</math>) to flow, and the pressure drop (<math>P_b-P_a</math>), all divided by the [[dynamic viscosity]] <math>\mu</math> (in SI units e.g. kg/(m·s) or Pa·s), and the length <math>L</math> the pressure drop is taking place over. The negative sign is needed because fluids flow from high pressure to low pressure. So if the change in pressure is negative (in the <math>x</math>-direction) then the flow will be positive (in the <math>x</math>-direction). Dividing both sides of the equation by the area and using more general notation leads to :<math>q=\frac{-\kappa}{\mu}\nabla P</math> where <math>q</math> is the [[flux]] (discharge per unit area, with units of length per time, m/s) and <math>\nabla P</math> is the [[pressure gradient]] vector. This value of flux, often referred to as the Darcy flux, is not the velocity which the water traveling through the pores is experiencing<ref>See {{cite journal | last=Stauffer | first=Philip H. | title=Flux Flummoxed: A Proposal for Consistent Usage | journal=Ground Water | year=2006 | volume=44 | issue=2 | pages=125–128 | doi=10.1111/j.1745-6584.2006.00197.x}} for a discussion of the many, sometimes confusing names given to (<math>q</math>) in the ground water literature.</ref>. The pore velocity (<math>v</math>) is related to the Darcy flux (<math>q</math>) by the [[porosity]] (<math>\phi</math>). The flux is divided by porosity to account for the fact that only a fraction of the total formation volume is available for flow. The pore velocity would be the velocity a conservative tracer would experience if carried by the fluid through the formation. :<math>v=\frac{q}{\phi}</math> ===In 3D=== In [[Three-dimensional space|three dimensions]], [[gravity]] must be accounted for, as the flow is not affected by the vertical pressure drop caused by gravity when assuming [[hydrostatic]] conditions. The solution is to subtract the gravitational pressure drop from the existing pressure drop in order to express the resulting flow, :<math>\boldsymbol{q} = \frac{-\boldsymbol{\kappa}}{\mu}(\nabla P - \rho g \hat{\boldsymbol{e}}_z)</math> where the flux <math>\boldsymbol{q}</math> is now a vector quantity, <math>\boldsymbol{\kappa}</math> is a tensor of permeability, <math>\nabla</math> is the [[gradient]] operator in 3D, <math>g</math> is the [[gravitational acceleration|acceleration due to gravity]], <math>\hat{\boldsymbol{e}}_z</math> is the [[unit vector]] in the vertical direction, pointing downwards and <math>\rho</math> is the [[density]]. Effects of [[anisotropy]] in three dimensions are addressed using a [[Symmetric matrix|symmetric]] second-order [[Permeability (fluid)#Tensor permeability|tensor of permeability]]: :<math>\boldsymbol{\kappa} = \begin{bmatrix} \kappa_{xx} & \kappa_{xy} & \kappa_{xz} \\ \kappa_{yx} & \kappa_{yy} & \kappa_{yz} \\ \kappa_{zx} & \kappa_{zy} & \kappa_{zz} \end{bmatrix}</math> where the magnitudes of permeability in the ''x'', ''y'', and ''z'' component directions are specified. Since this a symmetric matrix, there are ''at most'' six unique values. If the permeability is isotropic (equal magnitude in all directions), then the [[Identity matrix|diagonal values]] are equal, <math>\kappa_{xx} = \kappa_{yy} = \kappa_{zz} > 0\,</math>, while all other components are 0. The permeability tensor can be interpreted through an evaluation the relative magnitudes in each component. For example, rock with highly permeable vertical [[Fracture (geology)|fractures]] aligned in the ''x''-direction will have relatively higher values for <math>\kappa_{xx}\,</math> than other component values. ===Assumptions=== Darcy's law is a simple mathematical statement which neatly summarizes several familiar properties that [[groundwater]] flowing in [[aquifer]]s exhibits, including: * if there is no pressure gradient over a distance, no flow occurs (this is [[hydrostatics|hydrostatic]] conditions), * if there is a pressure gradient, flow will occur from high pressure towards low pressure (opposite the direction of increasing gradient - hence the negative sign in Darcy's law), * the greater the pressure gradient (through the same formation material), the greater the discharge rate, and * the discharge rate of fluid will often be different &mdash; through different formation materials (or even through the same material, in a different direction) &mdash; even if the same pressure gradient exists in both cases. A graphical illustration of the use of the steady-state [[groundwater flow equation]] (based on Darcy's law and the conservation of mass) is in the construction of [[flownet]]s, to quantify the amount of [[groundwater]] flowing under a [[dam]]. Darcy's law is only valid for slow, [[viscous]] flow; fortunately, most groundwater flow cases fall in this category. Typically any flow with a [[Reynolds number]] less than one is clearly laminar, and it would be valid to apply Darcy's law. Experimental tests have shown that for flow regimes with values of Reynolds number up to 10 may still be Darcian. Reynolds number (a dimensionless parameter) for porous media flow is typically expressed as : <math>Re = \frac{\rho v d_{30}}{\mu}</math> where <math>\rho</math> is the [[density]] of the [[fluid]] (units of mass per volume), <math>v</math> is the specific discharge (not the pore velocity &mdash; with units of length per time), <math>d_{30}</math> is a representative grain diameter for the porous medium (often taken as the 30% passing size from a [[Particle size|grain size]] analysis using sieves), and <math>\mu</math> is the dynamic viscosity of the fluid. == Derivation == Assuming stationary, creeping, incompressible flow, the Navier-Stokes equation becomes : <math> \mu\nabla^2 u_i +\rho g_i -\partial_i p=0</math>, where <math>\mu</math> is the viscosity, <math>u_i</math> is the velocity in the i direction, <math>g_i</math> is the gravity component in the i direction and p is the pressure. Assuming the viscous resisting force is proportional to the velocity, and opposite in direction. we may write <math>\frac{\mu\phi}{k_i}u_i+\rho g_i-\partial_i p</math>, where <math>\phi</math> is the porosity.This gives the velocity <math>u_i=-\frac{k_i}{\phi\mu}\left(\partial_i p-\rho g_i\right)</math>, which gives Darcy's law <math>q_i=-\frac{k_i}{\mu}\left(\partial_i p -\rho g_i\right)</math>. == Additional forms of Darcy's law == ===Time derivative of flux=== For very short time scales, a time derivative of flux may be added to Darcy's law, which results in valid solutions at very small times (in heat transfer, this is called the modified form of [[Fourier's law]]), :<math>\tau\frac{\partial q}{\partial t}+q=-K\nabla h,</math> where <math>\tau</math> is a very small time constant which causes this equation to reduce to the normal form of Darcy's law at "normal" times (>&nbsp;[[nanosecond]]s). The main reason for doing this is that the regular [[groundwater flow equation]] ([[diffusion equation]]) leads to [[Mathematical singularity|singularities]] at constant head boundaries at very small times. This form is more mathematically rigorous, but leads to a [[hyperbolic]] groundwater flow equation, which is more difficult to solve and is only useful at very small times, typically out of the realm of practical use. ===Brinkman term=== Another extension to the traditional form of Darcy's law is the Brinkman term, which is used to account for transitional flow between boundaries (introduced by Brinkman in 1947), :<math>\beta\nabla^{2}q+q=-K\nabla h,</math> where <math>\beta</math> is an [[effective viscosity]] term. This correction term accounts for flow through medium where the grains of the media are porous themselves, but is difficult to use, and is typically neglected. ===Multiphase flow=== For [[multiphase flow]], an approximation is to use Darcy's law for each phase, with permeability replaced by phase permeability, which is the permeability of the rock multiplied with [[relative permeability]]. This approximation is valid if the interfaces between the fluids remain static, which is not true in general, but it is still a reasonable model under steady-state conditions. Assuming that the flow of a phase in the presence of another phase can be viewed as single phase flow through a reduced pore network, we can add the subscript <math>i</math> for each phase to Darcy's law above written for Darcy flux, and obtain for each phase in multiphase flow :<math>q_i = -\frac{\kappa_i}{\mu_i} \nabla P_i</math> where <math>\kappa_i</math> is the ''phase permeability'' for phase <math>i</math>. From this we also define [[relative permeability]] <math>\kappa_{ri}</math> for phase <math>i</math> as :<math>\kappa_{ri} = \kappa_{i}/\kappa</math> where <math>\kappa</math> is the permeability for the porous medium, as in Darcy's law. ===Forchheimer equation for non-Darcy flow=== For a sufficiently high flow velocity, the flow is nonlinear, and Dupuit and Forchheimer has proposed to generalize the flow equation to :<math>-\nabla P = \frac{\mu}{\kappa} V + \beta \rho V^2</math> where <math>V</math> is the flow velocity and <math>\beta</math> is a factor to be experimentally deduced. ===In membrane operations=== In pressure-driven [[membrane operation]]s, Darcy's law is often used in the form, :<math>J = \frac{\Delta P - \Delta \Pi} {\mu (R_f + R_m)}</math> where, * <math>J</math> is the volumetric flux (<math>m.s^{-1}</math>), * <math>\Delta P</math> is the hydraulic pressure difference between the feed and permeate sides of the membrane (<math>Pa</math>), * <math>\Delta \Pi</math> is the osmotic pressure difference between the feed and permeate sides of the membrane (<math>Pa</math>), * <math>\mu</math> is the dynamic viscosity (<math>Pa.s</math>), * <math>R_f</math> is the fouling resistance (<math>m^{-1}</math>), and * <math>R_m</math> is the membrane resistance (<math>m^{-1}</math>) == References == <references/> ==See also== *The [[darcy]] unit of fluid [[Permeability (fluid)|permeability]] *[[Hydrogeology]] *[[Groundwater discharge]] *[[Groundwater flow equation]] *[[Groundwater energy balance]] *[[Richards equation]] *[[Darcy friction factor]] == External links == [http://www.calctool.org/CALC/eng/fluid/darcy Browser-based numerical calculator of permeability using Darcy's law.] [[Category:Equations of fluid dynamics]] [[Category:Geology]] [[Category:Hydraulics]] [[Category:Hydraulic engineering]] [[Category:Hydrology]] [[Category:Hydrogeology]] [[Category:Soil mechanics]] [[Category:Soil physics]] [[cs:Darcyho zákon]] [[de:Darcy-Gesetz]] [[eu:Darcyren Legea]] [[fr:Loi de Darcy]] [[ko:다르시의 법칙]] [[it:Legge di Darcy]] [[lt:Darsi dėsnis]] [[hu:Darcy-törvény]] [[nl:Wet van Darcy]] [[pt:Lei de Darcy]] [[ru:Закон Дарси]] [[sl:Darcyjev zakon]] [[sr:Дарсијев закон]] [[sv:Darcys lag]]