Past exam of the mathematics course of the University of Cambridge / 2015 / iii / Paper 73 / 2

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.