Stokes theorem in Geometric Algebra formalism.
Posted by peeterjoot on July 22, 2009
Obsolete with potential errors.
This post may be in error. I wrote this before understanding that the gradient used in Stokes Theorem must be projected onto the tangent space of the parameterized surface, as detailed in Alan MacDonald’s Vector and Geometric Calculus.
Relying on pictorial means and a brute force ugly comparison of left and right hand sides, a verification of Stokes theorem for the vector and bivector cases was performed (). This was more of a confirmation than a derivation, and the technique fails the transition to the trivector case. The trivector case is of particular interest in electromagnetism since that and a duality transformation provides a four-vector divergence theorem.
The fact that the pictorial means of defining the boundary surface doesn’t work well in four vector space is not the only unsatisfactory aspect of the previous treatment. The fact that a coordinate expansion of the hypervolume element and hypersurface element was performed in the LHS and RHS comparisons was required is particularly ugly. It is a lot of work and essentially has to be undone on the opposing side of the equation. Comparing to previous attempts to come to terms with Stokes theorem in () and () this more recent attempt at least avoids the requirement for a tensor expansion of the vector or bivector. It should be possible to build on this and minimize the amount of coordinate expansion required and go directly from the volume integral to the expression of the boundary surface.
Notation and Setup.
The desire is to relate the curl hypervolume integral to a hypersurface integral on the boundary
In order to put meaning to this statement the volume and surface elements need to be properly defined. In order that this be a scalar equation, the object in the integral is required to be of grade , and where is the dimension of the vector space that generates the object .
As evident in equation (2.1) a metric is required to define the dot product. If an affine non-metric formulation
of Stokes theorem is possible it will not be attempted here. A reciprocal basis pair will be utilized, defined by
Both of the sets and are taken to span the space, but are not required to be orthogonal. The notation is consistent with the Dirac reciprocal basis, and there will not be anything in this treatment that prohibits the Minkowski metric signature required for such a relativistic space.
Vector decomposition in terms of coordinates follows by taking dot products. We write
When working with a non-orthonormal basis, use of the reciprocal frame can be utilized to express the gradient.
This contains what may perhaps seem like an odd seeming mix of upper and lower indexes in this definition. This is how the gradient is defined in . Although it is possible to accept this definition and work with it, this form can be justified by require of the gradient consistency with the the definition of directional derivative. A definition of the directional derivative that works for single and multivector functions, in and other more general spaces is
Taylor expanding about in terms of coordinates we have
The lower index representation of the vector coordinates could also have been used, so using the directional derivative to imply a definition of the gradient, we have an additional alternate representation of the gradient
We define the hypervolume in terms of parametrized vector displacements . For the vector x we can form a pseudoscalar for the subspace spanned by this parametrization by wedging the displacements in each of the directions defined by variation of the parameters. For let
so the hypervolume element for the subspace in question is
This can be expanded explicitly in coordinates
Observe that when is also the dimension of the space, we can employ a pseudoscalar and can specify our volume element in terms of the Jacobian determinant.
However, we won’t have a requirement to express the Stokes result in terms of such Jacobians.
Expansion of the curl and volume element product
We are now prepared to go on to the meat of the issue. The first order of business is the expansion of the curl and volume element product
The wedge product within the scalar grade selection operator can be expanded in symmetric or antisymmetric sums, but this is a grade dependent operation. For odd grade blades (vector, trivector, …), and vector we have for the dot and wedge product respectively
Similarly for even grade blades we have
First treating the odd grade case for we have
Employing cyclic scalar reordering within the scalar product for the first term
The end result is
For even grade (and thus odd grade ) it is straightforward to show that (2.11) also holds.
Expanding the volume dot product
We want to expand the volume integral dot product
Picking will serve to illustrate the pattern, and the generalization (or degeneralization to lower grades) will be clear. We have
This avoids the requirement to do the entire Jacobian expansion of (2.9). The dot product of the differential displacement with can now be made explicit without as much mess.
We now have products of the form
Now we see that the differential form of (2.11) for this example is reduced to
While 2.11 was a statement of Stokes theorem in this Geometric Algebra formulation, it was really incomplete without this explicit expansion of . This expansion for the case serves to illustrate that we would write Stokes theorem as
Here the indexes have the range . This with the definitions 2.7, and 2.8 is really Stokes theorem in its full glory.
Observe that in this Geometric algebra form, the one forms are nothing more abstract that plain old vector differential elements. In the formalism of differential forms, this would be vectors, and would be a form. In a context where we are working with vectors, or blades already, the Geometric Algebra statement of the theorem avoids a requirement to translate to the language of forms.
With a statement of the general theorem complete, let’s return to our case where we can now integrate over each of the parameters. That is
This is precisely Stokes theorem for the trivector case and makes the enumeration of the boundary surfaces explicit. As derived there was no requirement for an orthonormal basis, nor a Euclidean metric, nor a parametrization along the basis directions. The only requirement of the parametrization is that the associated volume element is non-trivial (i.e. none of ).
For completeness, note that our boundary surface and associated Stokes statement for the bivector and vector cases is, by inspection respectively
These three expansions can be summarized by the original single statement of (2.1), which repeating for reference, is
Where it is implied that the blade is evaluated on the boundaries and dotted with the associated hypersurface boundary element. However, having expanded this we now have an explicit statement of exactly what that surface element is now for any desired parametrization.
Duality relations and special cases.
Some special (and more recognizable) cases of (2.1) are possible considering specific grades of , and in some cases employing duality relations.
curl surface integral
One important case is the vector result, which can be expressed in terms of the cross product.
Write . Then we have
This recovers the familiar cross product form of Stokes law.
3D divergence theorem
Duality applied to the bivector Stokes result provides the divergence theorem in . For bivector , let , , and . We then have
This recovers the divergence equation
4D divergence theorem
How about the four dimensional spacetime divergence? Write, express a trivector as a dual four-vector , and the four volume element . This gives
For the boundary volume integral write , for
So we have
the orientation of the fourspace volume element and the boundary normal is defined in terms of the parametrization, the duality relations and our explicit expansion of the 4D stokes boundary integral above.
4D divergence theorem, continued.
The basic idea of using duality to express the 4D divergence integral as a stokes boundary surface integral has been explored. Lets consider this in more detail picking a specific parametrization, namely rectangular four vector coordinates. For the volume element write
As seen previously (but not separately), the divergence can be expressed as the dual of the curl
So we have . Putting things together, and writing we have
It is straightforward to reduce each of these dot products. For example
The rest proceed the same and rather anticlimactically we end up coming full circle
This is however nothing more than the definition of the divergence itself and no need to resort to Stokes theorem is required. However, if we are integrating over a rectangle and perform each of the four integrals, we have (with ) from the dual Stokes equation the perhaps less obvious result
When stated this way one sees that this could have just as easily have followed directly from the left hand side. What’s the point then of the divergence theorem or Stokes theorem? I think that the value must really be the fact that the Stokes formulation naturally builds the volume element in a fashion independent of any specific parametrization. Here in rectangular coordinates the result seems obvious, but would the equivalent result seem obvious if non-rectangular spacetime coordinates were employed? Probably not.
 Peeter Joot. Stokes theorem applied to vector and bivector fields [online]. http://sites.google.com/site/peeterjoot/math2009/stokesGradeTwo.pdf.
 Peeter Joot. Stokes law in wedge product form [online]. http://sites.google.com/site/peeterjoot/geometric-algebra/vector_integral_relations.pdf.
 Peeter Joot. Stokes Law revisited with algebraic enumeration of boundary [online]. http://sites.google.com/site/peeterjoot/geometric-algebra/stokes_revisited.pdf.
 C. Doran and A.N. Lasenby. Geometric algebra for physicists. Cambridge University Press New York, Cambridge, UK, 1st edition, 2003.