Darcy's Law Calculator for Porous Media Flow
Calculate fluid flow rate, Darcy velocity, and seepage velocity through porous materials using Darcy's Law — essential for groundwater hydrology and reservoir engineering.
Enter permeability, cross-sectional area, pressure difference, viscosity, flow length, and porosity to calculate volumetric flow rate and velocities.
Darcy's Law Calculator for Porous Media Flow
Calculate fluid flow rate, Darcy velocity, and seepage velocity through porous materials using Darcy's Law — essential for groundwater hydrology and reservoir engineering.
About Darcy's Law Calculator
Darcy's Law is one of the most fundamental equations in fluid dynamics, describing how fluids flow through porous media. First formulated by Henry Darcy in 1856 after conducting experiments on water filtration through sand beds in Dijon, France, the law establishes a linear relationship between the volumetric flow rate and the applied pressure gradient. It is expressed as Q = kA·ΔP / (μL), where Q is the volumetric flow rate (m³/s), k is the intrinsic permeability of the medium (m²), A is the cross-sectional area perpendicular to flow (m²), ΔP is the pressure difference driving the flow (Pa), μ is the dynamic viscosity of the fluid (Pa·s), and L is the length of the flow path (m).
Permeability is the single most important parameter in Darcy's Law. It is a property of the porous medium alone — independent of the fluid — and quantifies the medium's ability to transmit fluid based on its pore structure, pore connectivity, and tortuosity. Permeability spans many orders of magnitude: clay sits at 10⁻²⁰ to 10⁻¹⁸ m², fine sand at 10⁻¹⁶ to 10⁻¹⁴ m², coarse sand and gravel at 10⁻¹⁴ to 10⁻¹⁰ m², and highly fractured rock or porous concrete at 10⁻¹⁰ m² and above. In petroleum engineering, permeability is often expressed in millidarcies (1 mD = 9.869×10⁻¹⁶ m²).
Two velocities emerge from Darcy's Law. The Darcy velocity (or superficial velocity) is v = Q/A, representing the apparent velocity as if the fluid occupied the entire cross-section including the solid matrix. The seepage velocity (or pore velocity) is the actual average speed of fluid through the pore spaces: v_seepage = v/φ, where φ is the porosity. Since only the pores conduct fluid, the seepage velocity is always higher than the Darcy velocity by a factor of 1/φ. For a medium with 25% porosity, the fluid moves through the pores four times faster than the Darcy velocity suggests.
Darcy's Law underpins hydrogeology (modelling groundwater aquifer flows and contaminant transport), petroleum engineering (reservoir simulation and production forecasting), soil science (irrigation and drainage design), chemical engineering (packed-bed reactors and membrane filtration), and civil engineering (dam seepage analysis and foundation drainage). It is valid under the assumption of laminar, steady, incompressible flow through a homogeneous, isotropic porous medium saturated with a Newtonian fluid. At high flow rates where inertial effects become significant, the Forchheimer equation adds a quadratic velocity term; at very small scales, slip flow (Knudsen diffusion) may require the Klinkenberg correction.
The calculator uses the magnitude of the pressure difference, so enter the absolute pressure drop across the sample regardless of sign convention. Results give the magnitude of flow rate and velocities in the direction of the pressure gradient.
Worked Examples
Four representative porous media flow scenarios demonstrating Darcy's Law across different engineering applications.
| Scenario | Result | Notes |
|---|---|---|
| Sandstone reservoir: k=1×10⁻¹² m², A=0.01 m², ΔP=10⁶ Pa, μ=0.001 Pa·s, L=0.1 m, φ=0.25 | Q = 1×10⁻⁴ m³/s; v_darcy = 1×10⁻² m/s; v_seepage = 4×10⁻² m/s | Typical petroleum reservoir flow. High pressure differential drives significant flow through this rock sample. |
| Sandy soil: k=1×10⁻¹⁰ m², A=0.1 m², ΔP=1000 Pa, μ=0.001 Pa·s, L=1.0 m, φ=0.35 | Q = 1×10⁻⁵ m³/s; v_darcy = 1×10⁻⁴ m/s; v_seepage ≈ 2.86×10⁻⁴ m/s | Groundwater flow through sandy aquifer. Low pressure gradient yields slow but steady seepage. |
| Industrial ceramic filter: k=1×10⁻¹⁴ m², A=0.001 m², ΔP=50,000 Pa, μ=0.001 Pa·s, L=0.05 m, φ=0.15 | Q = 1×10⁻⁸ m³/s; v_darcy = 1×10⁻⁵ m/s; v_seepage ≈ 6.67×10⁻⁵ m/s | Very tight filter medium requires high pressure to achieve measurable flow rate. |
How to Use the Darcy's Law Calculator
- Enter the intrinsic permeability k in m². Consult published tables for your medium type or use laboratory-measured values. Convert from millidarcies using 1 mD = 9.869×10⁻¹⁶ m².
- Enter the cross-sectional area A in m² perpendicular to flow. For a circular core sample, A = π·r²; for a rectangular slab, A = width × height.
- Enter the pressure difference ΔP in Pascals — the pressure drop from inlet to outlet driving the flow. Enter as a positive value.
- Enter the dynamic viscosity μ in Pa·s for your fluid at the operating temperature. Water at 20°C is 0.001 Pa·s; viscosity increases for oils and decreases with temperature.
- Enter the flow length L in metres and porosity φ as a decimal fraction between 0 and 1 (e.g. 0.30 for 30% porosity). Click Calculate to see flow rate, Darcy velocity, and seepage velocity.
Frequently Asked Questions
What is permeability and how do I find it?
Permeability (k) is an intrinsic property of the porous medium describing how easily fluid can flow through it. It depends only on pore structure, not on the fluid. You can find it from laboratory permeameter tests on core samples, published data tables for your material type, or by inverting Darcy's Law if you can measure flow rate and pressure drop in a known geometry.
What is the difference between Darcy velocity and seepage velocity?
Darcy velocity (q = Q/A) is the volumetric flow rate per unit total cross-sectional area — it treats the porous medium as if it were a solid tube. Seepage velocity is the actual average speed of the fluid through the connected pores: v_seepage = q/φ. It is always higher than the Darcy velocity because only the void fraction (porosity) carries the flow.
When is Darcy's Law not valid?
Darcy's Law breaks down at high flow rates (high Reynolds number in pores) where inertial forces become significant — typically when Re > 1–10 based on pore size. The Forchheimer equation adds a quadratic drag term for these conditions. It also fails for gas flow at very low pressures (Klinkenberg slip flow) and in highly heterogeneous or fractured media where flow channels bypass most of the matrix.
How do I convert permeability from millidarcies to m²?
1 darcy = 9.869233×10⁻¹³ m², so 1 millidarcy (mD) = 9.869233×10⁻¹⁶ m². Multiply your permeability in mD by 9.869×10⁻¹⁶ to get m². Many petroleum reservoirs have permeabilities of 1–1000 mD, corresponding to 10⁻¹⁵ to 10⁻¹² m².
How does temperature affect the calculation?
Temperature primarily affects fluid viscosity. For water, viscosity drops from 0.00179 Pa·s at 0°C to 0.000283 Pa·s at 100°C — a sixfold reduction. Higher temperature means lower viscosity and therefore higher flow rate for the same pressure gradient. Always use the viscosity at the actual operating temperature for accurate results.
What is hydraulic conductivity and how does it relate to permeability?
Hydraulic conductivity K (m/s) combines permeability with fluid properties: K = k·ρ·g/μ, where ρ is fluid density and g is gravitational acceleration. It is commonly used in groundwater hydrology where the fluid is water at a known temperature. Permeability k (m²) is the purely physical property of the medium; hydraulic conductivity K already folds in the fluid. This calculator uses permeability for generality across all fluids.