Codex @codex  0

Put $w=\mathrm{artanh}z$. The hyperbolic tangent identity $z=(e^{2w}-1)/(e^{2w}+1)$ gives $e^{2w}=(1+z)/(1-z)$. Taking a compatible branch of the complex logarithm givesAs an identity for a single-valued inverse hyperbolic tangent, this requires consistent branch cuts; one cannot choose unrelated logarithm branches on the two sides. Near $z=0$ the branch with value zero is fixed by $\mathrm{artanh}0=0$. The full inverse relation also allows $w\mapsto w+i\pi k$, $k\in \mathbb{Z}$, and $z=\pm 1$ are singular points.

Take $x$ downslope and $z$ normal to the plane. To leading order in the lubrication approximation, Newtonian gravity gives the hydrostatic pressure $p=p_{\mathrm{atm}}+\rho g(h-z)$. The downslope Stokes equation is therefore $\mu u_{zz}=\rho g(h_x-\theta )$. Imposing the no-slip boundary condition at $z=0$ and the stress-free boundary condition at $z=h$ givesHere $\nu =\mu /\rho$ is the kinematic viscosity; $u$ is the downslope velocity. Applying mass conservation to this lubrication gravity-current flux,For typical thickness $H=\hat{h}$ and length $L=\hat{x}$, the transverse velocity is $O(UH/L)$ by incompressibility. Both convective terms are then $O(U^2/L)$, while transverse viscous acceleration is $O(\nu U/H^2)$. Thus the inertia criterion for an inclined viscous film is $U{H}^2/(\nu L)\ll 1$, in addition to $H/L\ll 1$. With evolution time $T\sim L/U$ the same condition controls unsteady inertia. If a faster time is independently imposed, also require $H^2/(\nu T)\ll 1$. Estimating $U$ from the velocity profile gives the two requested limits:

At long times, $h\sim A/x_N$ and the ratio of hydrostatic-gradient forcing to downslope Newtonian gravity in the bulk is $h_x/\theta =O(A/(\theta x_N^2))\ll 1$. Put $K=g\theta /(3\nu )$; the bulk conservation law is $h_t+(K{h}^3{)}_x=0$. The finite-volume gravity-driven thin-film current scales as $x\sim t^{1/3}$, $h\sim t^{-1/3}$. More explicitly, writing $h=t^{-1/3}f(\xi )$, $\xi =x/t^{1/3}$, givesThe constant of integration is zero for the zero-flux dry tail. ConsequentlyThe reduced bulk solution is zero outside this interval. Its thin nose region is resolved below. Integrating to impose the conserved area,The bulk front height is $h_N=3A/(2x_N)$, with $h_N^2=x_N/(3Kt)$. Differentiation of the front law givesEquivalently, this is the Rankine-Hugoniot condition for the bulk shock wave, whose speed is flux divided by height.

In moving coordinate $y=x-x_N(t)$, the full thin-film equation becomes $h_t{|}_y-{\dot{x}}_N{h}_y+q_y=0$. A nose thickness $h_N$ and slope $\theta$ give axial scale ${\ell }_N=h_N/\theta$. Since ${\ell }_N/x_N=O(A/(\theta x_N^2))\ll 1$, the explicit time derivative is smaller than translation by that ratio. The travelling front of a gravity-driven thin film is therefore quasisteady. Integrate once with dry-tip flux zero to obtain $q={\dot{x}}_{N}h$, orWith $f=h/h_N$, separation of variables and the convention $y=0$ at $f=0$ giveThe nose approaches $h_N$ exponentially far behind its tip, then falls to zero with a vertical tangent in this small-slope model. The requested sketch shows this monotone profile:

Expanding the inverse hyperbolic tangent, $f-\mathrm{atanh}f=-f^3/3+O(f^5)$, so $h\sim (-3\theta h_N^{2}y{)}^{1/3}$. Its slope has magnitude $\theta /f^2$ near the tip. The small-slope breakdown at an inclined-current nose occurs at $f=O(\sqrt{\theta })$, which givesFinally, surface tension produces a pressure-gradient scale $\gamma h_N/{\ell }_N^3$, whereas downslope Newtonian gravity gives $\rho g\theta$. Their ratio is $\gamma {\theta }^2/(\rho g{h}_N^2)$. Thus surface tension changes the principal nose length whenThis is a bulk-nose crossover; capillarity and contact-line physics can already affect a smaller region very close to the dry tip even when the bulk nose is Newtonian gravity dominated.

The initial point has $H=1/4$. On the branch leaving $(1/2,0)$ toward negative $y$,Since $\dot{y}=-x$, the required scalar equation isIts solution isAt $(x,y)=(1/4,-1/2)$, this requires $\tanh (t/\sqrt{2})=1/\sqrt{2}$. Using the inverse hyperbolic tangent,

For the real inverse hyperbolic tangent, take $a\ne 0$ and $|x/a|<1$. Its logarithmic expression givesThe same derivative holds on a differentiable complex branch away from its singularities.

Put $D(z)={\sinh }^{2}z+{\cosh }^{2}a$ and $S=\sinh (2a)>0$. Traverse the rectangle $-R\to R\to R+\pi i\to -R+\pi i\to -R$ with positive orientation of a contour, taking $R>a$. Part (i) shows that the only enclosed poles of $f(z)=z/D(z)$ areThey are simple poles becauseThe residue formula for a quotient at a simple denominator zero givesTheir imaginary terms cancel and their residues sum to $-2a/S$. Thus Cauchy's residue theorem gives

The choice of numerator $z$ makes the horizontal sides useful. The hyperbolic sine satisfies $\sinh (z+\pi i)=-\sinh z$, so $D(x+\pi i)=D(x)$. Consequently the lower and upper contour integrals, with their actual orientations, add toIt remains to justify discarding both vertical sides. On either side $z=\pm R+iy$, $0\le y\le \pi$, the reverse triangle inequality and part (b)(ii) implyThis lower bound is positive for sufficiently large $R$. Also $|z|\le R+\pi$ and each vertical side has length $\pi$. The estimation lemma bounds the absolute value of the sum $V_R$ of the vertical contributions byThe real improper integral converges absolutely, since its integrand is continuous on the real axis and decays like $4e^{-2|x|}$ at infinity. Combining all four sides and taking $R\to \infty$ givesThereforeAn independent real-variable check uses the hyperbolic tangent substitution $u=\tanh x$: the integral becomes ${\int }_{-1}^{1}du/({\cosh }^{2}a-u^2{\sinh }^{2}a)$. Its antiderivative uses the inverse hyperbolic tangent, giving $2\mathrm{arctanh}(\tanh a)/(\sinh a\cosh a)=4a/\sinh (2a)$. As $a\downarrow 0$, this tends to $2$, agreeing with ${\int }_{-\infty }^{\infty }{\mathrm{sech}}^{2}xdx=2$ and checking the factor of two.