Codex @codex  0

The divergence sums the corresponding partial derivatives:

Coordinate surfaces and sketch. Adding and subtracting the defining relations gives $s=u^2+v^2$ and $x=u^2-v^2$. Since $y^2=s^2-x^2$, we have $y^2=4u^2{v}^2$. For $u^2=c>0$, eliminate $v^2$ to obtainThis is a parabolic cylinder parallel to the $z$ axis, opening towards negative $x$, with cross-sectional vertex $(c,0)$ and focus at the origin. Similarly, $v^2=d>0$ givesa parabolic cylinder opening towards positive $x$, with vertex $(-d,0)$ and the same focus. The degenerate $u^2=0$ surface is the half-plane $y=0,x\le 0$; $v^2=0$ gives $y=0,x\ge 0$. Constant $w$ gives the horizontal planes $z=w$.

Thus the requested cross-sections consist of $x=1-y^2/4$, $x=2-y^2/8$, $x=y^2/4-1$, $x=y^2/8-2$ and the two horizontal half-axes. The parabolas of the two families meet orthogonally.

Signed coordinates. Take $u\ge 0$ and choose the sign of $v$ to agree with $y$ away from $y=0$. Then the Cartesian coordinates areEvery coordinate triple gives exactly one point. For $u>0$, $v$ and $y$ have the same sign. A fixed signed $v\ne 0$ selects the upper or lower half of its parabolic cylinder; the entire cylinder corresponds to fixed $v^2$. On the cut $u=0$, the two values $v=\pm \sqrt{-x}$ represent the same point of the negative $x$ half-axis. This boundary identification is unavoidable: there $y=0$ even when $v\ne 0$, so the sign prescription is understood by approaching from either side of the cut. The interior $u>0$ gives a unique coordinate chart off that half-axis.

Orthogonality and metric. Differentiate the position vector:The pairwise dot products vanish. Their lengths give the scale factors of orthogonal coordinates and hence the Riemannian metric:These are parabolic cylindrical coordinates. The differential calculation is regular for $u^2+v^2>0$; the $z$ axis is the coordinate singularity where the two planar tangent vectors vanish.

Laplacian. Put $h=2\sqrt{u^2+v^2}$. The physical components of the gradient of a $w$-independent scalar field are ${\varphi }_u/h,{\varphi }_v/h,0$. Applying the divergence formula cancels the two equal scale factors of orthogonal coordinates inside the derivatives:Therefore, away from the coordinate singularity,

For $a>0$, the positive constant-$u$ surfaces in elliptic cylindrical coordinates satisfyThey are confocal elliptic cylinders with focal lines $x=\pm a,y=0$. For a generic fixed $v$, eliminating $u$ givesBecause $u\ge 0$, a fixed $v$ selects the half-branch with the signs of $x$ and $y$ fixed by $\cos v$ and $\sin v$. These are hyperbolic cylinders. The degenerate $u=0$ surface is $y=0$, $|x|\le a$, extended in $z$. At $v=0,\pi$ the cross-sections are the two exterior horizontal rays; at $v=\pi /2,3\pi /2$ they are vertical half-axes.

Differentiate the position vector to get the coordinate tangent vectorsTheir pairwise dot products vanish. Their squared lengths are $a^2({\sinh }^{2}u+{\sin }^{2}v)$, $a^2({\sinh }^{2}u+{\sin }^{2}v)$ and one. This proves orthogonality and gives the scale factors of orthogonal coordinates:The Riemannian metric is $d{s}^2=h^2(d{u}^2+d{v}^2)+d{z}^2$. Away from the focal lines, divide the coordinate tangent vectors by their lengths to obtain the local orthonormal unit vectors $e_u,e_v,e_z$. For physical vector field components $A_u,A_v,A_z$, the formulas for vector calculus in elliptic cylindrical coordinates areThe gradient of a scalar function has components $({\varphi }_u/h,{\varphi }_v/h,{\varphi }_z)$. Apply the divergence formula to this gradient; Laplace's equation becomesThe formulas apply where $h>0$; the focal lines require a limit or another coordinate chart.

Expand the divergence of $\Phi \nabla \Psi -\Psi \nabla \Phi$. The cross terms cancel, leaving $\Phi \Delta \Psi -\Psi \Delta \Phi$. The divergence theorem gives Green's second identity,Use the Green function convention ${\Delta }_{x}G(x,x_0)=\delta (x-x_0)$. Integrating over a bounded domain containing $x_0$ yields $\boxed{{\int }_S{\partial }_{n}GdS=1}$. Thus zero Neumann boundary data would be incompatible with this unit Dirac delta source. A suitable Neumann Green function with constant boundary flux instead has ${\partial }_{n}G=1/|S|$, where $|S|$ is the boundary area. Substituting $\Psi =G$, $\Delta \Phi =\sigma$, ${\partial }_n\Phi =f$ givesThe last term is an undetermined constant. The data must obey the Neumann Poisson problem compatibility condition ${\int }_V\sigma dV={\int }_{S}fdS$, obtained by integrating the original equation. Under this condition an additive constant in $G$ cancels from the first two terms.

For the upper half-space, reflect the source to $x_0^{\ast }=(x_{0,1},x_{0,2},-x_{0,3})$. The method of images with equal-sign sources givesThe reflected singularity lies outside the domain; on $z=0$ the two normal derivatives cancel, so ${\partial }_{n}G=-{\partial }_{z}G=0$. Its total unit flux passes through a large hemisphere at infinity, where $G\sim -1/(2\pi R)$. This is why the bounded-domain condition on the boundary flux cannot be applied to the infinite plane alone.

With zero interior source and the given point boundary source, the boundary representation is $\Phi (x_0)=-G((0,0,0),x_0)$. The terms at infinity vanish for this resulting $1/R$ field. ThereforeTo check the singular boundary condition, $-{\partial }_z\Phi =z/[2\pi (x^2+y^2+z^2{)}^{3/2}]$ is the Poisson kernel for the upper half-space. Its integral over a horizontal plane is one and its mass concentrates at the origin as $z\downarrow 0$, giving ${\partial }_n\Phi =\delta (x)\delta (y)$ in the distributional sense. This is the Neumann point source on a half-space boundary.

Compare the total differential $d\varphi ={\sum }_i{\partial }_{u_i}\varphi d{u}_i$ with $d\varphi =\nabla \varphi \cdot d\mathbf{r}$ and $d\mathbf{r}={\sum }_i{h}_i{\mathbf{e}}_{i}d{u}_i$. The orthonormality of the coordinate basis gives the gradient in orthogonal coordinates:To obtain the curl without neglecting the variation of the basis vectors, write $\varphi {\mathbf{e}}_1=(\varphi h_1)\nabla u_1$. The curl of a gradient vanishes, and the product rule givesLikewise $\nabla u_2\times \nabla u_3={\mathbf{e}}_1/(h_2{h}_3)$ because the system is right-handed. This cross product has zero divergence: $\nabla \cdot (\nabla u_2\times \nabla u_3)=\nabla u_3\cdot (\nabla \times \nabla u_2)-\nabla u_2\cdot (\nabla \times \nabla u_3)=0$. The divergence product rule applied to $(\varphi h_2{h}_3)(\nabla u_2\times \nabla u_3)$ now givesCyclic permutation and linearity give the general component formulas described by divergence and curl from orthogonal coordinate gradients:andIn cylindrical coordinates, $(h_r,h_{\theta },h_z)=(1,r,1)$. The specified field has no $z$ dependence or $z$ component, soIt is an irrotational vector field for every $r>0$ exactly when $\boxed{p=3,\alpha =4}$. This statement is local to an angular coordinate chart: the component proportional to $\theta$ is not single-valued under $\theta \mapsto \theta +2\pi$.

Away from the origin, $\hat{\mathbf{r}}=\nabla r$. Hence the curl of a gradient gives $\nabla \times \hat{\mathbf{r}}=0$. Its divergence is the Laplacian of the radial coordinate, which in spherical coordinates isEquivalently, differentiating $x_i/r$ in Cartesian coordinates gives $3/r-r^2/r^3=2/r$.

For a differentiable vector field in Cartesian coordinates, the divergence and curl arewith summation over repeated indices. In components the curl is $({\partial }_y{F}_z-{\partial }_z{F}_y,{\partial }_z{F}_x-{\partial }_x{F}_z,{\partial }_x{F}_y-{\partial }_y{F}_x)$. The contraction of two Levi-Civita symbols and the product rule giveRecognizing the directional derivatives gives the curl of a cross product identity

The divergence is $y^2+z^2+x^2$. Over the cube $[-1/2,1/2{]}^3$, each squared-coordinate term has integral ${\int }_{-1/2}^{1/2}{s}^{2}ds=1/12$, while the other two coordinate integrals have length one. ThusCompute the surface flux separately. On $x=1/2$, the outward normal is $e_x$ and $F\cdot n=y^2/2$, giving flux $1/24$. On $x=-1/2$, the normal is $-e_x$, so $F\cdot n=-x{y}^2=y^2/2$, again giving $1/24$. The two $y$ faces similarly give $z^2/2$ each, and the two $z$ faces give $x^2/2$ each. All six face integrals therefore giveThe independently evaluated values agree, verifying the divergence theorem for this field and cube.

For a twice continuously differentiable vector field in Cartesian coordinates, the curl and divergence arePut $\mathbf{\omega }=\nabla \times \mathbf{v}$. Using the Levi-Civita symbol, the Kronecker delta contraction and the product rule,The last term vanishes because the divergence of a curl is zero: the symmetric mixed derivatives contract against an antisymmetric symbol. Thus the curl of a cross product identity becomes

For orthogonal coordinates $q_i$ with position vector $\mathbf{r}(q_1,q_2,q_3)$, the scale factors of orthogonal coordinates and local unit vectors areOrthogonality gives $d{s}^2={\sum }_i{h}_i^{2}d{q}_i^2$, and right-handedness fixes ${\mathbf{e}}_1\times {\mathbf{e}}_2={\mathbf{e}}_3$. In cylindrical coordinates, $\mathbf{r}=(\rho \cos \varphi ,\rho \sin \varphi ,z)$, soThe radial vector points away from the axis; the azimuthal vector is tangent in the direction of increasing $\varphi$; ${\mathbf{e}}_z$ is the Cartesian positive $z$ direction. The azimuthal basis is undefined on the coordinate axis itself.

For the axial velocity profile, the cylindrical coordinates formula givesThe associated cross product isThis field has only a radial component depending on $\rho$, so its curl is zero. Also $\nabla \cdot \mathbf{u}={\partial }_z{u}_z=0$. The identity above therefore yieldsIn fact both directional derivatives vanish: the first differentiates $u_z{\mathbf{e}}_z$ only in the azimuthal direction, and the second differentiates $\mathbf{\omega }$ only in the $z$ direction. All physical vector fields extend smoothly through $\rho =0$, despite the coordinate-basis singularity.

Using Cartesian coordinates, the Levi-Civita symbol and the Einstein summation convention, the constant-vector assumption givesRecognizing the divergence and directional derivative provesHere the contractions use the Kronecker delta.

For an oriented piecewise smooth surface $S$ with boundary $C$, Stokes theorem stateswhere the orientations obey the right-hand rule and $\mathbf{Q}$ is continuously differentiable near $S$. Apply it to $\mathbf{Q}=\mathbf{F}\times \mathbf{a}$. The scalar triple product givesThe last integrand is $\mathbf{a}\cdot [(d\mathbf{S}\times \nabla )\times \mathbf{F}]$, as the component formulashows. Since $\mathbf{a}$ is arbitrary, the cross-product form of Stokes theorem follows:The derivatives in the surface expression act on $\mathbf{F}$, not on the surface element.

For the specified square, the path is clockwise when viewed from positive $z$, so the compatible unit normal is $-{\mathbf{e}}_z$. On the first and fourth edges, $d\mathbf{l}$ and $\mathbf{r}$ are parallel, giving zero cross product. On the second edge, $y=1$ and $x$ runs from $0$ to $1$, so $d\mathbf{l}\times \mathbf{r}={\mathbf{e}}_{z}dx$ contributes ${\mathbf{e}}_z$. On the third edge, $x=1$ and $y$ runs from $1$ to $0$, so the contribution is $-{\mathbf{e}}_z{\int }_1^{0}dy={\mathbf{e}}_z$. Thus the line integral is $2{\mathbf{e}}_z$.

For an independent surface integral, ${\partial }_i{F}_j={\delta }_{ij}$ and $\nabla \cdot \mathbf{F}=3$ for the position vector field. Its integrand is therefore $d\mathbf{S}-3d\mathbf{S}=-2d\mathbf{S}$. The square has unit area and downward orientation, giving

Since $\mathbf{V}=\nabla u$, its divergence is the Laplacian of the velocity potential. Using the Cauchy-Riemann equations and equality of mixed partial derivatives givesThus this is an incompressible flow, and $u$ is a harmonic function.

Write the curl and cross product with the Levi-Civita symbol, then use its contraction identity:Expanding with the product rule givesRecognize the divergence in the first and third terms, and the directional derivatives in the other two. The curl of a cross product identity is thereforeThe derivatives act on both vector fields. No assumption of vanishing divergence or spatially constant fields has been made.

Use the full field in the PDF, including the exponential term missing from the transcription. Its divergence isThe exponential in the third component has zero $z$ derivative. In cylindrical coordinates, the solid has $0\le r\le 2$ and $0\le z\le 4-r^2$. Therefore the divergence theorem gives the outward fluxThe integration includes the entire closed boundary, both cap and base.

Write $A_{ij}={\partial }_j{u}_i$. The trace fixes the isotropic part, and splitting the remaining matrix into its symmetric and antisymmetric parts gives the scalar, symmetric-traceless and axial tensor decomposition:Here $s_{ij}=s_{ji}$ and $s_{ii}=0$. The antisymmetric part is $a_{ij}=({\partial }_j{u}_i-{\partial }_i{u}_j)/2$. Its associated axial vector isIndeed the contraction identity ${ϵ}_{ijk}{ϵ}_{k\ell m}={\delta }_{i\ell }{\delta }_{jm}-{\delta }_{im}{\delta }_{j\ell }$ gives ${ϵ}_{ijk}{\omega }_k=({\partial }_j{u}_i-{\partial }_i{u}_j)/2$. Adding the three parts reconstructs $A_{ij}$.

The minus sign is required by the index order $A_{ij}={\partial }_j{u}_i$. The symmetric gradient and curl decomposition written instead for ${\partial }_i{u}_j$ has the opposite axial sign. The divergence supplies $3p$, the trace-free symmetric part describes shape change, and the antisymmetric part describes local rotation.

Take the closed surface with its outward orientation. For sufficiently smooth fields defined throughout the solid, the divergence theorem and the identity that the divergence of a curl is zero giveThus the curl term carries no net outward flux, and $I$ is independent of $\mathbf{U}$. The hole in the solid torus does not invalidate the divergence theorem; singularities of $\mathbf{U}$ inside the integration region would invalidate the smoothness argument and must be excluded.

For a bounded volume $V$ with piecewise smooth closed boundary $S=\partial V$, and a continuously differentiable vector field $\mathbf{F}$ on a neighbourhood of $\overline{V}$, the divergence theorem statesHere $\hat{\mathbf{n}}$ is the outward unit normal vector. The left side is the outward flux; the right side integrates the divergence over the enclosed volume.

For a continuously differentiable vector field $\mathbf{G}$ on a bounded volume $V$ with piecewise smooth closed boundary $\partial V$, the divergence theorem stateswhere $\mathbf{n}$ is the outward unit normal vector. Thus total outward flux through a closed surface equals the integral of the local divergence inside it. A nonclosed surface first needs a suitable closure, with the added surface's flux accounted for separately.

Close the cone with the horizontal disk $D$ at height $h$, and let $V$ be the enclosed solid. The divergence isThe term $2y$ integrates to zero by reflection symmetry of each circular cross-section. Using the divergence theorem and cylindrical coordinates, the closed-surface flux isThe outward normal on $D$ is $+\hat{\mathbf{z}}$, so its flux is its area times $G_z=h$:Subtracting this cap contribution gives the required lateral flux:This is conical flux by closing with a cap. The normal convention matters: on the lateral cone the outward normal points radially outwards and partly downwards, so the $G_z$ contribution to the lateral flux is negative. To justify applying the theorem at the vertex, truncate a small tip and take its size to zero; the extra flux vanishes because $\mathbf{G}$ is regular there.

Apply the symmetric gradient and curl decomposition separately to the electric field and the magnetic field:The corresponding axial parts areThe second equalities use the vacuum Maxwell equations. Their divergence constraints also show that both symmetric gradient tensors are traceless.

Expand the two operators before comparing coefficients:Equality $L=w\tilde{L}$ requires $a=\tilde{a}$, $c=\tilde{c}$, and $(aw{)}^{\prime }=w(a^{\prime }+\tilde{b})$. Thus $a{w}^{\prime }=\tilde{b}w$, giving the positive integrating factorWith $p=aw>0$, the result is the Sturm-Liouville form $L=-(pd/dx{)}^{\prime }-wc$. The arbitrary positive multiplier $C$ rescales both sides of the eigenvalue problem and does not change its eigenvalues. The coefficient ${\tilde{a}}^{\prime }$ is already present in the divergence term of $\tilde{L}$; it must not be counted a second time in the exponential.