## Exploring Stokes Theorem in tensor form.

Posted by peeterjoot on February 22, 2011

# 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.

See the post ‘stokes theorem in geometric algebra‘ [PDF], where this topic has been revisited with this in mind.

# Original Post:

[Click here for a PDF of this post with nicer formatting]

## Motivation.

I’ve worked through Stokes theorem concepts a couple times on my own now. One of the first times, I was trying to formulate this in a Geometric Algebra context. I had to resort to a tensor decomposition, and pictures, before ending back in the Geometric Algebra description. Later I figured out how to do it entirely with a Geometric Algebra description, and was able to eliminate reliance on the pictures that made the path to generalization to higher dimensional spaces unclear.

It’s my expectation that if one started with a tensor description, the proof entirely in tensor form would not be difficult. This is what I’d like to try this time. To start off, I’ll temporarily use the Geometric Algebra curl expression so I know what my tensor equation starting point will be, but once that starting point is found, we can work entirely in coordinate representation. For somebody who already knows that this is the starting point, all of this initial motivation can be skipped.

# Translating the exterior derivative to a coordinate representation.

Our starting point is a curl, dotted with a volume element of the same grade, so that the result is a scalar

Here is a blade of grade , and we wedge this with the gradient for the space

where we with with a basis (not necessarily orthonormal) , and the reciprocal frame for that basis defined by the relation

Our coordinates in these basis sets are

so that

The operator coordinates of the gradient are defined in the usual fashion

The volume element for the subspace that we are integrating over we will define in terms of an arbitrary parametrization

The subspace can be considered spanned by the differential elements in each of the respective curves where all but the th parameter are held constant.

We assume that the integral is being performed in a subspace for which none of these differential elements in that region are linearly dependent (i.e. our Jacobean determinant must be non-zero).

The magnitude of the wedge product of all such differential elements provides the volume of the parallelogram, or parallelepiped (or higher dimensional analogue), and is

The volume element is a oriented quantity, and may be adjusted with an arbitrary sign (or equivalently an arbitrary permutation of the differential elements in the wedge product), and we’ll see that it is convenient for the translation to tensor form, to express these in reversed order.

Let’s write

so that our volume element in coordinate form is

Our curl will also also be a grade blade. We write for the grade blade

where is antisymmetric (i.e. for a some set of vectors ).

With our gradient in coordinate form

the curl is then

The differential form for our integral can now be computed by expanding out the dot product. We want

Evaluation of the interior dot products introduces the intrinsic antisymmetry required for Stokes theorem. For example, with

Since our end result is a completely antisymmetric set of permutations of all the deltas

and the curl integral takes it’s coordinate form

One final contraction of the paired indexes gives us our Stokes integral in its coordinate representation

We now have a starting point that is free of any of the abstraction of Geometric Algebra or differential forms. We can identify the products of partials here as components of a scalar hypervolume element (possibly signed depending on the orientation of the parametrization)

This is also a specific computation recipe for these hypervolume components, something that may not be obvious when we allow for general metrics for the space. We are also allowing for non-orthonormal coordinate representations, and arbitrary parametrization of the subspace that we are integrating over (our integral need not have the same dimension as the underlying vector space).

Observe that when the number of parameters equals the dimension of the space, we can write out the antisymmetric term utilizing the determinant of the Jacobian matrix

When the dimension of the space is greater than the number of parameters for the integration hypervolume in question, the antisymmetric sum of partials is still the determinant of a Jacobian matrix

however, we will have one such Jacobian for each unique choice of indexes.

# The Stokes work starts here.

The task is to relate our integral to the boundary of this volume, coming up with an explicit recipe for the description of that bounding surface, and determining the exact form of the reduced rank integral. This job is essentially to reduce the ranks of the tensors that are being contracted in our Stokes integral. With the derivative applied to our rank antisymmetric tensor , we can apply the chain rule and examine the permutations so that this can be rewritten as a contraction of itself with a set of rank surface area elements.

Now, while the setup here has been completely general, this task is motivated by study of special relativity, where there is a requirement to work in a four dimensional space. Because of that explicit goal, I’m not going to attempt to formulate this in a completely abstract fashion. That task is really one of introducing sufficiently general notation. Instead, I’m going to proceed with a simpleton approach, and do this explicitly, and repeatedly for each of the rank 1, rank 2, and rank 3 tensor cases. It will be clear how this all generalizes by doing so, should one wish to work in still higher dimensional spaces.

## The rank 1 tensor case.

The equation we are working with for this vector case is

Expanding out the antisymmetric partials we have

with which we can reduce the integral to

Now, if it happens that

then each of the individual integrals in and can be carried out. In that case, without any real loss of generality we can designate the integration bounds over the unit parametrization space square , allowing this integral to be expressed as

It’s also fairly common to see used to designate evaluation of this first integral on the boundary, and using this we write

Also note that since we are summing over all , and have

we can write this summing over all unique pairs of instead, which eliminates a small bit of redundancy (especially once the dimension of the vector space gets higher)

In this form we have recovered the original geometric structure, with components of the curl multiplied by the component of the area element that shares the orientation and direction of that portion of the curl bivector.

This form of the result with evaluation at the boundaries in this form, assumed that was not a function of and was not a function of . When that is not the case, we appear to have a less pretty result

Can this be reduced any further in the general case? Having seen the statements of Stokes theorem in it’s differential forms formulation, I initially expected the answer was yes, and only when I got to evaluating my spacetime example below did I realize that the differentials displacements for the parallelogram that constituted the area element were functions of both parameters. Perhaps this detail is there in the differential forms version of the general Stokes theorem too, but is just hidden in a tricky fashion by the compact notation.

### Sanity check: case in rectangular coordinates.

For , and , we have for the LHS

Our RHS expands to

We have

The RHS is just a positively oriented line integral around the rectangle of integration

This special case is also recognizable as Green’s theorem, evident with the substitution , , which gives us

Strictly speaking, Green’s theorem is more general, since it applies to integration regions more general than rectangles, but that generalization can be arrived at easily enough, once the region is broken down into adjoining elementary regions.

### Sanity check: case in rectangular coordinates.

It is expected that we can recover the classical Kelvin-Stokes theorem if we use rectangular coordinates in . However, we see that we have to consider three different parametrizations. If one picks rectangular parametrizations in sequence, in each case holding the value of the additional coordinate fixed, we get three different independent Green’s function like relations

Note that we cannot just add these to form a complete integral since the curves are all have different orientations. To recover the Stokes theorem in rectangular coordinates, it appears that we’d have to consider a Riemann sum of triangular surface elements, and relate that to the loops over each of the surface elements. In that limiting argument, only the boundary of the complete surface would contribute to the RHS of the relation.

All that said, we shouldn’t actually have to go to all this work. Instead we can stick to a two variable parametrization of the surface, and use 3.30 directly.

### An illustration for a spacetime surface.

Suppose we have a particle trajectory defined by an active Lorentz transformation from an initial spacetime point

Let the Lorentz transformation be formed by a composition of boost and rotation

Different rates of evolution of and define different trajectories, and taken together we have a surface described by the two parameters

We can compute displacements along the trajectories formed by keeping either or fixed and varying the other. Those are

Writing the computation of the partials above yields

Different choices of the initial point yield different surfaces, but we can get the idea by picking a simple starting point leaving

We can now compute our Jacobian determinants

Those are

Using this, let’s see a specific 4D example in spacetime for the integral of the curl of some four vector , enumerating all the non-zero components of 3.31 for this particular spacetime surface

The LHS is thus found to be

On the RHS we have

Because of the complexity of the surface, only the second term on the RHS has the “evaluate on the boundary” characteristic that may have been expected from a Green’s theorem like line integral.

It is also worthwhile to point out that we have had to be very careful with upper and lower indexes all along (and have done so with the expectation that our application would include the special relativity case where our metric determinant is minus one.) Because we worked with upper indexes for the area element, we had to work with lower indexes for the four vector and the components of the gradient that we included in our curl evaluation.

## The rank 2 tensor case.

Let’s consider briefly the terms in the contraction sum

For any choice of a set of three distinct indexes ), we have ways of permuting those indexes in this sum

Observe that we have no sign alternation like we had in the vector (rank 1 tensor) case. That sign alternation in this summation expansion appears to occur only for odd grade tensors.

Returning to the problem, we wish to expand the determinant in order to apply a chain rule contraction as done in the rank-1 case. This can be done along any of rows or columns of the determinant, and we can write any of

This allows the contraction of the index , eliminating it from the result

Dividing out the common terms, we can summarize this result as

In general, as observed in the spacetime surface example above, the two index Jacobians can be functions of the integration variable first being eliminated. In the special cases where this is not the case (such as the case with rectangular coordinates), then we are left with just the evaluation of the tensor element on the boundaries of the respective integrals.

## The rank 3 tensor case.

The key step is once again just a determinant expansion

so that the sum can be reduced from a four index contraction to a 3 index contraction

That’s the essence of the theorem, but we can play the same combinatorial reduction games to reduce the built in redundancy in the result

## A note on Four diverence.

Our four divergence integral has the following form

We can relate this to the rank 3 Stokes theorem with a duality transformation, multiplying with a pseudoscalar

where can also be related back to the vector by the same sort of duality transformation

The divergence integral in terms of the rank 3 tensor is

and we are free to perform the same Stokes reduction of the integral. Of course, this is particularly simple in rectangular coordinates. I still have to think though one sublty that I feel may be important. We could have started off with an integral of the following form

and I think this differs from our starting point slightly because this has none of the antisymmetric structure of the signed 4 volume element that we have used. We do not take the absolute value of our Jacobians anywhere.

## Leave a Reply