Coordinate system changing
This notebook illustrates how to use the Atlas package to solve problems in elementary differential geometry. As an example, we discuss changing from Cartesian to bipolar coordinate systems.
What we do?
- We construct a Cartesian domain (C) with the standard flat metric and a bipolar domain (B).
- We map C into B and then calculate metric tensor field and connection induced on B by the mapping.
- We also obtain the Laplace operator, Killing vector fields, geodesic lines and tangent vector fields on the bipolar domain.
- We show by direct computation that the geometry on domains C and B is just the same namely flat geometry of the plane.
Examples:
| Domain[manifold] | manifold - string - a manifold name or a name of a manifold domain. |
| Metric[id→expr] | id - variable - metric identifier,
expr - expression - metric declaration. |
| Connection[id] | id - variable - connection identifier. |
Necessary functions.
Cartesian domain
First of all we have to describe the space we are working in. The Cartesian domain is 2-dimensional Euclidean (flat) space i.e. a plane with a Cartesian coordinate system. To define the space we declare
domain, forms, vectors, coframe, frame, flat metric and calculate
connection (it equals to zero of caurse).
Define the Euclidean space as a manifold:
| Out[56]= |  |
Declare 1-forms for the space coframe:
| Out[57]= |  |
Declare the vectors for the space frame:
| Out[58]= |  |
Declare the coframe on the space:
| Out[59]= |  |
Declare the frame on the space:
| Out[60]= |  |
Declare a flat metric on the space:
| Out[61]= |  |
Calculate the connection of the metric:
| Out[62]= |  |
| Curvature[id] | id - variable - curvature identifier |
Necessary functions.
Bipolar domain
| Out[63]= |  |
Declare 1-forms for the domain coframe:
| Out[64]= |  |
Declare the vectors for the domain frame:
| Out[65]= |  |
Declare the coframe on the domain:
| Out[66]= |  |
Declare the frame of the domain:
| Out[67]= |  |
Declare a mapping of the domain into
R2:
| Out[68]= |  |
| Out[69]= |  |
Now we can calculate the metric induced on the domain by the mapping:
| Out[70]= |  |
Calculate the connection:
| Out[71]= |  |
Out[72]//MatrixForm= |
| |  |
The connection is not equal to zero. Nevertheless the domain has flat metric. To verify this fact we compute curvature 2-forms:
| Out[73]= |  |
Out[74]//MatrixForm= |
| |  |
Functions[f1→f1[x1,x2,...,xn],f2→f2[y1,y2,...,ym],...,
fk→fk[z1,z2,...,zj]] | fk=fk[z1, z2, ..., zj] - equations where fk-function identifier and zj - variables. |
| div[expr] | expr - any vector expression. |
| grad[expr] | expr-any 0-form expression. |
Necessary functions.
Calculations in B domain
Declare abstract functions on B domain:
| Out[75]= |  |
Calculate the gradient of the function f:
| Out[76]= |  |
Divergence of the frame vector fields:
| Out[77]= |  |
| Out[78]= |  |
The divergences are not equal to zero. Thus the coordinate lines on the bipolar domain "have sources and drains" (see graph paper above). To calculate the Laplace operator one can use the grad and div operators:
| Out[79]= |  |
Verify that there are no Killing vector fields in frame vector fields:
| Out[80]= |  |
How do we calculate the Killing vector fields in a bipolar domain? To do this we can use condition
Lx(G)=0 for the vector field X. Construct the vector field X as a linear combination of the frame vector fields:
| Out[81]= |  |
Calculate the corresponding Lie derivative:
| Out[82]= |  |
We have to solve the corresponding PDE system (
Lx(G)=0) and find the functions f and h. But ...But ...
But it is obvious that the vector fields
j are Killing ones in the Cartesian domain. Thus these vectors are Killing ones on the B domain also. It means that we have to express fields
j by the frame vector fields
Wi on B.
First of all we find the expressions for corresponding 1-forms using restriction of the mapping
:
| Out[83]= |  |
| Out[84]= |  |
After that we use
Dual- operator to convert 1-forms into corresponding vector fields ("rising indexes"):
| Out[85]= |  |
| Out[86]= |  |
| Out[87]= |  |
| Out[88]= |  |
By the way there are "no sources and drains" for these vector fields (just straight lines; see graph paper above):
| Out[89]= |  |
| Out[90]= |  |
How to calculate the geodesic lines? To do this we can use the condition for tangent vectors X of a geodesic line.
Calculate the corresponding covariant derivative:
| Out[91]= |  |
We have to solve the corresponding PDE system (
cov(X,X)=0) and find functions f and h. But ...
But the obvious result is as follows:
| Out[92]= |  |
| Out[93]= |  |
Nothing surprising about it really. The vector fields E1, E2 are tangent fields of straight lines on the plane (just expressed by frame vectors of B domain). We are on the plane anyway (but with different coordinate system) thus the geometry on B is the same as on C.
- Lie derivative
- Exterior product
- Tensor product
- Hodge operator
- Covariant derivative
