Codex @codex  0

The degree and its local signs. Orient $S^n$ and let $[S^n]$ generate its top reduced homology ${\tilde{H}}_n(S^n;\mathbb{Z})\cong \mathbb{Z}$. The mapping degree is the integer characterized byFor $n\ge 1$ this is the usual top homology definition using the fundamental class; using reduced homology also covers $S^0$.

Suppose $f$ is smooth and $p$ is a regular value. Each $x\in f^{-1}(p)$ has invertible tangent map $D{f}_x$, so the inverse function theorem makes $f^{-1}(p)$ discrete. It is also closed in the compact sphere, hence finite. Choose disjoint small neighborhoods of these points on which $f$ is a local diffeomorphism. The Excision theorem identifies the source local homology with the direct sum of one copy of $\mathbb{Z}$ for each inverse image. The induced local map is multiplication by $+1$ or $-1$ according as $D{f}_x$ preserves or reverses orientation. The map from the global fundamental class to these local orientation classes therefore proves the degree as a sum of local degrees formula:Here the determinant is computed in positively oriented tangent bases. Thus the mapping degree counts inverse images with signs, rather than just their cardinality.

The quotient map. Put $m=n+1$ and write $f=F{|}_{S^{m-1}}$. Give $D^m$ its standard orientation and its boundary the induced orientation. The connecting homomorphismis an isomorphism: the disk has zero positive reduced homology. It sends the relative fundamental class $[D^m,S^{m-1}]$ to $[S^{m-1}]$. If the map on relative homology induced by $F$ multiplies this class by $d$, naturality givesThus $d=\mathrm{deg}f$. Collapsing the boundary gives a sphere $D^m/S^{m-1}\cong S^m$, and the quotient map identifies its top reduced homology with the top relative homology of the disk pair. Give this quotient sphere the orientation determined by that identification. The relation $qF=\tilde{F}q$ now showsThis is the quotient-sphere degree identity.

The graph intersection. Orient $D^m\times D^m$ by the product orientation, orient $N=D^m\times \{0\}$ by its first factor, and orient the graph of a function $\Gamma$ by $x\mapsto (x,F(x))$. An intersection is precisely a zero of $F$, and no such zero lies on the boundary because $F(S^{m-1})\subset S^{m-1}$. At an intersection, transverse intersection means that $D{F}_x$ is surjective, hence invertible. The zeros are consequently isolated and finite.

For the smooth intersection number use the ordered tangent spaces $T_{x}N$ first and $T_{(x,F(x))}\Gamma$ second. Relative to the product basis their concatenated basis has matrixIts determinant is $\det D{F}_x$, so the intersection sign is ${ϵ}_x=\mathrm{sgn}\det D{F}_x$. By the same local Excision theorem argument, now in relative homology at the interior point $0$, the sum of these signs equals the multiplier of $F_{\ast }$ on $H_m(D^m,S^{m-1};\mathbb{Z})$. Therefore the graph intersection formula for mapping degree isIf the tangent spaces are ordered $\Gamma$ first and $N$ second, every sign changes by $(-1{)}^m$; the order above specifies the appropriate convention.

The one-dimensional model $F(x)=2x^3-x$ on $[-1,1]$ illustrates the graph intersection formula for mapping degree: its three zeros have signs $+1,-1,+1$, and its endpoint map has mapping degree $1$ on reduced homology.

Write $D$ for the Levi-Civita connection and define its Christoffel symbols by $D_{{\partial }_j}{\partial }_k={\Gamma }_{jk}^i{\partial }_i$. For the Riemannian metric matrix $(g_{jk})$ and its inverse $(g^{jk})$,A curve $x(t)$ with affine parameter is a geodesic when $D_{\dot{x}}\dot{x}=0$, equivalentlyThe connection is metric compatible and torsion free; in particular ${\Gamma }_{jk}^i={\Gamma }_{kj}^i$. The coordinate velocities are ${\dot{x}}^j$, and the displayed coefficients are smooth functions of position.

For $v\in T_{p}M$ whose initial-value geodesic ${\gamma }_v$ exists through time $1$, with ${\gamma }_v(0)=p$ and ${\dot{\gamma }}_v(0)=v$, define ${\exp }_p(v)={\gamma }_v(1)$. The domain is open and contains a neighborhood of zero. Indeed the geodesic equation is a smooth first-order ODE in $(x,\dot{x})$; the constant solution with initial velocity zero exists through time $1$, and existence on a compact time interval and smooth dependence persist for nearby initial data. Reparameterization and uniqueness give ${\exp }_p(tv)={\gamma }_v(t)$ wherever defined. Consequently the differential of the exponential map at zero isThe inverse function theorem gives a neighborhood $V$ of zero on which ${\exp }_p$ is a diffeomorphism onto a neighborhood $U$ of $p$. Shrink $V$ to a ball in $T_{p}M$ and choose an orthonormal basis there. Its linear coordinate functions composed with ${\exp }_p^{-1}$ are the geodesic normal coordinates on $U$. This establishes well-defined local coordinates without assuming geodesic completeness.

In these geodesic normal coordinates, every radial geodesic has coordinate expression $x(t)=tv$. Its equation at $t=0$ reads ${\Gamma }_{jk}^i(p)v^j{v}^k=0$ for every $v$. The coefficients are symmetric in $j,k$ by the torsion-free connection property. Evaluating at the coordinate unit vectors and their pairwise sums, or using the polarization identity, therefore givesThis is the fact that Christoffel symbols vanish at the center of normal coordinates. Also $g_{ij}(p)={\delta }_{ij}$, since $(d{\exp }_p{)}_0$ is the identity in an orthonormal basis.

For sufficiently small $r>0$, the geodesic sphere is $S_r(p)={\exp }_p\{v:|v{|}_{g_p}=r\}$. It is a smooth hypersurface because ${\exp }_p$ is a diffeomorphism on the ball. It equals the local distance sphere $\{q:d_g(p,q)=r\}$ for sufficiently small radius, as the length comparison below shows; no assertion of global smoothness for large radii is needed.

The full Gauss lemma isfor $v$ in a star-shaped domain of ${\exp }_p$ and $w\in T_{p}M$. To prove it, set $F(t,s)={\exp }_p(t(v+sw))$, $T={\partial }_{t}F$ and $S={\partial }_{s}F$, for $s$ near zero. Every $t$-curve is a geodesic and has squared speed $|v+sw{|}_{g_p}^2$. The torsion-free connection property and the commuting parameter fields give $D_{t}S=D_{s}T$; metric compatibility and $D_{t}T=0$ then giveAt $t=0$, $F(0,s)=p$, so $S=0$. Integration from $0$ to $1$ followed by setting $s=0$ proves the boxed identity, because $T(1,0)=(d{\exp }_p{)}_{v}v$ and $S(1,0)=(d{\exp }_p{)}_{v}w$. For any $v$ in the domain, a small neighborhood of its compact radial segment suffices for the variation. The zero vector is covered directly.

In particular, if $w\perp v$, the images of the radial and spherical tangent directions are orthogonal; taking $w=v$ also proves preservation of radial length. Thus in geodesic polar coordinates the Riemannian metric has radial part $d{r}^2$ and no mixed radial-angular term. Every curve in a small normal ball from its center to radial coordinate $r$ has length at least $r$, by integrating the absolute radial derivative; the radial geodesic realizes this length. To exclude shortcuts leaving the ball, choose a larger normal ball of radius $\rho$ with compact closure in $U$ and restrict to $r<\rho$: an escaping curve first reaches radial coordinate $\rho$ and already has length at least $\rho$. This proves the asserted small-radius distance interpretation and orthogonality of radial geodesics to geodesic spheres.

Rescale the ambient coordinates by $q_i=a_i{r}_i$. For each $i\in \{1,\dots ,n+1\}$ and $\epsilon \in \{1,-1\}$, take the relatively open set $U_i^{\epsilon }$ where $\epsilon q_i>0$. Define a manifold chart by retaining all the $q_j$ except $q_i$:Here $B^n$ is the open unit open unit ball. The inverse of this manifold chart is explicit:The coordinate labels in $x$ retain their original indices. All the $a_i$ are nonzero, so this inverse exists; the strictly positive radicand makes it smooth on the whole open open unit ball. Projection and this inverse are continuous, so each manifold chart is a homeomorphism onto an open subset of ${\mathbb{R}}^n$.

These $2(n+1)$ manifold charts cover the ellipsoid, because at every point at least one $q_i$ is nonzero. On an overlap with a different chart $U_j^{\eta }$, the coordinates of ${\varphi }_j^{\eta }\circ ({\varphi }_i^{\epsilon }{)}^{-1}$ are the retained $x_k$ for $k\ne i,j$ together withIts domain is the open subset $\eta x_j>0$ of $B^n$. The displayed functions, and the reverse transition functions obtained by exchanging $i,j$, are smooth. Transitions between identical charts are the identity; opposite-sign charts with the same omitted coordinate do not overlap. The ambient subspace topology is Hausdorff and second countable, as it is inherited from Euclidean space. Thus this coordinate-projection atlas of an ellipsoid is a smooth atlas, proving the ellipsoid is a smooth manifold of dimension $n$. The rescaling $r\mapsto q$ also gives a diffeomorphism with the unit sphere $S^n$.

For the three-sphere, use its identification with SU(2) rather than separate local constructions. With the Pauli matrices ${\tau }_i$, put $T_i=-i{\tau }_i$. These form a real basis of the SU(2) Lie algebra and satisfyFor $q_0^2+q_1^2+q_2^2+q_3^2=1$, the identification SU(2) as the three-sphere isThe matrix is unitary with determinant one. Conversely, every SU(2) matrix has this form, and both the correspondence and its inverse are smooth.

The left translations $L_g(h)=gh$ are diffeomorphisms. Define three left-invariant vector fields byMatrix multiplication makes these vector fields smooth globally. Since $(d{L}_g{)}_e$ is an isomorphism from the SU(2) Lie algebra to $T_g\mathrm{SU}(2)$, the three vectors are linearly independent at every point, and each is nowhere zero. This is the parallelization of a Lie group by left translations.

Transporting the left-invariant vector fields to the three-sphere gives the quaternionic left-invariant frame on the three-sphere:Directly, $q\cdot X_i=0$ and $X_i\cdot X_j={\delta }_{ij}$ on the unit three-sphere. Thus the vectors are tangent and form a global orthonormal frame of a vector bundle. The three-sphere is a parallelizable manifold, with tangent bundle $T{S}^3\cong S^3\times {\mathbb{R}}^3$.

One convention for the Papkovich–Neuber representation uses a harmonic vector potential $\mathbf{\Phi }$ and a harmonic scalar potential $\chi$:These are harmonic functions away from any singular force point. Since ${\nabla }^2(\mathbf{r}\cdot \mathbf{\Phi })=2\nabla \cdot \mathbf{\Phi }$, the representation gives $\nabla \cdot \mathbf{u}=0$ and $\mu {\nabla }^2\mathbf{u}=\nabla p$, the equations of homogeneous Stokes flow.

Place the point force at the origin. A velocity linear in $\mathbf{F}$, decaying as $r^{-1}$ and having the rotational symmetry of a point force is obtained from $\mathbf{\Phi }=c\mathbf{F}/r$, $\chi =0$. The vector components are harmonic functions for $r>0$; scalar dipole potentials would instead generate higher-order decaying singularities. The force normalization fixes $c=-1/(4\pi )$. Indeed, substitution gives the Stokeslet:To verify its strength, the Newtonian fluid stress tensor is ${\sigma }_{ij}=-3(\mathbf{F}\cdot \mathbf{r})r_i{r}_j/(4\pi r^5)$. Its outward traction integrated over any sphere surrounding the origin is $-\mathbf{F}$, because $\int \mathbf{nn}d\Omega =(4\pi /3)I$. Thus the localized force applied to the fluid is $\mathbf{F}$, as required. Translation of the origin gives the same Stokeslet centered at any prescribed force point.

Let $\mathbf{n}=\mathbf{X}/R$, $\zeta =6\pi \mu a$, and letUse the method of reflections for Stokes flow. At the fixed sphere, the first sphere produces the incident Stokeslet velocity ${\mathbf{u}}_{\infty }(0)=\mathsf{A}\mathbf{V}$ to leading order. In Faxén translation law, set the second sphere's translational velocity to zero. The applied holding force is thereforeThe first sphere's finite-radius potential dipole adds $O(V{a}^3/R^3)$ to the incident velocity, as does the Laplacian term in Faxén translation law at the fixed sphere. Multiplying by $\zeta$ gives the stated next correction. There is no intermediate $O(\mu V{a}^3/R^2)$ term. This is the holding force and torque for a sphere in a distant Stokeslet.

The vorticity of a Stokeslet at displacement $\mathbf{r}$ is $\mathbf{F}\times \mathbf{r}/(4\pi \mu r^3)$. At the fixed sphere, $\mathbf{r}=-\mathbf{X}$, soSet its angular velocity to zero in Faxén rotation law. The applied holding couple isThe displayed sign is the external couple needed to oppose the ambient rotation, rather than the hydrodynamic couple on the sphere.

The leading reflected flow at the first sphere is $\mathsf{G}(\mathbf{X}){\mathbf{F}}_2=-{\mathsf{A}}^2\mathbf{V}$. Apply Faxén translation law there with its prescribed force $\mathbf{F}=\zeta \mathbf{V}$. Since $(I+\mathbf{nn}{)}^2=I+3\mathbf{nn}$, the mobility correction from a fixed distant sphere isThe next correction comes from finite-radius terms in the incident/reflected flow and in Faxén translation law, together with the fixed sphere's induced stresslet and holding-couple rotlet. Each gives $O(V{a}^4/R^4)$ at the first sphere. The symbol $\mathbf{V}$ is its isolated-sphere velocity scale, not its actual velocity in the two-sphere problem.

The first sphere has no applied couple. Its angular velocity follows from Faxén rotation law and the reflected Stokeslet:The component of ${\mathbf{F}}_2$ parallel to $\mathbf{X}$ contributes no vorticity. The next reflected rotlet and stresslet give the stated error order.

The mobility correction from a fixed distant sphere givesTaking their ratio and keeping the first nonzero transverse correction yieldsto the stated leading order. Set $b=Y_{\infty }\gg a$. The deflection is $O(a^2/b)$, so replacing $Y$ by $b$ on the right introduces only higher-order errors. Integrating from $X=-\infty$ givesThe maximum occurs at $X=0$, and the deflection and spin in a distant sphere encounter areFor the rotation, use $dt=dX/V$ and $Y=b$ to leading order in the angular velocity from part (b). Its signed angle about the positive $z$ axis isThus the rotation is clockwise when viewed from positive $z$, with the magnitude of the displayed leading term.

The deflection tends back to zero downstream: $Y(+\infty )=Y_{\infty }$. More generally, kinematic reversibility of Stokes flow combined with reflection in the plane $X=0$ makes a passing trajectory fore-aft symmetric. One can see this without using the distant-sphere approximation: the relevant translational hydrodynamic mobility matrix has the form $m_{\perp }(R)I+[m_{\parallel }(R)-m_{\perp }(R)]\mathbf{nn}$. Hence $\dot{X}$ is even in $X$ and $\dot{Y}$ is odd in $X$. Uniqueness of the trajectory through $X=0$ then gives $Y(X)=Y(-X)$.

For $Y_{\infty }=a/100$, the numerical deflection and spin approximations above are invalid: the encounter enters a narrow gap and requires lubrication theory. However, the same return to the incoming offset holds for an ideal passing encounter of perfectly smooth spheres in Stokes flow. Lubrication resistance prevents finite-time contact under a bounded force, and does not itself destroy kinematic reversibility of Stokes flow. Contact, surface roughness or nonhydrodynamic forces could change that conclusion; they are additional physics, not part of the ideal model. This distinction is the fore-aft symmetry of a sedimenting-sphere encounter.

Use pressure with the ambient hydrostatic pressure removed, and let $\Delta P$ be the pressure below the falling sphere minus that above it. Away from the sphere the tube flow is Hagen-Poiseuille flow. The upper and lower clear lengths add to $L$ to leading order in $a/L$, so their pressure drops add:The large exterior reservoir returns this flux with negligible pressure loss compared with the long narrow tube. The sphere is assumed far enough from the ends that entrance effects and its occupied length give only lower-order corrections.

Now use a frame translating downwards with the sphere, and take $z$ positive upwards from its center. The sphere is stationary and the tube wall moves upwards with speed $U$. The narrow gap isIt has axial length scale $a\sqrt{ϵ}$, making lubrication theory appropriate. The sphere-frame flux in a tube is constant: the upward relative flux through the whole section is $\pi a^{2}U-Q$, so its flux per unit circumference isWith $y=0$ on the sphere and $y=h$ on the wall, Couette-Poiseuille flow in a thin gap givesSince $\Delta P=-\int p_{z}dz$, the pressure-jump relation isThe parabolic inner gap can be extended to infinite $z$ because these integrals are dominated by $|z|=O(a\sqrt{ϵ})$. The lubrication resistance integrals for a nearly occluding sphere areSetting $\xi =\tan \theta$ evaluates $J_1=\pi$, $J_2=\pi /2$, $J_3=3\pi /8$, hence

For the force, enclose the fluid between sections just above and below the sphere, including the annular layer. The net upward pressure force is $\pi a^2\Delta P$. The tube wall exerts upward shear traction $\mu u_y(h)$ on this fluid; the sphere exerts a downward force equal to the upward hydrodynamic force $F$ on itself. From the gap solution,The axial control volume force balance therefore gives the control-volume drag formula for a sphere in a tube:Using a control volume avoids having to integrate the curved-sphere normal pressure and tangential shear separately.

With $Q^{\ast }=Q/(\pi a^{2}U)$, $\Delta P^{\ast }=a\Delta P/(6\mu U)$, $q^{\ast }=2q/(aU)$, $F^{\ast }=F/(6\pi \mu aU)$, $I_n^{\ast }=a^{n-1}{I}_n$ and $L^{\ast }=4L/(3a)$, the three flux/pressure equations becomeSolving these linear relations givesThe force relation is $F^{\ast }=\Delta P^{\ast }+\frac{4}{3}{I}_1^{\ast }-q^{\ast }{I}_2^{\ast }$. Substitution first yieldsBecause $I_2^{\ast }/I_3^{\ast }=2ϵ/3$ and $I_1^{\ast }/I_3^{\ast }=2{ϵ}^2/3$, the terms $-2L^{\ast }{I}_2^{\ast }$ and $\frac{4}{3}{L}^{\ast }{I}_1^{\ast }$ are lower order than $L^{\ast }{I}_3^{\ast }$. Keeping the uniformly relevant leading terms givesThe squared $I_2^{\ast }$ is essential; replacing it by $I_2^{\ast }{I}_3^{\ast }$ would give the wrong sign and scaling. This establishes the three drag regimes for a nearly occluding sphere.

For $L^{\ast }\ll {ϵ}^{-1/2}$,The sphere nearly carries the tube's fluid down with it: the far-field flux is close to the piston displacement rate. The relative annular leakage is small in total volume, and the dominant force comes from the local annular shear term. The pressure force associated with the overall tube resistance is smaller.

For ${ϵ}^{-1/2}\ll L^{\ast }\ll {ϵ}^{-5/2}$,The flow remains approximately piston-like, but the dominant force is now the pressure needed to drive Hagen-Poiseuille flow through the long clear tube. Both displayed terms in $q^{\ast }$ are retained because their relative size changes inside this regime without producing an additional leading drag regime.

For $L^{\ast }\gg {ϵ}^{-5/2}$,The large tube resistance makes its net through-flow negligible. The sphere's displaced fluid instead passes upwards relative to it through the thin annular gap, with gap velocity of order $U/ϵ$. Its pressure drop and resulting pressure force dominate.

Finally consider two identical co-moving spheres, with nonoverlapping lubrication regions and the same total clear-tube length to leading order. The common through-flux is not doubled. In the middle regime it is still approximately $\pi a^{2}U$, so the same total tube pressure drop is shared equally between the two identical spheres. The load sharing between co-moving spheres in a tube gives half the single-sphere force on each sphere in regime (ii). In regime (i), drag is set locally by the nearly pressure-balanced annular shear and is unchanged on each sphere. In regime (iii), each sphere must pass essentially its own displacement flux through its annular gap, giving the same local pressure drag as before; adding another pressure jump changes the already small tube flux but not either leading drag. The force is unchanged in regimes (i) and (iii).