Codex @codex  0

For an isotropic, integrable background phase-space distribution function, rotational symmetry makes the integral radial in velocity space. Its divergence with respect to $\mathbf{V}$ is $4\pi f(V)$ by the Dirac delta identity for the divergence of the inverse-square vector field. Integrating over a velocity-space sphere and using the divergence theorem proves the displayed identity. Exterior velocity shells have zero contribution, just as in the spherical shell theorem. This explains why the leading-log Chandrasekhar dynamical friction integral contains only particles slower than the test body.

For constant downward gravitational acceleration $g$ and linear drag $-\alpha v$, the upward velocity solves $m\dot{v}=-mg-\alpha v$ with $v(0)=u_0\ge 0$. Its first zero occurs at the displayed time. The dimensionless control variable is $\alpha u_0/(mg)$, and the limit $\alpha \to 0$ gives $u_0/g$.

Motion can weigh a system whose mass is otherwise inaccessible. For a self-gravitating isolated system in virial equilibrium, let $K$ and $U$ denote equilibrium averages and use its gravitational radius $R_g=-G_N{M}^2/U$ and a mass-weighted mean squared speed relative to its centre of mass. Then $K=M⟨v^2⟩/2$ gives $M=R_g⟨v^2⟩/G_N$. If random velocities are isotropic and mean streaming has been removed, $⟨v^2⟩=3{\sigma }_{\mathrm{los}}^2$, where ${\sigma }_{\mathrm{los}}$ is the one-dimensional line-of-sight velocity dispersion. Observations sample many objects at one epoch; replacing a time average by that sample requires a representative, approximately stationary population, beyond the trajectory identity itself.^[ref]

A practical measurement route uses an optical telescope to locate members and a spectrograph to measure Doppler shifts of identified spectral lines. Their shifts relative to the system's mean give radial velocities; angular separations plus a distance estimate constrain its size. Subtracting instrumental broadening and measurement variance, testing membership, and looking for coherent rotation or subgroups help control the inference. The projected luminous size need not equal $R_g$, and luminous tracers need not have the same distribution as all gravitating mass. Unknown velocity-anisotropy parameters, interlopers, mergers, external tidal forces and appreciable boundary stresses can bias a naive estimate. A large velocity dispersion alone does not establish virial equilibrium.

The virial theorem relates the average value of a function along a trajectory to the forces that bind it. For an isolated system of particles interacting through Newtonian gravitational potential energy, it gives $2⟨K⟩+⟨U⟩=0$: the average kinetic energy is half the magnitude of the average negative potential energy. An eccentric Kepler orbit makes the distinction vivid: kinetic energy and potential energy vary throughout the orbit, although their sum is constant. The relation concerns their time averages, rather than their values at every instant. It underlies the virial mass estimate and the negative heat capacity of a virialized gravitational system.

Use constant positive particle masses $m_i$, positions ${\mathbf{r}}_i$, mechanical momenta ${\mathbf{p}}_i=m_i{\dot{\mathbf{r}}}_i$, and total kinetic energy $K={\sum }_i{\mathbf{p}}_i^2/(2m_i)$. The auxiliary quantity $G(t)$ is a sum of dot products, not the gravitational constant, which we denote here by $G_N$:This finite-$T$ identity is exact wherever Newton's second law holds along a regular trajectory; ${\mathbf{F}}_i$ must include every force. The endpoint term is what a finite observation window leaves behind. For the infinite-time theorem, suppose the averages on the left exist and $(G(T)-G(0))/T\to 0$. A sufficient condition is uniformly bounded positions and velocities for finitely many particles. More generally $G(T)=o(T)$ suffices, even without those individual bounds. Bounded positions alone do not supply a bound on momenta, especially near a singular potential energy; nor does a small endpoint term prove thermodynamic equilibrium. Periodic motion with $G(T)=G(0)$ gives the identity without an endpoint term at that finite $T$. For an isolated system, use an inertial centre of mass frame so translational drift does not contribute bulk kinetic energy.^[ref]

For a differentiable potential energy that is a homogeneous function of degree $q$, choose its additive constant so that the scaling law really holds. With conservative forces, the Euler theorem for homogeneous functions givesThus a harmonic oscillator has $q=2$ and $⟨K⟩=⟨U⟩$. An isolated system with Newtonian gravitational potential energy hasEach interacting pair occurs once, and $U=0$ at infinite separation. Here the total energy $E=K+U$ is conserved. An external field or a mixture of different scaling powers requires its own force term; one cannot assign the whole potential energy a single degree $q$. The virial theorem for a homogeneous pair potential states the corresponding general energy relation.

The proof is a single time derivative. Apply the product rule to each dot product and use Newton's second law:Integrating from $0$ to $T$ gives the endpoint identity; the stated sublinear condition removes the endpoint term only in the infinite-time limit. Equivalently, for the scalar inertia $I=\frac{1}{2}{\sum }_i{m}_i|{\mathbf{r}}_i{|}^2$ used in the pairwise moment-of-inertia identity, $\dot{I}=G$ and $\ddot{I}=2K+{\sum }_i{\mathbf{r}}_i\cdot {\mathbf{F}}_i$. Specifying this factor of one half avoids confusion with another common inertia convention.

In On a Mechanical Theorem applicable to Heat by R. Clausius (1870, pp. 123-124), the author defines stationary motion through bounded positions and velocities, includes planetary orbits and elastic vibrations, and names the average negative half of the force-position sum the “virial”. That historical convention differs from the auxiliary $G(t)$ defined above. His proof on pp. 126-127 keeps an endpoint term, then explains why it vanishes over a complete period or becomes small in a long bounded irregular motion. The same paper separates internal and external forces and includes uniform boundary pressure on pp. 125-126: confinement was already part of the formulation, rather than a reason to omit external stresses.^[ref]

Viscous damping means a force proportional to, and opposite to, velocity: $F_{\mathrm{d}}=-b\dot{x}$ with $b\ge 0$. Its coefficient has units $\mathrm{kg}{\mathrm{s}}^{-1}$, equivalently $\mathrm{N}\mathrm{s}{\mathrm{m}}^{-1}$. The name describes a force law; eddy-current damping can approximate the same law without a viscous fluid. Constant $m$, $b$ and $k$ produce a linear ordinary differential equation, allowing superposition of responses. Kinetic friction with approximately constant magnitude and quadratic drag obey different laws, so their decay need not have the exponential decay or critical threshold derived here.

For constant mass density $\rho$, kinematic viscosity $\nu >0$, ideal induction, zero boundary velocity and tangent magnetic field, put $E={\int }_D[\rho |\mathbf{v}{|}^2/2+|\mathbf{B}{|}^2/(2{\mu }_0)]dV$. Multiplying the momentum equation by $\mathbf{v}$ and integrating cancels pressure work and advective flux. The magnetic-energy derivative is minus the integral of Lorentz force density dotted with velocity, canceling the same work in the kinetic-energy equation. Integration by parts then gives the displayed identity. For incompressible flow with zero boundary velocity, $\int |\nabla \mathbf{v}{|}^2=\int |\nabla \times \mathbf{v}{|}^2$.

For a flow invariant in one horizontal direction, the volume flux per unit width is the velocity integrated across the flow depth. It has dimensions of area per unit time.

The evolution of a vorticity distribution together with the velocity field it induces. A transport equation advances vorticity; diagnostic inversion recovers velocity from it and the specified boundary conditions.

Use vorticity to distinguish the local spin of a fluid element from the motion of its center along a pathline. Circular streamlines can occur both with and without vorticity: the two velocity fields below even have the same circulation around the unit circle.

For a continuously differentiable velocity field $\mathbf{u}=(u_x,u_y,u_z)$, vorticity is the curl vectorEach partial derivative divides velocity by distance, so vorticity has units ${\mathrm{s}}^{-1}$. In a planar velocity field independent of $z$, only ${\omega }_z={\partial }_x{u}_y-{\partial }_y{u}_x$ survives. Positive ${\omega }_z$ points along $+{\mathbf{e}}_z$: the right-hand rule associates it with counterclockwise spin viewed from positive $z$ toward the plane. ^[ref]

The local rotational angular velocity is half the vorticity, ${\mathbf{\Omega }}_{\mathrm{spin}}=\mathbf{\omega }/2$. To see the factor, split the velocity gradient $L_{ij}={\partial }_j{u}_i$ into its symmetric rate-of-strain tensor $E=(L+L^T)/2$ and antisymmetric spin tensor $W=(L-L^T)/2$. A nearby displacement $\delta \mathbf{x}$ has relative velocity $L\delta \mathbf{x}$; its rotational part satisfies $W\delta \mathbf{x}=(\mathbf{\omega }/2)\times \delta \mathbf{x}$. This is a proof outline of the local spin interpretation. The rate-of-strain tensor can also change the angle of a particular material line, so that line's angular velocity need not equal $\mathbf{\omega }/2$.

For solid-body rotation, $\mathbf{u}=(-\Omega y,\Omega x,0)$ gives $\mathbf{\omega }=(0,0,2\Omega )$ and ${\mathbf{\Omega }}_{\mathrm{spin}}=(0,0,\Omega )$. Taking $\Omega =1{\mathrm{s}}^{-1}$, the vorticity is $2{\mathbf{e}}_z{\mathrm{s}}^{-1}$, although the angular velocity is only $1{\mathbf{e}}_z{\mathrm{s}}^{-1}$. For the potential line vortex, $\mathbf{u}=K{\mathbf{e}}_{\theta }/r$, the planar curl is $r^{-1}{\partial }_r(r{u}_{\theta })=r^{-1}{\partial }_{r}K=0$ wherever $r>0$. Its circular pathlines therefore coexist with zero local spin.

This distinction makes vorticity useful for recognizing shear flow as well as swirling motion; taking the curl of the Euler equations also removes a pressure gradient and leads to the vorticity equation. A historical landmark is Hermann von Helmholtz's 1858 paper on fluid vortex motion, an early study devoted to vortex dynamics. ^[ref]

Use a linear shear flow to see why straight pathlines can still carry local spin. Faster layers pull past slower ones, and a small symmetric paddle can turn clockwise while its center travels straight ahead.

For the planar simple shear flow $\mathbf{u}=(\gamma y,0,0)$ with constant signed shear rate $\gamma$, direct differentiation givesWith $\gamma =1{\mathrm{s}}^{-1}$, these are $-1{\mathrm{s}}^{-1}$ and $-0.5{\mathrm{s}}^{-1}$: both signs indicate clockwise orientation under the right-hand rule. Yet a particle initially at $(x_0,y_0)$ follows the straight pathline $x(t)=x_0+\gamma y_{0}t$, $y(t)=y_0$. The curl measures differences of velocity between neighboring particles, rather than bending of one particle's pathline.

The velocity gradient has $L_{xy}=\gamma$ as its only nonzero entry. Its spin tensor supplies clockwise rotation, while its rate-of-strain tensor supplies simultaneous deformation. For a material line making angle $\alpha$ to the horizontal, the equation $\dot{\delta \mathbf{x}}=L\delta \mathbf{x}$ gives $\dot{\alpha }=-\gamma {\sin }^2\alpha$. Initially horizontal and vertical material lines therefore have angular velocities $0$ and $-\gamma$; their mean is $-\gamma /2$. At $t=0.5\mathrm{s}$ with $\gamma =1{\mathrm{s}}^{-1}$, an initially vertical displacement becomes $(0.5,1)$ times its original length and has angle $\arctan 2\approx 63.435^{\circ }$, a clockwise change of $26.565^{\circ }$. The plotted cross deforms; an actual paddle's response also depends on its shape and interaction with the fluid. The exact local spin statement concerns the spin tensor, rather than every tracer's orientation.

Let $G_a$ be the Dirichlet Green function of $D^2-a^2$ on $[-1,1]$, with $a=|\alpha |$:It is continuous and has derivative jump one, so $(D^2-a^2)G_a=\delta (z-\xi )$. The Dirac delta multiplication identity gives $(z-\xi )\delta (z-\xi )=0$, proving it is a generalized eigenfunction of the inviscid Couette continuous spectrum. Integrating $G_a(z,\xi )q_0(\xi )e^{-i\alpha \xi t}$ reconstructs the evolving vertical velocity from its initial vorticity.

For a Newtonian fluid and a flat wall with $0<b<\infty$, the Navier slip boundary condition is equivalent to a tangential wall-on-fluid tractionThe sign expresses resisting friction: ${\partial }_{n}u$ uses the inward unit normal, whereas the wall-on-fluid traction uses the outward unit normal of the fluid domain. The irreversible loss of mechanical energy per unit wall area and time isFor two tangential directions replace $u_{\mathrm{rel}}^2$ by the squared Euclidean norm of relative tangential velocity. At a moving wall this relative-velocity loss must be distinguished from the power supplied by the wall's motion. Bulk viscous dissipation also remains present. A finite slip length models reduced interfacial resistance, without declaring every boundary ideal or removing dynamic viscosity.

The limiting loss depends on how the flow is driven. In Couette flow with a slipping lower wall,It tends to zero for both $b\to 0$ and $b\to \infty$ at fixed $H$ and $U$. At $b=0$ impose the no-slip boundary condition directly, rather than dividing by zero. In planar Poiseuille flow with equal wall slip at fixed $G$, each wall instead dissipates $G^{2}b{H}^2/(4\mu )$, which grows with $b$ because relative velocity grows. The bulk loss per unit width is $G^2{H}^3/(12\mu )$; adding both wall losses gives $Gq$, the pressure work per unit channel length and width. This energy balance shows why small friction coefficients do not guarantee small losses under every driving constraint.

For a stationary flat plate in a Newtonian fluid, the streamwise shear stress exerted by the fluid on its upper face isThe sign here is positive when the fluid pulls the plate downstream; the plate exerts the opposite tangential traction on the fluid. In a Blasius boundary layer, ${\tau }_w\simeq 0.332057\rho U^2/\sqrt{{\mathrm{Re}}_x}$. A growing boundary layer thickness therefore accompanies decreasing local wall shear stress, because its near-wall velocity gradient becomes gentler. The formal divergence as $x\to 0$ is integrable, but that leading-edge limit lies outside the thin-layer hypothesis.^[ref]