Darcys Law Flow Rate Calculator - Geologic
Calculate groundwater flow rate, Darcy velocity, seepage velocity, and travel time using Darcy's Law for geologic and hydrogeologic applications.
Darcys Law Flow Rate Calculator - Geologic
Calculator
Adjust values & calculateEnter your values below. Every result is computed in your browser โ no data is sent to any server.
Formula: Q = K * A * (dH / dL); v_s = v_d / n
Worked example โ i: 0.015 | Q: 0.00375 m3/s (324 m3/day) | Darcy velocity: 6.5 m/day | Seepage velocity: 25.9 m/day
Formula
Q = K * A * (dH / dL); v_s = v_d / n
This Darcy computes results from your provided inputs using the calculator's underlying model.
Worked Examples
Example 1: Groundwater Flow Through a Sand Aquifer
Problem:A sand aquifer has hydraulic conductivity K = 0.005 m/s, a flow cross-section of 50 m2 and an effective porosity of 0.25. Two wells 200 m apart show a head difference of 3 m. Find the hydraulic gradient, the discharge, and the Darcy and seepage velocities.
Solution:Hydraulic gradient i = dh / L = 3 / 200 = 0.015 Discharge Q = K x A x i = 0.005 x 50 x 0.015 = 0.00375 m3/s Per day: 0.00375 x 86,400 = 324 m3/day Darcy (specific discharge) velocity q = K x i = 0.005 x 0.015 = 7.5e-5 m/s = 6.48 m/day Seepage (pore water) velocity v = q / n = 7.5e-5 / 0.25 = 3.0e-4 m/s = 25.9 m/day A tracer released at one well therefore reaches the other in about 200 / 25.9 = 7.7 days, four times faster than the Darcy velocity alone would suggest.
Result:i: 0.015 | Q: 0.00375 m3/s (324 m3/day) | Darcy velocity: 6.5 m/day | Seepage velocity: 25.9 m/day
Example 2: Seepage Through a Compacted Clay Dam Core
Problem:An earth dam core has K = 1e-7 m/s, a flow cross-section of 200 m2 and porosity 0.35. The reservoir stands 15 m above the downstream face across a 30 m core thickness. How much water is lost through the core?
Solution:Hydraulic gradient i = 15 / 30 = 0.5 Q = K x A x i = 1e-7 x 200 x 0.5 = 1.0e-5 m3/s Per day: 1.0e-5 x 86,400 = 0.864 m3/day, i.e. under one cubic metre a day Seepage velocity v = (K x i) / n = (1e-7 x 0.5) / 0.35 = 1.43e-7 m/s = 0.0123 m/day, so a water particle needs about 30 / 0.0123 = 2,400 days, roughly 6.6 years, to cross the core Even though the gradient here is 33 times steeper than in the sand aquifer, the conductivity is 50,000 times lower, and the discharge falls by a factor of 375.
Result:i: 0.5 | Q: 1.0e-5 m3/s (0.864 m3/day) | Travel time through core: about 6.6 years
Frequently Asked Questions
What is Darcy's Law and how does it apply to geology?
Darcy's Law describes the flow of fluid through a porous medium and is fundamental to groundwater hydrology and geologic engineering. The law states that the flow rate Q equals the hydraulic conductivity K times the cross-sectional area A times the hydraulic gradient (dH/dL). In geologic applications, it governs groundwater movement through aquifers, seepage through earth dams, contaminant transport in subsurface environments, and dewatering calculations for excavations.
What is the difference between Darcy velocity and seepage velocity?
Darcy velocity (also called specific discharge or flux) is the apparent flow velocity calculated as Q/A, assuming flow occurs through the entire cross-sectional area. Seepage velocity (also called pore velocity or actual velocity) is the true average velocity of water moving through the pore spaces, calculated by dividing the Darcy velocity by the effective porosity. Seepage velocity is always greater than Darcy velocity because water can only flow through the interconnected void spaces, not through solid particles.
When does Darcy's Law not apply?
Darcy's Law assumes laminar flow through a porous medium and breaks down under several conditions. At very high flow velocities in coarse gravel or rock fractures, turbulent flow occurs and the relationship between gradient and velocity becomes nonlinear. The Reynolds number threshold is typically around 1-10 for porous media. Darcy's Law also does not apply in unsaturated conditions without modification, in soils exhibiting non-Newtonian fluid behavior, or in very low permeability clays where electrochemical forces influence flow.
How is the hydraulic gradient calculated, and why is it dimensionless?
The hydraulic gradient i is the head loss divided by the flow path length, i = dH / dL. Both terms are lengths, so the ratio has no units. Hydraulic head itself combines elevation head and pressure head, so it must be read from piezometers screened in the same aquifer unit, not from ground surface elevations. Natural gradients are typically small, often between 0.0001 and 0.05, which is why groundwater moves in metres per year rather than metres per second.
What is the difference between hydraulic conductivity and intrinsic permeability?
Hydraulic conductivity K describes how easily a specific fluid moves through a specific medium, so it depends on both the pore geometry and on fluid density and viscosity. Intrinsic permeability k depends only on the rock or soil framework and has units of area, commonly reported in darcys or square metres. The two are linked by K = k * rho * g / mu. Petroleum and geothermal work prefers k because temperature and fluid changes alter viscosity strongly; groundwater work usually assumes cool fresh water and quotes K directly.
Why is effective porosity used for seepage velocity rather than total porosity?
Total porosity counts every void, including dead-end pores and the water held on grain surfaces by capillary and adsorptive forces. Only the interconnected, drainable fraction actually transmits flow, and that is the effective porosity. In clean sands the two values are close, roughly 0.20 to 0.35, but in clay a total porosity above 0.50 can pair with an effective porosity of only a few percent. Using total porosity would understate seepage velocity and make contaminant arrival estimates dangerously late.
How do you use Darcys Law Flow Rate Calculator - Geologic to estimate contaminant travel time?
Travel time follows from the seepage velocity: t = L / v_s, where v_s = K * i / n_e. This is the advective arrival time for a conservative, non-reactive tracer such as chloride. Dissolved organic compounds usually arrive later because sorption to organic carbon retards them by a factor R, so t_retarded = R * t. Hydrodynamic dispersion also spreads the front, so first detection typically precedes the calculated centre-of-mass arrival.
Why is horizontal hydraulic conductivity usually greater than vertical conductivity?
Sedimentary deposits are laid down in layers, so a stratified sequence behaves as an anisotropic unit. Horizontal flow is dominated by the most permeable bed, which acts like a parallel conduit, while vertical flow must cross every low-permeability layer in series and is throttled by the tightest one. The resulting anisotropy ratio Kh/Kv is commonly between 5 and 100 and can exceed 1000 in varved or clay-interbedded sequences. Using a single isotropic K therefore overestimates leakage across confining beds.
How does transmissivity relate to the hydraulic conductivity entered here?
Transmissivity T is the conductivity integrated over the saturated thickness b of the aquifer, T = K * b, with units of length squared per time such as m2/day. It measures how much water the whole aquifer thickness can transmit under a unit gradient, which is what pumping tests actually resolve. Hydraulic conductivity is the property of the material itself and is what Darcy's Law needs alongside a cross-sectional area. In an unconfined aquifer b is the water table thickness, so transmissivity falls as the water table is drawn down.
Which test should I use to measure hydraulic conductivity for a given site?
Constant-head permeameters suit clean sands and gravels in the laboratory, while falling-head permeameters handle silts and clays where flow is too slow to measure at constant head. In the field, a slug test gives a quick value for the small volume immediately around a single well, whereas a multi-well pumping test stresses a much larger volume and returns a bulk transmissivity that better represents formation-scale behaviour. Lab values from disturbed samples routinely underestimate field conductivity because they miss fractures and coarse lenses.
References
Background & Theory
History
Reviewed for accuracy by Daniel Agrici, Founder & Lead Developer ยท Editorial policy
Related Calculators
๐งฎGlacial Flow Velocity (glenโs Law) Calculator
Calculate glacial flow velocity (glenโs law) with inputs, formulas, and instant results.
๐งฎFault Slip Rate Calculator
Calculate fault slip rate with inputs, formulas, and instant results.
๐งฎGeothermal Heat Flow Calculator
Calculate geothermal heat flow with inputs, formulas, and instant results.
๐งฎHeat Flow from Gradient Calculator
Calculate heat flow from gradient and conductivity with inputs, formulas, and instant results.
๐งฎFlow Duration Curve Calculator
Calculate flow duration curve with inputs, formulas, and instant results.
๐งฎFlow Resistance Converter (nโf) Calculator
Calculate flow resistance converter (nโf) with inputs, formulas, and instant results.
๐งฎGroundwater Flow Velocity Calculator
Calculate groundwater flow velocity with inputs, formulas, and instant results.
๐งฎHydraulic Gradient Calculator (natural Flow)
Calculate hydraulic gradient (natural flow) with inputs, formulas, and instant results.