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.

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, <math>~\hat{g}_1 \cdot \hat{g}^1 = 1 \, ,</math>       and,       <math>~\hat{g}_2 \cdot \hat{g}^2 = 1 \, .</math> (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>

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>~(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>~\tan^{-1}\biggl[\frac{y}{x}\biggr]</math>	<math>~\frac{1}{\lambda_2}\biggl[\frac{x q^2 y p^2 z}{ \mathcal{D}}\biggr] </math>	<math>~\frac{\lambda_2}{x}</math>	<math>~\frac{\lambda_2}{q^2 y}</math>	<math>~-\frac{2\lambda_2}{p^2 z}</math>	<math>~\frac{x q^2 y p^2 z}{\mathcal{D}} \biggl(\frac{1}{x}\biggr)</math>	<math>~ \frac{x q^2 y p^2 z}{\mathcal{D}} \biggl(\frac{1}{q^2y}\biggr)</math>	<math>~\frac{x q^2 y p^2 z}{\mathcal{D}} \biggl(-\frac{2}{p^2z}\biggr)</math>
<math>~3</math>	<math>~z</math>	<math>~\frac{\ell_{3D}}{\mathcal{D}}</math>	<math>~-x(2 q^4y^2 + p^4z^2)</math>	<math>~q^2 y(p^4z^2 + 2x^2 )</math>	<math>~p^2z( x^2 - q^4y^2 )</math>	<math>~-\frac{\ell_{3D}}{\mathcal{D}}\biggl[ x(2 q^4y^2 + p^4z^2 ) \biggl]</math>	<math>~\frac{\ell_{3D}}{\mathcal{D}}\biggl[ q^2 y(p^4z^2 + 2x^2 ) \biggl]</math>	<math>~\frac{\ell_{3D}}{\mathcal{D}}\biggl[ p^2z( x^2 - q^4y^2 ) \biggl]</math>

<math>~\vec{g}_1</math>

<math>~=</math>

<math>~ \hat\imath h_1\biggl(\frac{\partial\lambda_1}{\partial x} \biggr) </math>

See Also

© 2014 - 2021 by Joel E. Tohline |   H_Book Home   |   YouTube   | Appendices: | Equations | Variables | References | Ramblings | Images | myphys.lsu | ADS | Recommended citation:   Tohline, Joel E. (2021), The Structure, Stability, & Dynamics of Self-Gravitating Fluids, a (MediaWiki-based) Vistrails.org publication, https://www.vistrails.org/index.php/User:Tohline/citation