Direct Flow Parametric Proof of Gauss' Divergence Theorem

Gauss' Divergence Theorem stands as a cornerstone of vector calculus, bridging the gap between the flux of a vector field through a closed surface and the divergence of that field within the volume enclosed by the surface. While there are several approaches to proving this theorem, including proofs based on the definition of divergence and decompositions into cuboids, a direct flow parametric proof offers profound insight into the geometric nature of the relationship. This method utilizes the power of differential calculus and parameter transformations to equate the flow through a boundary to the behavior of the field inside.

Theorem Statement

Let \( V \) be a compact volume in three-dimensional space bounded by a piecewise smooth closed surface \( S \) with outward unit normal vector \( \mathbf{n} \). Let \( \mathbf{F} \) be a continuously differentiable vector field defined on an open region containing \( V \). The theorem states:

The left-hand side represents the total flux (flow) of the field outward through the surface, while the right-hand side represents the total divergence (source rate) inside the volume.

Parametric Representation of the Volume

To construct a direct parametric proof, we first establish a mapping between a simpler reference domain and our complex volume. Consider a transformation \( \mathbf{r} \) that maps a unit cube \( C \) in the parameter space \( (u, v, w) \) to the volume \( V \) in physical space \( (x, y, z) \).

We assume this mapping is smooth and one-to-one on the interior of the domain. The fundamental tool for linking the parameter space to the physical space is the Jacobian matrix \( J \) and the Jacobian determinant \( J_{det} \). The position vectors tangent to the coordinate curves in the parameter domain are given by:

The volume element \( dV \) transforms according to the absolute value of the scalar triple product of these tangent vectors:

The Jacobian Identity

The crux of the direct parametric proof lies in a specific vector calculus identity involving the transformed field. While the divergence \( \nabla \cdot \mathbf{F} \) is a scalar defined in \( x,y,z \) space, we can relate it to the parameters via the following key identity:

This identity effectively distributes the divergence operator across the geometric stretching factors of the transformation. It asserts that the weighted divergence in the volume equals the sum of partial derivatives of the flux-like terms through the faces of the infinitesimal parameter cube.

Integration over the Parameter Volume

We now evaluate the right-hand side of the theorem: the volume integral of the divergence. Substituting the Jacobian identity, we rewrite the integral in parameter coordinates \( (u, v, w) \) over the unit cube \( C \):

Substituting the expanded identity from the previous section:

Applying the Fundamental Theorem of Calculus

The power of the parametric approach becomes apparent when we apply the Fundamental Theorem of Calculus to the triple integral. We can separate the integral into three parts and integrate each with respect to its specific variable.

Consider the first term:

Integrating with respect to \( u \) from 0 to 1 yields:

This expression represents the evaluation of the vector field flux through the faces of the cube where \( u \) is constant (u=0 and u=1).

Mapping to the Surface Integral

We must now interpret these boundary terms physically. The vector \( \mathbf{t}_v \times \mathbf{t}_w \) corresponds to the vector area element \( d\mathbf{S}_u \) on the surface patches mapped from the parameter lines of constant \( u \).

Specifically, on the face where \( u=1 \), the outward normal corresponds to \( \mathbf{t}_v \times \mathbf{t}_w \). On the face where \( u=0 \), the outward normal opposes the tangent orientation, corresponding to \( -(\mathbf{t}_v \times \mathbf{t}_w) \).

Therefore, the difference derived in the previous step is exactly the contribution of the flux through the two surfaces patches \( S_u^+ \) and \( S_u^- \) corresponding to the parameter boundaries: