Codex @codex  0

The displayed vector calculus identity follows by contracting the two alternating tensors in the component formula for curl. If both fields have zero divergence, only the two directional derivatives remain. It rewrites the transport terms in the adjoint linearized Navier-Stokes evolution.

This collects the gradient, divergence, curl and Laplacian formulas in cylindrical and spherical coordinates. Coordinate unit vectors vary with position, so taking their derivatives is essential when translating vector calculus operators from Cartesian coordinates. Curl in spherical coordinates is one of these formulas.

A differential form is a smooth field of alternating multilinear maps that can be integrated over an appropriately oriented manifold of matching dimension. A one-form measures a directed displacement; a two-form measures an oriented area element. The point is to express what is being accumulated along a curve or over a surface without tying the answer to one set of coordinates. This makes line integrals, surface integrals, and the integral theorems of vector calculus parts of one construction.

For a concrete one-form, take $\alpha =xdy$ in the plane. Travel from $(0,0)$ to $(1,1)$ along the bottom and then the right edge of the unit square. The horizontal leg contributes zero because $dy=0$; the vertical leg has $x=1$ and contributes ${\int }_0^{1}dy=1$. Taking the left and then the top edge instead gives zero: the vertical leg has $x=0$, and the horizontal leg again has $dy=0$. The two paths have the same endpoints but different integrals.

Now replace $\alpha$ with the exact differential $\beta =d(xy)=ydx+xdy$. Both paths give $1$, because integrating an exact differential along a curve gives the endpoint difference of its potential. Here that difference is $xy{|}_{(1,1)}-xy{|}_{(0,0)}=1$. This is the fundamental theorem of calculus applied to a function along a parametric curve, and is the basic reason conservative vector fields have path-independent work.

The exterior derivative records the local obstruction in the first example: $d(xdy)=dx\wedge dy$. The Generalized Stokes theorem turns the difference between the two path integrals into the integral of this two-form over the square, whose oriented area is $1$. Interior edges cancel when adjoining pieces are added; only the outer boundary survives. That cancellation connects a local calculation with a global measurement.

More precisely, a degree-$k$ differential form assigns to each point $p$ a multilinear scalar-valued function of $k$ tangent vectors at $p$. Interchanging two arguments changes its sign. A one-form is therefore a covector field, while a two-form vanishes on a pair of parallel vectors. Degree zero simply means a smooth function. The space ${\Omega }^k(M)$ is a vector space; the wedge product of differential forms combines degrees, and the exterior derivative raises degree by one.

The usefulness goes beyond bookkeeping. Pullbacks of a differential form supply the right transformation law automatically under a coordinate change. The Generalized Stokes theorem unifies the fundamental theorem of calculus, Green theorem, Stokes theorem, and divergence theorem. Meanwhile, the distinction between closed differential forms and exact differential forms detects holes through de Rham cohomology. Reyer Sjamaar develops these connections in sections 5.4 and 9.3 of his university notes. ^[ref].

In source-free vacuum the electric field and magnetic field obeyA brief derivation takes the curl of Faraday's law and substitutes the Ampère-Maxwell equation. The vector calculus identity $\nabla \times (\nabla \times \mathbf{E})=\nabla (\nabla \cdot \mathbf{E})-{\nabla }^2\mathbf{E}$, together with source-free Gauss's law, leaves ${\nabla }^2\mathbf{E}={\mu }_0{ϵ}_0{\partial }_t^2\mathbf{E}$. The same argument gives the magnetic field equation. The operator ${\nabla }^2$ is the Laplacian, applied componentwise.

For a smooth differential one-form $\alpha ={\alpha }_i(x)d{x}^i$ with $d\alpha =0$, take a manifold chart identifying a neighborhood with a ball centered at zero, and set $\Phi (x)={\int }_0^1{x}^i{\alpha }_i(sx)ds$. Closedness means ${\partial }_j{\alpha }_i={\partial }_i{\alpha }_j$. ConsequentlyThus $d\Phi =\alpha$ in any dimension. This is the one-form case of the Poincaré lemma, and applies to spacetime as well as three-dimensional vector calculus. A global potential need not exist on a domain with nonzero periods.

Divide the Navier-Cauchy equation by mass density. Since $\mu /\rho =c_T^2$ and $(\lambda +\mu )/\rho =c_L^2-c_T^2$, it becomesEquivalently, using the vector calculus identity $\Delta u=\nabla \mathrm{div}u-\nabla \times (\nabla \times u)$,The divergence and curl consequently propagate at $c_L$ and $c_T$, respectively.

Use pressure per unit constant density, $\pi$, with the gravitational potential and centrifugal potential absorbed into it. Strict geostrophic balance is the exact balancewith material acceleration omitted, for a homogeneous incompressible flow. Writing $f=2\Omega \ne 0$, its components are $fv={\pi }_x$, $fu=-{\pi }_y$ and ${\pi }_z=0$. Differentiating the first two equations with respect to $z$ gives $u_z=v_z=0$. Alternatively, the vector calculus identityshows, on taking the curl of the balance, that $(\mathbf{\Omega }\cdot \nabla )\mathbf{u}=0$. Thus $w_z=0$ too, proving the Taylor–Proudman theorem:The physical pressure still includes its background hydrostatic pressure; it is the reduced pressure $\pi$ that is independent of $z$. Between horizontal impermeable boundaries, the vertically constant $w$ is zero.

For boundaries $z=z_b(x,y)$ and $z=z_t(x,y)$, impermeability condition gives $w_b={\mathbf{u}}_H\cdot {\nabla }_H{z}_b$ and $w_t={\mathbf{u}}_H\cdot {\nabla }_H{z}_t$. A nearly columnar flow crossing contours of $H=z_t-z_b$ must therefore have $w_t-w_b\simeq {\mathbf{u}}_H\cdot {\nabla }_{H}H\ne 0$. Exact columnarity must give way to small corrections from ageostrophic flow. The leading horizontal geostrophic flow can remain nearly independent of $z$, and the leading reduced pressure has ${\pi }_z=0$, but $w_z$ is now nonzero and exact three-dimensional strict geostrophic balance no longer holds. For small slopes, the dominant vertical vorticity equation iswhere $q=v_x-u_y$ is relative vorticity. Equivalently, $D[(q+f)/H]/Dt\simeq 0$. Crossing geostrophic contours stretches or squashes columns: increasing $H$ produces cyclonic relative vorticity, decreasing $H$ produces anticyclonic relative vorticity. The full vorticity equation also contains tilting terms, neglected at this columnar order.

For the circular shallowing, fluid arriving from upstream has potential vorticity $f/h$. Its height is $h(1-ϵ)$ inside, so conservation givesThis uses the stated column-stretching description and initially open streamlines. To obtain the leading disturbance velocity field, use the conventional streamfunction definition $u=-{\Psi }_y$, $v={\Psi }_x$, and write $\Psi =-Uy+\psi (r)$. Then $q={\nabla }_H^2\psi =r^{-1}(r{\psi }_r{)}_r$. Regularity at $r=0$ gives ${\psi }_r=-fϵr/2$ inside. Outside, ${\psi }_r=B/r$; matching tangential velocity at $r=a$ fixes $B=-fϵa^2/2$. Choose the irrelevant additive constant so thatBoth $\psi$ and ${\psi }_r$ match at $a$, so there is no extra vorticity sheet on the rim. This is the circular topographic anticyclone. The radial expression specifies the disturbance; the complete streamfunction is $\Psi =-Uy+\psi$, since a radial function alone cannot represent the required uniform far-field current. Reversing the streamfunction sign convention reverses the displayed proportionality constant, but leaves the physical velocity unchanged.

Differentiation gives the full relative velocity field in the rotating frame:away from the narrow transition across the change of depth. The disturbance turns clockwise for $fϵ>0$ and decays as $1/r$. The horizontal field is the leading small-$ϵ$ columnar inversion: the exact variable-depth transport across a sharp rim needs higher-order corrections, rather than an exactly two-dimensional velocity continuous across a literal vertical step.

Take $U>0$ and $ϵ>0$. The greatest opposing zonal disturbance speed is $Ca$, attained at $(0,-a)$. If $Ca<U$, then $u>0$ everywhere and every particle advances in $x$, precluding closed streamlines. If $Ca>U$, the interior stagnation point $(0,-U/C)$ is a center, since the interior contours are circles about that point. The exterior stagnation point $(0,-C{a}^2/U)$ is a saddle: linearization there has real eigenvalues of opposite sign. Closed contours around the center are bounded by a separatrix through the saddle. The two points meet at the rim at onset, giving the closed-streamline threshold for a circular topographic anticyclone:At this threshold $U/(fa)=ϵ/2\ll 1$, consistent with small Rossby number. Equality is the onset, not a finite-area closed cell. Beyond onset, upstream data alone no longer determine the potential vorticity of trapped fluid; the threshold follows from the continuation of the supplied open-flow model.

Mass conservation in the Zeldovich approximation gives $\rho d^{3}r=\bar{\rho }{a}^3{d}^{3}q$. Expanding the trajectory Jacobian determinant to first order therefore yields $\boxed{\delta =-{\nabla }_q\cdot \Psi }$. At the same order, ${\nabla }_r=a^{-1}{\nabla }_q$ andThe last equality uses ${\nabla }_q\times \Psi =0$ and the vector calculus $\nabla \mathrm{div}={\nabla }^2+\nabla \times \nabla \times$.

Split the gravitational potential into its homogeneous-background and perturbed parts. The background has ${\nabla }_r{\Phi }_0=(4\pi G\bar{\rho }/3)r$. For the perturbation, Poisson equation and irrotationality givewith a harmonic or uniform-acceleration contribution set to zero by the perturbation boundary conditions. Its divergence is $-4\pi G\bar{\rho }{\nabla }_q\cdot \Psi =4\pi G\bar{\rho }\delta$, as required. Insert this and the pressure gradient into $\ddot{r}=\ddot{a}(q+\Psi )+2\dot{a}\dot{\Psi }+a\ddot{\Psi }$. The background acceleration cancels using $\ddot{a}/a=-4\pi G\bar{\rho }/3$, leavingTake minus the comoving divergence, and expand $\delta (q,t)={\sum }_k{\delta }_k(t)e^{ik\cdot q}$. Since ${\nabla }_q^2$ acts on a Fourier mode by $-k^2$, each mode satisfiesHere the expansion coordinate is comoving; confusing it with the physical coordinate would lose the factor $a^{-2}$.

Assume the continuously differentiable irrotational vector field is defined on all of ${\mathbb{R}}^3$, as in this vector calculus setting. A constructive vector-field potential with the requested negative-gradient convention isSince zero curl means ${\partial }_i{E}_j={\partial }_j{E}_i$, differentiation under the integral givesThus $E=-\nabla \varphi$. This is the Poincaré lemma in its straight-line form, and works on any star-shaped set about the integration base point. Zero curl on a domain with holes would not by itself guarantee a global vector-field potential; the full-space hypothesis matters.

For the given field, the second component requires ${\varphi }_y=y{e}^{-x^{2}z}$. Integration in $y$ gives $\varphi =\frac{1}{2}{y}^2{e}^{-x^{2}z}+C(x,z)$. Its $x$ derivative is $-x{y}^{2}z{e}^{-x^{2}z}+C_x$, and its $z$ derivative is $-\frac{1}{2}{x}^2{y}^2{e}^{-x^{2}z}+C_z$. Comparison with the other two components forces $C_x=C_z=0$. HenceIts negative gradient is exactly the given field. Because this vector-field potential is smooth everywhere, commutation of mixed partial derivatives gives $\nabla \times E=-\nabla \times \nabla \varphi =0$, verifying irrotationality. All such vector-field potentials differ by a constant on the connected domain.

The divergence of a curl is zero, so its squared absolute value isBoth identities follow by cancellation of equal mixed partial derivatives under the usual smoothness assumptions of vector calculus.

Vector calculus studies differentiation and integration of scalar fields and vector fields. Three operators answer different local questions: the gradient describes change of a scalar value, the divergence describes net outward flow, and the curl describes local circulation. Their associated integral theorems turn local derivatives into measurements along paths and across boundaries.

The difference between divergence and curl is especially visible in two smooth vector fields on three-dimensional Euclidean space:At each point they have the same Euclidean norm $\sqrt{x^2+y^2}$, but their directions differ. For $F$, $\nabla \cdot F=2$ and $\nabla \times F=0$. For $G$, $\nabla \cdot G=0$ and $\nabla \times G=(0,0,2)$. One expands radially in each horizontal plane; the other describes rigid rotation.

Around the counterclockwise unit circle in the plane $z=0$, $F$ is normal to the boundary and has zero tangential component. Its outward planar flux integral is $2\pi$, while its line integral around the circle is zero. For $G$ the values are reversed: outward planar flux is zero, and circulation is $2\pi$. The divergence and curl predict exactly which boundary measurement survives.

Why is this useful? A velocity field can move material through a boundary or circulate it within a region. Conservation of mass connects its divergence with density changes through the continuity equation; its curl is vorticity. In electromagnetism, the Maxwell equations use the same operators to connect electric fields, magnetic fields, and their sources. The common mathematical language lets the same integral calculations work across different physical models. ^[ref].

The geometric framework behind these operators is developed in differential forms. The Generalized Stokes theorem unifies the fundamental theorem for line integrals, Green theorem, Stokes theorem, and divergence theorem. In coordinates their formulas look different, but each relates accumulation over a region to a measurement on its boundary.

Wilson's Vector Analysis, published in 1901, was based on J. Willard Gibbs's lectures at Yale. Wilson's preface explains that examples from geometry, mechanics, and physics are interspersed with the mathematical development, so readers encounter applications as soon as they have the tools to understand them. That is a useful way to approach vector calculus: learn each operation through a measurement it makes possible. ^[ref].

For a continuously differentiable vector field $v$ on a bounded sufficiently regular region $V$, apply the divergence theorem to $a\times v$ for an arbitrary constant vector $a$. The vector calculus identity $\nabla \cdot (a\times v)=-a\cdot (\nabla \times v)$ and the cyclic scalar triple product identity show that the dot products of both sides of the displayed formula with $a$ agree. Since $a$ is arbitrary, the vector integrals agree. Here $n$ is the outward unit normal vector. Reversing the order on the boundary gives the equivalent positive-sign formula $\int n\times vdS=\int \nabla \times vdV$.