User Rating 0.0
Total Usage 0 times
Soil permeability coefficient
Soil Presets:
Area perpendicular to flow direction
h₁ − h₂ (upstream minus downstream)
Distance along flow direction
dimensionless
Fraction of interconnected pore space (0 to 1)
Is this tool helpful?

Your feedback helps us improve.

About

Darcy's Law quantifies fluid flow through porous media under a hydraulic gradient. The relationship Q = K × A × (Δh / L) governs subsurface hydrology, contaminant transport modeling, and aquifer yield estimation. The negative sign indicates flow occurs from high to low hydraulic head. Errors in hydraulic conductivity selection propagate directly into flow predictions - misidentifying clay (K 10−9 m/s) as sand (K 10−4 m/s) yields five orders of magnitude error.

This calculator computes volumetric discharge, specific discharge (Darcy velocity), and actual seepage velocity accounting for effective porosity. Darcy's Law assumes laminar flow (Reynolds number < 1 - 10), saturated conditions, and homogeneous isotropic media. The law breaks down in fractured rock, karst systems, and near pumping wells where turbulent flow dominates.

darcy's law groundwater flow hydraulic conductivity hydrogeology seepage velocity porous media aquifer

Formulas

Darcy's Law describes laminar flow of fluid through saturated porous media under a hydraulic gradient. The volumetric flow rate is proportional to the hydraulic conductivity, cross-sectional area, and hydraulic gradient.

Q = K × A × ΔhL

Where Q = volumetric flow rate (m³/s), K = hydraulic conductivity (m/s), A = cross-sectional area perpendicular to flow (), Δh = hydraulic head difference (m), and L = length of flow path (m). The negative sign indicates flow direction from high to low head.

The hydraulic gradient i is the dimensionless ratio:

i = ΔhL = h1 h2L

Specific discharge (Darcy velocity) represents apparent velocity through the total cross-section:

q = QA = K × i

Actual seepage velocity accounts for the fraction of area available for flow (effective porosity ne):

v = qne = QA × ne

The seepage velocity is always greater than Darcy velocity since ne < 1. For contaminant transport calculations, seepage velocity determines actual travel time through aquifers.

Reference Data

MaterialHydraulic Conductivity K (m/s)Hydraulic Conductivity K (m/day)Effective Porosity neClassification
Gravel (coarse)10−2 - 1864 - 86,4000.25 - 0.35Very High Permeability
Gravel (fine)10−3 - 10−286 - 8640.20 - 0.30High Permeability
Sand (coarse)10−4 - 10−38.6 - 860.25 - 0.35High Permeability
Sand (medium)10−5 - 10−40.86 - 8.60.25 - 0.32Moderate-High
Sand (fine)10−6 - 10−50.086 - 0.860.20 - 0.30Moderate Permeability
Silty Sand10−7 - 10−50.0086 - 0.860.15 - 0.25Low-Moderate
Silt10−8 - 10−60.00086 - 0.0860.01 - 0.20Low Permeability
Sandy Clay10−9 - 10−78.6×10−5 - 0.00860.05 - 0.15Very Low
Clay10−11 - 10−98.6×10−7 - 8.6×10−50.01 - 0.10Aquitard/Aquiclude
Glacial Till10−12 - 10−68.6×10−8 - 0.0860.05 - 0.20Variable
Sandstone10−10 - 10−68.6×10−6 - 0.0860.05 - 0.30Variable (fractured)
Limestone (unfractured)10−9 - 10−68.6×10−5 - 0.0860.01 - 0.20Low-Moderate
Limestone (karst)10−6 - 10−20.086 - 8640.05 - 0.50High (conduit flow)
Fractured Basalt10−8 - 10−40.00086 - 8.60.05 - 0.25Variable
Granite (unfractured)10−14 - 10−108.6×10−10 - 8.6×10−60.0001 - 0.01Impermeable
Shale10−13 - 10−98.6×10−9 - 8.6×10−50.005 - 0.05Aquiclude
Peat10−6 - 10−40.086 - 8.60.30 - 0.50Moderate
Loess10−7 - 10−50.0086 - 0.860.15 - 0.25Low-Moderate
Alluvium10−5 - 10−30.86 - 860.20 - 0.35High
Weathered Rock10−8 - 10−50.00086 - 0.860.10 - 0.30Variable

Frequently Asked Questions

Darcy velocity (specific discharge, q) is the volumetric flow rate divided by total cross-sectional area - it represents an apparent velocity as if water flowed through the entire section. Seepage velocity (v) is the actual average velocity of water particles moving through pore spaces, calculated by dividing Darcy velocity by effective porosity. Since effective porosity ne ranges from 0.01 to 0.35 for most materials, seepage velocity is 3 - 100 times greater than Darcy velocity. For contaminant transport modeling, seepage velocity determines actual travel time.
Darcy's Law assumes laminar flow where viscous forces dominate. It breaks down when Reynolds number exceeds approximately 1 - 10, occurring in: coarse gravel with high gradients, fractured rock with large apertures, karst conduits, near pumping wells at high extraction rates, and unsaturated conditions where two-phase flow occurs. In these cases, the Forchheimer equation adds a quadratic velocity term to account for inertial effects. Darcy's Law also assumes homogeneous isotropic media - anisotropic formations require a hydraulic conductivity tensor.
Hydraulic conductivity can be measured via: (1) laboratory permeameter tests on undisturbed samples - constant head for coarse materials (K > 10−6 m/s), falling head for fine materials; (2) field slug tests or pump tests for in-situ values; (3) empirical correlations using grain-size analysis (Hazen, Kozeny-Carman equations). Laboratory values typically underestimate field conductivity by factors of 2 - 10 due to preferential flow paths and macropores absent in small samples. For preliminary estimates, use the reference table values for your dominant soil type.
The negative sign indicates that flow occurs in the direction of decreasing hydraulic head (opposite to the gradient direction). In practical calculations where you define Δh as the head drop (high minus low) and want positive flow in the downstream direction, the sign is already accounted for. This calculator uses the convention where Δh = hupstream hdownstream (positive value), yielding positive flow rate without the explicit negative sign.
Hydraulic conductivity K is temperature-dependent because it incorporates fluid viscosity: K = kρg/μ, where k is intrinsic permeability (temperature-independent), ρ is fluid density, and μ is dynamic viscosity. Water viscosity at 5°C is approximately 1.5 times that at 25°C, so K decreases by 50% in cold groundwater. Published K values typically reference 20°C; apply correction factors for significantly different temperatures.
Yes, Darcy's Law applies to any flow direction. For vertical flow, include the gravitational component in the total head. The hydraulic head h = z + p/(ρg) combines elevation head and pressure head. For upward flow against gravity, the gradient is reduced; for downward flow, it increases. In vertical percolation through unsaturated zones (infiltration), use Richards' equation which extends Darcy's Law to variable saturation conditions.