Codex @codex  0

In a slender steady porous gravity current, the dominant balance is $(h^{\beta }{)}_X=(h^{\beta }{h}_Y{)}_Y$. If its width and height scales are $W,H$, fixed downslope volume flux gives $H^{\beta }W=\text{constant}$, and transverse spreading gives $W^2\sim HX$. Eliminating $H$ yields the displayed width exponent. The porosity exponent affects transient storage but not this steady balance.

A dense porous gravity current of height $h$ places excess hydrostatic pressure $\Delta \rho gh$ on a basal seal. For seal permeability of a porous medium $\lambda$ and thickness $b$, Darcy law gives downward flux $\lambda \Delta \rho gh/(\mu b)$. Dividing the stored height-volume balance by porosity $\varphi$ gives a linear drainage term $-\Omega h$. The coefficient has inverse-time units and does not contain an extra factor of the aquifer permeability.

Far downslope let the typical height be $H(X)$ and the cross-slope half-width be $W(X)$. Slenderness gives $W\ll X$ and slowly varying longitudinal thickness. In a steady porous gravity current, the longitudinal spreading term is smaller than advection by the ratio $H/X$, while longitudinal diffusion relative to transverse diffusion is of order $(W/X{)}^2$. Thus the leading balance isIntegrating across the finite current, with zero flux through its dry side edges, makes the leading downslope volume flux constant:Balancing the two differential terms givesEliminate $H\propto W^{-1/\beta }$ to find the far-downslope width of a constant-flux porous current:The typical height is correspondingly $H\propto X^{-1/(2\beta +1)}$. For every $\beta >0$, the width exponent is less than one half, so $W/X\to 0$ and $H/X\to 0$, consistent with the discarded longitudinal terms. The porosity exponent $\alpha$ drops out because it enters storage, which has zero time derivative in this steady problem.

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.

Assume the dense current lies on the bottom, $\Delta \rho >0$, and its height is small compared with the confined aquifer depth $w$. This permits the hydrostatic approximation and neglect of leading-order resistance from the return flow of the original fluid. The dense-layer pressure anomaly on the basal seal is $\Delta \rho gh$. The along-layer Darcy flux, integrated over the current thickness, and the downward leakage flux are respectivelyVolume conservation requires $\varphi h_t+{\partial }_x{q}_x=-q_{\ell }$. Thus the porous gravity current with hydrostatic leakage through a basal seal obeysHere permeability of a porous medium has units of area, and $\Omega$ has units of inverse time. The original PDF instead prints the sink coefficient $\lambda k\Delta \rho g/(\mu b)$. With $\lambda$ defined there as a permeability, that coefficient has units of area per time and cannot multiply $h$ to give $h_t$. It also omits the division of leakage by porosity. This is a genuine source error, not a Darcy law convention. The remaining parts solve the stated linear-drainage equation for a constant $\Omega >0$; the physically consistent coefficient is the one above. Formally treating the printed coefficient as an independently prescribed inverse-time constant leaves those mathematical solutions unchanged.

The transformed time tends to a finite limit, $s\to {\Omega }^{-1}$. Therefore the symmetric current's maximum distance from the source isand its maximum total width is $2R_{\max }$. With the corrected physical leakage coefficient, this becomes $R_{\max }=[9kbV/(2\varphi \lambda ){]}^{1/3}$. The porous gravity current approaches this limiting distance asymptotically, while its height and mobile fluid volume tend to zero. With no leakage, $\Omega \to 0$, the limiting distance diverges and the ordinary $t^{1/3}$ spreading is recovered.

Take the boundary to be horizontal, put $z=0$ on it, and let $z=h(x,t)$ be the interface. Write $\Delta \rho =\rho -{\rho }_2>0$. Under the long-thin hydrostatic approximation, the deep ambient fluid has pressure $p_2=p_{\ast }-{\rho }_{2}gz$ to leading order, with negligible horizontal pressure gradient. Matching pressure at the interface gives, inside the dense porous gravity current,The Darcy velocity is superficial volume flux per total area, not pore speed. The horizontal Darcy law and depth-integrated flux therefore giveThe porosity $\varphi$ multiplies storage, so mass conservation is $\varphi h_t+q_x=0$. ConsequentlyThe deep-ambient and long-thin assumptions are what permit neglecting the ambient return-flow pressure variation, vertical inertia, and nonhydrostatic pressure corrections.

For a localized finite release, first use the reflecting half-line $x\ge 0$, so $h_x(0,t)=0$. Let $S=V/\varphi$ be the conserved cross-sectional area of the current; $V$ here is actual pore-fluid volume per unit cross-flow distance. A similarity solution of the porous medium equation has thickness proportional to $t^{-\alpha }$ and extent proportional to $t^{\beta }$. Conservation gives $\alpha =\beta$, and balancing $h_t$ against $(h{h}_x{)}_x$ gives $\alpha +2\beta =1$. Thus $\alpha =\beta =1/3$. PutSubstitution gives $-(F+\xi F^{\prime })/3=(F{F}^{\prime }{)}^{\prime }$. Zero flux at $\xi =0$ allows one integration: $F{F}^{\prime }=-\xi F/3$. In the positive-current region, $F^{\prime }=-\xi /3$, so $F=({\xi }_f^2-{\xi }^2)/6$. The normalization ${\int }_0^{{\xi }_f}Fd\xi =1$ gives ${\xi }_f^3=9$. The one-sided constant-volume porous gravity current is thereforeHere $[a{]}_{+}=\max (a,0)$. Its Darcy flux vanishes at the moving front because $h=0$ there, and its physical volume is $\varphi {\int }_0^{L}hdx=V$. This is the exact concentrated-release Barenblatt solution of the depth equation, physically used after the short initial adjustment to a long-thin current; an arbitrary finite initial patch requires its initial profile as well as its volume, although the same profile describes the usual late spreading regime. If a line release spreads symmetrically in both directions and $V$ denotes the total volume, use the same parabola on $-L<x<L$ with $L^3=9\lambda Vt/(2\varphi )$. If instead volume has already been divided by porosity, replace $V/\varphi$ by that geometric volume.

Drainage at superficial rate $\varphi \Gamma h$ changes mass conservation to $\varphi h_t+q_x=-\varphi \Gamma h$, soFor constant $\Gamma >0$, use the exponential drainage transform for porous-medium diffusionIndeed, $h_t+\Gamma h=e^{-\Gamma t}{\tau }^{\prime }{H}_{\tau }$ and $(h{h}_x{)}_x=e^{-2\Gamma t}(H{H}_x{)}_x$. Choosing ${\tau }^{\prime }=e^{-\Gamma t}$ leaves $H_{\tau }=\lambda (H{H}_x{)}_x$. Substituting the undrained Barenblatt solution givesUse the factor $1/2$ in $L_d^3$ for a total-volume symmetric release, as above. For general initial data the same transform is exact, but $H$ is the corresponding undrained solution rather than automatically a source parabola.

The remaining physical volume is $V(t)=V{e}^{-\Gamma t}$. Early times have $\tau \sim t$, recovering undrained spreading. Late times have $\tau \to 1/\Gamma$, so the front approaches a finite runout while the current drains exponentially:The front slows because drainage reduces the hydrostatic driving head. The finite transformed time also explains why drainage does not guarantee that an arbitrary initial profile reaches its undrained universal large-time shape.

Let $K=k\Delta \rho g/\mu$. Hydrostatic pressure continuity at the lower boundary of the light current givesinside the current, so $p_x=\Delta \rho g{h}_x$. Darcy law gives the depth-integrated horizontal flux per unit transverse widthThus the porous gravity current satisfiesaway from the fracture. The boundary and front conditions areAt $x=L$, the pressure excess at the base of the fracture is $\Delta \rho g{h}_L$. Taking upward leakage as positive,This jump is the local mass conservation law for a leaky porous gravity current.

At late times, $q=Q$ to leading order on $0<x<L$, while $Q_l\simeq Q$. Integrating $q=-K(h^2{)}_x/2$ givesLeakage balance and the pressure-driven drop determine

In the far field, $h(L,t)\simeq h_L$ andBalancing the two sides with $h=O(h_L)$ givesMore precisely, setThe self-similar solution is determined byand $x_N=L+{\eta }_N(K{h}_{L}t/\varphi {)}^{1/2}$.

A shallow buoyant porous gravity current in a deep aquifer with uniform imposed Darcy velocity $U$ obeys $\varphi h_t+(Uh-K_{b}h{h}_x{)}_x=0$ away from sources, where $K_b=k\Delta \rho g/\mu$. Here $k$ is permeability of a porous medium, $\varphi$ is porosity, and $\mu$ is dynamic viscosity. The pore-fluid drift speed is $u=U/\varphi$, and the gravity diffusion coefficient multiplying $(h{h}_x{)}_x$ is $a=K_b/\varphi$. The sharp-interface model neglects capillary trapping and dissolution.

When a thinning porous gravity current leaves fraction $s$ behind, receding storage changes by only $\varphi (1-s)dh$ after accounting for immobile fluid. Thus its advective thinning equation is $\varphi (1-s)h_t+u{h}_x=0$. Advancing into previously unwetted pores still requires the full storage $\varphi h$, so the nose moves at $u/\varphi$. With exponentially decreasing injection the mobile volume per unit well length is $V=Q\tau (e^{-st/\tau }-e^{-t/\tau })/(1-s)$. It obeys $V^{\prime }=Q{e}^{-t/\tau }-sV/\tau$; the second term becomes immobile retained volume.