Past exam of the topics in fluid mechanics course of the University of Oxford / 2023 / 2 (Aquifer recharge and drainage)

Take $z$ vertically upward and the porous skeleton stationary. Darcy's law and incompressibility areHere $\mathbf{u}$ is the Darcy velocity, a volume flux per unit total cross-sectional area, including solid and void. If $\mathbf{v}$ is the liquid pore velocity and all the pore space is saturated and connected, thenThe porosity $\varphi$ therefore enters storage and the conversion to pore velocity; it does not multiply the Darcy velocity again when calculating discharge.

Let $H\ll L$ be the characteristic saturated thickness and let $\delta =H/L$. Scale pressure by $\rho gH$. The horizontal part of Darcy's law gives the scaleBy incompressibility, $w\sim \delta u\sim K{\delta }^2$. The vertical part of Darcy's law isIts left side is smaller than the scale $K$ by ${\delta }^2$. Therefore $p_z+\rho g=O(\rho g{\delta }^2)$, and the leading pressure is hydrostatic pressure:Here the water table has pressure $p_a$. This is the Dupuit approximation, valid away from the small region near an outlet where the surface slope need not remain small.

The hydrostatic pressure gives $p_x=\rho g{h}_x$, independent of depth. Thus Darcy's law gives $u=-K{h}_x$ and the horizontal volume flux per unit width isThe coefficient $K$ is the hydraulic conductivity, with units of velocity. For a column of unit horizontal area, the stored water volume is $\varphi h$, and rainfall supplies volume flux $R$. The local mass conservation equation is $\varphi h_t+q_x=R$. Therefore the Boussinesq equation for an unconfined aquifer isThis assumes rainfall is entirely converted into recharge and neglects unsaturated storage and losses before the water reaches the water table.

Set $x=L\tilde{x}$, $h=H\tilde{h}$, and $t=T\tilde{t}$. The Boussinesq equation for an unconfined aquifer becomesBalance rainfall recharge with lateral drainage, and then with storage:Dropping tildes, the dimensionless equation isThe horizontal Darcy velocity, vertical Darcy velocity, and discharge scales are respectively $\sqrt{KR}$, $R$, and $RL$. Thus the vertical thickness scale is set by rainfall balanced against hydraulic conductivity over the finite drainage length $L$, rather than by porosity. The latter determines the storage time $T$. The Dupuit approximation requires $H/L=\sqrt{R/K}\ll 1$ for these scales. The initially dry initial condition is $h(x,0)=0$.

Use the dimensionless equation from part (a). For the similarity ansatz $h=t^{a}f(\eta )$, $\eta =x/t^b$, its three terms have time factorsA nontrivial similarity solution retaining storage, lateral drainage and nonzero rainfall requires the exponents to agree:Unlike a source-free diffusion equation, the nonzero prescribed rainfall fixes an amplitude-versus-time scale directly, so no additional conserved quantity or boundary condition is required to determine these exponents. Degenerate spatially uniform solutions are an exception to this all-terms balance. The rain-fed aquifer similarity solution satisfiesIf the dimensional equation is used directly without scaling out coefficients, its similarity profile instead obeys $\varphi (f-\eta f^{\prime })=K(f{f}^{\prime }{)}^{\prime }+R$. The normalization above is the one in which the subsequent condition $f(\infty )=1$ applies.

The condition $f(0)=0$ means the river maintains the water table at the impermeable base level: $h(0,t)=0$. The liquid is exposed to atmospheric pressure there. It does not imply zero drainage, because $h_x$ may be singular at the outlet.

Far from the river, lateral drainage has not significantly affected storage. Rainfall simply fills pore volume, so the dimensional height is $h\sim Rt/\varphi$. In the dimensionless variables this is $h\sim t$, giving $f(\infty )=1$.

Although the physical aquifer is finite, these boundary conditions describe its early similarity solution before the drainage-affected region reaches $x=L$. In dimensional variables its horizontal extent is of order $\sqrt{KR}t/\varphi$, so the semi-infinite approximation requires $t\ll T=\varphi L/\sqrt{KR}$.

Introduce the positive outlet constantFor a draining similarity profile, $J$ is finite and positive. A regular linear profile would have $f{f}^{\prime }\to 0$ and $(f{f}^{\prime }{)}^{\prime }\ge 0$ at the origin, incompatible with the rain-fed aquifer similarity solution $(f{f}^{\prime }{)}^{\prime }(0)=-1$. Instead there is a square-root drainage boundary. At leading order, $(f^2{)}^{\prime }=2f{f}^{\prime }\to 2J$, givingMore explicitly, the rain-fed aquifer similarity solution implies $(f{f}^{\prime }{)}^{\prime }=-1+O(\sqrt{\eta })$, hence $f{f}^{\prime }=J-\eta +O({\eta }^{3/2})$ and $f^2=2J\eta -{\eta }^2+O({\eta }^{5/2})$. The height vanishes but its derivative diverges; their product remains finite.

Use $t_d,x_d,h_d$ temporarily for dimensional quantities. The nondimensionalization in part (a) givesTaking the drainage flux to be positive out of the left-hand end, $Q_0=-q(0,t_d)$, givesThe printed expression in terms of $f(0)$ and $f^{\prime }(0)$ must be understood as this limit; a literal product of their separate endpoint values is undefined. This is discharge per unit transverse width. Very close to the river, the divergent slope marks the limit of the Dupuit approximation, while the outer solution still supplies the leading discharge.

Return to dimensional variables. The river at base level imposes the drained boundary condition, while the right-hand end is impermeable:As long as the right-hand end remains wet, the latter is $h_x(L,t)=0$. The initial condition for the drought is the height profile inherited at $t=t^{\ast }$; stopping rainfall changes the forcing but does not instantaneously remove stored water.

With $R=0$, substitute $h=X(x)T(t)$ into the dimensional Boussinesq equation for an unconfined aquifer:For a positive draining profile, take $\lambda >0$. Normalize $T(t^{\ast })=1$, so $X$ has the dimensions of height. The temporal ordinary differential equation integrates toThe spatial ordinary differential equation and boundary conditions areHere $\lambda$ has units of inverse length with this normalization. The scale freedom between $X$ and $T$ is fixed by the chosen value $T(t^{\ast })=1$; the amplitude of $X$, and hence $\lambda$, must be matched to the available water.

This separable aquifer drawdown is exact if the inherited height profile has the separated spatial shape. For a general rainfall-generated initial condition, it describes the late drought regime; it cannot reproduce an arbitrary profile at the instant rainfall stops merely by choosing one temporal constant.

Let $J_X={\lim }_{x\downarrow 0}X(x)X^{\prime }(x)>0$. The spatial ordinary differential equation gives $(X{X}^{\prime }{)}^{\prime }=-\lambda X$, so its leading outlet behavior isAgain this is a square-root drainage boundary. For the separable aquifer drawdown, the dimensional drainage flux isEquivalently, with ${\tau }_d=\varphi /(K\lambda )$, it is $Q_0(t^{\ast })[1+(t-t^{\ast })/{\tau }_d{]}^{-2}$ for an exact separated profile. Integrating the spatial ordinary differential equation gives $J_X=\lambda {\int }_0^{L}Xdx$. This verifies the global mass conservation checkThus the late drainage flux follows an inverse-square power law, while stored water decays inversely with time.

The aquifer drainage hydrograph starts at zero, rises linearly in the early rainfall similarity solution, and approaches the total recharge $RL$ if rainfall lasts long enough for a finite-length steady state. Indeed, steady mass conservation gives $q_x=R$ and $q(L)=0$, hence $Q_0=RL$. At $t^{\ast }$, the outlet flux stays continuous for a regular inherited flux profile, then recedes as water drains; the eventual decay is proportional to $(t-t^{\ast }+\text{constant}{)}^{-2}$. An arbitrary transition requires the full time-dependent partial differential equation, so matching an early rainfall formula directly to a separated drought formula is only a schematic approximation.

The curve illustrates these regimes and a possible $t^{\ast }$. Its numerical shape and the relative duration of the regimes are schematic; the early rainfall and late drought formulas apply in their respective asymptotic ranges.