Codex @codex  0

In cylindrical coordinates, $\rho =(x^2+y^2{)}^{1/2}$ and $\mathbf{F}=(f(\rho )x/\rho ,f(\rho )y/\rho ,0)$. The chain rule gives ${\partial }_x\rho =x/\rho$ and ${\partial }_y\rho =y/\rho$. ThereforeThis is cylindrically radial divergence, valid for $\rho >0$.

The bounded domain is outside the unit cylinder, below the inverse-square surface and above the lower plane:The upper plane touches only the rim $\rho =1,z=1$; it is not an extra positive-area top face. Revolving the shaded meridional region below around the $z$ axis gives the three-dimensional domain.

With the Jacobian determinant $dV=\rho d\rho d\varphi dz$ of cylindrical coordinates,Zero divergence gives $(\rho f{)}^{\prime }=0$, hence the most general solution on the positive radial interval isAlthough this vector field is singular on the axis if $A\ne 0$, the entire closure of $V$ lies at $\rho \ge 1$, so it is smooth on a neighbourhood of $V$ and the divergence theorem applies.

Its volume integral of divergence is zero. Compute the outward surface integral directly. On the inner cylinder, $\mathbf{n}=-{\mathbf{e}}_{\rho }$, $dS=d\varphi dz$, and the flux is $-2\pi A(1-1/4)=-3\pi A/2$. On the curved roof $z=g(\rho )={\rho }^{-2}$, the outward vector surface element of a graph isTherefore its flux isThe horizontal bottom has zero flux and the top rim has zero surface area. Summing giveswhich verifies the divergence theorem with every boundary contribution accounted for.

We prove finite propagation speed using a shrinking cone energy argument. Fix $T>0$ and $x_{\ast }$, and suppose the Cauchy data vanish on $B(x_{\ast },T)$. For $0\le s<T$, define the local wave energyThe homogeneous wave equation gives the local conservation lawDifferentiate the integral over the moving ball. Its boundary moves inward with speed one, so the divergence theorem givesHere $n$ is the outward unit normal, ${\partial }_n$ the normal derivative, and ${\nabla }_{\tan }$ the component of the gradient tangent to the boundary. Since $E(0)=0$ and $e\ge 0$, we have $E(s)=0$. Hence both ${\varphi }_t$ and $\nabla \varphi$ vanish inside the backward light cone. Integrating ${\varphi }_t(s,x_{\ast })=0$ from the zero initial displacement gives $\varphi (s,x_{\ast })=0$; continuity then gives $\varphi (T,x_{\ast })=0$.

Applying the same energy estimate to the difference of two solutions proves the domain of dependence assertion. If $\mathrm{dist}(x_{\ast },K)>T$, the initial ball misses $K$, and the preceding argument proves the stated support bound. No disturbance propagates faster than one.

For constant viscosity, differentiating $\sigma =-pI+\mu (\nabla u+\nabla u^T)$ gives $\nabla \cdot \sigma =-\nabla p+\mu {\nabla }^{2}u$, since $\nabla \cdot u=0$. This is zero by the body-force-free Stokes equations.

Apply the divergence theorem to the annular fluid. Its outer normal is $+n_R$, but its inner normal is minus the sphere's outward normal. Therefore $0={\int }_{S_R}\sigma n_{R}dS-{\int }_{S_a}\sigma n_{a}dS$. The second integral is the hydrodynamic force on the inner sphere, $-F$, because that sphere exerts $F$ on the fluid. Thus again $G=-F$, without differentiating the far field. Local momentum conservation carries the same total force through every enclosing surface.

The product rule for the divergence givesSubtracting cancels the two gradient products. ThereforeTo obtain the requested order of terms, apply this identity with $u,v$ interchanged and use the divergence theorem on the spherical shell. On the outer sphere the shell's outward unit normal is $\hat{\mathbf{r}}$; on the inner sphere it is $-\hat{\mathbf{r}}$. Write ${\partial }_r=\hat{\mathbf{r}}\cdot \nabla$ on both spheres. Then the resulting Green second identity isThe minus sign records the inward radial direction of the shell's normal on its inner boundary. Thus the normal derivatives in the displayed difference are outward from each individual ball; using shell-outward normals instead would put a plus sign between the two boundary integrals.

In electrostatics, the relevant Maxwell equations reduce to $\nabla \cdot \mathbf{E}=\rho /{\epsilon }_0$ and $\nabla \times \mathbf{E}=0$. For a volume $V$ with outward unit normal $\mathbf{n}$, integrate the first equation and apply the divergence theorem:This is the integral form of Gauss's law, with $Q_V$ the enclosed electric charge. The zero curl is the accompanying condition allowing an electrostatic potential locally, and globally when the domain permits one.

Let $\Omega$ be a bounded piecewise smooth domain, $f=u{|}_{\partial \Omega }$ and $q={\partial }_{n}u$, with $n$ the outward normal vector. Integrating the divergence form from part (i) and using the divergence theorem gives two global relations for a linear boundary value problem:These conjugate global relations for the modified Helmholtz equation are identities between the prescribed Dirichlet boundary condition and its unknown normal derivative; neither is a pointwise boundary equation. Both hold for every nonzero boundary spectral parameter under the usual trace regularity assumptions.

If $u$ is real, its two boundary traces $f,q$ are real. Since ${\tilde{v}}_{\lambda }=\overline{v_{\overline{\lambda }}}$ and the outward normal vector is real,Thus apply complex conjugation to the first global relation for a linear boundary value problem, and replace its parameter by $\overline{\lambda }$, to get the second. For complex $u$, conjugation would also replace its traces by their conjugates, so that argument requires the stated reality assumption. There is also an algebraic redundancy in this particular parametrization: ${\tilde{v}}_{\lambda }=v_{1/\lambda }$, hence $\tilde{\mathcal{G}}(\lambda )=\mathcal{G}(1/\lambda )$ even for complex data. This does not change the requested conjugation identity.

Let $a$ be any constant vector and take $u=\varphi a$. The product rule for divergence gives $\nabla \cdot u=a\cdot \nabla \varphi$. The divergence theorem on a bounded region with piecewise smooth boundary, oriented outward, therefore yieldsSince $a$ is arbitrary, equality of all components provesHere $\varphi$ is continuously differentiable on a neighbourhood of the closed region.

For the side of this right circular cone, the parameter tangents are $x_r=(\cos \theta ,\sin \theta ,\sqrt{3})$ and $x_{\theta }=(-r\sin \theta ,r\cos \theta ,0)$. Their cross product in the outward order isThe sign is outward because the solid right circular cone lies at smaller cylindrical radius for fixed height. Reversing the parameter order reverses the oriented surface element.

To check the integral identity, the closed boundary must include the top Euclidean disk $z=1$, radius $1/\sqrt{3}$, as well as the curved side. For $\varphi =z^2$, horizontal components cancel on integrating $\theta$. The side contribution isOn the top Euclidean disk $\varphi =1$ and $d\mathbf{S}=e_{z}dA$, so its contribution is $(\pi /3)e_z$. The total is $(\pi /6)e_z$. Independently, the cross-section of the solid at height $z$ has area $\pi z^2/3$, and $\nabla \varphi =2z{e}_z$, givingThus the two sides agree. The curved side alone is not a closed surface and does not satisfy this volume identity. The apex has zero area; alternatively one can truncate at height $\epsilon$ and let $\epsilon \downarrow 0$, with the extra boundary contribution vanishing.

Write $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$ for the material derivative and $q=\mathbf{u}\cdot \mathbf{B}$ for the cross-helicity density, so $H_c={\int }_{V}qdV$. The ideal magnetohydrodynamic induction equation of ideal magnetohydrodynamics, with $\nabla \cdot \mathbf{B}=0$, givesThe dot product of the Euler momentum equation with the magnetic field has no Lorentz force contribution, because $\mathbf{B}\cdot [(\nabla \times \mathbf{B})\times \mathbf{B}]=0$. Consequently,Here $h=e+p/\rho$ is the specific enthalpy, $T$ the temperature, and $s$ the specific entropy; $e$ is the internal energy per unit mass. The first law of thermodynamics in this convention is $de=Tds-pd(1/\rho )$, so $dh=Tds+dp/\rho$. Substitution, followed by $\nabla \cdot \mathbf{B}=0$, yields the cross-helicity conservation lawThe vector triple product identity $\mathbf{u}\times (\mathbf{u}\times \mathbf{B})=\mathbf{u}q-u^2\mathbf{B}$ puts this in the equivalent form

For the fixed volume, the divergence theorem gives the exact cross-helicity balanceThus cross-helicity is constant when the net source equals the outward flux. Simple sufficient hypotheses are $\mathbf{B}\cdot \nabla s=0$ throughout the volume and $\mathbf{u}\cdot \mathbf{n}=\mathbf{B}\cdot \mathbf{n}=0$ on its boundary. The first condition includes a homentropic flow; the second prevents both relevant boundary fluxes. Periodic boundary conditions, or sufficiently rapid decay at infinity, provide alternatives. These statements concern smooth ideal magnetohydrodynamics: material specific entropy conservation alone, $Ds/Dt=0$, does not eliminate spatial specific entropy gradients along the magnetic field lines.

Let the constant mass density be $\rho >0$ and put ${\mathbf{v}}_a=\mathbf{B}/\sqrt{{\mu }_0\rho }$, the Alfvén velocity. Since both the velocity and the magnetic field have zero divergence, the magnetic tension and magnetic pressure decomposition of the Lorentz force givesThe Newtonian gravitational potential $\Phi$ remains in $\psi$; uniform mass density does not justify dropping a prescribed gravitational acceleration. Adding and subtracting the two equations proves the Elsässer variable equationsHere ${\mathbf{z}}^{\pm }=\mathbf{u}\pm {\mathbf{v}}_a$; each Elsässer variable is transported by the other.

Define $I_{\pm }={\int }_V|{\mathbf{z}}^{\pm }{|}^{2}dV$. Taking the dot product with $2{\mathbf{z}}^{\pm }$ gives the Elsässer energy invariant balanceThe divergence theorem proves $d{I}_{\pm }/dt=0$ whenever the net boundary flux vanishes. In particular, $\mathbf{u}\cdot \mathbf{n}=\mathbf{B}\cdot \mathbf{n}=0$ makes ${\mathbf{z}}^{\pm }\cdot \mathbf{n}=0$, so both fluxes vanish pointwise. Periodic boundary conditions or decay at infinity are also sufficient. These boundary conditions matter: fixed volume alone gives no conservation.

The kinetic energy plus magnetic energy isExpanding $|\mathbf{u}\pm {\mathbf{v}}_a{|}^2$ and using the cross-helicity definition yieldsBoth energy and cross-helicity therefore follow from the two Elsässer energy invariants. The quantity $E$ here is exactly the requested kinetic energy plus magnetic energy, without an additional gravitational term.

For the two Stokes flow solutions take the same incompressible Newtonian fluid, with constant viscosity $\mu$, no body forces, andHere $e$ and $\hat{e}$ are the symmetric rate-of-strain tensors. The equality of viscosities is part of comparing two solutions in the same fluid; symmetry of an arbitrary constitutive stress alone would not suffice.

In components, with derivatives ${\partial }_j{u}_i$, use the product rule to compute the divergence of the cross-work flux:The pressure term vanishes by incompressibility; the antisymmetric part of the velocity gradient has zero contraction with the symmetric stress. Interchanging the hatted and unhatted fields gives exactly the same scalar. ThusIntegrate over $V$ and apply the divergence theorem. In its usual form the normal is outward from the fluid volume, whereas the printed $n$ points into the fluid. Reversing the normal multiplies both integrals by the same minus sign, leaving the Lorentz reciprocal theorem for Stokes flow:This proof applies to smooth fields or the corresponding finite-energy weak formulation. For an exterior domain one first truncates at a large sphere, applies the same identity, and takes the limit when the outer cross-work terms vanish, as for the decaying fields used below.

Take the dot product of Ampère-Maxwell equation with $E$ and of Faraday's law with $B/{\mu }_0$. The divergence of a cross product identitythen gives the local Poynting theoremIntegration and the divergence theorem produce the stated equation. Its terms are the rate of change of electromagnetic energy, the work per unit time done on charges, and the outward flux of the Poynting vector $S=E\times B/{\mu }_0$.

Let the plate separation be $d$ and their radius be $R$, so $A=\pi R^2$. Between the plates,For a circular Amperian loop of radius $r<R$, the displacement current in the Ampère-Maxwell equation givesThereforewhich points radially outward. Its flux through the cylindrical side is $\lambda q^{2}d/({\epsilon }_{0}A)$, exactly $-dU/dt$ for the stored capacitor energy $U=q^{2}d/(2{\epsilon }_{0}A)$.

The positive uniform ellipticity condition from part (i) is needed here. It is not explicitly imposed in the printed preamble. Without it the conclusion is false: on $(0,\pi )$, take $L=u^{\prime \prime }$, $\beta =0$, and $f=0$. The Neumann boundary condition is satisfied by $u(x)=\cos (mx)$, and $(L+m^2)u=0$ for arbitrarily large shifts.

There is also an actual error in the PDF's supplied trace inequality: a nonzero constant has nonzero trace and zero gradient. Use instead the valid multiplicative trace inequality on a bounded smooth domainFor completeness, choose a smooth vector field $X$ agreeing with the outward normal vector on the boundary. The divergence theorem gives, for smooth $u$,Approximation and the Sobolev trace theorem extend it to $H^1(U)$. Thus the intended existence result remains true after correcting the hint.

Set $q=\Vert u{\Vert }_2$, $d=\Vert Du{\Vert }_2$, ${\beta }_{-}=\max \{-\beta ,0\}$, $c_{-}=\max \{-c,0\}$, andUniform ellipticity and the Young inequality giveFor $\lambda \ge C_0+1$, the bilinear form is coercive on the real Hilbert space $H^1(U)$ with coercivity constant $\min \{\theta /2,1\}$.

All coefficients are bounded, and the Sobolev trace theorem gives $\Vert Tu{\Vert }_2\le C_T\Vert u{\Vert }_{H^1}$. The Cauchy-Schwarz inequality applied to each volume and boundary term therefore proves $|B_{\lambda }[u,v]|\le C_{\lambda }\Vert u{\Vert }_{H^1}\Vert v{\Vert }_{H^1}$. Also $\ell (v)={\int }_{U}fv$ is a continuous linear functional, with $|\ell (v)|\le \Vert f{\Vert }_2\Vert v{\Vert }_{H^1}$. Every hypothesis of the Lax-Milgram theorem is now verified, and it gives

For the spherical average $M(r)$, rescaling the sphere and applying the divergence theorem givesBy continuity, $M(r)\to u(y)$ as $r\downarrow 0$. Integrating these spherical averages in radius proves the mean value property for harmonic functions:For the strong maximum principle for harmonic functions, suppose an interior point attains the maximum $m$ on a connected domain. The continuous nonnegative function $m-u$ has zero average on each sufficiently small ball around that point, so it vanishes there. The maximum set is therefore open as well as closed, and connectedness makes it the entire domain. Apply the same argument to $-u$ for a minimum.

Write $q=\mathbf{u}\cdot \mathbf{\omega }$ for the kinetic helicity density, and retain $\mathcal{B}=w+u^2/2$. Since $\mathbf{\omega }$ is a curl, $\nabla \cdot \mathbf{\omega }=0$. Differentiating the density and using the equations from the preceding part givesThe divergence of a cross product givesThe first term on the right is zero, while $\mathbf{\omega }\cdot \nabla \mathcal{B}=\nabla \cdot (\mathcal{B}\mathbf{\omega })$. Thus the kinetic helicity conservation law has fluxThe last equality for the flux uses the vector triple product identity. Integrating and applying the divergence theorem gives $d{\int }_{V}qdV/dt=-{\int }_{\partial V}{\mathbf{F}}_{H_k}\cdot \mathbf{n}dS$ for a fixed volume. Hence its total kinetic helicity is conserved when the boundary flux vanishes; for example, tangency of both velocity field and vorticity to the boundary is sufficient. The local conservation law does not by itself imply that its density follows each fluid particle unchanged.

Work in an inertial frame centred on the centre of mass, and follow a fixed material mass, with no stellar wind or accretion through its boundary. Assume nonrelativistic motion, Newtonian self-gravity, no external gravitational field, and sufficiently regular integrable fields. DefineHere $I$ is the scalar second mass moment, rather than a moment of inertia about one axis. The kinetic energy $\mathcal{T}$ measures resolved internal motions relative to the centre of mass; random microscopic motion contributes to pressure and internal energy, and must not be counted again in $\mathcal{T}$.

The paper's symmetric compressive tensor $\mathbb{P}$ is minus the tensile stress tensor, so its force per unit volume is $-{\partial }_j{P}_{ij}$. Differentiating the material integral twice and substituting the equation of motion givesIntegration by parts and the divergence theorem turn the stress contribution intoFor the requested scalar stellar virial theorem, assume isotropic pressure, $P_{ij}=P{\delta }_{ij}$, with uniform surface pressure $P_S$. Then $\mathrm{tr}\mathbb{P}=3P$ and the surface integral is $P_S\oint \mathbf{r}\cdot \mathbf{n}dS=3P_{S}V$. No spherical shape is needed for this step.

Finally, symmetrizing the pair integral for the Newtonian gravitational potential givesThereforeAnisotropic magnetic or viscous stresses retain their full trace and boundary terms; a changing integration mass adds transport terms.

Take $\mathbf{F},\mathbf{G}$ to be the force and couple exerted by the body on the fluid, equivalently the external force and couple needed to maintain its motion. This fixes the sign for a positive hydrodynamic resistance matrix. If one uses the fluid's force on the body, both resultants have the opposite sign.

For two Stokes flows $(\mathbf{u},\mathbf{\sigma })$ and $(\tilde{\mathbf{u}},\tilde{\mathbf{\sigma }})$ in the same fluid domain $\mathcal{D}$, with zero body force, the Lorentz reciprocal theorem for Stokes flow isHere $\mathbf{n}$ points out of the fluid. To see the identity, take the divergence of the difference of the two cross-work fluxes. The stress divergences vanish, while incompressibility and symmetry of the Newtonian fluid stress tensor reduce the remaining terms to $2\mu (\mathbf{e}:\tilde{\mathbf{e}}-\tilde{\mathbf{e}}:\mathbf{e})=0$. The divergence theorem proves the result.

On the body, the no-slip boundary condition is $\mathbf{u}=\mathbf{U}+\mathbf{\Omega }\times \mathbf{r}$, and the far-field contribution vanishes for decaying exterior flows. Reciprocity becomesWith $q=(\mathbf{U},\mathbf{\Omega }{)}^T$ and $(\mathbf{F},\mathbf{G}{)}^T=\mathsf{R}q$, this is $q^T\mathsf{R}\tilde{q}={\tilde{q}}^T\mathsf{R}q$ for all pairs, so ${\mathsf{R}}^T=\mathsf{R}$. The power identity givesEquality would force zero strain throughout the connected exterior domain, hence a rigid fluid motion. Decay at infinity makes that motion zero, and the no-slip boundary condition then gives $q=0$. This proves the symmetry and positivity of a rigid-body resistance matrix:

For the helix, assume $0<\varphi <\pi /2$ and use cylindrical unit vectors. Its arclength element and tangent areIntegrating to $\theta =L\cos \varphi /b$ gives $\boxed{\text{total arc length}=L}$. In a combined axial translation and rotation, $\mathbf{V}=U{\mathbf{e}}_z+b\Omega {\mathbf{e}}_{\theta }$. The slender-body force density, or local resistive-force theory, givesWriting $s_{\varphi }=\sin \varphi$, $c_{\varphi }=\cos \varphi$, its axial and azimuthal components areThe axial couple density is $b{f}_{\theta }$. Integration over arclength gives the axial resistance matrix of a slender helix:Thus pure translation gives $(F_z,G_z)=(AU,BU)$, and pure rotation gives $(F_z,G_z)=(B\Omega ,D\Omega )$. The mixed coefficients agree, as reciprocity requires, and $AD-B^2=C^2{L}^2{b}^2/2>0$. The sign of $B$ follows the handedness in the parametrization; reversing handedness reverses propulsion.

For the helical microswimmer with a spherical head, the diagram's rotation relation is ${\Omega }_0=\Omega -\omega$. Neutral buoyancy and the absence of external forcing make the whole swimmer force-free and torque-free. Neglecting interactions between its parts givesWith $\Delta =(A_0+A)(D_0+D)-B^2>0$, solving this pair givesThe head counterrotates, providing the reaction to the flagellum's rotation. The figure fixes this relative rotation convention; it is duplicated and corrupted in the local TeX.

For a very large head, the coefficient regime is $A_0\gg A$ and $D_0\gg D$. ConsequentlyThe large translational resistance $A_0$ makes the microswimmer slow even though the flagellum rotates almost at the motor rate.

For a very small head, $A_0\ll A$ and $D_0\ll D$. The appropriate denominator retains the translation-rotation coupling:Both tend to zero with $D_0\propto a^3$, while ${\Omega }_0\simeq -\omega$. The motor mainly rotates the low-resistance head, and produces little flagellar motion relative to the fluid.

In the intermediate coefficient regime $A_0\ll A$ and $D_0\gg D$, the head supplies a strong rotational reaction with little translational penalty. ThenThe head radius has dropped out. To find the optimal pitch of a helical microswimmer, put $z=\tan \varphi$. The dimensionless speed becomes $z/(z^2+2)$ and has derivative $(2-z^2)/(z^2+2{)}^2$. Hence, for $\omega >0$ with the chosen handedness,The angle is measured from the horizontal plane. An intermediate pitch combines axial and transverse tangent directions, allowing anisotropic drag to convert rotation into translation; either extreme removes that coupling.

The geometric head-size regimes quoted in the paper hold with the logarithmic slenderness factor treated as fixed. More precisely, $A_0\ll A$ requires $a\ll L(1+{\cos }^2\varphi )/(3|\log ϵ|)$, while $D_0\gg D$ requires $a^3\gg L{b}^2(1+{\sin }^2\varphi )/(4|\log ϵ|)$. These coefficient inequalities state the validity of the intermediate approximation when the logarithm is quantitatively important. The many-turn and local-drag assumptions of slender-body theory also remain in force.

Integrating the continuity equation over state space and applying the divergence theorem givesAssume the probability density function and its derivatives decay sufficiently rapidly at infinity, or impose a zero-flux boundary condition. The boundary integral then vanishes, so $\int PdV$ is constant. Since it is initially one, the distribution remains normalized for all time.

For each incompressible Stokes flow without body force, ${\partial }_j{\sigma }_{ij}=0$ andThe divergence theorem givesIncompressibility removes the pressure term, antisymmetric velocity-gradient parts do not contract with the symmetric stress, and the remaining integrand is $2\mu e_{ij}^{(1)}{e}_{ij}^{(2)}$, symmetric under $1\leftrightarrow 2$. This proves the Lorentz reciprocal theorem for Stokes flow:

By Linearity of Stokes flow, the drag on a translating body is linear in its velocity, so $F_i=A_{ij}{U}_j$. Apply the reciprocal theorem to two translation solutions with velocities ${\mathbf{U}}^{(1)}$ and ${\mathbf{U}}^{(2)}$. Their surface velocities are constant and the far-field contribution vanishes, yieldingSince the two vectors are arbitrary,Thus the hydrodynamic resistance matrix is symmetric and, for fixed fluid viscosity, depends only on body geometry.

A straight rod has fore-aft symmetry about its centre. Gravity acts through that centre, so the translational-rotational coupling and applied torque vanish. Kinematic reversibility of Stokes flow and Uniqueness of Stokes flow then rule out spontaneous rotation: the inclination $\theta$ remains constant.

Let the speed for perpendicular translation be $U_0$. The given factor of two makes the parallel speed $2U_0$ for the same force component. Resolving gravity parallel and perpendicular to the rod givesResolving these velocities horizontally and vertically givesTherefore

Using the Levi-Civita symbol and the epsilon-delta identity,ThereforeOne particularly simple vector potential for $u$ isDirect differentiation gives $\nabla \cdot A=0$ and $\nabla \times A=(y^2,z^2,x^2)=u$.

Parametrize the curved cone byThe outward oriented surface element isSubstitution and integration over $0\le z\le h$, $0\le \phi <2\pi$ give

Because $\nabla \cdot u=0$, the divergence theorem predicts that the two outward fluxes sum to zero. On the top disk $S_2$, $d\mathbf{S}=(0,0,1)dxdy$ and $u_z=x^2$, so, with $R=h\tan \alpha$,as predicted.