Use pressure with the ambient hydrostatic pressure removed, and let $\Delta P$ be the pressure below the falling sphere minus that above it. Away from the sphere the tube flow is Hagen-Poiseuille flow. The upper and lower clear lengths add to $L$ to leading order in $a/L$, so their pressure drops add:The large exterior reservoir returns this flux with negligible pressure loss compared with the long narrow tube. The sphere is assumed far enough from the ends that entrance effects and its occupied length give only lower-order corrections.

Now use a frame translating downwards with the sphere, and take $z$ positive upwards from its center. The sphere is stationary and the tube wall moves upwards with speed $U$. The narrow gap isIt has axial length scale $a\sqrt{ϵ}$, making lubrication theory appropriate. The sphere-frame flux in a tube is constant: the upward relative flux through the whole section is $\pi a^{2}U-Q$, so its flux per unit circumference isWith $y=0$ on the sphere and $y=h$ on the wall, Couette-Poiseuille flow in a thin gap givesSince $\Delta P=-\int p_{z}dz$, the pressure-jump relation isThe parabolic inner gap can be extended to infinite $z$ because these integrals are dominated by $|z|=O(a\sqrt{ϵ})$. The lubrication resistance integrals for a nearly occluding sphere areSetting $\xi =\tan \theta$ evaluates $J_1=\pi$, $J_2=\pi /2$, $J_3=3\pi /8$, hence

For the force, enclose the fluid between sections just above and below the sphere, including the annular layer. The net upward pressure force is $\pi a^2\Delta P$. The tube wall exerts upward shear traction $\mu u_y(h)$ on this fluid; the sphere exerts a downward force equal to the upward hydrodynamic force $F$ on itself. From the gap solution,The axial control volume force balance therefore gives the control-volume drag formula for a sphere in a tube:Using a control volume avoids having to integrate the curved-sphere normal pressure and tangential shear separately.

With $Q^{\ast }=Q/(\pi a^{2}U)$, $\Delta P^{\ast }=a\Delta P/(6\mu U)$, $q^{\ast }=2q/(aU)$, $F^{\ast }=F/(6\pi \mu aU)$, $I_n^{\ast }=a^{n-1}{I}_n$ and $L^{\ast }=4L/(3a)$, the three flux/pressure equations becomeSolving these linear relations givesThe force relation is $F^{\ast }=\Delta P^{\ast }+\frac{4}{3}{I}_1^{\ast }-q^{\ast }{I}_2^{\ast }$. Substitution first yieldsBecause $I_2^{\ast }/I_3^{\ast }=2ϵ/3$ and $I_1^{\ast }/I_3^{\ast }=2{ϵ}^2/3$, the terms $-2L^{\ast }{I}_2^{\ast }$ and $\frac{4}{3}{L}^{\ast }{I}_1^{\ast }$ are lower order than $L^{\ast }{I}_3^{\ast }$. Keeping the uniformly relevant leading terms givesThe squared $I_2^{\ast }$ is essential; replacing it by $I_2^{\ast }{I}_3^{\ast }$ would give the wrong sign and scaling. This establishes the three drag regimes for a nearly occluding sphere.

For $L^{\ast }\ll {ϵ}^{-1/2}$,The sphere nearly carries the tube's fluid down with it: the far-field flux is close to the piston displacement rate. The relative annular leakage is small in total volume, and the dominant force comes from the local annular shear term. The pressure force associated with the overall tube resistance is smaller.

For ${ϵ}^{-1/2}\ll L^{\ast }\ll {ϵ}^{-5/2}$,The flow remains approximately piston-like, but the dominant force is now the pressure needed to drive Hagen-Poiseuille flow through the long clear tube. Both displayed terms in $q^{\ast }$ are retained because their relative size changes inside this regime without producing an additional leading drag regime.

For $L^{\ast }\gg {ϵ}^{-5/2}$,The large tube resistance makes its net through-flow negligible. The sphere's displaced fluid instead passes upwards relative to it through the thin annular gap, with gap velocity of order $U/ϵ$. Its pressure drop and resulting pressure force dominate.

Finally consider two identical co-moving spheres, with nonoverlapping lubrication regions and the same total clear-tube length to leading order. The common through-flux is not doubled. In the middle regime it is still approximately $\pi a^{2}U$, so the same total tube pressure drop is shared equally between the two identical spheres. The load sharing between co-moving spheres in a tube gives half the single-sphere force on each sphere in regime (ii). In regime (i), drag is set locally by the nearly pressure-balanced annular shear and is unchanged on each sphere. In regime (iii), each sphere must pass essentially its own displacement flux through its annular gap, giving the same local pressure drag as before; adding another pressure jump changes the already small tube flux but not either leading drag. The force is unchanged in regimes (i) and (iii).

At leading order in the lubrication approximation, surface arclength is horizontal distance. An insoluble surfactant is transported by the surface velocity, without bulk exchange or diffusion. Conservation of insoluble surfactant on a moving interface therefore becomesThe local hydrostatic pressure and capillary normal stress giveSolve $\mu u_{zz}=p_x$ with no-slip boundary condition at $z=0$ and Marangoni stress $\mu u_z(h)={\gamma }_x=-A{C}_x$. The velocity, surface velocity and volume flux per unit width areMass conservation, $h_t+q_x=0$, yields the thin-film equations with insoluble surfactant:and the corresponding surfactant equation isThe hydrostatic and capillary terms thus affect both the film flux and the surface transport, with different coefficients.

For finite-mass Marangoni spreading on a liquid film, first neglect both pressure-gradient terms. If the spread has size $\ell$, then $C\sim M/\ell$ and $C_x\sim M/{\ell }^2$, giving surface speed $u_s\sim AM{h}_0/(\mu {\ell }^2)$. Equating $\ell /t$ to this speed givesSet $\ell =(AM{h}_{0}t/\mu {)}^{1/3}$ exactly as a similarity scale, and writeOn the positive half of the pool, substitution produces two ordinary differential equations:Symmetry fixes the flux integration constant in the second equation to zero, giving $H{\Gamma }^{\prime }=-\eta /3$ wherever $\Gamma >0$. Substitution into the first equation gives $\eta H^{\prime }=H$. Thus the linear similarity profiles for surfactant spreading have the formThe concentration vanishes at the moving front. The two integral constraints are total surfactant mass and conservation of the fluid volume relative to the undisturbed layer:The first gives $k={\eta }_N^2/3$; the second gives $k{\eta }_N=2$. Therefore ${\eta }_N^3=6$. In physical variables, the solution isOutside the pool the leading outer solution has $h=h_0$, $C=0$. In particular, $h(x_{N-},t)=2h_0$: the subscript means the one-sided limit just behind the front. As a consistency check, the surface speed there is $x_N/(3t)={\dot{x}}_N$ and the fluid jump satisfies the moving-front mass conservation condition. This ideal outer solution has a height jump from $2h_0$ to $h_0$.

That jump cannot persist in a physical interface with nonzero gravity or surface tension. Across a smoothing region of width $\Delta$, $h$ changes by $O(h_0)$. The Marangoni flux scales as $AM{h}_0^2/(\mu x_N^2)$. At $g=0$, its balance with the capillary flux ${\gamma }_0{h}_0^4/(\mu {\Delta }^3)$ gives capillary smoothing of a Marangoni front:Here the surface tension approaches ${\gamma }_0$ at the surfactant-free edge. The assumed concentration-gradient scale is $M/x_N^2$ even inside the smoothing region.

When gravity dominates capillarity, the hydrostatic flux is $\rho g{h}_0^4/(\mu \Delta )$. The same balance gives gravity smoothing of a Marangoni front:The gravity-dominated condition is ${\Delta }^2\gg {\gamma }_0/(\rho g)$, the square of the capillary length. The printed condition omits $\rho$; the expression above restores dimensional consistency.

This width grows faster than $x_N$. It becomes comparable with the pool size whenThe time is a scaling estimate, so numerical factors cannot be fixed by the width balance. For $t\gg t^{\ast }$, the sharp-front similarity profile no longer describes the film: hydrostatic pressure levels the layer over the whole pool, and its thickness tends towards $h_0$ with an increasingly small depression beneath the surfactant. The surface can still spread by Marangoni stress while a hydrostatic-pressure-driven return flow nearly cancels the net film flux. In this gravity-levelled surfactant film, $q\simeq 0$ gives $h_x\simeq -3A{C}_x/(2\rho g{h}_0)$, so the fractional height variation is of order $AM/(\rho g{h}_0^2{x}_N)\ll 1$.

Take $z$ positive upwards and let $f$ denote axial body force per unit volume. With inertia and surface tension omitted, the extensional equations for a slender Newtonian column areMass conservation gives the first equation. The second uses the Trouton ratio three: the leading radial velocity is $u_r=-r{w}_z/2$, the radial normal traction gives $p=p_{\mathrm{ext}}-\mu w_z$, and the axial stress is $-p_{\mathrm{ext}}+3\mu w_z$. An axial control volume balance then gives the displayed equation. The approximation requires a slender column and slowly varying radius.

Let $h=R-a$. Since the container base is closed, the total volume flux at every section is zero. The column carries flux $\pi a^{2}w$, so the annulus carries $-\pi a^{2}w$ and has mean speed of order $a|w|/h$. Its interfacial shear stress is therefore of order $\lambda \mu a|w|/h^2$. Transmitting this stress through the core produces an axial velocity variation of order $\lambda a^2|w|/h^2$. The plug-flow criterion for a column in a low-viscosity annulus is consequently

Use modified pressure $P=p_{\mathrm{ext}}+\rho gz$, removing the annular fluid's hydrostatic pressure. In a local planar gap, with $y=0$ at the column and $y=h$ at the container wall, lubrication theory givesThe required return flux per unit circumference is $q_g=-aw/2$ to leading order in $h/a$. Its pressure-driven component is larger than the Couette component by $a/h$, so the annular return flow around a slumping column givesThe annular shear contribution to the integrated axial force is smaller than this pressure-gradient term by $h/a$. Retaining the leading terms and substituting $f=-(\rho +\Delta \rho )g$ gives the dimensional equationsThe boundary conditions are $w(0,t)=0$ at the base and $w_z(L(t),t)=0$ at the free upper end; the latter expresses zero excess axial extensional stress.

For $\lambda =0$ and initially constant $a=a_0$, put $G=\Delta \rho g/(3\mu )$. Then $w_{zz}=G$, so

For $\lambda >0$, defineThis function records the leading thin-annulus resistance; differences from curvature and the Couette term enter at the already neglected relative order $h/a$. Choose the extensional screening length of a slumping column and the associated scales asHere $h_0=R-a_0$ denotes the initial annular gap in this question. With $Z=z/\hat{z}$, $T=t/\hat{t}$ and $W=w/\hat{w}$, the dimensional equations reduce to

Initially $\varphi ={\varphi }_0$ and ${\Lambda }_0=L_0/\hat{z}$, so $W_{ZZ}=1+W$, with $W(0)=0$ and $W_Z({\Lambda }_0)=0$. The initial velocity profile of an annularly confined column isThe column thickens everywhere below its stress-free top, while that initial thickening rate is zero at the top.

For a tall column, ${\Lambda }_0\gg 1$, the limiting forms away from exponentially small end corrections areThe bulk falls at approximately $-\hat{w}$, where excess weight balances annular hydraulic resistance. Extensional stresses adjust this velocity to zero in a basal boundary layer of thickness $\hat{z}$; almost all thickening occurs there.

For a short column on this scale, ${\Lambda }_0\ll 1$, while still slender, the limits areThe annular-resistance term $W$ is small, so excess weight is balanced predominantly by extensional viscous stress. The top speed is $-\hat{w}{\Lambda }_0^2/2$, recovering the zero-annular-viscosity result. Thus $\hat{z}$ is the vertical distance over which extensional stress communicates the basal constraint before annular resistance screens it. The sketches below use separate natural normalizations for the two limits.

Take $\alpha >0$; a real disturbance with negative wavenumber is recovered by complex conjugation. To obtain the pressure equation for a shear-flow normal mode, apply $i\alpha$ to the streamwise momentum equation and differentiate the normal momentum equation. The mass conservation equation cancels the viscous divergence and the terms proportional to $U-c$, leaving $2i\alpha U^{\prime }v={\alpha }^{2}p-p^{\prime \prime }$. Thus the fluid pressure equation and wall value areThe wall value follows from streamwise momentum with no-slip boundary condition. Normal momentum also gives $p^{\prime }(0)=v^{\prime \prime }(0)/R$.

In the lower viscous boundary layer, mass conservation gives $v=O(\alpha \delta )$ if $u=O(1)$. Since $U\sim y$, convection and shear have size $\alpha \delta$ when $c=O(\delta )$; matching these with transverse viscous diffusion gives $R^{-1}{\delta }^{-2}\sim \alpha \delta$. The pressure gradient balance then gives $p=O(\delta )$. Thus the lower scalings areThe neglected streamwise viscous term is smaller by $(\alpha \delta {)}^2$. In these variables the leading equations areDifferentiating the second and using the first yields the Airy reduction $\boxed{u_{YYY}=i(Y-C)u_Y}$.

Put $q=e^{i\pi /6}$, so $q^3=i$. The bounded matching branch of the Airy function gives $u_Y=A\mathrm{Ai}(q(Y-C))$; the other independent Airy solution grows along the ray and is excluded. Integrating with the two wall conditions giveswhile the wall momentum equation gives the constant fluid pressureDefine the convergent Airy displacement integralsThe Airy function equation and its decaying derivative imply $I_1-C{I}_0=iq{\mathrm{Ai}}^{\prime }(-qC)$. Consequently, as $Y\to \infty$,Both differentiated lower momentum and the constant of integration are satisfied; retaining that constant is essential for the final matching.

For the middle region introduce a perturbation stream function $\Phi$ with $u={\Phi }^{\prime }$ and $v=-i\alpha \Phi$. The leading inviscid flow streamwise equation isFor a shear flow, this is the integrated form of the Rayleigh equation for inviscid shear flow. The leading matched solution for fixed $y=O(1)$ is $\Phi \sim BU$, so $u\sim B{U}^{\prime }$ and $v\sim -i\alpha BU$. In the overlap with the lower viscous boundary layer, the first constant correction iswhich agrees with the integrated lower solution. For a regular smooth profile, the normal momentum equation gives $p^{\prime }\sim -{\alpha }^2(U-c)\Phi$. Thus fluid pressure variation across a fixed-width middle region is $O({\alpha }^{2}B)$; in the final distinguished scaling this is smaller than its $O(\delta )$ value, so $p\sim \delta P_0$ throughout that region. Since $U\to 1$, the middle velocity tends to $\boxed{v\sim -i\alpha B}$, which does not decay and requires another region.

The fluid pressure equation becomes $p^{\prime \prime }-{\alpha }^{2}p=0$ when $U\sim 1$. Therefore the upper scale is $\Delta =1/\alpha$, with $z=\alpha y$. The inviscid decaying upper solution isMatching velocity as $z\downarrow 0$ gives $P_{\infty }\sim \alpha B$; matching fluid pressure to the middle region gives $\delta P_0\sim \alpha B$. This is three-layer long-wave shear-flow matching, expressed as a matched asymptotic expansion.

Under $\alpha =k{R}^{-1/4}$ with fixed positive $k$, $\delta =k^{-1/3}{R}^{-1/4}$, and $\alpha /\delta =k^{4/3}$. Eliminating the nonzero amplitude $A$ gives the Airy wall-layer matching determinant and dispersion relationEquivalently, with $\mathcal{T}(C)={\int }_{-qC}^{\infty e^{i\pi /6}}\mathrm{Ai}(\zeta )d\zeta =q{I}_0(C)$,The determinant form does not divide by a possibly vanishing $I_0$. The wavespeed is $c=\delta C$; temporal growth has sign $\mathrm{Im}C$ for positive $\alpha$.

There is a genuine conflict in the printed assumptions: fixed-$k$ $\alpha =k{R}^{-1/4}$ cannot satisfy $R^{-1/7}\ll \alpha$ as $R\to \infty$. For order-one lower matching amplitudes, constant leading middle fluid pressure requires ${\alpha }^2\ll \delta$, equivalently $\alpha \ll R^{-1/7}$; this is the reverse inequality and is satisfied by the final distinguished scaling. Thinness also requires $\alpha R\gg 1$. The dispersion relation above is therefore derived under the final distinguished limit, for which the complete consistent window is $R^{-1}\ll \alpha \ll R^{-1/7}$. Both printed limits cannot be imposed simultaneously.