Let $\delta =H/L\ll 1$. Incompressibility makes the transverse velocity $O(\delta U)$. In the axial Navier-Stokes equation, transverse viscous stress scales as $\mu U/H^2$, longitudinal viscous stress as $\mu U/L^2$, and convective inertia as $\rho U^2/L$. Thus lubrication theory requiresUnsteady inertia is also negligible if $\rho H^2/(\mu T)\ll 1$; this coincides with the preceding condition when $T=L/U$. The transverse pressure variation is smaller than this pressure scale by $O({\delta }^2)$, and viscous normal stress at the free surface is similarly smaller. Consequently pressure is uniform across the film to leading order, while axial viscous shear balances its axial gradient. These estimates justify the lubrication approximation with small surface slopes and a slowly varying rigid wall.

Write the physical cylindrical-film thickness as $d=ϵ\bar{h}$, so that $\bar{h}=O(a)$, and use $\bar{z},\bar{t}$ for dimensional coordinates. For a surface of radius $a+d$, its outward mean curvature, with the convention giving cylindrical curvature $1/a$, isThe Young–Laplace equation gives $p=p_{\mathrm{ext}}+\gamma \kappa$. With local normal coordinate $y$ measured from the rigid cylinder, the axial velocity solves $\mu u_{yy}=p_{\bar{z}}$, with no-slip boundary condition at $y=0$ and stress-free boundary condition at $y=d$. ThereforeThe factors of circumference in axisymmetric mass conservation cancel at leading thin-film order, giving $d_{\bar{t}}+q_{\bar{z}}=0$. IntroduceThis produces the capillary instability of a thin film on a rigid cylinder governed byThe curvature term $h_z$ destabilizes long waves; the axial-curvature term $h_{zzz}$ smooths short waves.

Linearize about $h_0$ using $h=h_0+\eta e^{st+ikz}$ with small $\eta$. The linear equation is $s+h_0^3(-k^2+k^4)=0$, soThe growth rate is positive for $0<|k|<1$, zero at $k=0,\pm 1$, and negative for $|k|>1$. Its derivative is $2h_0^{3}k(1-2k^2)$, givingThe requested sketch is even in $k$:

For the Rayleigh–Plateau instability of an entire liquid cylinder in the purely viscous limit with negligible exterior viscosity, the fastest mode tends to infinite wavelength. The rigid-wall coating instead selects a finite wavelength. The difference comes from the flow resistance: a coating must shear across its thickness, so its long-wave growth vanishes as $k^2$; a free liquid cylinder can undergo nearly plug-like axial extension. For example, the long-wave free-cylinder balances ${\eta }_t+(a/2)u_z=0$ and $3\mu u_{zz}=-\gamma ({\eta }_z/a^2+{\eta }_{zzz})$ give $s\sim \gamma (1-a^2{k}^2)/(6\mu a)$, which does not vanish as $k\to 0$. This comparison is for the creeping-flow, negligible-exterior-viscosity limit; finite inertia or exterior viscosity changes the wavelength selection. The full free-cylinder dispersion relation is discussed in Plateau-Rayleigh Instability of a Cylinder of Viscous Liquid.

The collar is a quasisteady sliding collar on a coated cylinder. It must slowly grow as well as translate: an exact travelling wave would obey $q-ch=\text{constant}$, whereas its stationary uniform far films have $q=0$ and would require $-c{h}_{+}=-c{h}_{-}$. Thus the stipulated unequal films require the stated slow evolution of the collar amplitude. In its $h=O(1)$ core, $c\ll 1$ and the volume flux is small, so the leading pressure gradient vanishes: $h^{\prime }+h^{\prime \prime \prime }=0$. Integration gives a constant plus sine and cosine terms. Choose the collar centre at its maximum to remove the sine term, and parameterize its edge height as $B$. ThenThe core is a local capillary equilibrium. In the moving frame, swept-up volume minus deposited volume gives the growth rate $c(h_{+}-h_{-})$. Its leading excess volume is $2\pi A$, and for constant $B$,Thus the printed rate is the rate of collar volume, or an order-of-magnitude amplitude rate; interpreted as an exact equation for $\dot{A}$, it lacks $1/(2\pi )$. This slow change is smaller than the translational core terms by the small far-film thickness and does not alter the leading local matching below.

Near either edge $x_f=\pm \pi$, let $h_f$ denote the adjacent film thickness, $h_{-}$ at the rear and $h_{+}$ at the front. Across the short transition region the collar's slow growth is higher order, so local mass conservation integrates to $q=c(h-h_f)$. The axial third derivative dominates the first derivative in the capillary flux, givingFor collar-to-film capillary matching, chooseThe two sides become $c{h}_f{H}^3{H}^{\prime \prime \prime }$ and $c{h}_f(H-1)$ respectively, and henceThe resulting transition length is small compared with the core length: the matching relations below imply $h_f=O(c^{2/3}/A)$, so $h_f/c^{1/3}=O(c^{1/3}/A)$ for an order-one core amplitude.

To count admissible branches, put $H=1+\zeta$ and linearize. The characteristic equation is $r^3=1$, with rootsAt the rear, decay as $X\to -\infty$ permits only $\zeta \sim D{e}^X$. For the nontrivial branch growing into a positive collar, $D>0$; translation rescales $D$, leaving one orbit. This branch has positive $H^{\prime },H^{\prime \prime },H^{\prime \prime \prime }$. It exists for all forward $X$, since $0<H^{\prime \prime \prime }=(H-1)/H^3\le 4/27$ for $H\ge 1$, and once $H^{\prime \prime }>0$, $H$ grows at least quadratically. Then $H^{\prime \prime \prime }=O(X^{-4})$, so successive integrations give a quadratic asymptote with a finite constant term; translation removes the linear term. The constant solution is not a rear-to-collar match. On the downward branch, $H<1$ gives $H^{\prime \prime \prime }<0$ while $H^{\prime }$ and $H^{\prime \prime }$ are initially negative; the slope remains negative and $H$ reaches zero, so it cannot match a positive collar. Uniqueness here refers to the positive, nontrivial rear matching branch, modulo translations.

At the front, decay as $X\to +\infty$ permits the two real amplitudes inTranslation removes one degree of freedom, leaving one parameter $\lambda$. These oscillatory branches account for the small capillary ripples ahead of the collar. This local stable manifold count gives the family of front branches; admissible global branches must also remain positive and reach a parabolic core asymptote. It does not, by itself, prove the global uniqueness of the matching parameter supplied in the question.

Write the rear and front core asymptotes, after translation removes their linear terms, as $a_{-}{X}^2+b_{-}$ and $a_{+}(\lambda )X^2+b_{+}(\lambda )$. Expanding the outer collar near either edge giveswhereas the inner expansion givesMatching the constant and quadratic coefficients therefore yieldsEliminate $h_{-}/h_{+}$ to obtain $a_{+}(\lambda )b_{+}(\lambda )=a_{-}{b}_{-}$. Using the specified selected root ${\lambda }_0$ and writing $\alpha =a_{+}({\lambda }_0)$, $\beta =b_{+}({\lambda }_0)$,The physical branch has $\alpha >a_{-}$, as required by $h_{-}<h_{+}$. No numerical values of these universal matching coefficients are needed for the requested relations.

Finally, the energy source is decreasing capillary surface energy, not Newtonian gravity. This can be checked through the volume-constrained capillary energy of a cylindrical coating. The exact liquid volume per axial length is $\pi [(a+d{)}^2-a^2]$, and the free-surface area is $2\pi \int (a+d)\sqrt{1+d_{\bar{z}}^2}d\bar{z}$. Expanding the latter and eliminating $\int d$ with fixed exact volume gives the nonconstant quadratic partApart from its positive dimensional factor, this is the surface energy $\mathcal{E}=\frac{1}{2}\int (h_z^2-h^2)dz$. Its variational derivative is $-h-h_{zz}$, and the thin-film equation is a conserved gradient flow. With periodic or zero-energy-flux boundaries, integration by parts givesFor an isolated collar between unequal infinite films, use the corresponding local energy balance, including far-film fluxes and subtracting background energies; the unsubtracted total energy is infinite. Collecting film into the growing collar lowers surface energy, which supplies the viscous dissipation sustaining its slow motion.

Let $({\mathbf{u}}^{(1)},{\mathbf{\sigma }}^{(1)})$ and $({\mathbf{u}}^{(2)},{\mathbf{\sigma }}^{(2)})$ be two Stokes flows in the same fluid domain, with equal dynamic viscosity and zero body force. The Lorentz reciprocal theorem for Stokes flow isTo prove it, use the symmetric Newtonian fluid stress tensor ${\mathbf{\sigma }}^{(k)}=-p^{(k)}\mathbf{I}+2\mu {\mathbf{e}}^{(k)}$ and ${\partial }_j{\sigma }_{ij}^{(k)}=0$. Incompressibility removes the pressure contraction, while symmetry removes the antisymmetric velocity gradient. ThereforeThe divergence theorem proves the boundary identity. It applies to an exterior domain by a limit at infinity; the boundary integral there vanishes for decaying rigid-body Stokes flows.

Use a normal directed outward from the fluid, into the rigid body. Then $\mathbf{F}=\int \mathbf{\sigma }\mathbf{n}dS$ and $\mathbf{G}=\int \mathbf{r}\times (\mathbf{\sigma }\mathbf{n})dS$ are the force and couple exerted by the body on the fluid. If $\mathbf{q}=(\mathbf{U},\mathbf{\Omega })$ and $(\mathbf{F},\mathbf{G})=\mathsf{R}\mathbf{q}$, inserting the rigid boundary velocity $\mathbf{U}+\mathbf{\Omega }\times \mathbf{r}$ into the Lorentz reciprocal theorem for Stokes flow givesArbitrary rigid velocities imply $\mathsf{R}={\mathsf{R}}^T$. A separate integration by parts, using the Stokes equation, givesEquality requires zero rate-of-strain tensor everywhere. Such a velocity is a rigid motion in the connected exterior fluid; decay at infinity forces it to vanish, and the no-slip boundary condition then forces $\mathbf{U}=\mathbf{\Omega }=0$. Hence the resistance matrix is symmetric positive definite. Defining force as fluid-on-body would reverse its sign, so the positive convention matters.

For the rods, choose body axes along $OA$, $OB$, and their normal. On the first rod, $\mathbf{X}=s{\mathbf{e}}_x$ for $0\le s\le 2L$, and the slender-body force density isOn the second rod, $\mathbf{X}=s{\mathbf{e}}_y$, andIntegrate the force densities and the moments $\mathbf{X}\times \mathbf{f}$, using $\int ds=2L$, $\int sds=2L^2$, and $\int s^{2}ds=8L^3/3$. For example, $F_x=C(3L{U}_x-2L^2{\Omega }_z)$, $G_x=C(2L^2{U}_z+8L^3{\Omega }_x/3)$, and $G_z=C[-2L^2{U}_x+2L^2{U}_y+16L^3{\Omega }_z/3]$. In the ordering $(U_x,U_y,U_z,{\Omega }_x,{\Omega }_y,{\Omega }_z)$, the full right-angle two-rod resistance matrix isThe planar subsystem isInverting this symmetric matrix gives the hydrodynamic mobility matrixThe upper-right entry is positive, as printed in the PDF; the negative sign in the converted TeX would violate both this inversion and reciprocity.

Let laboratory $X$ point along the initial $OA$ and laboratory $Y$ point upward. At body angle $\theta$, ${\mathbf{e}}_x=(\cos \theta ,\sin \theta )$ and ${\mathbf{e}}_y=(-\sin \theta ,\cos \theta )$. The gravitational load, balanced by body-on-fluid resistance, has body componentsThe last expression follows by taking moments of the two endpoint weights about $O$. Substitution into the hydrodynamic mobility matrix givesFor $\lambda =1$, the initial condition $\theta =0$ persists, and $U_x=0$, $U_y=-mg/(CL)$: vertical descent without rotation. For $\lambda <1$, $\dot{\theta }>0$ throughout $0\le \theta <\pi /4$, and the zero of $\cos \theta -\sin \theta$ at $\pi /4$ is stable. A simple zero is approached only as $t\to \infty$, so $\theta$ increases monotonically to $\pi /4$. For $\lambda >1$, it decreases monotonically through $-\pi /4$ and $-\pi /2$ toward the next stable zero: $\theta \to -3\pi /4$, not $-\pi /4$.

The laboratory horizontal velocity of $O$ isFor $\lambda \ne 1$, divide by the angular equation to eliminate time:At both limiting orientations the sine and cosine are equal, so the sedimentation drift of a weighted two-rod body isFor $\lambda >1$ the initial drift is negative. It reaches $X-X(0)=2L(1-\sqrt{2})/3$ at $\theta =-\pi /4$, reverses, crosses its starting horizontal position at $\theta =-\pi /2$, and ends $2L/3$ to the right. At exactly $\lambda =1$ there is no horizontal displacement; the infinite-time displacement is discontinuous there because the rotation time diverges as $\lambda \to 1$.