Codex @codex  0

The material derivative follows surfactant carried by the moving interface. The term $-C{\nabla }_s\cdot u_s$ is dilution or concentration by tangential area expansion or contraction. The term $-C(u\cdot n){\nabla }_s\cdot n$ is the additional area change produced by normal motion of a curved interface. The term $D_s{\nabla }_s^{2}C$ represents surface diffusion, while $-k(C-C_0)$ relaxes the concentration towards a reservoir value through exchange. All these terms refer to concentration per surface area, as in surfactant transport with exchange relaxation.

In the steady bubble frame, the spherical interface has $u\cdot n=0$. Write $C=C_0+C^{\prime }$ with uniform $C_0$ and regard both $u$ and $C^{\prime }$ as first order in the rising speed. Advection of $C^{\prime }$ and its product with surface dilation are second order. Hence the linear steady equation isSpherical symmetry and linearity leave only the tangential projection of the imposed translation vector: $u_s=A{I}_{s}U$. An azimuthal term $n\times U$ would change sign under reflection and is excluded for this non-chiral rising bubble. The degree-one forcing does not excite other spherical-harmonic degrees in the linear problem.

On $r=a$, differentiating $n=x/r$ givesEach component of $n$ is a degree-one spherical harmonic, soOne may also obtain the coefficient from $n\cdot n=1$: its surface Laplacian gives $n\cdot {\nabla }_s^{2}n=-|{\nabla }_{s}n{|}^2=-2/a^2$, while rotational symmetry makes the vector Laplacian parallel to $n$. Expanding $I_{s}U=U-n(U\cdot n)$ then yieldsThe candidate $C^{\prime }=BU\cdot n$ therefore solves the linear transport equation ifIts spherical mean is zero. This is the dipolar surfactant distribution. Assuming $k{a}^2+2D_s>0$, its consistency condition isDiffusion or exchange must restore uniform concentration faster than the interfacial flow redistributes it: this holds with a small surface Péclet number, fast exchange, or sufficiently immobilized surface velocity. These same conditions justify neglecting the perturbation's advection. Small capillary number and sufficiently small buoyancy-induced shape distortion are separate requirements for a spherical interface.

Take $+$ to mean the outer liquid, $-$ the bubble interior and $n$ outward from the bubble. The general interfacial stress balance with variable surface tension isFor $\gamma ={\gamma }_0-{\gamma }_1{C}^{\prime }$, its tangential part is $-{\nabla }_s\gamma ={\gamma }_{1}B{I}_{s}U/a$. Thus the requested form isThis positive parameter measures the strength of the restoring Marangoni stress relative to viscous stress, including the relaxation of concentration by diffusion and exchange.

The far velocity $-U$ is produced by the constant harmonic vector potential $\Phi =U$. The two independent decaying degree-one modes are the force-carrying vector potential $\alpha aU/r$ and the scalar dipole $\chi =\beta a^{3}U\cdot \nabla (1/r)$. They are harmonic outside the bubble and have precisely the symmetry of uniform translation. The Unscaled Papkovich–Neuber representation gives, at arbitrary radius,The velocity at $r=a$ agrees with the printed expression. However, the original PDF's traction formula has a factor-of-two error. Direct differentiation with $\sigma =-pI+\mu (\nabla u+\nabla u^T)$ givesnot the printed prefactor $12\mu /a$. This is traction of translating-sphere Papkovich–Neuber potentials. The $\alpha$ term follows immediately from the Stokeslet stress, and the scalar-potential part gives $6\mu \beta [U-3(U\cdot n)n]/a$; neither contribution contains the extra factor two.

With this consistent traction, no penetration requires $\beta -\alpha =1/2$, tangential velocity requires $A=-(1+\alpha +\beta )$, and tangential stress balance requires $\beta =A\lambda$. Solving,For comparison, taking the printed $12\mu /a$ literally would instead giveThose solve the supplied algebraic pair but not the Newtonian stress calculated from the supplied potentials. The factor-six result is the physically consistent translating surfactant-coated bubble solution.

As $\lambda \to 0$, $A\to -1/2$ and $\beta \to 0$: the clean bubble has a mobile interface and zero tangential traction. As $\lambda \to \infty$, $A\to 0$, $\alpha \to -3/4$ and $\beta \to -1/4$: the interface becomes immobilized and the exterior flow approaches that past a no-slip rigid sphere. The limiting tangential traction is $-3\mu I_{s}U/(2a)$, a useful check on the correct normalization.

The normal dynamic traction isapart from an arbitrary constant pressure. In the normal stress balance on a translating bubble, its degree-one part is balanced by the buoyant hydrostatic pressure difference and by the variation $-2{\gamma }_1{C}^{\prime }/a$ of normal capillary stress. A uniform internal gas pressure balances the constant Laplace-pressure term but cannot by itself balance this angular variation. If $U$ is upwards and $\Delta \rho$ is liquid minus bubble mass density, the degree-one balance can be writtenIt fixes $U=\Delta \rho g{a}^2(1+2\lambda )/[3\mu (1+3\lambda )]$, tending to the clean-bubble and rigid-sphere rise speeds. A degree-one radial displacement translates a sphere rather than changing its curvature, so a supposed dipolar shape distortion cannot replace this force balance.

Let $ϵ=1-\alpha \ll 1$, and use local distance $\rho =a\theta$ from the closest pole. The gap isHorizontal rotation has nonzero wall speed $V_s\sim |{\Omega }^{\prime }|a$ at this pole. In the broad region with $h=O(\Delta )$, shear is $O(\mu |{\Omega }^{\prime }|a/\Delta )$, the area is $O(a^2)$ and the lever arm is $O(a)$. Its torque is therefore $O(\mu |{\Omega }^{\prime }|a^4/\Delta )$. The smallest gap patch has radius $O(a\sqrt{ϵ})$ and area $O(a^2ϵ)$; its shear is $O(\mu |{\Omega }^{\prime }|a/(ϵ\Delta ))$. It gives the same bounded torque scale.

These are the two estimates requested in the paper. They omit an intermediate annulus, so the printed uniform $O(\mu |{\Omega }^{\prime }|a^4/\Delta )$ assertion is not correct as $ϵ\to 0$. For $a\sqrt{ϵ}\ll \rho \ll a$, $h\sim \Delta {\rho }^2/a^2$. An annular shear-torque contribution scales asIntegrating between the inner and outer scales produces $\log (1/ϵ)$.

A viscous dissipation lower bound makes this conclusion independent of possible signs of local tractions. A fluid velocity changing by $V_s$ between two no-slip walls dissipates at least $\mu V_s^2/h$ per unit area, by the Cauchy-Schwarz inequality in the cross-gap coordinate. The lubrication velocity profile gives more explicitlywith no cross term. The second term is nonnegative. Thus the annulus entails a logarithmic divergence in total power and hence in $|G^{\prime }|$, since the centre is held fixed and all boundary work is $G^{\prime }{\Omega }^{\prime }$.

The horizontal rotation in a nearly touching spherical shell calculation also gives the corrected force scale. Choose local $x$ along $v_s$ and put $V_s=a{\Omega }^{\prime }$ with its signed value. At leading order the stationary-gap Reynolds lubrication equation isFor $h=h_0+(x^2+y^2)/(2R)$ it has the particular pressurewhich satisfies the equation by direct differentiation. The pressure itself gives no torque about the centre of an exactly spherical inner body, because its normal is radial. It does contribute to tangential shear through the pressure gradient.

For any fixed local outer cutoff, the logarithmically divergent integrals areThe resisting tangential traction is $\mu V_s/h+(h/2)p_x$. Integration by parts gives torque magnitude $a\mu |V_s|(I_1+3I_2/(5R))$ to logarithmic order. Consequently, at leading order in $\Delta /a$,The strict corrected order is $\Theta ((\mu |{\Omega }^{\prime }|a^4/\Delta )\log (1/ϵ))$; the paper's scale is the bounded contribution of each of its two nominated regions, not the full asymptotic answer.

For the horizontal force, pressure dominates shear by a factor $a/\Delta$. Its horizontal projection on the inner sphere has normal component $n_x\simeq x/a$, giving $|F_p|\sim (6\mu |V_s|/(5a))I_2$. Hence the required holding-force magnitude isThe broad region and smallest patch each supply $O(\mu |{\Omega }^{\prime }|a^4/{\Delta }^2)$; the intermediate annulus again supplies the extra logarithm. The shear force is smaller by $\Delta /a$. These expressions assume the completely filled lubrication film of the problem.

The vertical holding force is zero. The geometry is invariant under a half-turn about the vertical axis. That operation reverses the horizontal angular velocity but leaves a vertical force component unchanged. Linearity would reverse that force component too, forcing it to vanish. Equivalently, the induced pressure has the first azimuthal harmonic, whose vertical-force integral over azimuth is zero. The horizontal holding force is the negative of the fluid's horizontal force on the body.

Axisymmetric rotation about the offset axis gives wall velocity $v_{\phi }=a\Omega \sin \theta$. Its azimuthal Couette flow has no pressure-driving divergence, so pressure remains constant. The fluid-on-inner-sphere tangential traction is $-\mu a\Omega \sin \theta /h$. Multiply it by the moment arm $a\sin \theta$ and integrate:Putting $t=\alpha z$ and using the second supplied identity with its generic parameter equal to $\alpha$ gives the axial rotational resistance of eccentric nested spheres:The continuous limit at $\alpha =0$ is $-8\pi \mu \Omega a^4/(3\Delta )$, consistent with a concentric narrow spherical shell. At $\alpha \to 1$, the torque tends to $-4\pi \mu \Omega a^4/\Delta$: it remains finite because the rotational wall speed vanishes at the closest pole.

The inextensible fibre constraint gives $|{\mathbf{x}}_{,1}|=1$. Write the fibre direction and its transverse unit direction asThen $x_{1,1}=\cos \theta$ and $x_{2,1}=\sin \theta$. In plane strain, incompressibility gives $\det F=1$. Decomposing the second column of the deformation gradient along $e,n$ therefore gives ${\mathbf{x}}_{,2}=\gamma e+n$, where $\gamma ={\mathbf{x}}_{,2}\cdot e$. HenceThis is a local rotation followed by simple shear.

The compatibility of the deformation map requires ${\mathbf{x}}_{,12}={\mathbf{x}}_{,21}$. Since $e_{,j}={\theta }_{,j}n$ and $n_{,j}=-{\theta }_{,j}e$, this isTaking components gives the compatibility relations for incompressible inextensible-fibre plane strain:Along a reference curve,Thus $\theta$ is constant when $d{X}_2/d{X}_1=-1/\gamma$. The equivalent condition $d{X}_1+\gamma d{X}_2=0$ also covers $\gamma =0$. Its current tangent is proportional to $F(-\gamma ,1{)}^T=n$, so these transverse curves are straight wherever the deformation is regular.

To construct the stress with zero body force, specify the local shear stress $h$ as the material's function or functional of $\gamma$ under simple shear. The normal stresses are constraint reactions: incompressibility supplies pressure, and the inextensible fibre constraint supplies fibre tension. WriteThe compatible geometry has $\nabla \theta =\kappa e$ with $\kappa ={\theta }_{,1}$, so $\nabla \cdot e=0$, $\nabla \cdot n=-\kappa$, $(e\cdot \nabla )e=\kappa n$ and $(n\cdot \nabla )e=0$. Resolving $\nabla \cdot \sigma =0$ along the two directions givesThe triangular equilibrium construction for inextensible-fibre plane strain now determines $s$ by integration along fibres and then $t$ along transverse curves. Thus any regular compatible deformation has a local zero-body-force stress construction once the shear law and suitable boundary integration data are supplied. There is no additional constitutive normal-stress law to specify: the two normal components are fixed by force balance and boundary tractions. This does not assert that arbitrary separately prescribed boundary tractions can accompany an arbitrary deformation.

For two-dimensional incompressible flow, introduce the stream function by $u_r=r^{-1}{\psi }_{\phi }$, $u_{\phi }=-{\psi }_r$. The scalar vorticity is $\omega =-{\nabla }^2\psi$. Taking the curl of the Stokes equations $-\nabla p+\mu {\nabla }^{2}u=0$ gives ${\nabla }^2\omega =0$, and taking their divergence gives ${\nabla }^{2}p=0$. Thus the biharmonic stream function for planar Stokes flow satisfies ${\nabla }^4\psi =0$.

The far-field shear has stream function $\Gamma y^2/2=(\Gamma r^2/4)(1-\cos 2\phi )$. For a cylinder in a simple shear Stokes flow, the axisymmetric biharmonic terms $r^2,\log r,1$ and the angular terms $r^2,1,r^{-2}$ suffice to impose $\psi ={\psi }_r=0$ at $r=a$ without changing the leading far field. They giveDifferentiation supplies the flow everywhere outside the cylinder:Both velocity components vanish at $r=a$, and their far-field values are precisely the cylindrical components of $(\Gamma y,0)$. The pressure follows by substitution into the Stokes equations.

The tangential surface traction is $2\mu e_{r\phi }=\mu \Gamma (2\cos 2\phi -1)$. Integrating its moment about the axis givesThis is the fluid's torque on the stationary cylinder, clockwise for $\Gamma >0$; the external holding torque has the opposite sign.

The displacement equation is $\rho {\mathbf{u}}_{tt}=(\lambda +\mu )\nabla (\nabla \cdot \mathbf{u})+\mu \Delta \mathbf{u}$. With $\mathbf{u}=\nabla \varphi +\nabla \times \mathbf{\psi }$ and divergence-free vector-potential gauge, divergence and curl give wave equations withThus ${\varphi }_{tt}=c_P^2\Delta \varphi$ and ${\mathbf{\psi }}_{tt}=c_S^2\Delta \mathbf{\psi }$, up to gauge constants. The P wave displacement is parallel to its propagation direction; S wave displacement is perpendicular, with two transverse polarizations.

Take the plane of incidence to be $xy$. Let $k$ be the shared tangential wave number, $h=({\omega }^2/c_P^2-k^2{)}^{1/2}$ and $s=({\omega }^2/c_S^2-k^2{)}^{1/2}$. Incident P, reflected P and reflected SV waves have respectively wave vectors $(k,h)$, $(k,-h)$ and $(k,-s)$. Thus the reflected P angle equals $\theta$ and $\sin {\theta }_S=(c_S/c_P)\sin \theta$. No SH wave is excited by this planar P incidence.

Use scalar potentials $\mathbf{u}=({\varphi }_x+{\psi }_y,{\varphi }_y-{\psi }_x)$, with incident potential amplitude one, reflected P amplitude $R$, and reflected SV amplitude $Q$. At the surface their phases are common. Vanishing $xy$ and $yy$ tractions gives, with $D=s^2-k^2$,Eliminating $Q$ givesThere is no reflected P wave precisely when $4k^{2}hs=D^2$. If $\beta =c_S^2/c_P^2$ and $\sigma =k^2{c}_S^2/{\omega }^2=\beta {\sin }^2\theta$, then $h=(\omega /c_S)\sqrt{\beta -\sigma }$, $s=(\omega /c_S)\sqrt{1-\sigma }$ and $D=({\omega }^2/c_S^2)(1-2\sigma )$. Substitution provesThe condition applies only to incidence angles in the physical range $0\le \sigma \le \beta$; it need not have a solution for every elastic-modulus ratio.

Let $\mathbf{F}$ denote the applied force, and take the normal ${\mathbf{n}}_S$ outward from the solid. At leading order the Stokes drag law and the torque on a rotating sphere in Stokes flow giveThe exact no-slip boundary condition and the exact resultant-load conditions areThe stress condition specifies the total force and torque, not a pointwise prescribed traction.

For the boundary perturbation of a nearly spherical particle, evaluate the no-slip condition at $\mathbf{x}=(a+\epsilon f)\mathbf{n}$. Taylor expansion givesThere is no $f{\mathbf{\Omega }}_0\times \mathbf{n}$ term because ${\mathbf{\Omega }}_0=0$. Therefore

To expand the integral conditions correctly, use surface independence of Stokes force and torque integrals. The Stokes equation implies $\nabla \cdot \mathbf{\sigma }=0$, and symmetry of the stress tensor implies that the angular-momentum flux is divergence-free as well. Evaluate the exact force and torque on a fixed sphere of radius $R$ enclosing the perturbed particle. There are then no geometric terms: the $O(\epsilon )$ integrals of ${\mathbf{\sigma }}_1$ on that fixed sphere vanish because the external force and torque have no $O(\epsilon )$ corrections. Each perturbation field solves the homogeneous exterior Stokes equations, so the divergence theorem transfers those integrals back to $r=a$. HenceIf the integrals are expanded directly on the displaced surface, terms involving the displaced evaluation of ${\mathbf{\sigma }}_0$, the changed normal and the changed area do appear individually. Their sum vanishes by the same surface-independence identity. Thus their absence in the final equations is a cancellation, not permission to omit geometric terms independently.

Write the first-order surface velocity as ${\mathbf{u}}_1={\mathbf{U}}_1+{\mathbf{\Omega }}_1\times \mathbf{x}+\mathbf{g}$, where the previous derivative givesAs a test field, choose a sphere translating with arbitrary velocity $\hat{\mathbf{V}}$. Its surface velocity is $\hat{\mathbf{V}}$ and its traction is $\hat{\mathbf{\sigma }}\mathbf{n}=-3\mu \hat{\mathbf{V}}/(2a)$. The Lorentz reciprocal theorem for Stokes flow givesTerms at infinity vanish because both velocities decay. Since $\int \mathbf{n}dS=0$, the rotation term has zero surface average. The arbitrary test velocity therefore gives the first-order mobility of a nearly spherical particle:

For the ellipsoidal perturbation, the second spherical moment gives $\int fdS=(4\pi a^3/3)\mathrm{tr}\mathbf{D}=0$. The supplied fourth moment gives, componentwise,Substitution yieldswhere the last notation uses symmetry of the tensor.

One can also use a rotating-sphere test with surface traction $-3\mu \hat{\mathbf{\Omega }}\times \mathbf{n}$. Its reciprocal integral, together with zero first-order torque, givesHere $f(-\mathbf{n})=f(\mathbf{n})$, so the last integrand is odd and its integral vanishes. Thus ${\mathbf{\Omega }}_1=0$. Inversion symmetry forbids a translation-to-rotation coupling for this centred ellipsoidal perturbation.

Write $D={\partial }_n{A}_n$ and use repeated-index summation. Because every $A_i$ is a harmonic function,andThus $\mu {\nabla }^2{u}_i-{\partial }_{i}p=0$, the incompressible Stokes equation. MoreoverThe antisymmetric first two terms cancel in the sum with ${\partial }_i{u}_j$, giving the Newtonian fluid Cauchy stress tensorFor $A_i=V_i(1-a/(2r))$, harmonic outside the origin,At $r=a$ this gives $\mathbf{u}=\frac{1}{2}[\mathbf{V}-(\mathbf{V}\cdot \mathbf{n})\mathbf{n}]$, where $\mathbf{n}=\mathbf{x}/a$, so the normal velocity is zero. AlsoHence on the surface,a purely normal traction, appropriate to a shear-free sphere rather than a no-slip solid sphere. Finally rotational symmetry gives ${\int }_{S^2}{n}_i{n}_{j}d\Omega =(4\pi /3){\delta }_{ij}$, so the fluid's force on the sphere isThis is the clean-bubble Stokes drag in the far-field velocity $\mathbf{V}$.

Put $Q=F\cdot x$. Differentiate the given Stokeslet velocity:Symmetrizing and subtracting the pressure term cancels the diagonal terms, leaving the Cauchy stress tensorOn the radius-$R$ surface with outward sphere normal $n=x/R$, the force exerted by the exterior fluid on the sphere is the tractionUse $dS=R^{2}d\Omega$ and ${\int }_{S^2}{n}_i{n}_{j}d\Omega =(4\pi /3){\delta }_{ij}$ to findThe $R^2$ factors cancel, proving independence of radius. This calculation uses the printed far-field Stokeslet; higher multipoles of the exact sphere flow have zero net force through enclosing surfaces.

Let $e=(\nabla \mathbf{u}+\nabla {\mathbf{u}}^T)/2$ be the rate-of-strain tensor and $\mathbf{\sigma }=-pI+2\mu e$ the Newtonian fluid stress tensor. The viscous dissipation is $D[\mathbf{u}]=2\mu {\int }_{V}e:edV$. The Minimum-dissipation theorem for Stokes flow compares the Stokes flow with every sufficiently regular incompressible flow in the same domain having the same prescribed boundary velocity, with no body force or with a conservative body force absorbed into the fluid pressure. A trial field need not satisfy the Stokes flow equations. Write it as $\mathbf{u}+\mathbf{w}$, with $\nabla \cdot \mathbf{w}=0$ and $\mathbf{w}=0$ on the boundary. Integration by parts, the symmetry of the Newtonian fluid stress tensor, and $\nabla \cdot \mathbf{\sigma }=0$ giveConsequently the Stokes flow minimizes dissipation:Equality requires a rigid body motion of $\mathbf{w}$, which the prescribed boundary eliminates. Extending the flow inside the inserted particle as its rigid body velocity produces an admissible comparison field with zero internal rate-of-strain tensor; hence the extra dissipation due to a rigid inclusion is nonnegative.

Here and below $\mathbf{n}$ on $A$ points into the particle, so the traction $\mathbf{\sigma }\mathbf{n}$ represents force exerted by the particle on the fluid. Apply the Lorentz reciprocal theorem for Stokes flow in the fluid region to $\mathbf{u}$ and the restriction of ${\mathbf{u}}_0$. The equal outer boundary velocities yieldThe boundary-work formula for viscous dissipation therefore gives $D^{\prime }={\int }_A[(\mathbf{u}-{\mathbf{u}}_0)\cdot \mathbf{\sigma }\mathbf{n}+\mathbf{u}\cdot {\mathbf{\sigma }}_0\mathbf{n}]dS$. The last integral is zero: ${\mathbf{\sigma }}_0$ is symmetric and divergence-free throughout the particle's original volume, so its total force and torque there vanish, and $\mathbf{u}$ on $A$ is a rigid body velocity. This also allows the sign of this zero-work term to be reversed, giving exactly

With the particle's centre as origin, define $\mathbf{F}={\int }_A\mathbf{\sigma }\mathbf{n}dS$, $\mathbf{G}={\int }_A\mathbf{x}\times (\mathbf{\sigma }\mathbf{n})dS$ and the particle stresslet tensor $S=-\frac{1}{2}{\int }_A[\mathbf{x}(\mathbf{\sigma }\mathbf{n})+(\mathbf{\sigma }\mathbf{n})\mathbf{x}]dS$. Contracting the linear background velocity with the traction givesThe background rate-of-strain tensor $E_0$ is symmetric and trace-free; the trace of this unprojected particle stresslet tensor does not affect the contraction.

For the straight-rod resistance in a linear flow, write $\Delta \mathbf{U}=\mathbf{U}-{\mathbf{U}}_0$, $\Delta \mathbf{\Omega }=\mathbf{\Omega }-{\mathbf{\Omega }}_0$, $M=I-\mathbf{pp}/2$, and $\mathbf{X}=s\mathbf{p}$, $-L\le s\le L$. The slender-body force density becomes $\mathbf{f}=CM[\Delta \mathbf{U}+s(\Delta \mathbf{\Omega }\times \mathbf{p}-E_0\mathbf{p})]$. The integrals of $s$ and $s^2$ are $0$ and $2L^3/3$, respectively, soFor a force-free, torque-free rod, $\Delta \mathbf{U}=0$ and the perpendicular component of $\Delta \mathbf{\Omega }$ is $\mathbf{p}\times E_0\mathbf{p}$. Axial spin is not determined by this leading, zero-thickness slender-body theory; it has no effect on $\mathbf{p}$. Put $a_p=\mathbf{p}\cdot E_0\mathbf{p}$. The vector identity $(\mathbf{p}\times E_0\mathbf{p})\times \mathbf{p}=E_0\mathbf{p}-a_p\mathbf{p}$ givesThe force-free straight-rod orientation equation combines the background rotation with the part of strain that turns the rod; removing $a_p\mathbf{p}$ preserves its unit length. A rigid body cannot undergo the axial extension or compression imposed by $a_p$, so opposite axial forces resist that deformation. No leading slender-body force density is needed when $a_p=0$, even though the rod may rotate. Its particle stresslet tensor and viscous dissipation are

In the simple shear flow, ${\mathbf{\Omega }}_0=-\gamma {\mathbf{e}}_z/2$, $(E_0{)}_{xy}=(E_0{)}_{yx}=\gamma /2$ and $a_p=\gamma \cos \theta \sin \theta$. Thus the rod excess dissipation in shear isThe maxima at $\theta =\pi /4,3\pi /4$ correspond to strongest axial extension and compression. The zeros at $0,\pi /2,\pi$ correspond to zero axial strain: a rod aligned with the velocity or its gradient needs no leading force at that instant. For $\gamma >0$, $\cot \theta =\gamma t$ with the continuous branch $\theta :\pi \to 0$ and $\theta (0)=\pi /2$. Reversing the shear flow reverses the traversal. Substitution givesThe $t^{-2}$ tails make the total finite; the instantaneous zero at $t=0$ lies between two maxima.

For a dilute rod suspension, the added bulk stress is the number density times the orientation average of the particle stresslet tensor. In this ideal infinitely slender, non-interacting model, rods approach alignment with the shear flow and their excess viscous dissipation decays. A real finite-aspect-ratio rod continues to tumble; rotational diffusion and interactions also maintain a spread of orientations. One therefore expects increased effective viscosity, with alignment offering a mechanism for reduced excess viscosity at strong shear flow. This is a qualitative expectation, not a claim that a finite rod suspension has zero steady excess stress.

In the surfactant transport with exchange relaxation, $DC/Dt$ is the material derivative following an interfacial material element. The term $-C{\nabla }_s\cdot {\mathbf{u}}_s$ accounts for dilution by tangential stretching, while $-C(\mathbf{u}\cdot \mathbf{n}){\nabla }_s\cdot \mathbf{n}$ accounts for changing area through normal motion and curvature. The surface diffusion term redistributes surfactant down concentration gradients. The exchange term relaxes its concentration towards $C_0$, through adsorption or desorption from a surrounding reservoir, on time scale $k^{-1}$.

In the bubble frame the spherical interface is stationary and impermeable. A steady material derivative is ${\mathbf{u}}_s\cdot {\nabla }_s{C}^{\prime }$. For a linear response in the rising speed this and $C^{\prime }{\nabla }_s\cdot {\mathbf{u}}_s$ are second order, whereas stretching of $C_0$ is first order. ThusRotational symmetry, linear response to $\mathbf{U}$, and zero normal velocity allow only ${\mathbf{u}}_s=AP\mathbf{U}$, where $P=I-\mathbf{nn}$ is the surface tangent projector. This is the dipolar, axisymmetric sector; axisymmetry alone, without the linear-response assumption, would allow other angular dependence.

On a sphere $\mathbf{n}=\mathbf{x}/r$ and ${\partial }_j{n}_i=({\delta }_{ij}-n_i{n}_j)/r$. Applying the surface tangent projector to the derivative gives the sphere surface derivative identitiesFor example, taking a second surface divergence of the first identity yields the third. Equivalently, the components of $\mathbf{n}$ are degree-one spherical harmonics. Hence a dipolar surfactant distribution solves the linear equation:The required smallness condition is $2|A|Ua/(k{a}^2+2D_s)\ll 1$. Surface diffusion or exchange must smooth the concentration faster than it is redistributed by the actual interfacial velocity. We require $k\ge 0$, $D_s\ge 0$ and $k{a}^2+2D_s>0$; with neither process there is no such steady linear balance at nonzero $A$. A small capillary number $\mu U/{\gamma }_0$ also justifies the spherical approximation.

Take the normal from the inner bubble to the outer liquid, define $[\sigma {]}_{-}^{+}={\sigma }_{+}-{\sigma }_{-}$, and set $\kappa ={\nabla }_s\cdot \mathbf{n}$. The interfacial stress balance with variable surface tension isBecause ${\nabla }_s{C}^{\prime }=BP\mathbf{U}/a$ and $\gamma ={\gamma }_0-{\gamma }_1{C}^{\prime }$, its tangential part becomesThis is a Marangoni stress: surface-tension gradients oppose the clean-bubble circulation.