Calculate Darcy flux and volumetric flow

Separate the driving gradient from the medium’s response

A head difference provides a driving gradient, but the resulting flow also depends on the water-transmitting properties of the medium. This investigation assigns hydraulic conductivity K between 1 and 10 m/day. The selected homogeneous column has the same conductivity everywhere, a fixed length of 100 m and a bulk cross-section of 20 m².

Hydraulic conductivity depends on the porous material and fluid properties. It has dimensions of length per time but is not, by itself, the speed of a water parcel. Darcy’s law relates it to flow through a specified area under a head gradient. USGS explanation of conductivity and Darcy’s law.

The fictional setup assumes saturated, constant-density, one-dimensional laminar flow. Boundary heads are maintained. No transient adjustment, pumping well, variable layering or unsaturated behavior is calculated.

Calculate flux before multiplying by area

Let i be the signed driving head drop divided by length: i = (left head − right head)/100 m. Darcy flux q is K × i, with positive flow to the right. This is equivalent to q = −K dh/dx when dh/dx means the mathematical slope of head toward increasing x.

At left head 14 m and right head 10 m, i = 4/100 = 0.04. With K = 2 m/day, q = 2 × 0.04 = 0.08 m/day. The gradient is a ratio of lengths, so its units cancel. The resulting flux has length-per-time units.

Flux is volume per time per unit bulk cross-sectional area. Multiplying q by the full 20 m² gives the signed volumetric flow Q = 1.6 m³/day. Write the unit multiplication: m/day × m² = m³/day. Q answers how much water crosses the whole section per day; q normalizes that amount by bulk area.

Why flux is not pore-water speed

The bulk area includes solids as well as connected pore space. Dividing flow by that whole area produces a flux, not the average speed within the effective flowing pore space. The investigation’s effective porosity is 0.25, so its average pore velocity is q/0.25. At the default, that is 0.32 m/day, four times q.

Keep three labeled quantities in your notes: K = 2 m/day, q = 0.08 m/day and average pore velocity = 0.32 m/day. Their matching units do not make them interchangeable. Each has a different definition and role in the calculation.

Draw a bulk section as a rectangle and label an idealized quarter of it as effective pore area. This is an accounting sketch, not a realistic map of pore geometry. Use it to explain why the same total flow corresponds to a larger average velocity when normalized by the smaller effective area.

Test proportionality with one change at a time

Keep left head at 14 m and double K from 2 to 4 m/day. The driving gradient stays 0.04, q doubles to 0.16 m/day and Q doubles to 3.2 m³/day. The direction remains rightward. Capture both cases and identify which input was held constant.

Return K to 2 and change left head to 6 m. The gradient becomes −0.04, q becomes −0.08 m/day and Q becomes −1.6 m³/day. Negative signs encode right-to-left flow, not a negative amount of water physically present. Their magnitudes match the first case because the head-drop magnitude and conductivity match.

At equal boundary heads, q and Q are zero at every allowed positive K. This is a useful counterexample to the claim that high conductivity alone guarantees fast flow. The model needs a driving difference as well as a transmitting medium.

Transfer and evidence limits

For a paper-only case, take K = 3 m/day, driving gradient −0.02 and bulk cross-section 30 m². Calculate signed q and Q before checking. Then state what would be required to apply the result to a measured aquifer.

Check the Darcy-flow transfer

q = 3 × (−0.02) = −0.06 m/day and Q = 30 × (−0.06) = −1.8 m³/day. Flow is opposite the positive axis. Applying this relation to a real aquifer requires appropriate measurements and evidence that the geometry, conductivity, head gradient and flow assumptions are suitable; the arithmetic alone does not establish them.