2D Coordinate system changing - Cartesian to Elliptic
What we do?In this article we discuss coordinate system changing from Cartesian to Elliptic.
- We construct a Cartesian domain (C) with the standard flat metric and an Elliptic domain (E).
- We map C into E and then calculate metric tensor field and connection induced on E by the mapping.
- We also obtain the Laplace operator, Killing vector fields, geodesic lines and tangent vector fields on the elliptic domain.
- We show by straight computations that geometry on domains C and E is just the same namely flat geometry of the plane.
What we use?
Let's begin and have fun!
Load atlas package:
Redefine `atlas/simp` procedure to simplify the results:
`atlas/simp` := proc (a) factor(simplify(a, trig)) end proc;
Define the Euclidean space as a manifold:
Declare 1-forms for the space coframe:
Forms(e[k] = 1);
Declare the vectors for the space frame:
Declare the coframe on the space:
Coframe(e = d(x), e = d(y));
Declare the frame on the space:
Declare a flat metric on the space:
Metric(g = `&.`(d(x), d(x))+`&.`(d(y), d(y)));
Calculate the connection of the metric:
Elliptic domainThe Elliptic domain is a plane with 2-dimensional orthogonal coordinate system in which the coordinate lines are confocal ellipses and hyperbolae.
Define new domain:
Declare 1-forms for the domain coframe:
Forms(w[i] = 1);
Declare the vectors for the domain frame:
Declare the coframe on the domain:
Coframe(w = d(u), w = d(v));
Declare the frame of the domain:
Declare a mapping of the domain into :
Mapping(pi, E, C, x = cosh(u)*cos(v), y = sinh(u)*sin(v));
Now we can calculate the metric induced on the domain by the mapping:
Metric(G = `&/`(g, pi));
Calculate the connection:
Let us see the result:
Connection is not equal to zero. Nevertheless the domain has flat metric. To verify this fact we compute curvature 2-forms:
Calculations in E domainAnd now the most fun part!!! Let's answer the several questions:
Declare abstract functions on E domain:
Functions(f = f(u, v), h = h(u, v));
Calculate the gradient of the function:
Divergence of the frame vector fields:
'div(W)' = div(W);
'div(W)' = div(W);
The divergences are not equal zero.
To calculate the Laplace operator one can use the grad and div operators:
Verify that there are no Killing vector fields in frame vector fields:
'L[W[j]](G)' = L[W[j]](G);
To do this we can use condition for the vector field X.
Construct the vector field X as a linear combination of the frame vector fields:
X := f*W+h*W;
Calculate the corresponding Lie derivative:
But it is obvious that the vector fields are Killing ones in the Cartesian domain. Thus these vectors are Killing ones on the E domain also. It means that we have to express fields by frame vector fields on E.
First of all we find the expressions for corresponding 1-forms using restriction of the mapping:
e1 := `&/`(e, pi);
e2 := `&/`(e, pi);
After that we use dual operator to convert 1-forms into corresponding vector fields ("rising indexes"):
E1 := dual(e1);
E2 := dual(e2);
Verify the result:
'L['E1'](G)' = simplify(L[E1](G));
'L['E2'](G)' = simplify(L[E2](G));
By the way there are "no sources and drains" for these vector fields (just straight lines; see graph paper above):
'div(E1)' = div(E1);
'div(E2)' = div(E2);
To do this we can use the condition for tangent vectors X of a geodesic line.
Calculate the corresponding covariant derivative:
'cov[X](X)' = collect(cov[X](X), [W, W]);
We have to solve the corresponding PDE system () and find functions f and h. But ...
But the obvious result is as follows:
'cov['E1'](E1)' = simplify(cov[E1](E1));
'cov['E2'](E2)' = simplify(cov[E2](E2));
ConclusionsSo, what do we get?
We construct two domains (Cartesian and Elliptic) and map them one to another.
Calculate the Laplace operator, Killing vector and so on.
Also we show that geometry on the Cartesian and Elliptic domains is just the same namely flat geometry of the plane.
Although some formulas look quite big and impressive, but solving the problem take less than 10 minutes for me.