← GW Home
Chapter 2

Fundamentals of Groundwater Motion [8 Hours]

GW Motion & Darcy's Law ~12+ Past Questions

MOST IMPORTANT TOPICS & MARKS DISTRIBUTION

High-Yield Topics (by exam frequency):

  1. Darcy’s Law — derivation, validity, and limitations — Asked in 2078, 2073, 2072 Ashwin, 2072 Magh, 2071 [4–8 marks]
  2. REV and continuum approach — Asked in 2078, 2073, 2071 [4–6 marks]
  3. Hydraulic conductivity, permeability, transmissivity — definitions and relations — Asked in 2073, 2072 Ashwin [2–4 marks]
  4. Numericals: Darcy validity domain, seepage velocity — Asked in 2078, 2073, 2071 [4–8 marks]

Marks Each Time Asked:

ExamQuestionsTotal
2078 ChaitraQ2a [6], Q2b [6 num]12
2073 BhadraQ2a [4+2], Q2b [8 num]14
2072 AshwinQ1 [4], Q2 [8]12
2072 MaghQ2 [6]6
2071 BhadraQ2 [4], Q6 [4+4 num]12

CHAPTER ROADMAP

  1. Continuum approach and Representative Elementary Volume (REV)
  2. Darcy's experiment and one-dimensional groundwater flow
  3. Derivation of the three-dimensional groundwater flow equation
  4. Validity of Darcy's law and non-Darcian flow
  5. Hydraulic conductivity, permeability, and transmissivity
  6. Aquifer heterogeneity and anisotropy
  7. Past question bank for Chapter 2
  8. Worked numerical solutions
  9. Chapter summary and exam formula sheet

2.1 CONTINUUM APPROACH AND REPRESENTATIVE ELEMENTARY VOLUME

PAST EXAM QUESTIONS ADDRESSED IN THIS SECTION:

  1. Define Representative Elementary Volume and continuum approach of analysis. [Part of Q2a, 6 marks] (2078 Chaitra, Q2a)
  2. Explain the continuum and Representative Elementary Volume approach in groundwater flow analysis. [Part of Q2a, 4 marks] (2073 Bhadra, Q2a)
  3. What do you understand by Representative Elementary Volume? Explain. [4] (2071 Bhadra, Q2)
CORE IDEA: Groundwater actually moves through an intricate network of pores and fractures, but engineering analysis becomes possible only when the porous medium is treated as a continuum whose averaged properties can represent the microscopic void system.

At pore scale, water flow is highly irregular. Every grain, pore throat, and constriction changes the local direction and magnitude of velocity. If groundwater were analyzed pore by pore, the governing geometry would become too complex for routine engineering use. The continuum approach solves this difficulty by replacing the actual porous mass with an equivalent continuous material having averaged properties such as porosity, hydraulic conductivity, and storage.

The averaging is performed over a volume that is large compared with individual pores but still small compared with the total aquifer dimensions. This critical averaging scale is called the Representative Elementary Volume (REV). Inside the REV, the measured property stops fluctuating erratically and begins to stabilize around a meaningful average value.

DEFINITION — Continuum Approach:
The continuum approach is the method of treating a porous medium — which is actually a complex mixture of solid grains and void spaces — as if it were a single, continuous material with smoothly varying, spatially averaged properties (porosity, hydraulic conductivity, specific storage, etc.). Instead of tracking water flow through each individual pore and around each grain, the continuum approach replaces the discrete grain-pore system with an equivalent homogeneous medium at every point, so that the governing physics can be expressed as partial differential equations (Darcy's law, the groundwater flow equation). This replacement is valid only when the averaging is performed over a volume at least as large as the Representative Elementary Volume (REV).
DEFINITION — Representative Elementary Volume (REV):
The REV is the smallest volume of a porous medium over which a measurable property (such as porosity, hydraulic conductivity, or specific storage) yields a stable, reproducible average value that is independent of the sample size. Below this volume, the property fluctuates erratically due to pore-scale heterogeneity; above it, the value remains essentially constant (until formation-scale heterogeneity introduces new variability). The REV is the theoretical foundation that permits the replacement of the actual grain-and-pore system with a continuum amenable to partial differential equations (Darcy's law, flow equation).
1. Pore Scale Every grain changes flow path Too complex for engineering Average over REV 2. REV Scale Volume REV Model 3. Continuum h(x,y,z,t) q = −K∇h Smooth PDEs govern flow Pore complexity → REV averaging → Continuum PDEs (Darcy’s law, flow equation)
Figure 2.1: Groundwater modeling moves from pore-scale complexity to REV-based averaging and finally to a continuum description.
REV Threshold — Property Stabilization with Sample Size Sample volume → Measured property (n or K) Pore-scale fluctuations Stable — meaningful average Vᵣᵒᵥ n̄ Single pore or grain dominates Average independent of sample size Below REV: erratic | At/above REV: stable, reproducible average
Figure 2.2: The REV is reached when the averaged property becomes practically stable with increasing sample size.

Why the Continuum Approach Is Important

  1. It permits the use of differential equations for groundwater flow.
  2. It allows the definition of field-scale properties such as hydraulic conductivity and storativity.
  3. It links laboratory observations and pumping-test results to aquifer-scale predictions.
  4. It justifies the use of Darcy's law as a macroscopic law even though pore-scale velocities vary strongly.

Limitations of the REV/Continuum Assumption

  • It becomes weak in highly fractured rock where the flow is controlled by a few discrete fractures.
  • It may fail near cavities, karst conduits, or very coarse gravel where pore-scale effects are dominant.
  • It is scale-dependent: a volume acceptable as REV in one formation may be too small or too large in another.
  • It does not eliminate heterogeneity; it only smooths the material enough to define local averages.
Where the REV / Continuum Approach Fails Fractured Rock Flow only in discrete fractures Karst / Solution Channels Conduit/pipe flow, not porous Sample < REV Too small! Property fluctuates wildly Use discrete-fracture or dual-porosity models instead
Figure 2.2b: The REV/continuum approach fails in fractured rock, karst systems, and when the sample volume is smaller than the REV.

2.2 DARCY'S EXPERIMENT AND ONE-DIMENSIONAL GROUNDWATER FLOW

PAST EXAM QUESTIONS ADDRESSED IN THIS SECTION:

  1. What is Darcy's law for groundwater flow? Derive the one-dimensional form and explain its validity. [6] (2072 Magh, Q2)

In 1856, Henry Darcy performed controlled sand-column experiments and observed that the discharge through a porous medium is directly proportional to the cross-sectional area and to the head loss, and inversely proportional to the flow length. This gave the empirical law that forms the basis of groundwater hydraulics.

Darcy’s Sand-Column Experiment (1856) Saturated Sand Inflow Outflow Q h₁ h₂ Δh = h₁ − h₂ Q (flow direction) L Area = A Q = K · A · (h₁ − h₂) / L = K · A · i Q ∝ A, Q ∝ Δh, Q ∝ 1/L — Darcy’s empirical observation
Figure 2.3: Darcy’s original experiment — discharge through saturated sand is proportional to area and head loss, inversely proportional to length.
DARCY'S LAW: The discharge rate through a saturated porous medium is proportional to the hydraulic gradient and the cross-sectional area normal to flow.
$$ Q = K A \frac{h_1 - h_2}{L} = K A i $$

In sign-conscious form, groundwater moves from high head to low head. Therefore, if the positive coordinate axis is taken in the flow direction, Darcy's law is written as:

$$ q = \frac{Q}{A} = -K \frac{dh}{dl} $$

where q is the Darcy flux or specific discharge, K is hydraulic conductivity, and dh/dl is the hydraulic gradient. The negative sign indicates that head decreases in the direction of flow.

DARCY FLUX

Darcy flux is the discharge per total cross-sectional area of the porous medium. It does not represent the actual average velocity inside pores because part of the area is occupied by solids.

$$ q = \frac{Q}{A} $$

SEEPAGE VELOCITY

Actual average velocity through pore spaces is greater than Darcy flux and is found by dividing by effective porosity.

$$ v_s = \frac{q}{n_e} = \frac{K i}{n_e} $$
Darcy Flux (q) vs Seepage Velocity (vₛ) Darcy Flux: q = Q/A Total area A (solid + voids) Seepage Velocity: vₛ = q/nₑ Pore area only: Aᵥ = nₑA vₛ = q / nₑ = Ki / nₑ > q (since nₑ < 1)
Figure 2.3b: Darcy flux is computed over total cross-section; seepage velocity through pores is q/nₑ, always larger since nₑ < 1.
Symbol Meaning Typical unit
Q Total discharge m3/s or m3/day
A Area normal to flow m2
i Hydraulic gradient, \((h_1-h_2)/L\) Dimensionless
K Hydraulic conductivity m/s or m/day
q Darcy flux or specific discharge m/s
Darcy’s Law — Sign Convention h₁ (high) h₂ (low) dh/dx < 0 (head decreases) Flow direction → x L q = −K(dh/dx) = −K(h₂−h₁)/L = K(h₁−h₂)/L > 0 Negative sign makes q positive in the direction of decreasing head
Figure 2.3c: The negative sign in Darcy’s law ensures flow is in the direction of decreasing hydraulic head.

Derivation of the One-Dimensional Groundwater Flow Equation

Starting from Darcy's law, the one-dimensional groundwater flow equation is derived by applying the principle of conservation of mass (continuity) to a control volume within a saturated confined aquifer.

1D Control Volume for Flow Equation Derivation Saturated Confined Aquifer Control Volume qₓ qₓ + (∂qₓ/∂x)dx dx x Net inflow = −(∂qₓ/∂x)·dx = Storage change = Sₛ·(∂h/∂t)·dx
Figure 2.3d: One-dimensional control volume — inflow at x, outflow at x + dx; net inflow equals rate of storage change.

Step 1 — Darcy's law in 1D:

$$ q_x = -K \frac{\partial h}{\partial x} $$

Step 2 — Mass balance on the control volume:

Consider a small element of length \(dx\) and cross-sectional area \(A\). Using a first-order Taylor expansion, the flux leaving the element at \(x + dx\) is:

$$ q_{x+dx} = q_x + \frac{\partial q_x}{\partial x} \, dx $$

The net inflow per unit area into the element is:

$$ q_x - \left( q_x + \frac{\partial q_x}{\partial x} \, dx \right) = -\frac{\partial q_x}{\partial x} \, dx $$

Step 3 — Storage term:

In a confined aquifer, any change in hydraulic head causes a change in the volume of water stored, characterised by the specific storage \(S_s\). The rate of change of storage per unit volume is:

$$ \text{Storage change per unit volume} = S_s \frac{\partial h}{\partial t} $$

Step 4 — Combine mass balance and storage:

Equating net inflow to the rate of storage change for the element of volume \(A \cdot dx\):

$$ -\frac{\partial q_x}{\partial x} = S_s \frac{\partial h}{\partial t} $$

Step 5 — Substitute Darcy's law:

Replace \(q_x = -K \frac{\partial h}{\partial x}\) in the continuity equation:

$$ \frac{\partial}{\partial x}\left( K \frac{\partial h}{\partial x} \right) = S_s \frac{\partial h}{\partial t} $$

For a homogeneous aquifer where \(K\) is constant in space:

$$ K \frac{\partial^2 h}{\partial x^2} = S_s \frac{\partial h}{\partial t} \qquad \text{(1D Groundwater Flow Equation)} $$

For steady-state flow (\(\partial h/\partial t = 0\)), this reduces to:

$$ \frac{d^2 h}{d x^2} = 0 \qquad \text{(Laplace equation in 1D)} $$

This implies that the hydraulic head varies linearly with distance under steady one-dimensional flow — consistent with Darcy's original observation.

Validity of Darcy's Law

Darcy's law holds under the following conditions:

  1. Laminar flow: The pore-scale Reynolds number \(Re = q d_m / \nu\) must remain below approximately 1 to 10. Most natural groundwater systems satisfy this criterion.
  2. Viscous dominance: Fluid inertia must be negligible compared with viscous resistance. This is inherently tied to the Reynolds number condition.
  3. Saturated porous medium: The medium must be fully saturated (for unsaturated flow, a modified form using relative permeability is needed).
  4. Newtonian fluid: The flowing fluid must obey Newton's law of viscosity.

Darcy's law fails in the following situations:

  • Turbulent flow near pumping wells in coarse gravel packs.
  • Flow through karst conduits and large open fractures where channel-type flow dominates.
  • Extremely slow flow through dense clays where electrochemical and threshold-gradient effects become significant.
  • High-velocity flow through well screens and gravel envelopes.

2.3 DERIVATION OF THE THREE-DIMENSIONAL GROUNDWATER FLOW EQUATION

PAST EXAM QUESTIONS ADDRESSED IN THIS SECTION:

  1. Starting from the general expression of Darcy's law for groundwater flow, derive the three-dimensional form of groundwater flow. [6] (2072 Ashwin, Q2)

Natural groundwater flow is three-dimensional. Darcy's law applied independently along the three principal axes gives the flux components:

$$ q_x = -K_x \frac{\partial h}{\partial x}, \qquad q_y = -K_y \frac{\partial h}{\partial y}, \qquad q_z = -K_z \frac{\partial h}{\partial z} $$

In isotropic media, \(K_x = K_y = K_z = K\), and the vector form becomes:

$$ \vec{q} = -K \nabla h $$

Derivation Using Conservation of Mass

To derive the general 3D groundwater flow equation, consider a differential control volume \(dx \cdot dy \cdot dz\) within a saturated aquifer. We apply the principle of conservation of mass: the net rate of fluid inflow must equal the rate of change of fluid storage within the element.

3D Control Volume — Flux Balance for Flow Equation Derivation dx dy dz qₓ qᵧ qᵪ qₓ + (∂qₓ/∂x)dx at face x + dx qᵧ + (∂qᵧ/∂y)dy at face y + dy qᵪ + (∂qᵪ/∂z)dz at face z + dz Sₛ · (∂h/∂t) × dx·dy·dz x z y Inflow Outflow ∑ (Inflow − Outflow) in x, y, z = Sₛ · (∂h/∂t) · dx·dy·dz
Figure 2.4: Differential control volume (dx × dy × dz) in a saturated aquifer. Green arrows show inflow fluxes entering three faces; red arrows show outflow fluxes leaving opposite faces, each modified by a first-order Taylor expansion. The storage change inside the element equals the net inflow.

Step 1 — Net flux in the x-direction:

The flux entering through the face at \(x\) is \(q_x \cdot dy \cdot dz\). The flux leaving through the face at \(x + dx\) is \(\left(q_x + \frac{\partial q_x}{\partial x} dx\right) \cdot dy \cdot dz\). The net inflow in the x-direction is:

$$ \text{Net}_x = -\frac{\partial q_x}{\partial x} \, dx \, dy \, dz $$

Step 2 — Net flux in y and z directions:

Applying the same Taylor-expansion argument to the y and z faces:

$$ \text{Net}_y = -\frac{\partial q_y}{\partial y} \, dx \, dy \, dz, \qquad \text{Net}_z = -\frac{\partial q_z}{\partial z} \, dx \, dy \, dz $$

Step 3 — Total net inflow and storage:

The total net inflow into the element is the sum of contributions from all three directions. This must equal the rate of change of fluid storage within the element, expressed using the specific storage \(S_s\):

$$ -\left( \frac{\partial q_x}{\partial x} + \frac{\partial q_y}{\partial y} + \frac{\partial q_z}{\partial z} \right) dx \, dy \, dz = S_s \frac{\partial h}{\partial t} \, dx \, dy \, dz $$

Dividing both sides by \(dx \, dy \, dz\):

$$ -\frac{\partial q_x}{\partial x} - \frac{\partial q_y}{\partial y} - \frac{\partial q_z}{\partial z} = S_s \frac{\partial h}{\partial t} $$

Step 4 — Substitute 3D Darcy's law:

Replacing each flux component with its Darcy expression (\(q_x = -K_x \partial h / \partial x\), etc.):

$$ \frac{\partial}{\partial x}\left(K_x \frac{\partial h}{\partial x}\right) + \frac{\partial}{\partial y}\left(K_y \frac{\partial h}{\partial y}\right) + \frac{\partial}{\partial z}\left(K_z \frac{\partial h}{\partial z}\right) = S_s \frac{\partial h}{\partial t} $$

This is the general three-dimensional groundwater flow equation for a heterogeneous, anisotropic confined aquifer.

Important Special Cases

Case 1: Homogeneous, isotropic aquifer (\(K_x = K_y = K_z = K = \text{constant}\)):

$$ K \left( \frac{\partial^2 h}{\partial x^2} + \frac{\partial^2 h}{\partial y^2} + \frac{\partial^2 h}{\partial z^2} \right) = S_s \frac{\partial h}{\partial t} $$

or equivalently:

$$ \nabla^2 h = \frac{S_s}{K} \frac{\partial h}{\partial t} \qquad \text{(Diffusion equation for groundwater)}$$

Case 2: Steady-state flow (\(\partial h / \partial t = 0\)):

$$ \nabla^2 h = \frac{\partial^2 h}{\partial x^2} + \frac{\partial^2 h}{\partial y^2} + \frac{\partial^2 h}{\partial z^2} = 0 \qquad \text{(Laplace equation)} $$

Case 3: 2D horizontal flow in a confined aquifer (Dupuit assumption, integrating over thickness \(b\), \(T = Kb\), \(S = S_s b\)):

$$ T \left( \frac{\partial^2 h}{\partial x^2} + \frac{\partial^2 h}{\partial y^2} \right) = S \frac{\partial h}{\partial t} $$

EXAM NOTATION GUIDE

\(S_s\) = specific storage (per unit thickness), \(S = S_s b\) = storativity (dimensionless), \(T = Kb\) = transmissivity. The 3D equation uses \(S_s\); the 2D plan-view equation uses \(S\).

Hydraulic Head

The driving force in all these equations is the hydraulic head, which is the sum of elevation head and pressure head:

Components of Hydraulic Head Datum (z = 0) Aquifer P z Elevation head p/γ Pressure head h Total head h = z + p/γ (GW flows from higher h to lower h)
Figure 2.5: Hydraulic head = elevation head + pressure head. Groundwater flows from higher to lower total head.

2.4 VALIDITY OF DARCY'S LAW AND NON-DARCIAN FLOW

PAST EXAM QUESTIONS ADDRESSED IN THIS SECTION:

  1. Explain the importance, validity, and limitations of Darcy's law in groundwater motion. [Part of 6] (2078 Chaitra, Q2a)
  2. For a fully penetrating well in a confined aquifer, determine the valid domain around the well for which Darcy's law applies. [6] (2078 Chaitra, Q2b)
  3. A fully penetrating well is pumped at a constant rate from a confined aquifer. Determine the domain around the well for which Darcy's law is applicable. [8] (2073 Bhadra, Q2b)
  4. For a homogeneous isotropic confined aquifer, state whether Darcy's law is applicable and determine the average flow velocity in pores. [8] (2071 Bhadra, Q6)

Darcy's law is valid when flow through pores remains essentially laminar and viscous resistance dominates inertial effects. In groundwater engineering, this condition is commonly checked by a Reynolds number defined with mean grain diameter and either Darcy flux or seepage velocity, depending on the adopted convention.

$$ Re = \frac{q d_m}{\nu} \quad \text{or} \quad Re = \frac{v_s d_m}{\nu} $$

where \(d_m\) is average grain diameter and \(\nu\) is kinematic viscosity. Many groundwater exam problems state the limiting Reynolds number directly, so the accepted procedure is to use the form implied in class notes and solve for the region where \(Re \leq Re_{cr}\).

Darcy (Linear) vs Non-Darcian (Nonlinear) Flow Hydraulic gradient, i Specific discharge, q DARCY (linear) NON-DARCIAN Transition Re ≈ 1–10 slope = K Beyond Re ≈ 1–10, inertial effects cause deviation from Darcy linearity
Figure 2.6: Darcy’s law is valid only in the laminar range; at high gradients, inertial effects cause nonlinear (non-Darcian) behavior.
Darcy Validity Domain Around a Pumping Well Well DARCY VALID r ≥ rₘᵢₙ DARCY VALID r ≥ rₘᵢₙ Non-Darcy High velocity Re > Reᶜᵣ rₘᵢₙ rₘᵢₙ = Q·dₘ / (2π·b·ν·Reᶜᵣ) Darcy’s law valid for r ≥ rₘᵢₙ from the well center
Figure 2.7: Near a pumping well, radial velocity is highest; Darcy’s law is valid only beyond rₘᵢₙ where Re ≤ Reᶜᵣ.

STANDARD EXAM PROCEDURE FOR WELL-DOMAIN PROBLEMS

  1. Write radial Darcy flux for a fully penetrating confined aquifer: \( q_r = \frac{Q}{2 \pi r b} \).
  2. Use the Reynolds-number criterion prescribed in the problem.
  3. Set \( Re = Re_{cr} \) to obtain the boundary of validity.
  4. Solve for \( r = r_{min} \). Darcy flow is accepted for \( r \geq r_{min} \).
$$ Re = \frac{q_r d_m}{\nu} = \frac{Q d_m}{2 \pi r b \nu} $$
$$ r_{min} = \frac{Q d_m}{2 \pi b \nu Re_{cr}} $$

If class notes define Reynolds number with seepage velocity instead of Darcy flux, divide by effective porosity accordingly. The problem statement or the instructor's convention should control the final substitution.

Examples of Non-Darcian Flow in the Subsurface

  • Very coarse gravel packs near high-capacity wells.
  • Flow through large fractures and karst conduits.
  • Turbulent entrance losses in well screens and gravel envelopes.
  • Rapid infiltration through macropores, root channels, and fissured soils.

2.5 HYDRAULIC CONDUCTIVITY, PERMEABILITY, AND TRANSMISSIVITY

PAST EXAM QUESTIONS ADDRESSED IN THIS SECTION:

  1. Explain the continuum and REV approach in groundwater flow analysis. Define permeability and transmissibility. [4+2] (2073 Bhadra, Q2a)
  2. Define specific storage, specific yield, transmissibility, and anisotropy of an aquifer. [4] (2072 Ashwin, Q1)
IMPORTANT DISTINCTION: Permeability is a property of the porous medium alone, whereas hydraulic conductivity depends both on the medium and on the fluid flowing through it.

Hydraulic Conductivity, K

Hydraulic conductivity is the discharge that occurs through a unit area of aquifer under unit hydraulic gradient. It measures how easily a particular fluid can move through the porous formation. It depends on pore-size distribution, connectivity, and also on fluid density and viscosity.

Constant-Head Permeameter Overflow (constant level) Supply Soil Sample Collect Q Δh L K = QL / (AΔh)
Figure 2.5c: Constant-head permeameter — water flows through the soil sample under constant head difference. K is calculated from measured Q.

Intrinsic Permeability, k

Intrinsic permeability is a medium property independent of the fluid. It reflects the internal geometry of pores and is often expressed in m2 or in Darcy units in petroleum engineering.

Falling-Head Permeameter h₁ (at t = 0) h₂ (at t = t) falls a (pipe area) Soil Sample (A, L) Area = A L Outlet K = (aL)/(At) · ln(h₁/h₂)
Figure 2.5d: Falling-head permeameter — used for fine-grained soils. Head drops from h₁ to h₂ over time t.
$$ K = \frac{k \rho g}{\mu} = \frac{k g}{\nu} $$

where \(\rho\) is fluid density, \(g\) is gravitational acceleration, \(\mu\) is dynamic viscosity, and \(\nu\) is kinematic viscosity.

Transmissivity, T

Transmissivity is the integrated transmitting capacity of the full saturated thickness of an aquifer. It is especially useful in pumping-test and well-hydraulics calculations.

$$ T = K b $$

Specific Storage, Ss

Specific storage is the volume of water released from (or taken into) storage per unit volume of aquifer per unit decline (or rise) in hydraulic head, while the aquifer remains fully saturated. It arises from two mechanisms: compression of the aquifer skeleton and expansion of the water itself under pressure change.

$$ S_s = \rho g (\alpha + n \beta) $$

where \(\alpha\) = compressibility of the aquifer skeleton, \(n\) = porosity, and \(\beta\) = compressibility of water. Units: [m−1]. The storativity (or storage coefficient) of a confined aquifer is \(S = S_s b\), where \(b\) is the aquifer thickness (dimensionless).

Specific Yield, Sy

Specific yield is the volume of water released from storage per unit surface area of an unconfined aquifer per unit decline in the water table, due to gravity drainage of pores. It is the counterpart of storativity for unconfined aquifers.

$$ S_y = \frac{\text{Volume of water drained by gravity}}{\text{Total volume of aquifer}} $$

Typical values: coarse sand 0.25–0.35, fine sand 0.10–0.20, clay 0.01–0.05. Specific yield is always less than porosity because some water is retained by surface tension (specific retention: \(S_r = n - S_y\)).

Material Typical K (m/s) Groundwater implication
Clay 10-11 to 10-9 Acts as confining layer or aquiclude in practice
Silt 10-9 to 10-6 Slow drainage and leakage-dominated flow
Fine sand 10-6 to 10-4 Moderate aquifer behavior
Coarse sand 10-4 to 10-3 Productive shallow aquifer material
Gravel 10-3 to 10-1 Excellent water-bearing formation, though non-Darcy effects may appear near wells
Fractured rock Highly variable Controlled by fracture frequency and connectivity rather than matrix pores
Intrinsic Permeability → Hydraulic Conductivity → Transmissivity k Intrinsic Permeability Medium property only (pore geometry, grain size) Independent of fluid Units: [m²] × ρg/μ K Hydraulic Conductivity Medium + fluid property (K = kρg/μ = kg/ν) Flux per unit gradient Units: [m/s] × b T Transmissivity Whole aquifer thickness (T = Kb) Used in pumping tests Units: [m²/s] k (medium only) → K (medium + fluid) → T (K × saturated thickness b)
Figure 2.9: Conceptual distinction between intrinsic permeability (k), hydraulic conductivity (K), and transmissivity (T).

2.6 AQUIFER HETEROGENEITY AND ANISOTROPY

PAST EXAM QUESTIONS ADDRESSED IN THIS SECTION:

  1. Define specific storage, specific yield, transmissibility, and anisotropy of an aquifer. [4] (2072 Ashwin, Q1)
  2. Explain the continuum and REV approach in groundwater flow analysis. Define permeability and transmissibility. [Part of Q2a] (2073 Bhadra)

A groundwater formation is homogeneous if the hydraulic properties remain the same from place to place. It is heterogeneous if the properties vary spatially. A formation is isotropic if conductivity is the same in all directions at a point, and anisotropic if conductivity depends on direction.

Condition Meaning Groundwater example
Homogeneous Same K value everywhere in the region considered Uniform clean sand deposit
Heterogeneous K changes from point to point Alternating sand and silt lenses
Isotropic K is same in all directions at one point Randomly packed granular medium
Anisotropic K differs with direction at one point Layered alluvium where \(K_h > K_v\)
Anisotropy in Layered Deposits Layered Alluvial Deposit Sand (K₁ = high) Silty clay (K₂ = low) Sand (K₃ = high) Clay (K₄ = very low) Kₕ (large) Kᵥ (small) Conductivity Ellipse Kₕ Kᵥ Ellipse flattened → more anisotropic Always: Kₕ ≥ Kᵥ in stratified sediments Horizontal flow uses most permeable layer; vertical limited by least permeable
Figure 2.10: Anisotropy develops in layered deposits because horizontal flow encounters less resistance than vertical flow across bedding.

In many alluvial aquifers, the horizontal conductivity exceeds the vertical conductivity because sedimentary bedding creates easier pathways parallel to layering. This anisotropy strongly influences recharge, leakage, and pumping-test interpretation.

EXAM LANGUAGE TO REMEMBER

Heterogeneity refers to variation from point to point, whereas anisotropy refers to variation with direction at the same point. These two ideas are independent and can occur together.

Equivalent K for Layered Systems Flow PARALLEL to Layers (Kₕ) K₁, b₁ K₂, b₂ K₃, b₃ Kₕ = ∑(Kᵢbᵢ) / ∑bᵢ Weighted arithmetic mean Flow PERPENDICULAR (Kᵥ) K₁, b₁ K₂, b₂ K₃, b₃ Kᵥ = ∑bᵢ / ∑(bᵢ/Kᵢ) Weighted harmonic mean Always: Kₕ ≥ Kᵥ Arithmetic mean ≥ Harmonic mean — anisotropy ratio Kₕ/Kᵥ = 2 to 100+ in sedimentary deposits
Figure 2.10b: Equivalent K — arithmetic mean for parallel flow (Kₕ), harmonic mean for perpendicular flow (Kᵥ). Always Kₕ ≥ Kᵥ.

2.7 PAST QUESTION BANK FOR CHAPTER 2

The following list consolidates the chapter 2 past questions into theory and numerical categories. This is intended as a revision sheet that can be read independently from the detailed discussion above.

A. Theory and Definition Questions

2078 Chaitra Q2(a) Theory

Define Representative Elementary Volume and continuum approach of analysis. Explain the importance, validity, and limitations of Darcy's law in groundwater motion.

2073 Bhadra Q2(a) Theory

Explain the continuum and Representative Elementary Volume approach in groundwater flow analysis. Define permeability and transmissibility.

2072 Ashwin Q1 Definition

Define specific storage, specific yield, transmissibility, and anisotropy of an aquifer.

2072 Ashwin Q2 Derivation

Starting from the general expression of Darcy's law for groundwater flow, derive the three-dimensional form of groundwater flow.

2072 Magh Q2 Derivation

What is Darcy's law for groundwater flow? Derive the one-dimensional form and explain its validity.

2071 Bhadra Q2 Theory

What do you understand by Representative Elementary Volume? Explain.

B. Numerical and Applied Questions

2078 Chaitra Q2(b) Numerical

For a pumping well in a confined aquifer, determine the valid domain around the well for which Darcy's law applies using Reynolds number and kinematic viscosity data.

Aquifer thickness = 18 m
Pumping rate = 0.3 m^3/s
Average grain diameter = 1 mm
Reynolds limit = 8
Kinematic viscosity = 1 centistoke
2073 Bhadra Q2(b) Numerical

A fully penetrating well is pumped at a constant rate from a confined aquifer. Determine the domain around the well for which Darcy's law is applicable.

Pumping rate = 1000 m^3/hr
Aquifer thickness = 28 m
Average grain diameter = 1 mm
Reynolds limit = 6
Kinematic viscosity = 1 centistoke
2071 Bhadra Q6 Numerical

For a homogeneous isotropic confined aquifer with given porosity, permeability, piezometric heads, grain size, and viscosity, state whether Darcy's law is applicable and determine the average flow velocity in pores.

Aquifer thickness = 25 m
Effective porosity = 20%
Permeability = 15 m/day
Observation well spacing = 1100 m
Piezometric heads = 5.5 m and 3 m
Average grain diameter = 1 mm
Kinematic viscosity = 0.01 cm^2/s

C. Quick Revision Matrix

Theme Years asked What to prepare
REV and continuum 2078, 2073, 2071 Definition, necessity, and practical meaning in porous-media flow
Darcy law derivation 2072 Ashwin, 2072 Magh 1D and 3D forms, sign convention, and assumptions
Darcy validity 2078, 2073, 2071 Reynolds number, laminar range, radial-flow domain around a well
Hydraulic properties 2073, 2072 K, k, T, anisotropy, and unit conversions

2.8 WORKED NUMERICAL SOLUTIONS

SOLUTION 1: 2078 Chaitra Q2(b) — Darcy Validity Domain

Given: \(Q = 0.3\) m³/s, \(b = 18\) m, \(d_m = 1\) mm = \(0.001\) m, \(Re_{cr} = 8\), \(\nu = 1\) centistoke = \(1 \times 10^{-6}\) m²/s.

Required: Minimum radius \(r_{min}\) from the well beyond which Darcy's law is valid.

Step 1: For a fully penetrating well in a confined aquifer, the Darcy flux at radial distance \(r\) is:

$$ q_r = \frac{Q}{2\pi r b} $$

Step 2: Reynolds number defined using Darcy flux:

$$ Re = \frac{q_r \cdot d_m}{\nu} = \frac{Q \cdot d_m}{2\pi \cdot r \cdot b \cdot \nu} $$

Step 3: At the boundary of Darcy validity, \(Re = Re_{cr}\). Solve for \(r_{min}\):

$$ r_{min} = \frac{Q \cdot d_m}{2\pi \cdot b \cdot \nu \cdot Re_{cr}} $$

Step 4: Substitute values:

$$ r_{min} = \frac{0.3 \times 0.001}{2\pi \times 18 \times 1 \times 10^{-6} \times 8} $$
$$ r_{min} = \frac{3 \times 10^{-4}}{2\pi \times 144 \times 10^{-6}} = \frac{3 \times 10^{-4}}{9.0478 \times 10^{-4}} $$
$$ \boxed{r_{min} = 0.332 \text{ m} \approx 0.33 \text{ m}} $$

Conclusion: Darcy's law is valid for \(r \geq 0.33\) m from the well centre. Within 0.33 m, flow velocities are too high and non-Darcian effects dominate.

SOLUTION 2: 2073 Bhadra Q2(b) — Darcy Validity Domain

Given: \(Q = 1000\) m³/hr, \(b = 28\) m, \(d_m = 1\) mm = \(0.001\) m, \(Re_{cr} = 6\), \(\nu = 1\) centistoke = \(1 \times 10^{-6}\) m²/s.

Step 1: Convert units: \(Q = \frac{1000}{3600} = 0.2778\) m³/s.

Step 2: Apply the formula:

$$ r_{min} = \frac{Q \cdot d_m}{2\pi \cdot b \cdot \nu \cdot Re_{cr}} = \frac{0.2778 \times 0.001}{2\pi \times 28 \times 10^{-6} \times 6} $$
$$ r_{min} = \frac{2.778 \times 10^{-4}}{2\pi \times 168 \times 10^{-6}} = \frac{2.778 \times 10^{-4}}{1.0556 \times 10^{-3}} $$
$$ \boxed{r_{min} = 0.263 \text{ m} \approx 0.26 \text{ m}} $$

Conclusion: Darcy's law is applicable for \(r \geq 0.26\) m from the pumping well.

SOLUTION 3: 2071 Bhadra Q6 — Darcy Applicability and Pore Velocity

Given: \(b = 25\) m, \(n_e = 0.20\), \(K = 15\) m/day, \(L = 1100\) m, \(h_1 = 5.5\) m, \(h_2 = 3.0\) m, \(d_m = 1\) mm = \(0.001\) m, \(\nu = 0.01\) cm²/s = \(1 \times 10^{-6}\) m²/s.

Step 1: Hydraulic gradient

$$ i = \frac{h_1 - h_2}{L} = \frac{5.5 - 3.0}{1100} = \frac{2.5}{1100} = 2.273 \times 10^{-3} $$

Step 2: Darcy flux (specific discharge)

$$ q = K \cdot i = 15 \times 2.273 \times 10^{-3} = 0.0341 \text{ m/day} $$

Convert to m/s: \(q = \frac{0.0341}{86400} = 3.945 \times 10^{-7}\) m/s.

Step 3: Check Darcy validity (Reynolds number)

$$ Re = \frac{q \cdot d_m}{\nu} = \frac{3.945 \times 10^{-7} \times 0.001}{1 \times 10^{-6}} = \frac{3.945 \times 10^{-10}}{10^{-6}} $$
$$ \boxed{Re = 3.95 \times 10^{-4} \ll 1} $$

Conclusion: Since \(Re \approx 0.0004\), which is far below any commonly accepted limit (1–10), Darcy's law is fully applicable.

Step 4: Average pore velocity (seepage velocity)

$$ v_s = \frac{q}{n_e} = \frac{0.0341}{0.20} $$
$$ \boxed{v_s = 0.170 \text{ m/day} = 1.973 \times 10^{-6} \text{ m/s}} $$

The actual groundwater velocity in the pores is about 5 times the Darcy flux, and the flow is deeply in the laminar regime.

CHAPTER SUMMARY AND EXAM FORMULA SHEET

Concept Key relation / statement
Continuum assumption Porous medium is treated as continuous once a representative elementary volume is identified.
Darcy law \( Q = KAi \), or \( q = -K \frac{dh}{dl} \)
Seepage velocity \( v_s = q / n_e \)
3D Darcy law \( q_x = -K_x \frac{\partial h}{\partial x}, q_y = -K_y \frac{\partial h}{\partial y}, q_z = -K_z \frac{\partial h}{\partial z} \)
Vector form \( \vec{q} = -K \nabla h \) for isotropic media
1D flow equation \( K \frac{\partial^2 h}{\partial x^2} = S_s \frac{\partial h}{\partial t} \) (homogeneous); steady-state: \( \frac{d^2h}{dx^2} = 0 \)
3D flow equation (general) \( \frac{\partial}{\partial x}\!\left(K_x \frac{\partial h}{\partial x}\right) + \frac{\partial}{\partial y}\!\left(K_y \frac{\partial h}{\partial y}\right) + \frac{\partial}{\partial z}\!\left(K_z \frac{\partial h}{\partial z}\right) = S_s \frac{\partial h}{\partial t} \)
3D flow equation (isotropic) \( \nabla^2 h = \frac{S_s}{K} \frac{\partial h}{\partial t} \); steady-state: \( \nabla^2 h = 0 \) (Laplace equation)
Hydraulic head \( h = z + p/\gamma \)
Intrinsic permeability relation \( K = k \rho g / \mu = k g / \nu \)
Transmissivity \( T = K b \)
Specific storage \( S_s = \rho g (\alpha + n\beta) \); storativity \( S = S_s b \)
Specific yield \( S_y = V_{drained} / V_{total} \); always \( S_y < n \); \( S_r = n - S_y \)
Darcy validity test \( Re = q d_m / \nu \) or \( Re = v_s d_m / \nu \), with the adopted class convention
Radial confined-flow criterion \( q_r = Q / (2 \pi r b) \), so \( r_{min} = Q d_m / (2 \pi b \nu Re_{cr}) \) when Reynolds number is defined using Darcy flux
Anisotropy Directional conductivity varies at a point, commonly \(K_h > K_v\) in layered alluvium

LAST-MINUTE EXAM FOCUS

  1. Memorize the physical meaning of REV, not just its wording.
  2. Be able to derive 1D and 3D groundwater flow equations from mass conservation + Darcy substitution.
  3. Keep the distinction between Darcy flux and seepage velocity clear.
  4. Practice Reynolds-number based validity problems for radial well flow.
  5. Do not mix up hydraulic conductivity, intrinsic permeability, and transmissivity.