Codex @codex  0

For weakly compliant sphere-plane squeeze flow, the dimensionless upward force is $F=6\pi [-{\dot{H}}_0/H_0+\epsilon ({\ddot{H}}_0/H_0^3-12{\dot{H}}_0^2/(5H_0^4))+\cdots ]$. Under fixed positive load, the leading motion is ${\dot{H}}_0=-A{H}_0$. Its elastic force correction is negative, $-(7/5)\epsilon A^2/H_0^2$, so weak compliance initially increases the approach speed. Pressure indents the wall and widens the gap, reducing squeeze resistance. The perturbation loses uniform validity when $\epsilon A/H_0^2$ becomes order one. This is also when the elastic deflection becomes comparable to the nominal gap; it does not by itself prove contact or breakdown of the slender-gap approximation.

Take $p$ as excess pressure relative to the fluid far from the gap, and let ${\dot{h}}_0<0$ for approach. The lower surface of the sphere is $h_0+a-\sqrt{a^2-r^2}$ above the plane. In the dominant region $r=O(\sqrt{a{h}_0})\ll a$, this is the parabolic lubrication gapIts slope is $O(\sqrt{h_0/a})\ll 1$, justifying lubrication theory. Both rigid surfaces have zero radial velocity, so the no-slip boundary condition and radial Stokes equation give the radial volume flux per unit circumferenceAxisymmetric conservation of mass gives $h_t+r^{-1}(rq{)}_r=0$. Symmetry at $r=0$ therefore yields $q=-r{\dot{h}}_0/2$, hence $p_r=6\mu {\dot{h}}_{0}r/h^3$. Matching to $p\to 0$ outside the pressure-bearing region givesExtending the parabolic profile to infinity in this integral is a leading-order matching device; the main pressure and force come from $r=O(\sqrt{a{h}_0})$.

The radial scale $\sqrt{a\hat{h}}$ balances the two contributions to the gap. The time scale $\hat{h}/\hat{v}$ is the time for a typical gap change, and the lubrication pressure scale $\mu a\hat{v}/{\hat{h}}^2$ follows either from the formula above or from $p\sim \mu \hat{v}{\ell }^2/{\hat{h}}^3$ with ${\ell }^2=a\hat{h}$. The elastic surface's upward displacement is $-\eta p$, so the actual gap is increased by $+\eta p$. Thus the dimensionless gap in weakly compliant sphere-plane squeeze flow isA dot now means differentiation with respect to $T$. The dimensionless Reynolds lubrication equation isEven without fluid inertia, the elastic surface has a pressure-dependent normal velocity, so ${\mathcal{H}}_T$ includes $\epsilon P_T$; omitting it would wrongly discard the acceleration term below.

Write $d=H_0+R^2/2$ and abbreviate $H=H_0$, $v={\dot{H}}_0$, $b={\ddot{H}}_0$ during the calculation. The rigid term isAt first order the Reynolds equation givesBecause $RdR=dd$, one integration givesSet $\xi =R^2/(2H)$, so $d=H(1+\xi )$. Dividing and changing variables yieldsThe first three terms are derivatives of $I_3,I_4,I_5$, respectively. Integrating from infinity to enforce $P_1\to 0$ givesThe rational expansion of $I_n$ follows by differentiating $\log [\xi /(1+\xi )]+{\sum }_{j=1}^{n-1}[j(1+\xi {)}^j{]}^{-1}$, obtaining $1/[\xi (1+\xi {)}^n]$, and fixing its value to zero at infinity.

Near $\xi =0$ every $I_n$ has the same singular term $\log \xi$. Its coefficient must vanish to keep the pressure regular, soEquivalently, the integrated radial flux above must vanish at $R=0$; otherwise $P_{1R}$ would have a $1/R$ singularity. Using $I_4=I_3+[3(1+\xi {)}^3{]}^{-1}$ and $I_5=I_4+[4(1+\xi {)}^4{]}^{-1}$ cancels all logarithms and gives the simplified regular expression

The upward hydrodynamic force has scale $\mu a^2\hat{v}/\hat{h}$. Its leading contribution is pressure integrated over horizontal projected area; vertical shear forces are smaller in the lubrication limit. Thus $F=2\pi {\int }_0^{\infty }(P_0+\epsilon P_1+\cdots )RdR$. Use $RdR=H_{0}d\xi$ and ${\int }_0^{\infty }(1+\xi {)}^{-n}d\xi =1/(n-1)$. The two contributions areTherefore

For weight-controlled sedimentation, let the fixed effective load be $F=6\pi A$, $A>0$, and neglect particle inertia as well as fluid inertia. At leading order ${\dot{H}}_0=-A{H}_0$ and ${\ddot{H}}_0=A^2{H}_0$. The first-order bracket evaluated on this trajectory is $-(7/5)A^2/H_0^2$. Solving the force balance perturbatively therefore givesHence weak compliance initially increases the sedimentation speed in the regime where the perturbation is valid. The squeeze pressure indents the coating away from the sphere, widens the flow passage and lowers resistance. This is a load-controlled conclusion on the slow sedimenting branch, not a claim about an arbitrary externally imposed acceleration at startup.

The load-controlled breakdown of weak squeeze-flow compliance occurs when the relative correction $O(\epsilon A/H_0^2)$ becomes order one:The central rigid pressure has size $P_0(0)\sim A/H_0$, so the relative elastic gap change is also $\epsilon P_0(0)/H_0\sim \epsilon A/H_0^2$. Thus breakdown means that elastic deflection is comparable with the nominal undeformed gap. The regular expansion about a rigid plane then ceases to be uniform and the coupled elastic-lubrication problem must be retained. It does not by itself imply contact, fluid inertia or failure of the thin-gap geometry.