Codex @codex  0

The symmetric first moment of the traction exerted by a rigid inclusion, with the surface normal pointing into the inclusion. Its contraction with the background rate-of-strain tensor is the strain-related contribution to extra dissipation due to a rigid inclusion. A trace-free version is equivalent for an incompressible flow. This tensor is distinct from the axisymmetric force-dipole flow used as a far-field singularity.

For a Newtonian fluid, write $\mathbf{\sigma }=-pI+2\mu e$, where $e=(\nabla u+\nabla u^T)/2$ is the rate-of-strain tensor. The Stokes equation without body force gives ${\partial }_j{\sigma }_{ij}=0$. Since the stress tensor is symmetric, its contraction with the antisymmetric part of the velocity gradient is zero. ThusThe divergence theorem proves the boundary-work identityFor an incompressible flow, $\mathrm{tr}e=0$, so $D=2\mu {\int }_{V}e:edV\ge 0$ is precisely the viscous dissipation.

Apply this identity to the particle-free volume and then to the fluid outside the added particle. On the rigid particle, the velocity has the form $V+\Omega \times x$. Its contribution to the boundary work is $V\cdot F+\Omega \cdot T$, which vanishes because the particle is force-free and torque-free. The outer velocity is unchanged. Subtracting the two boundary-work identities therefore givesTo move this expression to the particle surface, apply the Lorentz reciprocal theorem to $(u_0,{\sigma }_0)$ and $(u^{\prime },{\sigma }^{\prime })$ in the punctured fluid domain. For completeness, its local proof isOn the outer boundary $u^{\prime }=0$. Let $N=-n$ be the particle-outward normal, opposite to the fluid-outward normal there. Integrating the preceding divergence identity givesThis expression includes the removal of the original fluid inside the particle; that removal must not be counted separately again.

For a sphere in a uniform straining Stokes flow, $E$ is symmetric and trace free, and $u_0=Ex$. Symmetry gives zero rigid translation and rotation. Put $H=x\cdot Ex$ and $r=|x|$. Use the Unscaled Papkovich–Neuber representation printed on the paper, with the decaying harmonic potentialsTheir harmonicity follows from $x_i/r^3=-{\partial }_i(1/r)$ and from $H/r^5$ being a trace-free second derivative of $1/r$. Direct differentiation givesAt $r=a$, the no-slip boundary condition requires $u^{\prime }=-Ex$. Matching the independent $Ex$ and $Hx$ terms fixes $B=-1/2$ and $C=5/6$. The disturbance velocity and pressure are consequentlyThe pressure follows from $p^{\prime }=2\mu \nabla \cdot {\Phi }^{\prime }$. The disturbance decays, the total velocity vanishes at the sphere, and the harmonic representation verifies the Stokes equation. The boundary traction has zero total force by odd symmetry, and zero total torque because $E$ is symmetric, confirming the assumed force-free and torque-free motion.

At the sphere $u^{\prime }=-u_0$, so the extra dissipation due to a rigid inclusion becomes ${\int }_A{u}_0\cdot ({\sigma }_0+{\sigma }^{\prime })NdS$. Contracting the supplied total stress tensor with $N$ cancels its terms proportional to $(N\cdot EN)N$, leaving $5\mu EN$. ThusUsing ${\int }_A{N}_i{N}_{j}dS=4\pi a^2{\delta }_{ij}/3$ proves the single-sphere dissipation increaseOne may perform the comparison in a large bounded vessel with prescribed straining velocity and then take its size to infinity; the given near-sphere approximation makes this limit compatible with the fixed-outer-velocity definition of $D^{\prime }$.

The sphere number density is $n_s=\varphi /(4\pi a^3/3)=3\varphi /(4\pi a^3)$. Neglecting hydrodynamic interactions, the additional viscous dissipation per unit volume is $n_s{D}^{\prime }=5\mu \varphi E:E$. Adding the particle-free value $2\mu E:E$ givesThis is the Einstein viscosity formula for a dilute suspension, to first order in the particle volume fraction.

For a fixed fluid domain with prescribed velocity on its whole boundary, the Minimum-dissipation theorem for Stokes flow compares the actual incompressible Stokes flow $u$ with every admissible incompressible velocity field $v$ having the same boundary values. The trial fields need not themselves satisfy the momentum equation. With $e(u)=(\nabla u+\nabla u^T)/2$ and constant dynamic viscosity $\mu$, define the viscous dissipationAssume no nonconservative body force, or absorb a conservative force into the pressure. Set $w=v-u$, so $\nabla \cdot w=0$ and $w=0$ on the boundary. The Stokes stress tensor $\sigma =-pI+2\mu e(u)$ obeys $\nabla \cdot \sigma =0$. Because pressure has zero contraction with $e(w)$,Expanding the square provesEquality requires $e(w)=0$, so $w$ is a rigid motion; the homogeneous boundary values force it to vanish. The theorem concerns fixed boundary velocities, not fixed applied forces or arbitrary comparisons between different fluid domains.

For the concentric-sphere motion, seek $u=F(r)\mathbf{\Omega }\times x$. Incompressibility is automatic. The choice $F=A+B/r^3$ gives a harmonic velocity: the linear term is harmonic, and each component of $x/r^3$ is a derivative of the harmonic function $1/r$. The no-slip conditions $F(a)=1$, $F(b)=0$ giveThus the Rotational Stokes flow between concentric spheres isIt satisfies the Stokes equation with constant pressure, and the uniqueness implied by the theorem identifies it as the required flow.

In the Unscaled Papkovich–Neuber representation, take $\Phi =-u/2$ and $\chi =0$. Both potentials are harmonic and $x\cdot \Phi =0$, so the representation gives $u=-2\Phi$. Alsoand therefore $p=0$, after choosing the arbitrary pressure constant.

This pressure conclusion also follows from symmetry and Linearity of Stokes flow. Pressure is a scalar linear in the axial vector $\mathbf{\Omega }$, and isotropy leaves only a possible term proportional to $\mathbf{\Omega }\cdot x$. That is a pseudoscalar, forbidden by reflection symmetry of the concentric spherical domain. Hence the pressure perturbation vanishes. Equivalently, axisymmetry and reflection with flow reversal make the motion purely azimuthal, while its meridional momentum balance has no pressure gradient.

Let $v=\mathbf{\Omega }\times x$ in the following stress calculation. The symmetric velocity gradient isThe full stress tensor, in the zero-pressure gauge, is thereforeFor $n=x/r$, the traction on a radial surface is $\sigma n=-3\mu B(\mathbf{\Omega }\times n)/r^3$. This completely specifies the stress field; with the polar axis along $\mathbf{\Omega }$, its only nonzero spherical shear component is ${\sigma }_{r\varphi }=-3\mu B\Omega \sin \theta /r^3$.

The fluid couple on the inner sphere is the integral of $x\times \sigma n$ with the normal pointing from the solid into the fluid. Using $\int n_i{n}_{j}d\omega =(4\pi /3){\delta }_{ij}$,Hence the sustaining applied Torque in rotational Stokes flow between concentric spheres isFor $a\ll b$, it approaches the isolated rotating sphere in Stokes flow result $8\pi \mu a^3\mathbf{\Omega }$. For a narrow gap $d=b-a\ll a$, it becomes $8\pi \mu a^4\mathbf{\Omega }/(3d)$: the divergent resistance is that of a thin Couette shear layer, with speed scale $a\Omega$ and gradient scale $a\Omega /d$.

Now add the force-free, couple-free rigid particles and denote the resulting velocity by $u_p$. Extend $u_p$ through each particle as its rigid translational and rotational velocity. No slip makes this extension continuous across every particle surface; it is incompressible and has zero rate-of-strain tensor inside each particle. It therefore gives an admissible trial velocity in the original particle-free annulus, with exactly the same sphere boundary velocities. Its full-annulus dissipation is the actual fluid dissipation ${\mathcal{D}}_p$, since the added interior integrals are zero. The Minimum-dissipation theorem for Stokes flow consequently gives ${\mathcal{D}}_p\ge {\mathcal{D}}_0$.

The boundary-work identity follows by integrating $\nabla \cdot (\sigma u)$: dissipation equals the mechanical power supplied at all rigid boundaries. The outer sphere is stationary; the inner sphere has no translation, so any force needed to hold its centre does no work; and each added particle contributes zero power because its net force and couple are zero. Thus ${\mathcal{D}}_p={\mathbf{G}}_p\cdot \mathbf{\Omega }$, while ${\mathcal{D}}_0={\mathbf{G}}_0\cdot \mathbf{\Omega }$. ThereforeFor $\mathbf{\Omega }\ne 0$ and at least one particle of positive volume in the open annulus, the inequality is strict. Equality would make the extended velocity identical to the particle-free solution, yet that solution has nonzero strain on every open region away from the rotation axis and cannot be rigid throughout any particle interior. Hence the applied couple component along the imposed angular velocity increases. This is the fixed-velocity extra dissipation due to a rigid inclusion argument; it does not assert an increase in every transverse torque component.

Let downward velocity be positive in this part, and let $W>0$ be the dense sphere's weight minus buoyancy. For the original centreline configuration, write its six-component translation-and-rotation velocity as $q$, and its diagonal positive hydrodynamic resistance matrix as $R_0$. Its single-sphere dissipation isWith only axial force and no torque, the original velocity is $V_0=W/R_{zz}$.

Now add the neutral sphere. In its actual Stokes flow, extend the velocity rigidly through the new sphere's interior. This extension is a valid incompressible trial field for the domain without that sphere, with exactly the dense sphere's actual six boundary velocities; its strain and dissipation inside the inserted ball are zero. The Minimum-dissipation theorem for Stokes flow and extra dissipation due to a rigid inclusion therefore giveThe neutral sphere has zero net external force and torque, and the dense sphere has only the axial force $W$. The dissipation equals total external power, so $D=W{V}_z$. Positivity gives $V_z>0$, and the comparison yields $V_z\le W/R_{zz}$ without assuming that either sphere translates only vertically or that either angular velocity vanishes.

For a finite nonoverlapping placement, the first inequality is strict. Equality would make the extended field the old Stokes solution and thus make that solution rigid throughout the added ball. Analyticity and connectedness then force its strain to vanish throughout the original fluid domain; the stationary tube wall and quiescent far field would force the whole motion to vanish, contradicting $V_z>0$. Hence the force-free inclusion reduces axial mobility of a centred settling sphere:The difference tends to zero as the added sphere recedes infinitely far away. The lateral and rotational contributions are nonnegative in the resistance comparison, which is precisely why they cannot invalidate the bound.

For an incompressible flow, write the Newtonian fluid stress tensor as $\sigma =-pI+2\mu e$, where $e=(\nabla u+\nabla u^T)/2$ is the rate-of-strain tensor. The Stokes equation gives ${\partial }_j{\sigma }_{ij}=0$. Symmetry of the Cauchy stress tensor therefore givesIntegrating and applying the divergence theorem expresses the viscous dissipation as boundary power:At a moving rigid body, the surface velocity is $U+\Omega \times (x-x_c)$. Its boundary-power contribution is $U\cdot F+\Omega \cdot G$, with the force and torque evaluated using the fluid's outward normal vector. Both resultants vanish for a force-free, torque-free inclusion. Thus the new viscous dissipation comes entirely from the unchanged outer boundary velocity. Subtracting the particle-free boundary power givesApply the Lorentz reciprocal theorem to $(u_0,{\sigma }_0)$ and $(u^{\prime },{\sigma }^{\prime })$ in the fluid outside the inclusion. On the outer boundary $u^{\prime }=0$, soWriting $n_p=-n$ for the particle's outward normal vector gives the required extra dissipation due to a rigid inclusion:This subtraction already accounts for the fluid volume displaced by the inclusion; it is not just the integral of the disturbance's local viscous dissipation.

For the sphere, the ambient rate-of-strain tensor $E$ is symmetric and trace free. The sphere in a uniform straining Stokes flow has zero translational and rotational velocity: inversion symmetry eliminates its force, and symmetry of $E$ eliminates its torque. To derive the disturbance, put $r=|x|$, $Q=x\cdot Ex$ and seekThe incompressible flow condition and Stokes equation reduce toA decaying family satisfying these equations is $f=A{r}^{-5}$, $g=B{r}^{-5}-5A{r}^{-7}/2$, $j=2B{r}^{-5}$. The no-slip boundary condition requires $f(a)=-1$ and $g(a)=0$, so $A=-a^5$ and $B=-5a^3/2$. ConsequentlyThe total velocity is zero at $r=a$, and the disturbance decays as $r^{-2}$. The pressure constant has been set to zero. These boundary and far-field conditions, together with Uniqueness of Stokes flow, establish the solution. There is no rotational background in a pure strain flow.

On the sphere $u^{\prime }=-u_0=-aE{n}_p$. Contracting the supplied Newtonian fluid stress tensor with $n_p$ makes its two terms proportional to $(n_p\cdot E{n}_p)n_p$ cancel, leaving $\sigma n_p=5\mu E{n}_p$. Hence the extra dissipation due to a rigid inclusion isHere ${\int }_{S^2}{n}_j{n}_{k}d\Omega =4\pi {\delta }_{jk}/3$. Taking the large outer boundary limit only after using the fixed-boundary power identity avoids replacing that prescribed boundary by an uncontrolled boundary at infinity.

If $N_v$ denotes the number of spheres per unit volume, $N_v(4\pi a^3/3)=\varphi$. ThusAdd this to the particle-free viscous dissipation density $2\mu E:E$. Comparing with $2{\mu }_{\mathrm{eff}}E:E$ gives the Einstein viscosity formula for a dilute suspension:to first order in $\varphi$. The hydrodynamic interactions neglected here contribute beyond that dilute order.

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.