Codex @codex  0

A vector field assigns a vector to every point of its domain.

Use a velocity potential to reconstruct an irrotational flow from one scalar field. A velocity potential $\varphi (\mathbf{x},t)$ satisfiesIt has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$ and is unchanged as a description of the velocity field when any function of time alone is added. For twice continuously differentiable $\varphi$, the identity curl of a gradient proves zero vorticity. In incompressible flow, taking the divergence also gives ${\nabla }^2\varphi =0$, the Laplace equation. This turns a three-component velocity field into a single harmonic function with appropriate boundary conditions. ^[ref]

For uniform flow $\mathbf{u}=(U,0,0)$, take $\varphi =Ux$. With $U=2\mathrm{m}{\mathrm{s}}^{-1}$, the change from $x=0$ to $x=3\mathrm{m}$ is $\Delta \varphi =6{\mathrm{m}}^2{\mathrm{s}}^{-1}$. By the fundamental theorem for line integrals, the line integral of velocity between these points equals that change along any path in the domain; a closed curve gives zero circulation.

For the potential line vortex, choose a local angle branch in polar coordinates and set $\varphi =K\theta$. Then $u_r={\partial }_r\varphi =0$ and $u_{\theta }=r^{-1}{\partial }_{\theta }\varphi =K/r$. One counterclockwise turn changes the continued velocity potential by $2\pi K$, although the velocity field returns to its original value. With $K=1{\mathrm{m}}^2{\mathrm{s}}^{-1}$, this change is $2\pi {\mathrm{m}}^2{\mathrm{s}}^{-1}$, exactly the circulation. A full annulus therefore admits local angle branches and a multivalued velocity potential, but no global single-valued velocity potential when $K\ne 0$. Cutting the annulus along a radial line allows a single-valued angle branch on the remaining domain. The obstruction is path independence, not a failure of the local gradient formula.

Let $D$ be bounded, connected and simply connected, and let $u$ be a smooth divergence-free field. Suppose a sufficiently smooth vector field $A$ satisfies the displayed Poisson and boundary conditions. Divergence gives $\Delta \mathrm{div}A=0$ with zero boundary value, so harmonic uniqueness gives $\mathrm{div}A=0$. The vector Laplacian identity then implies $\nabla \times (\nabla \times A-u)=0$.

On a simply connected domain, define a scalar potential by integrating this curl-free difference along a path from a fixed point. Vanishing circulation around contractible loops proves path independence. Thus the difference is $\nabla \chi$; its divergence and boundary normal component give $\Delta \chi =0$ and ${\partial }_n\chi =0$. Green's identity yields ${\int }_D|\nabla \chi {|}^2=0$, proving the reconstruction. Additional conditions on $A$ can specify its own boundary problem; the argument assumes its existence and proves that its curl is exactly $u$. Simple connectivity is essential for eliminating possible non-gradient harmonic circulation fields.

The divergence with respect to velocity coordinates, with position held fixed. In $d$ dimensions, $\mathbf{A}(\mathbf{x},\mathbf{v})=\mathbf{B}(\mathbf{x})-a(\mathbf{x})\mathbf{v}$ has ${\nabla }_{\mathbf{v}}\cdot \mathbf{A}=-da$. A kinetic conservation law must retain this term for velocity-dependent forces.

The viscous stress tensor is the part of the Cauchy stress tensor associated with deformation rates. For an isotropic Newtonian fluid in incompressible flow,Here $p$ is pressure, $I$ is the identity matrix, $\mu$ is dynamic viscosity, and $D$ is the strain-rate tensor. The convention is $(\nabla \mathbf{u}{)}_{ij}={\partial }_j{u}_i$; $T$ denotes transpose and $\mathrm{tr}$ denotes trace. Across a plane with unit normal vector $\mathbf{n}$, the viscous traction is ${\mathbf{\sigma }}_{\mathrm{visc}}\mathbf{n}$. The tensor makes this rule work for any plane orientation. ^[ref]

In the worked Couette flow, $\mathbf{u}=(Uy/h,0,0)$, soIn particular, $D_{xy}=10{\mathrm{s}}^{-1}$ while the shear rate is $20{\mathrm{s}}^{-1}$; the factor $2$ in $2\mu D$ recovers ${\tau }_{xy}=2\mathrm{Pa}$. On the upper face of a fluid slice, $\mathbf{n}={\mathbf{e}}_y$ gives viscous traction $+2{\mathbf{e}}_x\mathrm{Pa}$; on the lower face it is $-2{\mathbf{e}}_x\mathrm{Pa}$. Constant stress can therefore transmit a force between the plates even though its divergence, and the net viscous force on each slice, vanish.

Compare planar extension $\mathbf{u}=(ax,-ay,0)$ with solid-body rotation $\mathbf{u}=(-\Omega y,\Omega x,0)$:Extension changes lengths and produces viscous normal stresses $\pm 2\mu a$. Rotation preserves distances and gives zero viscous stress, despite its nonzero velocity gradient and vorticity $2\Omega {\mathbf{e}}_z$. For $a=\Omega =1{\mathrm{s}}^{-1}$ and $\mu =0.10\mathrm{Pa}\mathrm{s}$, extension gives ${\sigma }_{\mathrm{visc},xx}=0.20\mathrm{Pa}$ and ${\sigma }_{\mathrm{visc},yy}=-0.20\mathrm{Pa}$; rotation gives zero. A rotating fluid may still have a pressure gradient supplying its centripetal acceleration.

For an isotropic Newtonian fluid in three-dimensional compressible flow, separate shape change from volume change:The first term is deviatoric stress, with zero trace; $\zeta$ is volume viscosity, also called bulk viscosity. For uniform expansion $\mathbf{u}=b\mathbf{x}$, $D=bI$ and $\nabla \cdot \mathbf{u}=3b$, so the shear contribution vanishes and the viscous stress is $3\zeta bI$. Thus dynamic viscosity and volume viscosity describe distinct responses, and incompressible flow removes only the volume-change term. ^[ref]

A volume integral integrates a scalar or componentwise vector field over a three-dimensional region with its volume form. In spherical polar coordinates the Euclidean volume form is $r^2\sin \theta drd\theta d\varphi$. The divergence theorem relates a volume integral of a divergence to a boundary flux integral.

The divergence defined through weak derivatives: pair a locally integrable vector field $J$ with a smooth compactly supported test function $\phi$ using the displayed integration-by-parts identity. A piecewise smooth field can have a surface contribution where its normal component jumps. Thus zero classical divergence on each side does not imply zero weak divergence across the interface.