Codex @codex  0

Integrate the axial ordinary differential equation twice: $u=-G{y}^2/(2\mu )+C_{1}y+C_2$. The no-slip boundary condition at zero gives $C_2=0$, and the one at $h$ gives $C_1=U/h+Gh/(2\mu )$. HenceThe first term is Couette flow; the second is Poiseuille flow driven in positive $x$ when $G>0$. Their addition is valid because the reduced Stokes equation is linear.

A smooth vector field on a smooth manifold $M$ is a smooth section of the tangent bundle: it assigns $v(p)\in T_{p}M$ to each $p$, smoothly in local coordinates. In a coordinate chart it has the form $v={\sum }_i{a}_i\partial /\partial x_i$, with smooth coefficients $a_i$. Equivalently, it acts on smooth functions as a derivation, $v(fh)=v(f)h+fv(h)$.

For a Lie group $G$, write $L_a(h)=ah$ for Left translation on a Lie group. A left-invariant vector field satisfiesThus it is determined by its value at the identity. Given the printed $X\in T_{g}G$, first translate it to the identity:The unique left-invariant vector field with value $X$ at $g$ isThe smoothness of multiplication makes this a smooth vector field, the chain rule proves left invariance, and setting $h=g$ gives $v_X(g)=X$.

We prove the completeness of left-invariant vector fields: this left-invariant vector field is a complete vector field. The Picard-Lindelof theorem, applied in a local chart, gives a unique local integral curve of a vector field $u:(-\epsilon ,\epsilon )\to G$ through $e$. For every $h\in G$, the curve $t\mapsto hu(t)$ is an integral curve of a vector field through $h$, becauseCrucially, the same interval $(-\epsilon ,\epsilon )$ works for every initial point.

Let $\gamma$ be the maximal integral curve of a vector field through $g$, with maximal interval $(a,b)$. If $b<\infty$, choose $t_0\in (a,b)$ with $b-t_0<\epsilon /2$. The curveexists on $(t_0-\epsilon ,t_0+\epsilon )$ and agrees with $\gamma$ on the overlap by uniqueness of the local ordinary differential equation. It extends $\gamma$ past $b$, a contradiction. The same argument at a finite $a$ excludes that possibility. Thus $(a,b)=\mathbb{R}$, and uniqueness on overlapping intervals gives uniqueness on all of $\mathbb{R}$.

After establishing completeness, uniqueness also gives the one-parameter subgroup law for the global curve through $e$:Both sides, as curves in $t$, are integral curves of a vector field through $u(s)$ at $t=0$. Defining the Exponential map of a Lie group by $u(t)=\exp (t\xi )$, the answer is

The identity component $G^{\circ }$ is an open normal subgroup, and every open identity neighbourhood generates it. Let $G^{\circ }$ be the connected component of $e$. The product of connected spaces is connected, so the image of $G^{\circ }\times G^{\circ }$ under multiplication is connected and contains $e$. It is therefore contained in $G^{\circ }$. Inversion has the same property. Thus $G^{\circ }$ is a subgroup.

A smooth manifold is locally connected. In particular, a coordinate neighbourhood of $e$ can be chosen homeomorphic to an open ball, so there is a connected open neighbourhood $O$ of $e$ contained in $G^{\circ }$. Its translates $hO$, $h\in G^{\circ }$, are open and lie in $G^{\circ }$, and cover $G^{\circ }$. Therefore the identity component of a Lie group is open in $G$. Conjugation by any $a\in G$ is a homeomorphism fixing $e$, so it maps $G^{\circ }$ into itself; conjugation by $a^{-1}$ gives the reverse inclusion. Hence $G^{\circ }$ is a normal subgroup.

If $U$ is an open neighbourhood of $e$ in $G^{\circ }$, let $H=⟨U⟩$, allowing inverses in the meaning of generated subgroup. For each $h\in H$, $hU$ is open and lies in $H$, so $H={\bigcup }_{h\in H}hU$ is open in $G^{\circ }$. Every other left coset of $H$ is also open. Thus $H$ is both open and closed in the connected space $G^{\circ }$. It is nonempty, so $H=G^{\circ }$.

A quadratic form on $V={\mathbb{R}}^n$ is a function $Q(v)=B(v,v)$ for a symmetric bilinear form $B$. Equivalently, $Q(tv)=t^{2}Q(v)$ and the polarizationis a bilinear map. In coordinates there is a unique real symmetric matrix $A$ with $Q(v)=v^{T}Av$. No assumption of nondegeneracy or positive definiteness is needed.

An element $g\in {\mathrm{GL}}_n(\mathbb{R})$ stabilizes $Q$ when $Q(gv)=Q(v)$ for every $v$. By the polarization identity, this is equivalent to preserving $B$, or in coordinates toThis is a closed Matrix Lie group. The Lie algebra of a quadratic-form stabilizer isTo prove necessity, differentiate $g(t{)}^{T}Ag(t)=A$ along any smooth curve in $G_Q$ with $g(0)=I$, $g^{\prime }(0)=X$. The derivative at zero is $X^{T}A+AX$.

To prove sufficiency, suppose $X^{T}A+AX=0$ and consider the matrix exponential $g(t)=e^{tX}$. ThenIts value at zero is $A$, so $e^{tX}\in G_Q$ for every real $t$. This curve has derivative $X$ at zero, proving the claimed tangent space description even for a degenerate quadratic form.

Equivalently, the condition is $B(Xu,v)+B(u,Xv)=0$ for all $u,v$. For a nondegenerate quadratic form it is the corresponding Special orthogonal Lie algebra. To see the degenerate case explicitly, choose a basis withWriting $X=\left(\begin{array}{cc}P& R\\ S& T\end{array}\right)$, the condition becomeswhile $S,T$ are arbitrary. Thus the nullspace of the quadratic form is invariant, but arbitrary infinitesimal maps into it are allowed. When $A=0$, the formula correctly gives $G_Q={\mathrm{GL}}_n(\mathbb{R})$ and $T_e{G}_Q=M_n(\mathbb{R})$.

Along the downward gradient flow, the gradient-flow dissipation identity, the chain rule and the defining property of the Riemannian gradient giveA smooth vector field on a compact manifold is complete, so the trajectory exists for every $t\ge 0$. In particular, its function value never increases.

If the trajectory is still in $C_{[a,b]}$ at time $t$, its monotone function value has stayed in $[a,b]$ throughout $[0,t]$. ConsequentlyThus a uniform exit bound isThe trajectory is outside the closed band for every $t>T$. Once it leaves below level $a$, monotonicity prevents its return. The argument also covers the zero-width band $a=b$.

For the total occupation-time bound for a gradient flow, let $m={\min }_{M}f$ and $F={\max }_{M}f$, which exist by compactness. For every $R>0$, the same gradient flow identity givesTaking $R\to \infty$, by monotone convergence theorem, gives the total occupation-time boundThis allows arbitrarily many visits to $C_{\epsilon }$; it bounds their total duration. Neither of these two estimates needs the Morse function hypothesis.

For exponential convergence of a Morse gradient flow, first observe that a limiting point $z=x_{\infty }$ is a critical point. Indeed, writing ${\varphi }_s$ for the flow and using continuous dependence on initial data,Differentiation at $s=0$ gives $\nabla f(z)=0$.

We use the Morse lemma: near a nondegenerate critical point of Morse index $\lambda$, there are coordinates $(u,v)\in {\mathbb{R}}^{\lambda }\times {\mathbb{R}}^{n-\lambda }$ in whichChoose a Riemannian metric equal to the Euclidean metric on a smaller ball in these coordinates. Such a global Riemannian metric exists by a partition of unity: patch this metric to any background metric outside a slightly larger ball. One can do this simultaneously at every critical point, since a Morse function on a compact manifold has only finitely many critical points. The metric is chosen before forming its gradient flow; changing the metric can change trajectories.

For this metric, any trajectory converging to $z$ eventually remains in the smaller coordinate ball. From some time $t_0$ onward its equations are exactlyConvergence forces $u(t_0)=0$, and thenThe straight radial segment stays in the coordinate ball, so its length bounds the Riemannian distance from above:Increase $C$ beyond both $e^{2t_0}|v(t_0)|$ and ${\max }_{0\le t\le t_0}{e}^{2t}d({\gamma }_x(t),z)$. Then the requested strict estimate holds for all $t\ge 0$:This proof uses the Morse lemma, existence and uniqueness for smooth ordinary differential equations, continuous dependence on initial data, and the existence of a partition of unity. The expanding coordinates explain why convergence to a saddle occurs only along its stable manifold.

For a counterexample, use an arc coordinate $x$ around a point of the circle, its flat local metric, and a globally smooth function that equals $x^4$ on this arc. Such a function can be made by a smooth cutoff, agreeing with a positive constant outside a larger arc. For $x_0>0$ sufficiently small the downward gradient flow remains in the arc and isIt converges to the degenerate critical point $0$, and its Riemannian distance from $0$ is $x(t)$ for sufficiently large $t$. For every $k>0$, $e^{kt}x(t)\to \infty$, so no exponential bound is possible. This counterexample is on a compact manifold as well.

Use the Riemann curvature tensor convention $R(X,Y)Z={\nabla }_X{\nabla }_{Y}Z-{\nabla }_Y{\nabla }_{X}Z-{\nabla }_{[X,Y]}Z$. A Jacobi field is a smooth vector field $J$ along $\gamma$ satisfyingwhere $D_t$ is the covariant derivative along a curve for the Levi-Civita connection. Differentiating a geodesic variation gives this equation, and conversely every Jacobi field arises from a geodesic variation by varying its initial point and velocity.

Choose a parallel frame $E_1(t),\dots ,E_n(t)$ along $\gamma$ and write $J={\sum }_a{y}_a{E}_a$. The Jacobi equation becomes the linear systemThe existence and uniqueness theorem for linear ordinary differential equations says that each pair $(y(0),y^{\prime }(0))\in {\mathbb{R}}^n\oplus {\mathbb{R}}^n$ determines exactly one solution on $[0,1]$. Addition and scalar multiplication preserve the equation. Thus evaluation of initial position and covariant velocity is a linear isomorphismand the dimension is

For the convexity assertion, use the convex-normal meaning of a geodesically convex open set: its points are joined by a unique geodesic within the set, and the joining geodesic depends smoothly on its endpoints. Equivalently, the appropriate star-shaped restriction of ${\exp }_p$ is a diffeomorphism onto the set. The convex-normal-neighbourhood theorem supplies such sets around every point. Mere existence of some minimizing geodesic, without this uniqueness and normality, would not imply the assertion.

Let $p=\gamma (0)$ and $v=\dot{\gamma }(0)$. For a Jacobi field with $J(0)=0$ and $D_{t}J(0)=w$, differentiating ${\exp }_p(t(v+sw))$ gives the standard differential-of-the-exponential map identityIn a convex normal neighbourhood, $(d{\exp }_p{)}_v$ is invertible. Therefore a Jacobi field vanishing at both endpoints has $w=0$ and hence is identically zero. The difference of two fields with the same endpoint values consequently vanishes. More precisely, the mapis an isomorphism: it is injective and both spaces have dimension $2n$. The constant-geodesic case follows directly from $D_t^{2}J=0$. This uses the convex-normal-neighbourhood theorem, the differential-of-the-exponential map identity, and linear ordinary differential equation uniqueness.

For the special unitary group example, use Jacobi fields from conjugation at a central endpoint. WriteThe endpoints of $\gamma (t)=\exp (tA)$ are $I$ and $zI$. The latter is central. For every $X\in \mathfrak{su}(3)$ considerConjugation is an isometry for a bi-invariant Riemannian metric, and geodesics of a bi-invariant metric are one-parameter subgroups. Thus $F_X$ is a geodesic variation. Its Jacobi field iswhich vanishes at $t=0$ and $t=1$ because both endpoints are central. Its initial covariant derivative is $[X,A]$.

The linear map $X\mapsto J_X$ has kernel equal to the centralizer of an element of a Lie algebra of $A$. Indeed, $J_X\equiv 0$ implies $[X,A]=0$ on differentiating at $t=0$, and that commutation conversely implies $J_X\equiv 0$. The repeated first two eigenvalues giveAs ${\dim }_{\mathbb{R}}\mathfrak{su}(3)=8$, rank-nullity theorem givesExplicit independent generators are the fields belonging to $E_{13}-E_{31}$, $i(E_{13}+E_{31})$, $E_{23}-E_{32}$ and $i(E_{23}+E_{32})$. Their initial derivatives $[X,A]$ are independent because the eigenvalue differences in the $13$ and $23$ entries are $2\pi i$. This also demonstrates directly why the endpoint-value conclusion fails along this geodesic.

The deterministic SIR model follows the ordinary differential equationsUse $\beta \ge 0$ and $\gamma >0$. The state space is $S,I,R\ge 0$, $S+I+R=N+1$; summing the three equations proves conservation of this population. Infection is supported only where $S>0$ and $I>0$, and recovery only where $I>0$. The rates vanish at the corresponding boundaries, preventing negative compartment sizes.

This convention has an unnormalized population infection rate $\beta SI$, as required by the printed threshold. A convention with infection rate $\beta SI/(N+1)$ would use a differently scaled parameter.

Orient both spheres in the standard way and normalize their area forms to total area $4\pi$. The degree of a map between oriented manifolds can be obtained in two ways. For a regular value $y$, use the degree as a sum of local degrees:Each inverse image is isolated by the inverse function theorem, and compactness makes the set finite. The determinant is computed in consistently oriented local coordinates. A second method is spherical degree by area pullback:where the last formula represents $g$ as a unit vector in ${\mathbb{R}}^3$. The pullback of a differential form already contains the signed Jacobian determinant; no extra $\sin \theta$ is to be inserted in that last coordinate expression.

To relate the methods, replace $\omega$ by a smooth top-degree differential form with the same total integral supported in a small neighbourhood of a regular value. Two such top-degree forms with equal integral differ by an exact differential form on $S^2$, by its top-degree de Rham cohomology. Their pullbacks therefore have the same integral by Stokes theorem. Over the chosen neighbourhood, $g$ splits into local inverse branches; the change of variables formula makes the contribution of each branch its orientation sign times $4\pi$. Their sum is precisely the first formula. Thus the area integral is an integer and agrees with the signed inverse-image count.

For a nonconstant rational map, first use common-factor reduction of a rational map so $p$ and $q$ are coprime. Write $k=\max (\mathrm{deg}p,\mathrm{deg}q)$ for these reduced polynomials. A generic finite target value $w$ has inverse images at the roots of $p-wq$: avoiding exceptional values makes its degree $k$ and its roots simple. The fundamental theorem of algebra supplies $k$ roots. A holomorphic map has positive real Jacobian determinant $|R^{\prime }{|}^2$ at a regular point, so every local sign is $+1$. HenceThe source leaves coprimality implicit. In an unreduced representation the answer is $\max (\mathrm{deg}p,\mathrm{deg}q)-\mathrm{deg}\gcd (p,q)$, including degree zero for a constant reduced map. For example $(z^2-1)/(z-1)$ extends to $z+1$ and has degree one, although the unreduced maximum degree is two. Exceptional inverse images at infinity or multiple roots do not change the degree of a rational map of the Riemann sphere.

For the rational map approximation for Skyrmions, use stereographic projection $z=\tan (\theta /2)e^{i\phi }$ and the unit target vectorCombine this rational map with a radial profile to form a special unitary group field:where $\mathbf{\sigma }$ are the Pauli matrices. The endpoint values make $U(0)=-1$ independent of angle and $U(\infty )=1$. Appropriate radial behaviour gives an admissible finite-energy field configuration. With $L_i=U^{†}{\partial }_{i}U$, choose the topological baryon number in the Skyrme model conventionSeparating the radial and angular factors givesThus the degree of a rational map of the Riemann sphere supplies the Skyrmion charge.

In conventional dimensionless massless Skyrme model units, its static energy reduces towith the angular Jacobian of a rational mapThe Cauchy-Schwarz inequality gives $\mathcal{I}\ge k^2$. These formulas follow from the radial strain $|f^{\prime }|$ and the two equal angular strains $\sin f\sqrt{J_R}/r$: the quadratic energy sums their squares and the quartic Skyrme term sums their pairwise products of squares. Minimize the angular integral in the rational map approximation over degree-$k$ maps, then minimize the remaining radial energy with the stated endpoints. This replaces a three-dimensional field minimization by finitely many map coefficients and an ordinary differential equation for $f$.

The method constructs a charge-$k$ variational approximation, with topology built in and with rotational symmetry of a rational map translated into combined spatial and isospin rotations. It is efficient for identifying shapes and providing initial data for unrestricted numerical relaxation. Its restrictions are equally concrete: it uses one radial profile and a holomorphic angular map independent of radius, so it cannot represent arbitrary radial-angular correlations, separated clusters, or all deformations. Apart from the degree-one Skyrmion hedgehog ansatz, it generally does not solve the full field equation exactly. Massive-pion terms can be included in the radial functional but do not remove these restrictions, and multi-shell or unrestricted fields may be needed for larger charges. Approximate energy minima and a final collective-coordinate quantization are distinct steps.

Take $z$ positive downwards from the roof. Write $T_0=T_m-m{C}_0$, $C_B=(T_m-T_B)/m$, $\Delta T=T_0-T_B=m(C_B-C_0)>0$ and $\Delta C=C_B-C_0$. The cooled magma follows the decreasing liquidus from $(C_0,T_0)$ towards $(C_B,T_B)$; the initially superheated liquid starts at $(C_0,T_{\infty })$ above that liquidus. The solid has zero solute concentration. At each temperature in the mushy layer, a horizontal tie line joins the pure solid to the residual liquid on the liquidus. The lever rule gives the liquid fraction from the bulk concentration. Since $T_B>T_E$, this path stops before the eutectic temperature and does not cross the eutectic composition.

The porosity $\varphi$ denotes the liquid fraction, so $1-\varphi$ is the solid fraction. Equal phase mass densities remove volume-change flow; the stable density stratification suppresses convection. In the stagnant mushy-layer model, define bulk specific enthalpy and bulk solute concentration byThe arbitrary reference of specific enthalpy has been omitted. The two conservation laws in $0<z<h(t)$ areHere $\kappa =k_T/(\rho c_p)$ is the thermal diffusivity, with $k_T$ the thermal conductivity. The positive $L{\varphi }_t$ term records increased specific enthalpy on melting; during freezing it supplies latent heat. In the liquid $z>h(t)$, the heat equation and solute diffusion equation arewhere $D$ is the solutal diffusivity. Neglecting solute diffusion inside the mushy layer does not permit dropping $D{C}_{zz}$ in the liquid at this stage.

Use $T(0,t)=T_B$, liquid $T\to T_{\infty }$ and $C\to C_0$ as $z\to \infty$, and initially $h=0$, $T=T_{\infty }$, $C=C_0$ for $z>0$. These are the half-space boundary conditions appropriate before cooling reaches the bottom of the chamber. For a closed chamber of finite depth $H$, replace the far-field conditions by the specified bottom conditions, for example $T_z=C_z=0$ at an insulated impermeable floor; the half-space similarity solution then ceases to apply when the thermal penetration depth approaches $H$.

At $z=h(t)$, let $C_i$ and $T_i=T_m-m{C}_i$ be the common liquid concentration and temperature, and let ${\varphi }_i$ be the limiting porosity on the mush side. Continuity of $T,C$, the Rankine-Hugoniot condition for solute, and the interfacial conservation of energy giveThe first condition has the sign of solute rejection into the liquid. For the sharp-edge stagnant mushy-layer model with nonzero liquid solutal diffusivity, close the free-boundary problem by marginal equilibrium at a mush-liquid boundary:This prevents constitutional supercooling immediately ahead of the mushy layer. A small solid fraction jump at its edge is allowed: imposing ${\varphi }_i=1$ as well as these finite-$D$ conditions would generally overdetermine the reduced model. At the roof, $C=C_B$ follows from the liquidus, and no additional solute-flux boundary condition is needed in the diffusion-free mush.

To obtain the large-concentration similarity solution, putThe positive superheat specifies the physically intended liquid initial state. On a similarity solution, $C_i,{\varphi }_i$ are constant, so every newly incorporated level leaves the same constant $B={\varphi }_i{C}_i$. Write $C_i=C_0+\Delta C\delta$ and $B=\Delta C(\mathcal{C}+a)$. Then throughout the mushy layerFor $a,\theta =O(1)$, $\mathcal{C}\gg 1$, and $\mathcal{S}/\mathcal{C}=O(1)$, the effective heat capacity is constant to relative error $O({\mathcal{C}}^{-1})$. Thus setIntegration of the similarity solution ordinary differential equation gives the leading analytical profilesThe temperature is $T=T_0+\Delta T\theta$. In the mushy layer, the concentration and porosity are the displayed $C,\varphi$; in the liquid,These retain the liquid solutal boundary layer. For completeness, finite-$r$ values of $\delta ,a,\lambda$ are determined analytically, to the same large-$\mathcal{C}$ order, by settingand solvingThe solute jump gives ${\varphi }_i=1-b/(\mathcal{C}+\delta )$ exactly; replacing $\mathcal{C}$ by $\mathcal{C}+\delta$ on the right of the last equation retains that jump's next-order factor, but does not make the linearized thermal profiles exact at finite $\mathcal{C}$. The branch has $0\le \delta <1$, $0<\varphi \le 1$. All these profiles are leading asymptotic expressions, not exact solutions of the nonlinear finite-$\mathcal{C}$ heat equation.

For $r\to 0$ at fixed positive $\lambda$, the complementary error function asymptotic gives $K(\lambda /\sqrt{r})\sim \sqrt{\pi }\lambda /\sqrt{r}$, hence $\delta =O(r)$, $b=O(r)$, $a\to 0$, and ${\varphi }_i\to 1$. The result simplifies toIn particular, $\varphi (0,t)=C_0/C_B=\mathcal{C}/(\mathcal{C}+1)$ and $\varphi (h^{-},t)=1$. There is no finite interfacial latent heat jump in this limit: the latent heat is released continuously throughout the mushy layer. Matching the two temperature gradients givesThe printed final equation is missing the square on $\Omega$ in the exponential. The corrected exponent follows directly from differentiating $\mathrm{erf}(\Omega \eta )$ and is consistent with ${\Omega }^2=1+\mathcal{S}/\mathcal{C}$. The positive root is unique: the mush-side gradient decreases from infinity to zero with $\lambda$, whereas ${\theta }_{\infty }K(\lambda )$ increases. The zero-superheat limit is singular and is not described by a finite positive $\lambda$ in this two-region construction.

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.

Use the time convention $e^{-i\omega t}$. The free-space Helmholtz equation, with $\psi =e^{ikx}E$, gives the exact reduced equationFor a forward plane wave at angle $\alpha$, $E_x=O(k{\alpha }^{2}E)$ and $E_{zz}=O(k^2{\alpha }^{2}E)$, whereas $E_{xx}=O(k^2{\alpha }^{4}E)$. Thus the paraxial approximation neglects $E_{xx}$ relative to $2ik{E}_x$. The reduced field satisfies the parabolic wave equationWith the Fourier transform convention used below and its inverse, this becomes an ordinary differential equation for each transverse wavenumber:Consequently the initial-value solution isEquivalently, evaluating the oscillatory Gaussian integral gives the one-dimensional transverse Fresnel propagation formulaThe square-root branch has $\sqrt{1/i}=e^{-i\pi /4}$. The formula is an oscillatory integral for general data, or the Fresnel propagator acting on square-integrable functions. It approaches the initial field as $x\downarrow 0$. A plane wave has transverse wavenumber $\nu =k\sin \alpha$; the reduced axial phase $-{\nu }^{2}x/(2k)$ agrees with $k(\cos \alpha -1)x$ through order ${\alpha }^2$. This also checks the sign. The paraxial approximation applies to $|\nu |\ll k$.