Curvilinear coordinates

Short description: Coordinate system whose directions vary in space

"Lamé coefficients" redirects here. It is not to be confused with Lamé parameters (solid mechanics).

Curvilinear (top), affine (right), and Cartesian (left) coordinates in two-dimensional space

In geometry, curvilinear coordinates are a coordinate system for Euclidean space in which the coordinate lines may be curved. These coordinates may be derived from a set of Cartesian coordinates by using a transformation that is locally invertible (a one-to-one map) at each point. This means that one can convert a point given in a Cartesian coordinate system to its curvilinear coordinates and back. The name curvilinear coordinates, coined by the French mathematician Lamé, derives from the fact that the coordinate surfaces of the curvilinear systems are curved.

Well-known examples of curvilinear coordinate systems in three-dimensional Euclidean space (R3) are cylindrical and spherical coordinates. A Cartesian coordinate surface in this space is a coordinate plane; for example z = 0 defines the x-y plane. In the same space, the coordinate surface r = 1 in spherical coordinates is the surface of a unit sphere, which is curved. The formalism of curvilinear coordinates provides a unified and general description of the standard coordinate systems.

Curvilinear coordinates are often used to define the location or distribution of physical quantities which may be, for example, scalars, vectors, or tensors. Mathematical expressions involving these quantities in vector calculus and tensor analysis (such as the gradient, divergence, curl, and Laplacian) can be transformed from one coordinate system to another, according to transformation rules for scalars, vectors, and tensors. Such expressions then become valid for any curvilinear coordinate system.

A curvilinear coordinate system may be simpler to use than the Cartesian coordinate system for some applications. The motion of particles under the influence of central forces is usually easier to solve in spherical coordinates than in Cartesian coordinates; this is true of many physical problems with spherical symmetry defined in R3. Equations with boundary conditions that follow coordinate surfaces for a particular curvilinear coordinate system may be easier to solve in that system. While one might describe the motion of a particle in a rectangular box using Cartesian coordinates, it is easier to describe the motion in a sphere with spherical coordinates. Spherical coordinates are the most common curvilinear coordinate systems and are used in Earth sciences, cartography, quantum mechanics, relativity, and engineering.

## Contents

* 1 Orthogonal curvilinear coordinates in 3 dimensions
* 1.1 Coordinates, basis, and vectors
* 2 Vector calculus
* 2.1 Differential elements
* 3 Covariant and contravariant bases
* 4 Integration
* 4.1 Constructing a covariant basis in one dimension
* 4.2 Constructing a covariant basis in three dimensions
* 4.3 Jacobian of the transformation
* 5 Generalization to n dimensions
* 6 Transformation of coordinates
* 7 Vector and tensor algebra in three-dimensional curvilinear coordinates
* 8 Tensors in curvilinear coordinates
* 8.1 The metric tensor in orthogonal curvilinear coordinates
* 8.1.1 Relation to Lamé coefficients
* 8.1.2 Example: Polar coordinates
* 8.2 The alternating tensor
* 8.3 Christoffel symbols
* 8.4 Vector operations
* 9 Vector and tensor calculus in three-dimensional curvilinear coordinates
* 9.1 Geometric elements
* 9.2 Integration
* 9.3 Differentiation
* 10 Fictitious forces in general curvilinear coordinates
* 11 See also
* 12 References
* 13 Further reading
* 14 External links

## Orthogonal curvilinear coordinates in 3 dimensions

### Coordinates, basis, and vectors

Fig. 1 - Coordinate surfaces, coordinate lines, and coordinate axes of general curvilinear coordinates. Fig. 2 - Coordinate surfaces, coordinate lines, and coordinate axes of spherical coordinates. Surfaces: r \- spheres, θ - cones, Φ \- half-planes; Lines: r \- straight beams, θ - vertical semicircles, Φ \- horizontal circles; Axes: r \- straight beams, θ - tangents to vertical semicircles, Φ \- tangents to horizontal circles

For now, consider 3-D space. A point P in 3-D space (or its position vector r) can be defined using Cartesian coordinates (x, y, z) [equivalently written (x1, x2, x3)], by 𝐫=x𝐞x+y𝐞y+z𝐞z, where ex, ey, ez are the standard basis vectors.

It can also be defined by its curvilinear coordinates (q1, q2, q3) if this triplet of numbers defines a single point in an unambiguous way. The relation between the coordinates is then given by the invertible transformation functions:

x=f1(q1,q2,q3),y=f2(q1,q2,q3),z=f3(q1,q2,q3)
q1=g1(x,y,z),q2=g2(x,y,z),q3=g3(x,y,z)

The surfaces q1 = constant, q2 = constant, q3 = constant are called the coordinate surfaces; and the space curves formed by their intersection in pairs are called the coordinate curves. The coordinate axes are determined by the tangents to the coordinate curves at the intersection of three surfaces. They are not in general fixed directions in space, which happens to be the case for simple Cartesian coordinates, and thus there is generally no natural global basis for curvilinear coordinates.

In the Cartesian system, the standard basis vectors can be derived from the derivative of the location of point P with respect to the local coordinate

𝐞x=∂𝐫∂x;𝐞y=∂𝐫∂y;𝐞z=∂𝐫∂z.

Applying the same derivatives to the curvilinear system locally at point P defines the natural basis vectors:

𝐡1=∂𝐫∂q1;𝐡2=∂𝐫∂q2;𝐡3=∂𝐫∂q3.

Such a basis, whose vectors change their direction and/or magnitude from point to point is called a local basis. All bases associated with curvilinear coordinates are necessarily local. Basis vectors that are the same at all points are global bases, and can be associated only with linear or affine coordinate systems.

For this article e is reserved for the standard basis (Cartesian) and h or b is for the curvilinear basis.

These may not have unit length, and may also not be orthogonal. In the case that they are orthogonal at all points where the derivatives are well-defined, we define the Lamé coefficients (after Gabriel Lamé) by

h1=|𝐡1|;h2=|𝐡2|;h3=|𝐡3|

and the curvilinear orthonormal basis vectors by

𝐛1=𝐡1h1;𝐛2=𝐡2h2;𝐛3=𝐡3h3.

These basis vectors may well depend upon the position of P; it is therefore necessary that they are not assumed to be constant over a region. (They technically form a basis for the tangent space of ℝ3 at P, and so are local to P.)

In general, curvilinear coordinates allow the natural basis vectors hi not all mutually perpendicular to each other, and not required to be of unit length: they can be of arbitrary magnitude and direction. The use of an orthogonal basis makes vector manipulations simpler than for non-orthogonal. However, some areas of physics and engineering, particularly fluid mechanics and continuum mechanics, require non-orthogonal bases to describe deformations and fluid transport to account for complicated directional dependences of physical quantities. A discussion of the general case appears later on this page.

## Vector calculus

See also: Differential geometry

### Differential elements

In orthogonal curvilinear coordinates, since the total differential change in r is

d𝐫=∂𝐫∂q1dq1+∂𝐫∂q2dq2+∂𝐫∂q3dq3=h1dq1𝐛1+h2dq2𝐛2+h3dq3𝐛3

so scale factors are hi=|∂𝐫∂qi|

In non-orthogonal coordinates the length of d𝐫=dq1𝐡1+dq2𝐡2+dq3𝐡3 is the positive square root of d𝐫⋅d𝐫=dqidqj𝐡i⋅𝐡j (with Einstein summation convention). The six independent scalar products gij=hi.hj of the natural basis vectors generalize the three scale factors defined above for orthogonal coordinates. The nine gij are the components of the metric tensor, which has only three non zero components in orthogonal coordinates: g11=h1h1, g22=h2h2, g33=h3h3.

## Covariant and contravariant bases

Main articles: Covariance and contravariance of vectors and Raising and lowering indices

A vector v (red) represented by • a vector basis (yellow, left: e1, e2, e3), tangent vectors to coordinate curves (black) and • a covector basis or cobasis (blue, right: e1, e2, e3), normal vectors to coordinate surfaces (grey) in general (not necessarily orthogonal) curvilinear coordinates (q1, q2, q3). The basis and cobasis do not coincide unless the coordinate system is orthogonal.[1]

Spatial gradients, distances, time derivatives and scale factors are interrelated within a coordinate system by two groups of basis vectors:

1. basis vectors that are locally tangent to their associated coordinate pathline: 𝐛i=∂𝐫∂qi are contravariant vectors (denoted by lowered indices), and
2. basis vectors that are locally normal to the isosurface created by the other coordinates: 𝐛i=∇qi are covariant vectors (denoted by raised indices), ∇ is the del operator.

Note that, because of Einstein's summation convention, the position of the indices of the vectors is the opposite of that of the coordinates.

Consequently, a general curvilinear coordinate system has two sets of basis vectors for every point: {b1, b2, b3} is the contravariant basis, and {b1, b2, b3} is the covariant (a.k.a. reciprocal) basis. The covariant and contravariant basis vectors types have identical direction for orthogonal curvilinear coordinate systems, but as usual have inverted units with respect to each other.

Note the following important equality: 𝐛i⋅𝐛j=δji wherein δji denotes the generalized Kronecker delta.

Proof

In the Cartesian coordinate system (𝐞x,𝐞y,𝐞z), we can write the dot product as:

𝐛i⋅𝐛j=(∂x∂qi,∂y∂qi,∂z∂qi)⋅(∂qj∂x,∂qj∂y,∂qj∂z)=∂x∂qi∂qj∂x+∂y∂qi∂qj∂y+∂z∂qi∂qj∂z

Consider an infinitesimal displacement d𝐫=dx⋅𝐞x+dy⋅𝐞y+dz⋅𝐞z. Let dq1, dq2 and dq3 denote the corresponding infinitesimal changes in curvilinear coordinates q1, q2 and q3 respectively.

By the chain rule, dq1 can be expressed as:

dq1=∂q1∂xdx+∂q1∂ydy+∂q1∂zdz=∂q1∂xdx+∂q1∂y(∂y∂q1dq1+∂y∂q2dq2+∂y∂q3dq3)+∂q1∂z(∂z∂q1dq1+∂z∂q2dq2+∂z∂q3dq3)

If the displacement dr is such that dq2 = dq3 = 0, i.e. the position vector r moves by an infinitesimal amount along the coordinate axis q2=const and q3=const, then:

dq1=∂q1∂xdx+∂q1∂y∂y∂q1dq1+∂q1∂z∂z∂q1dq1

Dividing by dq1, and taking the limit dq1 → 0:

1=∂q1∂x∂x∂q1+∂q1∂y∂y∂q1+∂q1∂z∂z∂q1=∂x∂q1∂q1∂x+∂y∂q1∂q1∂y+∂z∂q1∂q1∂z

or equivalently:

𝐛1⋅𝐛1=1

Now if the displacement dr is such that dq1=dq3=0, i.e. the position vector r moves by an infinitesimal amount along the coordinate axis q1=const and q3=const, then:

0=∂q1∂xdx+∂q1∂y∂y∂q2dq2+∂q1∂z∂z∂q2dq2

Dividing by dq2, and taking the limit dq2 → 0:

0=∂q1∂x∂x∂q2+∂q1∂y∂y∂q2+∂q1∂z∂z∂q2=∂x∂q2∂q1∂x+∂y∂q2∂q1∂y+∂z∂q2∂q1∂z

or equivalently:

𝐛2⋅𝐛1=0

And so forth for the other dot products.

Alternative Proof:

δjidqj=dqi=∇qi⋅d𝐫=𝐛i⋅∂𝐫∂qjdqj=𝐛i⋅𝐛jdqj

and the Einstein summation convention is implied.

A vector v can be specified in terms of either basis, i.e.,

𝐯=v1𝐛1+v2𝐛2+v3𝐛3=v1𝐛1+v2𝐛2+v3𝐛3

Using the Einstein summation convention, the basis vectors relate to the components by[2]: 30–32

𝐯⋅𝐛i=vk𝐛k⋅𝐛i=vkδki=vi
𝐯⋅𝐛i=vk𝐛k⋅𝐛i=vkδik=vi

and

𝐯⋅𝐛i=vk𝐛k⋅𝐛i=gkivk
𝐯⋅𝐛i=vk𝐛k⋅𝐛i=gkivk

where g is the metric tensor (see below).

A vector can be specified with covariant coordinates (lowered indices, written vk) or contravariant coordinates (raised indices, written vk). From the above vector sums, it can be seen that contravariant coordinates are associated with covariant basis vectors, and covariant coordinates are associated with contravariant basis vectors.

A key feature of the representation of vectors and tensors in terms of indexed components and basis vectors is invariance in the sense that vector components which transform in a covariant manner (or contravariant manner) are paired with basis vectors that transform in a contravariant manner (or covariant manner).

## Integration

Main article: Covariant transformation

### Constructing a covariant basis in one dimension

Fig. 3 – Transformation of local covariant basis in the case of general curvilinear coordinates

Consider the one-dimensional curve shown in Fig. 3. At point P, taken as an origin, x is one of the Cartesian coordinates, and q1 is one of the curvilinear coordinates. The local (non-unit) basis vector is b1 (notated h1 above, with b reserved for unit vectors) and it is built on the q1 axis which is a tangent to that coordinate line at the point P. The axis q1 and thus the vector b1 form an angle α with the Cartesian x axis and the Cartesian basis vector e1.

It can be seen from triangle PAB that

cos⁡α=|𝐞1||𝐛1|⇒|𝐞1|=|𝐛1|cos⁡α

where |e1|, |b1| are the magnitudes of the two basis vectors, i.e., the scalar intercepts PB and PA. PA is also the projection of b1 on the x axis.

However, this method for basis vector transformations using directional cosines is inapplicable to curvilinear coordinates for the following reasons:

1. By increasing the distance from P, the angle between the curved line q1 and Cartesian axis x increasingly deviates from α.
2. At the distance PB the true angle is that which the tangent at point C forms with the x axis and the latter angle is clearly different from α.

The angles that the q1 line and that axis form with the x axis become closer in value the closer one moves towards point P and become exactly equal at P.

Let point E be located very close to P, so close that the distance PE is infinitesimally small. Then PE measured on the q1 axis almost coincides with PE measured on the q1 line. At the same time, the ratio PD/PE (PD being the projection of PE on the x axis) becomes almost exactly equal to cos⁡α.

Let the infinitesimally small intercepts PD and PE be labelled, respectively, as dx and dq1. Then

cos⁡α=dxdq1=|𝐞1||𝐛1|.

Thus, the directional cosines can be substituted in transformations with the more exact ratios between infinitesimally small coordinate intercepts. It follows that the component (projection) of b1 on the x axis is

p1=𝐛1⋅𝐞1|𝐞1|=|𝐛1||𝐞1||𝐞1|cos⁡α=|𝐛1|dxdq1⇒p1|𝐛1|=dxdq1.

If qi = qi(x1, x2, x3) and xi = xi(q1, q2, q3) are smooth (continuously differentiable) functions the transformation ratios can be written as ∂qi∂xj and ∂xi∂qj. That is, those ratios are partial derivatives of coordinates belonging to one system with respect to coordinates belonging to the other system.

### Constructing a covariant basis in three dimensions

Doing the same for the coordinates in the other 2 dimensions, b1 can be expressed as:

𝐛1=p1𝐞1+p2𝐞2+p3𝐞3=∂x1∂q1𝐞1+∂x2∂q1𝐞2+∂x3∂q1𝐞3

Similar equations hold for b2 and b3 so that the standard basis {e1, e2, e3} is transformed to a local (ordered and normalised) basis {b1, b2, b3} by the following system of equations:

𝐛1=∂x1∂q1𝐞1+∂x2∂q1𝐞2+∂x3∂q1𝐞3𝐛2=∂x1∂q2𝐞1+∂x2∂q2𝐞2+∂x3∂q2𝐞3𝐛3=∂x1∂q3𝐞1+∂x2∂q3𝐞2+∂x3∂q3𝐞3

By analogous reasoning, one can obtain the inverse transformation from local basis to standard basis:

𝐞1=∂q1∂x1𝐛1+∂q2∂x1𝐛2+∂q3∂x1𝐛3𝐞2=∂q1∂x2𝐛1+∂q2∂x2𝐛2+∂q3∂x2𝐛3𝐞3=∂q1∂x3𝐛1+∂q2∂x3𝐛2+∂q3∂x3𝐛3

### Jacobian of the transformation

The above systems of linear equations can be written in matrix form using the Einstein summation convention as

∂xi∂qk𝐞i=𝐛k,∂qi∂xk𝐛i=𝐞k.

This coefficient matrix of the linear system is the Jacobian matrix (and its inverse) of the transformation. These are the equations that can be used to transform a Cartesian basis into a curvilinear basis, and vice versa.

In three dimensions, the expanded forms of these matrices are

𝐉=[∂x1∂q1∂x1∂q2∂x1∂q3∂x2∂q1∂x2∂q2∂x2∂q3∂x3∂q1∂x3∂q2∂x3∂q3],𝐉−1=[∂q1∂x1∂q1∂x2∂q1∂x3∂q2∂x1∂q2∂x2∂q2∂x3∂q3∂x1∂q3∂x2∂q3∂x3]

In the inverse transformation (second equation system), the unknowns are the curvilinear basis vectors. For any specific location there can only exist one and only one set of basis vectors (else the basis is not well defined at that point). This condition is satisfied if and only if the equation system has a single solution. In linear algebra, a linear equation system has a single solution (non-trivial) only if the determinant of its system matrix is non-zero:

det⁡(𝐉−1)≠0

which shows the rationale behind the above requirement concerning the inverse Jacobian determinant.

## Generalization to n dimensions

The formalism extends to any finite dimension as follows.

Consider the real Euclidean n-dimensional space, that is Rn = R × R × ... × R (n times) where R is the set of real numbers and × denotes the Cartesian product, which is a vector space.

The coordinates of this space can be denoted by: x = (x1, x2,...,xn). Since this is a vector (an element of the vector space), it can be written as:

𝐱=∑i=1nxi𝐞i

where e1 = (1,0,0...,0), e2 = (0,1,0...,0), e3 = (0,0,1...,0),...,en = (0,0,0...,1) is the standard basis set of vectors for the space Rn, and i = 1, 2,...n is an index labelling components. Each vector has exactly one component in each dimension (or "axis") and they are mutually orthogonal (perpendicular) and normalized (has unit magnitude).

More generally, we can define basis vectors bi so that they depend on q = (q1, q2,...,qn), i.e. they change from point to point: bi = bi(q). In which case to define the same point x in terms of this alternative basis: the coordinates with respect to this basis vi also necessarily depend on x also, that is vi = vi(x). Then a vector v in this space, with respect to these alternative coordinates and basis vectors, can be expanded as a linear combination in this basis (which simply means to multiply each basis vector ei by a number vi – scalar multiplication):

𝐯=∑j=1nv¯j𝐛j=∑j=1nv¯j(𝐪)𝐛j(𝐪)

The vector sum that describes v in the new basis is composed of different vectors, although the sum itself remains the same.

## Transformation of coordinates

From a more general and abstract perspective, a curvilinear coordinate system is simply a coordinate patch on the differentiable manifold En (n-dimensional Euclidean space) that is diffeomorphic to the Cartesian coordinate patch on the manifold.[3] Two diffeomorphic coordinate patches on a differential manifold need not overlap differentiably. With this simple definition of a curvilinear coordinate system, all the results that follow below are simply applications of standard theorems in differential topology.

The transformation functions are such that there's a one-to-one relationship between points in the "old" and "new" coordinates, that is, those functions are bijections, and fulfil the following requirements within their domains:

1. They are smooth functions: qi = qi(x)
2. The inverse Jacobian determinant

J−1=|∂q1∂x1∂q1∂x2⋯∂q1∂xn∂q2∂x1∂q2∂x2⋯∂q2∂xn⋮⋮⋱⋮∂qn∂x1∂qn∂x2⋯∂qn∂xn|≠0
is not zero; meaning the transformation is invertible: xi(q) according to the inverse function theorem. The condition that the Jacobian determinant is not zero reflects the fact that three surfaces from different families intersect in one and only one point and thus determine the position of this point in a unique way.[4]

## Vector and tensor algebra in three-dimensional curvilinear coordinates

Note: the Einstein summation convention of summing on repeated indices is used below.

Elementary vector and tensor algebra in curvilinear coordinates is used in some of the older scientific literature in mechanics and physics and can be indispensable to understanding work from the early and mid-1900s, for example the text by Green and Zerna.[5] Some useful relations in the algebra of vectors and second-order tensors in curvilinear coordinates are given in this section. The notation and contents are primarily from Ogden,[6] Naghdi,[7] Simmonds,[2] Green and Zerna,[5] Basar and Weichert,[8] and Ciarlet.[9]

## Tensors in curvilinear coordinates

Main article: Tensors in curvilinear coordinates

A second-order tensor can be expressed as

S=Sij𝐛i⊗𝐛j=Sij𝐛i⊗𝐛j=Sij𝐛i⊗𝐛j=Sij𝐛i⊗𝐛j

where ⊗ denotes the tensor product. The components Sij are called the contravariant components, Si j the mixed right-covariant components, Si j the mixed left-covariant components, and Sij the covariant components of the second-order tensor. The components of the second-order tensor are related by

Sij=gikSkj=gjkSik=gikgjℓSkℓ

### The metric tensor in orthogonal curvilinear coordinates

Main article: Metric tensor

At each point, one can construct a small line element dx, so the square of the length of the line element is the scalar product dx • dx and is called the metric of the space, given by:

d𝐱⋅d𝐱=∂xi∂qj∂xi∂qkdqjdqk.

The following portion of the above equation

∂xk∂qi∂xk∂qj=gij(qi,qj)=𝐛i⋅𝐛j

is a symmetric tensor called the fundamental (or metric) tensor of the Euclidean space in curvilinear coordinates.

Indices can be raised and lowered by the metric:

vi=gikvk

#### Relation to Lamé coefficients

Defining the scale factors hi by

hihj=gij=𝐛i⋅𝐛j⇒hi=gii=|𝐛i|=|∂𝐱∂qi|

gives a relation between the metric tensor and the Lamé coefficients, and

gij=∂𝐱∂qi⋅∂𝐱∂qj=(hki𝐞k)⋅(hmj𝐞m)=hkihkj

where hij are the Lamé coefficients. For an orthogonal basis we also have:

g=g11g22g33=h12h22h32⇒g=h1h2h3=J

#### Example: Polar coordinates

If we consider polar coordinates for R2,

(x,y)=(rcos⁡θ,rsin⁡θ)

(r, θ) are the curvilinear coordinates, and the Jacobian determinant of the transformation (r,θ) → (r cos θ, r sin θ) is r.

The orthogonal basis vectors are br = (cos θ, sin θ), bθ = (−r sin θ, r cos θ). The scale factors are hr = 1 and hθ= r. The fundamental tensor is g11 =1, g22 =r2, g12 = g21 =0.

### The alternating tensor

In an orthonormal right-handed basis, the third-order alternating tensor is defined as

ℰ=εijk𝐞i⊗𝐞j⊗𝐞k

In a general curvilinear basis the same tensor may be expressed as

ℰ=ℰijk𝐛i⊗𝐛j⊗𝐛k=ℰijk𝐛i⊗𝐛j⊗𝐛k

It can also be shown that

ℰijk=1Jεijk=1+gεijk

### Christoffel symbols

Christoffel symbols of the first kind Γkij

𝐛i,j=∂𝐛i∂qj=𝐛kΓkij⇒𝐛k⋅𝐛i,j=Γkij

where the comma denotes a partial derivative (see Ricci calculus). To express Γkij in terms of gij,

gij,k=(𝐛i⋅𝐛j),k=𝐛i,k⋅𝐛j+𝐛i⋅𝐛j,k=Γjik+Γijkgik,j=(𝐛i⋅𝐛k),j=𝐛i,j⋅𝐛k+𝐛i⋅𝐛k,j=Γkij+Γikjgjk,i=(𝐛j⋅𝐛k),i=𝐛j,i⋅𝐛k+𝐛j⋅𝐛k,i=Γkji+Γjki

Since

𝐛i,j=𝐛j,i⇒Γkij=Γkji

using these to rearrange the above relations gives

Γkij=12(gik,j+gjk,i−gij,k)=12[(𝐛i⋅𝐛k),j+(𝐛j⋅𝐛k),i−(𝐛i⋅𝐛j),k]

Christoffel symbols of the second kind Γkji

Γkij=gklΓlij=Γkji,∂𝐛i∂qj=𝐛kΓkij

This implies that

Γkij=∂𝐛i∂qj⋅𝐛k=−𝐛i⋅∂𝐛k∂qj since ∂∂qj(𝐛i⋅𝐛k)=0.

Other relations that follow are

∂𝐛i∂qj=−Γijk𝐛k,∇𝐛i=Γkij𝐛k⊗𝐛j,∇𝐛i=−Γijk𝐛k⊗𝐛j

### Vector operations

1. Dot product:

The scalar product of two vectors in curvilinear coordinates is[2]: 32

𝐮⋅𝐯=uivi=uivi=gijuivj=gijuivj
2. Cross product:

The cross product of two vectors is given by[2]: 32–34

𝐮×𝐯=ϵijkujvk𝐞i

where ϵijk is the permutation symbol and 𝐞i is a Cartesian basis vector. In curvilinear coordinates, the equivalent expression is

𝐮×𝐯=[(𝐛m×𝐛n)⋅𝐛s]umvn𝐛s=ℰsmnumvn𝐛s
where ℰijk is the third-order alternating tensor.

## Vector and tensor calculus in three-dimensional curvilinear coordinates

Note: the Einstein summation convention of summing on repeated indices is used below.

Adjustments need to be made in the calculation of line, surface and volume integrals. For simplicity, the following restricts to three dimensions and orthogonal curvilinear coordinates. However, the same arguments apply for n-dimensional spaces. When the coordinate system is not orthogonal, there are some additional terms in the expressions.

Simmonds,[2] in his book on tensor analysis, quotes Albert Einstein saying[10]

> The magic of this theory will hardly fail to impose itself on anybody who has truly understood it; it represents a genuine triumph of the method of absolute differential calculus, founded by Gauss, Riemann, Ricci, and Levi-Civita.

Vector and tensor calculus in general curvilinear coordinates is used in tensor analysis on four-dimensional curvilinear manifolds in general relativity,[11] in the mechanics of curved shells,[9] in examining the invariance properties of Maxwell's equations which has been of interest in metamaterials[12][13] and in many other fields.

Some useful relations in the calculus of vectors and second-order tensors in curvilinear coordinates are given in this section. The notation and contents are primarily from Ogden,[14] Simmonds,[2] Green and Zerna,[5] Basar and Weichert,[8] and Ciarlet.[9]

Let φ = φ(x) be a well defined scalar field and v = v(x) a well-defined vector field, and λ1, λ2... be parameters of the coordinates

### Geometric elements

1. Tangent vector: If x(λ) parametrizes a curve C in Cartesian coordinates, then

∂𝐱∂λ=∂𝐱∂qi∂qi∂λ=(hki∂qi∂λ)𝐛k

is a tangent vector to C in curvilinear coordinates (using the chain rule). Using the definition of the Lamé coefficients, and that for the metric gij = 0 when i ≠ j, the magnitude is:

|∂𝐱∂λ|=hkihkj∂qi∂λ∂qj∂λ=gij∂qi∂λ∂qj∂λ=hi2(∂qi∂λ)2
2. Tangent plane element: If x(λ1, λ2) parametrizes a surface S in Cartesian coordinates, then the following cross product of tangent vectors is a normal vector to S with the magnitude of infinitesimal plane element, in curvilinear coordinates. Using the above result,

∂𝐱∂λ1×∂𝐱∂λ2=(∂𝐱∂qi∂qi∂λ1)×(∂𝐱∂qj∂qj∂λ2)=ℰkmp(hki∂qi∂λ1)(hmj∂qj∂λ2)𝐛p

where ℰ is the permutation symbol. In determinant form:

∂𝐱∂λ1×∂𝐱∂λ2=|𝐞1𝐞2𝐞3h1i∂qi∂λ1h2i∂qi∂λ1h3i∂qi∂λ1h1j∂qj∂λ2h2j∂qj∂λ2h3j∂qj∂λ2|

### Integration

Operator  | Scalar field  | Vector field

|
|

Line integral | ∫Cφ(𝐱)ds=∫abφ(𝐱(λ))|∂𝐱∂λ|dλ | ∫C𝐯(𝐱)⋅d𝐬=∫ab𝐯(𝐱(λ))⋅(∂𝐱∂λ)dλ
Surface integral | ∫Sφ(𝐱)dS=∬Tφ(𝐱(λ1,λ2))|∂𝐱∂λ1×∂𝐱∂λ2|dλ1dλ2 | ∫S𝐯(𝐱)⋅dS=∬T𝐯(𝐱(λ1,λ2))⋅(∂𝐱∂λ1×∂𝐱∂λ2)dλ1dλ2
Volume integral | ∭Vφ(x,y,z)dV=∭Vχ(q1,q2,q3)Jdq1dq2dq3 | ∭V𝐮(x,y,z)dV=∭V𝐯(q1,q2,q3)Jdq1dq2dq3

### Differentiation

The expressions for the gradient, divergence, and Laplacian can be directly extended to n-dimensions, however the curl is only defined in 3D.

The vector field bi is tangent to the qi coordinate curve and forms a natural basis at each point on the curve. This basis, as discussed at the beginning of this article, is also called the covariant curvilinear basis. We can also define a reciprocal basis, or contravariant curvilinear basis, bi. All the algebraic relations between the basis vectors, as discussed in the section on tensor algebra, apply for the natural basis and its reciprocal at each point x.

Operator  | Scalar field  | Vector field  | 2nd order tensor field

|
|
|

Gradient | ∇φ=1hi∂φ∂qi𝐛i | ∇𝐯=1hi2∂𝐯∂qi⊗𝐛i | ∇S=∂S∂qi⊗𝐛i
Divergence | N/A  | ∇⋅𝐯=1∏jhj∂∂qi(vi∏j≠ihj) | (∇⋅S)⋅𝐚=∇⋅(S⋅𝐚)

where a is an arbitrary constant vector. In curvilinear coordinates,

∇⋅S=[∂Sij∂qk−ΓkilSlj−ΓkjlSil]gik𝐛j

Laplacian  | ∇2φ=1∏jhj∂∂qi(∏jhjhi2∂φ∂qi) |

∇2𝐯≡∇∇⋅𝐯−∇×∇×𝐯 =𝐱^∇2vx+𝐲^∇2vy+𝐳^∇2vz (First equality in 3D only; second equality in Cartesian components only)

Curl | N/A  | For vector fields in 3D only,

∇×𝐯=1h1h2h3𝐞iϵijkhi∂(hkvk)∂qj

where ϵijk is the Levi-Civita symbol.

| See Curl of a tensor field

## Fictitious forces in general curvilinear coordinates

By definition, if a particle with no forces acting on it has its position expressed in an inertial coordinate system, (x1, x2, x3, t), then there it will have no acceleration (d2xj/dt2 = 0).[15] In this context, a coordinate system can fail to be "inertial" either due to non-straight time axis or non-straight space axes (or both). In other words, the basis vectors of the coordinates may vary in time at fixed positions, or they may vary with position at fixed times, or both. When equations of motion are expressed in terms of any non-inertial coordinate system (in this sense), extra terms appear, called Christoffel symbols. Strictly speaking, these terms represent components of the absolute acceleration (in classical mechanics), but we may also choose to continue to regard d2xj/dt2 as the acceleration (as if the coordinates were inertial) and treat the extra terms as if they were forces, in which case they are called fictitious forces.[16] The component of any such fictitious force normal to the path of the particle and in the plane of the path's curvature is then called centrifugal force.[17]

This more general context makes clear the correspondence between the concepts of centrifugal force in rotating coordinate systems and in stationary curvilinear coordinate systems. (Both of these concepts appear frequently in the literature.[18][19][20]) For a simple example, consider a particle of mass m moving in a circle of radius r with angular speed w relative to a system of polar coordinates rotating with angular speed W. The radial equation of motion is mr” = Fr + mr(w + W)2. Thus the centrifugal force is mr times the square of the absolute rotational speed A = w + W of the particle. If we choose a coordinate system rotating at the speed of the particle, then W = A and w = 0, in which case the centrifugal force is mrA2, whereas if we choose a stationary coordinate system we have W = 0 and w = A, in which case the centrifugal force is again mrA2. The reason for this equality of results is that in both cases the basis vectors at the particle's location are changing in time in exactly the same way. Hence these are really just two different ways of describing exactly the same thing, one description being in terms of rotating coordinates and the other being in terms of stationary curvilinear coordinates, both of which are non-inertial according to the more abstract meaning of that term.

When describing general motion, the actual forces acting on a particle are often referred to the instantaneous osculating circle tangent to the path of motion, and this circle in the general case is not centered at a fixed location, and so the decomposition into centrifugal and Coriolis components is constantly changing. This is true regardless of whether the motion is described in terms of stationary or rotating coordinates.

## See also

* Curvilinear motion
* Covariance and contravariance
* Introduction to the mathematics of general relativity
* Special cases:
* Orthogonal coordinates
* Skew coordinates
* Tensors in curvilinear coordinates
* Frenet–Serret formulas
* Covariant derivative
* Tensor derivative (continuum mechanics)
* Curvilinear perspective
* Del in cylindrical and spherical coordinates

## References

1. ↑ J.A. Wheeler; C. Misner; K.S. Thorne (1973). Gravitation. W.H. Freeman & Co. ISBN 0-7167-0344-0.
2. ↑ 2.0 2.1 2.2 2.3 2.4 2.5 Simmonds, J. G. (1994). A brief on tensor analysis. Springer. ISBN 0-387-90639-8.
3. ↑ Boothby, W. M. (2002). An Introduction to Differential Manifolds and Riemannian Geometry (revised ed.). New York, NY: Academic Press.
4. ↑ McConnell, A. J. (1957). Application of Tensor Analysis. New York, NY: Dover Publications, Inc.. Ch. 9, sec. 1. ISBN 0-486-60373-3. https://archive.org/details/applicationoften0000mcco.
5. ↑ 5.0 5.1 5.2 Green, A. E.; Zerna, W. (1968). Theoretical Elasticity. Oxford University Press. ISBN 0-19-853486-8.
6. ↑ Ogden, R. W. (2000). Nonlinear elastic deformations. Dover.
7. ↑ Naghdi, P. M. (1972). "Theory of shells and plates". in S. Flügge. Handbook of Physics. VIa/2. pp. 425–640.
8. ↑ 8.0 8.1 Basar, Y.; Weichert, D. (2000). Numerical continuum mechanics of solids: fundamental concepts and perspectives. Springer.
9. ↑ 9.0 9.1 9.2 Ciarlet, P. G. (2000). Theory of Shells. 1. Elsevier Science.
10. ↑ Einstein, A. (1915). "Contribution to the Theory of General Relativity". in Laczos, C.. The Einstein Decade. p. 213. ISBN 0-521-38105-3.
11. ↑ Misner, C. W.; Thorne, K. S.; Wheeler, J. A. (1973). Gravitation. W. H. Freeman and Co.. ISBN 0-7167-0344-0.
12. ↑ Greenleaf, A.; Lassas, M.; Uhlmann, G. (2003). "Anisotropic conductivities that cannot be detected by EIT". Physiological Measurement 24 (2): 413–419. doi:10.1088/0967-3334/24/2/353. PMID 12812426.
13. ↑ Leonhardt, U.; Philbin, T.G. (2006). "General relativity in electrical engineering". New Journal of Physics 8 (10): 247. doi:10.1088/1367-2630/8/10/247. Bibcode: 2006NJPh....8..247L.
14. ↑ Ogden
15. ↑ Friedman, Michael (1989). The Foundations of Space–Time Theories. Princeton University Press. ISBN 0-691-07239-6.
16. ↑ Stommel, Henry M.; Moore, Dennis W. (1989). An Introduction to the Coriolis Force. Columbia University Press. ISBN 0-231-06636-8. https://archive.org/details/introductiontoco0000stom.
17. ↑ Beer; Johnston (1972). Statics and Dynamics (2nd ed.). McGraw–Hill. p. 485. ISBN 0-07-736650-6.
18. ↑ Hildebrand, Francis B. (1992). Methods of Applied Mathematics. Dover. p. 156. ISBN 0-13-579201-0. https://archive.org/details/methodsofapplied00hild/page/156.
19. ↑ McQuarrie, Donald Allan (2000). Statistical Mechanics. University Science Books. ISBN 0-06-044366-9. https://archive.org/details/statisticalmecha00mcqu_0.
20. ↑ Weber, Hans-Jurgen; Arfken, George Brown (2004). Essential Mathematical Methods for Physicists. Academic Press. p. 843. ISBN 0-12-059877-9.

## Further reading

* Spiegel, M. R. (1959). Vector Analysis. New York: Schaum's Outline Series. ISBN 0-07-084378-3.
* Arfken, George (1995). Mathematical Methods for Physicists. Academic Press. ISBN 0-12-059877-9.

## External links

* Planetmath.org Derivation of Unit vectors in curvilinear coordinates
* MathWorld's page on Curvilinear Coordinates
* Prof. R. Brannon's E-Book on Curvilinear Coordinates
* Wikiversity:Introduction to Elasticity/Tensors#The divergence of a tensor field – Wikiversity, Introduction to Elasticity/Tensors.

* v
* t
* e

Orthogonal coordinate systems

Two dimensional|

* Cartesian
* Polar
* Parabolic
* Bipolar
* Elliptic

Three dimensional|

* Cartesian coordinate system
* Cylindrical coordinate system
* Spherical coordinate system
* Parabolic cylindrical coordinates
* Paraboloidal coordinates
* Oblate spheroidal coordinates
* Prolate spheroidal coordinates
* Ellipsoidal coordinates
* Elliptic cylindrical coordinates
* Toroidal coordinates
* Bispherical coordinates
* Bipolar cylindrical coordinates
* Conical coordinates
* 6-sphere coordinates
* Flat-ring cyclide coordinates
* Flat-disk cyclide coordinates
* Bi-cyclide coordinates
* Cap-cyclide coordinates
* Concave bi-sinusoidal single-centered coordinates
* Concave bi-sinusoidal double-centered coordinates
* Convex inverted-sinusoidal spherically aligned coordinates
* Quasi-random-intersection cartesian coordinates

0.00

(0 votes)

Original source: https://en.wikipedia.org/wiki/Curvilinear coordinates. Read more

|

Retrieved from "https://handwiki.org/wiki/index.php?title=Curvilinear_coordinates&oldid=4298660"
*[v]: View this template
*[t]: Discuss this template
*[e]: Edit this template