Past exam of the mathematics course of the University of Cambridge / 2006 / iii / Paper 83 / 2 / Solution

Use constant porosity $\varphi$, isotropic rock permeability of a porous medium $K$, common liquid dynamic viscosity $\mu$ and gravitational acceleration $g$. Treat the two liquids as separated by a sharp interface, with negligible capillary pressure, molecular mixing, ambient-current resistance and inertia. At distances beyond the initial source scale the current is shallow relative to its horizontal length, so its pressure is hydrostatic. The initial release is symmetric about $x=0$ and has actual liquid volume $V$ per unit length of source. Consequently its geometric thickness integral is $V/\varphi$, not $V$.

Below the interface the excess pressure relative to the ambient hydrostatic field is $p^{\prime }=\Delta \rho g(h-z)$, so Darcy law gives horizontal discharge $q_x=-K\Delta \rho g{h}_x/\mu$. The vertically integrated discharge is $Q_x=-K\Delta \rho gh{h}_x/\mu$. Take the specified downward leakage $v=\gamma h$ to mean Darcy flux per bulk horizontal area. Mass conservation then givesIf $v$ is instead defined as interstitial seepage speed, the actual leakage discharge is $\varphi v$ and the corresponding coefficient is $\lambda =\gamma$. Stating the velocity convention removes that otherwise ambiguous porosity factor. For a shale of thickness $d_s$ and permeability $K_s$, a simple hydrostatic leakage model has $\gamma =K_s\Delta \rho g/(\mu d_s)$, if the reservoir below it has the same ambient pressure reference and no capillary entry threshold.

Let $R(t)$ be the half-width of the support; the full end-to-end length will be $L(t)=2R(t)$. Appropriate conditions are symmetry $h_x(0,t)=0$, $h(\pm R,t)=0$, no flux beyond the support, and the moving-front condition $\dot{R}=-\beta h_x(R^{-},t)$. Initially $\varphi h(x,t)dx$ converges weakly to $V{\delta }_0$, expressing the concentrated line release. The equations describe the spreading porous gravity current, not the singular immediate release itself.

With $\lambda =0$, conserved volume and balance of $h_t$ against $\beta (h{h}_x{)}_x$ give width proportional to $t^{1/3}$ and height proportional to $t^{-1/3}$. Write $h=t^{-1/3}f(\eta )$, $\eta =x/t^{1/3}$. Substitution givesIntegrating from the symmetry axis gives $\beta f{f}^{\prime }=-\eta f/3$. Inside the support $f>0$, so $f^{\prime }=-\eta /(3\beta )$ and $f=A-{\eta }^2/(6\beta )$. Thus the compact Barenblatt solution isThe positive-part notation makes $h=0$ outside the current. Volume fixes its remaining constant:If “length” denotes the distance from the source to one front, it is $R$ rather than $L$; both conventions are given explicitly. Differentiating gives $\dot{R}=R/(3t)=-\beta h_x(R^{-},t)$, so the front condition holds. Using $V/2$ on each reflecting half-line gives the same normalization.

For positive leakage define actual retained volume $\mathcal{V}(t)=\varphi \int hdx$. Since the edge thickness and boundary flux vanish, integration of the equation yieldsPut $h=e^{-\lambda t}H(x,\tau )$, where $\tau =(1-e^{-\lambda t})/\lambda$. Because $\dot{\tau }=e^{-\lambda t}$, substitution cancels the leakage term and leaves $H_{\tau }=\beta (H{H}_x{)}_x$. This derives the exponential drainage transform for porous-medium diffusion. Its conserved geometric volume is still $V/\varphi$, soThe full length is again $L(t)=2R(t)$. In particular,At early times $\tau \sim t$, recovering impermeable-bed spreading. At late times, drainage reduces the current thickness and therefore its hydrostatic driving pressure; the effective spreading time saturates at $1/\lambda$. The current approaches a finite horizontal reach while its retained volume and thickness vanish exponentially. This does not mean a positive-volume pool remains at a finite final depth. For small positive $\lambda$ the length scale diverges as ${\lambda }^{-1/3}$, consistently with unbounded spreading when leakage is absent.