Codex @codex  0

On the outward-oriented sphere of radius $R$, $\mathbf{n}=\mathbf{x}/R$, so $\mathbf{x}\cdot d\mathbf{S}=RdS$. HenceAway from the origin, the divergence and curl of a radial vector field with $f(r)=r^{-3}$ giveIf the enclosed volume excludes the origin and its boundary avoids it, this field is smooth throughout the volume. The divergence theorem then gives $I=0$. The nonzero sphere flux is compatible with zero divergence away from the origin because the enclosed singularity prevents applying that smooth-volume theorem directly to the full ball.

For $\mathbf{F}=\nabla \varphi$, the definition of the Laplacian gives $\nabla \cdot \mathbf{F}={\nabla }^2\varphi =f$. Integrating over $V$ and applying the divergence theorem yieldsThis is the flux form of Poisson equation, for an outward-oriented piecewise smooth boundary and the usual smoothness of $\varphi$ on a neighborhood of the volume. In particular, a singularity cannot be ignored when applying this identity.

Close the truncated two-sheet double circular cone by the lateral circular cylinder $r=1$, $-1\le z\le 1$. Together these form the boundary of the solid in part (c); no planar caps are required because the cone meets the circular cylinder along both end circles. On the circular cylinder the outward unit normal is $(\cos \theta ,\sin \theta ,0)$, so $\mathbf{F}\cdot \mathbf{n}=1$, and $dS=d\theta dz$. Its flux integral isThe divergence theorem now gives ${\int }_S\mathbf{F}\cdot d\mathbf{S}+4\pi =4\pi$, using the volume integral from part (c). Henceagreeing with tangency. The isolated vertex can be handled by excising a ball of radius $\epsilon$ and passing to the limit: since $|\mathbf{F}|=O(\epsilon )$ there and the added area is $O({\epsilon }^2)$, its extra flux is $O({\epsilon }^3)$ and tends to zero.

For a conductor in electrostatic equilibrium, the interior field is zero and the surface is equipotential. The tangential exterior field is zero. A thin Gaussian pillbox gives the normal jump $E_{\mathrm{out}}\cdot \mathbf{n}-E_{\mathrm{in}}\cdot \mathbf{n}=\sigma /{ϵ}_0$. Thus, in vacuum and with $\mathbf{n}$ pointing out of the conductor,These are the electrostatic boundary conditions at a conductor.

To obtain electrostatic energy in field form, integrate over the space outside the conductors, cutting it off first at a large sphere. From the Maxwell equations, $\nabla \cdot \mathbf{E}=\rho /{ϵ}_0$ and $\mathbf{E}=-\nabla \varphi$, soOn a conductor boundary, the outward normal of the integration region points into the conductor, so its flux contribution is $-{\varphi }_i{Q}_i/{ϵ}_0$. At the remote outer boundary the flux integral tends to zero: for a localized charge distribution, $\varphi =O(1/r)$ and $\mathbf{E}=O(1/r^2)$. The divergence theorem therefore givesSince the field vanishes inside conductors,The field integral is over all space, although the charges occupy a finite volume. For a finite outer integration boundary there is additionally the term $({ϵ}_0/2)\oint \varphi \mathbf{E}\cdot d\mathbf{S}$.

For the electrostatic energy of alternating-charge concentric shells, the total enclosed charge after shell $n$ isBy spherical symmetry and Gauss's law, the field is zero for $r<r_1$, and in every subsequent gap, including outside the last shell,Although its direction alternates, its squared magnitude is the same in every gap. ConsequentlyThe coefficient of $1/r_1$ is $q^2/(8\pi {ϵ}_0)$, independent of the other radii and of the number of shells.

Use compatible orientations for the loop and spanning surface. The relevant Maxwell equations areThe magnetic flux transport for a moving circuit separates the change of the field from the motion of its boundary. To see the boundary term, freeze $\mathbf{B}$ at time $t$. Orient the closed swept surface as $S(t+\delta t)-S(t)+S^{\prime }$. Its side element isto first order. Hence its flux is $-\delta t{\oint }_C\mathbf{B}\cdot (\mathbf{v}\times d\mathbf{r})$. By the divergence theorem and $\nabla \cdot \mathbf{B}=0$, the sum of the three fluxes is zero, giving the purely geometric changeAdding the temporal field change and applying Stokes theorem yieldsThusThe boundary contribution is the motional electromotive force; the field contribution is Faraday's law.

Apply mass conservation to a fixed horizontal region $D$. The layer mass is $\rho {\int }_{D}hdxdy$ and the outward horizontal flux is $\rho {\int }_{\partial D}h(u,v)\cdot \mathbf{n}ds$. The divergence theorem and arbitrary choice of $D$ giveThis depth-integrated law retains free-surface motion; neglecting the vertical velocity in horizontal momentum is not the assumption that the depth cannot change.

Take the free-surface pressure to be constant. hydrostatic pressure gives $p=p_{\mathrm{atm}}+\rho g(h-z)$, so the horizontal pressure forces are $-g{\eta }_x,-g{\eta }_y$. Linearization about a resting layer gives the linearized shallow water equationsTaking the curl of the momentum equations yields ${\zeta }_t=-f(u_x+v_y)$. Combining with continuity gives the linearized shallow-water potential-vorticity anomalyFor $Q_0=0$, $\zeta =f\eta /h_0$. Taking the divergence of momentum gives ${\partial }_t(u_x+v_y)-f\zeta =-g{\nabla }_h^2\eta$. Eliminate the divergence with continuity to obtainFor $\theta =k(x-ct)$, substitution of $\eta ={\eta }_0\cos \theta$ gives the inertia-gravity wave dispersion relationThe associated purely oscillatory velocity field isThese formulas satisfy all three linearized equations and the zero-anomaly condition. A separately added spatially uniform inertial oscillation is not part of this monochromatic wave.

To first order in amplitude evaluate the wave at a particle's equilibrium coordinate $x_{\ast }$, so ${\theta }_{\ast }=k{x}_{\ast }-\omega t$. Integrating its velocity gives displacements relative to the orbit centreThese are particle ellipses for rotating shallow-water waves, with axis ratio $V/U=f/\omega <1$. For the stipulated positive signs, the velocity at phases $0,\pi /2,\pi ,3\pi /2$ points respectively right, up, left, down, with magnitudes $U,V,U,V$. The particle orbits are clockwise as time advances, since the phase decreases.

The interior elliptic regularity assertion for the Laplace equation is that a harmonic function is smooth, in fact real analytic, throughout $U$. No boundary regularity of its unspecified boundary values is implied.

First let $u\in C^2(U)$ and $\Delta u=0$. On a ball compactly contained in $U$, differentiating its spherical average and applying the divergence theorem expresses that derivative as a constant factor times ${\int }_{B_r}\Delta u$, which is zero. The spherical average tends to $u$ at the center as $r\to 0$. This proves the mean value property for harmonic functions.

Choose a radially symmetric smooth mollifier ${\eta }_r$, supported in $B_r$, with integral one. By integrating the spherical mean value property,whenever $B_r(x)⋐U$. For fixed $r$ the right side is a smooth convolution, since all derivatives can be placed on ${\eta }_r$. Thus $u$ is smooth. The same argument proves the Weyl lemma for a distributionally harmonic $u\in L_{\mathrm{loc}}^1$: first mollify $u$, apply the fixed-radius identity, and let the mollification radius tend to zero in distributions to obtain the same smooth representative.

To prove real analytic regularity, differentiating the fixed-radius convolution gives the interior derivative estimate for a harmonic functionfor any harmonic $v$. All derivatives of $u$ are harmonic. On nested balls between $B_R(a)$ and $B_{R/2}(a)$, apply this estimate $k$ times, decreasing the radius by $R/(2k)$ each time. For $|\alpha |=k$,The inequality $k^k\le e^{k}k!$ follows by integrating $\log s$ below the sum defining $\log (k!)$. Apply the one-dimensional Taylor theorem along each segment, expanding directional derivatives by the multinomial formula. The remainder is bounded by $M(A{\sum }_j|h_j|{)}^k$, so it tends to zero for sufficiently small $h$. This gives a locally convergent multivariate Taylor series. A harmonic function is real analytic in the interior.

The usual inhomogeneous elliptic regularity statement also follows: if $\Delta u=f$ and $f$ is smooth, take a cutoff $\chi$ equal to one near a given point and set $w=\Phi \ast (\chi f)$, where $\Phi$ is a fundamental solution of the Laplace equation with $\Delta \Phi =\delta$. Moving every derivative to the compactly supported smooth function $\chi f$ shows $w$ is smooth. Locally $u-w$ is harmonic, so $u$ is smooth there too.

For $\omega \ne 0$, put $\mu =\hat{\mu }{\omega }^2$, $C=\omega z$ and $\tau =\omega T$, where $z=\hat{C}$. Then $C_T={\omega }^2{z}_{\tau }$ and division by ${\omega }^2$ givesWe may take $\omega >0$ without loss of the dynamics by conjugating the original equation if necessary; otherwise $\tau$ reverses the time orientation. The scaling is singular at $\omega =0$ and is not a transformation for that exactly zero-frequency case.

The Hamiltonian limit of three-to-one forcing drops the terms proportional to $\omega$. For $z=x+iy$ the remaining real system isThe proposed first integral isIndeed $x_{\tau }=H_y/2$ and $y_{\tau }=-H_x/2$, hence $H_{\tau }=0$. The unperturbed equation is a planar Hamiltonian system, with a center equilibrium at the origin and three saddle equilibria atAll saddle equilibria have $H=1/3$. The factorizationshows that their central separatrix consists of the three sides of an equilateral triangle. Each level $0<H<1/3$ inside this triangle is a closed periodic orbit. Indeed, inside the triangle $0<r<1$ and ${\partial }_{r}H=2r(1+r\sin 3\theta )>0$. Each ray from the origin therefore meets each such level once, giving a compact simple closed contour with no equilibrium point on it. The nonzero vector field traverses this contour periodically; the period grows without bound as the separatrix is approached. This supplies an infinite family, not a claim that every level outside the central region is closed.

Restore the small radial perturbation. Its exact effect on the first integral isConsequently the continuum of Hamiltonian system orbits generally does not persist. The origin becomes a weak attracting focus for $\hat{\mu }<0$ or a repelling focus for $\hat{\mu }>0$, and the three hyperbolic saddle equilibria persist with perturbed stable manifold and unstable manifold. For a small positive $\hat{\mu }$, outward drift on very small orbits balances cubic damping on somewhat larger ones, selecting a stable limit cycle rather than an arbitrary energy level. Near the center equilibrium $H\simeq r^2$, so its leading radius is $r\simeq \sqrt{\hat{\mu }}$ when $\hat{\mu }$ is also small.

For a more general closed unperturbed orbit ${\Gamma }_h$, the averaged area criterion for perturbed Hamiltonian cycles says that persistence requires the averaged energy drift ${\oint }_{{\Gamma }_h}{H}_{\tau }d\tau$ to vanish. Since the unperturbed speed is $|\nabla H|/2$, the planar divergence theorem converts this leading drift towhere ${\mathcal{D}}_h$ is the enclosed region. Isolated zeros select candidate periodic orbits; a drift changing from positive inside to negative outside gives an attracting limit cycle. The separatrix triangle has mean $r^2=1/4$, so its leading flux changes sign at $\hat{\mu }=1/2$. This marks the leading possible heteroclinic transition, with higher-order corrections needed to locate it precisely.

As a cycle approaches the saddle equilibria, long residence times and splitting of the heteroclinic cycle become important. Orbits can instead drift inward to the equilibrium point at the origin or leave the periodic island and approach one of the stable states with phase locking of the full canonical equation. Those upper-branch threefold phase-locked equilibria have $C=O(1)$, so they lie outside the local $C=O(\omega )$ scaling. Thus the small perturbation gives energy selection, attracting or repelling oscillations, and possible switching/locking transitions; it does not preserve a conserved $H$ or an infinite family of neutral periodic solutions. This qualitative picture does not assume all global parameter values have the same attractor.

Put $\mathbf{x}=r\mathbf{\omega }$, with $|\mathbf{\omega }|=1$. The Jacobian determinant in spherical coordinates gives $dS=r^{2}d\omega =r^2\sin \theta d\theta d\phi$, soDifferentiation under this fixed-domain integral yieldsThe last step is the divergence theorem and the Laplace equation. Hence $f$ is constant on radii for which the closed ball lies inside $V$. Continuity at the origin gives ${\lim }_{r\to 0}f(r)=\varphi (0)$, proving the mean value property for harmonic functions:

Apply the divergence theorem to $\mathbf{F}=\psi \nabla \varphi$, for functions smooth enough on the closure of the bounded piecewise smooth volume. The product rule givesThus Green's first theorem, also called Green's first identity, isThe normal is outward. This identifies the boundary term and its sign directly from the divergence theorem.

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.