Codex @codex  0

An instantaneous point release of actual volume $V$ per unit width in a porous gravity current with background flow has the Barenblatt solution $h=[R(t{)}^2-(x-ut{)}^2{]}_{+}/(6at)$, where $R^3=9aVt/(2\varphi )$. Its centre translates at pore speed $u=U/\varphi$ and its radius grows as $t^{1/3}$. The conserved physical volume is $\varphi \int hdx=V$.

A localized release of actual fluid volume $V$ per unit width on a reflecting half-line $x\ge 0$ has the Barenblatt solution $h=[L^2-x^2{]}_{+}/(6\lambda t)$, where $L^3=9\lambda Vt/\varphi$ and $\lambda =k\Delta \rho g/(\varphi \mu )$. The mass conservation normalization is $\varphi {\int }_0^{L}hdx=V$. Zero slope at the reflecting boundary gives zero Darcy flux. For symmetric spreading on the full line, the same formula has $L^3=9\lambda Vt/(2\varphi )$ when $V$ denotes the total released volume. A finite initial patch is not determined by volume alone; this exact source profile has a concentrated initial release.

Let $z$ measure height above the substrate, and write $\mu =\rho \nu$. In the lubrication approximation, the horizontal scale $\ell$ is much larger than the typical height $H_0$. Vertical momentum balance is hydrostatic to leading order, while the horizontal balance retains the dominant vertical shear:The no-slip boundary condition gives $u(0)=0$; negligible air viscosity and constant surface tension give $u_z(h)=0$. Integrating twice and then through the depth yieldsThe free-surface kinematic boundary condition and incompressibility imply $h_t+q_x=0$. ThereforeThe cross-sectional area $A=\int hdx$ is conserved when the endpoint volume flux vanishes.

The relevant small groups follow from the same balances. Put ${ϵ}_h=H_0/\ell$ and $U_{\ast }=g{H}_0^3/(3\nu \ell )$. Small ${ϵ}_h$ suppresses horizontal viscous derivatives and makes vertical viscous corrections to hydrostatic pressure smaller by $O({ϵ}_h^2)$. The ratio of horizontal inertia to vertical viscous stress isThus a small full-length Reynolds number is sufficient but not necessary: thin-layer geometry reduces the inertial ratio. The unsteady viscous ratio $H_0^2/(\nu t_{\ast })$ must likewise be small; on the spreading timescale $t_{\ast }=\ell /U_{\ast }$ it equals the displayed group. If surface tension is retained, its pressure-gradient correction relative to gravity isThe Bond number here uses the horizontal scale, because the capillary pressure is $O(\gamma H_0/{\ell }^2)$ whereas the hydrostatic pressure is $O(\rho g{H}_0)$. We also neglect air shear and evaporation and assume a Newtonian fluid with constant kinematic viscosity.

To find the similarity solution, put $K=g/(3\nu )$. Conservation of $A$ gives $h\sim A/\ell$, and balancing the time derivative against the nonlinear diffusion gives ${\ell }^5\sim K{A}^{3}t$. A normalized ansatz isSubstitution gives $(f^3{f}^{\prime }{)}^{\prime }=-(f+\xi f^{\prime })/5$. For the symmetric zero-flux profile, one integration gives $f^3{f}^{\prime }=-\xi f/5$, and on its positive support a second integration givesThus the solution is a compactly supported Barenblatt solution of the porous medium equation. Let $\ell (t)$ denote its half-width. In dimensional variables,Here the positive-part notation sets the height to zero outside $|x-x_c|<\ell$. The area normalization fixes the coefficient, not merely the exponent. With $x-x_c=\ell \sin \theta$,Writing the remaining integral as $I$, the final widths areThe height falls like $t^{-1/5}$ and the width grows like $t^{1/5}$. The volume flux tends to zero at the front, even though the profile slope diverges; this gives a valid compactly supported weak solution of the reduced equation. The diverging edge slope means the lubrication approximation is not uniform at a microscopic contact line. Surface-tension or contact-line physics supplies an inner regularization without changing the gravity-dominated bulk similarity law in its validity regime. This is an exact source-type solution for $t>0$; an arbitrary finite initial drop approaches an appropriate late-time profile rather than being exactly represented by this singular initial profile.

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.

For a concentrated line release, take $V$ to be actual fluid volume per unit source length and allow symmetric spreading to both sides of $x=0$. PutSubstitution in the porous medium equation gives the exponential drainage transform for porous-medium diffusion:The Barenblatt solution has conserved area ${\int }_{-\infty }^{\infty }Hdx=V/\varphi$. A quadratic cap satisfies this equation:Indeed, inside the current $H_x=-x/(3Ds)$ and $R\propto s^{1/3}$; differentiating verifies $H_s=D(H{H}_x{)}_x$. Its area is $\int Hdx=2R^3/(9Ds)=V/\varphi$, fixing the prefactor. The flux vanishes at the two moving noses because $H=0$ there. Hence the remaining physical volume and occupied width areFor a one-sided current on a reflecting half-line, replace the $2$ in $R^3$ by $1$ when $V$ denotes the total release into that side. Volume alone does not select this exact solution for an arbitrary finite initial patch; it describes an impulsive line release. Its initial central height is singular, so a finite injection region and the constraint $h\ll w$ regularize the near-source early-time behavior.

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.

Use the positive transport coefficient $D$ from the preceding solution. The porous medium equation here is $h_t=D(h^2{h}_x{)}_x=(D/3)(h^3{)}_{xx}$. A planar pulse conserves the water cross-sectional area $\mathcal{V}=\int hdx$. If the supplied lake volume is a three-dimensional volume $V$, introduce the constant out-of-plane width $B_y$ and take $\mathcal{V}=V/B_y$; alternatively $V$ may be understood as volume per unit span. A two-dimensional model cannot determine an absolute extent from an unspecified three-dimensional volume alone.

For a symmetric localized release, write $h=t^{-a}f(x/t^b)$. Area conservation requires $a=b$, and balancing the PDE powers gives $a+1=3a+2b$, hence $a=b=1/4$. With $\eta =x/t^{1/4}$, the profile equation isIntegrate once, using symmetry and zero central flux: $D{f}^2{f}^{\prime }=-\eta f/4$. Inside the wet region this gives $f^2=C-{\eta }^2/(4D)$. The dry continuation is zero. Write the result asThe normalization follows from ${\int }_{-L}^{L}hdx=\pi HL/2=\mathcal{V}$. This cubic diffusion pulse has finite support $-L<x<L$, total extent $2L$, and spreading speed $\dot{L}=L/(4t)$. Although $h_x$ is singular at the ideal nose, the flux vanishes there and $u=-Dh{h}_x=x/(4t)$ tends to the finite front speed.

For a one-sided pulse on $0<x<L$ with a reflecting boundary at zero and the same area $\mathcal{V}$, replace $\mathcal{V}$ in the displayed full-line formulas by $2\mathcal{V}$. The power laws are unchanged. These are source-type Barenblatt solutions for an instantaneous localized release; a finite initial lake footprint approaches the profile at long times and may require a virtual time origin, rather than matching this singular initial condition exactly.

Take $\Delta \rho >0$ and $\tau >0$, with $x$ pointing along the imposed shear stress, and put $G=\Delta \rho g$. Relative to the upper fluid's hydrostatic reference, the hydrostatic pressure in the current is $p=G(h-z)$. Thus the horizontal pressure gradient is $G\nabla h$. The lubrication theory momentum balance, lower no-slip boundary condition, and imposed upper shear stress giveIts depth-integrated volume flux isDefine $D=G/(3\mu )$ and $A=\tau /(2\mu )$. The continuity equation gives the shear-driven viscous gravity current equationThe two terms have opposite roles: imposed shear stress carries fluid downstream, whereas the hydrostatic pressure gradient spreads it from thick regions towards thin ones.

For the approximation, a representative horizontal velocity is $U\sim \tau H/\mu +G{H}^3/(\mu L)$, with vertical velocity $UH/L$. Small slopes require $H/L\ll 1$. Horizontal fluid inertia relative to vertical viscous resistance is ${\rho }_{c}U{H}^2/(\mu L)$, where ${\rho }_c=\rho +\Delta \rho$. Hence sufficient small parameters, expressed without an unknown velocity, areFor the hydrostatic pressure and normal-stress approximation also require $\mu U/(GHL)\ll 1$. Its gravity-driven part is already $(H/L{)}^2$; its shear-driven part addsThis condition controls viscous normal stress, vertical viscous corrections and the normal projection of shear at a slightly tilted interface relative to $GH$. With it and the previous inertia bounds, vertical inertia is small as well. Time variations here are on the transport timescale $L/U$; separately imposed rapid forcing would need its own unsteady inertia bound. A small ordinary Reynolds number is a stronger sufficient restriction, but the reduced inertia ratios above are what the thin-layer momentum balance directly requires.

In the steady far-downstream regime, cross-stream derivatives dominate gravity-driven spreading. Dropping the smaller downstream gravity flux leavesWrite $Y=y_N(x)$ and let $h_c$ be a typical central thickness. Balance gives $Y^2\sim (D/A)h_c^{2}x$, while conservation of the source volume flux gives $Q\sim A{h}_c^{2}Y$. Therefore the downstream similarity of a shear-driven gravity current hasThe omitted downstream gravity flux relative to the imposed shear flux is $D{h}_c^2/(Ax)\sim (Y/x{)}^2\ll 1$, so this approximation is self-consistent in the stated regime.

For the full profile put $w=h^2$. Then $h^3{h}_y=(w^2{)}_y/4$ andThis is a porous medium equation with downstream distance as its evolution coordinate. Conservation of $\int wdy=Q/A$ and a similarity solution $w=x^{-1/3}f(\eta )$, $\eta =y/x^{1/3}$, givewhere symmetry gives zero integration constant at $\eta =0$. Inside the positive region, $f^{\prime }=-2A\eta /(3D)$, so $f=(A/(3D))({\eta }_N^2-{\eta }^2)$. Requiring the total downstream volume flux to be $Q$ fixes the coefficient:Consequentlyfor $|y|<y_N$, with $h=0$ outside. The squared height is a parabolic Barenblatt solution; the height cross-section is a semicircular profile after rescaling its axes. The cross-stream volume flux vanishes at the edges even though the height slope becomes singular there. The profile describes the outer lubrication theory region, not a resolved microscopic front.

Near the point source the two horizontal dimensions are comparable, say $\ell$. Upstream spreading arrests where outward gravity-driven volume flux balances downstream shear-driven volume flux:Eliminating $h_{\ast }$ givesThis is a scaling estimate, not a determination of a numerical prefactor. The associated depth is $h_{\ast }\sim (\mu Q/(\Delta \rho g){)}^{1/4}$. The same horizontal scale is obtained by setting $y_N\sim x$ in the downstream similarity solution, confirming where that solution fails.

For the line source the steady volume flux per unit width is $q=A{h}^2-D{h}^3{h}_x$. On the downstream constant-height branch, $q=Q_{2d}$ gives $h_d^2=Q_{2d}/A$. Upstream there is no net flux through the finite nose, so $q=0$. Thus $h{h}_x=A/D$, and continuity of height at the source givesThe line-source upstream reach under imposed shear is thereforeIn the upstream part, shear-driven and gravity-driven volume fluxes cancel. The volume flux per unit width jumps by exactly $Q_{2d}$ at the source. As in the point-source profile, the square-root nose is a formal outer solution with an unresolved steep front; its finite upstream reach follows from flux balance, not from assuming a small slope all the way to the nose.

Choose the positive $x$ direction along the background flow, so $U>0$; reverse $x$ if necessary. Let $z$ measure depth below the impermeable roof, and set $\Delta \rho ={\rho }_w-{\rho }_c>0$. The deep ambient water has $p_w=p_0(x)+{\rho }_{w}gz$. Continuity of pressure at $z=h$ and hydrostatic pressure inside the buoyant CO2 giveHere $k$ is permeability of a porous medium and $\mu$ is dynamic viscosity. Since $h/H\ll 1$, the ambient return-flow correction to the imposed pressure gradient is higher order. The small aspect ratio supplies the hydrostatic pressure approximation; capillary trapping, dissolution, compressibility and dispersion are omitted in this sharp-interface model.

The depth-integrated Darcy flux is $q=Uh-K_{b}h{h}_x$. The actual stored fluid volume is $\varphi h$ per unit horizontal area, where $\varphi$ is the aquifer porosity. Thus the porous gravity current with background flow satisfieswith $u=U/\varphi$, $a=K_b/\varphi$ and $J=Q/\varphi$. Both $Q$ and the later $V$ denote actual injected volume per unit transverse width. Away from injection, the equation is an advected porous medium equation. At each finite front impose $h=0$ and vanishing outgoing volume flux per unit width. At the injection point $h$ is continuous and $q(0^{+})-q(0^{-})=Q$. For an initially empty aquifer, $\varphi \int hdx=Qt$. The advective speed is $U/\varphi$, not $U$.

For a constant-flux porous gravity current, write characteristic depth $H_c$ and extent $L_c$ temporarily. At early times volume conservation gives $H_c{L}_c\sim Jt$; gravity spreading gives $L_c^2/t\sim a{H}_c$. HenceBackground advection becomes comparable when $L_c\sim ut$, so the transition scales areFor $t\ll t_c$, the leading similarity solution is symmetric:For $0<\eta <{\eta }_N$, its model ordinary differential equation and boundary conditions areChoose the nonnegative finite-front branch, reflect evenly in $x$, and require $F{F}^{\prime }\to 0$ at the front. The factor $1/2$ divides the injection equally between the two halves. Integration of the ordinary differential equation gives $2{\int }_0^{{\eta }_N}Fd\eta =1$; the limiting front slope is $F^{\prime }({\eta }_N^{-})=-2{\eta }_N/3$. This boundary-value problem determines the initially symmetric shape without assuming an unjustified elementary formula. Independent numerical integration gives ${\eta }_N\simeq 1.17619$ in this normalization.

For $t\gg t_c$, upstream storage approaches a finite steady value. With no upstream leakage, $q=0$ for $x<0$, giving $h_x=U/K_b$ wherever $h>0$. Downstream the uniform interior carries $q=Q$, soThe upstream volume flux vanishes even though local Darcy velocity contributions from background flow and buoyancy cancel. The upstream stored volume is $\varphi h_{\mathrm{int}}|x_{-}|/2$, a finite constant. The downstream plateau therefore extends a distance $ut+O({\ell }_c)$, and its front advances at $u$ to leading order.

The diffusive nose of an advected porous gravity current has widthIn a frame translating at $u$, balance the time derivative against the nonlinear gravity diffusion. With $\eta =(x-x_c(t))/\sqrt{a{h}_{\mathrm{int}}t}$ and $x_c(t)=ut+O({\ell }_c)$, the leading transition is $h=h_{\mathrm{int}}F(\eta )$, withIts finite front has $F\sim ({\eta }_N/2)({\eta }_N-\eta )$ as $\eta \uparrow {\eta }_N$. Independent numerical integration of this boundary-value problem gives ${\eta }_N\simeq 1.23849$, so a more accurate leading nose position is $x_N\simeq x_c+1.23849\sqrt{a{h}_{\mathrm{int}}t}$ and its local shape is $F\simeq 0.619247(1.23849-\eta )$. The entire transition has no elementary closed form. A useful explicit, approximate shape is the volume-preserving linear rampThis gives $x_N\simeq x_c+\sqrt{a{h}_{\mathrm{int}}t}$. Its ramp width is $2\sqrt{a{h}_{\mathrm{int}}t}$, and its edge obeys the exact kinematic boundary condition ${\dot{x}}_N=u-a{h}_x(x_N^{-},t)$ to the displayed order. It also has the same area as a sharp step at $x_c$. This is an approximate ramp closure, not an exact solution of the nonlinear nose ordinary differential equation. A constant phase shift $x_c=ut-|x_{-}|/2$ balances the finite upstream storage within this ramp approximation; the full solution's subleading matching may alter constant offsets.

For a rapid injection of total volume $V$, remove the source for $t>0$ and introduce $\xi =x-ut$. The advected constant-volume porous gravity current obeys $h_t=a(h{h}_{\xi }{)}_{\xi }$ and $\int hd\xi =M=V/\varphi$. A similarity solution has depth proportional to $t^{-1/3}$ and radius proportional to $t^{1/3}$. Direct substitution gives the compactly supported Barenblatt solutionIndeed ${\int }_{-R}^R(R^2-{\xi }^2)/(6at)d\xi =2R^3/(9at)=V/\varphi$. Its centre translates at the pore speed $u$, its two edges are $ut\pm R(t)$, and its maximum depth decays as $t^{-1/3}$. The expression is exact for a point release within this reduced model; finite-width initial data approach this similarity solution at late times. If $U=0$, the constant-rate current remains in the symmetric spreading regime and the constant-volume solution has a stationary centre.

In a cross-section normal to the cylinder axes, measure $y$ from the line of contact. The two circular boundaries have the parabolic approximationfor their separation. The contact point of the meniscus is at $y=w+o(w)$, so the leading cross-sectional area isThe area of the small meniscus cap is $O(w^4/a^2)$ and is lower order. At its upper end the gap has width $h(w)=w^2/a$. A tangent semicircle therefore has radius $w^2/(2a)$ and curvature

The Young–Laplace equation makes the liquid pressure, relative to the nearly uniform gas pressure,At a fixed $y$, lubrication theory gives a planar Poiseuille flow through a gap of width $h(y)$, with axial flux per unit $y$Integrating across the cusp,The continuity equation $A_t+Q_z=0$ now gives

Set $q=w^3$. Since $w^4{w}_z=\frac{1}{5}(q^{5/3}{)}_z$, this is the porous medium equationConservation of the fixed volumeand dimensional analysis give the self-similar solution $q=t^{-3/8}F(z/t^{3/8})$. One integration of the resulting ordinary differential equation yields the compactly supported Barenblatt solutionEquivalently,Its tip is $z_N=(C/B{)}^{1/2}{t}^{3/8}$. Usingin the volume constraint gives $C^2=8aV\sqrt{B}/\pi$. Therefore

For vertical cylinders at equilibrium, hydrostatic pressure gives $p-p_{\mathrm{atm}}=-\rho gz$. Balancing this with the capillary pressure $-2\gamma a/w^2$ gives the large-height profile

When sliding dominates and the whole cap is below the snowline, its thickness obeysWith $\alpha =0$, this is a one-dimensional porous medium equation. Its mass-preserving Barenblatt solution isWithin its support,Now set $h=h_{ss}+ct$. Since $h_{ss}$ solves the unforced equation,Matching the uniform ablation term $-\alpha$ givesThe similarity solution has a virtual origin. Writing $s=t+t_0$ gives the exact decaying familyIf its initial centre thickness is $h^{\ast }$, then $C=(h^{\ast }+3\alpha t_0/4)t_0^{1/3}$ and the exact extinction time isIts value depends on the initial cap width through the virtual age $t_0$, but its dimensional scale is unambiguouslyThe additive melting correction alone has the corresponding time $4h^{\ast }/(3\alpha )$.