Introduction to Tensor Calculus

A.V.Smirnov
c Draft date November 6, 2003
Contents
1 Coordinates and Tensors 4

2 Cartesian Tensors 7
2.1 Tensor Notation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.2 Tensor Derivatives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

3 Curvilinear coordinates 15
3.1 Tensor invariance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
3.2 Covariant differentiation . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.3 Orthogonal coordinates . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
3.3.1 Unit vectors and stretching factors . . . . . . . . . . . . . . . . . 21
3.3.2 Physical components of tensors . . . . . . . . . . . . . . . . . . 24

4 Problems 27

A Solutions to problems 30

Bibliography 45

Index 46

2
Preface

This material offers a short introduction to tensor calculus. It is directed toward

students of continuum mechanics and engineers. The emphasis is made on ten-
sor notation and invariant forms. A knowledge of calculus is assumed. A more
complete coverage of tensor calculus can be found in [1, 2].

Nomenclature

A B A is defined as B, or A is equivalent to B
Ai Bi ∑3i Ai Bi . Note: Ai Bi A j B j

Ȧ partial derivative over time: ∂∂tA

Ai
 partial derivative over xi : ∂∂ xAi
V control volume
t time 

xi i-th component of a coordinate (i=0,1,2), or xi xuz

  

RHS Right-hand-side
LHS Left-hand-side
PDE Partial differential equation
.. Continued list of items

3
 There are two aspects of tensors that are of practical and fundamental im-
portance: tensor notation and tensor invariance. Tensor notation is of great prac-
tical importance, since it simplifies handling of complex equation systems. The
idea of tensor invariance is of both practical and fundamental importance, since it
provides a powerful apparatus to describe non-Euclidean spaces in general and
curvilinear coordinate systems in particular.
A definition of a tensor is given in Section 1. Section 2 deals with an im-
portant class of Cartesian tensors, and describes the rules of tensor notation.
Section 3 provides a brief introduction to general curvilinear coordinates, invari-
ant forms and the rules of covariant differentiation.

1 Coordinates and Tensors

Consider a space of real numbers of dimension n, Rn , and a single real time,

t. Continuum properties in this space can be described by arrays of different
dimensions, m, such as scalars (m 0), vectors (m 1), matrices (m 2), and

general multi-dimensional arrays. In this space we shall introduce a coordinate

system, xi i 1

n , as a way of assigning n real numbers1 for every point of space


There can be a variety of possible coordinate systems. A general transformation

rule between the coordinate systems is


x̃i

x̃i x1  xn  (1)

Consider a small displacement dxi . Then it can be transformed from coordi-

nate system xi to a new coordinate system x̃i using the partial differentiation rules
applied to (1):

∂ x̃i j
d x̃i dx (2)
∂xj

This transformation rule2 can be generalized to a set of vectors that we shall call
contravariant vectors:

∂x̃i j
Ãi A (3)
∂xj

1
Super-indexes denote components of a vector (i  1   n) and not the power exponent, for the reason
explained later (Definition 1.1)
2 The repeated indexes imply summation (See. Proposition 21)

4
That is, a contravariant vector is defined as a vector which transforms to a new
coordinate system according to (3). We can also introduce the transformation
matrix as:

∂ x̃i
aij (4)
∂xj
With which (3) can be rewritten as:

Ai

aij A j (5)

Transformation rule (3) will not apply to all the vectors in our space. For
example, a partial derivative ∂ ∂ xi will transform as:

∂ ∂ ∂xj ∂xj ∂
(6)
∂ x̃i ∂ x̃i ∂ x j ∂ x̃i ∂ x j

that is, the transformation coefficients are the other way up compared to (2). Now
we can generalize this transformation rule, so that each vector that transforms
according to (6) will be called a Contravariant vector:

∂xj
Ãi Aj (7)
∂x̃i

This provides the reason for using lower and upper indexes in a general
tensor notation.

Definition 1.1 Tensor

Tensor of order m is a set of nm numbers identified by m integer indexes. For
example, a 3rd order tensor A can be denoted as Ai jk and an m-order tensor can
be denoted as Ai1

im . Each index of a tensor changes between 1 and n. For ex-
ample, in a 3-dimensional space (n=3) a second order tensor will be represented
by 32 9 components.

Each index of a tensor should comply to one of the two transformation rules:
(3) or (7). An index that complies to the rule (7) is called a covariant index and is
denoted as a sub-index, and an index complying to the transformation rule (3) is
called a contravariant index and is denoted as a super-index.

Each index of a tensor can be covariant or a contravariant, thus tensor Akij

is a 2-covariant, 1-contravariant tensor of third order.

5
 Tensors are usually functions of space and time:


A i1

im

Ai1

im x1   xn t  

which defines a tensor field, i.e. for every point xi and time t there are a set of mn
nubers Ai1

im .

Remark 1.2 Tensor character of coordinate vectors

Note, that the coordinates xi are not tensors, since generally, they are not
transformed as (5). Transformation law for the coordinates is actually given by (1).
Nevertheless, we shall use the upper (contravariant) indexes for the coordinates.

Definition 1.3 Kronecker delta tensor

Second order delta tensor, δi j is defined as

i

j δi j

1
i


j δi j

0 (8)

From this definition and since coordinates xi are independent of each other
it follows that:

∂ xi
δi j (9)
∂xj

Corollary 1.4 Delta product

From the definition (1.3) and the summation convention (21), follows that

δi j A j

Ai (10)

Assume that there exists the transformation inverse to (5), which we call bij :

dxi

bij d x̃ j (11)

Then by analogy to (4) bij can be defined as:

∂ xi
bij (12)
∂ x̃ j

6
 j
From this relation and the independence of coordinates (9) it follows that aij bk

j
bij ak

δik , namely:

j ∂ x̃i ∂ x j
aij bk
∂ x j ∂ x̃k

∂ x j ∂ x̃i ∂ x̃i
δik (13)
∂ x j ∂ x̃k ∂ x̃k

2 Cartesian Tensors

Cartesian tensors are a sub-set of general tensors for which the transformation
matrix (4) satisfies the following relation:

∂ x̃k ∂ x̃k
aki akj δi j (14)
∂ xi ∂ x j

For Cartesian tensors we have

∂ x̃i ∂ xk
(15)
∂ xk ∂ x̃i

(see Problem 4.3), which means that both (5) and (6) are transformed with the
same matrix aik . This in turn means that the difference between the covariant and
contravariant indexes vanishes for the Cartesian tensors. Considering this we
shall only use the sub-indexes whenever we deal with Cartesian tensors.

2.1 Tensor Notation

Tensor notation simplifies writing complex equations involving multi-dimensional

objects. This notation is based on a set of tensor rules. The rules introduced in
this section represent a complete set of rules for Cartesian tensors and will be
extended in the case of general tensors (Sec.3). The importance of tensor rules
is given by the following general remark:

Remark 2.1 Tensor rules Tensor rules guarantee that if an expression follows
these rules it represents a tensor according to Definition 1.1.

7
 Thus, following tensor rules, one can build tensor expressions that will pre-
serve tensor properties of coordinate transformations (Definition 1.1) and coordi-
nate invariance (Section 3).
Tensor rules are based on the following definitions and propositions.

Definition 2.2 Tensor terms

A tensor term is a product of tensors.

For example:

Ai jk B jkC pq Eq Fp (16)

Definition 2.3 Tensor expression

Tensor expression is a sum of tensor terms. For example:

Ai jk B jk Ci D pq Eq Fp (17)

Generally the terms in the expression may come with plus or minus sign.

Proposition 2.4 Allowed operations

The only allowed algebraic operations in tensor expressions are the addi-
tion, subtraction and multiplication. Divisions are only allowed for constants, like
1 C. If a tensor index appears in a denominator, such term should be redefined,
so as not to have tensor indexes in a denominator. For example, 1 Ai should be
redefined as: Bi 1 Ai .

Definition 2.5 Tensor equality

Tensor equality is an equality of two tensor expressions.

For example:

Ai j B j

Cikp Dk E p E jC jki Bk (18)

8
Definition 2.6 Free indexes
A free index is any index that occurs only once in a tensor term. For exam-
ple, index i is a free index in the term (16).

Proposition 2.7 Free index restriction

Every term in a tensor equality should have the same set of free indexes.

For example, if index i is a free index in any term of tensor equality, such as
(18), it should be the free index in all other terms. For example

Ai j B j

C jD j

is not a valid tensor equality since index i is a free index in the term on the
RHS but not in the LHS.

Definition 2.8 Rank of a term

A rank of a tensor term is equal to the number of its free indexes.

For example, the rank of the term Ai jk B jCk is equal to 1.

It follows from (2.7) that ranks of all the terms in a valid tensor expression
should be the same. Note, that the difference between the order and the rank is
that the order is equal to the number of indexes of a tensor, and the rank is equal
to the number of free indexes in a tensor term.

Proposition 2.9 Renaming of free indexes

Any free index in a tensor expression can be named by any symbol as long
as this symbol does not already occur in the tensor expression.

For example, the equality

Ai j B j

Ci D j E j (19)

is equivalent to

Ak j B j

Ck D j E j (20)

Here we replaced the free index i with k.

9
Definition 2.10 Dummy indexes
A dummy index is any index that occurs twice in a tensor term.

For example, indexes j k p q in (16) are dummy indexes.

  

Proposition 2.11 Summation rule

Any dummy index implies summation, i.e.

n
Ai Bi

∑ Ai Bi (21)
i

Proposition 2.12 Summation rule exception If there should be no summation

over the repeated indices, it can be indicated by enclosing such indices in paren-
theses.

For example, expression:

C i A i Bj





Di j

does not imply summation over i.

Corollary 2.13 Scalar product


A scalar product notation from vector algebra: A B  is expressed in tensor


notation as Ai Bi .

The scalar product operation is also called a contraction of indexes.

Proposition 2.14 Dummy index restriction

No index can occur more than twice in any tensor term.

Remark 2.15 Repeated indexes

In case if an index occurs more than twice in a term this term should be
redefined so as not to contain more than two occurrences of the same index. For
example, term Aik B jkCk should be rewritten as Aik D jk , where D jk is defined as
D jk B j k C k with no summation over k in the last term.

10
Proposition 2.16 Renaming of dummy indexes
Any dummy index in a tensor term can be renamed to any symbol as long
as this symbol does not already occur in this term.

For example, term Ai Bi is equivalent to A j B j , and so are terms Ai jk B jCk and

Aipq B pCq .

Remark 2.17 Renaming rules

Note that while the dummy index renaming rule (2.16) is applied to each
tensor term separately, the free index naming rule (2.9) should apply to the whole
tensor expression. For example, the equality (19) above

Ai j B j

Ci D j E j

can also be rewritten as

Akp B p

Ck D j E j (22)

without changing its meaning.

(See Problem 4.1).

Definition 2.18 Permutation tensor

The components of a third order permutation tensor εi jk are defined to be
equal to 0 when any index is equal to any other index; equal to 1 when the set of
indexes can be obtained by cyclic permutation of 123; and -1 when the indexes
can be obtained by cyclic permutation from 132. In a mathematical language it
can be expressed as:

i j i k j k εi jk 0


i jk PG 123  εi jk 1




i jk PG 132  εi jk

1 (23)
 
where PG abc  is a permutation group of a triple of indexes abc, i.e. PG abc 


abc bca cab . For example, the permutation group of 123 will consist of three
  

combinations: 123, 231 and 312, and the permutation group of 123 consists of
132, 321 and 213.

11
Corollary 2.19 Permutation of the permutation tensor indexes
From the definition of the permutation tensor it follows that the permutation
of any of its two indexes changes its sign:

εi jk
 εik j (24)

A tensor with this property is called skew-symmetric.

Corollary 2.20 Vector product

A vector product (cross-product) of two vectors in vector notation is ex-
pressed as

A

B 

C (25)

which in tensor notation can be expressed as

Ai

εi jk B jCk (26)

Remark 2.21 Cross product

Tensor expression (26) is more accurate than its vector counterpart (25),
since it explicitly shows how to compute each component of a vector product.

Theorem 2.22 Symmetric identity

If Ai j is a symmetric tensor, then the following identity is true:

εi jk A jk

0 (27)

Proof:
From the symmetry of Ai j we have:

εi jk A jk

εi jk Ak j (28)

Let’s rename index j into k and k into j in the RHS of this expression, according
to rule (2.16):

εi jk Ak j

εik j A jk

12
 Using (24) we finally obtain:

εik j A jk
 εi jk A jk

Comparing the RHS of this expression to the LHS of (28) we have:

εi jk A jk
 εi jk A jk

from which we conclude that (27) is true.

Theorem 2.23 Tensor identity

The following tensor identity is true:

εi jk εipq

δ j p δkq  δ jq δkp (29)

Proof
This identity can be proved by examining the components of equality (29)
component-by-component.

Corollary 2.24 Vector identity

Using the tensor identity (29) it is possible to prove the following important
vector identity:

 

 

A 

B 

C

B A C  C A B  (30)

See Problem 4.4.

2.2 Tensor Derivatives

For Cartesian tensors derivatives introduce the following notation.

Definition 2.25 Time derivative of a tensor

A partial derivative of a tensor over time is designated as

∂A
Ȧ
∂t

13
Definition 2.26 Spatial derivative of a tensor
A partial derivative of a tensor A over one or its spacial components is de-
noted as A i :


∂A
Ai (31)
∂ xi


that is, the index of the spatial component that the derivation is done over is
delimited by a comma (’,’) from other indexes. For example, Ai j k is a derivative of


a second order tensor Ai j .

Definition 2.27 Nabla

Nabla operator acting on a tensor A is defined as

∇i A Ai
 (32)

Even though the notation in (31) is sufficient to define the derivative, In some
instances it is convenient to introduce the nabla operator as defined above.

Remark 2.28 Tensor derivative

In a more general context of non-Cartesian tensors the coordinate indepen-
dent derivative will have a different form from (31). See the chapter on covariant
differentiation in [1].

Remark 2.29 Rank of a tensor derivative

The derivative of a zero order tensor (scalar) as given by (31) forms a first
order tensor (vector). Generally, the derivative of an m-order tensor forms an m 1
order tensor. However, if the derivation index is a dummy index, then the rank of
the derivative will be lower than that of the original tensor. For example, the rank
of the derivative Ai j j is one, since there is only one free index in this term.


Remark 2.30 Gradient

Expression (31) represents a gradient, which in a vector notation is ∇ A:

∇A  Ai 

14
Corollary 2.31 Derivative of a coordinate
From (9) it follows that:

xi 

j

δi j (33)

In particular, the following identity is true:

xi i


x1 1
 x2 2
 x3 3 

1 1 1

3 (34)

Remark 2.32 Divergence operator

A divergence operator in a vector notation is represented in a tensor notation
as Ai i : 


∇ A 

 Ai i

Remark 2.33 Laplace operator

The Laplace operator in vector notation is represented in tensor notation as
A ii :


∆A  A ii 

Remark 2.34 Tensor notation

Examples (2.30), (2.32) and (2.33) clearly show that tensor notation is more
concise and accurate than vector notation, since it explicitly shows how each
component should be computed. It is also more general since it covers cases
that don’t have representation in vector notation, for example: Aik k j . 

3 Curvilinear coordinates

In this section 3 we introduce the idea of tensor invariance and introduce the rules
for constructing invariant forms.
3
In this section we reinstall the difference between covariant and contravariant indexes.

15
3.1 Tensor invariance

The distance between the material points in a Cartesian coordinate system is

computed as dl 2 dxi dxi . The metric tensor, gi j is introduced to generalize the

notion of distance (39) to curvilinear coordinates.

Definition 3.1 Metric Tensor

The distance element in curvilinear coordinate system is computed as:

dl 2

gi j dxi dx j (35)

where gi j is called the metric tensor.

Thus, if we know the metric tensor in a given curvilinear coordinate system then
the distance element is computed by (35). The metric tensor is defined as a
tensor since we need to preserve the invariance of distance in different coordinate
systems, that is, the distance should be independent of the coordinate system,
thus:

dl 2

gi j dxi dx j

g̃i j d x̃i d x̃ j (36)

The metric tensor is symmetric, which can be shown by rewriting (35) as

follows:

gi j dxi dx j

gi j dx j dxi

g ji dxi dx j

where we first swapped places of dxi and dx j , and then renamed index i into j
and j into i. We can rewrite the equality above as:

gi j dxi dx j  g ji dxi dx j


gi j  g ji  dxi dx j

0

Since the equality above should hold for any dxi dx j , we get:

gi j

g ji (37)

The metric tensor is also called the fundamental tensor. The inverse of the
metric tensor is also called the conjugate metric tensor, gi j , which satisfies the
relation:

16
 gik gk j

δi j (38)

Let xi be a Cartesian coordinate system, and x̃ j - the new curvilinear coordi-

nate system. Both systems are related by transformation rules (5) and (11). Then
from (36) we get:

∂ xi j ∂ xi k ∂ xi ∂ xi j k
dl 2 dxi dxi d x̃ k d x̃ d x̃ d x̃ (39)
∂x̃ j ∂x̃ ∂x̃ j ∂x̃k

When we transform from a Cartesian to curvilinear coordinates the metric

tensor in curvilinear coordinate system, g̃i j can be determined by comparing re-
lations (39) and (35):

∂ xk ∂ xk
g̃i j (40)
∂ x̃i ∂ x̃ j

Using (38) we can also find its inverse as:

∂x̃i ∂x̃ j
g̃i j (41)
∂ xk ∂ xk

Using these expression one can compute gi j and gi j in various curvilinear coordi-
nate systems (see Problem 4.6).

Definition 3.2 Conjugate tensors

For each index of a tensor we introduce the conjugate tensor where this
index is transfered to its counterpart (covariant/contravariant) using the relations:

Ai

gi j A j (42)
Ai

gi j A j (43)

Conjugate tensor is also called the associate tensor. Relations (42), (43)
are also called as operations of raising/lowering of indexes.

Remark 3.3 Tensor invariance

Since the transformation rules defined by (1.1) have a simple multiplicative
character, any tensor expression should retain it’s original form under transforma-
tion into a new coordinate system. Thus if an expression is given in a tensor form
it will be invariant under coordinate transformations.

17
 Not all the expressions constructed from tensor terms in curvilinear coordi-
nates will be tensors themselves. For example, if vectors Ai and Bi are tensors,
then Ai Bi is not generally a tensor4 . However, if we consider the same operation
on a contravariant tensor Ai and a covariant tenso Bi then the product will form an
invariant:

Āi B̄i

Ai Bi (44)

Thus in curvilinear coordinates we have to refine the definition of the scalar

product (Corollary 2.13) or the index contraction operation to make it invariant
(Problem 4.12).

Definition 3.4 Invariant Scalar Product

The invariant form of the scalar product between two covariant vectors Ai
and Bi is gi j Ai B j . Similarly, the invariant form of a scalar product between two
contravariant vectors Ai and Bi is gi j Ai B j , where gi j is the metric tensor (40) and
gi j is its conjugate (38).

Corollary 3.5 Two forms of a scalar product

According to (42), (43) the scalar product can be represented by two invari-
ant forms: Ai Bi and Ai Bi . It can be easily shown that these two forms have the
same values (see Problem 4.12).

Corollary 3.6 Rules of invariant expressions

To build invariant tensor expressions we add two more rules to Cartesian
tensor rules outlined in Section 2.1:

1. Each free index should keep its vertical position in every term, i.e. if the
index is covariant in one term it should be covariant in every other term, and
vise versa.
2. Every pair of dummy indexes should be complementary, that is one should
be covariant, and another contravariant.

For example, a Cartesian formulation of a momentum equation for an in-

compressible viscous fluid is

u̇i uk ui k





Pi
ρ


ντik k


4 For Cartesian tensors any product of tensors will always be a tensor, but this is not so for general tensors

18
 The invariant form of this equation is:

u̇i uk ui k 




Pi
ρ


ντki k 
(45)

where the rising of indexes was done using relation (42): uk

gk j u j , and τki

g k j τi j .

3.2 Covariant differentiation

A simple scalar value, S, is invariant under coordinate transformations. A partial

derivative of an invariant is a first order covariant tensor (vector):

∂S
Ai Si
∂ xi


However, a partial derivative of a tensor of the order one and greater is not
generally an invariant under coordinate transformations of type (7) and (3).
In curvilinear coordinate system we should use more complex differentiation
rules to preserve the invariance of the derivative. These rules are called the rules
of covariant differentiation and they guarantee that the derivative itself is a tensor.
According to these rules the derivatives for covariant and contravariant indices
will be slightly different. They are expressed as follows:

Ai 

j
∂ Ai
∂xj
 k
ij

Ak (46)
∂ Ai
Ai j i
Ak (47)
∂xj kj


where the contstruct k is defined as

ij

k
ij



1 kl
2
g 
∂ gil
∂xj
∂ g jl
∂ xi
 ∂ gi j
∂ xl 

and is also known in tensor calculus as Christoffel’s symbol of the second kind
[1]. Tensor gi j represents the inverse of the metric tensor gi j (38). As can be seen
differentiation of a single component of a vector will involve all other components
of this vector.

19
 In differentiating higher order tensors each index should be treated inde-
pendently. Thus differentiating a second order tensor, Ai j , should be performed
as:

Ai j k 

∂ Ai j
∂ xk
 

m
ik  Am j  

m
jk  Aim

and as can be seen also involves all the components of this tensor. Likewise for
the contravariant second order tensor Ai j we have:

ij ∂ Ai j i j
Am j Aim
 

Ak (48)
∂ xk mk mk
 

And for a general n-covariant, m-contravariant tensor we have:

j

j ∂ j1

jm
Ai11

inmp A
∂ x p i1

in k


j1 q j2

jm 

jm j1

jm 1 q
qp Ai1

in
   

qp Ai1

in


q j1

jm q j

j
Ai11

inm 1 q
 

i1 p Aqi2

in    

in p  (49)

Despite their seeming complexity, the relations of covariant differentiation

can be easily implemented algorithmically and used in numerical solutions on
arbitrary curved computational grids (Problem 4.8).

Remark 3.7 Rules of invariant expressions

As was pointed out in Corollary 3.6, the rules to build invariant expressions
involve raising or lowering indexes (42), (43). However, since we did not intro-
duce the notation for contravariant derivative, the only way to raise the index of a
covariant derivative, say A i , it to use the relation (42) directly, that is: gi j A j .
 

For example, we can re-formulate the momentum equation (45) in terms of

contravariant free index i as:

u̇ i k i
u uk 

 gik P k
ρ


ντikk 
(50)

where the index of the pressure term was raised by means of (42).
Using the invariance of the scalar product one can construct two important
differential operators in curvilinear coordinates: divergence of a vector div A Aii 

(51) and Laplacian, ∆ A gik A ki (55). 

20
Definition 3.8 Divergence
Divergence of a vector is defined as Aii : 

div A A ii 
(51)

From this definition and the rule of covariant differentiation (47) we have:

∂ Ai
A ii i
Ak


(52)
∂ xi ki 

this can be shown [2] to be equal to:

i ∂ Ai 1 ∂
Ai 

g Ai
∂ xi g ∂x i


1 ∂
gAi 

(53)
g ∂ xi

where g is the determinant of the metric tensor gi j .

The divergence of a covariant vector Ai is defined as a divergence of its
conjugate contravariant tensor (42):

A ii 

gi j A j i  (54)

Definition 3.9 Laplacian

A Laplace operator or a Laplacian of a scalar A is defined as

∆A gik A ki 
(55)

The definitions (3.8), (3.9) of differential operators are invariant under coor-
dinate transformations. They can be programmed using a symbolic manipulation
packages and used to derive expressions in different curvilinear coordinate sys-
tems (Problem 4.9).

3.3 Orthogonal coordinates

3.3.1 Unit vectors and stretching factors

The coordinate system is orthogonal if the tangential vectors to coordinate lines

are orthogonal at every point.

21
 Consider three unit vectors, ai bi ci , each directed along one of the coordi-
 

nate axis (tangential unit vectors), that is:

ai a1 0 0


   (56)
bi 0 b2 0


   (57)
ci 0 0 c3


   (58)

The condition of orthogonality means that the scalar product between any
two of these unit vectors should be zero. According to the definition of a scalar
product (Definition 3.4) it should be written in form (44), that is, a scalar product
between vectors ai and bi can be written as: ai bi or ai bi . Let’s use the first form for
definiteness. Then, applying the operation of rising indexes (42), we can express
the scalar product in contravariant components only:

0

a i bi

gi j ai b j

1 1 2
g11 a 0 g12 a b g13 00
g21 a2 b1 g22 0b2 g23 00
g31 a3 0 g32 0b2 g33 00


g12 g21  a1 b2 2g12 a1 b2 0


(59)

where we used the symmetry of gi j , (37). Since vectors a1 and b2 were chosen to
be non-zero, we have: g12 0. Applying the same reasoning for scalar products of

other vectors, we conclude that the metric tensor has only diagonal components
non-zero5 :

gi j

δi j g ii

(60)

Let’s introduce stretching factors, hi , as the square roots of these diagonal com-
ponents of gi j :

  
1 2 1 2 1 2
h1 g11  ; h2 g22  ; h3 g33  ; (61)

Now, consider the scalar product of each of the unit vectors (56)-(58) with
itself. Since all vectors are unit, the scalar product of each with itself should be
one:
5 We use parenthesis to preclude summation (Proposition 2.12)

22
 ai ai

bi bi

ci ci

1

Or, expressed in contravariant components only the condition of unity is:

gi j ai a j

gi j bi b j

g i j ci c j

1

Now, consider the first term above and substitute the components of a from (56).
The only non-zero term will be:


2
g11 a1 a1

h1  a1  2


1

and consequently:

a1



1
h1
(62)

where the negative solution identifies a vector directed into the opposite direction,
and we can neglect it for definiteness. Applying the same reasoning for each of
the tree unit vectors ai bi ci , we can rewrite (56), (57) and (58) as:
 

1
ai


00    (63)
h1

1
bi


0 0   (64)
h2

1
ci


00    (65)
h3

which means that the components of unit vectors in a curved space should be
scaled with coefficients hi . It follows from this that the expression for the element
of length in curvilinear coordinates, (35), can be written as:


dl 2

gi j d x̃i d x̃ j

h2i d x̃i  2
(66)

Similarly, we introduce the hi coefficients for the conjugate metric tensor

(38):


gi j

δi j h i

 2
(67)

Combining the latter with (38), we obtain: δi j h i h i

δi j , from which it follows that

23
 hi



1 hi

(68)

3.3.2 Physical components of tensors

Consider a direction in space determined by a unit vector ei . Then the physical

component of a vector Ai in the direction ei is given by a scalar product between
Ai and ei (Definition 3.4), namely:


A e

g i j Ai e j

According to Corollary 3.5 the above can also be rewritten as:


A e

A i ei

Ai e i (69)

Suppose the unit vector is directed along one of the axis: ei e1 0 0 . From


  

(63) it follows that:

e1

1 h1

where h1 is defined by (61). Thus according to (69) the physical component of

vector Ai in direction 1 in orthogonal coordinate system is equal to:


A 1

A 1 h1

or, repeating the argument for other components, we have for the physical com-
ponents of a covariant vector:

A1 h 1 A2 h 2 A3 h 3
  (70)

Following the same reasoning, for the contravariant vector Ai , we have:

h 1 A1 h 2 A2 h 3 A3 

General rules of covariant differentiation introduced in (Sec.3.2) simplify

considerably in orthogonal coordinate systems. In particular, we can define the
nabla operator by the physical components of a covariant vector composed of
partial differentials:

24
 1 ∂
∇i (71)
h i ∂ xi

where the parentheses indicate that there’s no summation with respect to index i.
In orthogonal coordinate system the general expressions for divergence (53)
and Laplacian (55)) operators can be expressed in terms of stretching factors only
[3]:

1 ∂ H
Aii Ai (72)
H ∂ xi hi


1 ∂ H A
∆A
H ∂ xi h i ∂ xi

n
H ∏ hi
i 1

Important examples of orthogonal coordinate systems are spherical and cylindri-

cal coordinate systems. Consider the example of a cylindrical coordinate system:
rθ l :
 

xi

x1 x2 x3 and x̃i
  

  

x1

r cos θ
x2

r sin θ
x3 l

According to (40) only few components of the metric tensor will survive
(Problem 4.5). Then we can compute nabla, divergence and Laplacian oper-
ators according to (71), (52) and (55), or using simplified relations (72)-(73):

∂ 1 ∂ ∂
∇ 

∂r r ∂θ ∂z
 

∂ A1 1 ∂ A2 ∂ A3 1
div A A1
∂ x̃1 x̃1 ∂ x̃2 ∂ x̃3 x̃1

∂ Ar 1 ∂ Aθ ∂ Az 1
Ar
∂r r ∂θ ∂z r

Note, that instead of using the contravariant components as implied by the gen-
eral definition of the divergence operator (51) we are using the covariant compo-
nents as dictated by relation (70). The expression of the Laplacian becomes:

25
 ∂2 A 1 ∂2 A ∂2 A 1 ∂A
∆A   
∂ x̃1  2 x̃21 ∂ x̃2  2 ∂ x̃3  2 x̃1 ∂ x̃1

∂2 A 1 ∂2 A ∂2 A 1 ∂A
  
∂ r 2 r2 ∂ θ  2 ∂ z 2 r ∂r

(see Problems 4.9,4.10).

The advantages of the tensor approach are that it can be used for any type
of curvilinear coordinate transformations, not necessarily analytically defined, like
cylindrical (84) or spherical. Another advantage is that the equations above can
be easily produced automatically using symbolic manipulation packages, such
as Mathematica (wolfram.com) (Problems 4.6,4.7,4.9). For further reading see
[1, 2].

26
4 Problems

Problem 4.1 Check tensor expressions for consistency

Check if the following Cartesian tensor expressions violate tensor rules:

Ai j B jk B pqCq Dk

0

E pqi Fk jC pk B p j D jq Gq

Fkp

Ei jk A j Bk  D i j Ai B j

Fi j G jk Hk j

Problem 4.2 Construct tensor expression

Construct a valid Cartesian tensor expression, consisting of three terms,
each including some of the four tensors: Ai jk Bi j Ci Di j . Term 1 should include   

tensors A B C only, term 2 tensor B C D and term 3 tensors C D A. The ex-

     

pression should have 2 free indexes, which should always come first among the
indexes of a tensor. The free indexes should be at A and B in the first term, at
B and C in the second term and C and D in the last term. How many different
tensor expressions can be constructed?

Problem 4.3 Cartesian identity

Prove identity (15)

Problem 4.4 Vector identity

Using tensor identity (29):

εi jk εipq

δ j p δkq  δ jq δkp

prove vector identity (30):

 

 

A 

B 

C

B A C  C A B 

Problem 4.5 Metric tensor in cylindrical coordinates

r θ l (84) is given by the following trans-


Cylindrical coordinate system yi

  

formation rules to a Cartesian coordinate system, xi xyz :

  

27
 x

r cos θ
y

r sin θ
z l

Obtain the components of the metric tensor (40) gi j and its inverse gi j (38)
in cylindrical coordinates.

Problem 4.6 Metric tensor in curvilinear coordinates

Using Mathematica Compute the metric tensor, g, (40) and its conjugate, ĝ,
(38) in spherical coordinate system (r φ θ):  

x

r sin θ cos φ
y

r sin θ sin φ
z r cos θ


(73)

Problem 4.7 Christoffel’s symbols with Mathematica

Using the Mathematica package, write the routines for computing Christof-
fel’s symbols.

Problem 4.8 Covariant differentiation with Mathematica

Using the Mathematica package, and the routines developed in Problem 4.7
write the routines for covariant differentiation of tensors up to second order.

Problem 4.9 Divergence of a vector in curvilinear coordinates

Using the Mathematica package and the solution of Problem 4.8, write the
routines for computing divergence of a vector in curvilinear coordinates.

Problem 4.10 Laplacian in curvilinear coordinates

Using the Mathematica package and the solution of Problem 4.8, write the
routines for computing the Laplacian in curvilinear coordinates.

Problem 4.11 Invariant expressions

Check if any of these tensor expressions are invariant, and correct them if
not:

28
 j
Ai BijkCt k


Dt (74)

ij
A jk BipqC kq  Fk j Gkp H j

jq
H k Ak j Cti B pit q
 (75)

p j
E i Bikp DkqC jq

Dki G p i  (76)

Problem 4.12 Contraction invariance

Prove that Ai Bi is an invariant and Ai Bi is not.

29
A Solutions to problems

Problem 4.1: Check tensor expressions

Check if the following Cartesian tensor expressions violate tensor rules:

Ai j B jk B pqCq Dk

0

Answer: term (1): ik = free, term (2): pk=free

E pqi Fk jC pk B p j D jq Gq

Fkp

Answer: (1): ijq=free (2): p=free (3): kp=free

Ei jk A j Bk  D i j Ai B j

Fi j G jk Hk j

Answer: (1): i=free (2): none, (3): i=free, j = tripple occurrence

Problem 4.2: Construct tensor expression

Construct a valid Cartesian tensor expression, consisting of three terms,
each including some of the four tensors: Ai jk Bi j Ci Di j . Term 1 should include
  

tensors A B C only, term 2 tensor B C D and term 3 tensors C D A. The ex-

     

pression should have 2 free indexes, which should always come first among the
indexes of a tensor. The free indexes should be at A and B in the first term, at
B and C in the second term and C and D in the last term. How many different
tensor expressions can be constructed?
Solution
One possibility is:

Aipk B jkC p BiqC pC j D pq Ci D j p A pqq

0

Since there are four locations for dummy indexes in each term, there could
be three different combinations of dummies in each term. Thus, the total number
of different expression is 33 27

30
Problem 4.3: Cartesian identity
Prove identity (15).
Proof
Integrating (5) in the case of constant transformation marix coefficients, we
have:

x̃i

aik xk bi (77)

where the transformation matrix is given by (4):

∂ x̃i
aik (78)
∂ xk
By the definition of the Cartesian coordinates (79) we have:

∂ x̃k ∂ x̃k
aki akj δi j (79)
∂ xi ∂ x j

Let’s multiply the transformation rule (77) by aij . Then we get:

aij x̃i

aij aik xk aij bi

δ jk xk aij bi

xj aij bi

Differentiation this over x̃i , we have:

∂xj
aij
∂x̃i

Now rename index j into k:

∂ xk
aik
∂x̃i

Comparing this with (78), we have

∂x̃i ∂xj
∂xj ∂x̃i

which proves (15).

31
Problem 4.4: Tensor identity
Prove identity (30).
Proof

 

 

A 

B 

C

B A C  C A B  (80)

Proof
Applying (26) twice to the RHS of (80), we have:



A B C  

εi jk A j εkpq B pCq


εi jk εkpq A j B pCq

From (24) it follows that εi jk
 εik j

εki j . Then we have:

εi jk εkpq A j B pCq

εki j εkpq A j B pCq (81)

Now rename the dummy indexes: k ii  j j  k, so that the expression

looks like one in (29):


εi jk εipq  Ak B pCq

δ j p δkq δ jq δkp  Ak B pCq 


δ j p B p δkq AkCq δ jqCq δkp Ak B p (82)

Using (10), and since A j

B j is the same as Ai

Bi the latter can be rewrit-
ten as:

B j AqCq  C j A pB p (83)

which is the same as

 

B A C  C A B 

32
Problem 4.5: Metric tensor in cylindrical coordinates.
Cylindrical coordinate system x̃i r θ l (84) is given by the following trans-


  

formation rules to a Cartesian coordinate system, xi



xyz :

  

x

r cos θ
y

r sin θ
z l

Obtain the components of the metric tensor (40) gi j and its inverse gi j (38)
in cylindrical coordinates.
Solution:
First compute the derivatives of xi x y z with respect to x̃i rθ l :
 

     

∂x1 ∂x
xr cos θ
∂x̃1 ∂r

∂x2 ∂y
yr sin θ
∂x̃1 ∂r

∂x1
∂x̃ 2


∂x
∂θ
xθ 


r sin θ
∂x2 ∂y
yθ r cos θ
∂x̃2 ∂ θ

∂x3 ∂ z
zl 1 (84)
∂x̃3 ∂ z

Then the components of the metric tensor are:

grr

xr xr yr yr 1

gθ θ

xθ xθ yθ yθ r 2

gzz 1

grr 1

1
gθ θ
r2

zz
g 1

33
Problem 4.6: Metric tensor in curvilinear coordinates
Using Mathematica, write a procedure to compute metric tensor in curvilin-
ear coordinate system, and use it to obtain the components of metric tensor, g,
(40) and its conjugate, ĝ, (38) in spherical coordinate system (r φ θ):
 

x

r sin θ cos φ
y

r sin θ sin φ
z r cos θ


(85)

Solution with Mathematica

NX = 3

(* Curvilinear cooridnate system *)

Y = Array[,NX] (* Spherical coordinate system *)
Y[[1]] = r; (* radius *)
Y[[2]] = th; (* angle theta *)
Y[[3]] = phi; (* angle phi *)

(* Cartesian coordinate system *)

X = Array[,NX]
X[[1]] = r Sin[th] Cos[phi];
X[[2]] = r Sin[th] Sin[phi];
X[[3]] = r Cos[th];

(* Compute the Jacobian: dXi/dYj *)

J = Array[,{NX,NX}]
Do[
J [[i,j]] = D[X[[i]],Y[[j]]],
{j,1,NX},{i,1,NX}
]

(* Covariant Metric tensor *)

g = Array[,{NX,NX}] (* covariant *)
Do[
g [[i,j]] = Sum[J[[k,i]] J[[k,j]],{k,NX}],
{j,1,NX},{i,1,NX}
];
g=Simplify[g]

34
(* Contravariant metric tensor *)
g1 =Array[,{NX,NX}]
g1=Inverse[g]

With the result:


0 r2 0 0 0 r2 sin θ  2
   

g

100
           


  

2 

csc θ  2
ĝ 100 0r 0 00
r2
           

Problem 4.7: Christoffel’s symbols with Mathematica

Using the Mathematica package, write the routines to compute Christoffel’s
symbols
Solution

(************* File g.m *************

The metric tensor

and Christoffel symbols

*************************************)
DIM = 3
(*
The metric tensor
*)
g = Array[,{DIM,DIM}] (* covariant *)
g1 =Array[,{DIM,DIM}] (* contravariant *)
Do[
g [[i,j]] = 0;
g1[[i,j]] = 0
,
{j,1,DIM},{i,1,DIM}
]
(*

35
 Cylindrical coordinates
*)
Z=Array[,DIM]
Z[[1]] = r
Z[[2]] = th
Z[[3]] = z
g [[1,1]] = 1
g [[2,2]] = rˆ2
g [[3,3]] = 1
g1[[1,1]] = 1
g1[[2,2]] = 1/rˆ2
g1[[3,3]] = 1
(*
Christoffel symbols of the first and second type
*)
Cr1 = Array[,{DIM,DIM,DIM}]
Cr2 = Array[,{DIM,DIM,DIM}]
Do[
Cr1[[i,j,k]] = 1/2
(
D[ g [[i,k]], Z[[j]] ]
+ D[ g [[j,k]], Z[[i]] ]
- D[ g [[i,j]], Z[[k]] ]
),
{k,DIM},{j,DIM},{i,DIM}
]
Do[
Cr2[[l,i,j]] =
Sum[
g1[[l,k]] Cr1[[i,j,k]],
{k,DIM}
],
{j,DIM},{i,DIM},{l,DIM}
]

Problem 4.8: Covariant differentiation with Mathematica

Using the Mathematica package, write the routines for covariant differentia-
tion of tensors up to second order.

36
 solution

(************** File D.m *******************

Rules of covariant differentiation

********************************************)
(*
B.Spain
Tensor Calculus, 1965
Eq.(22.2)
*)
D1[N_,A_,k_,X_,j_]:=
(*
Computes covariant derivative
of a mixed tensor of second order
with index k - covariant (upper)
*)
Module[
{i,s},
s = Sum[Cr2[[k,i,j]] A[[i]],{i,N}];
D[A[[k]],X[[j]]] + s
]
Dl1[N_,A_,l_,X_,t_]:=
(*
Computes covariant derivative
of a mixed tensor of second order
with index l - covariant (lower)
*)
Module[
{s,r},
s =Sum[Cr2[[r,l,t]] A[[r]],{r,N}];
D[A[[l]],X[[t]]] - s
]
D1l1[N_,A_,m_,l_,X_,t_]:=
(*
Computes covariant derivative
of a mixed tensor of second order
with index m - contravariant (upper) and
index l - covariant (lower)

37
*)
Module[
{s1,s2,r},
s1 =Sum[Cr2[[m,r,t]] A[[r,l]],{r,N}];
s2 =Sum[Cr2[[r,l,t]] A[[m,r]],{r,N}];
D[A[[m,l]],X[[t]]] + s1 - s2
]
D2[N_,A_,i_,j_,X_,n_]:=
(*
Computes covariant derivative
of second order tensor with
both m and l contravariant (upper)
indexes
B.Spain
Tensor Calculus, 1965
Eq.(23.3)
*)
Module[
{s1,s2,k},
s1 =Sum[Cr2[[i,k,n]] A[[k,j]],{k,N}];
s2 =Sum[Cr2[[j,k,n]] A[[i,k]],{k,N}];
D[A[[i,j]],X[[n]]] + s1 + s2
]
D2l1[N_,A_,i_,j_,k_,X_,n_]:=
(*
Computes covariant derivative
of third order tensor with
i and j contravariant (upper)
and k contravariant (lower)
indexes
B.Spain
Tensor Calculus, 1965
Eq.(23.3)
*)
Module[
{s1,s2,s3,m},
s1 =Sum[Cr2[[i,m,n]] A[[m,j,k]],{m,N}];
s2 =Sum[Cr2[[j,m,n]] A[[i,m,k]],{m,N}];
s3 =Sum[Cr2[[m,k,n]] A[[i,j,m]],{m,N}];
D[A[[i,j,k]],X[[n]]] + s1 + s2 - s3
]

38
D4l1[N_,A_,i1_,i2_,i3_,i4_,i5,X_,i6_]:=
(*
Computes covariant derivative
of 5 order tensor with
4 first indexes contravariant (upper)
and the last one contravariant (lower)
B.Spain
Tensor Calculus, 1965
Eq.(23.3)
*)
Module[
{k,s1,s2,s3,s4,s5},
s1= Sum[Cr2[[i1,k,n]] A[[k,i2,i3,i4,i5]],{k,N}];
s2= Sum[Cr2[[i2,k,n]] A[[i1,k,i3,i4,i5]],{k,N}];
s3= Sum[Cr2[[i3,k,n]] A[[i1,i2,k,i4,i5]],{k,N}];
s4= Sum[Cr2[[i4,k,n]] A[[i1,i2,i3,k,i5]],{k,N}];
s5=-Sum[Cr2[[k,i5,n]] A[[i1,i2,i3,i4,k]],{k,N}];
D[A[[i1,i2,i3,i4,i5]],X[[i6]]]+s1+s2+s3+s4+s5
]

Problem 4.9: Divergence of a vector in curvilinear coordinates

Using the Mathematica package and the solution of Problem 4.8, write the
routines for computing divergence of a vector in curvilinear coordinates.
Solution
Using the algorithms of covariant differentiation developed in Problem 4.8
we have:

<<"./g.m" (* The g-tensor and Christoffel symbols *)

<<"./D.m" (* Rules of covariant differentiation *)

(* The original coordinates: *)

NX = DIM
X = Array[,NX]

(* Variables: *)
NV = DIM
U = Array[,NV]

39
(* New coordinate system *)
Y = Array[,NX]
Y[[1]] = r;
Y[[2]] = th;
Y[[3]] = z;
X[[1]] = r Cos[th];
X[[2]] = r Sin[th];
X[[3]] = z;

(* Compute the Jacobian *)

J = Array[,{DIM,DIM}]
Do[
J [[i,j]] = D[X[[i]],Y[[j]]],
{j,1,DIM},{i,1,DIM}
]
J1=Simplify[Inverse[J]]

(* Derivatives of a vector *)

V0 = Array[,NX]
V0[[1]] = Vr[r,th,z];
V0[[2]] = Vt[r,th,z];
V0[[3]] = Vz[r,th,z];

(*
Rescaling for physical
(dimensionally correct) coordinates
(\cite[5.102-5.110]{SyScTC69})
*)
V = Array[,NX]
Do[
V[[i]] = PowerExpand[V0[[i]]/g[[i,i]]ˆ(1/2)],
{i,1,NX}
]

(*
Transform vectors
as first order contravariant tensors
*)
U = Array[,NX]

40
SetAttributes[RV1,HoldAll]
RV1[NX,V,U]
(*
Compute first covariant derivatives
of vectors
*)
DV = Array[,{NX,NX}];
Do[
DV[[i,j]] = D1[NX,V,i,Y,j],
{j,1,NX},{i,1,NX}
]
(* Divergence *)
div=0
Do[
div=div+DV[[i,i]],
{i,NX}
]
div0 = div/.th->0

Problem 4.10: Laplacian in curvilinear coordinates

Using the Mathematica package, write the routines for computing Laplacian
in curvilinear coordinates.
solution
Using the algorithms of covariant differentiation developed in Problem 4.8
we have:

<<"./g.m" (* The g-tensor and Christoffel symbols *)

<<"./D.m" (* Rules of covariant differentiation *)

(* The original coordinates: *)

NX = DIM
X = Array[,NX]

(* Variables: *)
NV = DIM
U = Array[,NV]

41
(* New coordinate system *)
Y = Array[,NX]
Y[[1]] = r;
Y[[2]] = th;
Y[[3]] = z;
X[[1]] = r Cos[th];
X[[2]] = r Sin[th];
X[[3]] = z;

(* Compute the Jacobian *)

J = Array[,{DIM,DIM}]
Do[
J [[i,j]] = D[X[[i]],Y[[j]]],
{j,1,DIM},{i,1,DIM}
]
J1=Simplify[Inverse[J]]

(* Derivative of a scalar *)

DP = Array[,NX];
Do[
DP[[i]] = D[p[r,th,z],Y[[i]]],
{i,1,NX}
]
DDP = Array[,{NX,NX}];
Do[
DDP[[i,j]] = Dl1[NX,DP,i,Y,j],
{i,1,NX},{j,1,NX}
]
DDQ = Array[,{NX,NX}];
Do[
DDQ[[i,j]] = Sum[DDP[[k,l]] J1[[k,i]] J1[[l,j]],{k,NX},{l,NX}],
{i,1,NX},{j,1,NX}
]

(* Laplacian *)
(*** lap=lap+Sum[g[[i,j]]*Dl1[NX,DS,j,Y,i],{i,1,NX},{j,1,NX}],*)
lap=Sum[DDQ[[i,i]],{i,NX}]
lap0=lap/.th->0

42
Problem 4.11: Invariant expressions
Check if any of these tensor expressions are invariant, and correct them if
not:

j
Ai BijkCt k 


Dt (86)

ij
A jk BipqC kq  Fk j Gkp H j

jq
H k Ak j Cti B pit q  (87)

p j
E i Bikp DkqC jq

Dki G p i  (88)

Answers:
A corrected form of (86) is:

j
Ai BikjCt k 


Dt

Equality (88) requires no corrections. A corrected form of (88) is:

jq
Ei Bikp D pkqC j

Dik G p i 

Since there are two combinations for an invariant combination of dummy

indexes (Corollary 3.5), there can be several different invariant expressions.

Problem 4.12: Contraction invariance

Prove that Ai Bi

Ai Bi and both are invariant, while Ai Bi is not.
Proof
Using the operation of rising/lowering indexes (42), (43), we have

Ai Bi

gi j A j gik Bk

gi j gik A j Bk

δ jk A j Bk

A jB j

which proves that both forms have the same values. If we now consider the first
form then:

∂x̄i j ∂ xk
Āi B̄i A Bk δ jk A j Bk A jB j Ai Bi
∂ x j ∂x̄i

43
which proves the point.
Consider now Ai Bi :

∂ x j ∂ xk
Āi B̄i Aj Bk
∂x̄i ∂x̄i

which can not be reduced further and, therefore is not invariant, since it has a
different form from the LHS.

44
References
[1] Barry Spain. Tensor Calculus. Oliver and Boyd, 1965.

[2] J.L. Synge and A. Schild. Tensor Calculus. Dover Publications, 1969.

[3] P. Morse and H. Feshbach. Methods of Theoretical Physics. McGraw-Hill,

New York, 1953.

45
Index
Associate tensor, 17 Order of a tensor, 5, 9
Orthogonal coordinate system, 24
Cartesian Tensors, 7 Orthogonal coordinates, 21
Christoffel’s symbol, 19
Conjugate metric tensor, 16 Permutation tensor, 11
Conjugate tensor, 17, 21 Physical component, 24
Contraction of indexes, 10
Contraction operation, 18 Raising indices, 17
Contravariant index, 5 Rank of a tensor derivative, 14
Contravariant tensor, 20 Rank of a term, 9
Contravariant vector, 5 Renaming indexes, 11
Contravariant vectors, 4 Renaming of dummy indexes, 10
Coordinate system, 4 Scalar product, 10, 18
Covariant differentiation, 19 Skew-symmetric tensor, 12
Covariant index, 5 Spatial derivative of a tensor, 14
Cross product, 12 Stretching factors, 22, 25
Cylindrical coordinates, 25
Tensor, 5
Divergence operator, 15, 21, 25 Tensor derivative, 14
Dummy index restriction, 10 Tensor equality, 8
Dummy indexes, 9 Tensor expression, 8
Free indexes, 8 Tensor identity, 13
Fundamental tensor, 16 Tensor notation, 4, 7
Tensor rules, 7
Gradient, 14 Tensor terms, 8
Time derivative of a tensor, 13
Invariance, 4, 15, 16 Transformation matrix, 5, 7
Invariant, 17–19 Transformation rule, 4
Invariant forms, 15, 18
Vector product, 12
Kronecker delta tensor, 6

Laplacian, 15, 21, 25

Lowering indices, 17

Metric tensor, 16, 18

Momentum equation, 18

Nabla, 14, 24, 25

46