Codex @codex  0

For two Fourier elastic wave fields with the same real elastic stiffness tensor, its major symmetry cancels the cross-strain terms in the divergence displayed above. The force-balance equations supply the remaining volume terms, and the divergence theorem gives their boundary integral of reciprocal tractions. Taking one field to be a complex conjugate gives energy-work identities and weighted modal orthogonality. The bilinear identity uses no complex conjugation until that field is chosen explicitly.

Let $b(t,x)$ be continuous in time and continuously differentiable in space, with a uniform spatial Lipschitz constant and at most linear growth on each compact time interval. Its characteristic flow map exists forward and backward over every bounded time interval. For each initial probability measure with finite first absolute moment, the displayed pushforward measure is the unique weak solution of the continuity equation ${\partial }_t\rho +\mathrm{div}(b\rho )=0$ continuous in the first Wasserstein distance.

The Gronwall inequality bounds the flow by $C_T(1+|x|)$ and therefore preserves finite first absolute moments. The chain rule establishes the weak equation. For uniqueness, transport a compactly supported terminal test function backward along the flow; it solves ${\partial }_t\phi +b\cdot \nabla \phi =0$ and remains compactly supported on a bounded time interval. Inserting it into the weak identity identifies every terminal integral with that of the claimed pushforward measure. A zero divergence makes the corresponding scalar transport equation and the continuity equation coincide.

Charge conservation says that the total signed electric charge of an isolated collection stays constant. More generally, a collection closed to electric charge exchange has constant total electric charge, even if energy crosses its boundary. This is the charge bookkeeping relevant to a closed system or an isolated system; isolation from energy exchange is not required. Positive and negative electric charge count with opposite signs. Creating or removing an electron and a Positron together is compatible with charge conservation, because their combined electric charge is zero. ^[ref].

Imagine a fixed box initially containing electric charge $Q(0)=10\mu \mathrm{C}$. Stipulate incoming conventional current $I_{\mathrm{in}}=3\mu \mathrm{A}$ and outgoing conventional current $I_{\mathrm{out}}=1\mu \mathrm{A}$ throughout $0\le t\le 2\mathrm{s}$, with no other boundary transport. Since one ampere is one coulomb per second, the stored electric charge grows at $2\mu C/s$:If instead the outgoing electric current also equals $3\mu \mathrm{A}$, the stored electric charge remains $10\mu \mathrm{C}$ while electric charge passes through. These are stipulated transport examples; constant terminal electric currents alone do not guarantee a time-independent interior charge density.

Locally, charge conservation is the charge continuity equation, an instance of the general continuity equation:Here $\rho (\mathbf{r},t)$ is signed volume charge density and $\mathbf{J}(\mathbf{r},t)$ is current density. A positive divergence of current density means net outward transport from a small neighborhood, so its charge density decreases. For a fixed control volume $V$ with closed surface $\partial V$ and outward unit normal $\mathbf{n}$, the equivalent integrated balance isThe surface integral is signed: incoming positive electric charge contributes negatively to outward electric current, and outgoing positive electric charge contributes positively. Integrate the charge continuity equation over fixed $V$ and apply the divergence theorem to obtain this balance. Zero net boundary transport gives constant $Q_V$; the global and local statements describe the same conservation law at different scales. ^[ref].

In a steady transport state, charge density and current density are time-independent, so the charge continuity equation reduces to $\nabla \cdot \mathbf{J}=0$. At an electrical circuit junction, choose each terminal electric current $I_k$ positive outward and assume the terminals account for all boundary transport. ThenKirchhoff's current law, ${\sum }_k{I}_k=0$, follows when accumulation of electric charge at the junction is negligible. During storage, retain the same balance with its nonzero derivative; charge conservation continues to hold. This makes charge conservation useful for checking electrical circuit models and for tracking the charging of a capacitor. ^[ref].

The beautiful consistency with Maxwell equations is that the electric field and magnetic field already enforce this transport balance. For sufficiently smooth fields in vacuum SI, take the divergence of the Ampère-Maxwell equation, use that the divergence of a curl is zero, and substitute Gauss's law, $\nabla \cdot \mathbf{E}=\rho /{ϵ}_0$. Constant vacuum permittivity ${ϵ}_0$ givesThe displacement current term supplies the part needed when electric charge accumulates. Historically, Maxwell added that term to extend the steady-current equation; see the existing charge conservation from Maxwell equations for this connection. ^[ref].

The local form of conservation of electric charge relates charge density to current density. Taking the divergence of the Ampère-Maxwell equation and substituting Gauss's law proves the equation, using that the divergence of a curl is zero.

For a closed polyhedron, let $A_F$ be a face area and ${\mathbf{n}}_F$ its outward unit normal vector. The sum of the face area vectors is zero. To prove it, fix any constant vector $\mathbf{q}$ and apply the divergence theorem to the constant vector field $\mathbf{q}$: its divergence vanishes, while its outward flux is $\mathbf{q}\cdot {\sum }_F{A}_F{\mathbf{n}}_F$. This is zero for every $\mathbf{q}$, so the sum vanishes. For a tetrahedron, the same identity follows directly by expanding the four oriented face cross products; multiplying every face vector by two leaves the identity intact.

In Cartesian coordinates on four-dimensional Minkowski spacetime, all Conformal Killing vector fields have the formTheir conformal Killing factor is $c+2b\cdot x$. There are four translation generators, six Lorentz transformation generators, one dilation generator, and four special conformal transformation generators. These form a fifteen-dimensional Lie algebra under the Lie bracket of vector fields.

For completeness, differentiate the Conformal Killing equation in three index orders to obtainIts divergence gives $2{\partial }_{\mu }{\partial }_{\nu }\lambda =-{\eta }_{\mu \nu }\square \lambda$; its tensor contraction gives $\square \lambda =0$. Hence $\lambda$ is an affine function. Subtracting the displayed quadratic and dilation terms leaves a Killing vector field whose second partial derivatives vanish, proving completeness of the list. The vector fields are globally smooth, but some special conformal transformations have incomplete flows on affine Minkowski spacetime; the dilation flow $x\mapsto e^{s}x$ is global.

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.

For a vector field directed radially from a fixed axis and depending only on distance $\rho$ from that axis, cylindrical coordinates give $\nabla \cdot [f(\rho ){\mathbf{e}}_{\rho }]=f\prime +f/\rho$. Zero divergence gives $f=A/\rho$ on any positive radial interval. Unless $A=0$, this field is singular on the axis, so applications of the divergence theorem must use a domain avoiding it or account separately for that singularity.

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.

Diffusion is net transport associated with random microscopic motion. A localized patch of a dilute solute can spread through a stationary solvent: individual molecules move in both directions, but more cross from the crowded side of an imaginary surface than from the less crowded side. This imbalance gives a net diffusive flux down the concentration gradient; the molecules do not need to sense which direction has lower concentration. At uniform concentration, microscopic exchanges continue even though the net diffusive flux vanishes. ^[ref]

The distinction from advection is easiest to see in a dye patch. Uniform bulk motion carries its centre downstream; diffusion broadens it relative to that centre. Both can happen together. A moving or spreading patch by itself does not identify the mechanism: a nonuniform velocity field can also stretch a patch. ^[ref]

For a dilute nonreacting solute in a uniform medium at fixed temperature, Fick's first law models the relative diffusive flux as $\mathbf{j}=-D\nabla C$, with scalar diffusion coefficient $D>0$ of units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. Here $C$ is the amount of substance per unit volume, so $\mathbf{j}$ has units $\mathrm{mol}{\mathrm{m}}^{-2}{\mathrm{s}}^{-1}$. With carrier velocity field $\mathbf{u}$, the total flux is $C\mathbf{u}+\mathbf{j}$. The continuity equation, with no sources or chemical reactions, givesFor constant $D$ and incompressible flow, this reduces to $C_t+\mathbf{u}\cdot \nabla C=D{\nabla }^{2}C$, where ${\nabla }^2$ is the Laplacian. In a stationary carrier it is the diffusion equation. These assumptions matter: variable $D$ must remain inside the divergence, and changing sources require extra terms. ^[ref]

The useful surprise is the common mathematics: thermal conduction smooths temperature through the analogous heat equation, while molecular diffusion smooths composition. This analogy already motivated the original formulation of Fick's laws. ^[ref] However, the thermal diffusivity and the solutal diffusivity need not be equal. A diffusive time scale $t_D\sim L^2/D$ means that crossing ten times the distance requires roughly one hundred times the time. This is why molecular diffusion is especially effective across thin layers, and why chaotic advection can assist mixing by making those layers thinner. Physical diffusion is a transport mechanism, distinct from the broader probabilistic class of Markov diffusion.

In vacuum, displacement current is the term ${\mathbf{J}}_D={ϵ}_0{\partial }_t\mathbf{E}$ in the Ampère-Maxwell equation. It has the units of current density, but need not describe electric charge crossing the region. In Maxwell equations in matter, the corresponding term is ${\partial }_t\mathbf{D}$, with $\mathbf{D}$ the electric displacement field.

Consider a slowly charging parallel-plate capacitor with vacuum gap, plate area $A$, and charge $Q(t)$ on its positive plate. Neglect fringing and assume the charging timescale is long compared with the light-travel time across the apparatus. The nearly uniform gap electric field has magnitude $E=Q/({ϵ}_{0}A)$, so the displacement current through the gap isFor $A=0.010{\mathrm{m}}^2$ and charging electric current $I=1.0\mathrm{mA}$, this gives $dE/dt\simeq 1.13\times 10^{10}\mathrm{V}{\mathrm{m}}^{-1}{\mathrm{s}}^{-1}$ and $J_D=0.10A/{\mathrm{m}}^2$. No conducting electric current crosses the gap. A spanning surface cutting the wire and another passing between the plates nevertheless give the same magnetic circulation when the full Ampère-Maxwell equation is used. ^[ref].

The local consistency check is short. Without displacement current, taking the divergence of $\nabla \times \mathbf{B}={\mu }_0\mathbf{J}$ would force $\nabla \cdot \mathbf{J}=0$, allowing no local accumulation of electric charge. With displacement current, the divergence of a curl is zero and Gauss's law giveThis is precisely the charge continuity equation. Electric charge conservation fixes the required divergence of the added term; it does not by itself exclude other divergence-free additions. The specific Maxwell equations also yield the electromagnetic wave equation, linking this consistency correction to propagation at the speed of light.

For a right-handed regular system of orthogonal coordinates, the gradient in orthogonal coordinates gives $\nabla u_i={\mathbf{e}}_i/h_i$. Thus $\varphi {\mathbf{e}}_1=(\varphi h_1)\nabla u_1$, and the curl product rule givesAlso $\nabla u_2\times \nabla u_3={\mathbf{e}}_1/(h_2{h}_3)$ and the divergence of this cross product is zero, by the divergence-of-a-cross-product identity and the vanishing curl of each gradient. Applying the divergence product rule to $\varphi {\mathbf{e}}_1=(\varphi h_2{h}_3)(\nabla u_2\times \nabla u_3)$ givesCyclic permutation and addition give the full component formulas for any vector field. This derivation accounts for the varying coordinate basis without differentiating its unit vectors explicitly. The factors $h_i>0$ are scale factors of orthogonal coordinates; the chart must avoid coordinate singularities.

For $r=|x|>0$, the gradient satisfies ${\partial }_{i}r=x_i/r$. The product rule and Einstein notation then give the displayed divergence formula in dimension $n$. In dimension three, its derivative matrix is symmetric, so the curl vanishes. Zero divergence forces $f(r)=C{r}^{-n}$ on each positive radial interval. This formula does not assert regularity at the origin.

A differential operator is in divergence form when its highest derivatives occur inside a divergence, for example $Lu=\nabla \cdot (A\nabla u)$ or $(p(x)u^{\prime }{)}^{\prime }$. This form makes integration by parts express the operator through a lower-order bilinear form and boundary flux.

For a twice continuously differentiable vector field in Cartesian coordinates, expand the curl and then its divergence. Every mixed partial derivative occurs twice with opposite signs, proving the identity. In electromagnetism this forces the right side of the Ampère-Maxwell equation to be divergence-free, yielding local charge conservation.

Let $V$ be a bounded region in Euclidean space with piecewise smooth boundary and outward unit normal vector $n$. For a continuously differentiable vector field $F$ on a neighborhood of its closure,The total outward flux equals the accumulated divergence inside. Contributions across shared internal faces cancel when a region is divided into small pieces.

For $F=(x,y,z)$ on a ball of radius $R$, the divergence is $3$, giving $3(4\pi R^3/3)=4\pi R^3$. On the boundary, $F\cdot n=R$, so the surface integral is also $R(4\pi R^2)=\boxed{4\pi R^3}$. This agreement turns the derivative into a geometric flux measurement. The required regularity matters: a singular field at the center cannot be included without checking its contribution separately.

Double Compton scattering is $e^{-}+\gamma \to e^{-}+\gamma +\gamma$, together with its inverse. The forward reaction has one incoming and two outgoing photons, so it changes their number and can source the divergence of the photon number current. Ordinary Compton scattering has one incoming and one outgoing photon and conserves photon number. Forward and inverse rates can balance without either reaction individually conserving it.

A finite electric dipole consists of equal and opposite electric charges separated by a nonzero displacement. It has charge neutrality but a nonzero electric dipole moment $\mathbf{p}=q\mathbf{d}$. An ideal point electric dipole is the mathematical limit $d=|\mathbf{d}|\to 0$ with $\mathbf{p}$ fixed, so $q$ must increase as $1/d$. Shrinking the separation while holding $q$ fixed instead makes the electric dipole moment vanish. The ideal electric dipole describes an exterior electric field; it does not reproduce the two near-charge singularities or the internal structure of a molecule. ^[ref].

Place the pair's midpoint at the origin and write $k_e=1/(4\pi {ϵ}_0)$, where ${ϵ}_0$ is vacuum permittivity. By superposition principle and Coulomb's law, its exact electric potential away from the two point charges isFor $r=|\mathbf{r}|\gg d$, retain the first nonzero term of a Taylor series. The electric potential and its negative gradient becomeHere $\hat{\mathbf{r}}=\mathbf{r}/r$ and the products in parentheses are dot products. These are exact exterior formulas for the ideal electric dipole and distant approximations for a finite pair. The electric field on the outer axis points along $\mathbf{p}$ with magnitude $2k_{e}p/r^3$; on the perpendicular bisector it points opposite $\mathbf{p}$ with magnitude $k_{e}p/r^3$. This directional structure is what a scalar total electric charge misses. ^[ref].

In a prescribed uniform external electric field ${\mathbf{E}}_{\mathrm{ext}}$, a rigid pair has opposite forces ${\mathbf{F}}_{\pm }=\pm q{\mathbf{E}}_{\mathrm{ext}}$. Its net force is zero, but summing their moments about the midpoint gives the torqueThe cross product makes the torque perpendicular to both the electric dipole moment and the external electric field. If $\theta$ is the angle from ${\mathbf{E}}_{\mathrm{ext}}$ to $\mathbf{p}$ in their plane, the signed turning torque is ${\tau }_{\theta }=-p{E}_{\mathrm{ext}}\sin \theta$. It tends to reduce a small positive $\theta$. ^[ref].

The external interaction potential energy, with zero chosen at $\theta =\pi /2$, isParallel alignment minimizes this potential energy; antiparallel alignment maximizes it even though the torque is zero in both positions. This omits the pair's constant internal binding potential energy and assumes a fixed permanent electric dipole moment. ^[ref].

A nonuniform external electric field can exert a net force on a neutral electric dipole. Before taking the small-pair approximation, the exact external force is $\mathbf{F}=q[{\mathbf{E}}_{\mathrm{ext}}(\mathbf{R}+\mathbf{d}/2)-{\mathbf{E}}_{\mathrm{ext}}(\mathbf{R}-\mathbf{d}/2)]$, where $\mathbf{R}$ is the midpoint. A first-order Taylor series gives, for a fixed electric dipole moment and an external electric field varying slowly across $d$,in electrostatics. The last equality uses the vanishing curl of the external electric field, with $\mathbf{p}$ held fixed while differentiating. Thus uniform external electric fields turn a misaligned electric dipole, whereas gradients can also move its midpoint. For example, the local source-free electric field ${\mathbf{E}}_{\mathrm{ext}}=Gx\hat{\mathbf{x}}+(E_0-Gy)\hat{\mathbf{y}}$ gives $\mathbf{F}=(G{p}_x,-G{p}_y,0)$, so reversing the electric dipole moment reverses the net force. This stipulated local model has zero divergence and zero curl.

A motivating experiment is Filsinger and colleagues' study of electrostatic deflection of polar molecules. Benzonitrile (${\mathrm{C}}_7{\mathrm{H}}_5\mathrm{N}$) in helium expanded from a pulsed valve into vacuum; two 1-mm skimmers narrowed the beam before a 15-cm deflector. A grounded trough and a high-voltage rod produced a vertical electric field gradient. Scanning a focused laser vertically downstream and recording benzonitrile ions from photoionization measured the beam's relative density. With voltage off, the profile spanned about 2 mm; applying voltage shifted and broadened it upwards. At 10 kV, many molecules lay outside the original profile. ^[ref].

Interpretation uses a quantum state's Stark effect energy $W(E)$: ${\mu }_{\mathrm{eff}}=-dW/dE$ gives $F_y={\mu }_{\mathrm{eff}}{\partial }_{y}E$. States whose energy falls with increasing electric field seek stronger fields. The voltage-off profile controls the mechanical aperture. Inferring state populations needs calculated trajectories, speed distributions and a rotational temperature model; profile displacement alone does not measure a permanent electric dipole moment uniquely. This is an electric Stark effect experiment, with photoionization occurring after deflection. ^[ref].

For the symmetric finite pair, the first omitted electric multipole expansion term after the electric dipole is third order, so typical relative errors scale as $(d/r{)}^2$. This stronger error estimate need not hold for an arbitrary neutral molecule with a nonzero electric quadrupole moment. A concrete check for a pair on the $z$ axis compares the exact and ideal electric fields at $r>d/2$ on the positive axis and at $r>0$ on the perpendicular bisector:Both perpendicular-bisector components are negative. At $r=5d$, the ratios are approximately $1.0203$ and $0.9852$: errors of $+2.03\%$ and $-1.48\%$ relative to the ideal electric field. The condition $r\gg d$ concerns observing the pair's own electric field; the separate condition for its external force approximation concerns variation of the applied electric field across the pair.

Electric flux measures how much electric field crosses a chosen oriented surface, with a sign distinguishing the two crossing directions. A strong electric field running along a surface contributes nothing to its electric flux: only the normal component counts. This lets Gauss's law turn a potentially complicated electric field into one number recording the net enclosed electric charge.

For a surface $S$ with chosen unit normal $\mathbf{n}$, the electric flux is the surface integralHere $dA$ is positive surface area, while $d\mathbf{S}$ is the vector area element. The dot product makes electric flux a signed scalar, rather than a vector like the electric field. Its SI unit is $\mathrm{N}{\mathrm{m}}^2/C$: electric field measured using the newton per coulomb, multiplied by surface area in square metres. Reversing the orientation of a surface replaces $\mathbf{n}$ by $-\mathbf{n}$ and reverses the electric flux. Calling it electric flux does not imply that electric charge is moving through the surface. ^[ref].

For an open patch in a plane, with surface area $A$ and a constant electric field, the unit normal is constant and the surface integral simplifies toThe angle $\theta$ is between the electric field and the chosen unit normal, not between the electric field and the plane. The cosine describes the signed projected surface area seen by the electric field. Thus electric flux is positive for $0\le \theta <90^{\circ }$, zero at $90^{\circ }$, and negative for $90^{\circ }<\theta \le 180^{\circ }$. ^[ref].

For example, take $E=100N/C$ and $A=0.02{\mathrm{m}}^2$, so $EA=2\mathrm{N}{\mathrm{m}}^2/C$. Three choices of angle giveKeeping the $60^{\circ }$ patch fixed and reversing its unit normal gives $120^{\circ }$ and changes $+1$ to $-1\mathrm{N}{\mathrm{m}}^2/C$. The electric field and physical patch have not changed. If the electric field is instead tangent to the patch, its dot product with the unit normal vanishes: electric flux is zero even though $E=100N/C$.

For a closed surface enclosing a volume $V$, choose the outward unit normal. Gauss's law givesHere ${\epsilon }_0$ is the vacuum permittivity, and $Q_{\mathrm{enclosed}}$ is the signed total electric charge inside. The electric field in the surface integral includes all sources, including external electric charges. Assume no ideal point charge lies on the integration surface. Zero enclosed electric charge implies zero total electric flux; it does not imply zero electric field everywhere on the closed surface or inside it. Positive and negative contributions to the surface integral can cancel. ^[ref].

To see the cancellation with sources inside, place point charges $+q$ at $(-a,0,0)$ and $-q$ at $(a,0,0)$, where $q>0$, $a=0.30\mathrm{m}$, and both lie inside the sphere of radius $R=1\mathrm{m}$ centered at the origin. This separated pair has the electric dipole geometry. Its enclosed electric charge is $q-q=0$, so its total outward electric flux is zero. At the point $P=(0,0,R)$ on the sphere, Coulomb's law and the superposition principle give the individual electric fieldsThe $z$ components cancel and the positive $x$ components add:For the stipulated approximation $k=8.99\times 10^9\mathrm{N}{\mathrm{m}}^2/{\mathrm{C}}^2$ and $q=1\mathrm{nC}$, this is about $4.740N/C$ along positive $x$. At $P$, the outward unit normal is $\hat{\mathbf{z}}$, so this nonzero electric field is tangent to the sphere. At other points of the sphere, its normal component can be positive or negative. The source-field formula applies away from the two point charges, where the electric field has zero divergence; it is singular at the sources. ^[ref].

The diagram shows streamlines of the three-dimensional electric field in the section $y=0$. The circle is a cross-section of the sphere, not the integration surface: Gauss's law uses the entire three-dimensional surface integral. A two-dimensional line integral of outward field around that circle would be a different quantity, with units $\mathrm{N}m/C$ rather than the electric flux units $\mathrm{N}{\mathrm{m}}^2/C$.

In a linear isotropic nondispersive medium with time-independent permittivity $ϵ$ and magnetic permeability $\mu$, the stored electromagnetic energy per volume is $u=ϵ|\mathbf{E}{|}^2/2+|\mathbf{B}{|}^2/(2\mu )$. For a short derivation, the Maxwell equations in matter and the divergence identity for a cross product give $\nabla \cdot (\mathbf{E}\times \mathbf{H})=-\mathbf{H}\cdot {\partial }_t\mathbf{B}-\mathbf{E}\cdot {\partial }_t\mathbf{D}-{\mathbf{J}}_f\cdot \mathbf{E}$. The constitutive relations turn the first two terms into $-{\partial }_{t}u$, proving the Poynting theorem. Choose $u=0$ when the electric field and magnetic field both vanish. Dispersive or time-dependent material response requires a different accounting of stored energy.