User:Tohline/Appendix/Mathematics/ScaleFactors
Scale Factors for Orthogonal Curvilinear Coordinate Systems
Here we lean heavily on the class notes and associated references that have been provided by P. A. Kelly in a collection titled, Mechanics Lecture Notes: An Introduction to Solid Mechanics, as they appeared online in early 2021. See especially the subsection of Part III in which the properties of Vectors and Tensors are discussed.
| Tiled Menu | Tables of Content | Banner Video | Tohline Home Page | |
Getting Started
Following Kelly, we will use <math>~\hat{e}_i</math> and <math>~x_i</math> when referencing, respectively, the three (i = 1,3) basis vectors and coordinate "curves" of the Cartesian coordinate system; and we will use <math>~\hat{g}_i</math> and <math>~\Theta_i</math> when referencing, respectively, the three (i = 1,3) basis vectors and coordinate curves of some other, curvilinear coordinate system.
2D Oblique Coordinate System Example
Consider a vector, <math>~\vec{v}</math>, which in Cartesian coordinates is described by the expression,
<math>~\vec{v}</math> |
<math>~=</math> |
<math>~ \hat{e}_1 v_x + \hat{e}_2 v_y \, . </math> |
Referencing Figure 1.16.4 of Kelly's Part III, we appreciate that in a two-dimensional (2D) oblique coordinate system where <math>~\alpha</math> is the (less than 90°) angle between the two basis vectors, the same vector will be represented by the expression,
<math>~\vec{v}</math> |
<math>~=</math> |
<math>~ \hat{g}_1 v^1 + \hat{g}_2 v^2 \, . </math> |
The angle between <math>~\hat{g}_2</math> and <math>~\hat{e}_2</math> is, (π/2 - α), so we appreciate that,
<math>~v_y</math> |
<math>~=</math> |
<math>~v^2\cos\biggl(\frac{\pi}{2} - \alpha \biggr) = v^2 \sin\alpha</math> |
<math>~\Rightarrow~~~v^2</math> |
<math>~=</math> |
<math>~\frac{v_y}{\sin\alpha} \, .</math> |
Next, from a visual inspection of the figure, we appreciate that <math>~v_x</math> is longer than <math>~v^1</math> by the amount, <math>~v^2\cos\alpha</math>; that is,
<math>~v_x</math> |
<math>~=</math> |
<math>~v^1 + v^2\cos\alpha = v_1 + \frac{v_y}{\tan\alpha}</math> |
<math>~\Rightarrow ~~~ v^1</math> |
<math>~=</math> |
<math>~v_x - \frac{v_y}{\tan\alpha} \, .</math> |
(These are the same pair of transformation relations that appear as Eq. (1.16.3) of Kelly's Part III.)
Covarient: The set of basis vectors, <math>~\hat{g}_1</math> and <math>~\hat{g}_2</math> (note the subscript indices), that are aligned with the coordinate directions, <math>~\Theta_1</math> and <math>~\Theta_2</math>, are generically referred to as covariant base vectors. Contravarient: A second set of vectors, which will be termed contravariant base vectors, <math>~\hat{g}^1</math> and <math>~\hat{g}^2</math> (denoted by superscript indices), will be aligned with a new set of coordinate directions, <math>~\Theta^1</math> and <math>~\Theta^2</math>. This new set of base vectors is defined as follows (see Fig. 1.15.5 of Kelly's Part III): the base vector <math>~\hat{g}^1</math> is perpendicular to <math>~\hat{g}_1</math> — that is, <math>~\hat{g}^1 \cdot \hat{g}_2 = 0</math> — and the base vector <math>~\hat{g}^2</math> is perpendicular to <math>~\hat{g}_2</math> — that is, <math>~\hat{g}_1 \cdot \hat{g}^2 = 0</math>. Further, we ensure that,
(Verbatim from p. 138 of Kelly's Part III) A good trick for remembering which are the covariant and which are the contravariant is that the third letter of the word tells us whether the word is associated with subscripts or with superscripts. In "covariant," the "v" is pointing down, so we use subscripts; for "contravariant," the "n" is (with a bit of imagination) pointing up, so we use superscripts. |
Continuing with our 2D oblique coordinate system example and appreciating that Kelly has chosen to align the <math>~\hat{g}_1</math> basis vector with the <math>~\hat{e}_1</math> (Cartesian) basis vector, we see that the transformation between the two sets of covariant basis vectors is given by the relations,
<math>~\hat{g}_1 = \hat{e}_1 \, ,</math> |
and, |
<math>~\hat{g}_2 = \hat{e}_1\cos\alpha + \hat{e}_2\sin\alpha \, .</math> |
These conditions lead to the following complementary set of contravariant basis vectors:
<math>~\hat{g}^1 = \hat{e}_1 - \hat{e}_2 \biggl( \frac{1}{\tan\alpha}\biggr) \, ,</math> |
and, |
<math>~\hat{g}^2 = \hat{e}_2 \biggl( \frac{1}{\sin\alpha} \biggr) \, .</math> |
Note that, as defined herein, the magnitude (i.e., scalar lengths) of these contravariant basis vectors is not unity; they are, instead,
<math>~|\hat{g}^1| \equiv \biggl[ \hat{g}^1 \cdot \hat{g}^1 \biggr]^{1 / 2}</math> |
<math>~=</math> |
<math>~ \biggl\{ \biggl[\hat{e}_1 - \hat{e}_2 \biggl( \frac{1}{\tan\alpha}\biggr)\biggr] \cdot \biggl[\hat{e}_1 - \hat{e}_2 \biggl( \frac{1}{\tan\alpha}\biggr)\biggr] \biggr\}^{1 / 2} </math> |
|
<math>~=</math> |
<math>~ \biggl\{ 1 + \frac{1}{\tan^2\alpha} \biggr\}^{1 / 2} = \frac{1}{\sin\alpha} \, ; </math> |
<math>~|\hat{g}^2| \equiv \biggl[ \hat{g}^2 \cdot \hat{g}^2 \biggr]^{1 / 2}</math> |
<math>~=</math> |
<math>~ \frac{1}{\sin\alpha} \, . </math> |
Some Useful Relations
Examples
Once expressions for the nine separate direction cosines are known for a system of orthogonal coordinates, then the following hold:
<math>~\mathbf{\hat{g}}_n</math> |
<math>~=</math> |
<math>~\hat\imath \gamma_{n1} + \hat\jmath \gamma_{n2} + \hat{k} \gamma_{n3} \, ;</math> |
and,
<math>~\hat\imath</math> |
<math>~=</math> |
<math>~ \mathbf{\hat{g}}_1 \gamma_{11} + \mathbf{\hat{g}}_2 \gamma_{21} + \mathbf{\hat{g}}_3 \gamma_{31} \, , </math> |
<math>~\hat\jmath</math> |
<math>~=</math> |
<math>~ \mathbf{\hat{g}}_1 \gamma_{12} + \mathbf{\hat{g}}_2 \gamma_{22} + \mathbf{\hat{g}}_3 \gamma_{32} \, , </math> |
<math>~\hat{k}</math> |
<math>~=</math> |
<math>~ \mathbf{\hat{g}}_1 \gamma_{13} + \mathbf{\hat{g}}_2 \gamma_{23} + \mathbf{\hat{g}}_3 \gamma_{33} \, . </math> |
Hence, the position vector is given by the expression,
<math>~\mathbf{\vec{x}} = \hat\imath x + \hat\jmath y + \hat{k}z</math> |
<math>~=</math> |
<math>~ \mathbf{\hat{g}}_1 (\gamma_{11} x + \gamma_{12} y + \gamma_{13} z) + \mathbf{\hat{g}}_2 (\gamma_{21} x + \gamma_{22} y + \gamma_{23} z) + \mathbf{\hat{g}}_3 (\gamma_{31} x + \gamma_{32} y + \gamma_{33} z) </math> |
Cylindrical Coordinates
This is drawn principally from Example #1 (starting on p. 148) of Kelly.
<math>~\lambda_1</math> |
<math>~=</math> |
<math>~ \varpi \equiv \biggl[(x^1)^2 + (x^2)^2 \bigg]^{1 / 2} \, , </math> |
<math>~\lambda_2</math> |
<math>~\equiv</math> |
<math>~ \tan^{-1}\biggl[\frac{x^2}{x^1} \biggr] \, , </math> |
<math>~\lambda_3</math> |
<math>~\equiv</math> |
<math>~ x^3 \, . </math> |
Direction Cosine Components for Cylindrical Coordinates | ||||||||
<math>~n</math> | <math>~\lambda_n</math> | <math>~h_n</math> | <math>~\frac{\partial \lambda_n}{\partial x}</math> | <math>~\frac{\partial \lambda_n}{\partial y}</math> | <math>~\frac{\partial \lambda_n}{\partial z}</math> | <math>~\gamma_{n1}</math> | <math>~\gamma_{n2}</math> | <math>~\gamma_{n3}</math> |
<math>~1</math> | <math>~\varpi \equiv (x^2 + y^2 )^{1 / 2} </math> | <math>~1</math> | <math>~\frac{x}{\lambda_1}</math> | <math>~\frac{y}{\lambda_1}</math> | <math>~0</math> | <math>~\frac{x}{\lambda_1}</math> | <math>~\frac{y}{\lambda_1}</math> | <math>~0</math> |
<math>~2</math> | <math>~\varphi \equiv \tan^{-1}\biggl[\frac{y}{x}\biggr]</math> | <math>~\lambda_1</math> | <math>~- \frac{y}{\varpi^2}</math> | <math>~\frac{x}{\varpi^2}</math> | <math>~0</math> | <math>~- \frac{y}{\varpi}</math> | <math>~\frac{x}{\varpi}</math> | <math>~0</math> |
<math>~3</math> | <math>~z</math> | <math>~1</math> | <math>~0</math> | <math>~0</math> | <math>~1</math> | <math>~0</math> | <math>~0</math> | <math>~1</math> |
<math>~\mathbf{ \hat{g}}_1</math> |
<math>~=</math> |
<math>~ \hat\imath \gamma_{11} + \hat\jmath \gamma_{12} + \hat{k} \gamma_{13} = \hat\imath \biggl( \frac{x}{\varpi} \biggr) + \hat\jmath \biggl( \frac{y}{\varpi} \biggr) = \hat\imath \cos\varphi + \hat\jmath \sin\varphi \, , </math> |
<math>~\mathbf{ \hat{g}}_2</math> |
<math>~=</math> |
<math>~ \hat\imath \gamma_{21} + \hat\jmath \gamma_{22} + \hat{k} \gamma_{23} = - \hat\imath \biggl( \frac{y}{\varpi} \biggr) + \hat\jmath \biggl( \frac{x}{\varpi} \biggr) = \hat\imath \sin\varphi + \hat\jmath \cos\varphi \, , </math> |
<math>~\mathbf{ \hat{g}}_3</math> |
<math>~=</math> |
<math>~ \hat\imath \gamma_{31} + \hat\jmath \gamma_{32} + \hat{k} \gamma_{33} = \hat{k} \, . </math> |
And the position vector is,
<math>~\mathbf{\vec{x}} </math> |
<math>~=</math> |
<math>~ \mathbf{\hat{g}}_1 (\gamma_{11} x + \gamma_{12} y + \gamma_{13} z) + \mathbf{\hat{g}}_2 (\gamma_{21} x + \gamma_{22} y + \gamma_{23} z) + \mathbf{\hat{g}}_3 (\gamma_{31} x + \gamma_{32} y + \gamma_{33} z) </math> |
|
<math>~=</math> |
<math>~ \mathbf{\hat{g}}_1 \biggl[\frac{x^2}{\varpi} + \frac{y^2}{\varpi} \biggr] + \mathbf{\hat{g}}_2 \biggl[-\frac{xy}{\varpi} + \frac{xy}{\varpi} \biggr] + \mathbf{\hat{g}}_3 z </math> |
|
<math>~=</math> |
<math>~ \mathbf{\hat{g}}_1 \varpi + \mathbf{\hat{g}}_3 z \, . </math> |
The line element is,
<math>~ds^2</math> |
<math>~=</math> |
<math>~ \sum_{i=1}^3 h_i^2 d\lambda_i^2 = d\lambda_1^2 + \varpi^2 d\lambda_2^2 + d\lambda_3^2 \, . </math> |
In terms of Cartesian basis vector, this is,
<math>~ds^2</math> |
<math>~=</math> |
<math>~ h_1^2 \biggl[ \biggl( \frac{\partial \lambda_1}{\partial x}\biggr)dx + \biggl( \frac{\partial \lambda_1}{\partial y}\biggr)dy + \biggl( \frac{\partial \lambda_1}{\partial z}\biggr)dz \biggl]^2 + h_2^2 \biggl[ \biggl( \frac{\partial \lambda_2}{\partial x}\biggr)dx + \biggl( \frac{\partial \lambda_2}{\partial y}\biggr)dy + \biggl( \frac{\partial \lambda_2}{\partial z}\biggr)dz \biggl]^2 + h_3^2 \biggl[ \biggl( \frac{\partial \lambda_3}{\partial x}\biggr)dx + \biggl( \frac{\partial \lambda_3}{\partial y}\biggr)dy + \biggl( \frac{\partial \lambda_3}{\partial z}\biggr)dz \biggl]^2 </math> |
|
<math>~=</math> |
<math>~ \biggl[ \gamma_{11} dx + \gamma_{12} dy + \gamma_{13} dz \biggl]^2 + \biggl[ \gamma_{21} dx + \gamma_{22} dy + \gamma_{23} dz \biggl]^2 + \biggl[ \gamma_{31} dx + \gamma_{32} dy + \gamma_{33} dz \biggl]^2 </math> |
|
<math>~=</math> |
<math>~ \biggl[ \biggl(\frac{x}{\varpi}\biggr) dx + \biggl(\frac{y}{\varpi}\biggr) dy \biggl]^2 + \biggl[- \biggl(\frac{y}{\varpi}\biggr) dx + \biggl(\frac{x}{\varpi}\biggr) dy \biggl]^2 + dz^2 </math> |
|
<math>~=</math> |
<math>~ \biggl(\frac{x}{\varpi}\biggr)^2 dx^2 + \biggl(\frac{y}{\varpi}\biggr)^2 dy^2 + 2 \biggl(\frac{x}{\varpi}\biggr)\biggl(\frac{y}{\varpi}\biggr)dxdy ~-~ 2\biggl(\frac{y}{\varpi}\biggr)\biggl(\frac{x}{\varpi}\biggr) dx dy + \biggl(\frac{y}{\varpi}\biggr)^2 dx^2 + \biggl(\frac{x}{\varpi}\biggr)^2 dy^2 + dz^2 </math> |
|
<math>~=</math> |
<math>~ \biggl[ \biggl(\frac{x}{\varpi}\biggr)^2+ \biggl(\frac{y}{\varpi}\biggr)^2 \biggr] dx^2 + \biggl[ \biggl(\frac{y}{\varpi}\biggr)^2 + \biggl(\frac{x}{\varpi}\biggr)^2 \biggr] dy^2 + dz^2 </math> |
|
<math>~=</math> |
<math>~ dx^2 + dy^2 + dz^2 \, . </math> Yes! |
T10 Coordinates
Position Vector
Pulling from our accompanying Table of Direction Cosine Components for T10 Coordinates, the position vector is given by the expression,
<math>~\mathbf{\vec{x}} </math> |
<math>~=</math> |
<math>~ \mathbf{\hat{g}}_1 (\gamma_{11} x + \gamma_{12} y + \gamma_{13} z) + \mathbf{\hat{g}}_2 (\gamma_{41} x + \gamma_{42} y + \gamma_{43} z) + \mathbf{\hat{g}}_3 (\gamma_{51} x + \gamma_{52} y + \gamma_{53} z) </math> |
|
<math>~=</math> |
<math>~ \mathbf{\hat{g}}_1 \biggl[ \lambda_1^2 \ell_{3D} \biggr] + \mathbf{\hat{g}}_2 \biggl[ 1 + \frac{1}{q^2} - \frac{2}{p^2} \biggr]\frac{xq^2y p^2z}{\mathcal{D}} + \mathbf{\hat{g}}_3 \biggl[ -x^2(2q^4y^2 + p^4z^2) + q^2y^2(p^4z^2 + 2x^2) + p^2z^2(x^2-q^4y^2) \biggr] \frac{\ell_{3D}}{\mathcal{D}} </math> |
Line Element
In terms of Cartesian basis vector, the line element is,
<math>~ds^2</math> |
<math>~=</math> |
<math>~ h_1^2 \biggl[ \biggl( \frac{\partial \lambda_1}{\partial x}\biggr)dx + \biggl( \frac{\partial \lambda_1}{\partial y}\biggr)dy + \biggl( \frac{\partial \lambda_1}{\partial z}\biggr)dz \biggl]^2 + h_2^2 \biggl[ \biggl( \frac{\partial \lambda_4}{\partial x}\biggr)dx + \biggl( \frac{\partial \lambda_4}{\partial y}\biggr)dy + \biggl( \frac{\partial \lambda_4}{\partial z}\biggr)dz \biggl]^2 + h_3^2 \biggl[ \biggl( \frac{\partial \lambda_5}{\partial x}\biggr)dx + \biggl( \frac{\partial \lambda_5}{\partial y}\biggr)dy + \biggl( \frac{\partial \lambda_5}{\partial z}\biggr)dz \biggl]^2 </math> |
|
<math>~=</math> |
<math>~ \biggl[ \gamma_{11} dx + \gamma_{12} dy + \gamma_{13} dz \biggl]^2 + \biggl[ \gamma_{41} dx + \gamma_{42} dy + \gamma_{43} dz \biggl]^2 + \biggl[ \gamma_{51} dx + \gamma_{52} dy + \gamma_{53} dz \biggl]^2 </math> |
|
<math>~=</math> |
<math>~ \biggl[ (x) dx + (q^2 y) dy + (p^2 z) dz \biggl]^2 \ell_{3D}^2 + \biggl[ (q^2y p^2z) dx +( x p^2z) dy - ( 2xq^2y) dz \biggl]^2 \frac{1}{\mathcal{D}^2} + \biggl[ -x(2q^4y^2 + p^4z^2) dx + q^2y(p^4z^2 + 2x^2) dy + p^2z(x^2 - q^4y^2) dz \biggl]^2 \frac{\ell_{3D}^2}{\mathcal{D}^2} </math> |
<math>~\Rightarrow ~~~ \biggl( \frac{\mathcal{D}}{\ell_{3D}}\biggr)^2 ds^2</math> |
<math>~=</math> |
<math>~ \biggl[ (x) dx + (q^2 y) dy + (p^2 z) dz \biggl]^2 \mathcal{D}^2 + \biggl[ (q^2y p^2z) dx +( x p^2z) dy - ( 2xq^2y) dz \biggl]^2 (\ell_{3D})^{-2} + \biggl[ -x(2q^4y^2 + p^4z^2) dx + q^2y(p^4z^2 + 2x^2) dy + p^2z(x^2 - q^4y^2) dz \biggl]^2 \, . </math> |
Notice that the coefficient that corresponds to each term is given by the following expressions:
<math>~[dx^2]:</math> |
|
<math>~ x^2 \mathcal{D}^2 + q^4y^2 p^4z^2 \ell_{3D}^{-2} + x^2(2q^4y^2 + p^4z^2)^2 </math> |
|
<math>~=</math> |
<math>~ x^2 (q^4y^2p^4z^2 + x^2p^4z^2 + 4x^2q^4y^2) + q^4y^2 p^4z^2 (x^2 + q^4 y^2 + p^4z^2) + x^2(4q^8y^4 + 4q^4y^2p^4z^2 + p^8z^4) </math> |
|
<math>~=</math> |
<math>~ 6x^2 q^4y^2p^4z^2 + x^4(q^4y^2 + p^4z^2) +q^8z^4(4x^2 + p^4z^2) + p^8z^2(x^2 + q^4y^2) </math> |
|
<math>~=</math> |
<math>~ \biggl( \frac{\mathcal{D}}{\ell_{3D}}\biggr)^2 \, . </math> |
<math>~[dx\cdot dy]:</math> |
|
<math>~ 2xq^2y \mathcal{D}^2 +2q^2yp^2z \cdot xp^2z (\ell_{3D})^{-2} - 2x(2q^4y^2 + p^4z^2)q^2y(p^4z^2 + 2x^2) </math> |
|
<math>~=</math> |
<math>~ 2xq^2y[\mathcal{D}^2 + p^4z^2(\ell_{3D})^{-2} - (2q^4y^2 + p^4z^2)(p^4z^2 + 2x^2)] </math> |
|
<math>~=</math> |
<math>~ 2xq^2y[q^4y^2p^4z^2 + x^2p^4z^2 + 4x^2q^4y^2 + p^4z^2(x^2 + q^4 y^2 + p^4 z^2) - 2q^4y^2 (p^4z^2 + 2x^2) - p^4z^2(p^4z^2 + 2x^2)] </math> |
|
<math>~=</math> |
<math>~ 2xq^2y[ ~0~] = 0\, . </math> |
Volume Element
|
|
|
<math>~dV = h_1 h_4 h_5 ~d\lambda_1 d\lambda_4 d\lambda_5</math> |
<math>~=</math> |
<math>~ \biggl[ \gamma_{11} dx + \gamma_{14} dy + \gamma_{15} dz \biggl] \biggl[ \gamma_{41} dx + \gamma_{44} dy + \gamma_{45} dz \biggl] \biggl[ \gamma_{51} dx + \gamma_{54} dy + \gamma_{55} dz \biggl] </math> |
|
<math>~=</math> |
<math>~ \biggl[ \gamma_{11} dx + \gamma_{14} dy + \gamma_{15} dz \biggl] \biggl[ dx^2 \gamma_{41}\gamma_{51} + dy^2 \gamma_{44}\gamma_{54} + dz^2 \gamma_{45}\gamma_{55} + dx dy ( \gamma_{41}\gamma_{54} + \gamma_{44}\gamma_{51} ) + dx dz ( \gamma_{41}\gamma_{55} + \gamma_{45}\gamma_{51} ) + dy dz ( \gamma_{44}\gamma_{55} + \gamma_{45}\gamma_{54} ) \biggr] </math> |
|
<math>~=</math> |
<math>~ dx^3 \gamma_{11} \gamma_{41}\gamma_{51} + dx dy^2 \gamma_{11} \gamma_{44}\gamma_{54} + dx dz^2 \gamma_{11} \gamma_{45}\gamma_{55} + dx^2 dy \gamma_{11} ( \gamma_{41}\gamma_{54} + \gamma_{44}\gamma_{51} ) + dx^2 dz \gamma_{11} ( \gamma_{41}\gamma_{55} + \gamma_{45}\gamma_{51} ) + dx dy dz \gamma_{11} ( \gamma_{44}\gamma_{55} + \gamma_{45}\gamma_{54} ) \biggr] </math> |
|
|
<math>~ + dx^2 dy \gamma_{14} \gamma_{41}\gamma_{51} + dy^3 \gamma_{14} \gamma_{44}\gamma_{54} + dy dz^2 \gamma_{14} \gamma_{45}\gamma_{55} + dx dy^2 \gamma_{14} ( \gamma_{41}\gamma_{54} + \gamma_{44}\gamma_{51} ) + dx dy dz \gamma_{14} ( \gamma_{41}\gamma_{55} + \gamma_{45}\gamma_{51} ) + dy^2 dz \gamma_{14} ( \gamma_{44}\gamma_{55} + \gamma_{45}\gamma_{54} ) \biggr] </math> |
|
|
<math>~ + dx^2 dz \gamma_{15} \gamma_{41}\gamma_{51} + dy^2 dz \gamma_{15} \gamma_{44}\gamma_{54} + dz^3 \gamma_{15} \gamma_{45}\gamma_{55} + dx dy dz \gamma_{15} ( \gamma_{41}\gamma_{54} + \gamma_{44}\gamma_{51} ) + dx dz^2 \gamma_{15} ( \gamma_{41}\gamma_{55} + \gamma_{45}\gamma_{51} ) + dy dz^2 \gamma_{15} ( \gamma_{44}\gamma_{55} + \gamma_{45}\gamma_{54} ) \biggr] </math> |
See Also
© 2014 - 2021 by Joel E. Tohline |