Calculus/Points, paths, surfaces, and volumes

This chapter will provide an intuitive interpretation of vector calculus using simple concepts such as multi-points, multi-paths, multi-surfaces, and multi-volumes. Scalar fields will not be simply treated as a function $f : ℝ^{3} \to ℝ$ that returns a number given an input point, and vector fields will not be simply treated as a function $𝐅 : ℝ^{3} \to ℝ^{3}$ that returns a vector given an input point.

Basic structures

The basic structures are multi-points, multi-paths, multi-surfaces, and multi-volumes.

Multi-points

A point $𝐪_{0}$ is an arbitrary position. A "multi-point" is a set of point/weight pairs: $𝐐 = {(𝐪_{1}, w_{1}), (𝐪_{2}, w_{2}), . . ., (𝐪_{k}, w_{k})}$ where $w_{i}$ is the "weight" that is assigned to point $𝐪_{i}$ . Given two point/weight pairs $(𝐪, w_{1})$ and $(𝐪, w_{2})$ that cover the same point $𝐪$ , the weights add up to get $(𝐪, w_{1} + w_{2})$ which replaces $(𝐪, w_{1})$ and $(𝐪, w_{2})$ . Any pair $(𝐪, 0)$ is removed. $𝐐$ can consist of infinitely many points, and each point may have an infinitesimal weight.

An arbitrary point $𝐪_{0}$ can be described by the scalar field $δ_{0} (𝐪; 𝐪_{0}) = {\begin{matrix} + \infty^{3} & (𝐪 = 𝐪_{0}) \\ 0 & (𝐪 \neq 𝐪_{0}) \end{matrix}$ . This is the "Dirac delta function" centered on point $𝐪_{0}$ . The $+ \infty^{3}$ is the inverse of an infinitely small volume that wraps point $𝐪_{0}$ . To further explain this, let $ω_{0} (𝐪_{0}, v)$ be a tiny volume with volume $v$ that wraps point $𝐪_{0}$ . $δ_{0} (𝐪; 𝐪_{0})$ can be approximated by $Δ_{0} (𝐪; 𝐪_{0}, v) = {\begin{matrix} 1 / v & (𝐪 \in ω_{0} (𝐪_{0}, v)) \\ 0 & (𝐪 \notin ω_{0} (𝐪_{0}, v)) \end{matrix}$ . A mass of 1 is being crammed into $ω_{0} (𝐪_{0}, v)$ yielding an infinitely high density. Since $δ_{0} (𝐪; 𝐪_{0})$ is essentially a density function, it brings with it the units $[{length}^{- 3}]$ .

Multi-point $𝐐 = {(𝐪_{1}, w_{1}), (𝐪_{2}, w_{2}), . . ., (𝐪_{k}, w_{k})}$ can be described by the scalar field $δ_{0} (𝐪; 𝐐) = \sum_{i = 1}^{k} w_{i} δ_{0} (𝐪; 𝐪_{i})$ . If $𝐐$ consists of infinitely many points with each point having infinitesimal weight, then $δ_{0} (𝐪; 𝐐)$ is a density function.

In the image below, the multi-point in the left panel is converted to the scalar field in the center panel by averaging the point weight over each cell. The volume of each cell should be infinitesimal. The multi-point in the right panel corresponds to the same scalar field, and is in a more canonical form where oppositely weighted points have cancelled out.

The image below shows how a continuous scalar field $ρ : ℝ^{3} \to ℝ$ can be generated as a collection of points. Consider position $𝐪_{0}$ and the infinitesimal volume $ω_{0} (𝐪_{0}, v)$ with volume $v$ . The total point weight contained by $ω_{0} (𝐪_{0}, v)$ is $\underset{𝐪 \in ω_{0} (𝐪_{0}, v)}{∭} ρ (𝐪) d V \approx v \cdot ρ (𝐪_{0})$ . This weight of $v \cdot ρ (𝐪_{0})$ is then split up over an arbitrarily large number of points that are scattered over the volume $ω_{0} (𝐪_{0}, v)$ .

In summary, a multi-point is denoted by a scalar field that quantifies the density at each point, and any scalar field that quantifies density at each point is best interpreted as a multi-point.

Multi-paths

A simple path (also called a simple curve) $C$ is an oriented continuous curve that extends from a starting point $C (0)$ to an ending point $C (1)$ . Intermediate points are indexed by $t \in [0, 1]$ and are denoted by $C (t)$ . A simple path should be continuous (no breaks), and may intersect or retrace itself. A "multi-path" is a set of simple-path/weight pairs: $𝐂 = {(C_{1}, w_{1}), (C_{2}, w_{2}), . . ., (C_{k}, w_{k})}$ where $w_{i}$ is the weight that is assigned to path $C_{i}$ . Given two path/weight pairs $(C, w_{1})$ and $(C, w_{2})$ that cover the same path $C$ , the weights add up to get $(C, w_{1} + w_{2})$ which replaces $(C, w_{1})$ and $(C, w_{2})$ . Any pair $(C, 0)$ is removed. In addition given two path/weight pairs $(C_{1}, w)$ and $(C_{2}, w)$ with the same weight $w$ and $C_{1} (1) = C_{2} (0)$ , then $C_{1}$ and $C_{2}$ can be linked end-to-end to get the pair $(C_{1} + C_{2}, w)$ which replaces $(C_{1}, w)$ and $(C_{2}, w)$ . Assigning a path a negative weight effectively reverses its orientation: if $- C$ denotes path $C$ with the opposite orientation, then $(C, - w)$ is equivalent to $(- C, w)$ . $𝐂$ can consist of infinitely many paths, and each path may have an infinitesimal weight.

An arbitrary curve $C$ can be described by the vector field $δ_{1} (𝐪; C) = {\begin{matrix} (+ \infty^{2}) \hat{𝐧} (𝐪; C) & (𝐪 \in C) \\ 𝟎 & (𝐪 \notin C) \end{matrix}$ . This is the "Dirac delta function" for the curve $C$ . $\hat{𝐧} (𝐪; C)$ is the unit length tangent vector to path $C$ at point $𝐪 \in C$ . $\hat{𝐧} (𝐪; C) = 𝟎$ if $𝐪 \notin C$ . If there are multiple tangent vectors due to $C$ intersecting itself, then $\hat{𝐧} (𝐪; C)$ is the sum of these tangent vectors. The $+ \infty^{2}$ is the inverse of the cross-sectional area of an infinitely thin tube that encloses $C$ . To further explain this, let $ω_{1} (C, a)$ be a thin tube with cross-sectional area $a$ that encloses $C$ . $δ_{1} (𝐪; C)$ can be approximated by $Δ_{1} (𝐪; C, a) = {\begin{matrix} (1 / a) {\hat{𝐧}}_{*} (𝐪; C, a) & (𝐪 \in ω_{1} (C, a)) \\ 𝟎 & (𝐪 \notin ω_{1} (C, a)) \end{matrix}$ . ${\hat{𝐧}}_{*} (𝐪; C, a)$ is the generalization of $\hat{𝐧} (𝐪; C)$ to the tube $ω_{1} (C, a)$ . A path weight of 1 is being crammed into the cross-sectional area of $ω_{1} (C, a)$ yielding an infinitely high path density. Since $δ_{1} (𝐪; C)$ is essentially a density over area, it brings with it the units $[{length}^{- 2}]$ .

The image to the right gives a depiction of the Dirac delta function for a simple curve. The vector field $δ_{1} (𝐪; C)$ is $𝟎$ everywhere outside of an infinitely thin tube that encloses the path. Inside the tube, the vectors are parallel to the path, and have a magnitude equal to the inverse of the cross-sectional area. The Dirac delta function is the limit as the tube becomes infinitely thin.

Multi-path $𝐂 = {(C_{1}, w_{1}), (C_{2}, w_{2}), . . ., (C_{k}, w_{k})}$ can be described by the vector field $δ_{1} (𝐪; 𝐂) = \sum_{i = 1}^{k} w_{i} δ_{1} (𝐪; C_{i})$ . If $𝐂$ consists of infinitely many paths with each path having infinitesimal weight, then $δ_{1} (𝐪; 𝐂)$ is a flow density function.

In the image below, the multi-path in the left panel is converted to the vector field in the center panel by computing the total displacement in each cell and averaging over the volume. The volume of each cell should be infinitesimal. The multi-path in the right panel corresponds to the same vector field, and is in a more canonical form where the individual paths do not cross each other.

In summary, a multi-path is denoted by a vector field that quantifies the path/flow density at each point, and any vector field that quantifies a flow density at each point (such as current density) is best interpreted as a multi-path. (Flow density is a vector that points in the net direction of a flow, and has a length equal to the flow rate per unit area through a surface that is perpendicular to the net flow.)

Multi-surfaces

A simple surface $σ$ is an oriented continuous surface. A simple surface may intersect or fold back on itself. A "multi-surface" is a set of simple-surface/weight pairs: $𝐒 = {(σ_{1}, w_{1}), (σ_{2}, w_{2}), . . ., (σ_{k}, w_{k})}$ where $w_{i}$ is the weight that is assigned to surface $σ_{i}$ . Given two surface/weight pairs $(σ, w_{1})$ and $(σ, w_{2})$ that cover the same surface $σ$ , the weights add up to get $(σ, w_{1} + w_{2})$ which replaces $(σ, w_{1})$ and $(σ, w_{2})$ . Any pair $(σ, 0)$ is removed. In addition given two surface/weight pairs $(σ_{1}, w)$ and $(σ_{2}, w)$ with the same weight $w$ , then $σ_{1}$ and $σ_{2}$ can be combined to get the pair $(σ_{1} + σ_{2}, w)$ which replaces $(σ_{1}, w)$ and $(σ_{2}, w)$ . Assigning a surface a negative weight effectively reverses its orientation: if $- σ$ denotes surface $σ$ with the opposite orientation, then $(σ, - w)$ is equivalent to $(- σ, w)$ . $𝐒$ can consist of infinitely many surfaces, and each surface may have an infinitesimal weight.

An arbitrary surface $σ$ can be described by the vector field $δ_{2} (𝐪; σ) = {\begin{matrix} (+ \infty) \hat{𝐧} (𝐪; σ) & (𝐪 \in σ) \\ 𝟎 & (𝐪 \notin σ) \end{matrix}$ . This is the "Dirac delta function" for the surface $σ$ . $\hat{𝐧} (𝐪; σ)$ is the unit length normal vector to surface $σ$ at point $𝐪 \in σ$ . $\hat{𝐧} (𝐪; σ) = 𝟎$ if $𝐪 \notin σ$ . If there are multiple normal vectors due to $σ$ intersecting itself, then $\hat{𝐧} (𝐪; σ)$ is the sum of these normal vectors. The $+ \infty$ is the inverse of the thickness of an infinitely thin membrane that encloses $σ$ . To further explain this, let $ω_{2} (σ, t)$ be a thin membrane with thickness $t$ that encloses $σ$ . $δ_{2} (𝐪; σ)$ can be approximated by $Δ_{2} (𝐪; σ, t) = {\begin{matrix} (1 / t) {\hat{𝐧}}_{*} (𝐪; σ, t) & (𝐪 \in ω_{2} (σ, t)) \\ 𝟎 & (𝐪 \notin ω_{2} (σ, t)) \end{matrix}$ . ${\hat{𝐧}}_{*} (𝐪; σ, t)$ is the generalization of $\hat{𝐧} (𝐪; σ)$ to the membrane $ω_{2} (σ, t)$ . A surface weight of 1 is being sandwiched into the thickness of $ω_{2} (σ, t)$ yielding an infinitely high surface density. Since $δ_{2} (𝐪; σ)$ is essentially a density over length, it brings with it the units $[{length}^{- 1}]$ .

Multi-surface $𝐒 = {(σ_{1}, w_{1}), (σ_{2}, w_{2}), . . ., (σ_{k}, w_{k})}$ can be described by the vector field $δ_{2} (𝐪; 𝐒) = \sum_{i = 1}^{k} w_{i} δ_{2} (𝐪; σ_{i})$ . If $𝐒$ consists of infinitely many surfaces with each surface having infinitesimal weight, then $δ_{2} (𝐪; 𝐒)$ is a rate-of-gain function.

In the image below, the multi-surface in the left panel is converted to the vector field in the center panel by computing the total surface in each cell and averaging over the volume. The volume of each cell should be infinitesimal. The multi-surface in the right panel corresponds to the same vector field, and is in a more canonical form where the individual surfaces do not intersect each other.

In summary, a multi-surface is denoted by a vector field that quantifies the rate of gain at each point. To describe the rate-of-gain, imagine that passing through a surface in the preferred direction gives "energy". The rate of gain is a vector that points in the direction that yields the greatest rate of energy increase per unit length, and has a length equal to the maximum rate of energy increase per unit length. Any vector field that quantifies a rate of gain at each point (such as a force field) is best interpreted as a multi-surface.

Multi-volumes

A volume $Ω$ is an arbitrary region of space. A "multi-volume" is a set of volume/weight pairs: $𝐔 = {(Ω_{1}, w_{1}), (Ω_{2}, w_{2}), . . ., (Ω_{k}, w_{k})}$ where $w_{i}$ is the "weight" that is assigned to volume $Ω_{i}$ . Given two volume/weight pairs $(Ω, w_{1})$ and $(Ω, w_{2})$ that cover the same volume $Ω$ , the weights add up to get $(Ω, w_{1} + w_{2})$ which replaces $(Ω, w_{1})$ and $(Ω, w_{2})$ . Any pair $(Ω, 0)$ is removed. In addition given two volume/weight pairs $(Ω_{1}, w)$ and $(Ω_{2}, w)$ with the same weight $w$ and $Ω_{1} \cap Ω_{2} = \emptyset$ , then $Ω_{1}$ and $Ω_{2}$ can be combined to get the pair $(Ω_{1} \cup Ω_{2}, w)$ which replaces $(Ω_{1}, w)$ and $(Ω_{2}, w)$ . $𝐔$ can consist of infinitely many volumes, and each volume may have an infinitesimal weight.

An arbitrary volume $Ω$ can be described by the scalar field $δ_{3} (𝐪; Ω) = {\begin{matrix} 1 & (𝐪 \in Ω) \\ 0 & (𝐪 \notin Ω) \end{matrix}$ . This is the "Dirac delta function" analog for volumes, and is essentially an indicator function that indicates whether or not a position is contained by $Ω$ or not, 1 being yes and 0 being no. Since $δ_{3} (𝐪; Ω)$ is simply an indicator function, it brings with it no units (it is dimensionless).

Multi-volume $𝐔 = {(Ω_{1}, w_{1}), (Ω_{2}, w_{2}), . . ., (Ω_{k}, w_{k})}$ can be described by the scalar field $δ_{3} (𝐪; 𝐔) = \sum_{i = 1}^{k} w_{i} δ_{3} (𝐪; Ω_{i})$ . If $𝐔$ consists of infinitely many volumes with each volume having infinitesimal weight, then $δ_{3} (𝐪; 𝐔)$ is a potential function.

In the image below, the multi-volume in the left panel is converted to the scalar field in the center panel by averaging the volume weight in each cell. The volume of each cell should be infinitesimal. The multi-volume in the right panel corresponds to the same scalar field, and is in a more canonical form where oppositely weighted volumes have cancelled out, and the remaining volume has diffused to fill each cell.

In summary, a multi-volume is denoted by a scalar field that quantifies a potential at each point, and any scalar field that quantifies a potential at each point is best interpreted as a multi-volume.

At infinity

An important requirement is that all multi-points, multi-paths, multi-surfaces, and multi-volumes not extend to infinity. All structures can extend to an arbitrarily large range, as long as this range is not unbounded. Allowing the structures to extend to infinity will cause problems in the later discussions.

Totals

These sections will describe the total weight of multi-points, the total displacement of multi-paths, the total surface of multi-surfaces, and the total volumes of multi-volumes.

Total point weight

Given a multi-point $𝐐 = {(𝐪_{1}, w_{1}), (𝐪_{2}, w_{2}), . . ., (𝐪_{k}, w_{k})}$ , the total point weight is clearly $\sum_{i = 1}^{k} w_{i}$ . Given a scalar field $ρ$ that denotes a multi-point, the total weight of $ρ$ is $\underset{𝐪 \in ℝ^{3}}{∭} ρ (𝐪) d V$ . Given a simple point $𝐪_{0}$ , the total weight is 1 so $\underset{𝐪 \in ℝ^{3}}{∭} δ_{0} (𝐪; 𝐪_{0}) d V = 1$ .

Total displacement

Given a simple path $C$ that starts at point $C (0)$ and ends at point $C (1)$ , the total displacement generated by $C$ is $\int_{𝐪 \in C} d 𝐪 = C (1) - C (0)$ . This displacement is solely dependent on the endpoints as indicated by the top image to the right.

The displacement generated by a closed loop is $𝟎$ .

Given a multi-path $𝐂 = {(C_{1}, w_{1}), (C_{2}, w_{2}), . . ., (C_{k}, w_{k})}$ , the total displacement generated by $𝐂$ is $\sum_{i = 1}^{k} w_{i} \int_{𝐪 \in C_{i}} d 𝐪 = \sum_{i = 1}^{k} w_{i} (C_{i} (1) - C_{i} (0))$ .

Given a vector field $𝐉$ that denotes a multi-path, the total displacement generated by $𝐉$ is $\underset{𝐪 \in ℝ^{3}}{∭} 𝐉 (𝐪) d V$ . Since the displacement generated by a simple path $C$ is $\int_{𝐪 \in C} d 𝐪 = C (1) - C (0)$ , it is the case that $\underset{𝐪 \in ℝ^{3}}{∭} δ_{1} (𝐪; C) d V = \int_{𝐪 \in C} d 𝐪 = C (1) - C (0)$ .

One important observation from $\int_{𝐪 \in C} d 𝐪 = \underset{𝐪 \in ℝ^{3}}{∭} δ_{1} (𝐪; C) d V$ is that given a path integral over path $C$ , the differential $d 𝐪$ is equal to $δ_{1} (𝐪; C) d V$ in a volume integral: $\int_{𝐪 \in C} f (𝐪, d 𝐪) = \underset{𝐪 \in ℝ^{3}}{∭} f (𝐪, δ_{1} (𝐪; C) d V)$ provided that function $f$ is linear in the second parameter. In the lower image to the right, the displacement differential $d 𝐪 = \hat{𝐧} \cdot Δ l$ is equated to the volume differential $(\frac{\hat{𝐧}}{Δ A}) d V = δ_{1} (𝐪; C) d V$ by diffusing the path over an infintely thin cross-sectional area $Δ A$ . The integrand at points outside of the infinitely thin tube is 0: for all points $𝐪 \notin C$ , $f (𝐪, δ_{1} (𝐪; C) d V) = f (𝐪, 𝟎) = 0$ .

Total surface vector

Given an arbitrary oriented surface $σ$ , its "counter-clockwise boundary", denoted by $\partial σ$ , is the boundary of $σ$ whose orientation is determined in the following manner: Looking at $σ$ so that the preferred direction through $σ$ is oriented towards the viewer, the boundary $\partial σ$ wraps $σ$ in a counter-clockwise direction.

Given a flat surface as shown in the image to the right, this surface can be quantified by the "surface vector" which is a vector that is perpendicular (normal) to the surface in the preferred orientation, and has a length equal to the area of the surface. In the image to the right, a flat surface has an area of $A$ and is oriented to be perpendicular to unit-length normal vector $\hat{𝐧}$ . The "surface vector" of this surface is $A \cdot \hat{𝐧}$ .

Given a non-flat surface $σ$ , the total surface vector of $σ$ is computed by summing the surface vectors of each infinitesimal portion of $σ$ . The total surface vector is $𝐒 = \iint_{𝐪 \in σ} d 𝐒$ .

In a manner similar to how the total displacement of a path is solely a function of the endpoints, the total surface vector of a surface is solely a function of its counter-clockwise boundary. This is not intuitive, and will be explained in greater detail below using two approaches:

Generalizing from surfaces in 2D space

Below are shown two images related to surface vectors in 2D space. The image to the left shows surface vectors in 2D space. In 2 dimensions, surfaces are called 1D surfaces and are similar to paths. The boundary of a 1D surface consists of 2 points. The surface vector of a 1D surface segment is a 90 degree rotation of the segment and is oriented in the direction of the surface's orientation. The total surface vector of a 1D surface is the sum of all of the surface vectors of the individual components. For each component of the surface, the surface vector is a 90 degree rotation of the displacement that traverses the component, so the total surface vector is a 90 degree rotation of the displacement between the points that form the boundary of the surface. This proves that in two dimensions, the total surface vector depends only on the boundary of the 1D surface.

The image to the right extrudes the 1D surfaces in two dimensional space into 2D "ribbons" in 3 dimensional space. At the top a closed "ribbon" is shown. This "ribbon" is a surface that is always parallel to the vertical dimension, and whose boundary forms two identical loops that are vertically displaced from each other. The boundary loops are also perpendicular to the vertical dimension. The ribbon itself is partitioned into tiny rectangles whose height is equal to that of the ribbon. To the bottom left, a view of the same ribbon from the top down is shown. It can be seen that the length of the each surface vector is proportional to the length of the corresponding rectangular segment, since the heights are all uniform. To the bottom right, by rotating the surface vectors 90 degrees around the vertical dimension, the surface vectors now sum to $𝟎$ , so the sum of the unrotated surface vectors is also $𝟎$ .

The fact that the total surface vector of a closed ribbon is $𝟎$ means that if relief is added to a surface without changing its boundaries, the total surface vector is conserved. The two left images below give examples of distorting the interior of a surface by hammering in relief. The vertical surfaces introduced by the relief are ribbons which contribute $𝟎$ to the total surface vector, while the horizontal surfaces are simply displaced vertically be the relief. The rightmost image below shows how the total surface vector is preserved if the "texture" of the surface at infinitesimal scales is converted from "steps" (a union of horizontal and vertical surfaces) to "smooth slopes" and vice versa. The surface formed from the red and green planes is a step, while the surface formed from the blue plane is a slope. These two surfaces can be seen to have equal total surface vectors from the right-angled triangle at the right side of the image.

Generalizing from displacement vectors

The total displacement along a simple oriented curve can be used to compute the net displacement in a specific direction. Given a simple oriented curve $C$ and an oriented straight line with the direction indicated by normal vector $\hat{𝐧}$ , the total displacement $Δ 𝐪$ along $C$ can be used to compute the net displacement in the direction indicated by the line. This displacement is $\hat{𝐧} \cdot Δ 𝐪$ , and depends only on the endpoints of the curve.

In a direct analogy, given a simple oriented surface $σ$ with counter-clockwise boundary $\partial σ$ , and an oriented flat plane whose surface normal is $\hat{𝐧}$ , a quantity of interest is the total signed area of $σ$ perpendicularly projected onto the plane. The signed area that is projected by a flat infinitesimal portion of $σ$ with surface vector $d 𝐒$ is $\hat{𝐧} \cdot d 𝐒$ , and the total signed area is $\iint_{𝐪 \in σ} \hat{𝐧} \cdot d 𝐒 = \hat{𝐧} \cdot \iint_{𝐪 \in σ} d 𝐒 = \hat{𝐧} \cdot 𝐒$ where $𝐒$ is the total surface vector of $σ$ .

The total signed projected area $\hat{𝐧} \cdot 𝐒$ onto the plane is purely a function of the boundary $\partial σ$ , and does not depend on how $σ$ fills its boundary $\partial σ$ . This is much more obvious and clearer than the claim that the total surface vector $𝐒$ is only a function of $\partial σ$ : the area enclosed by a boundary in 2D space is purely a function of that boundary. Since the projected area is signed, "upside down" surfaces project negative area, and folds and overhangs cancel each other out.

Since $\hat{𝐧} \cdot 𝐒$ is purely a function of $\partial σ$ for all choices of plane orientation $\hat{𝐧}$ , then the total surface vector $𝐒$ is purely a function of $\partial σ$ .

Summary

The total surface vector generated by a closed surface is $𝟎$ .

Given a multi-surface $𝐒 = {(σ_{1}, w_{1}), (σ_{2}, w_{2}), . . ., (σ_{k}, w_{k})}$ the total surface vector generated by $𝐒$ is $\sum_{i = 1}^{k} w_{i} \iint_{𝐪 \in σ_{i}} d 𝐒$ .

Given a vector field $𝐅$ that denotes a multi-surface, the total surface vector generated by $𝐅$ is $\underset{𝐪 \in ℝ^{3}}{∭} 𝐅 (𝐪) d V$ . Since the surface vector generated by simple surface $σ$ is $\iint_{𝐪 \in σ} d 𝐒$ , it is the case that $\underset{𝐪 \in ℝ^{3}}{∭} δ_{2} (𝐪; σ) d V = \iint_{𝐪 \in σ} d 𝐒$ . One important observation is that given a surface integral over $σ$ , the differential $d 𝐒$ is equal to $δ_{2} (𝐪; σ) d V$ in a volume integral: $\iint_{𝐪 \in σ} f (𝐪, d 𝐒) = \underset{𝐪 \in ℝ^{3}}{∭} f (𝐪, δ_{2} (𝐪; σ) d V)$ provided that function $f$ is a linear in the second parameter.

Total volume

Consider a multi-volume $𝐔 = {(Ω_{1}, w_{1}), (Ω_{2}, w_{2}), . . ., (Ω_{k}, w_{k})}$ , where the volumes of $Ω_{1}, Ω_{2}, . . ., Ω_{k}$ are respectively $V_{1}, V_{2}, . . ., V_{k}$ , then the total volume of $𝐔$ is $\sum_{i = 1}^{k} w_{i} V_{i}$ . Each volume $V_{i}$ can be computed by $V_{i} = \underset{𝐪 \in Ω_{i}}{∭} d V = \underset{𝐪 \in ℝ^{3}}{∭} δ_{3} (𝐪; Ω_{i}) d V$ . The total volume of $𝐔$ is $V = \sum_{i = 1}^{k} w_{i} V_{i} = \sum_{i = 1}^{k} w_{i} \underset{𝐪 \in Ω_{i}}{∭} d V = \sum_{i = 1}^{k} w_{i} \underset{𝐪 \in ℝ^{3}}{∭} δ_{3} (𝐪; Ω_{i}) d V$ $= \underset{𝐪 \in ℝ^{3}}{∭} (\sum_{i = 1}^{k} w_{i} δ_{3} (𝐪; Ω_{i})) d V = \underset{𝐪 \in ℝ^{3}}{∭} δ_{3} (𝐪; 𝐔) d V$ .

If a multi-volume $𝐔$ can be denoted by scalar field $U$ , then the volume of $𝐔$ is $\underset{𝐪 \in ℝ^{3}}{∭} U (𝐪) d V$ .

Given an arbitrary volume $Ω$ , a volume integral over $Ω$ can be converted to a volume integral over $ℝ^{3}$ by replacing the differential $d V$ with $δ_{3} (𝐪; Ω) d V$ :

$\underset{𝐪 \in Ω}{∭} f (𝐪, d V) = \underset{𝐪 \in ℝ^{3}}{∭} f (𝐪, δ_{3} (𝐪; Ω) d V)$ provided that $f$ is linear in the second parameter.

Intersections

The union of two multi-points denoted by scalar fields $ρ_{1}$ and $ρ_{2}$ is simply $ρ_{1} + ρ_{2}$ , and the same is true for the union of two multi-paths, the union of two multi-surfaces, and the union of two multi-volumes. The union of two structures with different types, such as a multi-point with a multi-path, is forbidden however.

Unions
structure	multi-point $ρ_{2}$	multi-path $𝐉_{2}$	multi-surface $𝐅_{2}$	multi-volume $U_{2}$
multi-point $ρ_{1}$	multi-point $ρ_{1} + ρ_{2}$	n/a	n/a	n/a
multi-path $𝐉_{1}$	n/a	multi-path $𝐉_{1} + 𝐉_{2}$	n/a	n/a
multi-surface $𝐅_{1}$	n/a	n/a	multi-surface $𝐅_{1} + 𝐅_{2}$	n/a
multi-volume $U_{1}$	n/a	n/a	n/a	multi-volume $U_{1} + U_{2}$

The intersection on the other hand, is less trivial and can occur between structures of different types.

Point-Volume intersections

When a point $𝐪$ with weight $w_{1}$ intersects a volume $Ω$ with weight $w_{2}$ , then the intersection is point $𝐪$ with weight $w_{1} w_{2}$ . Given a multi-point and a multi-volume, the intersection is the sum of the pair-wise intersections of each simple point with each simple volume. The image below gives an example of the intersection of a multi-point with a multi-volume.

Given a multi-point with scalar field $ρ$ , and a multi-volume with scalar field $U$ , then the intersection is a multi-point with scalar field $ρ U$ .

The total intersection between a multi-point $ρ$ and a multi-volume $U$ is $\underset{𝐪 \in ℝ^{3}}{∭} ρ (𝐪) U (𝐪) d V$ .

If $ρ$ denotes a simple point $𝐪_{0}$ , then the total intersection is $\underset{𝐪 \in ℝ^{3}}{∭} δ_{0} (𝐪; 𝐪_{0}) U (𝐪) d V = U (𝐪_{0})$ .

If $U$ denotes a simple volume $Ω$ , then the total intersection is $\underset{𝐪 \in ℝ^{3}}{∭} ρ (𝐪) δ_{3} (𝐪; Ω) d V = \underset{𝐪 \in Ω}{∭} ρ (𝐪) d V$ .

Path-Surface intersections

When a path $C$ with weight $w_{1}$ intersects a surface $σ$ with weight $w_{2}$ at point $𝐪$ , then the intersection is point $𝐪$ with weight $\pm w_{1} w_{2}$ . The weight is $+ w_{1} w_{2}$ if $C$ passes through $σ$ in the direction in which $σ$ is oriented. The weight is $- w_{1} w_{2}$ if $C$ passes through $σ$ opposite the direction in which $σ$ is oriented. Given a multi-path and a multi-surface, the intersection is the sum of the pair-wise intersections of each simple path with each simple surface. The images below give examples of the intersections of a multi-path with a multi-surface.

In the image above to the far right, the multi-path is denoted by a vector field which has the value $𝐅$ inside the blue tube, and is $𝟎$ everywhere else. The multi-surface is denoted by a vector field which has the value $𝐆$ among the red sheets, and is $𝟎$ everywhere else. The total path weight in the blue tube is $| 𝐅 | Δ A$ . The total surface weight in the red sheets is $| 𝐆 | Δ t$ . The total weight of all the intersection points is $(| 𝐅 | Δ A) (| 𝐆 | Δ t) = | 𝐅 | | 𝐆 | Δ A Δ t$ . The volume that the intersection points are evenly spread out in is $Δ A Δ t / \cos θ$ . The intersection point density is $\frac{| 𝐅 | | 𝐆 | Δ A Δ t}{Δ A Δ t / \cos θ} = | 𝐅 | | 𝐆 | \cos θ = 𝐅 \cdot 𝐆$ .

Given a multi-path with vector field $𝐉$ , and a multi-surface with vector field $𝐅$ , then the intersection is a multi-point with scalar field $𝐉 \cdot 𝐅$ .

The total intersection between a multi-path $𝐉$ and a multi-surface $𝐅$ is $\underset{𝐪 \in ℝ^{3}}{∭} (𝐉 (𝐪) \cdot 𝐅 (𝐪)) d V$ .

If $𝐉$ is a simple path $C$ , then the total intersection is $\underset{𝐪 \in ℝ^{3}}{∭} (δ_{1} (𝐪; C) \cdot 𝐅 (𝐪)) d V = \int_{𝐪 \in C} 𝐅 (𝐪) \cdot d 𝐪$ .

If $𝐅$ is a simple surface $σ$ , then the total intersection is $\underset{𝐪 \in ℝ^{3}}{∭} (𝐉 (𝐪) \cdot δ_{2} (𝐪; σ)) d V = \iint_{𝐪 \in σ} 𝐉 (𝐪) \cdot d 𝐒$ .

Path-Volume intersections

When a path $C$ with weight $w_{1}$ intersects a volume $Ω$ with weight $w_{2}$ , then the intersection is path $C \cap Ω$ with weight $w_{1} w_{2}$ . Given a multi-path and a multi-volume, the intersection is the sum of the pair-wise intersections of each simple path with each simple volume. The image below gives an example of the intersection of a multi-path with a multi-volume.

Given a multi-path with vector field $𝐉$ , and a multi-volume with scalar field $U$ , then the intersection is a multi-path with vector field $𝐉 U$ .

The total intersection between a multi-path $𝐉$ and a multi-volume $U$ is $\underset{𝐪 \in ℝ^{3}}{∭} 𝐉 (𝐪) U (𝐪) d V$ .

If $𝐉$ denotes a simple path $C$ , then the total intersection is $\underset{𝐪 \in ℝ^{3}}{∭} δ_{1} (𝐪; C) U (𝐪) d V = \int_{𝐪 \in C} U (𝐪) d 𝐪$ .

If $U$ denotes a simple volume $Ω$ , then the total intersection is $\underset{𝐪 \in ℝ^{3}}{∭} 𝐉 (𝐪) δ_{3} (𝐪; Ω) d V = \underset{𝐪 \in Ω}{∭} 𝐉 (𝐪) d V$ .

Surface-Surface intersections

When a surface $σ_{1}$ with weight $w_{1}$ intersects a surface $σ_{2}$ with weight $w_{2}$ , then the intersection is the path $σ_{1} \cap σ_{2}$ with weight $w_{1} w_{2}$ . The orientation given to path $σ_{1} \cap σ_{2}$ is defined as follows: viewing the intersection where the surface normal vectors of $σ_{1}$ and $σ_{2}$ are oriented towards the viewer, the intersection path has $σ_{1}$ to its right, and $σ_{2}$ to its left. Put another way, the intersection path is oriented according to the "right-hand rule" where the surface normals of $σ_{1}$ are the "x" direction, and the surface normals of $σ_{2}$ are the "y" direction. The images below give examples of the intersections of a multi-surface with a multi-surface.

In the image above to the right, the first multi-surface is denoted by a vector field that has the value $𝐅$ among the blue sheets, and is $𝟎$ everywhere else. The second multi-surface is denoted by a vector field that has the value $𝐆$ among the red sheets, and is $𝟎$ everywhere else. The total surface weight in the blue sheets is $| 𝐅 | Δ t_{1}$ , and the total surface weight in the red sheets is $| 𝐆 | Δ t_{2}$ . The total weight of all the intersection paths is $(| 𝐅 | Δ t_{1}) (| 𝐆 | Δ t_{2}) = | 𝐅 | | 𝐆 | Δ t_{1} Δ t_{2}$ . The cross-sectional area that the intersection paths are evenly spread out over is $Δ t_{1} Δ t_{2} / \sin θ$ . The intersection path density is $\frac{| 𝐅 | | 𝐆 | Δ t_{1} Δ t_{2}}{Δ t_{1} Δ t_{2} / \sin θ} = | 𝐅 | | 𝐆 | \sin θ = | 𝐅 \times 𝐆 |$ . Lastly, it should be noted that the intersection paths are oriented out of the screen as per the right-hand rule.

Given a multi-surface with vector field $𝐅_{1}$ , and a multi-surface with vector field $𝐅_{2}$ , then the intersection is the multi-path with vector field $𝐅_{1} \times 𝐅_{2}$ .

The total intersection between multi-surface $𝐅_{1}$ and multi-surface $𝐅_{2}$ is $\underset{𝐪 \in ℝ^{3}}{∭} (𝐅_{1} (𝐪) \times 𝐅_{2} (𝐪)) d V$ .

If $𝐅_{2}$ denotes a simple surface $σ$ , then the total intersection is $\underset{𝐪 \in ℝ^{3}}{∭} (𝐅_{1} (𝐪) \times δ_{2} (𝐪; σ)) d V = \iint_{𝐪 \in σ} 𝐅_{1} (𝐪) \times d 𝐒$ .

Surface-Volume intersections

When a surface $σ$ with weight $w_{1}$ intersects a volume $Ω$ with weight $w_{2}$ , then the intersection is surface $σ \cap Ω$ with weight $w_{1} w_{2}$ . Given a multi-surface and a multi-volume, the intersection is the sum of the pair-wise intersections of each simple surface with each simple volume. The image below gives an example of the intersection of a multi-surface with a multi-volume.

Given a multi-surface with vector field $𝐅$ , and a multi-volume with scalar field $U$ , then the intersection is a multi-surface with vector field $𝐅 U$ .

The total intersection between a multi-surface $𝐅$ and a multi-volume $U$ is $\underset{𝐪 \in ℝ^{3}}{∭} 𝐅 (𝐪) U (𝐪) d V$ .

If $𝐅$ denotes a simple surface $σ$ , then the total intersection is $\underset{𝐪 \in ℝ^{3}}{∭} δ_{2} (𝐪; σ) U (𝐪) d V = \iint_{𝐪 \in σ} U (𝐪) d 𝐒$ .

If $U$ denotes a simple volume $Ω$ , then the total intersection is $\underset{𝐪 \in ℝ^{3}}{∭} 𝐅 (𝐪) δ_{3} (𝐪; Ω) d V = \underset{𝐪 \in Ω}{∭} 𝐅 (𝐪) d V$ .

Volume-Volume intersections

When a volume $Ω_{1}$ with weight $w_{1}$ intersects a volume $Ω_{2}$ with weight $w_{2}$ , then the intersection is the volume $Ω_{1} \cap Ω_{2}$ with weight $w_{1} w_{2}$ . Given two multi-volumes, the intersection is the sum of the pair-wise intersections of each simple volume from the first multi-volume with each simple volume from the second multi-volume. The image below gives an example of the intersection between two multi-volumes.

Given a multi-volume with scalar field $U_{1}$ , and a multi-volume with scalar field $U_{2}$ , then the intersection is a multi-volume with scalar field $U_{1} U_{2}$ .

The total intersection between multi-volume $U_{1}$ and multi-volume $U_{2}$ is $\underset{𝐪 \in ℝ^{3}}{∭} U_{1} (𝐪) U_{2} (𝐪) d V$ .

If $U_{2}$ denotes a simple volume $Ω$ , then the total intersection is $\underset{𝐪 \in ℝ^{3}}{∭} U_{1} (𝐪) δ_{3} (𝐪; Ω) d V = \underset{𝐪 \in Ω}{∭} U_{1} (𝐪) d V$ .

Other intersections

Other types of intersections, such as Point-Point intersections, Point-Path intersections, Point-Surface intersections, and Path-Path intersections, are not considered since these kinds of intersections can occur only by design. For example, the probability that two randomly chosen points will intersect each other is 0, but if a point and a volume are randomly chosen, then the probability of the point landing in the volume is nonzero. Given two unrelated points, these two points will never land on each other, since a prior relationship has to exist for the points to coincide. Below is summarized the various types of intersections:

Intersections
structure	multi-point $ρ_{2}$	multi-path $𝐉_{2}$	multi-surface $𝐅_{2}$	multi-volume $U_{2}$
multi-point $ρ_{1}$	n/a	n/a	n/a	multi-point $ρ_{1} U_{2}$
multi-path $𝐉_{1}$	n/a	n/a	multi-point $𝐉_{1} \cdot 𝐅_{2}$	multi-path $𝐉_{1} U_{2}$
multi-surface $𝐅_{1}$	n/a	multi-point $𝐅_{1} \cdot 𝐉_{2}$	multi-path $𝐅_{1} \times 𝐅_{2}$	multi-surface $𝐅_{1} U_{2}$
multi-volume $U_{1}$	multi-point $U_{1} ρ_{2}$	multi-path $U_{1} 𝐉_{2}$	multi-surface $U_{1} 𝐅_{2}$	multi-volume $U_{1} U_{2}$

Boundaries

The endpoints of paths

Given a simple path $C$ that starts at point $C (0)$ and ends at point $C (1)$ , the "endpoints" of $C$ is the multi-point ${(C (0), + 1), (C (1), - 1)}$ that consists of the starting point with a weight of +1, and the final point with a weight of -1. While $C$ is denoted by the vector field $δ_{1} (𝐪; C)$ , the endpoints are denoted by the scalar field $δ_{0} (𝐪; C (0)) - δ_{0} (𝐪; C (1))$ . The image below gives several examples of simple paths and their associated endpoints.

Given a multi-path $𝐂 = {(C_{1}, w_{1}), (C_{2}, w_{2}), . . ., (C_{k}, w_{k})}$ , the endpoints of $𝐂$ is the multi-point ${(C_{1} (0), + 1), (C_{1} (1), - 1), (C_{2} (0), + 1), (C_{2} (1), - 1), . . ., (C_{k} (0), + 1), (C_{k} (1), - 1)}$ .

Given a vector field $𝐉$ that denotes a multi-path, the multi-point that denotes the endpoints of $𝐉$ is denoted by scalar field $\nabla \cdot 𝐉$ . The scalar field $\nabla \cdot 𝐉$ evaluated at point $𝐪$ is denoted by $\nabla \cdot 𝐉 (𝐪)$ , $(\nabla \cdot 𝐉) (𝐪)$ or $\nabla \cdot 𝐉 |_{𝐪}$ .

The requirement that no path extends to infinity means that every starting point is paired with a final point, and therefore the total weight of all of the endpoints together is 0: $\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \cdot 𝐉 (𝐪)) d V = 0$ .

The similarity of the notation $\nabla \cdot 𝐉$ to the intersection of multi-path $𝐉$ with multi-surface $𝐅$ , denoted by $𝐅 \cdot 𝐉$ , makes sense if $\nabla$ is interpreted as the "surface of reality". A starting point forms when a path pokes into reality, and a final point forms when a path pokes out of reality.

In the image to the right, a depiction of the "surface of reality" interpretation of $\nabla$ is shown. On the right is a simple path $𝐅$ , along with its endpoints $\nabla \cdot 𝐅$ . On the left $𝐅_{ext}$ is an extension of $𝐅$ that is behind the "veil" of surface $𝐆_{\nabla}$ . $𝐅_{ext}$ pokes out of and into $𝐆_{\nabla}$ at points consistent with the endpoints of $𝐅$ : i.e. $𝐆_{\nabla} \cdot 𝐅_{ext} = \nabla \cdot 𝐅$ .

The counter-clockwise oriented boundaries of surfaces

Given an oriented surface $σ$ , the "counter-clockwise oriented boundary" of $σ$ is a path $\partial σ$ that traces the boundary of $σ$ in a counter-clockwise direction. The counter-clockwise direction is better described as follows: While located on the boundary, the counter-clockwise direction is the "forwards" direction when the surface normal vectors point "up" and the surface itself is on the "left". The image below gives several examples of oriented surfaces and their counter-clockwise boundaries. Note in particular the 4th panel that shows that the orientation around a hole in the surface appears to be clockwise.

Given a multi-surface $𝐒 = {(σ_{1}, w_{1}), (σ_{2}, w_{2}), . . ., (σ_{k}, w_{k})}$ , the counter-clockwise boundary of $𝐒$ is the multi-path ${(\partial σ_{1}, w_{1}), (\partial σ_{2}, w_{2}), . . ., (\partial σ_{k}, w_{k})}$ .

Given a vector field $𝐅$ that denotes a multi-surface, the multi-path that denotes the counter-clockwise boundary of $𝐅$ is denoted by vector field $\nabla \times 𝐅$ . The vector field $\nabla \times 𝐅$ evaluated at point $𝐪$ is denoted by $\nabla \times 𝐅 (𝐪)$ , $(\nabla \times 𝐅) (𝐪)$ , or $\nabla \times 𝐅 |_{𝐪}$ .

The requirement that no surface weight extends to infinity means that all counter-clockwise boundaries form closed loops, and therefore the total displacement of the total counter-clockwise boundary is $𝟎$ : $\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \times 𝐅 (𝐪)) d V = 𝟎$ .

It is also important to note that the counter-clockwise boundary has no endpoints: $\nabla \cdot (\nabla \times 𝐅) = 0$ .

The similarity of the notation $\nabla \times 𝐅_{2}$ to the intersection of multi-surface $𝐅_{1}$ with multi-surface $𝐅_{2}$ , denoted by $𝐅_{1} \times 𝐅_{2}$ , again makes sense if $\nabla$ is interpreted as the "surface of reality". An edge is formed when a surface "slices" into reality.

In the image to the right, a depiction of the "surface of reality" interpretation of $\nabla$ is shown. On the right is a simple surface $𝐅$ , along with its counter-clockwise boundary $\nabla \times 𝐅$ . On the left $𝐅_{ext}$ is an extension of $𝐅$ that is behind the "veil" of surface $𝐆_{\nabla}$ . $𝐅_{ext}$ slices into $𝐆_{\nabla}$ at curves consistent with the boundary of $𝐅$ : i.e. $𝐆_{\nabla} \times 𝐅_{ext} = \nabla \times 𝐅$ .

The inwards-oriented surfaces of volumes

Given a volume $Ω$ , the "inwards oriented surface" of $Ω$ is a surface $\partial Ω$ that wraps the volume $Ω$ with the surface normals pointing inwards. The image below gives several examples of volumes and their inwards oriented surfaces.

Given a multi-volume $𝐔 = {(Ω_{1}, w_{1}), (Ω_{2}, w_{2}), . . ., (Ω_{k}, w_{k})}$ , the inwards oriented surface of $𝐔$ is the multi-surface ${(\partial Ω_{1}, w_{1}), (\partial Ω_{2}, w_{2}), . . ., (\partial Ω_{k}, w_{k})}$ .

Given a scalar field $U$ that denotes a multi-volume, the multi-surface that denotes the inwards oriented surface of $U$ is denoted by vector field $\nabla U$ . The vector field $\nabla U$ evaluated at point $𝐪$ is denoted by $\nabla U (𝐪)$ , $(\nabla U) (𝐪)$ , or $\nabla U |_{𝐪}$ .

The requirement that no volume weight extends to infinity means that all inwards oriented surfaces form closed surfaces, and therefore the total surface vector of the total inwards oriented surface is $𝟎$ : $\underset{𝐪 \in ℝ^{3}}{∭} (\nabla U (𝐪)) d V = 𝟎$ .

It is also important to note that the inwards oriented surface has no boundary: $\nabla \times (\nabla U) = 𝟎$ .

The similarity of the notation $\nabla U$ to the intersection of multi-surface $𝐅$ with multi-volume $U$ , denoted by $𝐅 U$ , again makes sense if $\nabla$ is interpreted as the "surface of reality". A surface is formed from the surface of reality when the volume "pushes" into reality.

In the image to the right, a depiction of the "surface of reality" interpretation of $\nabla$ is shown. The image is a 2D cross-section for simplicity. On the right is a simple volume $U$ , along with its inwards oriented surface $\nabla U$ . On the left $U_{ext}$ is an extension of $U$ that is behind the "veil" of surface $𝐆_{\nabla}$ . $U_{ext}$ pushes through $𝐆_{\nabla}$ at surfaces consistent with the surface of $U$ : i.e. $𝐆_{\nabla} U_{ext} = \nabla U$ .

Closed loops and closed surfaces

A simple path is "closed" or a "loop" if its starting and final points are the same, so the total endpoints is 0 since the weights of the starting and final points cancel out. More generally, a multi-path $𝐉$ is "closed" or a "multi-loop" if $\nabla \cdot 𝐉 = 0$ . As previously noted, the counter-clockwise boundary of a surface is closed.

A simple surface is "closed" or a "bubble" if it has no boundary. More generally, a multi-surface $𝐅$ is "closed" or a "multi-bubble" if $\nabla \times 𝐅 = 𝟎$ . As previously noted, the inwards oriented surface of a volume is closed.

It is clear that the total displacement present in a closed multi-path is $𝟎$ : $\nabla \cdot 𝐉 = 0 ⟹ \underset{𝐪 \in ℝ^{3}}{∭} 𝐉 d V = 𝟎$ , and it is also clear that the total surface vector of a closed multi-surface is also $𝟎$ : $\nabla \times 𝐅 = 𝟎 ⟹ \underset{𝐪 \in ℝ^{3}}{∭} 𝐅 d V = 𝟎$ .

Given both a simple loop and a simple bubble, the total point weight of all intersection points is 0: every time the loop enters the bubble, it must also leave the bubble, and the weights of these two intersection points cancel out. More generally, given a closed multi-path $𝐉$ and a closed multi-surface $𝐅$ , then the total intersection point weight is 0: $(\nabla \cdot 𝐉 = 0 and \nabla \times 𝐅 = 𝟎) ⟹ \underset{𝐪 \in ℝ^{3}}{∭} (𝐉 \cdot 𝐅) d V = 0$ .

The above identity gives rise to the following observations:

The total intersection point weight of a multi-loop and a multi-surface is purely a function of the multi-loop and the multi-surface's counter-clockwise boundary: the interior of the multi-surface does not matter. If $\nabla \cdot 𝐉 = 0$ and $\nabla \times 𝐅_{1} = \nabla \times 𝐅_{2}$ , then $\underset{𝐪 \in ℝ^{3}}{∭} (𝐉 \cdot 𝐅_{1}) d V = \underset{𝐪 \in ℝ^{3}}{∭} (𝐉 \cdot 𝐅_{2}) d V$ .
The total intersection point weight of a multi-path and a multi-bubble is purely a function of the multi-bubble and the multi-path's endpoints: the interior of the multi-path does not matter. If $\nabla \times 𝐅 = 𝟎$ and $\nabla \cdot 𝐉_{1} = \nabla \cdot 𝐉_{2}$ , then $\underset{𝐪 \in ℝ^{3}}{∭} (𝐉_{1} \cdot 𝐅) d V = \underset{𝐪 \in ℝ^{3}}{∭} (𝐉_{2} \cdot 𝐅) d V$ .

The inwards oriented surface of a volume is closed. Conversely, given a closed surface, there exists a volume that "fills" the surface. More generally, given a multi-bubble $𝐅$ , there exists a multi-volume $U$ for which $𝐅$ is the inwards oriented multi-surface of $U$ : $\nabla \times 𝐅 = 𝟎 ⟹ \exists U : \nabla U = 𝐅$ . This multi-volume is referred to as the "scalar potential" of $𝐅$ . The requirement that volumes cannot extend to infinity means that $U$ is unique.

The counter-clockwise oriented boundary of a surface is closed. Conversely, given a loop, there exists a surface that "fills" the loop. More generally, given a multi-loop $𝐉$ , there exists a multi-surface $𝐅$ for which $𝐉$ is the counter-clockwise boundary of $𝐅$ : $\nabla \cdot 𝐉 = 0 ⟹ \exists 𝐅 : \nabla \times 𝐅 = 𝐉$ . This multi-surface is referred to as the "vector potential" of $𝐉$ . Even with the requirement that surfaces cannot extend to infinity, $𝐅$ is not unique without additional restrictions.

Coordinate Systems

This section will describe how to compute various quantities such as intersections, endpoints, boundaries, and surfaces given a curvilinear coordinate system.

Let the curvilinear coordinate system be arbitrary. Let the 3 coordinates that index all points be $c_{1}, c_{2}, c_{3}$ . Coordinates will be denoted by the triple $(c_{1}, c_{2}, c_{3})$ .

The following notation will be used in the following discussions:

Given an arbitrary expression $f : {1, 2, 3} \to ℝ$ that assigns a real number to each index $i = 1, 2, 3$ , then $(i; f (i))$ will denote the triple $(f (1), f (2), f (3))$ .

Given index variables $i, j \in {1, 2, 3}$ , the expression $𝟏 (i = j)$ equals 1 if $i = j$ and 0 if otherwise.

Given an arbitrary expression $f : {1, 2, 3} \to ℝ$ that assigns a real number to each index $i = 1, 2, 3$ , then $\sum_{i} f (i)$ will denote the sum $f (1) + f (2) + f (3)$ .

Given an index variable $i \in {1, 2, 3}$ , $i + 1$ will rotate $i$ forwards by 1, and $i + 2$ will rotate $i$ forwards by 2. In essence, $i + 1 = {\begin{matrix} i + 1 & (i = 1, 2) \\ 1 & (i = 3) \end{matrix}$ and $i + 2 = {\begin{matrix} 3 & (i = 1) \\ i - 1 & (i = 2, 3) \end{matrix}$ .

Start with an arbitrary coordinate $(c'_{1}, c'_{2}, c'_{3}) = (j; c'_{j})$ and infinitesimal differences $Δ c_{1}$ , $Δ c_{2}$ , and $Δ c_{3}$ . The following 3 paths, 3 surfaces, and volume will be associated with point $(j; c'_{j})$ :

For each $i \in {1, 2, 3}$ there exists an infinitely short path $C_{i} ((j; c'_{j}))$ starting from point $(j; c'_{j})$ and ending on point $(j; c'_{j} + Δ c_{i} 𝟏 (j = i))$ along the curve defined by $c'_{i} \leq c_{i} < c'_{i} + Δ c_{i}$ , $c_{i + 1} = c'_{i + 1}$ and $c_{i + 2} = c'_{i + 2}$ . The displacement covered by $C_{i} ((j; c'_{j}))$ is approximately $Δ c_{i} \cdot l_{i} ((j; c'_{j})) \cdot {\hat{𝐚}}_{i} ((j; c'_{j}))$ where ${\hat{𝐚}}_{i} ((j; c'_{j}))$ is a unit length vector that is parallel to the displacement between points $(j; c'_{j})$ and $(j; c'_{j} + Δ c_{i} 𝟏 (j = i))$ , and $Δ c_{i} \cdot l_{i} ((j; c'_{j}))$ is the length of the displacement. Note that the length of the displacement is proportional to $Δ c_{i}$ , with $l_{i} ((j; c'_{j}))$ being the constant of proportionality. The set of vectors ${{\hat{𝐚}}_{1} ((j; c'_{j})), {\hat{𝐚}}_{2} ((j; c'_{j})), {\hat{𝐚}}_{3} ((j; c'_{j}))}$ is the set of displacement basis vectors.

For each $i \in {1, 2, 3}$ there exists an infinitely small surface $σ_{i} ((j; c'_{j}))$ that is defined by the following: $c_{i} = c'_{i}$ , $c'_{i + 1} \leq c_{i + 1} < c'_{i + 1} + Δ c_{i + 1}$ , and $c'_{i + 2} \leq c_{i + 2} < c'_{i + 2} + Δ c_{i + 2}$ . The orientation of $σ_{i} ((j; c'_{j}))$ is in the direction of increasing $c_{i}$ . The surface vector of $σ_{i} ((j; c'_{j}))$ is approximately $Δ c_{i + 1} Δ c_{i + 2} \cdot A_{i} ((j; c'_{j})) \cdot {\hat{𝐚}}^{i} ((j; c'_{j}))$ where ${\hat{𝐚}}^{i} ((j; c'_{j}))$ is a unit length vector that is perpendicular to $σ_{i} ((j; c'_{j}))$ , and $Δ c_{i + 1} Δ c_{i + 2} \cdot A_{i} ((j; c'_{j}))$ is the area of $σ_{i} ((j; c'_{j}))$ . Note that the area of $σ_{i} ((j; c'_{j}))$ is proportional to $Δ c_{i + 1} Δ c_{i + 2}$ , with $A_{i} ((j; c'_{j}))$ being the constant of proportionality. The set of vectors ${{\hat{𝐚}}^{1} ((j; c'_{j})), {\hat{𝐚}}^{2} ((j; c'_{j})), {\hat{𝐚}}^{3} ((j; c'_{j}))}$ is the set of surface basis vectors.

There is an infinitely small volume $Ω ((j; c'_{j}))$ defined by $c'_{1} \leq c_{1} < c'_{1} + Δ c_{1}$ , $c'_{2} \leq c_{2} < c'_{2} + Δ c_{2}$ , and $c'_{3} \leq c_{3} < c'_{3} + Δ c_{3}$ . $Ω ((j; c'_{j}))$ has a shape that is approximately that of a parallelepiped. The volume of $Ω ((j; c'_{j}))$ is approximately $Δ c_{1} Δ c_{2} Δ c_{3} \cdot V ((j; c'_{j}))$ . Note that the volume of $Ω ((j; c'_{j}))$ is proportional to $Δ c_{1} Δ c_{2} Δ c_{3}$ , with $V ((j; c'_{j}))$ being the constant of proportionality.

It is important to note that:

$(i; c_{i}) \in Ω ((j; c'_{j}))$ if and only if $c'_{1} \leq c_{1} < c'_{1} + Δ c_{1}$ , $c'_{2} \leq c_{2} < c'_{2} + Δ c_{2}$ , and $c'_{3} \leq c_{3} < c'_{3} + Δ c_{3}$ (note the strictness of the upper bounds).
For all $i \in {1, 2, 3}$ , $C_{i} ((j; c_{j})) \subseteq Ω ((j; c'_{j}))$ if and only if $c_{i} = c'_{i}$ , $c'_{i + 1} \leq c_{i + 1} < c'_{i + 1} + Δ c_{i + 1}$ , and $c'_{i + 2} \leq c_{i + 2} < c'_{i + 2} + Δ c_{i + 2}$ (note the strictness of the upper bounds).
For all $i \in {1, 2, 3}$ , $σ_{i} ((j; c_{j})) \subseteq Ω ((j; c'_{j}))$ if and only if $c'_{i} \leq c_{i} < c'_{i} + Δ c_{i}$ (note the strictness of the upper bound), $c_{i + 1} = c'_{i + 1}$ , and $c_{i + 2} = c'_{i + 2}$ .

Converting between multi-points, multi-paths, multi-surfaces, and multi-volumes and their respective scalar fields and vector fields proceeds as follows:

This conversion is performed by subdividing space into discrete volumes or cells. Infinitesimal differences $Δ c_{1}$ , $Δ c_{2}$ , and $Δ c_{3}$ are chosen, and a lattice consisting of the points $(j; k_{j} Δ c_{j})$ where $(j; k_{j})$ is an arbitrary triple of integers is generated. The cell indexed by $(j; k_{j})$ consists of the point $(j; k_{j} Δ c_{j})$ , the paths $C_{i} ((j; k_{j} Δ c_{j}))$ for each $i \in {1, 2, 3}$ , the surfaces $σ_{i} ((j; k_{j} Δ c_{j}))$ for each $i \in {1, 2, 3}$ , and the volume $Ω ((j; k_{j} Δ c_{j}))$ . All points $(j; c_{j})$ where $k_{i} Δ c_{i} \leq c_{i} < (k_{i} + 1) Δ c_{i}$ for all $i \in {1, 2, 3}$ "belong" to the cell indexed by $(j; k_{j})$ (note that the upper bounds are excluded). Given an arbitrary point $(j; c_{j})$ , the cell that contains $(j; c_{j})$ is indexed by $(j; k_{j}) = (j; ⌊ \frac{c_{j}}{Δ c_{j}} ⌋)$ . The point $(j; c'_{j}) = (j; k_{j} Δ c_{j})$ is the vertex that the cell is associated with.

A multi-point, multi-path, multi-surface, or multi-volume is converted to a scalar field or vector field by computing the total point weight, displacement, surface vector, or volume contained by each cell and then averaging over the cell's volume.

A scalar-field $ρ$ is converted to a multi-point by doing the following for each cell $(j; k_{j})$ . First compute the total point weight contained inside the cell: $\underset{𝐪 \in Ω ((j; k_{j} Δ c_{j}))}{∭} ρ (𝐪) d V \approx ρ ((j; k_{j} Δ c_{j})) V ((j; k_{j} Δ c_{j})) Δ c_{1} Δ c_{2} Δ c_{3}$ . Next assign this weight to the point $(j; k_{j} Δ c_{j})$ .

A vector-field $𝐉 = \sum_{i} J_{i} {\hat{𝐚}}_{i}$ is converted to a multi-path by doing the following for each cell $(j; k_{j})$ . First compute the total displacement contained inside the cell: $\underset{𝐪 \in Ω ((j; k_{j} Δ c_{j}))}{∭} 𝐉 (𝐪) d V \approx (\sum_{i} J_{i} ((j; k_{j} Δ c_{j})) {\hat{𝐚}}_{i} ((j; k_{j} Δ c_{j}))) V ((j; k_{j} Δ c_{j})) Δ c_{1} Δ c_{2} Δ c_{3}$ . Next separate this total displacement into components according to the basis ${\hat{𝐚}}_{1}$ , ${\hat{𝐚}}_{2}$ , and ${\hat{𝐚}}_{3}$ : for each $i \in {1, 2, 3}$ the coefficient of ${\hat{𝐚}}_{i}$ is $\underset{𝐪 \in Ω ((j; k_{j} Δ c_{j}))}{∭} J_{i} (𝐪) d V \approx J_{i} ((j; k_{j} Δ c_{j})) V ((j; k_{j} Δ c_{j})) Δ c_{1} Δ c_{2} Δ c_{3}$ . Next for each $i \in {1, 2, 3}$ , divide the coefficient of ${\hat{𝐚}}_{i}$ by the length of $C_{i} ((j; k_{j} Δ c_{j}))$ , which results in approximately $J_{i} ((j; k_{j} Δ c_{j})) \frac{V ((j; k_{j} Δ c_{j}))}{l_{i} ((j; k_{j} Δ c_{j}))} Δ c_{i + 1} Δ c_{i + 2}$ , and assign this weight to $C_{i} ((j; k_{j} Δ c_{j}))$ .

A vector-field $𝐅 = \sum_{i} F_{i} {\hat{𝐚}}^{i}$ is converted to a multi-surface by doing the following for each cell $(j; k_{j})$ . First compute the total surface vector contained inside the cell: $\underset{𝐪 \in Ω ((j; k_{j} Δ c_{j}))}{∭} 𝐅 (𝐪) d V \approx (\sum_{i} F_{i} ((j; k_{j} Δ c_{j})) {\hat{𝐚}}^{i} ((j; k_{j} Δ c_{j}))) V ((j; k_{j} Δ c_{j})) Δ c_{1} Δ c_{2} Δ c_{3}$ . Next separate this total surface vector into components according to the basis ${\hat{𝐚}}^{1}$ , ${\hat{𝐚}}^{2}$ , and ${\hat{𝐚}}^{3}$ : for each $i \in {1, 2, 3}$ the coefficient of ${\hat{𝐚}}^{i}$ is $\underset{𝐪 \in Ω ((j; k_{j} Δ c_{j}))}{∭} F_{i} (𝐪) d V \approx F_{i} ((j; k_{j} Δ c_{j})) V ((j; k_{j} Δ c_{j})) Δ c_{1} Δ c_{2} Δ c_{3}$ . Next for each $i \in {1, 2, 3}$ , divide the coefficient of ${\hat{𝐚}}^{i}$ by the area of $σ_{i} ((j; k_{j} Δ c_{j}))$ , which results in approximately $F_{i} ((j; k_{j} Δ c_{j})) \frac{V ((j; k_{j} Δ c_{j}))}{A_{i} ((j; k_{j} Δ c_{j}))} Δ c_{i}$ , and assign this weight to $σ_{i} ((j; k_{j} Δ c_{j}))$ .

A scalar-field $U$ is converted to a multi-volume by doing the following for each cell $(j; k_{j})$ . First compute the total volume contained inside the cell: $\underset{𝐪 \in Ω ((j; k_{j} Δ c_{j}))}{∭} U (𝐪) d V \approx U ((j; k_{j} Δ c_{j})) V ((j; k_{j} Δ c_{j})) Δ c_{1} Δ c_{2} Δ c_{3}$ . Next divide this weight by the volume of $Ω ((j; k_{j} Δ c_{j}))$ , which results in approximately $U ((j; k_{j} Δ c_{j}))$ , and assign this weight to $Ω ((j; k_{j} Δ c_{j}))$ .

Computing various intersections

Computing the intersection of any structure with a multi-volume is trivial matter: Simply multiply the scalar of vector field by the scalar field that denotes the multi-volume. When both structures are denoted by vector fields however, computing the intersection is far less trivial.

Computing path-surface intersections

To save space, the notation $(j; c_{j})$ and $(j; k_{j} Δ c_{j})$ will be omitted from the various terms.

Given a multi-path $𝐂$ denoted by vector field $𝐉 = \sum_{i} J_{i} {\hat{𝐚}}_{i}$ , and a multi-surface $𝐒$ denoted by vector field $𝐅 = \sum_{i} F_{i} {\hat{𝐚}}^{i}$ , the scalar field that denotes the intersection can be computed as follows:

The following computations applies to each cell:

For each $i \in {1, 2, 3}$ , the weight assigned to $C_{i}$ by $𝐂$ is computed as follows: $Δ c_{1} Δ c_{2} Δ c_{3} \cdot V \cdot J_{i}$ is the ${\hat{𝐚}}_{i}$ component of the total displacement contained by the current cell. Computing the weight assigned to $C_{i}$ requires that this displacement be spread over the length of $C_{i}$ : $\frac{Δ c_{1} Δ c_{2} Δ c_{3} \cdot V \cdot J_{i}}{Δ c_{i} \cdot l_{i}} = \frac{V}{l_{i}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot J_{i}$ .

For each $i \in {1, 2, 3}$ , the weight assigned to $σ_{i}$ by $𝐒$ is computed as follows: $Δ c_{1} Δ c_{2} Δ c_{3} \cdot V \cdot F_{i}$ is the ${\hat{𝐚}}^{i}$ component of the total surface vector contained by the current cell. Computing the weight assigned to $σ_{i}$ requires that this surface vector be spread over the area of $σ_{i}$ : $\frac{Δ c_{1} Δ c_{2} Δ c_{3} \cdot V \cdot F_{i}}{Δ c_{i + 1} Δ c_{i + 2} \cdot A_{i}} = \frac{V}{A_{i}} \cdot Δ c_{i} \cdot F_{i}$ .

The intersection between $C_{i}$ and $σ_{i}$ is the current lattice point with weight $(\frac{V}{l_{i}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot J_{i}) (\frac{V}{A_{i}} \cdot Δ c_{i} \cdot F_{i})$ $= \frac{V^{2}}{l_{i} A_{i}} \cdot Δ c_{1} Δ c_{2} Δ c_{3} \cdot J_{i} F_{i}$ .

Aside from the intersections between $C_{i}$ and $σ_{i}$ for each cell and $i \in {1, 2, 3}$ , no other intersections occur. The total weight of the intersection at the vertex of the current cell is $\sum_{i} \frac{V^{2}}{l_{i} A_{i}} \cdot Δ c_{1} Δ c_{2} Δ c_{3} \cdot J_{i} F_{i}$ $= V^{2} \cdot Δ c_{1} Δ c_{2} Δ c_{3} \sum_{i} \frac{1}{l_{i} A_{i}} \cdot J_{i} F_{i}$ .

The value of $𝐉 \cdot 𝐅$ at the current cell is approximately $\frac{1}{Δ c_{1} Δ c_{2} Δ c_{3} \cdot V} \cdot V^{2} \cdot Δ c_{1} Δ c_{2} Δ c_{3} \sum_{i} \frac{1}{l_{i} A_{i}} \cdot J_{i} F_{i}$ $= V \sum_{i} \frac{1}{l_{i} A_{i}} \cdot J_{i} F_{i}$ . The coefficient of $\frac{1}{Δ c_{1} Δ c_{2} Δ c_{3} \cdot V}$ exists to spread the point weight over the current cell.

Therefore $𝐉 \cdot 𝐅 = V \sum_{i} \frac{1}{l_{i} A_{i}} \cdot J_{i} F_{i}$ . Note that $𝐉$ is expressed using the displacement basis vectors, while $𝐅$ is expressed using the surface basis vectors.

Computing surface-surface intersections

To save space, the notation $(j; c_{j})$ and $(j; k_{j} Δ c_{j})$ will be omitted from the various terms.

Given a multi-surface $𝐒_{1}$ denoted by vector field $𝐅 = \sum_{i} F_{i} {\hat{𝐚}}^{i}$ , and a multi-surface $𝐒_{2}$ denoted by vector field $𝐆 = \sum_{i} G_{i} {\hat{𝐚}}^{i}$ , the vector field that denotes the intersection can be computed as follows:

The following computations applies to each cell:

For each $i \in {1, 2, 3}$ , the weight assigned to $σ_{i}$ by $𝐒_{1}$ is computed as follows: $Δ c_{1} Δ c_{2} Δ c_{3} \cdot V \cdot F_{i}$ is the ${\hat{𝐚}}^{i}$ component of the total surface vector contained by the current cell. Computing the weight assigned to $σ_{i}$ requires that this surface vector be spread over the area of $σ_{i}$ : $\frac{Δ c_{1} Δ c_{2} Δ c_{3} \cdot V \cdot F_{i}}{Δ c_{i + 1} Δ c_{i + 2} \cdot A_{i}} = \frac{V}{A_{i}} \cdot Δ c_{i} \cdot F_{i}$ . Similarly, the weight assigned to $σ_{i}$ by $𝐒_{2}$ is $\frac{V}{A_{i}} \cdot Δ c_{i} \cdot G_{i}$ .

The intersection between $σ_{i + 1}$ and $σ_{i + 2}$ is path $C_{i}$ with weight $(\frac{V}{A_{i + 1}} \cdot Δ c_{i + 1} \cdot F_{i + 1}) (\frac{V}{A_{i + 2}} \cdot Δ c_{i + 2} \cdot G_{i + 2})$ $= \frac{V^{2}}{A_{i + 1} A_{i + 2}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot F_{i + 1} G_{i + 2}$ . Conversely, the intersection between $σ_{i + 2}$ and $σ_{i + 1}$ is path $C_{i}$ with weight $- \frac{V^{2}}{A_{i + 2} A_{i + 1}} \cdot Δ c_{i + 2} Δ c_{i + 1} \cdot F_{i + 2} G_{i + 1}$ .

Aside from the intersections between $σ_{i + 1}$ and $σ_{i + 2}$ , and the intersections between $σ_{i + 2}$ and $σ_{i + 1}$ , for each cell and $i \in {1, 2, 3}$ , no other intersections occur. For each $i \in {1, 2, 3}$ , the total weight assigned to $C_{i}$ is $\frac{V^{2}}{A_{i + 1} A_{i + 2}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot F_{i + 1} G_{i + 2} - \frac{V^{2}}{A_{i + 2} A_{i + 1}} \cdot Δ c_{i + 2} Δ c_{i + 1} \cdot F_{i + 2} G_{i + 1}$ $= \frac{V^{2}}{A_{i + 1} A_{i + 2}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot (F_{i + 1} G_{i + 2} - F_{i + 2} G_{i + 1})$ .

The value of $𝐅 \times 𝐆$ at the current cell is approximately $\sum_{i} \frac{l_{i} \cdot Δ c_{i} \cdot {\hat{𝐚}}_{i}}{V \cdot Δ c_{1} Δ c_{2} Δ c_{3}} \cdot \frac{V^{2}}{A_{i + 1} A_{i + 2}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot (F_{i + 1} G_{i + 2} - F_{i + 2} G_{i + 1})$ $= \sum_{i} \frac{V l_{i}}{A_{i + 1} A_{i + 2}} \cdot (F_{i + 1} G_{i + 2} - F_{i + 2} G_{i + 1}) {\hat{𝐚}}_{i}$ . The coefficient of $\frac{l_{i} \cdot Δ c_{i} \cdot {\hat{𝐚}}_{i}}{V \cdot Δ c_{1} Δ c_{2} Δ c_{3}}$ exists to spread the displacement of each path over the current cell.

Therefore $𝐅 \times 𝐆 = \sum_{i} \frac{V l_{i}}{A_{i + 1} A_{i + 2}} \cdot (F_{i + 1} G_{i + 2} - F_{i + 2} G_{i + 1}) {\hat{𝐚}}_{i}$ . Note that both $𝐅$ and $𝐆$ are both expressed using the surface basis vectors, but $𝐅 \times 𝐆$ is using the displacement basis vectors.

Computing the endpoints of paths

To save space, the notation $(j; k_{j})$ , $(j; c_{j})$ , and $(j; k_{j} Δ c_{j})$ will be omitted from the various terms. However, given a quantity $Q$ and an arbitrary $i \in {1, 2, 3}$ , the notation $[Q]_{- i}$ will denote the quantity in the adjacent cell by moving back one step along the dimension indexed by $i$ . This cell will be referred to as the $- i$ neighbor of the current cell.

Given a multi-path $𝐂$ denoted by vector field $𝐉 = \sum_{i} J_{i} {\hat{𝐚}}_{i}$ , the scalar field that denotes the endpoints can be computed as follows:

The following computations apply to each cell:

For each $i \in {1, 2, 3}$ , the weight assigned to $C_{i}$ by $𝐂$ is computed as follows: $Δ c_{1} Δ c_{2} Δ c_{3} \cdot V \cdot J_{i}$ is the ${\hat{𝐚}}_{i}$ component of the total displacement contained by the current cell. Computing the weight assigned to $C_{i}$ requires that this displacement be spread over the length of $C_{i}$ : $\frac{Δ c_{1} Δ c_{2} Δ c_{3} \cdot V \cdot J_{i}}{Δ c_{i} \cdot l_{i}} = \frac{V}{l_{i}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot J_{i}$ .

For each $i \in {1, 2, 3}$ , path $C_{i}$ contributes a weight of $+ \frac{V}{l_{i}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot J_{i}$ to the lattice point of the current cell, and path $[C_{i}]_{- i}$ contributes a weight of $- {[\frac{V}{l_{i}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot J_{i}]}_{- i}$ to the lattice point of the current cell.

The total weight assigned to the lattice point of the current cell is $\sum_{i} (+ \frac{V}{l_{i}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot J_{i} - {[\frac{V}{l_{i}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot J_{i}]}_{- i})$ $\approx \sum_{i} Δ c_{i} \frac{\partial}{\partial c_{i}} (\frac{V}{l_{i}} \cdot Δ c_{i + 1} Δ c_{i + 2} \cdot J_{i})$ $= Δ c_{1} Δ c_{2} Δ c_{3} \sum_{i} \frac{\partial}{\partial c_{i}} (\frac{V}{l_{i}} \cdot J_{i})$ .

Spreading the weight assigned to the current lattice point over the volume of the current cell gives: $\nabla \cdot 𝐉 = \frac{1}{Δ c_{1} Δ c_{2} Δ c_{3} \cdot V} \cdot Δ c_{1} Δ c_{2} Δ c_{3} \sum_{i} \frac{\partial}{\partial c_{i}} (\frac{V}{l_{i}} \cdot J_{i})$ $= \frac{1}{V} \sum_{i} \frac{\partial}{\partial c_{i}} (\frac{V}{l_{i}} \cdot J_{i})$ .

Therefore: $\nabla \cdot 𝐉 = \frac{1}{V} \sum_{i} \frac{\partial}{\partial c_{i}} (\frac{V}{l_{i}} \cdot J_{i})$ . Note that $𝐉$ is expressed using the displacement basis vectors.

Computing the counter-clockwise boundaries of surfaces

To save space, the notation $(j; k_{j})$ , $(j; c_{j})$ , and $(j; k_{j} Δ c_{j})$ will be omitted from the various terms. However, given a quantity $Q$ and an arbitrary $i \in {1, 2, 3}$ , the notation $[Q]_{- i}$ will denote the quantity in the adjacent cell by moving back one step along the dimension indexed by $i$ . This cell will be referred to as the $- i$ neighbor of the current cell.

Given a multi-surface $𝐒$ denoted by vector field $𝐅 = \sum_{i} F_{i} {\hat{𝐚}}^{i}$ , the vector field that denotes the counter-clockwise boundary can be computed as follows:

The following computations apply to each cell:

For each $i \in {1, 2, 3}$ , the weight assigned to $σ_{i}$ by $𝐒$ is computed as follows: $Δ c_{1} Δ c_{2} Δ c_{3} \cdot V \cdot F_{i}$ is the ${\hat{𝐚}}^{i}$ component of the total surface vector contained by the current cell. Computing the weight assigned to $σ_{i}$ requires that this surface vector be spread over the area of $σ_{i}$ : $\frac{Δ c_{1} Δ c_{2} Δ c_{3} \cdot V \cdot F_{i}}{Δ c_{i + 1} Δ c_{i + 2} \cdot A_{i}} = \frac{V}{A_{i}} \cdot Δ c_{i} \cdot F_{i}$ .

For each $i \in {1, 2, 3}$ , surfaces that contain path $C_{i}$ as part of their boundary include $σ_{i + 1}$ , $[σ_{i + 1}]_{- (i + 2)}$ , $σ_{i + 2}$ , and $[σ_{i + 2}]_{- (i + 1)}$ . $C_{i}$ receives a mass of $- \frac{V}{A_{i + 1}} \cdot Δ c_{i + 1} \cdot F_{i + 1}$ from $σ_{i + 1}$ ; a mass of $+ {[\frac{V}{A_{i + 1}} \cdot Δ c_{i + 1} \cdot F_{i + 1}]}_{- (i + 2)}$ from $[σ_{i + 1}]_{- (i + 2)}$ ; a mass of $+ \frac{V}{A_{i + 2}} \cdot Δ c_{i + 2} \cdot F_{i + 2}$ from $σ_{i + 2}$ ; and a mass of $- {[\frac{V}{A_{i + 2}} \cdot Δ c_{i + 2} \cdot F_{i + 2}]}_{- (i + 1)}$ from $[σ_{i + 2}]_{- (i + 1)}$ . The total mass assigned to $C_{i}$ is $- \frac{V}{A_{i + 1}} \cdot Δ c_{i + 1} \cdot F_{i + 1} + {[\frac{V}{A_{i + 1}} \cdot Δ c_{i + 1} \cdot F_{i + 1}]}_{- (i + 2)} + \frac{V}{A_{i + 2}} \cdot Δ c_{i + 2} \cdot F_{i + 2} - {[\frac{V}{A_{i + 2}} \cdot Δ c_{i + 2} \cdot F_{i + 2}]}_{- (i + 1)}$ $\approx - Δ c_{i + 2} \frac{\partial}{\partial c_{i + 2}} (\frac{V}{A_{i + 1}} \cdot Δ c_{i + 1} \cdot F_{i + 1}) + Δ c_{i + 1} \cdot \frac{\partial}{\partial c_{i + 1}} (\frac{V}{A_{i + 2}} \cdot Δ c_{i + 2} \cdot F_{i + 2})$ $= Δ c_{i + 1} Δ c_{i + 2} (\frac{\partial}{\partial c_{i + 1}} (\frac{V}{A_{i + 2}} \cdot F_{i + 2}) - \frac{\partial}{\partial c_{i + 2}} (\frac{V}{A_{i + 1}} \cdot F_{i + 1}))$ .

Spreading the displacement generated by each $C_{i}$ over the volume of the current cell gives: $\nabla \times 𝐅 = \sum_{i} \frac{Δ c_{i} \cdot l_{i} \cdot {\hat{𝐚}}_{i}}{Δ c_{1} Δ c_{2} Δ c_{3} \cdot V} \cdot Δ c_{i + 1} Δ c_{i + 2} (\frac{\partial}{\partial c_{i + 1}} (\frac{V}{A_{i + 2}} \cdot F_{i + 2}) - \frac{\partial}{\partial c_{i + 2}} (\frac{V}{A_{i + 1}} \cdot F_{i + 1}))$ $= \sum_{i} \frac{l_{i}}{V} (\frac{\partial}{\partial c_{i + 1}} (\frac{V}{A_{i + 2}} \cdot F_{i + 2}) - \frac{\partial}{\partial c_{i + 2}} (\frac{V}{A_{i + 1}} \cdot F_{i + 1})) {\hat{𝐚}}_{i}$ .

Therefore: $\nabla \times 𝐅 = \sum_{i} \frac{l_{i}}{V} (\frac{\partial}{\partial c_{i + 1}} (\frac{V}{A_{i + 2}} \cdot F_{i + 2}) - \frac{\partial}{\partial c_{i + 2}} (\frac{V}{A_{i + 1}} \cdot F_{i + 1})) {\hat{𝐚}}_{i}$ . Note that $𝐅$ is expressed using the surface basis vectors, but $\nabla \times 𝐅$ is using the displacement basis vectors.

Computing the inwards-oriented surfaces of volumes

To save space, the notation $(j; k_{j})$ , $(j; c_{j})$ , and $(j; k_{j} Δ c_{j})$ will be omitted from the various terms. However, given a quantity $Q$ and an arbitrary $i \in {1, 2, 3}$ , the notation $[Q]_{- i}$ will denote the quantity in the adjacent cell by moving back one step along the dimension indexed by $i$ . This cell will be referred to as the $- i$ neighbor of the current cell.

Given a multi-volume $𝐔$ denoted by scalar field $U$ , the vector field that denotes the inwards-oriented surface can be computed as follows:

The following computations apply to each cell:

The cell's volume $Ω$ has the weight $U$ .

For each $i \in {1, 2, 3}$ , surface $σ_{i}$ receives a weight of $U$ from the current cell, and a weight of $- [U]_{- i}$ from the $- i$ neighbor of the current cell. The total weight is simply $U - [U]_{- i} \approx Δ c_{i} \frac{\partial U}{\partial c_{i}}$ . Spreading the surface vector generated by each $σ_{i}$ over the volume of the current cell gives: $\nabla U = \sum_{i} \frac{Δ c_{i + 1} Δ c_{i + 2} \cdot A_{i} \cdot {\hat{𝐚}}^{i}}{Δ c_{1} Δ c_{2} Δ c_{3} \cdot V} \cdot Δ c_{i} \frac{\partial U}{\partial c_{i}}$ $= \sum_{i} \frac{A_{i}}{V} \frac{\partial U}{\partial c_{i}} {\hat{𝐚}}^{i}$ .

Therefore: $\nabla U = \sum_{i} \frac{A_{i}}{V} \frac{\partial U}{\partial c_{i}} {\hat{𝐚}}^{i}$ . Note that $\nabla U$ uses the surface basis vectors.

Summary

Given multi-path $𝐉 = \sum_{i} J_{i} {\hat{𝐚}}_{i}$ and multi-surface $𝐅 = \sum_{i} F_{i} {\hat{𝐚}}^{i}$ , the intersection of $𝐉$ with $𝐅$ is multi-point $𝐉 \cdot 𝐅 = \sum_{i} \frac{V}{l_{i} A_{i}} J_{i} F_{i}$ .
Given multi-surfaces $𝐅 = \sum_{i} F_{i} {\hat{𝐚}}^{i}$ and $𝐆 = \sum_{i} G_{i} {\hat{𝐚}}^{i}$ , the intersection of $𝐅$ with $𝐆$ is multi-path $𝐅 \times 𝐆 = \sum_{i} \frac{l_{i} V}{A_{i + 1} A_{i + 2}} (F_{i + 1} G_{i + 2} - F_{i + 2} G_{i + 1}) {\hat{𝐚}}_{i}$ .
Given multi-path $𝐉 = \sum_{i} J_{i} {\hat{𝐚}}_{i}$ , the endpoints of $𝐉$ is multi-point $\nabla \cdot 𝐉 = \sum_{i} \frac{1}{V} \frac{\partial}{\partial c_{i}} (\frac{V}{l_{i}} J_{i})$ .
Given multi-surface $𝐅 = \sum_{i} F_{i} {\hat{𝐚}}^{i}$ , the counter-clockwise boundary of $𝐅$ is multi-path $\nabla \times 𝐅 = \sum_{i} \frac{l_{i}}{V} (\frac{\partial}{\partial c_{i + 1}} (\frac{V}{A_{i + 2}} F_{i + 2}) - \frac{\partial}{\partial c_{i + 2}} (\frac{V}{A_{i + 1}} F_{i + 1})) {\hat{𝐚}}_{i}$ .
Given multi-volume $U$ , the inwards-oriented surface of $U$ is multi-surface $\nabla U = \sum_{i} \frac{A_{i}}{V} \frac{\partial U}{\partial c_{i}} {\hat{𝐚}}^{i}$ .

Orthogonal coordinate systems

In the special case where the displacement basis vectors ${{\hat{𝐚}}_{1}, {\hat{𝐚}}_{2}, {\hat{𝐚}}_{3}}$ are all mutually orthogonal (perpendicular), then:

The surface basis vectors are identical to the displacement basis vectors: $\forall i \in {1, 2, 3} : {\hat{𝐚}}^{i} = {\hat{𝐚}}_{i}$ .
For each $i \in {1, 2, 3}$ , $A_{i} = l_{i + 1} l_{i + 2}$ .
$V = l_{1} l_{2} l_{3}$ .

The above formulas simplify to:

$(\sum_{i} J_{i} {\hat{𝐚}}_{i}) \cdot (\sum_{i} F_{i} {\hat{𝐚}}^{i}) = \sum_{i} J_{i} F_{i}$ .
$(\sum_{i} F_{i} {\hat{𝐚}}^{i}) \times (\sum_{i} G_{i} {\hat{𝐚}}^{i}) = \sum_{i} (F_{i + 1} G_{i + 2} - F_{i + 2} G_{i + 1}) {\hat{𝐚}}_{i}$ .
$\nabla \cdot (\sum_{i} J_{i} {\hat{𝐚}}_{i}) = \sum_{i} \frac{1}{l_{1} l_{2} l_{3}} \frac{\partial}{\partial c_{i}} (l_{i + 1} l_{i + 2} J_{i})$ .
$\nabla \times (\sum_{i} F_{i} {\hat{𝐚}}^{i}) = \sum_{i} \frac{1}{l_{i + 1} l_{i + 2}} (\frac{\partial}{\partial c_{i + 1}} (l_{i + 2} F_{i + 2}) - \frac{\partial}{\partial c_{i + 2}} (l_{i + 1} F_{i + 1})) {\hat{𝐚}}_{i}$ .
$\nabla U = \sum_{i} \frac{1}{l_{i}} \frac{\partial U}{\partial c_{i}} {\hat{𝐚}}^{i}$ .

For Cartesian coordinates, $c_{1} = x$ , $c_{2} = y$ , $c_{3} = z$ , and ${\hat{𝐚}}^{1} = {\hat{𝐚}}_{1} = \hat{𝐱}$ , ${\hat{𝐚}}^{2} = {\hat{𝐚}}_{2} = \hat{𝐲}$ , ${\hat{𝐚}}^{3} = {\hat{𝐚}}_{3} = \hat{𝐳}$ , and $l_{1} = 1$ , $l_{2} = 1$ , $l_{3} = 1$ . Therefore:

$(J_{x} \hat{𝐱} + J_{y} \hat{𝐲} + J_{z} \hat{𝐳}) \cdot (F_{x} \hat{𝐱} + F_{y} \hat{𝐲} + F_{z} \hat{𝐳}) = J_{x} F_{x} + J_{y} F_{y} + J_{z} F_{z}$ .
$(F_{x} \hat{𝐱} + F_{y} \hat{𝐲} + F_{z} \hat{𝐳}) \times (G_{x} \hat{𝐱} + G_{y} \hat{𝐲} + G_{z} \hat{𝐳}) = (F_{y} G_{z} - F_{z} G_{y}) \hat{𝐱} + (F_{z} G_{x} - F_{x} G_{z}) \hat{𝐲} + (F_{x} G_{y} - F_{y} G_{x}) \hat{𝐳}$ .
$\nabla \cdot (J_{x} \hat{𝐱} + J_{y} \hat{𝐲} + J_{z} \hat{𝐳}) = \frac{\partial J_{x}}{\partial x} + \frac{\partial J_{y}}{\partial y} + \frac{\partial J_{z}}{\partial z}$ .
$\nabla \times (F_{x} \hat{𝐱} + F_{y} \hat{𝐲} + F_{z} \hat{𝐳}) = (\frac{\partial F_{z}}{\partial y} - \frac{\partial F_{y}}{\partial z}) \hat{𝐱} + (\frac{\partial F_{x}}{\partial z} - \frac{\partial F_{z}}{\partial x}) \hat{𝐲} + (\frac{\partial F_{y}}{\partial x} - \frac{\partial F_{x}}{\partial y}) \hat{𝐳}$ .
$\nabla U = \frac{\partial U}{\partial x} \hat{𝐱} + \frac{\partial U}{\partial y} \hat{𝐲} + \frac{\partial U}{\partial z} \hat{𝐳}$ .

For cylindrical coordinates, $c_{1} = ρ$ , $c_{2} = ϕ$ , $c_{3} = z$ , and ${\hat{𝐚}}^{1} = {\hat{𝐚}}_{1} = \hat{ρ}$ , ${\hat{𝐚}}^{2} = {\hat{𝐚}}_{2} = \hat{ϕ}$ , ${\hat{𝐚}}^{3} = {\hat{𝐚}}_{3} = \hat{𝐳}$ , and $l_{1} = 1$ , $l_{2} = ρ$ , $l_{3} = 1$ . Therefore:

$(J_{ρ} \hat{ρ} + J_{ϕ} \hat{ϕ} + J_{z} \hat{𝐳}) \cdot (F_{ρ} \hat{ρ} + F_{ϕ} \hat{ϕ} + F_{z} \hat{𝐳}) = J_{ρ} F_{ρ} + J_{ϕ} F_{ϕ} + J_{z} F_{z}$ .
$(F_{ρ} \hat{ρ} + F_{ϕ} \hat{ϕ} + F_{z} \hat{𝐳}) \times (F_{ρ} \hat{ρ} + F_{ϕ} \hat{ϕ} + F_{z} \hat{𝐳}) = (F_{ϕ} G_{z} - F_{z} G_{ϕ}) \hat{ρ} + (F_{z} G_{ρ} - F_{ρ} G_{z}) \hat{ϕ} + (F_{ρ} G_{ϕ} - F_{ϕ} G_{ρ}) \hat{𝐳}$ .
$\nabla \cdot (J_{ρ} \hat{ρ} + J_{ϕ} \hat{ϕ} + J_{z} \hat{𝐳}) = \frac{1}{ρ} (\frac{\partial}{\partial ρ} (ρ J_{ρ}) + \frac{\partial F_{ϕ}}{\partial ϕ} + \frac{\partial}{\partial z} (ρ F_{z}))$ .
$\nabla \times (F_{ρ} \hat{ρ} + F_{ϕ} \hat{ϕ} + F_{z} \hat{𝐳}) = \frac{1}{ρ} (\frac{\partial F_{z}}{\partial ϕ} - \frac{\partial}{\partial z} (ρ F_{ϕ})) \hat{ρ} + (\frac{\partial F_{ρ}}{\partial z} - \frac{\partial F_{z}}{\partial ρ}) \hat{ϕ} + \frac{1}{ρ} (\frac{\partial}{\partial ρ} (ρ F_{ϕ}) - \frac{\partial F_{ρ}}{\partial ϕ}) \hat{𝐳}$ .
$\nabla U = \frac{\partial U}{\partial ρ} \hat{ρ} + \frac{1}{ρ} \frac{\partial U}{\partial ϕ} \hat{ϕ} + \frac{\partial U}{\partial z} \hat{𝐳}$ .

For spherical coordinates, $c_{1} = r$ , $c_{2} = θ$ , $c_{3} = ϕ$ , and ${\hat{𝐚}}^{1} = {\hat{𝐚}}_{1} = \hat{𝐫}$ , ${\hat{𝐚}}^{2} = {\hat{𝐚}}_{2} = \hat{θ}$ , ${\hat{𝐚}}^{3} = {\hat{𝐚}}_{3} = \hat{ϕ}$ , and $l_{1} = 1$ , $l_{2} = r$ , $l_{3} = r \sin θ$ . Therefore:

$(J_{r} \hat{𝐫} + J_{θ} \hat{θ} + J_{ϕ} \hat{ϕ}) \cdot (F_{r} \hat{𝐫} + F_{θ} \hat{θ} + F_{ϕ} \hat{ϕ}) = J_{r} F_{r} + J_{θ} F_{θ} + J_{ϕ} F_{ϕ}$ .
$(F_{r} \hat{𝐫} + F_{θ} \hat{θ} + F_{ϕ} \hat{ϕ}) \times (G_{r} \hat{𝐫} + G_{θ} \hat{θ} + G_{ϕ} \hat{ϕ}) = (F_{θ} G_{ϕ} - F_{ϕ} G_{θ}) \hat{𝐫} + (F_{ϕ} G_{r} - F_{r} G_{ϕ}) \hat{θ} + (F_{r} G_{θ} - F_{θ} G_{r}) \hat{ϕ}$ .
$\nabla \cdot (J_{r} \hat{𝐫} + J_{θ} \hat{θ} + J_{ϕ} \hat{ϕ}) = \frac{1}{r^{2} \sin θ} (\frac{\partial}{\partial r} (r^{2} \sin θ J_{r}) + \frac{\partial}{\partial θ} (r \sin θ J_{θ}) + \frac{\partial}{\partial ϕ} (r F_{ϕ}))$ .
$\nabla \times (F_{r} \hat{𝐫} + F_{θ} \hat{θ} + F_{ϕ} \hat{ϕ}) = \frac{1}{r^{2} \sin θ} (\frac{\partial}{\partial θ} (r \sin θ F_{ϕ}) - \frac{\partial}{\partial ϕ} (r F_{θ})) \hat{𝐫} + \frac{1}{r \sin θ} (\frac{\partial F_{r}}{\partial ϕ} - \frac{\partial}{\partial r} (r \sin θ F_{ϕ})) \hat{θ} + \frac{1}{r} (\frac{\partial}{\partial r} (r F_{θ}) - \frac{\partial F_{r}}{\partial θ}) \hat{ϕ}$ .
$\nabla U = \frac{\partial U}{\partial r} \hat{𝐫} + \frac{1}{r} \frac{\partial U}{\partial θ} \hat{θ} + \frac{1}{r \sin θ} \frac{\partial U}{\partial ϕ} \hat{ϕ}$ .

Intersection Boundaries

The endpoints of intersections

Many identities related to vector calculus can be derived from examining the endpoints of path-volume intersections and surface-surface intersections.

The endpoints of path-volume intersections

Start with a multi-path $𝐂$ , denoted by vector field $𝐉$ , and a multi-volume $𝐔$ , denoted by scalar field $U$ . The intersection $𝐂 \cap 𝐔$ is denoted by vector field $𝐉 U$ .

Any time a path $C$ with weight $w_{1}$ starts in a volume $Ω$ with weight $w_{2}$ , the intersection $𝐂 \cap 𝐔$ gains an endpoint at the starting point of $C$ with weight $w_{1} w_{2}$ . Any time a path $C$ with weight $w_{1}$ finishes in a volume $Ω$ with weight $w_{2}$ , the intersection $𝐂 \cap 𝐔$ gains an endpoint at the finishing point of $C$ with weight $- w_{1} w_{2}$ . The endpoints for $𝐂 \cap 𝐔$ that are generated when paths from $𝐂$ start or finish in volumes from $𝐔$ is the intersection of the endpoints of $𝐂$ with multi-volume $𝐔$ . This contributes the term $(\nabla \cdot 𝐉) U$ to $\nabla \cdot (𝐉 U)$ .

Any time a path $C$ with weight $w_{1}$ enters a volume $Ω$ with weight $w_{2}$ , the intersection $𝐂 \cap 𝐔$ gains an endpoint at the point of entry with weight $w_{1} w_{2}$ . Any time a path $C$ with weight $w_{1}$ leaves a volume $Ω$ with weight $w_{2}$ , the intersection $𝐂 \cap 𝐔$ gains an endpoint at the point of exit with weight $- w_{1} w_{2}$ . The endpoints for $𝐂 \cap 𝐔$ that are generated when paths from $𝐂$ enter or exit volumes from $𝐔$ is the intersection of multi-path $𝐂$ with the inwards oriented multi-surface of $𝐔$ . This contributes the term $𝐉 \cdot (\nabla U)$ to $\nabla \cdot (𝐉 U)$ .

The total endpoints of $𝐂 \cap 𝐔$ are: $\nabla \cdot (𝐉 U) = (\nabla \cdot 𝐉) U + 𝐉 \cdot (\nabla U)$ . In essence, the endpoints of $𝐂 \cap 𝐔$ are the endpoints of $𝐂$ that are contained in $𝐔$ , plus the points at which paths from $𝐂$ enter or exit volumes from $𝐔$ . This is depicted in the images on the right.

From the identity $\nabla \cdot (𝐉 U) = (\nabla \cdot 𝐉) U + 𝐉 \cdot (\nabla U)$ , counting the total point weight gives: $\underset{𝐪 \in ℝ^{3}}{∭} \nabla \cdot (𝐉 (𝐪) U (𝐪)) d V = \underset{𝐪 \in ℝ^{3}}{∭} (\nabla \cdot 𝐉 (𝐪)) U (𝐪) d V + \underset{𝐪 \in ℝ^{3}}{∭} 𝐉 (𝐪) \cdot (\nabla U (𝐪)) d V$ . For the endpoints of a multi-path, every starting point must be paired with a finishing point so the total point weight of the endpoints of a multi-path is 0. $\underset{𝐪 \in ℝ^{3}}{∭} \nabla \cdot (𝐉 (𝐪) U (𝐪)) d V = 0$ so hence $\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \cdot 𝐉 (𝐪)) U (𝐪) d V = - \underset{𝐪 \in ℝ^{3}}{∭} 𝐉 (𝐪) \cdot (\nabla U (𝐪)) d V$ . The total intersection between the endpoints of multi-path $𝐂$ and multi-volume $𝐔$ is the negative of the total intersection between $𝐂$ and the inwards oriented surface of $𝐔$ .

If $𝐉$ denotes a simple path $C$ that starts at point $𝐪_{0}$ and ends at point $𝐪_{1}$ , then the above integral identity becomes:

$\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \cdot δ_{1} (𝐪; C)) U (𝐪) d V = - \underset{𝐪 \in ℝ^{3}}{∭} δ_{1} (𝐪; C) \cdot (\nabla U (𝐪)) d V$ $⟺ \underset{𝐪 \in ℝ^{3}}{∭} (δ_{0} (𝐪; 𝐪_{0}) - δ_{0} (𝐪; 𝐪_{1})) U (𝐪) d V = - \int_{𝐪 \in C} (\nabla U (𝐪)) \cdot d 𝐪$ $⟺ U (𝐪_{0}) - U (𝐪_{1}) = - \int_{𝐪 \in C} (\nabla U (𝐪)) \cdot d 𝐪$ $⟺ \int_{𝐪 \in C} (\nabla U (𝐪)) \cdot d 𝐪 = U (𝐪_{1}) - U (𝐪_{0})$ This is known as the gradient theorem.

If $U$ denotes a simple volume $Ω$ with a outwards oriented surface $σ$ , then the integral identity becomes:

$\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \cdot 𝐉 (𝐪)) δ_{3} (𝐪; Ω) d V = - \underset{𝐪 \in ℝ^{3}}{∭} 𝐉 (𝐪) \cdot (\nabla δ_{3} (𝐪; Ω)) d V$ $⟺ \underset{𝐪 \in Ω}{∭} (\nabla \cdot 𝐉 (𝐪)) d V = - \underset{𝐪 \in ℝ^{3}}{∭} 𝐉 (𝐪) \cdot (- δ_{2} (𝐪; σ)) d V$ $⟺ \underset{𝐪 \in Ω}{∭} (\nabla \cdot 𝐉 (𝐪)) d V = \iint_{𝐪 \in σ} 𝐉 (𝐪) \cdot d 𝐒$ This is known as Gauss's divergence theorem.

In summary:

Given a multi-path denoted by vector field $𝐉$ , and a multi-volume denoted by scalar field $U$ , then the endpoints of the intersection are: $\nabla \cdot (𝐉 U) = (\nabla \cdot 𝐉) U + 𝐉 \cdot (\nabla U)$ .
Given a multi-path denoted by vector field $𝐉$ , and a multi-volume denoted by scalar field $U$ , then $\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \cdot 𝐉 (𝐪)) U (𝐪) d V = - \underset{𝐪 \in ℝ^{3}}{∭} 𝐉 (𝐪) \cdot (\nabla U (𝐪)) d V$ .
Given a simple path $C$ that starts at point $𝐪_{0}$ and ends at point $𝐪_{1}$ , and a multi-volume denoted by scalar field $U$ , then $\int_{𝐪 \in C} (\nabla U (𝐪)) \cdot d 𝐪 = U (𝐪_{1}) - U (𝐪_{0})$ . This is the gradient theorem.
Given a multi-path denoted by vector field $𝐉$ , and a simple volume $Ω$ with outwards oriented surface $σ$ , then $\underset{𝐪 \in Ω}{∭} (\nabla \cdot 𝐉 (𝐪)) d V = \iint_{𝐪 \in σ} 𝐉 (𝐪) \cdot d 𝐒$ . This is Gauss's divergence theorem.

The endpoints of surface-surface intersections

Start with multi-surface $𝐒_{1}$ , denoted by vector field $𝐅_{1}$ , and a second multi-surface $𝐒_{2}$ , denoted by vector field $𝐅_{2}$ . The intersection $𝐒_{1} \cap 𝐒_{2}$ is denoted by vector field $𝐅_{1} \times 𝐅_{2}$ .

Consider a surface $σ_{1}$ with weight $w_{1}$ from $𝐒_{1}$ , and a surface $σ_{2}$ with weight $w_{2}$ from $𝐒_{2}$ . Let $\partial σ_{1}$ denote the counter-clockwise boundary of $σ_{1}$ , and let $\partial σ_{2}$ denote the counter-clockwise boundary of $σ_{2}$ . There are 4 scenarios regarding the endpoints of $σ_{1} \cap σ_{2}$ :

When the $\partial σ_{1}$ intersects $σ_{2}$ in the preferred direction, the intersection point $\partial σ_{1} \cap σ_{2}$ has a weight of $+ w_{1} w_{2}$ , and an endpoint with weight $+ w_{1} w_{2}$ (starting point) for $σ_{1} \cap σ_{2}$ forms at $\partial σ_{1} \cap σ_{2}$ .
When the $\partial σ_{1}$ intersects $σ_{2}$ in the opposite direction, the intersection point $\partial σ_{1} \cap σ_{2}$ has a weight of $- w_{1} w_{2}$ , and an endpoint with weight $- w_{1} w_{2}$ (finishing point) for $σ_{1} \cap σ_{2}$ forms at $\partial σ_{1} \cap σ_{2}$ .
When the $\partial σ_{2}$ intersects $σ_{1}$ in the preferred direction, the intersection point $σ_{1} \cap \partial σ_{2}$ has a weight of $+ w_{1} w_{2}$ , and an endpoint with weight $- w_{1} w_{2}$ (finishing point) for $σ_{1} \cap σ_{2}$ forms at $σ_{1} \cap \partial σ_{2}$ .
When the $\partial σ_{2}$ intersects $σ_{1}$ in the opposite direction, the intersection point $σ_{1} \cap \partial σ_{2}$ has a weight of $- w_{1} w_{2}$ , and an endpoint with weight $+ w_{1} w_{2}$ (starting point) for $σ_{1} \cap σ_{2}$ forms at $σ_{1} \cap \partial σ_{2}$ .

It can be seen that the intersection $\partial σ_{1} \cap σ_{2}$ forms endpoints for $σ_{1} \cap σ_{2}$ with the correct polarity, and that the intersection $σ_{1} \cap \partial σ_{2}$ forms endpoints for $σ_{1} \cap σ_{2}$ with the opposite polarity. This can be observed in the image of the right. This implies that the endpoints of $𝐒_{1} \cap 𝐒_{2}$ are: $\nabla \cdot (𝐅_{1} \times 𝐅_{2}) = (\nabla \times 𝐅_{1}) \cdot 𝐅_{2} - 𝐅_{1} \cdot (\nabla \times 𝐅_{2})$ .

From the identity $\nabla \cdot (𝐅_{1} \times 𝐅_{2}) = (\nabla \times 𝐅_{1}) \cdot 𝐅_{2} - 𝐅_{1} \cdot (\nabla \times 𝐅_{2})$ , counting the total point weight gives: $\underset{𝐪 \in ℝ^{3}}{∭} \nabla \cdot (𝐅_{1} (𝐪) \times 𝐅_{2} (𝐪)) d V = \underset{𝐪 \in ℝ^{3}}{∭} (\nabla \times 𝐅_{1} (𝐪)) \cdot 𝐅_{2} (𝐪) d V - \underset{𝐪 \in ℝ^{3}}{∭} 𝐅_{1} (𝐪) \cdot (\nabla \times 𝐅_{2} (𝐪)) d V$ . For the endpoints of a multi-path, every starting point must be paired with a finishing point so the total point weight of the endpoints of a multi-path is 0. $\underset{𝐪 \in ℝ^{3}}{∭} \nabla \cdot (𝐅_{1} (𝐪) \times 𝐅_{2} (𝐪)) d V = 0$ so hence $\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \times 𝐅_{1} (𝐪)) \cdot 𝐅_{2} (𝐪) d V = \underset{𝐪 \in ℝ^{3}}{∭} 𝐅_{1} (𝐪) \cdot (\nabla \times 𝐅_{2} (𝐪)) d V$ . The total intersection of the counter-clockwise boundary of multi-surface $𝐒_{1}$ with multi-surface $𝐒_{2}$ is the total intersection of the counter-clockwise boundary of $𝐒_{2}$ with $𝐒_{1}$ . This is illustrated by the image on the right.

If $𝐅_{2}$ denotes a simple surface $σ$ with a counter-clockwise boundary $\partial σ$ , then the above integral identity becomes:

$\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \times 𝐅_{1} (𝐪)) \cdot δ_{2} (𝐪; σ) d V = \underset{𝐪 \in ℝ^{3}}{∭} 𝐅_{1} (𝐪) \cdot (\nabla \times δ_{2} (𝐪; σ)) d V$ $⟺ \iint_{𝐪 \in σ} (\nabla \times 𝐅_{1} (𝐪)) \cdot d 𝐒 = \underset{𝐪 \in ℝ^{3}}{∭} 𝐅_{1} (𝐪) \cdot δ_{1} (𝐪; \partial σ) d V$ $⟺ \iint_{𝐪 \in σ} (\nabla \times 𝐅_{1} (𝐪)) \cdot d 𝐒 = \int_{𝐪 \in \partial σ} 𝐅_{1} (𝐪) \cdot d 𝐪$

This is known as Stokes' theorem.

In summary:

Given two multi-surfaces denoted by vector fields $𝐅_{1}$ and $𝐅_{2}$ , then the endpoints of the intersection are: $\nabla \cdot (𝐅_{1} \times 𝐅_{2}) = (\nabla \times 𝐅_{1}) \cdot 𝐅_{2} - 𝐅_{1} \cdot (\nabla \times 𝐅_{2})$ .
Given two multi-surfaces denoted by vector fields $𝐅_{1}$ and $𝐅_{2}$ , then $\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \times 𝐅_{1} (𝐪)) \cdot 𝐅_{2} (𝐪) d V = \underset{𝐪 \in ℝ^{3}}{∭} 𝐅_{1} (𝐪) \cdot (\nabla \times 𝐅_{2} (𝐪)) d V$ .
Given a multi-surface denoted by vector field $𝐅_{1}$ and a simple surface $σ$ with counter-clockwise oriented boundary $\partial σ$ , then $\iint_{𝐪 \in σ} (\nabla \times 𝐅_{1} (𝐪)) \cdot d 𝐒 = \int_{𝐪 \in \partial σ} 𝐅_{1} (𝐪) \cdot d 𝐪$ . This is Stokes' theorem.

The boundaries of intersections

In addition to the identities derived from examining the endpoints of intersections, some more identities can be derived by examining the counter-clockwise boundaries of surfaces that result from intersections.

The counter-clockwise boundary of surface-volume intersections

Start with a multi-surface $𝐒$ , denoted by vector field $𝐅$ , and a multi-volume $𝐔$ , denoted by scalar field $U$ . The intersection $𝐒 \cap U$ is denoted by vector field $𝐅 U$ .

Consider a surface $σ$ with weight $w_{1}$ from $𝐒$ , and a volume $Ω$ with weight $w_{2}$ from $𝐔$ . Let $\partial σ$ denote the counter-clockwise boundary of $σ$ , and let $\partial Ω$ denote the inwards oriented surface of $Ω$ . There are two sources for the counter-clockwise boundary of $σ \cap Ω$ . Any time $\partial σ$ intersects $Ω$ , the intersection $\partial σ \cap Ω$ contributes to the boundary of $σ \cap Ω$ . When $\partial σ$ leaves $Ω$ , the boundary of $σ \cap Ω$ cannot follow, and instead must trace along the surface of $Ω$ while remaining in the surface $σ$ as indicated in the image to the right. The boundary of the total intersection $𝐒 \cap 𝐔$ , denoted by $\nabla \times (𝐅 U)$ , consists of two parts: the intersection of the boundary of $𝐒$ with $𝐔$ , denoted by $(\nabla \times 𝐅) U$ , and the intersection of the inwards-oriented surface of $𝐔$ with $𝐒$ , denoted by $(\nabla U) \times 𝐅 = - 𝐅 \times (\nabla U)$ . Therefore: $\nabla \times (𝐅 U) = (\nabla \times 𝐅) U + (\nabla U) \times 𝐅 = (\nabla \times 𝐅) U - 𝐅 \times (\nabla U)$ .

From the identity $\nabla \times (𝐅 U) = (\nabla \times 𝐅) U - 𝐅 \times (\nabla U)$ , computing the total displacement gives: $\underset{𝐪 \in ℝ^{3}}{∭} \nabla \times (𝐅 (𝐪) U (𝐪)) d V = \underset{𝐪 \in ℝ^{3}}{∭} (\nabla \times 𝐅 (𝐪)) U (𝐪) d V - \underset{𝐪 \in ℝ^{3}}{∭} 𝐅 (𝐪) \times (\nabla U (𝐪)) d V$ . The counter-clockwise boundary of a multi-surface is a closed multi-loop, and the total displacement generated by a loop is $𝟎$ . $\underset{𝐪 \in ℝ^{3}}{∭} \nabla \times (𝐅 (𝐪) U (𝐪)) d V = 𝟎$ so hence $\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \times 𝐅 (𝐪)) U (𝐪) d V = \underset{𝐪 \in ℝ^{3}}{∭} 𝐅 (𝐪) \times (\nabla U (𝐪)) d V$ . The total intersection of a the boundary of multi-surface $𝐒$ with multi-volume $𝐔$ is the total intersection of $𝐒$ with the surface of $𝐔$ .

If $𝐅$ denotes a simple surface $σ$ with counter-clockwise boundary $\partial σ$ , then the above integral identity becomes:

$\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \times δ_{2} (𝐪; σ)) U (𝐪) d V = \underset{𝐪 \in ℝ^{3}}{∭} δ_{2} (𝐪; σ) \times (\nabla U (𝐪)) d V$ $⟺ \underset{𝐪 \in ℝ^{3}}{∭} δ_{1} (𝐪; \partial σ) U (𝐪) d V = \iint_{𝐪 \in σ} d 𝐒 \times (\nabla U (𝐪))$ $⟺ \int_{𝐪 \in \partial σ} U (𝐪) d 𝐪 = - \iint_{𝐪 \in σ} (\nabla U (𝐪)) \times d 𝐒$

If $Ω$ denotes a simple volume with outwards oriented surface $σ$ , then the integral identity becomes:

$\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \times 𝐅 (𝐪)) δ_{3} (𝐪; Ω) d V = \underset{𝐪 \in ℝ^{3}}{∭} 𝐅 (𝐪) \times (\nabla δ_{3} (𝐪; Ω)) d V$ $⟺ \underset{𝐪 \in Ω}{∭} (\nabla \times 𝐅 (𝐪)) d V = \underset{𝐪 \in ℝ^{3}}{∭} 𝐅 (𝐪) \times (- δ_{2} (𝐪; σ)) d V$ $⟺ \underset{𝐪 \in Ω}{∭} (\nabla \times 𝐅 (𝐪)) d V = - \iint_{𝐪 \in σ} 𝐅 (𝐪) \times d 𝐒$

In summary:

Given a multi-surface denoted by vector field $𝐅$ , and a multi-volume denoted by scalar field $U$ , then the counter-clockwise boundary of the intersection is: $\nabla \times (𝐅 U) = (\nabla \times 𝐅) U - 𝐅 \times (\nabla U)$ .
Given a multi-surface denoted by vector field $𝐅$ , and a multi-volume denoted by scalar field $U$ , then $\underset{𝐪 \in ℝ^{3}}{∭} (\nabla \times 𝐅 (𝐪)) U (𝐪) d V = \underset{𝐪 \in ℝ^{3}}{∭} 𝐅 (𝐪) \times (\nabla U (𝐪)) d V$
Given a simple surface $σ$ with counter-clockwise boundary $\partial σ$ , and a multi-volume denoted by scalar field $U$ , then $\int_{𝐪 \in \partial σ} U (𝐪) d 𝐪 = - \iint_{𝐪 \in σ} (\nabla U (𝐪)) \times d 𝐒$ .
Given a multi-surface denoted by vector field $𝐅$ , and a simple volume $Ω$ with outwards oriented surface $σ$ , then $\underset{𝐪 \in Ω}{∭} (\nabla \times 𝐅 (𝐪)) d V = - \iint_{𝐪 \in σ} 𝐅 (𝐪) \times d 𝐒$ .

The surfaces of intersections

Some more identities can be derived by examining the surfaces of volumes that result from intersections.

The inwards-oriented surface of volume-volume intersections

Start with a multi-volume $𝐔_{1}$ , denoted by scalar field $U_{1}$ , and a second multi-volume $𝐔_{2}$ , denoted by scalar field $U_{2}$ . The intersection $𝐔_{1} \cap 𝐔_{2}$ is denoted by the scalar field $U_{1} U_{2}$ .

Consider a volume $Ω_{1}$ with weight $w_{1}$ from $𝐔_{1}$ , and a volume $Ω_{2}$ with weight $w_{2}$ from $𝐔_{2}$ . Let $σ_{1}$ denote the inwards-oriented surface of $Ω_{1}$ , and let $σ_{2}$ denote the inwards-oriented surface of $Ω_{2}$ . There are two parts to the inwards-oriented surface of the intersection $Ω_{1} \cap Ω_{2}$ , as shown in the image to the right. Part of the surface of $Ω_{1} \cap Ω_{2}$ consists of the portion of $σ_{2}$ that is contained by $Ω_{1}$ , which contributes the term $U_{1} (\nabla U_{2})$ to $\nabla (U_{1} U_{2})$ . The other part of the surface of $Ω_{1} \cap Ω_{2}$ consists of the portion of $σ_{1}$ that is contained by $Ω_{2}$ , which contributes the term $(\nabla U_{1}) U_{2}$ to $\nabla (U_{1} U_{2})$ . Therefore the total surface of $𝐔_{1} \cap 𝐔_{2}$ is $\nabla (U_{1} U_{2}) = U_{1} (\nabla U_{2}) + (\nabla U_{1}) U_{2}$ .

From the identity $\nabla (U_{1} U_{2}) = U_{1} (\nabla U_{2}) + (\nabla U_{1}) U_{2}$ , computing the total surface vector gives: $\underset{𝐪 \in ℝ^{3}}{∭} \nabla (U_{1} (𝐪) U_{2} (𝐪)) d V = \underset{𝐪 \in ℝ^{3}}{∭} U_{1} (𝐪) (\nabla U_{2} (𝐪)) d V + \underset{𝐪 \in ℝ^{3}}{∭} (\nabla U_{1} (𝐪)) U_{2} (𝐪) d V$ . The inwards-oriented surface of a multi-volume is a closed multi-surface, and the total surface vector of a closed surface is $𝟎$ . $\underset{𝐪 \in ℝ^{3}}{∭} \nabla (U_{1} (𝐪) U_{2} (𝐪)) d V = 𝟎$ so hence $\underset{𝐪 \in ℝ^{3}}{∭} U_{1} (𝐪) (\nabla U_{2} (𝐪)) d V = - \underset{𝐪 \in ℝ^{3}}{∭} (\nabla U_{1} (𝐪)) U_{2} (𝐪) d V$ . The total surface vector of the intersection of multi-volume $𝐔_{1}$ with the inwards oriented surface of multi-volume $𝐔_{2}$ is the opposite of the total surface vector of the intersection of the inwards oriented surface of $𝐔_{1}$ with $𝐔_{2}$ .

If $U_{1}$ denotes a simple volume $Ω$ with outwards oriented surface $σ$ , then the above integral identity becomes: $\underset{𝐪 \in ℝ^{3}}{∭} δ_{3} (𝐪; Ω) (\nabla U_{2} (𝐪)) d V = - \underset{𝐪 \in ℝ^{3}}{∭} (\nabla δ_{3} (𝐪; Ω)) U_{2} (𝐪) d V$ $⟺ \underset{𝐪 \in Ω}{∭} (\nabla U_{2} (𝐪)) d V = - \underset{𝐪 \in ℝ^{3}}{∭} (- δ_{2} (𝐪; σ)) U_{2} (𝐪) d V$ $⟺ \underset{𝐪 \in Ω}{∭} (\nabla U_{2} (𝐪)) d V = \iint_{𝐪 \in σ} U_{2} (𝐪) d 𝐒$

In summary:

Given two multi-volumes denoted by scalar fields $U_{1}$ and $U_{2}$ , then the inwards-oriented surface of the intersection is: $\nabla (U_{1} U_{2}) = U_{1} (\nabla U_{2}) + (\nabla U_{1}) U_{2}$ .
Given two multi-volumes denoted by scalar fields $U_{1}$ and $U_{2}$ , then $\underset{𝐪 \in ℝ^{3}}{∭} U_{1} (𝐪) (\nabla U_{2} (𝐪)) d V = - \underset{𝐪 \in ℝ^{3}}{∭} (\nabla U_{1} (𝐪)) U_{2} (𝐪) d V$ .
Given a simple volume $Ω$ with outwards oriented surface $σ$ and a multi-volume denoted by scalar field $U_{2}$ , then $\underset{𝐪 \in Ω}{∭} (\nabla U_{2} (𝐪)) d V = \iint_{𝐪 \in σ} U_{2} (𝐪) d 𝐒$ .

Summary

The tables below summarizes the results of the previous sections:

The endpoints, boundaries, and surfaces of intersections
structure 1	structure 2	intersection	endpoints, boundary, or surface
multi-path $𝐉_{1}$	multi-volume $U_{2}$	multi-path $𝐉_{1} U_{2}$	multi-point $\nabla \cdot (𝐉_{1} U_{2}) = (\nabla \cdot 𝐉_{1}) U_{2} + 𝐉_{1} \cdot (\nabla U_{2})$
multi-surface $𝐅_{1}$	multi-surface $𝐅_{2}$	multi-path $𝐅_{1} \times 𝐅_{2}$	multi-point $\nabla \cdot (𝐅_{1} \times 𝐅_{2}) = (\nabla \times 𝐅_{1}) \cdot 𝐅_{2} - 𝐅_{1} \cdot (\nabla \times 𝐅_{2})$
multi-surface $𝐅_{1}$	multi-volume $U_{2}$	multi-surface $𝐅_{1} U_{2}$	multi-path $\nabla \times (𝐅_{1} U_{2}) = (\nabla \times 𝐅_{1}) U_{2} - 𝐅_{1} \times (\nabla U_{2})$
multi-volume $U_{1}$	multi-volume $U_{2}$	multi-volume $U_{1} U_{2}$	multi-surface $\nabla (U_{1} U_{2}) = (\nabla U_{1}) U_{2} + U_{1} (\nabla U_{2})$

Integral identities
Simple structure	Multi-structure	Integral identity	Identity name
simple path $C$ , with starting point $C (0)$ and final point $C (1)$	multi-volume $U$	$\int_{𝐪 \in C} (\nabla U) \cdot d 𝐪 = U (C (1)) - U (C (0))$	the gradient theorem
simple volume $Ω$ , with outwards-oriented surface $σ$	multi-path $𝐉$	$\underset{𝐪 \in Ω}{∭} (\nabla \cdot 𝐉) d V = \iint_{𝐪 \in σ} 𝐉 \cdot d 𝐒$	Gauss's divergence theorem
simple surface $σ$ with counter-clockwise oriented boundary $\partial σ$	multi-surface $𝐅$	$\iint_{𝐪 \in σ} (\nabla \times 𝐅) \cdot d 𝐒 = \int_{𝐪 \in \partial σ} 𝐅 \cdot d 𝐪$	Stokes' theorem
simple surface $σ$ with counter-clockwise oriented boundary $\partial σ$	multi-volume $U$	$\iint_{𝐪 \in σ} (\nabla U) \times d 𝐒 = - \int_{𝐪 \in \partial σ} U d 𝐪$	unnamed
simple volume $Ω$ , with outwards-oriented surface $σ$	multi-surface $𝐅$	$\underset{𝐪 \in Ω}{∭} (\nabla \times 𝐅) d V = - \iint_{𝐪 \in σ} 𝐅 \times d 𝐒$	unnamed
simple volume $Ω$ , with outwards-oriented surface $σ$	multi-volume $U$	$\underset{𝐪 \in Ω}{∭} (\nabla U) d V = \iint_{𝐪 \in σ} U d 𝐒$	unnamed

Multi-path and Multi-surface duality

Template:BookCat

Calculus/Points, paths, surfaces, and volumes

Basic structures

Multi-points

Multi-paths

Multi-surfaces

Multi-volumes

At infinity

Totals

Total point weight

Total displacement

Total surface vector

Generalizing from surfaces in 2D space

Generalizing from displacement vectors

Summary

Total volume

Intersections

Point-Volume intersections

Path-Surface intersections

Path-Volume intersections

Surface-Surface intersections

Surface-Volume intersections

Volume-Volume intersections

Other intersections

Boundaries

The endpoints of paths

The counter-clockwise oriented boundaries of surfaces

The inwards-oriented surfaces of volumes

Closed loops and closed surfaces

Coordinate Systems

Computing various intersections

Computing path-surface intersections

Computing surface-surface intersections

Computing the endpoints of paths

Computing the counter-clockwise boundaries of surfaces

Computing the inwards-oriented surfaces of volumes

Summary

Orthogonal coordinate systems

Intersection Boundaries

The endpoints of intersections

The endpoints of path-volume intersections

The endpoints of surface-surface intersections

The boundaries of intersections

The counter-clockwise boundary of surface-volume intersections

The surfaces of intersections

The inwards-oriented surface of volume-volume intersections

Summary

Multi-path and Multi-surface duality

Navigation menu

Search