Coordinate Systems and Vector Fields

Total Page:16

File Type:pdf, Size:1020Kb

Coordinate Systems and Vector Fields Vector calculus in curvilinear coordinates D. Jaksch Goals: • Understand the difference between coordinates and vector components • Understand the implications of basis vectors in curvilinear coordinates not being constant • Learn how to use curvilinear coordinate systems in vector calculus Coordinate systems and vector fields Coordinate systems A point in coordinate space r is often represented as r = (x; y; z)T with x, y, and z the distances along the three coordinate axes. We can equally introduce cylindrical polar coordinates which we will use here as the prime example for curvilinear coordinate systems. They are defined through the relations x = ρ cos('); y = ρ sin('); and z = z: The point in space is now written as r = (ρ cos('); ρ sin('); z)T where the meaning of the components is still that they give the distances along the three coordinate axes. Vector fields A vector field assigns a vector to each point r and is usually denoted as F(r) or simply F. The vector field is often defined through components Fi(r) which are the projections of the vector onto the three coordinate axes. For instance F = (−y; x; 0)T =px2 + y2 assigns vectors as indicated in figure 1a). Using cylindrical polar coordinates this vector field is given by F = (− sin('); cos('); 0)T . Here we have rewritten the vector field in different coordinates but not changed the meaning of its components. a) b) c) 푒̂휑 퐹푥 푒휌̂ 휑 퐹푦 푒̂푦 퐹 z y z y z y 푒̂푥 푒̂휑 퐹푥 퐹휑 퐹푦 푒̂푦 x x x 휌 푒̂푥 푒̂ Figure 1: (a) Vector field F = (−y; x; 0)T =px2 + y2 (red arrows) and showing cartesian (b) and cylindrical polar (c) components for two of its vectors. Basis vectors Each point in coordinate space has a vector space associated with it where the vectors of vector fields live. A coordinate system fxig allows us to define bases for all of these vector spaces through @r ei = : @xi ^ ^ ^ For cartesian coordinates the normalized basis vectors are e^x = i, e^y = j, and e^z = k pointing along the three coordinate axes. They are orthogonal, normalized and constant, i.e. their direction does not change with the point r1. Next we calculate basis vectors for a curvilinear coordinate systems using again cylindrical polar coordinates. They are given by @r 1 @r e^ = = cos(')e^ + sin(')e^ ; e^ = = − sin(')e^ + cos(')e^ ; and e^ = e^ : ρ @ρ x y ' ρ @' x y z z The 1/ρ in the definition of e^' is required for the vector to be properly normalized to 1. These basis vectors are mutually orthogonal and normalized. However, they are not constant, their direction changes with position. By inverting this set of linear equations for the basis vectors we find e^x = cos(')e^ρ − sin(')e^' ; e^y = sin(')e^ρ + cos(')e^' ; and e^z = e^z : T We can now write a vector field F(r) = (Fx(r);Fy(r);Fz(r))c ≡ Fx(r)e^x + Fy(r)e^y + Fz(r)e^z in new components using the above relations between the basis vectors for different coordinate systems. This is still T the same vector field but now written as F(r) = (Fρ(r);F'(r);Fz(r))p ≡ Fρ(r)e^ρ + F'(r)e^' + Fz(r)e^z. The new components are projections on the basis vectors e^ρ, e^', and e^z. We use indices c and p to make explicit which components are used. In detail, the change of components is carried out by 0 1 Fx F = @ Fy A = Fxe^x + Fye^y + Fze^z = Fx[cos(')e^ρ − sin(')e^'] F z c 0 1 Fx cos(') + Fy sin(') +Fy[sin(')e^ρ + cos(')e^'] + Fze^z = @ −Fx sin(') + Fy cos(') A F z p T T Cartesian and cylindrical polar components for the vector field F = (− sin('); cos('); 0)c = (0; 1; 0)p are shown T T in figures 1b) and 1c) respectively. Produce a similar figure for the vector field F = r = (x; y; z)c = (ρ, 0; z)p . The most commonly used coordinate systems produce orthogonal right handed bases. This means that scalar and vector products between the basis vectors obey the familiar relations e^i · e^j = δij and e^i × e^j = e^kijk, where δij is the Kronecker delta and ijk the totally antisymmetric tensor. For cylindrical polar coordinates (ρ, '; z) we explicitly have e^ρ · e^ρ = e^' · e^' = e^z · e^z = 1 and e^ρ · e^' = e^' · e^z = e^z · e^ρ = 0 ; e^ρ × e^ρ = e^' × e^' = e^z × e^z = 0; e^ρ × e^' = e^z ; e^' × e^z = e^ρ and e^z × e^ρ = e^' ; Vector Calculus Partial derivatives The partial derivatives with respect to the coordinates are found using the chain rule2 @ρ = cos(')@x + sin(')@y ;@' = −ρ sin(')@x + ρ cos(')@y ; and @z = @z ; and sin(') cos(') @ = cos(')@ − @ ;@ = sin(')@ + @ ; and @ = @ : x ρ ρ ' y ρ ρ ' z z r-operator We now rewrite the r-operator by using these relations 0 1 0 1 @x @ρ e^' 1 r = @ = e^ @ + e^ @ + e^ @ = e^ @ + @ + e^ @ = @' @ y A x x y y z z ρ ρ ρ ' z z @ ρ A @ @ z c z p 1This might seem obvious but needs to be revisited in relativity and has far reaching consequences for the physical nature of space. 2 The calculation is identical to working out the basis vectors but replacing @r=@i ! @i. The partial derivatives can thus directly be read off from the relations between the basis vectors. Gradient The gradient of a scalar field U in cylindrical polar coordinates is now given by 0 @U 1 @ρ 1 @U e^' gradU = rU = @ ρ @' A = e^ρ@ρU + @'U + e^z@zU: @U ρ @z p The expression for the r-operator in cylindrical polar components is thus indirectly given on the data-sheet3. There you will find an expression for rU and the del-operator is found by simply leaving out U in this expression. Divergence When working out the divergence we need to properly take into account that the basis vectors are not constant in general curvilinear coordinates. For cylindrical polar coordinates we have two nonzero derivatives @'e^' = − cos(')e^x − sin(')e^y = −e^ρ and @'e^ρ = − sin(')e^x + cos(')e^y = e^' : P The divergence will thus in general not be given by r · F(r) = i @iFi(r) which is only true for an orthogonal coordinate system whose basis vectors are constant in space. Using the product rule we find e^ e^ r · F = [e^ @ + ' @ + e^ @ ] · [F e^ + F e^ + F e^ ] = @ F + ' · [@ (F e^ ) + @ (F e^ )] + @ F ρ ρ ρ ' z z ρ ρ ' ' z z ρ ρ ρ ' ρ ρ ' ' ' z z e^ @ (ρF ) @ F = @ F + ' · [e^ @ F + F e^ + e^ @ F − e^ F ] + @ F = ρ ρ + ' ' + @ F : ρ ρ ρ ρ ' ρ ρ ' ' ' ' ρ ' z z ρ ρ z z Rotation The calculation is similar to working out the divergence e^ r × F = [e^ @ + ' @ + e^ @ ] × [F e^ + F e^ + F e^ ] = e^ × [e^ @ F + e^ @ F ] + ρ ρ ρ ' z z ρ ρ ' ' z z ρ ' ρ ' z ρ z e^ ' × [e^ @ F − e^ F + e^ @ F ] + e^ × [e^ @ F + e^ @ F ] ρ ρ ' ρ ρ ' z ' z z ρ z ρ ' z ' @ F @ (ρF ) @ F = e^ ' z − @ F + e^ [−@ F + @ F ] + e^ ρ ' − ' ρ : ρ ρ z ' ' ρ z z ρ z ρ ρ Laplace operator We obtain the Laplace operator by replacing F ! r in the expression for the divergence 1 1 ∆ = r · r = @ (ρ∂ ) + @2 + @2 : ρ ρ ρ ρ2 ' z Example A velocity field describing a vortex is given by u = Ae^'/ρ with A a constant. Its divergence and rotation are given by e^ Ae^ e^ Ae^ A r · u = e^ @ + ' @ + e^ @ · ' = 0 ; r × u = e^ @ + ' @ + e^ @ × ' = (e^ − e^ ) = 0 : ρ ρ ρ ' z z ρ ρ ρ ρ ' z z ρ ρ2 z z The acceleration of a parcel of fluid is given by (this term appears in the Navier Stokes equation) A Ae^ A (u · r)u = @ ' = − e^ : ρ2 ' ρ ρ3 ρ Work out these expressions using cartesian components. 3Also spherical polar coordinates can be found on the data sheet. Summary Cylindrical polar coordinates (ρ, '; z) • Relation to cartesian coordinates (care is required to obtain ' in the correct quadrant) x = ρ cos(') ; y = ρ sin(') and z = z p ρ = x2 + y2 ;' = arctan(y=x) and z = z • Basis vectors e^ρ = @ρr, ρe^' = @'r, e^z = @zr e^ρ = cos(')e^x + sin(')e^y ; e^' = − sin(')e^x + cos(')e^y and e^z = e^z e^x = cos(')e^ρ − sin(')e^' ; e^y = sin(')e^ρ + cos(')e^' and e^z = e^z • Non-zero derivatives of basis vectors @'e^ρ = e^' ;@'e^' = −e^ρ • r-operator e^ r = e^ @ + ' @ + e^ @ ρ ρ ρ ' z z • Divergence @ (ρF ) @ F r · F = ρ ρ + ' ' + @ F ρ ρ z z • Rotation @ F @ (ρF ) @ F r × F(r) = e^ ' z − @ F + e^ [−@ F + @ F ] + e^ ρ ' − ' ρ ρ ρ z ' ' ρ z z ρ z ρ ρ • Laplace 1 1 ∆ = @ (ρ∂ ) + @2 + @2 ρ ρ ρ ρ2 ' z Spherical polar coordinates (r; θ; ') • Relation to cartesian and cylindrical coordinates (care is required to obtain ' in the correct quadrant) x = r cos(') sin(θ) ; y = r sin(') sin(θ) ; z = r cos(θ) and ρ = r sin(θ) p r = x2 + y2 + z2 ;' = arctan(y=x) and θ = arccos(z=r) = arctan(ρ/z) • Basis vectors e^r = @rr, re^θ = @θr, r sin(θ)e^' = @'r e^r = cos(θ)e^z + sin(θ)e^ρ ; e^θ = − sin(θ)e^z + cos(θ)e^ρ and e^' = e^' e^z = cos(θ)e^r − sin(θ)e^θ ; e^ρ = sin(θ)e^r + cos(θ)e^θ and e^' = e^' e^x = sin(θ) cos(')e^r + cos(θ) cos(')e^θ − sin(')e^' and e^y = sin(θ) sin(')e^r + cos(θ) sin(')e^θ + cos(')e^' • Non-zero derivatives of basis vectors @θe^r = e^θ ;@θe^θ = −e^r ;@'e^r = sin(θ)e^' ;@'e^θ = cos(θ)e^' and @'e^' = −e^ρ • r-operator e^ e^ r = e^ @ + θ @ + ' @ r r r θ r sin(θ) ' • Divergence @ (r2F ) @ (sin(θ)F ) @ F r · F = r r + θ θ + ' ' r2 r sin(θ) r sin(θ) • Rotation e^ e^ @ F e^ r × F = r [@ (F sin(θ)) − @ F ] + θ −@ (rF ) + ' r + ' [@ (rF ) − @ F ] r sin(θ) θ ' ' θ r r ' sin(θ) r r θ θ r • Laplace 1 2 1 1 2 ∆ = @r(r @r) + @θ(sin(θ)@θ) + @ r2 r2 sin(θ) r2 sin2(θ) '.
Recommended publications
  • A.7 Orthogonal Curvilinear Coordinates
    ~NSORS Orthogonal Curvilinear Coordinates 569 )osition ated by converting its components (but not the unit dyads) to spherical coordinates, and 1 r+r', integrating each over the two spherical angles (see Section A.7). The off-diagonal terms in Eq. (A.6-13) vanish, again due to the symmetry. (A.6-6) A.7 ORTHOGONAL CURVILINEAR COORDINATES (A.6-7) Enormous simplificatons are achieved in solving a partial differential equation if all boundaries in the problem correspond to coordinate surfaces, which are surfaces gener­ (A.6-8) ated by holding one coordinate constant and varying the other two. Accordingly, many special coordinate systems have been devised to solve problems in particular geometries. to func­ The most useful of these systems are orthogonal; that is, at any point in space the he usual vectors aligned with the three coordinate directions are mutually perpendicular. In gen­ st calcu­ eral, the variation of a single coordinate will generate a curve in space, rather than a :lrSe and straight line; hence the term curvilinear. In this section a general discussion of orthogo­ nal curvilinear systems is given first, and then the relationships for cylindrical and spher­ ical coordinates are derived as special cases. The presentation here closely follows that in Hildebrand (1976). ial prop­ Base Vectors ve obtain Let (Ul, U2' U3) represent the three coordinates in a general, curvilinear system, and let (A.6-9) ei be the unit vector that points in the direction of increasing ui• A curve produced by varying U;, with uj (j =1= i) held constant, will be referred to as a "u; curve." Although the base vectors are each of constant (unit) magnitude, the fact that a U; curve is not gener­ (A.6-1O) ally a straight line means that their direction is variable.
    [Show full text]
  • Area and Volume in Curvilinear Coordinates 1 Overview
    Area and Volume in curvilinear coordinates 1 Overview There are several reasons why we need to have a way of measuring ares and volumes in general relativity. the gravitational field of a point mass m will be singular at the location of the mass, so we typically compute the field of a mass density; i.e., the mass contained in a unit volume. If we are dealing with the Gauss law of electrodynamics (and there will be a similar law for gravity) then we need to compute the flux through an area. How should areas and volumes be defined in general metrics? Let us start in flat 3-dimensional space with Cartesian coordinates. Consider two vectors V;~ W~ . These vectors describe the sides of a parallelogram, whose area can be written as A~ = V~ × W~ (1) The cross-product has magnitude jA~j = jV~ jjW~ j sin θ (2) and is seen to equal the area of the parallelogram. The direction of the cross product also serves a useful purpose: it gives the normal direction to the area spanned by the vectors. Note that there are two normal directions to the area { pointing above the plane and below the plane { and the `right hand rule' for cross products is a convention that picks out one of these over the other. Thus there is a choice of convention in the choice of normal; we will see this fact again later. The cross product can be written in terms of a determinant 0 x^ y^ z^ 1 V~ × W~ = @ V x V y V z A (3) W x W y W z A determinant is computed by computing a product of numbers, taking one number from each row, while making sure that we also take only one number from each column.
    [Show full text]
  • Curvilinear Coordinate Systems
    Appendix A Curvilinear coordinate systems Results in the main text are given in one of the three most frequently used coordinate systems: Cartesian, cylindrical or spherical. Here, we provide the material necessary for formulation of the elasticity problems in an arbitrary curvilinear orthogonal coor- dinate system. We then specify the elasticity equations for four coordinate systems: elliptic cylindrical, bipolar cylindrical, toroidal and ellipsoidal. These coordinate systems (and a number of other, more complex systems) are discussed in detail in the book of Blokh (1964). The elasticity equations in general curvilinear coordinates are given in terms of either tensorial or physical components. Therefore, explanation of the relationship between these components and the necessary background are given in the text to follow. Curvilinear coordinates ~i (i = 1,2,3) can be specified by expressing them ~i = ~i(Xl,X2,X3) in terms of Cartesian coordinates Xi (i = 1,2,3). It is as- sumed that these functions have continuous first partial derivatives and Jacobian J = la~i/axjl =I- 0 everywhere. Surfaces ~i = const (i = 1,2,3) are called coordinate surfaces and pairs of these surfaces intersects along coordinate curves. If the coordinate surfaces intersect at an angle :rr /2, the curvilinear coordinate system is called orthogonal. The usual summation convention (summation over repeated indices from 1 to 3) is assumed, unless noted otherwise. Metric tensor In Cartesian coordinates Xi the square of the distance ds between two neighboring points in space is determined by ds 2 = dxi dXi Since we have ds2 = gkm d~k d~m where functions aXi aXi gkm = a~k a~m constitute components of the metric tensor.
    [Show full text]
  • Appendix a Orthogonal Curvilinear Coordinates
    Appendix A Orthogonal Curvilinear Coordinates Given a symmetry for a system under study, the calculations can be simplified by choosing, instead of a Cartesian coordinate system, another set of coordinates which takes advantage of that symmetry. For example calculations in spherical coordinates result easier for systems with spherical symmetry. In this chapter we will write the general form of the differential operators used in electrodynamics and then give their expressions in spherical and cylindrical coordi- nates.1 A.1 Orthogonal curvilinear coordinates A system of coordinates u1, u2, u3, can be defined so that the Cartesian coordinates x, y and z are known functions of the new coordinates: x = x(u1, u2, u3) y = y(u1, u2, u3) z = z(u1, u2, u3) (A.1) Systems of orthogonal curvilinear coordinates are defined as systems for which locally, nearby each point P(u1, u2, u3), the surfaces u1 = const, u2 = const, u3 = const are mutually orthogonal. An elementary cube, bounded by the surfaces u1 = const, u2 = const, u3 = const, as shown in Fig. A.1, will have its edges with lengths h1du1, h2du2, h3du3 where h1, h2, h3 are in general functions of u1, u2, u3. The length ds of the line- element OG, one diagonal of the cube, in Cartesian coordinates is: ds = (dx)2 + (dy)2 + (dz)2 1The orthogonal coordinates are presented with more details in J.A. Stratton, Electromagnetic The- ory, McGraw-Hill, 1941, where many other useful coordinate systems (elliptic, parabolic, bipolar, spheroidal, paraboloidal, ellipsoidal) are given. © Springer International Publishing Switzerland 2016 189 F. Lacava, Classical Electrodynamics, Undergraduate Lecture Notes in Physics, DOI 10.1007/978-3-319-39474-9 190 Appendix A: Orthogonal Curvilinear Coordinates Fig.
    [Show full text]
  • Curvilinear Coordinates
    UNM SUPPLEMENTAL BOOK DRAFT June 2004 Curvilinear Analysis in a Euclidean Space Presented in a framework and notation customized for students and professionals who are already familiar with Cartesian analysis in ordinary 3D physical engineering space. Rebecca M. Brannon Written by Rebecca Moss Brannon of Albuquerque NM, USA, in connection with adjunct teaching at the University of New Mexico. This document is the intellectual property of Rebecca Brannon. Copyright is reserved. June 2004 Table of contents PREFACE ................................................................................................................. iv Introduction ............................................................................................................. 1 Vector and Tensor Notation ........................................................................................................... 5 Homogeneous coordinates ............................................................................................................. 10 Curvilinear coordinates .................................................................................................................. 10 Difference between Affine (non-metric) and Metric spaces .......................................................... 11 Dual bases for irregular bases .............................................................................. 11 Modified summation convention ................................................................................................... 15 Important notation
    [Show full text]
  • Appendix a Relations Between Covariant and Contravariant Bases
    Appendix A Relations Between Covariant and Contravariant Bases The contravariant basis vector gk of the curvilinear coordinate of uk at the point P is perpendicular to the covariant bases gi and gj, as shown in Fig. A.1.This contravariant basis gk can be defined as or or a gk g  g ¼  ðA:1Þ i j oui ou j where a is the scalar factor; gk is the contravariant basis of the curvilinear coordinate of uk. Multiplying Eq. (A.1) by the covariant basis gk, the scalar factor a results in k k ðgi  gjÞ: gk ¼ aðg : gkÞ¼ad ¼ a ÂÃk ðA:2Þ ) a ¼ðgi  gjÞ : gk gi; gj; gk The scalar triple product of the covariant bases can be written as pffiffiffi a ¼ ½¼ðg1; g2; g3 g1  g2Þ : g3 ¼ g ¼ J ðA:3Þ where Jacobian J is the determinant of the covariant basis tensor G. The direction of the cross product vector in Eq. (A.1) is opposite if the dummy indices are interchanged with each other in Einstein summation convention. Therefore, the Levi-Civita permutation symbols (pseudo-tensor components) can be used in expression of the contravariant basis. ffiffiffi p k k g g ¼ J g ¼ðgi  gjÞ¼Àðgj  giÞ eijkðgi  gjÞ eijkðgi  gjÞ ðA:4Þ ) gk ¼ pffiffiffi ¼ g J where the Levi-Civita permutation symbols are defined by 8 <> þ1ifði; j; kÞ is an even permutation; eijk ¼ > À1ifði; j; kÞ is an odd permutation; : A:5 0ifi ¼ j; or i ¼ k; or j ¼ k ð Þ 1 , e ¼ ði À jÞÁðj À kÞÁðk À iÞ for i; j; k ¼ 1; 2; 3 ijk 2 H.
    [Show full text]
  • Tensor Analysis and Curvilinear Coordinates.Pdf
    Tensor Analysis and Curvilinear Coordinates Phil Lucht Rimrock Digital Technology, Salt Lake City, Utah 84103 last update: May 19, 2016 Maple code is available upon request. Comments and errata are welcome. The material in this document is copyrighted by the author. The graphics look ratty in Windows Adobe PDF viewers when not scaled up, but look just fine in this excellent freeware viewer: http://www.tracker-software.com/pdf-xchange-products-comparison-chart . The table of contents has live links. Most PDF viewers provide these links as bookmarks on the left. Overview and Summary.........................................................................................................................7 1. The Transformation F: invertibility, coordinate lines, and level surfaces..................................12 Example 1: Polar coordinates (N=2)..................................................................................................13 Example 2: Spherical coordinates (N=3)...........................................................................................14 Cartesian Space and Quasi-Cartesian Space.......................................................................................15 Pictures A,B,C and D..........................................................................................................................16 Coordinate Lines.................................................................................................................................16 Example 1: Polar coordinates, coordinate lines
    [Show full text]
  • A Tensor Operations in Orthogonal Curvilinear Coordinate Systems
    A Tensor Operations in Orthogonal Curvilinear Coordinate Systems A.1 Change of Coordinate System The vector and tensor operations we have discussed in the foregoing chapters were performed solely in rectangular coordinate system. It should be pointed out that we were dealing with quantities such as velocity, acceleration, and pressure gradient that are independent of any coordinate system within a certain frame of reference. In this connection it is necessary to distinguish between a coordinate system and a frame of reference. The following example should clarify this distinction. In an absolute frame of reference, the flow velocity vector may be described by the rectangular Cartesian coordinate xi: (A.1) It may also be described by a cylindrical coordinate system, which is a non-Cartesian coordinate system: (A.2) or generally by any other non-Cartesian or curvilinear coordinate ȟi that describes the flow channel geometry: (A.3) By changing the coordinate system, the flow velocity vector will not change. It remains invariant under any transformation of coordinates. This is true for any other quantities such as acceleration, force, pressure or temperature gradient. The concept of invariance, however, is generally no longer valid if we change the frame of reference. For example, if the flow particles leave the absolute frame of reference and enter the relative frame of reference, for example a moving or rotating frame, its velocity will experience a change. In this Chapter, we will pursue the concept of quantity invariance and discuss the fundamentals that are needed for coordinate transformation. A.2 Co- and Contravariant Base Vectors, Metric Coefficients As we saw in the previous chapter, a vector quantity is described in Cartesian coordinate system xi by its components: M.T.
    [Show full text]
  • Orthogonal Curvilinear Coordinates 28.3   
    ® Orthogonal Curvilinear Coordinates 28.3 Introduction The derivatives div, grad and curl from Section 28.2 can be carried out using coordinate systems other than the rectangular Cartesian coordinates. This Section shows how to calculate these derivatives in other coordinate systems. Two coordinate systems - cylindrical polar coordinates and spherical polar coordinates - will be illustrated. • be able to find the gradient, divergence and Prerequisites curl of a field in Cartesian coordinates Before starting this Section you should ... • be familiar with polar coordinates • find the divergence, gradient or curl of a Learning Outcomes vector or scalar field expressed in terms of On completion you should be able to ... orthogonal curvilinear coordinates HELM (2008): 37 Section 28.3: Orthogonal Curvilinear Coordinates 1. Orthogonal curvilinear coordinates The results shown in Section 28.2 have been given in terms of the familiar Cartesian (x, y, z) co- ordinate system. However, other coordinate systems can be used to better describe some physical situations. A set of coordinates u = u(x, y, z), v = v(x, y, z) and w = w(x, y, z) where the direc- tions at any point indicated by u, v and w are orthogonal (perpendicular) to each other is referred to as a set of orthogonal curvilinear coordinates. With each coordinate is associated a scale factor q ∂x 2 ∂y 2 ∂z 2 hu, hv or hw respectively where hu = ∂u + ∂u + ∂u (with similar expressions for hv and hw). The scale factor gives a measure of how a change in the coordinate changes the position of a point. Two commonly-used sets of orthogonal curvilinear coordinates are cylindrical polar coordinates and spherical polar coordinates.
    [Show full text]
  • Appendix A: Relations Between Covariant and Contravariant Bases
    Appendix A: Relations Between Covariant and Contravariant Bases The contravariant basis vector gk of the curvilinear coordinate of uk at the point P is perpendicular to the covariant bases gi and gj, as shown in Fig. A.1. This contravariant basis gk can be defined as ∂r ∂r α gk g  g ¼  ðA:1Þ i j ∂ui ∂uj where α is the scalar factor; gk is the contravariant basis of the curvilinear coordi- nate of uk. g3 u3 g3 x3 g1 P g e 2 3 g2 u2 u1 g 0 e 1 e1 2 x2 x1 Fig. A.1 Covariant and contravariant bases of curvilinear coordinates © Springer-Verlag Berlin Heidelberg 2017 313 H. Nguyen-Scha¨fer, J.-P. Schmidt, Tensor Analysis and Elementary Differential Geometry for Physicists and Engineers, Mathematical Engineering, DOI 10.1007/978-3-662-48497-5 314 Appendix A: Relations Between Covariant and Contravariant Bases Multiplying Eq. (A.1) by the covariant basis gk, the scalar factor α results in ÀÁ  Á ¼ α k Á ¼ αδk ¼ α gi gj gk g gk k hi ðA:2Þ ) α ¼  Á ; ; gi gj gk gi gj gk The scalar triple product of the covariant bases can be written as pffiffiffi α ¼ ½¼; ; ðÞÁ ¼ ¼ ð : Þ g1 g2 g3 g1 g2 g3 g J A 3 where Jacobian J is the determinant of the covariant basis tensor G. The direction of the cross-product vector in Eq. (A.1) is opposite if the dummy indices are interchanged with each other in Einstein summation convention. There- fore, the Levi-Civita permutation symbols (pseudo-tensor components) can be used in expression of the contravariant basis.
    [Show full text]
  • Class: M. Sc. Subject: Mathematical Methods for Physics
    Class: M. Sc. Subject: Mathematical Methods for Physics I have taken all course materials for Mathematical Methods for Physics from G. B. Arfken and H.J. Weber, Mathematical Methods for Physicists. Students can download this book form given web address; Web address: https://b-ok.cc/book/446170/6e805d All topics have been taken from Chapter-2 & Chapter-4 from above said book ( https://b- ok.cc/book/446170/6e805d ). I am sending pdf file of Chapter 2 and Chapter 4 . Name of Topics Vector Analysis in Curved Coordinates and Tensors, Orthogonal Coordinates in R3, Differential Vector Operators, Special Coordinate Systems, Circular Cylinder Coordinates, Spherical and Polar Coordinates. Tensor Analysis, Contraction, Direct Product, Quotient Rule, Pseudotensors, Dual Tensors, General Tensors, Tensor Derivative Operators. Group Theory, Introduction, Generators of Continuous Group, Orbital Angular Momentum, Angular Momentum Coupling, Homogeneous Lorentz Group, Discrete Groups CHAPTER 2 VECTOR ANALYSIS IN CURVED COORDINATES AND TENSORS In Chapter 1 we restricted ourselves almost completely to rectangular or Cartesian coordi- nate systems. A Cartesian coordinate system offers the unique advantage that all three unit vectors, x, y, and z, are constant in direction as well as in magnitude. We did introduce the radial distanceˆ ˆ r, butˆ even this was treated as a function of x, y, and z. Unfortunately, not all physical problems are well adapted to a solution in Cartesian coordinates. For instance, if we have a central force problem, F rF(r), such as gravitational or electrostatic force, Cartesian coordinates may be unusually= inappropriate.ˆ Such a problem demands the use of a coordinate system in which the radial distance is taken to be one of the coordinates, that is, spherical polar coordinates.
    [Show full text]
  • Gradient, Divergence and Curl in Curvilinear Coordinates
    Gradient, Divergence and Curl in Curvilinear Coordinates Although cartesian orthogonal coordinates are very intuitive and easy to use, it is often found more convenient to work with other coordinate systems. Being able to change all variables and expression involved in a given problem, when a different coordinate system is chosen, is one of those skills a physicist, and even more a theoretical physicist, needs to possess. In this lecture a general method to express any variable and expression in an arbitrary curvilinear coordinate system will be introduced and explained. We will be mainly interested to find out gen- eral expressions for the gradient, the divergence and the curl of scalar and vector fields. Specific applications to the widely used cylindrical and spherical systems will conclude this lecture. 1 The concept of orthogonal curvilinear coordinates The cartesian orthogonal coordinate system is very intuitive and easy to handle. Once an origin has been fixed in space and three orthogonal scaled axis are anchored to this origin, any point in space is uniquely determined by three real numbers, its cartesian coordinates. A curvilinear coordinate system can be defined starting from the orthogonal cartesian one. If x, y, z are the cartesian coordinates, the curvilinear ones, u, v, w, can be expressed as smooth functions of x, y, z, according to: u = u(x; y; z) v = v(x; y; z) (1) w = w(x; y; z) These functions can be inverted to give x, y, z-dependency on u, v, w: x = x(u; v; w) y = y(u; v; w) (2) z = z(u; v; w) There are infinitely many curvilinear systems that can be defined using equations (1) and (2).
    [Show full text]