Physics[ThreePlusOne][D3_] - indexed covariant differential operator with respect to curvilinear space coordinates acting on a 3D hypersurface
|
Calling Sequence
|
|
D3_[mu](A)
D3_[mu](A, [X])
D3_(A)
D3_(A, [X])
|
|
Parameters
|
|
mu
|
-
|
a name representing an integer number between 1 and the spacetime dimension or a space index running from 1 to the dimension - 1, can also be the number itself or preceded by ~ to represent a contravariant index
|
A
|
-
|
any mathematical expression or relation between expressions, or a set or list of them, or an rtable
|
[X]
|
-
|
a list of differentiation variables to which the index j refers
|
|
|
|
|
Description
|
|
•
|
As is the case of all the tensors of the ThreePlusOne package, D3_[mu] is a 4D spacetime tensor with special meaning (acting) on a 3D hypersurface specified by the values of the Lapse and Shift. The D3_[mu] command is a computational representation for the covariant differential operator in that hypersurface. Its space components also form a tensor in the 3D hypersurface, so D3_ can be used with 3D space indices, for instances as in D_[i](A[j]).
|
•
|
On the screen, D3_, is displayed as the corresponding D_ command, but in black , instead of in blue .
|
•
|
Note that the space components of - say - D3_[mu](A[nu](x, y, z, t)) are not necessarily the same as the space components of the 4D tensor D_[mu](A[nu](x, y, z, t)). Instead, the relationship between these two is given by
|
|
where is the gamma3_, a 4D tensor, the spatial metric induced in the 3D hypersurface by the 4D g_ metric, also a projection operator that projects 4D tensors into the 3D hypersurface, resulting in purely spatial 4D tensors, all of whose components are equal to zero when their indices are contravariant and any of them has a timelike value (i.e. contravariant 0).
|
•
|
D3_ can be used as well without an index, as in D3_(A) displayed as , in which case it represents the total differential in curvilinear coordinates, and the output comes automatically expanded as D3_[mu](A) * D3_(X[~mu]), where X[~mu] is the corresponding spacetime vector.
|
•
|
The %D3_ command is the inert form of D3_, so it represents the same mathematical operation but without performing it. To perform the operation, use value.
|
|
|
Examples
|
|
>
|
|
| |
| |
| |
| |
| |
| |
| (1) |
>
|
|
| (2) |
Set up a coordinate system to work with - the first one to be set is automatically taken as the differentiation variables for D_ and so also for D3_
>
|
|
| |
When Physics is initialized, the default spacetime metric is of Minkowski type. You can see the metric querying Setup, as in Setup(metric);, or directly entering the metric g_ with no indices
| (4) |
Check the nonzero components of Christoffel: since this is a Minkowski spacetime, there are none, then
The nonzero components of Christoffel3 are also none
Note the difference in color: the 4D tensors are displayed in blue while the corresponding tensors of ThreePlusOne are displayed in black.
Because the Christoffel symbols are all equal to zero, the covariant derivative is equal to the standard derivative expressed using the d_ operator, both for D_ and D3_
| (7) |
Set the scenario as a Schwarzschild spacetime in spherical coordinates; you can do this entering Setup(metric = Schwarzschild) or in the simpler form taking advantage of abbreviations and directly using the spacetime metric g_ command
| |
| |
| |
| |
| |
| (8) |
The covariant derivative of a scalar is always equal to the derivative
The covariant differential of a scalar function
| (10) |
For illustration purposes Define an arbitrary 3D space tensor ; to indicate that this is a 3D tensor, define it with its space indices explicit, for example
| |
| (11) |
Use a compact display for this function to avoid redundant repeated display of its functionality plus having derivatives displayed indexed by the differentiation variables
| (12) |
Now we have, for the total 3D covariant differential and the 3D covariant derivative, respectively,
| (13) |
| (14) |
This D3_ covariant derivative is expressed in terms of the Christoffel3 symbols in the same way the covariant derivative D_ is expressed in terms of the 4D Christoffel. To see that you can use expand or convert to d_
| (16) |
The generalized divergence of a vector in curvilinear coordinates
| (18) |
| (19) |
To compute with a representation for D3_ without actually performing the operation, use the inert form %D3_. To afterwards perform the operation use value. For example, when the space part of the 3+1 metric gamma3_ is diagonal, its covariant derivative is equal to 0
The general form of such a 3+1 metric happens when the space part of the 4D metric g_ is diagonal. For example, using arbitrary functions of X for the other components of the 4D metric, we have
>
|
|
| |
| |
| |
| |
| |
| (22) |
Set the spacetime metric g_ by giving the square of the spacetime interval such that
>
|
|
| (24) |
>
|
|
| |
| |
| |
| |
| |
| |
| |
| (25) |
The divergence is now not zero
Compute for instance its first component, for
>
|
|
|
|
See Also
|
|
Christoffel, CompactDisplay, Coordinates, D_, Define, g_, Physics, Physics conventions, Physics examples, Physics Updates, Tensors - a complete guide, Mini-Course Computer Algebra for Physicists, Setup, ThreePlusOne:-gamma3_, value
|
|
References
|
|
|
[1] Landau, L.D., and Lifshitz, E.M. The Classical Theory of Fields, Course of Theoretical Physics Volume 2, fourth revised English edition. Elsevier, 1975.
|
|
[2] Alcubierre, M., Introduction to 3+1 Numerical Relativity, International Series of Monographs on Physics 140, Oxford University Press, 2008.
|
|
[3] Baumgarte, T.W., Shapiro, S.L., Numerical Relativity, Solving Einstein's Equations on a Computer, Cambridge University Press, 2010.
|
|
|
Compatibility
|
|
•
|
The Physics[ThreePlusOne][D3_] command was introduced in Maple 2017.
|
|
|