Ehrenfest equations: Difference between revisions
en>Saehry m Moved from "Thermodynamics" category to subcategory. |
|||
Line 1: | Line 1: | ||
A system of '''skew coordinates''' is a [[coordinate system]] where the coordinate surfaces are not [[orthogonal]],<ref>[http://mathworld.wolfram.com/SkewCoordinateSystem.html Skew Coordinate System] at [[Mathworld]]</ref> in contrast to [[orthogonal coordinates]]. | |||
Skew coordinates tend to be more complicated to work with compared to orthogonal coordinates since the [[metric tensor]] will have nonzero off-diagonal components, preventing many drastic simplifications in formulas for [[tensor algebra]] and [[tensor calculus]]. The nonzero off-diagonal components of the metric tensor are a direct result of the non-orthogonality of the basis vectors of the coordinates, since by definition:<ref name=p13> | |||
{{cite book | |||
| last = Lebedev | |||
| first = Leonid P. | |||
| authorlink = | |||
| coauthors = | |||
| title = Tensor Analysis | |||
| publisher = World Scientific | |||
| year = 2003 | |||
| location = | |||
| pages = 13 | |||
| url = | |||
| doi = | |||
| id = | |||
| isbn = 981-238-360-3}} | |||
</ref> | |||
:<math>g_{i j} = \mathbf e_i \cdot \mathbf e_j</math> | |||
where <math>g_{i j}</math> is the metric tensor and <math>\mathbf e_i</math> the (covariant) [[basis vector]]s. | |||
These coordinate systems can be useful if the geometry of a problem fits well into a skewed system. For example, solving [[Laplace's equation]] in a [[parallelogram]] will be easiest when done in appropriately skewed coordinates. | |||
==Cartesian coordinates with one skewed axis== | |||
[[Image:SkewCartesianSystem.svg|thumb|right|A Cartesian coordinate system where the ''x'' axis has been bent toward the ''z'' axis.]] | |||
The simplest 3D case of a skew coordinate system is a [[Cartesian coordinates|Cartesian]] one where one of the axes (say the ''x'' axis) has been bent by some angle <math>\phi</math>, staying orthogonal to one of the remaining two axes. For this example, the ''x'' axis of a Cartesian coordinate has been bent toward the ''z'' axis by <math>\phi</math>, remaining orthogonal to the ''y'' axis. | |||
===Algebra and useful quantities=== | |||
Let <math>\mathbf e_1</math>, <math>\mathbf e_2</math>, and <math>\mathbf e_3</math> respectively be unit vectors along the <math>x</math>, <math>y</math>, and <math>z</math> axes. These represent the [[Covariance and contravariance of vectors|covariant]] basis; computing their dot products gives the following components of the [[metric tensor]]: | |||
:<math>g_{11} = g_{22} = g_{33} = 1 \quad ; \quad g_{12} = g_{23} = 0 \quad ; \quad g_{13} = \cos\left(\frac \pi 2 - \phi\right) = \sin(\phi)</math> | |||
:<math>\sqrt{g} = \mathbf e_1 \cdot (\mathbf e_2 \times \mathbf e_3) = \cos(\phi)</math> | |||
which are quantities that will be useful later on. | |||
The contravariant basis is given by<ref name=p13/> | |||
:<math>\mathbf e^1 = \frac{\mathbf e_2 \times \mathbf e_3}{\sqrt{g}} = \frac{\mathbf e_2 \times \mathbf e_3}{\cos(\phi)}</math> | |||
:<math>\mathbf e^2 = \frac{\mathbf e_3 \times \mathbf e_1}{\sqrt{g}} = \mathbf e_2</math> | |||
:<math>\mathbf e^3 = \frac{\mathbf e_1 \times \mathbf e_2}{\sqrt{g}} = \frac{\mathbf e_1 \times \mathbf e_2}{\cos(\phi)}</math> | |||
The contravariant basis isn't a very convenient one to use, however it shows up in definitions so must be considered. We'll favor writing quantities with respect to the covariant basis. | |||
Since the basis vectors are all constant, vector addition and subtraction will simply be familiar component-wise adding and subtraction. Now, let | |||
:<math>\mathbf a = \sum_i a^i \mathbf e_i \quad \mbox{and} \quad \mathbf b = \sum_i b^i \mathbf e_i</math> | |||
where the sums indicate summation over all values of the index (in this case, ''i'' = 1, 2, 3). The [[Covariance and contravariance of vectors|contravariant and covariant]] components of these vectors may be related by | |||
:<math>a^i = \sum_j a_j g^{ij}\,</math> | |||
so that, explicitly, | |||
:<math>a^1 = \frac{a_1 - \sin(\phi) a_3}{\cos(\phi)^2},\,</math> | |||
:<math>a^2 = a_2,\,</math> | |||
:<math>a^3 = \frac{-\sin(\phi) a_1 + a_3}{\cos(\phi)^2}.\,</math> | |||
The [[dot product]] in terms of contravariant components is then | |||
:<math>\mathbf a \cdot \mathbf b = \sum_i a^i b_i = a^1 b^1 + a^2 b^2 + a^3 b^3 + \sin(\phi) (a^1 b^3 + a^3 b^1).</math> | |||
===Calculus=== | |||
By definition,<ref>{{cite book | |||
| last = Lebedev | |||
| first = Leonid P. | |||
| authorlink = | |||
| coauthors = | |||
| title = Tensor Analysis | |||
| publisher = World Scientific | |||
| year = 2003 | |||
| location = | |||
| pages = 63 | |||
| url = | |||
| doi = | |||
| id = | |||
| isbn = 981-238-360-3}} | |||
</ref> the [[gradient]] of a scalar function ''f'' is | |||
:<math>\nabla f = \sum_i \mathbf e^i \frac{\partial f}{\partial q^i} = \frac{\partial f}{\partial x} \mathbf e^1 + \frac{\partial f}{\partial y} \mathbf e^2 + \frac{\partial f}{\partial z} \mathbf e^3</math> | |||
where <math>q_i</math> are the coordinates ''x'', ''y'', ''z'' indexed. Recognizing this as a vector written in terms of the contravariant basis, it may be rewritten: | |||
:<math>\nabla f = | |||
\frac{\frac{\partial f}{\partial x} - \sin(\phi) \frac{\partial f}{\partial z}}{\cos(\phi)^2} \mathbf e_1 + | |||
\frac{\partial f}{\partial y} \mathbf e_2 + | |||
\frac{-\sin(\phi) \frac{\partial f}{\partial x} + \frac{\partial f}{\partial z}}{\cos(\phi)^2} \mathbf e_3.</math> | |||
The [[divergence]] of a vector <math>\mathbf a</math> is | |||
:<math>\nabla \cdot \mathbf a = \frac{1}{\sqrt{g}} \sum_i \frac{\partial}{\partial q^i}\left(\sqrt{g} a^i\right) = \frac{\partial a^1}{\partial x} + \frac{\partial a^2}{\partial y} + \frac{\partial a^3}{\partial z}.</math> | |||
and of a tensor <math>\mathbf A</math> | |||
:<math>\nabla \cdot \mathbf A = \frac{1}{\sqrt{g}} \sum_{i, j} \frac{\partial}{\partial q^i}\left(\sqrt{g} a^{ij} \mathbf e_j\right) = | |||
\sum_{i, j} \mathbf e_j \frac{\partial a^{ij}}{\partial q^i}.</math> | |||
The [[Laplacian]] of ''f'' is | |||
:<math>\nabla^2 f = \nabla \cdot \nabla f = | |||
\frac{1}{\cos(\phi)^2}\left(\frac{\partial^2 f}{\partial x^2} + \frac{\partial^2 f}{\partial z^2} - 2 \sin(\phi) \frac{\partial^2 f}{\partial x \partial z}\right) + \frac{\partial^2 f}{\partial y^2}</math> | |||
and, since the covariant basis is normal and constant, the [[vector Laplacian]] is the same as the componentwise Laplacian of a vector written in terms of the covariant basis. | |||
While both the dot product and gradient are somewhat messy in that they have extra terms (compared to a Cartesian system) the [[advection|advection operator]] which combines a dot product with a gradient turns out very simple: | |||
:<math>(\mathbf a \cdot \nabla) = \left(\sum_i a^i e_i\right) \cdot \left(\sum_i \frac{\partial}{\partial q^i} \mathbf e^i\right) = \left(\sum_i a^i \frac{\partial}{\partial q^i}\right) | |||
</math> | |||
which may be applied to both scalar functions and vector functions, componentwise when expressed in the covariant basis. | |||
Finally, the [[curl (mathematics)|curl]] of a vector is | |||
:<math>\nabla \times \mathbf a = \sum_{i, j, k} \mathbf e_k \epsilon^{ijk} \frac{\partial a_i}{\partial q^i} = </math> | |||
::<math>\frac{1}{\cos(\phi)}\left( | |||
\left(\sin(\phi) \frac{\partial a^1}{\partial y} + \frac{\partial a^3}{\partial y} - \frac{\partial a^2}{\partial z}\right) \mathbf e_1 + | |||
\left(\frac{\partial a^1}{\partial z} + \sin(\phi) \left(\frac{\partial a^3}{\partial z} - \frac{\partial a^1}{\partial x}\right) - \frac{\partial a^3}{\partial x}\right) \mathbf e_2 + | |||
\left(\frac{\partial a^2}{\partial x} - \frac{\partial a^1}{\partial y} - \sin(\phi) \frac{\partial a^3}{\partial y}\right) \mathbf e_3 | |||
\right).</math> | |||
==References== | |||
{{Reflist}} | |||
[[Category:Coordinate systems]] |
Latest revision as of 11:21, 18 May 2013
A system of skew coordinates is a coordinate system where the coordinate surfaces are not orthogonal,[1] in contrast to orthogonal coordinates.
Skew coordinates tend to be more complicated to work with compared to orthogonal coordinates since the metric tensor will have nonzero off-diagonal components, preventing many drastic simplifications in formulas for tensor algebra and tensor calculus. The nonzero off-diagonal components of the metric tensor are a direct result of the non-orthogonality of the basis vectors of the coordinates, since by definition:[2]
where is the metric tensor and the (covariant) basis vectors.
These coordinate systems can be useful if the geometry of a problem fits well into a skewed system. For example, solving Laplace's equation in a parallelogram will be easiest when done in appropriately skewed coordinates.
Cartesian coordinates with one skewed axis
The simplest 3D case of a skew coordinate system is a Cartesian one where one of the axes (say the x axis) has been bent by some angle , staying orthogonal to one of the remaining two axes. For this example, the x axis of a Cartesian coordinate has been bent toward the z axis by , remaining orthogonal to the y axis.
Algebra and useful quantities
Let , , and respectively be unit vectors along the , , and axes. These represent the covariant basis; computing their dot products gives the following components of the metric tensor:
which are quantities that will be useful later on.
The contravariant basis is given by[2]
The contravariant basis isn't a very convenient one to use, however it shows up in definitions so must be considered. We'll favor writing quantities with respect to the covariant basis.
Since the basis vectors are all constant, vector addition and subtraction will simply be familiar component-wise adding and subtraction. Now, let
where the sums indicate summation over all values of the index (in this case, i = 1, 2, 3). The contravariant and covariant components of these vectors may be related by
so that, explicitly,
The dot product in terms of contravariant components is then
Calculus
By definition,[3] the gradient of a scalar function f is
where are the coordinates x, y, z indexed. Recognizing this as a vector written in terms of the contravariant basis, it may be rewritten:
The divergence of a vector is
The Laplacian of f is
and, since the covariant basis is normal and constant, the vector Laplacian is the same as the componentwise Laplacian of a vector written in terms of the covariant basis.
While both the dot product and gradient are somewhat messy in that they have extra terms (compared to a Cartesian system) the advection operator which combines a dot product with a gradient turns out very simple:
which may be applied to both scalar functions and vector functions, componentwise when expressed in the covariant basis.
Finally, the curl of a vector is
References
43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.
- ↑ Skew Coordinate System at Mathworld
- ↑ 2.0 2.1
20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
My blog: http://www.primaboinca.com/view_profile.php?userid=5889534 - ↑ 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
My blog: http://www.primaboinca.com/view_profile.php?userid=5889534