Darcy's Law Flow Calculator Back
Porous Media Flow

Darcy's Law Flow Calculator

Select porous medium and fluid, set the pressure gradient and cross-section area, and instantly compute Darcy velocity, volumetric flow rate, hydraulic conductivity and seepage Reynolds number for groundwater work.

Parameters
Porous Medium
Permeability k (m², log slider)
Pressure gradient ΔP/L (Pa/m)
Pa/m
Cross-section area A (m²)
Sample length L (m)
m
Fluid
Results
Results
Darcy velocity q (m/s)
Flow rate Q (m³/s)
Flow rate (L/min)
Hydraulic cond. K (m/s)
Darcy

Porous medium cross-section with REV concept and seepage flow arrows

Theory & Key Formulas
Darcy's law: $q = \dfrac{k}{\mu}\dfrac{\Delta P}{L}$
Volumetric flow: $Q = q \cdot A$
Hydraulic conductivity: $K = \dfrac{k\,\rho g}{\mu}$
Darcy Re: $Re_D = \dfrac{q\sqrt{k}}{\nu}$

Forchheimer equation:
$\dfrac{\Delta P}{L}= \dfrac{\mu}{k}q + \beta\rho q^2$
$\beta \approx 1.75/(\sqrt{k}\,\phi^{1.5})$

What is Darcy's Law?

🙋
What exactly is Darcy's Law? I see it's about flow through sand, but why is it so important in engineering?
🎓
Basically, it's the fundamental equation that describes how fluids like water or oil slowly seep through porous materials like soil, sand, or rock. In practice, it tells us that the flow rate depends on the material's permeability and the pressure pushing the fluid. Try moving the "Permeability k" slider in the simulator above—you'll see how drastically flow changes between gravel and clay.
🙋
Wait, really? So the "Darcy velocity" it calculates isn't the actual fluid speed? That seems confusing.
🎓
Great observation! You're right. The Darcy velocity, often called the "specific discharge" (q), is a volumetric flow rate per area. The actual average fluid speed in the pores is faster, because the fluid only flows through the pore spaces, not the solid grains. For instance, in an aquifer, water might appear to move 1 m/day using Darcy's law, but it's actually zipping through the pores at maybe 3 m/day. The simulator visualizes this with the REV (Representative Elementary Volume) concept.
🙋
So when does Darcy's Law break down? If I crank the pressure gradient way up for water in gravel, will it still be accurate?
🎓
That's the key limitation! Darcy's law assumes slow, laminar flow. At very high flow rates—like in very coarse gravel or near a pumping well—inertial forces become significant and flow becomes non-linear. This is described by the Forchheimer equation. You can explore this boundary by selecting "Gravel" and setting a huge pressure gradient. The flow visualization will hint at this transition, and we check it using the Darcy Reynolds number: $Re_D = q\sqrt{k}/\nu$. If $Re_D \gt 1-10$, Darcy's law may not hold.

Physical Model & Key Equations

The core of the simulator is Darcy's empirical law, which states that the volumetric flow rate through a porous medium is proportional to the cross-sectional area and the pressure gradient, and inversely proportional to the fluid's viscosity.

$$ q = -\frac{k}{\mu}\frac{\Delta P}{L}$$

Here, $q$ is the Darcy velocity or specific discharge [m/s], $k$ is the intrinsic permeability of the medium [m²], $\mu$ is the dynamic viscosity of the fluid [Pa·s], and $\Delta P / L$ is the pressure gradient [Pa/m]. The negative sign indicates flow from high to low pressure.

To get the total volumetric flow rate, we multiply the Darcy velocity by the total cross-sectional area perpendicular to the flow. This is what the simulator calculates as Q.

$$ Q = q \cdot A = -\frac{k A}{\mu}\frac{\Delta P}{L} $$

$Q$ is the total volumetric flow rate [m³/s], and $A$ is the cross-sectional area [m²]. In hydrogeology, a related parameter called Hydraulic Conductivity ( $K$ ) is often used for water: $K = k \rho g / \mu$, which combines medium and fluid properties.

Frequently Asked Questions

Darcy velocity is the apparent velocity obtained by dividing the flow rate by the total cross-sectional area of the porous medium. The actual velocity through the interstitial spaces (interstitial velocity) is the Darcy velocity divided by the porosity. The simulator displays the Darcy velocity.
It allows for a quick comparison of flow rate differences due to variations in permeability and porosity among different porous media (e.g., sand, clay, sandstone). This is useful for selecting the optimal medium during design or for intuitively understanding the impact of changes in pressure gradient.
It is an indicator for determining whether the flow within the porous medium is laminar or turbulent. Generally, inertial effects become non-negligible when Re > 10. Since the simulator calculates it automatically, please use it to verify the applicable range.
Selecting a preset such as water or oil from the fluid selection dropdown will automatically set the dynamic viscosity and density. To input a custom fluid, you can directly edit each value to reflect the desired physical properties.

Real-World Applications

Groundwater Resource Management: Hydrogeologists use Darcy's Law to model aquifer recharge, map contaminant plume migration, and design sustainable pumping wells. For instance, calculating how fast a chemical spill might reach a drinking water well is a direct application of this tool's parameters.

Oil and Gas Reservoir Engineering: Predicting the flow of crude oil or natural gas through underground rock formations is fundamental to estimating well production rates and planning enhanced recovery methods. The choice of fluid (oil vs. water) in the simulator shows how viscosity drastically impacts recoverable reserves.

Geotechnical and Civil Engineering: Designing effective drainage systems for earth dams, road embankments, or building foundations requires accurate seepage analysis to prevent structural failure. The pressure gradient and material permeability are key design inputs.

Filter and Membrane Design: Industrial filters for water treatment or air purification operate on Darcy flow principles. Engineers select materials with specific permeabilities (like the "Sand" or "Silt" options) to achieve the desired flow rate and particle removal efficiency.

Common Misconceptions and Points to Note

When you start using this tool, there are a few points where newcomers, especially those tasked with calculations from the field, often get tripped up. First, "the Darcy flux q is not the actual average fluid velocity." Darcy flux is the apparent velocity obtained by dividing by the total cross-sectional area of the porous medium. Actual water flows through the complex pathways between sand grains, so its real velocity is much faster. For example, in sandstone with a porosity of 0.3 (30%), the actual average interstitial velocity would be the Darcy flux divided by 0.3.

Next, be careful not to mix up the units and order of magnitude of permeability k. While m² is the standard in textbooks, "Darcy (D)" or "milliDarcy (mD)" are commonly used in the field. One Darcy is an extremely small value, approximately 9.87×10⁻¹³ m². Good sandstone might be 100-1000 mD, while tight shale can be below 0.001 mD. Don't be alarmed if you see exponential notation like "1e-12" on the slider.

Finally, watch out for the pitfall that "the Forchheimer correction coefficient β is not independent of k." While you can set them independently in the tool, in reality, k and β are related based on the porous structure (grain size and shape). In the absence of experimental data, you might estimate it using an empirical rule like "β ≈ 1.8 / (k^0.5)". Since including the correction can dramatically change the results, always verify the source and justification for any coefficient value you use.

How to Use

  1. Select porous medium type (sand, gravel, or clay) using the material dropdown to set initial hydraulic conductivity K
  2. Set pressure gradient (ΔP/L) in kPa/m using the slider or numeric input—typical range 0.1 to 100 kPa/m for groundwater flow
  3. Specify cross-sectional area in m² and flow path length in meters
  4. The calculator computes Darcy velocity q (m/s), volumetric flow rate Q (m³/s and L/min), and Reynolds number to verify laminar flow validity

Worked Example

Medium-grained sand with K=0.001 m/s, pressure gradient ΔP/L=5 kPa/m, area A=0.5 m², flow length L=10 m. Darcy velocity q = K × (ΔP/L)/ρg = 0.001 × (5000/10)/(9.81×1000) ≈ 5.1×10⁻⁵ m/s. Flow rate Q = q × A = 5.1×10⁻⁵ × 0.5 = 2.55×10⁻⁵ m³/s or 1.53 L/min. Reynolds number Re ≈ 0.008, confirming laminar seepage regime.

Practical Notes

  1. Keep Reynolds number below 1–10 to ensure Darcy's law validity; higher Re indicates turbulent flow requiring non-linear models
  2. Hydraulic conductivity varies with grain size: fine sand K≈10⁻⁵ m/s, gravel K≈10⁻² m/s, clay K≈10⁻⁸ m/s
  3. For contaminant transport modeling, use computed seepage velocity to estimate pollutant arrival time at monitoring wells