Криволинейные координаты могут быть сформулированы в тензорном исчислении с важными приложениями в физика и инженерия, особенно для описания перемещения физических величин и деформации материи в механике жидкости и механике сплошных сред.
Элементарная векторная и тензорная алгебра в криволинейных координатах используется в некоторой старой научной литературе по механике и физика и может быть незаменим для понимания работ начала и середины 1900-х годов, например текста Грина и Зерны. Некоторые полезные rel В этом разделе приведены элементы алгебры векторов и тензоров второго порядка в криволинейных координатах. Обозначения и содержание в основном взяты из Огдена, Нагди, Симмондса, Грина и Зерны, Басара и Вейхерта и Чиарлета.
Рассмотрим две системы координат с координатными переменными (Z 1, Z 2, Z 3) {\ displaystyle (Z ^ {1}, Z ^ {2}, Z ^ {3})}и (Z 1 ´, Z 2 ´, Z 3 ´) {\ displaystyle (Z ^ {\ Acustric {1}}, Z ^ {\ sharp {2}}, Z ^ {\ sharp {3}})}, который мы будем представлять вкратце, просто Z i {\ displaystyle Z ^ {i}}и Z i ´ {\ displaystyle Z ^ {\ строго {i}}}соответственно и всегда предполагаем, что наш индекс i {\ displaystyle i}проходит от 1 до 3. Мы будем предполагать,что эти системы координат встроены в трехмерное евклидово пространство. Могут использоваться координаты Z i {\ displaystyle Z ^ {i}}и Z i ´ {\ displaystyle Z ^ {\ строго {i}}}. чтобы объяснить друг друга, потому что, перемещаясь по линии координат в одной системе координат, мы можем использовать другую для описания нашего положения. Таким образом, координаты Z i {\ displaystyle Z ^ {i}}и Z i ´ {\ displaystyle Z ^ {\ строго {i}}}являются функциями друг друга
Z я = fi (Z 1 ´, Z 2 ´, Z 3 ´) {\ displaystyle Z ^ {i} = f ^ {i} (Z ^ {\ строго {1}}, Z ^ {\ острый {2}}, Z ^ {\ острый {3}})}для i = 1, 2, 3 {\ displaystyle i = 1,2,3}
, который можно записать как
Z i = Z i (Z 1 ´, Z 2 ´, Z 3 ´) = Z i (Z i ´) {\ displaystyle Z ^ {i} = Z ^ {i} (Z ^ {\ строго {1}}, Z ^ {\ строго {2}}, Z ^ {\ строго {3}}) = Z ^ {i} (Z ^ {\ строго {i}})}для i ´, i = 1, 2, 3 {\ displaystyle {\ строго {i}}, i = 1,2,3}
Эти три уравнения вместе также называются координатами преобразование из Z i ´ {\ displaystyle Z ^ {\ строго {i}}}в Z i {\ displaystyle Z ^ {i}}. мы обозначим это преобразование как T {\ displaystyle T}. Поэтому мы представим преобразование из системы координат с координатными переменными Z i ´ {\ displaystyle Z ^ {\ строго {i}}}в систему координат с координатами Z i { \ displaystyle Z ^ {i}}как:
Z = T (z ´) {\ displaystyle Z = T ({\ строго {z}})}
Аналогичным образом мы можем представить Z я ´ {\ displaystyle Z ^ {\ строго {i}}}как функция от Z i {\ displaystyle Z ^ {i}}следующим образом :
Z я ´ = gi ´ (Z 1, Z 2, Z 3) {\ displaystyle Z ^ {\ строго {i}} = g ^ {\ строго {i}} (Z ^ {1}, Z ^ {2}, Z ^ {3})}для i ´ = 1, 2, 3 {\ displaystyle {\ строго {i}} = 1,2,3}
аналогичным образом мы можем записать свободные уравнения болеекомпактно как
Z i ´ = Z i ´ (Z 1, Z 2, Z 3) = Z i ´ (Z i) {\ displaystyle Z ^ {\ строго {i}} = Z ^ {\ строго {i}} (Z ^ {1}, Z ^ {2}, Z ^ {3}) = Z ^ {\ строго {i}} (Z ^ {i})}для i ´, i = 1, 2, 3 {\ displaystyle {\ строго {i}}, i = 1,2,3}
Эти три уравнения вместе также называются преобразованием координат из Z i {\ displaystyle Z ^ {i}}до Z i ´ curvilinear coordinatesLet (b1, b2, b3) be an arbitrary basis for three-dimensional Euclidean space. In general, the basis vectors are neither unit vectors nor mutually orthogonal. However, they are required to be linearly independent. Then a vector vcan be expressed asv = v k b k {\displaystyle \mathbf {v} =v^{k}\,\mathbf {b} _{k}}The components v are the contravariantcomponents of the vector v. The reciprocal basis(b, b, b) is defined by the relationb i ⋅ b j = δ j i {\displaystyle \mathbf {b} ^{i}\cdot \mathbf {b} _{j}=\delta _{j}^{i}}where δjis the Kronecker delta. The vector vcan also be expressed in terms of the reciprocal basis:v = v k b k {\displaystyle \mathbf {v} =v_{k}~\mathbf {b} ^{k}}The components vkare the covariantcomponents of the vector v {\displaystyle \mathbf {v} }.Second-order tensors in curvilinear coordinatesA second-order tensor can be expressed asS = S i j b i ⊗ b j = S j i b i ⊗ b j = S i j b i ⊗ b j = S i j b i ⊗ b j {\displaystyle {\boldsymbol {S}}=S^{ij}~\mathbf {b} _{i}\otimes \mathbf {b} _{j}=S_{~j}^{i}~\mathbf {b} _{i}\otimes \mathbf {b} ^{j}=S_{i}^{~j}~\mathbf {b} ^{i}\otimes \mathbf {b} _{j}=S_{ij}~\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}}The components S are called the contravariantcomponents, Sjthe mixed right-covariantcomponents, Sithe mixed left-covariantcomponents, and Sijthe covariantcomponents of the second-order tensor.Metric tensor and relations between componentsThe quantities gij, g are defined asg i j = b i ⋅ b j = g j i ; g i j = b i ⋅ b j = g j i {\displaystyle g_{ij}=\mathbf {b} _{i}\cdot \mathbf {b} _{j}=g_{ji}~;~~g^{ij}=\mathbf {b} ^{i}\cdot \mathbf {b} ^{j}=g^{ji}}From the above equations we havev i = g i k v k ; v i = g i k v k ; b i = g i j b j ; b i = g i j b j {\displaystyle v^{i}=g^{ik}~v_{k}~;~~v_{i}=g_{ik}~v^{k}~;~~\mathbf {b} ^{i}=g^{ij}~\mathbf {b} _{j}~;~~\mathbf {b} _{i}=g_{ij}~\mathbf {b} ^{j}}The components of a vector are related byv ⋅ b i = v k b k ⋅ b i = v k δ k i = v i {\displaystyle \mathbf {v} \cdot \mathbf {b} ^{i}=v^{k}~\mathbf {b} _{k}\cdot \mathbf {b} ^{i}=v^{k}~\delta _{k}^{i}=v^{i}}v ⋅ b i = v k b k ⋅ b i = v k δ i k = v i {\displaystyle \mathbf {v} \cdot \mathbf {b} _{i}=v_{k}~\mathbf {b} ^{k}\cdot \mathbf {b} _{i}=v_{k}~\delta _{i}^{k}=v_{i}}Also,v ⋅ b i = v k b k ⋅ b i = g k i v k {\displaystyle \mathbf {v} \cdot \mathbf {b} _{i}=v^{k}~\mathbf {b} _{k}\cdot \mathbf {b} _{i}=g_{ki}~v^{k}}v ⋅ b i = v k b k ⋅ b i = g k i v k {\displaystyle \mathbf {v} \cdot \mathbf {b} ^{i}=v_{k}~\mathbf {b} ^{k}\cdot \mathbf {b} ^{i}=g^{ki}~v_{k}}The components of the second-order tensor are related byS i j = g i k S k j = g j k S k i = g i k g j l S k l {\displaystyleS^{ij}=g^{ik}~S_{k}^{~j}=g^{jk}~S_{~k}^{i}=g^{ik}~g^{jl}~S_{kl}}The alternating tensorIn an orthonormal right-handed basis, the third-order alternating tensor is defined asE = ε i j k e i ⊗ e j ⊗ e k {\displaystyle {\boldsymbol {\mathcal {E}}}=\varepsilon _{ijk}~\mathbf {e} ^{i}\otimes \mathbf {e} ^{j}\otimes \mathbf {e} ^{k}}In a general curvilinear basis the same tensor may be expressed asE = E i j k b i ⊗ b j ⊗ b k = E i j k b i ⊗ b j ⊗ b k {\displaystyle {\boldsymbol {\mathcal {E}}}={\mathcal {E}}_{ijk}~\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}\otimes \mathbf {b} ^{k}={\mathcal {E}}^{ijk}~\mathbf {b} _{i}\otimes \mathbf {b} _{j}\otimes \mathbf {b} _{k}}It can be shown thatE i j k = [ b i, b j, b k ] = ( b i × b j) ⋅ b k ; E i j k = [ b i, b j, b k ] {\displaystyle {\mathcal {E}}_{ijk}=\left[\mathbf {b} _{i},\mathbf {b} _{j},\mathbf {b} _{k}\right]=(\mathbf {b} _{i}\times \mathbf {b} _{j})\cdot \mathbf {b} _{k}~;~~{\mathcal {E}}^{ijk}=\left[\mathbf {b} ^{i},\mathbf {b} ^{j},\mathbf {b} ^{k}\right]}Now,b i × b j = J ε i j p b p = g ε i j p b p {\displaystyle \mathbf {b} _{i}\times \mathbf {b} _{j}=J~\varepsilon _{ijp}~\mathbf {b} ^{p}={\sqrt {g}}~\varepsilon _{ijp}~\mathbf {b} ^{p}}Hence,E i j k = J ε i j k = g ε i j k {\displaystyle {\mathcal {E}}_{ijk}=J~\varepsilon _{ijk}={\sqrt {g}}~\varepsilon _{ijk}}Similarly, we can show thatE i j k = 1 J ε i j k = 1 g ε i j k {\displaystyle {\mathcal {E}}^{ijk}={\cfrac {1}{J}}~\varepsilon ^{ijk}={\cfrac {1}{\sqrt {g}}}~\varepsilon ^{ijk}}Vector operationsIdentity mapThe identity map Idefined by I ⋅ v = v {\displaystyle \mathbf {I} \cdot \mathbf {v} =\mathbf {v} }can be shown to be:I = g i j b i ⊗ b j = g i j b i ⊗ b j = b i ⊗ b i = b i ⊗ b i {\displaystyle \mathbf {I} =g^{ij}\mathbf {b} _{i}\otimes \mathbf {b} _{j}=g_{ij}\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}=\mathbf {b} _{i}\otimes \mathbf {b} ^{i}=\mathbf {b} ^{i}\otimes \mathbf {b} _{i}}Scalar (dot) productThe scalar product of two vectors in curvilinear coordinates isu ⋅ v = u i v i = u i v i = g i j u i v j = g i j u i v j {\displaystyle \mathbf {u} \cdot \mathbf {v} =u^{i}v_{i}=u_{i}v^{i}=g_{ij}u^{i}v^{j}=g^{ij}u_{i}v_{j}}Vector (cross) productThe cross product of two vectors is given by:u × v = ε i j k u j v k e i {\displaystyle \mathbf {u} \times \mathbf {v} =\varepsilon _{ijk}u_{j}v_{k}\mathbf {e} _{i}}where εijkis the permutation symbol and eiis a Cartesian basis vector. In curvilinear coordinates, the equivalent expression is:u × v = [ ( b m × b n) ⋅ b s ] u m v n b s = E s m n u m v n b s {\displaystyle \mathbf {u} \times \mathbf {v} =[(\mathbf {b} _{m}\times \mathbf {b} _{n})\cdot \mathbf {b} _{s}]u^{m}v^{n}\mathbf {b} ^{s}={\mathcal {E}}_{smn}u^{m}v^{n}\mathbf {b} ^{s}}where E i j k{\displaystyle {\mathcal {E}}_{ijk}}is the third-order alternating tensor. The cross product of two vectors is given by:u × v = ε i j k u ^ j v ^ k e i {\displaystyle \mathbf {u} \times \mathbf {v} =\varepsilon _{ijk}{\hat {u}}_{j}{\hat {v}}_{k}\mathbf {e} _{i}}where εijkis the permutation symbol and e i {\displaystyle \mathbf {e} _{i}}is a Cartesian basis vector. Therefore,e p × e q = ε i p q e i {\displaystyle \mathbf {e} _{p}\times \mathbf {e} _{q}=\varepsilon _{ipq}\mathbf {e} _{i}}andb m × b n = ∂ x ∂ q m × ∂ x ∂ q n = ∂ ( x p e p) ∂ q m × ∂ ( x q e q) ∂ q n = ∂ x p ∂ q m ∂ x q ∂ q n e p × e q = ε i p q ∂ x p ∂ q m ∂ x q ∂ q n e i. {\displaystyle \mathbf {b} _{m}\times \mathbf {b} _{n}={\frac {\partial \mathbf {x} }{\partial q^{m}}}\times {\frac {\partial \mathbf {x} }{\partial q^{n}}}={\frac {\partial (x_{p}\mathbf {e} _{p})}{\partial q^{m}}}\times {\frac {\partial (x_{q}\mathbf {e} _{q})}{\partial q^{n}}}={\frac {\partial x_{p}}{\partial q^{m}}}{\frac {\partial x_{q}}{\partial q^{n}}}\mathbf {e} _{p}\times \mathbf {e} _{q}=\varepsilon _{ipq}{\frac {\partial x_{p}}{\partial q^{m}}}{\frac {\partial x_{q}}{\partial q^{n}}}\mathbf {e} _{i}.}Hence,( b m × b n) ⋅ b s = ε i p q ∂ x p ∂ q m ∂ x q ∂ q n ∂ x i ∂ q s {\displaystyle (\mathbf {b} _{m}\times \mathbf {b} _{n})\cdot \mathbf {b} _{s}=\varepsilon _{ipq}{\frac {\partial x_{p}}{\partial q^{m}}}{\frac {\partial x_{q}}{\partial q^{n}}}{\frac {\partial x_{i}}{\partial q^{s}}}}Returning to the vector product andusing the relations:u ^ j = ∂ x j ∂ q m u m, v ^ k = ∂ x k ∂ q n v n, e i = ∂ x i ∂ q s b s, {\displaystyle {\hat {u}}_{j}={\frac {\partial x_{j}}{\partial q^{m}}}u^{m},\quad {\hat {v}}_{k}={\frac {\partial x_{k}}{\partial q^{n}}}v^{n},\quad \mathbf {e} _{i}={\frac {\partial x_{i}}{\partial q^{s}}}\mathbf {b} ^{s},}gives us:u × v = ε i j k u ^ j v ^ k e i = ε i j k ∂ x j ∂ q m ∂ x k ∂ q n ∂ x i ∂ q s u m v n b s = [ ( b m × b n) ⋅ b s ] u m v n b s = E s m n u m v n b s {\displaystyle \mathbf {u} \times \mathbf {v} =\varepsilon _{ijk}{\hat {u}}_{j}{\hat {v}}_{k}\mathbf {e} _{i}=\varepsilon _{ijk}{\frac {\partial x_{j}}{\partial q^{m}}}{\frac {\partial x_{k}}{\partial q^{n}}}{\frac {\partial x_{i}}{\partial q^{s}}}u^{m}v^{n}\mathbf {b} ^{s}=[(\mathbf {b} _{m}\times \mathbf {b} _{n})\cdot \mathbf {b} _{s}]u^{m}v^{n}\mathbf {b} ^{s}={\mathcal {E}}_{smn}u^{m}v^{n}\mathbf {b} ^{s}}Tensor operationsIdentity mapThe identity map I {\displaystyle {\mathsf {I}}}defined by I ⋅ v = v {\displaystyle {\mathsf {I}}\cdot \mathbf {v} =\mathbf {v} }can be shown to beI = g i j b i ⊗ b j = g i j b i ⊗ b j = b i ⊗ b i = b i ⊗ b i {\displaystyle {\mathsf {I}}=g^{ij}\mathbf {b} _{i}\otimes \mathbf {b} _{j}=g_{ij}\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}=\mathbf {b} _{i}\otimes \mathbf {b} ^{i}=\mathbf {b} ^{i}\otimes\mathbf {b} _{i}}Action of a second-order tensor on a vectorThe action v = S ⋅ u {\displaystyle \mathbf {v} ={\boldsymbol {S}}\cdot \mathbf {u} }can be expressed in curvilinear coordinates asv i b i = S i j u j b i = S j i u j b i ; v i b i = S i j u i b i = S i j u j b i {\displaystyle v^{i}\mathbf {b} _{i}=S^{ij}u_{j}\mathbf {b} _{i}=S_{j}^{i}u^{j}\mathbf {b} _{i};\qquad v_{i}\mathbf {b} ^{i}=S_{ij}u^{i}\mathbf {b} ^{i}=S_{i}^{j}u_{j}\mathbf {b} ^{i}}Inner product of two second-order tensorsThe inner product of two second-order tEnsors U = S ⋅ T {\ displaystyle {\ boldsymbol {U}} = {\ boldsymbol {S}} \ cdot {\ boldsymbol {T}}}может быть выражено в криволинейных координатах какU ijbi ⊗ bj = S ik T. Дж К б я ⊗ б J знак равно S я. к T kjbi ⊗ bj {\ displaystyle U_ {ij} \ mathbf {b} ^ {i} \ otimes \ mathbf {b} ^ {j} = S_ {ik} T _ {. j} ^ {k} \ mathbf {b } ^ {i} \ otimes \ mathbf {b} ^ {j} = S_ {i} ^ {. k} T_ {kj} \ mathbf {b} ^ {i} \ otimes \ mathbf {b} ^ {j} }В качестве альтернативыU = S ij T. н м г дж м б я ⊗ б н знак равно S. м я т. nmbi ⊗ bn = S ij T jnbi ⊗ bn {\ displaystyle {\ boldsymbol {U}} = S ^ {ij} T _ {. n} ^ {m} g_ {jm} \ mathbf {b} _ {i} \ otimes \ mathbf {b} ^ {n} = S _ {. m} ^ {i} T _ {. n} ^ {m} \ mathbf {b} _ {i} \ otimes \ mathbf {b} ^ {n} = S ^ {ij} T_ {jn} \ mathbf {b} _ {i} \ otimes \ mathbf {b} ^ {n}}Определитель тензора второго порядкаЕсли S {\ displaystyle {\ boldsymbol {S}}}- тензор второго порядка, тогда определитель определяется соотношением[S ⋅ u, S ⋅ v, S ⋅ w ] = det S [u, v, w] {\ displaystyle \ left [{\ boldsymbol {S}} \ cdot \ mathbf {u}, {\ boldsymbol {S}} \ cdot \ mathbf {v}, {\ boldsymbol {S}} \ cdot \ mathbf {w} \ right] = \ det {\ boldsymbol {S}} \ left [\ mathbf {u}, \ mathbf {v}, \ mathbf {w} \ right]}где u, v, w {\ displaystyle \ mathbf {u}, \ mathbf {v}, \ mathbf {w}}- произвольные векторы, а[u, v, w ]: = u ⋅ (v × w). {\ displaystyle \ left [\ mathbf {u}, \ mathbf {v}, \ mathbf {w} \ right]: = \ mathbf {u} \ cdot (\ mathbf {v} \ times \ mathbf {w}). }Отношения между криволинейными и декартовыми базисными векторамиПусть (e1, e2, e3) - обычные декартовы базисные векторы для интересующего евклидова пространства, и пустьbi = F ⋅ ei {\ displaystyle \ mathbf {b} _ {i} = {\ boldsymbol {F}} \ cdot \ mathbf {e} _ {i}}, где Fi- тензор преобразования второго порядка, который отображает eiв bi. Тогдаb i ⊗ e i = (F ⋅ e i)⊗ e i = F ⋅ (e i ⊗ e i) = F. {\ displaystyle \ mathbf {b} _ {i} \ otimes \ mathbf {e} _ {i} = ({\ boldsymbol {F}} \ cdot \ mathbf {e} _ {i}) \ otimes \ mathbf {e } _ {i} = {\ boldsymbol {F}} \ cdot (\ mathbf {e} _ {i} \ otimes \ mathbf {e} _ {i}) = {\ boldsymbol {F}} ~.}Из этого соотношения мы можем показать, чтоbi = F - T ⋅ ei; g i j = [F - 1 ⋅ F - T] i j; gij = [gij] - 1 = [FT ⋅ F] ij {\ displaystyle \ mathbf {b} ^ {i} = {\ boldsymbol {F}} ^ {- {\ rm {T}}} \ cdot \ mathbf { e} ^ {i} ~; ~~ g ^ {ij} = [{\ boldsymbol {F}} ^ {- {\ rm {1}}} \ cdot {\ boldsymbol {F}} ^ {- {\ rm {T}}}] _ {ij} ~; ~~ g_ {ij} = [g ^ {ij}] ^ {- 1} = [{\ bol dsymbol {F}} ^ {\ rm {T}} \ cdot {\ boldsymbol {F}}] _ {ij}}Пусть J: = det F {\ displaystyle J: = \ det {\ boldsymbol {F}}}будет Якобиан преобразования. Тогда, согласно определению определителя,[b 1, b 2, b 3] = det F [e 1, e 2, e 3]. {\ displaystyle \ left [\ mathbf {b} _ {1}, \ mathbf {b} _ {2}, \ mathbf {b} _ {3} \ right] = \ det {\ boldsymbol {F}} \ left [\ mathbf {e} _ {1}, \ mathbf {e} _ {2}, \ mathbf {e} _ {3} \ right] ~.}Поскольку[e 1, e 2, е 3] = 1 {\ displaystyle \ left [\ mathbf {e} _ {1}, \ mathbf {e} _ {2}, \ mathbf {e} _ {3} \ right] = 1}мы иметьJ = det F = [b 1, b 2, b 3] = b 1 ⋅ (b 2 × b 3) {\ displaystyle J = \ det {\ boldsymbol {F}} = \ left [\ mathbf {b} _ {1}, \ mathbf {b} _ {2}, \ mathbf {b} _ {3} \ right] = \ mathbf {b} _ {1} \ cdot (\ mathbf {b} _ { 2} \ times \ mathbf {b} _ {3})}С помощью приведенных выше соотношений можно получить ряд интересных результатов. Сначала рассмотримg: = det [gij] {\ displaystyle g: = \ det [g_ {ij}]}Затемg = det [FT] ⋅ det [ F] = J ⋅ J = J 2 {\ displaystyle g = \ det [{\ boldsymbol {F}} ^ {\ rm {T}}] \ cdot \ det [{\ boldsymbol {F}}] = J \ cdot J = J ^ {2}}Аналогичным образом мы можем показать, чтоdet [gij] = 1 J 2 {\ displaystyle \ det [g ^ {ij}] = {\ cfrac {1} {J ^ {2}}}}Следовательно, используя тот факт, что [gij] = [gij] - 1 {\ displaystyle [g ^ {ij}] = [g_ {ij}] ^ {- 1}},∂ g ∂ gij знак равно 2 J ∂ J ∂ gij = ggij {\ displaystyle {\ cfrac {\ partial g} {\ partial g_ {ij}}} = 2 ~ J ~ {\ cfrac {\ partial J} {\ partial g_ {ij }}} = g ~ g ^ {ij}}Еще одно интересное соотношение выводится ниже. Напомним, чтоbi ⋅ bj = δ ji ⇒ b 1 ⋅ b 1 = 1, b 1 ⋅ b 2 bf {e} _ {i} = \ mathbf {b} _ {j} ~ {\ cfrac {\ partial q ^ {j}} {\ partial x_ {i}}}}Из этих результатов имеемek ⋅ bi = ∂ xk ∂ qi ⇒ ∂ xk ∂ qibi = ek ⋅ (bi ⊗ bi) = эк {\ displaystyle \ mathbf {e} ^ {k} \ cdot \ mathbf {b} _ {i} = {\ frac {\ partial x_ {k}} {\ partial q ^ {i}}} \ quad \ Rightarrow \ quad {\ frac {\ partial x_ {k}} {\ partial q ^ {i}}} ~ \ mathbf {b} ^ {i} = \ mathbf {e} ^ {k} \ cdot (\ mathbf {b } _ {i} \ otimes \ mathbf {b} ^ {i}) = \ mathbf {e} ^ {k}}иbk = ∂ qk ∂ xiei {\ displaystyle \ mathbf {b} ^ {k} = {\ frac {\ partial q ^ {k}} {\ partial x_ {i}}} ~ \ mathbf {e} ^ {i}}Векторное и тензорное исчисление в трехмерных криволинейных координатахПримечание: соглашение Эйнштейна о суммировании суммирования по повторяющимся индексам используется ниже.Симмондс в своей книге по тензорному анализу цитирует Альберта Эйнштейна высказывание Магия этой теории вряд ли перестанет навязываться любому, кто действительно понимает od это; он представляет собой подлинный триумф метода абсолютного дифференциального исчисления, основанного Гауссом, Риманом, Риччи и Леви-Чивита. Векторное и тензорное исчисление в общих криволинейных координатах используется в тензорном анализе четырехмерных криволинейных многообразий в общей теории относительности, в механике криволинейных оболочки, при исследовании свойств инвариантности уравнений Максвелла, которые представляли интерес в метаматериалах и во многих других областях. В этом разделе приведены некоторые полезные соотношения в исчислении векторов и тензоров второго порядка в криволинейных координатах. ион. Обозначения и содержание в основном взяты из Огдена, Симмондса, Грина и Зерны, Басара и Вейхерта и Чиарлета.Основные определенияПусть положение точки в пространстве характеризуется тремя координатными переменными (q 1, q 2, q 3) {\ displaystyle (q ^ {1}, q ^ {2}, q ^ {3})}. Координатная кривая q представляет кривая, на которой q, q постоянны. Пусть x будет вектором положения точки относительно некоторого начала. Тогда, предполагая, что такое отображение и обратное ему отображение существуют и непрерывны, мы можем записатьx = φ (q 1, q 2, q 3); qi знак равно ψ я (Икс) знак равно [φ - 1 (Икс)] я {\ Displaystyle \ mathbf {x} = {\ boldsymbol {\ varphi}} (q ^ {1}, q ^ {2}, q ^ { 3}) ~; ~~ q ^ {i} = \ psi ^ {i} (\ mathbf {x}) = [{\ boldsymbol {\ varphi}} ^ {- 1} (\ mathbf {x})] ^ {i}}Поля ψ (x ) называются функциями криволинейных координат в криволинейной системе координат ψ(x) = φ(x). Координатные кривые q определяются однопараметрическим семейством функций, задаваемымxi (α) = φ (α, qj, qk), i ≠ j ≠ k {\ displaystyle \ mathbf {x} _ {i} (\ alpha) = {\ boldsymbol {\ varphi}} (\ alpha, q ^ {j}, q ^ {k}) ~, ~~ i \ neq j \ neq k}с фиксированными q, q.Касательный вектор к координатным кривымКасательный вектор к кривой xiв точке xi(α) (или к координатной кривой q я в точке x ) равноdxid α ≡ ∂ x ∂ qi {\ displaystyle {\ cfrac {\ rm {{d} \ mathbf {x} _ {i}} } {\ rm {{d} \ alpha}}} \ Equiv {\ cfrac {\ partial \ mathbf {x}} {\ partial q ^ {i}}}}ГрадиентСкалярное полеПусть f (x ) будет скалярным полем в пространстве. Тогдаf (x) = f [φ (q 1, q 2, q 3)] = f φ (q 1, q 2, q 3) {\ displaystyle f (\ mathbf {x}) = f [ {\ boldsymbol {\ varphi}} (q ^ {1}, q ^ {2}, q ^ {3})] = f _ {\ varphi} (q ^ {1}, q ^ {2}, q ^ { 3})}Градиент поля f определяется как[∇ f (x)] ⋅ c = dd α f (x + α c) | α знак равно 0 {\ Displaystyle [{\ boldsymbol {\ nabla}} е (\ mathbf {x})] \ cdot \ mat hbf {c} = {\ cfrac {\ rm {d}} {\ rm {{d} \ alpha}}} f (\ mathbf {x} + \ alpha \ mathbf {c}) {\ biggr |} _ { \ alpha = 0}}где c - произвольный постоянный вектор. Если мы определим компоненты c из c таковы, чтоqi + α ci = ψ i (x + α c) {\ displaystyle q ^ {i} + \ alpha ~ c ^ {i} = \ psi ^ {i} (\ mathbf {x} + \ alpha ~ \ mathbf {c})}, затем[∇ f (x)] ⋅ c = dd α f φ (q 1 + α c 1, q 2 + α c 2, q 3 + α c 3) | α знак равно 0 знак равно ∂ е φ ∂ qici = ∂ е ∂ qici {\ displaystyle [{\ boldsymbol {\ nabla}} f (\ mathbf {x})] \ cdot \ mathbf {c} = {\ cfrac {\ rm { d}} {\ rm {{d} \ alpha}}} f _ {\ varphi} (q ^ {1} + \ alpha ~ c ^ {1}, q ^ {2} + \ alpha ~ c ^ {2}, q ^ {3} + \ alpha ~ c ^ {3}) {\ biggr |} _ {\ alpha = 0} = {\ cfrac {\ partial f _ {\ varphi}} {\ partial q ^ {i}} } ~ c ^ {i} = {\ cfrac {\ partial f} {\ partial q ^ {i}}} ~ c ^ {i}}Если мы положим f (x) = ψ i ( x) {\ displaystyle f (\ mathbf {x}) = \ psi ^ {i} (\ mathbf {x})}, то, поскольку qi = ψ i (x) {\ displaystyle q ^ {i} = \ psi ^ {i} (\ mathbf {x})}, имеем[∇ ψ i (x)] ⋅ c = ∂ ψ i ∂ qjcj = ci {\ displaystyle [{\ boldsymbol {\ nabla}} \ psi ^ {i} (\ mathbf {x})] \ cdot \ mathbf {c} = {\ cfrac {\ partial \ psi ^ {i}} {\ partial q ^ {j}}} ~ c ^ {j} = c ^ {i}}, который обеспечивает средства извлечения контравариантной компоненты вектора c. Если biявляется ковариантным (или естественным) базисом в точке, и если b является контравариантным (или обратным) базисом в этой точке, то[∇ f (x)] ⋅ c = ∂ f ∂ qi ci знак равно (∂ е ∂ qibi) (cibi) ⇒ ∇ е (x) = ∂ f ∂ qibi {\ displaystyle [{\ boldsymbol {\ nabla}} f (\ mathbf {x})] \ cdot \ mathbf {c} = {\ cfrac {\ partial f} {\ partial q ^ {i}}} ~ c ^ {i} = \ left ({\ cfrac {\ partial f} {\ partial q ^ {i}}} ~ \ mathbf {b} ^ {i} \ right) \ left (c ^ {i} ~ \ mathbf {b} _ {i} \ right) \ quad \ Rightarrow \ quad {\ boldsymbol {\ nabla}} f (\ mathbf { x}) = {\ cfrac {\ partial f} {\ partial q ^ {i}}} ~ \ mathbf {b} ^ {i}}Краткое обоснование этого выбора базиса дается в следующем разделе..Векторное полеАналогичный процесс можно использовать для получения градиента a vector field f(x). The gradient is given by[ ∇ f ( x) ] ⋅ c = ∂ f ∂ q i c i {\displaystyle [{\boldsymbol {\nabla }}\mathbf {f} (\mathbf {x})]\cdot \mathbf {c} ={\cfrac {\partial \mathbf {f} }{\partial q^{i}}}~c^{i}}If we consider the gradient of the position vector field r(x) = x, then we can show thatc = ∂ x ∂ q i c i = b i ( x) c i ; b i ( x) := ∂ x ∂ q i{\displaystyle \mathbf {c} ={\cfrac {\partial \mathbf {x} }{\partial q^{i}}}~c^{i}=\mathbf {b} _{i}(\mathbf {x})~c^{i}~;~~\mathbf {b} _{i}(\mathbf {x}):={\cfrac {\partial \mathbf {x} }{\partial q^{i}}}}The vector field biis tangent to the q coordinate curve and forms a natural basisat each point on the curve. This basis, as discussed at the beginning of this article, is also called the covariantcurvilinear basis. We can also define a reciprocal basis, or contravariantcurvilinear basis, b. 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. Since cis arbitrary, we can write∇ f ( x) = ∂ f ∂ q i ⊗ b i {\displaystyle {\boldsymbol {\nabla }}\mathbf {f} (\mathbf {x})={\cfrac {\partial \mathbf {f} }{\partial q^{i}}}\otimes \mathbf {b} ^{i}}Note that the contravariant basis vector bis perpendicular to the surface of constant ψ and is given byb i = ∇ ψ i {\displaystyle \mathbf {b} ^{i}={\boldsymbol {\nabla }}\psi ^{i}}Christoffel symbols of the first kindThe Christoffel symbols of the first kind are defined asb i, j = ∂ b i ∂ q j := Γ i j k b k ⇒ b i, j ⋅ b l = Γ i j l {\displaystyle \mathbf {b} _{i,j}={\frac {\partial \mathbf {b} _{i}}{\partial q^{j}}}:=\Gamma _{ijk}~\mathbf {b} ^{k}\quad \Rightarrow \quad \mathbf {b} _{i,j}\cdot \mathbf {b} _{l}=\Gamma _{ijl}}To express Γijkin terms of gijwe note thatg ij, k = ( b i ⋅ b j), k = b i, k ⋅ b j + b i ⋅ b j, k = Γ i k j + Γ j k i g i k, j = ( b i ⋅ b k), j = b i, j ⋅ b k + b i ⋅ b k, j = Γ i j k + Γ k j i g j k, i = ( b j ⋅ b k), i = b j, i ⋅ b k + b j ⋅ b k, i = Γ j i k + Γ k i j {\displaystyle {\begin{aligned}g_{ij,k}=(\mathbf {b} _{i}\cdot \mathbf {b} _{j})_{,k}=\mathbf {b} _{i,k}\cdot \mathbf {b} _{j}+\mathbf {b} _{i}\cdot \math bf {b} _{j,k}=\Gamma _{ikj}+\Gamma _{jki}\\g_{ik,j}=(\mathbf {b} _{i}\cdot \mathbf {b} _{k})_{,j}=\mathbf {b} _{i,j}\cdot \mathbf {b} _{k}+\mathbf {b} _{i}\cdot \mathbf {b} _{k,j}=\Gamma _{ijk}+\Gamma _{kji}\\g_{jk,i}=(\mathbf {b} _{j}\cdot \mathbf {b} _{k})_{,i}=\mathbf {b} _{j,i}\cdot \mathbf {b} _{k}+\mathbf {b} _{j}\cdot \mathbf {b} _{k,i}=\Gamma _{jik}+\Gamma _{kij}\end{aligned}}}Since b i,j= bj,iwe have Γijk= Γjik. Using these to rearrange the above relations givesΓ i j k = 1 2 ( g i k, j + g j k, i − g i j, k) = 1 2 [ ( b i ⋅ b k), j + ( b j ⋅ b k), i − ( b i ⋅ b j), k ] {\displaystyle \Gamma _{ijk}={\frac {1}{2}}(g_{ik,j}+g_{jk,i}-g_{ij,k})={\frac {1}{2}}[(\mathbf {b} _{i}\cdot \mathbf {b} _{k})_{,j}+(\mathbf {b} _{j}\cdot \mathbf {b} _{k})_{,i}-(\mathbf {b} _{i}\cdot \mathbf {b} _{j})_{,k}]}Christoffel symbols of the second kindThe Christoffel symbols of the second kind are defined asΓ i j k = Γ j i k {\displaystyle \Gamma _{ij}^{k}=\Gamma _{ji}^{k}}in which∂ b i ∂ q j = Γ i j k b k {\displaystyle {\cfrac {\partial \mathbf {b} _{i}}{\partial q^{j}}}=\Gamma _{ij}^{k}~\mathbf {b} _{k}}This implies thatΓ i j k = ∂ b i ∂ q j ⋅ b k = − b i ⋅ ∂ b k ∂ q j {\displaystyle \Gamma _{ij}^{k}={\cfrac {\partial \mathbf {b} _{i}}{\partial q^{j}}}\cdot \mathbf {b} ^{k}=-\mathbf {b} _{i}\cdot {\cfrac {\partial \mathbf {b} ^{k}}{\partial q^{j}}}}Other relations that follow are∂ b i ∂ q j = − Γ j k i b k ; ∇ b i = Γ i j k b k ⊗ b j ; ∇ b i = − Γ j k i b k ⊗ b j{\ displaystyle {\ cfrac {\ partial \ mathbf {b} ^ {i}} {\ partial q ^ {j}}} = - \ Gamma _ {jk} ^ {i} ~ \ mathbf {b} ^ {k } ~; ~~ {\ boldsymbol {\ nabla}} \ mathbf {b} _ {i} = \ Gamma _ {ij} ^ {k} ~ \ mathbf {b} _ {k} \ otimes \ mathbf {b} ^ {j} ~; ~~ {\ boldsymbol {\ nabla}} \ mathbf {b} ^ {i} = - \ Gamma _ {jk} ^ {i} ~ \ mathbf {b} ^ {k} \ otimes \ mathbf {b} ^ {j}}Еще одно особенно полезное соотношение, которое показывает, что символ Кристоффеля зависит от l y от метрического тензора и его производных, этоΓ ijk = gkm 2 (∂ gmi ∂ qj + ∂ gmj ∂ qi - ∂ gij ∂ qm) {\ displaystyle \ Gamma _ {ij} ^ {k} = {\ frac {g ^ {km}} {2}} \ left ({\ frac {\ partial g_ {mi} } {\ partial q ^ {j}}} + {\ frac {\ partial g_ {mj}} {\ partial q ^ {i}}} - {\ frac {\ partial g_ {ij}} {\ partial q ^ {m}}} \ right)}Явное выражение для градиента векторного поляСледующие выражения для градиента векторного поля в криволинейных координатах весьма полезны.∇ v = [∂ vi ∂ qk + Γ lkivl] bi ⊗ bk = [∂ vi ∂ qk - Γ kilvl] bi ⊗ bk {\ displaystyle {\ begin {align} {\ boldsymbol {\ nabla}} \ mathbf {v} = \ left [{\ cfrac {\ partial v ^ {i}} {\ partial q ^ {k}}} + \ Gamma _ {lk} ^ {i} ~ v ^ {l} \ right] ~ \ mathbf {b} _ {i} \ otimes \ mathbf {b} ^ {k} \\ [8pt] = \ left [{\ cfrac {\ partial v_ {i}} {\ partial q ^ {k} }} - \ Gamma _ {ki} ^ {l} ~ v_ {l} \ right] ~ \ mathbf {b} ^ {i} \ otimes \ mathbf {b} ^ {k} \ end {align}}}Представление физического векторного поляВекторное поле v может быть представлено какv = vibi = v ^ ib ^ i {\ displaystyle \ mathbf {v} = v_ { i} ~ \ mathbf {b} ^ {i} = {\ hat {v}} _ {i} ~ {\ hat {\ mathbf {b}}} ^ {i}}где vi { \ displaystyle v_ {i}}- ковариантные компоненты поля, v ^ i {\ displaystyle {\ hat {v}} _ {i}}- физические компоненты и (без суммирования )b ^ i = bigii {\ displaystyle {\ hat {\ mathbf {b}}} ^ {i} = {\ cfrac {\ m athbf {b} ^ {i }} {\ sqrt {g ^ {ii}}}}}- нормализованный контравариантный базис v эктор.Тензорное поле второго порядкаГрадиент тензорного поля второго порядка может аналогичным образом быть выражено как∇ S = ∂ S ∂ qi ⊗ bi {\ displaystyle {\ boldsymbol {\ nabla}} {\ boldsymbol {S}} = {\ cfrac {\ partial {\ boldsymbol {S}}} {\ частичное q ^ {i}}} \ otimes \ mathbf {b} ^ {i}}Явные выражения для градиентаЕсли мы рассмотрим выражение для тензора в терминах контравариантного базиса, то∇ S = ∂ ∂ qk [S ijbi ⊗ bj] ⊗ bk = [∂ S ij ∂ qk - Γ kil S lj - Γ kjl S il] bi ⊗ bj ⊗ bk {\ displaystyle {\ boldsymbol {\ nabla} } {\ boldsymbol {S}} = {\ cfrac {\ partial} {\ partial q ^ {k}}} [S_ {ij} ~ \ mathbf {b} ^ {i} \ otimes \ mathbf {b} ^ { j}] \ otimes \ mathbf {b} ^ {k} = \ left [{\ cfrac {\ partial S_ {ij}} {\ partial q ^ {k}}} - \ Gamma _ {ki} ^ {l} ~ S_ {lj} - \ Gamma _ {kj} ^ {l} ~ S_ {il} \ right] ~ \ mathbf {b} ^ {i} \ otimes \ mathbf {b} ^ {j} \ otimes \ mathbf { b} ^ {k}}Мы также можем написать∇ S = [∂ S ij ∂ qk + Γ kli S lj + Γ klj S il] bi ⊗ bj ⊗ bk = [∂ S ji ∂ qk + Γ kli S jl - Γ kjl S li] bi ⊗ bj ⊗ bk = [∂ S ij ∂ qk - Γ ikl S lj + Γ klj S il] bi ⊗ bj ⊗ b к {\ displaystyle {\ begin {align} {\ boldsymbol {\ nabla}} {\ boldsymbol {S}} = \ left [{\ cfrac {\ partial S ^ {ij}} {\ partial q ^ {k} }} + \ Gamma _ {kl} ^ {i} ~ S ^ {lj} + \ Gamma _ {kl} ^ {j} ~ S ^ {il} \ right] ~ \ mathbf {b} _ {i} \ otimes \ mathbf {b} _ {j} \ otimes \ mathbf {b} ^ {k} \\ [8pt] = \ left [{\ cfrac {\ partial S_ {~ j} ^ {i}} {\ partial q ^ {k}}} + \ Gamma _ {kl} ^ {i} ~ S_ {~ j} ^ {l} - \ Gamma _ {kj} ^ {l} ~ S_ {~ l} ^ {i} \ right] ~ \ mathbf {b} _ {i} \ otimes \ mathbf {b} ^ {j} \ otimes \ mathbf {b} ^ {k} \\ [8pt] = \ left [{\ cfrac {\ partial S_ {i} ^ {~ j}} {\ partial q ^ {k}}} - \ Gamma _ {ik} ^ {l} ~ S_ {l} ^ {~ j} + \ Gamma _ {kl} ^ { j} ~ S_ {i} ^ {~ l} \ right] ~ \ mathbf {b} ^ {i} \ otimes \ mathbf {b} _ {j} \ otimes \ mathbf {b} ^ {k} \ end { выровнено}}}Представление физического тензорного поля второго порядкаФизические компоненты тензорного поля второго порядка могут быть получены с использованием нормированного контравариантного базиса, т. е.S = S ijbi ⊗ bj = S ^ ijb ^ i ⊗ b ^ j {\ displaystyle {\ boldsymbol {S}} = S_ {ij} ~ \ mathbf {b} ^ {i} \ otimes \ mathbf {b} ^ {j} = { \ hat {S}} _ {ij} ~ {\ hat {\ mathbf {b}}} ^ {i} \ otimes {\ hat {\ mathbf {b}}} ^ {j}}, где штрихованные базисные векторы нормализованы. Это означает, что (опять же без суммирования)S ^ ij = S ijgiigjj {\ displaystyle {\ hat {S}} _ {ij} = S_ {ij} ~ {\ sqrt {g ^ {ii} ~ g ^ { jj}}}}ДивергенцияВекторное полеДивергенция векторного поля (v {\ displaystyle \ mathbf {v}}) определяется какdiv v = ∇ ⋅ v = tr (∇ v) {\ displaystyle \ operatorname {div} ~ \ mathbf {v} = {\ boldsymbol {\ nabla}} \ cdot \ mathbf {v} = {\ text {tr}} ({\ boldsymbol {\ nabla}} \ mathbf {v})}В терминах компонентов относительно криволинейного базиса∇ ⋅ v = ∂ vi ∂ qi + Γ ℓ iiv ℓ = [∂ vi ∂ qj - Γ ji ℓ v ℓ] gij {\ displaystyle {\ boldsymbol {\ nabla}} \ cdot \ mathbf {v} = {\ cfrac {\ partial v ^ { i}} {\ partial q ^ {i}}} + \ Gamma _ {\ ell i} ^ {i} ~ v ^ {\ ell} = \ left [{\ cfrac {\ partial v_ {i}} {\ частичное q ^ {j}}} - \ Gamma _ {ji} ^ {\ ell} ~ v _ {\ ell} \ right] ~ g ^ {ij}}Альтернативным уравнением для расходимости векторного поля является часто используемый. Чтобы вывести это соотношение, напомним, что∇ ⋅ v = ∂ vi ∂ qi + Γ ℓ iiv ℓ {\ displaystyle {\ boldsymbol {\ nabla}} \ cdot \ mathbf {v} = {\ frac {\ partial v ^ { i}} {\ partial q ^ {i}}} + \ Gamma _ {\ ell i} ^ {i} ~ v ^ {\ ell}}ТеперьΓ ℓ ii = Γ i ℓ i знак равно gmi 2 [∂ gim ∂ q ℓ + ∂ g ℓ m ∂ qi - ∂ gil ∂ qm] {\ displaystyle \ Gamma _ {\ ell i} ^ {i} = \ Gamma _ {i \ ell} ^ {i} = {\ cfrac {g ^ {mi}} {2}} \ left [{\ frac {\ partial g_ {im}} {\ partial q ^ {\ ell}}} + {\ frac {\ partial g _ {\ ell m}} {\ partial q ^ {i}}} - {\ frac {\ partial g_ {il}} {\ partial q ^ {m}}} \ right]}Отмечая, что из-за симметрии из g {\ displaystyle {\ boldsymbol {g}}},gmi ∂ g ℓ m ∂ qi = gmi ∂ gi ℓ ∂ qm {\ displaystyle g ^ {mi} ~ {\ frac {\ partial g_ { \ ell m}} {\ partial q ^ {i}}} = g ^ {mi} ~ {\ frac {\ partial g_ {i \ ell}} {\ partial q ^ {m}}}}мы иметь∇ ⋅ v = ∂ vi ∂ qi + gmi 2 ∂ gim ∂ q ℓ v ℓ {\ displaystyle {\ boldsymbol {\ nabla}} \ cdot \ mathbf {v} = {\ frac {\ partial v ^ { i}} {\ partial q ^ {i}}} + {\ cfrac {g ^ {mi}} {2}} ~ {\ frac {\ partial g_ {im}} {\ partial q ^ {\ ell}}}~v^{\ell }}Recall that if [gij] is the matrix whose components are gij, then the inverse of the matrix is [ g i j ] − 1 = [ g i j ] {\displaystyle [g_{ij}]^{-1}=[g^{ij}]}. The inverse of the matrix is given by[ g i j ] = [ g i j ] − 1 = A i j g ; g := det ( [ g i j ]) = det g {\displaystyle [g^{ij}]=[g_{ij}]^{-1}={\cfrac {A^{ij}}{g}}~;~~g:=\det( [g_{ij}])=\det {\boldsymbol {g}}}where A are the Cofactor matrix of the components gij. From matrix algebra we haveg = det ( [ g i j ]) = ∑ i g i j A i j ⇒ ∂ g ∂ g i j = A i j {\displaystyle g=\det([g_{ij}])=\sum _{i}g_{ij}~A^{ij}\quad \Rightarrow \quad {\frac {\partial g}{\partial g_{ij}}}=A^{ij}}Hence,[ g i j ] = 1 g ∂ g ∂ g i j {\displaystyle [g^{ij}]={\cfrac {1}{g}}~{\frac {\partial g}{\partial g_{ij}}}}Plugging this relation into the expression for the divergence gives∇ ⋅ v = ∂ v i ∂ q i + 1 2 g ∂ g ∂ g m i ∂ g i m ∂ q ℓ v ℓ = ∂ v i ∂ q i + 1 2 g ∂ g ∂ q ℓ v ℓ {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\frac {\partial v^{i}}{\partial q^{i}}}+{\cfrac {1}{2g}}~{\frac {\partial g}{\partial g_{mi}}}~{\frac {\partial g_{im}}{\partial q^{\ell }}}~v^{\ell }={\frac {\partial v^{i}}{\partial q^{i}}}+{\cfrac {1}{2g}}~{\frac {\partial g}{\partial q^{\ell }}}~v^{\ell }}A little manipulation leads to the more compact form∇ ⋅ v = 1 g ∂ ∂ q i ( v i g) {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\cfrac {1}{\sqrt {g}}}~{\frac {\partial }{\partial q^{i}}}(v^{i}~{\sqrt {g}})}Second-order tensor fieldThe divergence of a second-order tensor field is defined using( ∇ ⋅ S) ⋅ a = ∇ ⋅ ( S ⋅ a) {\displaystyle ({\boldsymbol {\nabla }}\cdot {\boldsymbol {S}})\cdot \mathbf{a} ={\boldsymbol {\nabla }}\cdot ({\boldsymbol {S}}\cdot \mathbf {a})}where ais an arbitrary constant vector. In curvilinear coordinates,∇ ⋅ S = [ ∂ S i j ∂ q k − Γ k i l Slj - Γ kjl S il] gikbj = [∂ S ij ∂ qi + Γ ili S lj + Γ ilj S il] bj = [∂ S ji ∂ qi + Γ ili S jl - Γ ijl S li] bj = [∂ S ij ∂ qk - Γ ikl S lj + Γ klj S il] gikbj {\ displaystyle {\ begin {align} {\ boldsymbol {\ nabla}} \ cdot {\ boldsymbol {S}} = \ left [{\ cfrac { \ partial S_ {ij}} {\ partial q ^ {k}}} - \ Gamma _ {ki} ^ {l} ~ S_ {lj} - \ Gamma _ {kj} ^ {l} ~ S_ {il} \ right] ~ g ^ {ik} ~ \ math bf {b} ^ {j} \\ [8pt] = \ left [{\ cfrac {\ partial S ^ {ij}} {\ partial q ^ {i}} } + \ Gamma _ {il} ^ {i} ~ S ^ {lj} + \ Gamma _ {il} ^ {j} ~ S ^ {il} \ right] ~ \ mathbf {b} _ {j} \\ [8pt] = \ left [{\ cfrac {\ partial S_ {~ j} ^ {i}} {\ partial q ^ {i}}} + \ Gamma _ {il} ^ {i} ~ S_ {~ j } ^ {l} - \ Gamma _ {ij} ^ {l} ~ S_ {~ l} ^ {i} \ right] ~ \ mathbf {b} ^ {j} \\ [8pt] = \ left [{ \ cfrac {\ partial S_ {i} ^ {~ j}} {\ partial q ^ {k}}} - \ Gamma _ {ik} ^ {l} ~ S_ {l} ^ {~ j} + \ Gamma _ {kl} ^ {j} ~ S_ {i} ^ {~ l} \ right] ~ g ^ {ik} ~ \ mathbf {b} _ {j} \ end {align}}}лапласианСкалярное полеЛапласиан скалярного поля φ (x ) определяется как∇ 2 φ: = ∇ ⋅ (∇ φ) {\ displaystyle \ na bla ^ {2} \ varphi: = {\ boldsymbol {\ nabla}} \ cdot ({\ boldsymbol {\ nabla}} \ varphi)}Использование альтернативного выражения для дивергенции векторного поля дает нам∇ 2 φ знак равно 1 г ∂ ∂ qi ([∇ φ] ig) {\ displaystyle \ nabla ^ {2} \ varphi = {\ cfrac {1} {\ sqrt {g}}} ~ {\ frac {\ partial } {\ partial q ^ {i}}} ([{\ boldsymbol {\ nabla}} \ varphi] ^ {i} ~ {\ sqrt {g}})}Теперь∇ φ = ∂ φ ∂ qlbl знак равно gli ∂ φ ∂ qlbi ⇒ [∇ φ] я = gli ∂ φ ∂ ql {\ displaystyle {\ boldsymbol {\ nabla}} \ varphi = {\ frac {\ partial \ varphi} {\ partial q ^ {l }}} ~ \ mathbf {b} ^ {l} = g ^ {li} ~ {\ frac {\ partial \ varphi} {\ partial q ^ {l}}} ~ \ mathbf {b} _ {i} \ quad \ Rightarrow \ quad [{\ boldsymbol {\ nabla}} \ varphi] ^ {i} = g ^ {li} ~ {\ frac {\ partial \ varphi} {\ partial q ^ {l}}}}Следовательно,∇ 2 φ = 1 г ∂ ∂ qi (gli ∂ φ ∂ qlg) {\ displaystyle \ nabla ^ {2} \ varphi = {\ cfrac {1} {\ sqrt {g}}} ~ { \ frac {\ partial} {\ partial q ^ {i}}} \ left (g ^ {li} ~ {\ frac {\ partial \ varphi} {\ partial q ^ {l}}} ~ {\ sqrt {g }} \ right)}Ротор векторного поля The curl of a vector field vin covariant curvilinear coordinates can be written as∇ × v = E r s t v s | r b t {\displaystyle {\boldsymbol {\nabla }}\times \mathbf {v} ={\mathcal {E}}^{rst}v_{s|r}~\mathbf {b} _{t}}wherev s | r = v s, r − Γ s r i v i {\displaystyle v_{s|r}=v_{s,r}-\Gamma _{sr}^{i}~v_{i}}Orthogonal curvilinear coordinatesAssume, for thepurposes of this section, that the curvilinear coordinate system is orthogonal, i.e.,b i ⋅ b j = { g i i if i = j 0 if i ≠ j, {\displaystyle \mathbf {b} _{i}\cdot \mathbf {b} _{j}={\begin{cases}g_{ii}{\text{if }}i=j\\0{\text{if }}i\neq j,\end{cases}}}or equivalently,b i ⋅ b j = { g i i if i = j 0 if i ≠ j, {\displaystyle \mathbf {b} ^{i}\cdot \mathbf {b}^{j}={\begin{cases}g^{ii}{\text{if }}i=j\\0{\text{if }}i\neq j,\end{cases}}}where g i i = g i i − 1 {\displaystyle g^{ii}=g_{ii}^{-1}}. As before, b i, b j {\displaystyle \mathbf {b} _{i},\mathbf {b} _{j}}are covariant basis vectors and b, bare contravariant basis vectors. Also, let (e, e, e) be a background, fixed, Cartesian basis. A list of orthogonal curvilinear coordinates is given below.Metric tensor in orthogonal curvilinear coordinatesLet r(x) be the position vector of the point xwith respect to the origin of the coordinate system. The notation can be simplified by noting that x= r(x). At each point we can construct a small line element dx. The square of the length of the line element is the scalar product dx• dxand is called the metric of the space. Recall that the space of interest is assumed to be Euclidean when we talk of curvilinear coordinates. Let us express the position vector in terms of the background, fixed, Cartesian basis, i.e.,x = ∑ i = 1 3 x i e i {\displaystyle \mathbf {x} =\sum _{i=1}^{3}x_{i}~\mathbf {e} _{i}} Using the chain rule, we can then express dxin terms of three-dimensional orthogonal curvilinear coordinates (q, q, q) asd x = ∑ i = 1 3 ∑ j = 1 3 ( ∂ x i ∂ q j e i) d q j {\displaystyle \mathrm {d} \mathbf {x} =\sum _{i=1}^{3}\sum _{j=1}^{3}\left({\cfrac {\partial x_{i}}{\partial q^{j}}}~\mathbf {e} _{i}\right)\mathrm {d} q^{j}}Therefore, the metric is given byd x ⋅ d x = ∑ i = 1 3 ∑ j = 1 3 ∑ k = 1 3 ∂ x i ∂ q j ∂ x i ∂ q k d q j d q k {\displaystyle \mathrm {d} \mathbf {x} \cdot \mathrm {d} \mathbf {x} =\sum _{i=1}^{3}\sum _{j=1}^{3}\sum _{k=1}^{3}{\cfrac {\partial x_{i}}{\partial q^{j}}}~{\cfrac {\partial x_{i}}{\partial q^{k}}}~\mathrm {d} q^{j}~\mathrm {d} q^{k}}The symmetric quantityg i j ( q i, q j) = ∑ k = 1 3 ∂ x k ∂ q i ∂ x k ∂ q j = b i ⋅ b j {\displaystyle g_{ij}(q^{i},q^{j})=\sum _{k=1}^{3}{\cfrac {\partial x_{k}}{\partial q^{i}}}~{\cfrac {\partial x_{k}}{\partial q^{j}}}=\mathbf {b} _{i}\cdot \mathbf {b} _{j}}is called the fundamental (or metric) tensor of the Euclidean space in curvilinear coordinates. Note also thatg i j = ∂ x ∂ q i ⋅ ∂ x ∂ q j = ( ∑ k h k i e k) ⋅ ( ∑ m h m j e m) = ∑ k h k i h k j {\displaystyle g_{ij}={\cfrac {\partial \mathbf {x} }{\partial q^{i}}}\cdot {\cfrac {\partial \mathbf {x} }{\partial q^{j}}}=\left(\sum _{k}h_{ki}~\mathbf {e} _{k}\right)\cdot \left(\sum _{m}h_{mj}~\mathbf {e} _{m}\right)=\sum _{k}h_{ki}~h_{kj}}where hijare the Lamé coefficients. If we define the scale factors, hi, usingb i ⋅ b i = g i i = ∑ k h k i 2 =: h i 2 ⇒ | ∂ x ∂ q i | = | b i | = g i i = h i {\displaystyle \mathbf {b} _{i}\cdot \mathbf {b} _{i}=g_{ii}=\sum _{k}h_{ki}^{2}=:h_{i}^{2}\quad \Rightarrow \quad\left|{\cfrac {\partial \mathbf {x} }{\partial q^{i}}}\right|=\left|\mathbf {b} _{i}\right|={\sqrt {g_{ii}}}=h_{i}}we get a relation between the fundamental tensor and the Lamé coefficients. Пример: полярные координатыЕсли мы рассмотрим полярные координаты для R, обратите внимание, что(x, y) = (r cos θ, r sin θ) {\ displaystyle ( x, y) = (r \ cos \ theta, r \ sin \ theta)}(r, θ) - криволинейные координаты, а определитель Якоби преобразования (r, θ) → (r cos θ, r sin θ) равно r. Ортогональные базисные векторы равны br= (cos θ, sin θ), bθ= (−r sin θ, r cos θ). Нормализованные базисные векторы: er= (cos θ, sin θ), eθ= (−sin θ, cos θ), а масштабные коэффициенты - h r = 1 и h θ. = р. Фундаментальный тензор равен g 11 = 1, g 22 = r, g 12 = g 21 = 0.Линейные и поверхностные интегралыЕсли мы хотим использовать криволинейные координаты для вычислений векторного исчисления, необходимо внести корректировки в вычисление линейных, поверхностных и объемных интегралов. Для простоты мы снова ограничим обсуждение тремя измерениями и ортогональными криволинейными координатами. Однако те же аргументы применимы к n {\ displaystyle n}-мерным задачам, хотя в выражениях есть некоторые дополнительные члены, когда система координат не ортогональна.Линейные интегралыОбычно при вычислении линейных интегралов нас интересует вычисление∫ C f d s = ∫ a b f (x (t)) | ∂ x ∂ t | dt {\ displaystyle \ int _ {C} f \, ds = \ int _ {a} ^ {b} f (\ mathbf {x} (t)) \ left | {\ partial \ mathbf {x} \ over \ partial t} \ right | \; dt}где x (t) параметризует C в декартовых координатах. В криволинейных координатах член| ∂ x ∂ t | = | ∑ i = 1 3 ∂ x ∂ q i ∂ q i ∂ t | {\ displaystyle \ left | {\ partial \ mathbf {x} \ over \ partial t} \ right | = \ left | \ sum _ {i = 1} ^ {3} {\ partial \ mathbf {x} \ over \ partial q ^ {i}} {\ partial q ^ {i} \ over \ partial t} \ right |}по правилу цепочки . И из определения коэффициентов Ламе,∂ x ∂ q i = ∑ k h k i e k {\ displaystyle {\ partial \ mathbf {x} \ over \ partial q ^ {i}} = \ sum _ {k} h _{ki}~\mathbf {e} _{k}}and thus| ∂ x ∂ t | = | ∑ k ( ∑ i h k i ∂ q i ∂ t) e k | = ∑ i ∑ j ∑ k h k i h k j ∂ q i ∂ t ∂ q j ∂ t = ∑ i ∑ j g i j ∂ q i ∂ t ∂ q j ∂ t {\displaystyle {\begin{aligned}\left|{\partial \mathbf {x} \over \partial t}\right|=\left|\sum _{k}\left(\sum _{i}h_{ki}~{\cfrac {\partial q^{i}}{\partial t}}\right)\mathbf {e} _{k}\right|\\[8pt]={\sqrt {\sum _{i}\sum _{j}\sum_{k}h_{ki}~h_{kj}{\cfrac {\partial q^{i}}{\partial t}}{\cfrac {\partial q^{j}}{\partial t}}}}={\sqrt {\sum _{i}\sum _{j}g_{ij}~{\cfrac {\partial q^{i}}{\partial t}}{\cfrac {\partial q^{j}}{\partial t}}}}\end{aligned}}}Now, since g i j = 0 {\displaystyle g_{ij}=0}when i ≠ j {\displaystyle i\neq j}, we have| ∂ x ∂ t | = ∑ i g i i ( ∂ q i ∂ t) 2 = ∑ ih i 2 ( ∂ q i ∂ t) 2 {\displaystyle \left|{\partial \mathbf {x} \over \partial t}\right|={\sqrt {\sum _{i}g_{ii}~\left({\cfrac {\partial q^{i}}{\partial t}}\right)^{2}}}={\sqrt {\sum _{i}h_{i}^{2}~\left({\cfrac {\partial q^{i}}{\partial t}}\right)^{2}}}}and we can proceed normally.Surface integralsLikewise, if we are interested in a surface integral, the relevant calculation, with the parameterization of the surface in Cartesian coordinates is:∫ S f d S = ∬ T f ( x ( s, t)) | ∂ x ∂ s × ∂ x ∂ t | d s d t {\displaystyle \int _{S}f\,dS=\iint _{T}f(\mathbf {x} (s,t))\left|{\partial \mathbf {x} \over \partial s}\times {\partial \mathbf {x} \over \partial t}\right|\,ds\,dt}Again, in curvilinear coordinates, we have| ∂ x ∂ s × ∂ x ∂ t | = | ( ∑ i ∂ x ∂ q i ∂ q i ∂ s) × ( ∑ j ∂ x ∂ q j ∂ q j ∂ t) | {\displaystyle \left|{\partial \mathbf {x} \over \partial s}\times {\partial \mathbf {x} \over \partial t}\right|=\left|\left(\sum _{i}{\partial \mathbf {x} \over \partial q^{i}}{\partial q^{i} \over \partial s}\right)\times \left(\sum _{j}{\partial \mathbf {x} \over \partial q^{j}}{\partial q^{j} \over \partial t}\right)\right|}и снова используем определение криволинейных координат, чтобы получить∂ x ∂ q i ∂ q i ∂ s = ∑ k (∑ i = 1 3 h k i ∂ q i ∂ s) e k; ∂ x ∂ qj ∂ qj ∂ T знак равно ∑ м (∑ j = 1 3 hmj ∂ qj ∂ t) em {\ displaystyle {\ partial \ mathbf {x} \ over \ partial q ^ {i}} {\ partial q ^ {i} \ over \ partial s} = \ sum _ {k} \ left (\ sum _ {i = 1} ^ {3} h_ {ki} ~ {\ partial q ^ {i} \ over \ partial s} \ right) \ mathbf {e} _ {k} ~; ~~ {\ partial \ mathbf {x} \ over \ partial q ^ {j}} {\ partial q ^ {j} \ over \ partial t} = \ сумма _ {m} \ left (\ sum _ {j = 1} ^ {3} h_ {mj} ~ {\ partial q ^ {j} \ over \ partial t} \ right) \ mathbf {e} _ {m }}Следовательно,| ∂ x ∂ s × ∂ x ∂ t | = | ∑ k ∑ m (∑ i = 1 3 h k i ∂ q i ∂ s) (∑ j = 1 3 h m j ∂ q j ∂ t) e k × e m | = | ∑ p ∑ k ∑ m E k m p (∑ i = 1 3 h k i ∂ q i ∂ s) (∑ j = 1 3 h m j ∂ q j ∂ t) e p | {\ displaystyle {\ begin {align} \ left | {\ partial \ mathbf {x} \ over \ partial s} \ times {\ partial \ mathbf {x} \ over \ partial t} \ right | = \ left | \ sum _ {k} \ sum _ {m} \ left (\ sum _ {i = 1} ^ {3} h_ {ki} ~ {\ partial q ^ {i} \ over \ partial s} \ right) \ left (\ sum _ {j = 1} ^ {3} h_ {mj} ~ {\ partial q ^ {j} \ over \ partial t} \ right) \ mathbf {e} _ {k} \ times \ mathbf { e} _ {m} \ right | \\ [8pt] = \ left | \ sum _ {p} \ sum _ {k} \ sum _ {m} {\ mathcal {E}} _ {kmp} \ left (\ sum _ {i = 1} ^ {3} h_ {ki} ~ {\ partial q ^ {i} \ over \ partial s} \ right) \ left (\ sum _ {j = 1} ^ {3} h_ {mj} ~ {\ partial q ^ {j} \ over \ partial t} \ right) \ mathbf {e} _ {p} \ right | \ end {align}}}где E { \ displaystyle {\ mathcal {E}}}- это символ перестановки . В детерминантной форме перекрестное произведение с точки зрения криволинейных координат будет:| e 1 e 2 e 3 ∑ ih 1 i ∂ qi ∂ s ∑ ih 2 i ∂ qi ∂ s ∑ ih 3 i ∂ qi ∂ s ∑ jh 1 j ∂ qj ∂ t ∑ jh 2 j ∂ qj ∂ t ∑ jh 3 j ∂ qj ∂ t | {\ displaystyle {\ begin {vmatrix} \ mathbf {e} _ {1} \ mathbf {e} _ {2} \ mathbf {e} _ {3} \\ \\\ сумма _ {i} h_ {1i} {\ partial q ^ {i} \ over \ partial s} \ sum _ {i} h_ {2i} {\ partial q ^ {i} \ over \ partial s} \ sum _ {i} h_ {3i} {\ partial q ^ {i} \ over \ partial s} \\ \\\ сумма _ {j} h_ {1j} {\ partial q ^ {j} \ over \ partial t} \ sum _ {j} h_ {2j} {\ partial q ^ {j} \over \partial t}\sum _{j}h_{3j}{\partial q^{j} \over \partial t}\end{vmatrix}}}Grad, curl, div, LaplacianIn orthogonal curvilinear coordinates of 3 dimensions, whereb i = ∑ k g i k b k ; g i i = 1 g i i = 1 h i 2 {\displaystyle \mathbf {b} ^{i}=\sum _{k}g^{ik}~\mathbf {b} _{k}~;~~g^{ii}={\cfrac {1}{g_{ii}}}={\cfrac {1}{h_{i}^{2}}}}one can express the gradient of a scalar or vector field as∇ φ = ∑ i ∂ φ ∂ q i b i = ∑ i ∑ j ∂ φ ∂ q i g i j b j = ∑ i 1 h i 2 ∂ f ∂ q i b i ; ∇ v = ∑ i 1 h i 2 ∂ v ∂ q i ⊗ b i {\displaystyle \nabla \varphi =\sum _{i}{\partial \varphi \over \partial q^{i}}~\mathbf {b} ^{i}=\sum _{i}\sum _{j}{\partial \varphi \over \partial q^{i}}~g^{ij}~\mathbf {b} _{j}=\sum _{i}{\cfrac {1}{h_{i}^{2}}}~{\partial f \over \partial q^{i}}~\mathbf {b} _{i}~;~~\nabla \mathbf {v} =\sum _{i}{\cfrac {1}{h_{i}^{2}}}~{\partial \mathbf {v} \over \partial q^{i}}\otimes \mathbf {b} _{i}}For an orthogonal basisg = g 11 g 22 g 33 = h 1 2 h 2 2 h 3 2 ⇒ g = h 1 h 2 h 3 {\displaystyle g=g_{11}~g_{22}~g_{33}=h_{1}^{2}~h_{2}^{2}~h_{3}^{2}\quad \Rightarrow \quad {\sqrt {g}}=h_{1}h_{2}h_{3}}The divergence of a vector field can then be written as∇ ⋅ v = 1 h 1 h 2 h 3 ∂ ∂ q i ( h 1 h 2 h 3 v i) {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\cfrac {1}{h_{1}h_{2}h_{3}}}~{\frac {\partial }{\partial q^{i}}}(h_{1}h_{2}h_{3}~v^{i})}Also,v i = g i k v k ⇒ v 1 = g 11 v 1 = v 1 h 1 2 ; v 2 = g 22 v 2 = v 2 h 2 2 ; v 3 = g 33 v 3 = v 3 h 3 2 {\displaystyle v^{i}=g^{ik}~v_{k}\quad \Rightarrow v^{1}=g^{11}~v_{1}={\cfrac {v_{1}}{h_{1}^{2}}}~;~~v^{2}=g^{22}~v_{2}={\cfrac {v_{2}}{h_{2}^{2}}}~;~~v^{3}=g^{33}~v_{3}={\cfrac {v_{3}}{h_{3}^{2}}}}Therefore,∇ ⋅ v = 1 h 1 h 2 h 3 ∑ i ∂ ∂ q i ( h 1 h 2 h 3 h i 2 v i) {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\cfrac {1}{h_{1}h_{2}h_{3}}}~\sum _{i}{\frac {\partial }{\partial q^{i}}}\left({\cfrac {h _{1}h_{2}h_{3}}{h_{i}^{2}}}~v_{i}\right)}We can get an expression for the Laplacian in a similar manner by noting thatg l i ∂ φ ∂ q l = { g 11 ∂ φ ∂ q 1, g 22 ∂ φ ∂ q 2, g 33 ∂ φ ∂ q 3 } = { 1 h 1 2 ∂ φ ∂ q 1, 1 h 2 2 ∂ φ ∂ q 2, 1 h 3 2 ∂ φ ∂ q 3 } {\displaystyle g^{li}~{\frac {\partial \varphi }{\partial q^{l}}}=\left\{g^{11}~{\frac {\partial \varphi }{\partialq^{1}}},g^{22}~{\frac {\partial \varphi }{\partial q^{2}}},g^{33}~{\frac {\partial \varphi }{\partial q^{3}}}\right\}=\left\{{\cfrac {1}{h_{1}^{2}}}~{\frac {\partial \varphi }{\partial q^{1}}},{\cfrac {1}{h_{2}^{2}}}~{\frac {\partial \varphi }{\partial q^{2}}},{\cfrac {1}{h_{3}^{2}}}~{\frac {\partial \varphi }{\partial q^{3}}}\right\}}Then we have∇ 2 φ = 1 h 1 h 2 h 3 ∑ i ∂ ∂ q i ( h 1 h 2 h 3 h i 2∂ φ ∂ q i) {\displaystyle \nabla ^{2}\varphi ={\cfrac {1}{h_{1}h_{2}h_{3}}}~\sum _{i}{\frac {\partial }{\partial q^{i}}}\left({\cfrac {h_{1}h_{2}h_{3}}{h_{i}^{2}}}~{\frac {\partial \varphi }{\partial q^{i}}}\right)}The expressions for the gradient, divergence, and Laplacian can be directly extended to n-dimensions. The curl of a vector field is given by∇ × v = 1 h 1 h 2 h 3 ∑ i = 1 n e i ∑ j k ε i j k h i ∂ ( h k v k) ∂ q j {\displaystyle \nabla \times \mathbf {v} ={\frac {1}{h_{1}h_{2}h_{3}}}\sum _{i=1}^{n}\mathbf {e} _{i}\sum _{jk}\varepsilon _{ijk}h_{i}{\frac {\partial (h_{k}v_{k})}{\partial q^{j}}}}where εijkis the Levi-Civita symbol.Example: Cylindrical polar coordinatesFor cylindrical coordinates we have( x 1, x 2, x 3) = x = φ ( q 1, q 2, q 3) = φ ( r, θ, z) = { r cos θ, r sin θ, z } {\displaystyle (x_{1},x_{2},x_{3})=\mathbf {x} ={\boldsymbol {\varphi }}(q^{1},q^{2},q^{3})={\boldsymbol {\varphi }}(r,\theta,z)=\{r\cos \theta,r\sin \theta,z\}}and{ ψ 1 ( x), ψ 2 ( x), ψ 3 ( x) } = ( q 1, q 2, q 3) ≡ ( r, θ, z) = { x 1 2 + x 2 2, tan − 1 (Икс 2 / Икс 1), Икс 3} {\ Displaystyle \ {\ psi ^ {1} (\ mathbf {x}), \ psi ^ {2} (\ mathbf {x}), \ psi ^ {3 } (\ mathbf {x}) \} = (q ^ {1}, q ^ {2}, q ^ {3}) \ эквив (r, \ theta, z) = \ {{\ sqrt {x_ {1 } ^ {2} + x_ {2} ^ {2}}}, \ tan ^ {- 1} (x_ {2} / x_ {1}), x_ {3} \}}где0 < r < ∞, 0 < θ < 2 π, − ∞ < z < ∞ {\displaystyle 0Тогда ковариантный и контравариантный базисные векторы равныb 1 = er = b 1 b 2 = re θ = r 2 b 2 b 3 = ez = b 3 {\ displaystyle {\ begin {align} \ mathbf {b} _ {1} = \ mathbf {e} _ {r} = \ mathbf {b} ^ {1} \\\ mathbf {b} _ {2} = r ~ \ mathbf {e} _ {\ theta} = r ^ {2} ~ \ mathbf {b} ^ {2} \\\ mathbf {b} _ {3} = \ mathbf {e} _ {z} = \ mathbf {b} ^ {3} \ end { выровнено}}}где er, e θ, ez {\ displaystyle \ mathbf {e} _ {r}, \ mathbf {e} _ {\ theta}, \ mathbf {e} _ {z} }- единичные векторы в направлениях r, θ, z {\ displaystyle r, \ theta, z}. Обратите внимание, что компоненты метрического тензора таковы, чтоg i j = g i j = 0 (i ≠ j); g 11 = 1, g 22 = 1 r, g 33 = 1 {\ displaystyle g ^ {ij} = g_ {ij} = 0 (i \ neq j) ~; ~~ {\ sqrt {g ^ {11}} } = 1, ~ {\ sqrt {g ^ {22}}} = {\ cfrac {1} {r}}, ~ {\ sqrt {g ^ {33}}} = 1}, что показывает, что базис ортогонален. Ненулевые компоненты символа Кристоффеля второго типа равныΓ 12 2 = Γ 21 2 = 1 r; Γ 22 1 = - р {\ displaystyle \ Gamma _ {12} ^ {2} = \ Gamma _ {21} ^ {2} = {\ cfrac {1} {r}} ~; ~~ \ Gamma _ {22 } ^ {1} = - r}Представление физического векторного поляНормализованные контрвариантные базисные векторы в цилиндрических полярных координатах:b ^ 1 = er; Ь ^ 2 = е θ; b ^ 3 = ez {\ displaystyle {\ hat {\ mathbf {b}}} ^ {1} = \ mathbf {e} _ {r} ~; ~~ {\ hat {\ mathbf {b}}} ^ { 2} = \ mathbf {e} _ {\ theta} ~; ~~ {\ hat {\ mathbf {b}}} ^ {3} = \ mathbf {e} _ {z}}и физические компоненты вектора v равны(v ^ 1, v ^ 2, v ^ 3) = (v 1, v 2 / r, v 3) =: (vr, v θ, vz) {\ displaystyle ({\ hat {v}} _ {1}, {\ hat {v}} _ {2}, {\ hat {v}} _ {3}) = (v_ {1}, v_ {2 } / r, v_ {3}) = :( v_ {r}, v _ {\ theta}, v_ {z})}Градиент скалярного поляГрадиент скалярного поля, f (x ), в цилиндрических координатах может теперь вычисляется из общего выражения в криволинейных координатах и имеет вид∇ f = ∂ f ∂ rer + 1 r ∂ f ∂ θ e θ + ∂ f ∂ zez {\ displaystyle {\ boldsymbol {\ nabla}} f = {\ cfrac {\ partial f} {\ partial r}} ~ \ mathbf {e} _ {r} + {\ cfrac {1} {r}} ~ {\ cfrac {\ partial f} {\ partial \ theta }} ~ \ mathbf {e} _ {\ theta} + {\ cfrac {\ partial f} {\ partial z}} ~ \ mathbf {e} _ {z}}Градиент вектора или поляАналогично, градиент векторного поля, v(x), в цилиндрических координатах может быть показан как∇ v = ∂ vr ∂ rer ⊗ er + 1 r (∂ vr ∂ θ - v θ) er ⊗ e θ + ∂ vr ∂ zer ⊗ ez + ∂ v θ ∂ re θ ⊗ er + 1 r (∂ v θ ∂ θ + vr) e θ ⊗ e θ + ∂ v θ ∂ ze θ ⊗ ez + ∂ vz ∂ rez ⊗ er + 1 р ∂ vz ∂ θ ez ⊗ e θ + ∂ vz ∂ zez ⊗ ez {\ displaystyle {\ begin {align} {\ boldsymbol {\ nabla}} \ mathbf {v} = {\ cfrac {\ частичный v_ {r}} {\ partial r}} ~ \ mathbf {e} _ {r} \ otimes \ mathbf {e} _ {r} + {\ cfrac {1} {r}} \ left ({\ cfrac {\ partial v_ {r}} {\ partial \ theta}} - v _ {\ theta} \ right) ~ \ mathbf {e} _ {r} \ otimes \ mathbf {e} _ {\ theta} + {\ cfrac {\ partial v_ {r}} {\ partial z}} ~ \ mathbf {e} _ {r} \ otimes \ mathbf {e} _ {z} \\ [8pt] + {\ cfrac {\ partial v _ {\ theta}} {\ partial r}} ~ \ mathbf {e} _ {\ theta} \ otimes \ mathbf {e} _ {r} + {\ cfrac { 1} {r}} \ left ({\ cfrac {\ partial v _ {\ theta}} {\ partial \ theta}} + v_ {r} \ right) ~ \ mathbf {e} _ {\ theta} \ otimes \ mathbf {e} _ {\ theta} + {\ cfrac {\ partial v _ {\ theta}} {\ partial z}} ~ \ mathbf {e} _ {\ theta} \ otimes \ mathbf {e} _ {z} \\ [8pt] + {\ cfrac {\ partial v_ {z}} {\ partial r}} ~ \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {r} + {\ cfrac { 1} {r}} {\ cfrac {\ partial v_ {z}} {\ partial \ theta}} ~ \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {\ theta} + {\ cfrac {\ partial v_ {z}} {\ partial z}} ~ \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {z} \ end {align}}}Дивергенция векторного поляИспользуя уравнение для расходимости векторного поля в криволинейных координатах, можно показать, что расхождение в цилиндрических координатах составляет∇ ⋅ v = ∂ vr ∂ r + 1 r (∂ v θ ∂ θ + vr) + ∂ vz ∂ Z {\ displaystyle {\ begin {align} {\ bo л dsymbol {\ nabla}} \ cdot \ mathbf {v} = {\ cfrac {\ partial v_ {r}} {\ partial r}} + {\ cfrac {1} {r}} \ left ({\ cfrac { \ partial v _ {\ theta}} {\ partial \ theta}} + v_ {r} \ right) + {\ cfrac {\ partial v_ {z}} {\ partial z}} \ end {align}}}Лапласиан скалярного поляЛапласиан легче вычислить, если учесть, что ∇ 2 f = ∇ ⋅ ∇ f {\ displaystyle {\ boldsymbol {\ nabla}} ^ {2} f = { \ boldsymbol {\ nabla}} \ cdot {\ boldsymbol {\ nabla}} f}. В цилиндрических полярных координатахv = ∇ f = [vrv θ vz] = [∂ f ∂ r 1 r ∂ f ∂ θ ∂ f ∂ z] {\ displaystyle \ mathbf {v} = {\ boldsymbol {\ nabla} } f = \ left [v_ {r} ~~ v _ {\ theta} ~~ v_ {z} \ right] = \ left [{\ cfrac {\ partial f} {\ partial r}} ~~ {\ cfrac { 1} {r}} {\ cfrac {\ partial f} {\ partial \ theta}} ~~ {\ cfrac {\ partial f} {\ partial z}} \ right]}Следовательно,∇ ⋅ v = ∇ 2 f = ∂ 2 f ∂ r 2 + 1 r (1 r ∂ 2 f ∂ θ 2 + ∂ f ∂ r) + ∂ 2 f ∂ z 2 = 1 r [∂ ∂ r (r ∂ f ∂ r)] + 1 р 2 ∂ 2 е ∂ θ 2 + ∂ 2 е ∂ Z 2 {\ displaystyle {\ boldsymbol {\ nabla}} \ cdot \ mathbf {v} = {\ boldsymbol {\ nabla}} ^ {2 } f = {\ cfrac {\ partial ^ {2} f} {\ partial r ^ {2}}} + {\ cfrac {1} {r}} \ left ({\ cfrac {1} {r}} { \ cfrac {\ partial ^ {2} f} {\ partial \ theta ^ {2}}} + {\ cfrac {\ partial f} {\ partial r}} \ right) + {\ cfrac {\ partial ^ {2 } f} {\ partial z ^ {2}}} = {\ cfrac {1} {r}} \ left [{\ cfrac {\ partial} {\ partial r}} \ left (r {\ cfrac {\ partial f} {\ partial r}} \ right) \ right] + {\ cfrac {1} {r ^ {2}}} {\ cfrac {\ partial ^ {2} f} {\ partial \ theta ^ {2} }} + {\ cfrac {\ partial ^ {2} f} {\ partial z ^ {2}}} }Представление физического тензорного поля второго порядкаФизические компоненты тензорного поля второго порядка - это те, которые получаются, когда тензор выражается в терминах нормированного контравариантного базиса. В цилиндрических полярных координатах эти компоненты следующие:S ^ 11 = S 11 =: S rr, S ^ 12 = S 12 r =: S r θ, S ^ 13 = S 13 =: S rz S ^ 21 = S 21 r =: S θ r, S ^ 22 = S 22 r 2 =: S θ θ, S ^ 23 = S 23 r =: S θ z S ^ 31 = S 31 =: S zr, S ^ 32 = S 32 р знак равно: S z θ, S ^ 33 = S 33 =: S zz {\ displaystyle {\ begin {align} {\ hat {S}} _ {11} = S_ {11} =: S_ {rr}, {\ hat {S}} _ {12} = {\ frac {S_ {12}} {r}} =: S_ {r \ theta}, {\ hat {S}} _ {13} = S_ {13} =: S_ {rz} \\ [6pt] {\ hat {S}} _ {21} = { \ frac {S_ {21}} {r}} =: S _ {\ theta r}, {\ hat {S}} _ {22} = {\ frac {S_ {22}} {r ^ {2} }} =: S _ {\ theta \ theta}, {\ hat {S}} _ {23} = {\ frac {S_ {23}} {r}} =: S _ {\ theta z} \\ [ 6pt] {\ hat {S}} _ {31} = S_ {31} =: S_ {zr}, {\ hat {S}} _ {32} = {\ frac {S_ {32}} { r}} =: S_ {z \ theta}, {\ hat {S}} _ {33} = S_ {33} =: S_ {zz} \ end {align}}}Градиент секунды тензорное поле -порядкаИспользуя приведенные выше определения, мы можем показать, что градиент тензорного поля второго порядка в цилиндрических полярных координатах может быть выражен как∇ S = ∂ S rr ∂ rer ⊗ er ⊗ er + 1 r [∂ S rr ∂ θ - (S θ r + S r θ)] er ⊗ er ⊗ e θ + ∂ S rr ∂ zer ⊗ er ⊗ ez + ∂ S r θ ∂ rer ⊗ e θ ⊗ er + 1 r [∂ S r θ ∂ θ + (S rr - S θ θ)] er ⊗ e θ ⊗ e θ + ∂ S r θ ∂ zer ⊗ e θ ez + ∂ S rz ∂ rer ⊗ ez er + 1 r [∂ S rz ∂ θ - S θ z] er ⊗ ez ⊗ e θ + ∂ S rz ∂ zer ⊗ ez ⊗ ez + ∂ S θ r ∂ re θ ⊗ er ⊗ er + 1 r [∂ S θ r ∂ θ + (S rr - S θ θ)] e θ ⊗ er ⊗ e θ + ∂ S θ r ∂ ze θ ⊗ er ⊗ ez + ∂ S θ θ ∂ re θ ⊗ e θ ⊗ er + 1 r [∂ S θ θ ∂ θ + (S r θ + S θ r)] e θ ⊗ e θ ⊗ e θ + ∂ S θ θ ∂ ze θ ⊗ e θ ⊗ ez + ∂ S θ z ∂ re θ ⊗ ez ⊗ er + 1 r [∂ S θ z ∂ θ + S rz] e θ ⊗ ez ⊗ e θ + ∂ S θ z ∂ ze θ ⊗ ez ⊗ ez + ∂ S zr ∂ rez ⊗ er ⊗ er + 1 r [∂ S zr ∂ θ - S z θ] ez ⊗ er ⊗ e θ + ∂ S zr ∂ zez ⊗ er ⊗ ez + ∂ S z θ ∂ rez ⊗ e θ ⊗ er + 1 r [∂ S z θ ∂ θ + S zr] ez ⊗ e θ ⊗ e θ + ∂ S z θ ∂ zez ⊗ e θ ⊗ ez + ∂ S zz ∂ rez ⊗ ez ⊗ er + 1 r ∂ S zz ∂ θ ez ⊗ ez ⊗ e θ + ∂ S zz ∂ zez ⊗ ez ⊗ ez {\ displaystyle {\ begin {align} {\ boldsymbol {\ nabla}} {\ boldsymbol {S}} = {\ frac {\ partial S_ {rr}} {\ partial r}} ~ \ mathbf {e} _ {r} \ otimes \ mathbf {e} _ {r} \ otimes \ mathbf {e} _ {r} + {\ cfrac {1} {r}} \ left [{\ frac { \ partial S_ {rr}} {\ partial \ theta}} - (S _ {\ theta r} + S_ {r \ theta}) \ right] ~ \ mathbf {e} _ {r} \ otimes \ mathbf {e} _ {r} \ otimes \ mathbf {e} _ {\ theta} + {\ frac {\ partial S_ {rr}} {\ partial z}}~\mathbf {e} _{r}\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{z}\\[8pt]+{\frac {\partial S_{r\theta }}{\partial r}}~\mathbf {e} _{r}\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{r\theta }}{\partial \theta }}+(S_{rr}-S_{\theta \theta })\right]~\mathbf {e} _{r}\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }+{\frac{\partial S_{r\theta }}{\partial z}}~\mathbf {e} _{r}\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{z}\\[8pt]+{\frac {\partial S_{rz}}{\partial r}}~\mathbf {e} _{r}\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{rz}}{\partial \theta }}-S_{\theta z}\right]~\mathbf {e} _{r}\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{rz}}{\partial z}}~\mathbf {e} _{r}\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{z}\\[8pt]+{\frac {\partial S_{\theta r}}{\partial r}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{\theta r}}{\partial \theta }}+(S_{rr}-S_{\theta \theta })\right]~\mathbf {e} _{\theta }\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{\theta r}}{\partial z}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{z}\\[8pt]+{\frac {\partial S_{\theta \theta }}{\partial r}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{\theta \theta }}{\partial \theta }}+(S_{r\theta }+S_{\theta r})\right]~\mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{\theta \theta }}{\partial z}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{z}\\[8pt]+{\frac {\partial S_{\theta z}}{\partial r}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{\theta z}}{\partial \theta }}+S_{rz}\right]~\mathbf {e} _{\theta }\otimes \mathbf {e} _{z}\otimes\ mathbf {e} _ {\ theta} + {\ frac {\ partial S _ {\ theta z}} {\ partial z}} ~ \ mathbf {e} _ {\ theta} \ otimes \ mathbf {e} _ { z} \ otimes \ mathbf {e} _ {z} \\ [8pt] + {\ frac {\ partial S_ {zr}} {\ partial r}} ~ \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {r} \ otimes \ mathbf {e} _ {r} + {\ cfrac {1} {r}} \ left [{\ frac {\ partial S_ {zr}} {\ partial \ theta} } -S_ {z \ theta} \ right] ~ \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {r} \ otimes \ mathbf {e} _ {\ theta} + {\ frac {\ частичный S_ {zr}} {\ partial z}} ~ \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {r} \ otimes \ mathbf {e} _ {z} \\ [8pt] + {\ frac {\ partial S_ {z \ theta}} {\ partial r}} ~ \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {\ theta} \ otimes \ mathbf {e} _ {r} + {\ cfrac {1} {r}} \ left [{\ frac {\ partial S_ {z \ theta}} {\ partial \ theta}} + S_ {zr} \ right] ~ \ mathbf {e } _ {z} \ otimes \ mathbf {e} _ {\ theta} \ otimes \ mathbf {e} _ {\ theta} + {\ frac {\ partial S_ {z \ theta}} {\ partial z}} ~ \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {\ theta} \ otimes \ mathbf {e} _ {z} \\ [8pt] + {\ frac {\ partial S_ {zz}} {\ partial r}} ~ \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {z} \ oti mes \ mathbf {e} _ {r} + {\ cfrac {1} {r}} ~ {\ frac {\ partial S_ {zz}} {\ partial \ theta}} ~ \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {\ theta} + {\ frac {\ partial S_ {zz}} {\ partial z}} ~ \ mathbf {e} _ {z } \ otimes \ mathbf {e} _ {z} \ otimes \ mathbf {e} _ {z} \ end {align}}}Дивергенция тензорного поля второго порядкаДивергенция тензорное поле второго порядка в цилиндрических полярных координатах может быть получено из выражения для градиента путем сбора членов, в которых скалярное произведение двух внешних векторов в диадических произведениях отлично от нуля. Следовательно,∇ ⋅ S = ∂ S rr ∂ rer + ∂ S r θ ∂ re θ + ∂ S rz ∂ rez + 1 r [∂ S r θ ∂ θ + (S rr - S θ θ)] er + 1 r [∂ S θ θ ∂ θ + (S r θ + S θ r)] e θ + 1 r [∂ S θ z ∂ θ + S rz] ez + ∂ S zr ∂ zer + ∂ S z θ ∂ ze θ + ∂ S zz ∂ zez {\ displaystyle {\ begin {align} {\ boldsymbol {\ nabla}} \ cdot {\ boldsymbol {S}} = {\ frac {\ partial S_ {rr}} {\ partial r }} ~ \ mathbf {e} _ {r} + {\ frac {\ partial S_ {r \ theta}} {\ partial r}} ~ \ m athbf {e} _{\theta }+{\frac {\partial S_{rz}}{\partial r}}~\mathbf {e} _{z}\\[8pt]+{\cfrac {1}{r}}\left[{\frac {\partial S_{r\theta }}{\partial \theta }}+(S_{rr}-S_{\theta \theta })\right]~\mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{\theta \theta }}{\partial \theta }}+(S_{r\theta }+S_{\theta r})\right]~\mathbf {e} _{\theta }+{\cfrac {1}{r}}\left[{\frac {\partialS_{\theta z}}{\partial \theta }}+S_{rz}\right]~\mathbf {e} _{z}\\[8pt]+{\frac {\partial S_{zr}}{\partial z}}~\mathbf {e} _{r}+{\frac {\partial S_{z\theta }}{\partial z}}~\mathbf {e} _{\theta }+{\frac {\partial S_{zz}}{\partial z}}~\mathbf {e} _{z}\end{aligned}}}See alsoCovariance and contravariance Basic introduction to the mathematics of curved spacetime Orthogonal coordinates Frenet–Serret formulas Covariant derivative Tensor derivative (continuum mechanics) Curvilinear perspective Del in cylindrical and spherical coordinates ReferencesNotesFurther readingExternal linksDerivation of Unit Vectors in Curvilinear Coordinates MathWorld's page on Curvilinear Coordinates Prof. R. Brannon's E-Book on Curvilinear Coordinates
Let (b1, b2, b3) be an arbitrary basis for three-dimensional Euclidean space. In general, the basis vectors are neither unit vectors nor mutually orthogonal. However, they are required to be linearly independent. Then a vector vcan be expressed as
The components v are the contravariantcomponents of the vector v.
The reciprocal basis(b, b, b) is defined by the relation
where δjis the Kronecker delta.
The vector vcan also be expressed in terms of the reciprocal basis:
The components vkare the covariantcomponents of the vector v {\displaystyle \mathbf {v} }.
A second-order tensor can be expressed as
The components S are called the contravariantcomponents, Sjthe mixed right-covariantcomponents, Sithe mixed left-covariantcomponents, and Sijthe covariantcomponents of the second-order tensor.
The quantities gij, g are defined as
From the above equations we have
The components of a vector are related by
Also,
The components of the second-order tensor are related by
In an orthonormal right-handed basis, the third-order alternating tensor is defined as
In a general curvilinear basis the same tensor may be expressed as
It can be shown that
Now,
Hence,
Similarly, we can show that
The identity map Idefined by I ⋅ v = v {\displaystyle \mathbf {I} \cdot \mathbf {v} =\mathbf {v} }can be shown to be:
The scalar product of two vectors in curvilinear coordinates is
The cross product of two vectors is given by:
where εijkis the permutation symbol and eiis a Cartesian basis vector. In curvilinear coordinates, the equivalent expression is:
where E i j k{\displaystyle {\mathcal {E}}_{ijk}}is the third-order alternating tensor. The cross product of two vectors is given by:
where εijkis the permutation symbol and e i {\displaystyle \mathbf {e} _{i}}is a Cartesian basis vector. Therefore,
and
Returning to the vector product andusing the relations:
gives us:
The identity map I {\displaystyle {\mathsf {I}}}defined by I ⋅ v = v {\displaystyle {\mathsf {I}}\cdot \mathbf {v} =\mathbf {v} }can be shown to be
The action v = S ⋅ u {\displaystyle \mathbf {v} ={\boldsymbol {S}}\cdot \mathbf {u} }can be expressed in curvilinear coordinates as
The inner product of two second-order tEnsors U = S ⋅ T {\ displaystyle {\ boldsymbol {U}} = {\ boldsymbol {S}} \ cdot {\ boldsymbol {T}}}может быть выражено в криволинейных координатах как
В качестве альтернативы
Если S {\ displaystyle {\ boldsymbol {S}}}- тензор второго порядка, тогда определитель определяется соотношением
где u, v, w {\ displaystyle \ mathbf {u}, \ mathbf {v}, \ mathbf {w}}- произвольные векторы, а
Пусть (e1, e2, e3) - обычные декартовы базисные векторы для интересующего евклидова пространства, и пусть
, где Fi- тензор преобразования второго порядка, который отображает eiв bi. Тогда
Из этого соотношения мы можем показать, что
Пусть J: = det F {\ displaystyle J: = \ det {\ boldsymbol {F}}}будет Якобиан преобразования. Тогда, согласно определению определителя,
Поскольку
мы иметь
С помощью приведенных выше соотношений можно получить ряд интересных результатов.
Сначала рассмотрим
Затем
Аналогичным образом мы можем показать, что
Следовательно, используя тот факт, что [gij] = [gij] - 1 {\ displaystyle [g ^ {ij}] = [g_ {ij}] ^ {- 1}},
Еще одно интересное соотношение выводится ниже. Напомним, что
Из этих результатов имеем
и
Симмондс в своей книге по тензорному анализу цитирует Альберта Эйнштейна высказывание
Магия этой теории вряд ли перестанет навязываться любому, кто действительно понимает od это; он представляет собой подлинный триумф метода абсолютного дифференциального исчисления, основанного Гауссом, Риманом, Риччи и Леви-Чивита.
Векторное и тензорное исчисление в общих криволинейных координатах используется в тензорном анализе четырехмерных криволинейных многообразий в общей теории относительности, в механике криволинейных оболочки, при исследовании свойств инвариантности уравнений Максвелла, которые представляли интерес в метаматериалах и во многих других областях.
В этом разделе приведены некоторые полезные соотношения в исчислении векторов и тензоров второго порядка в криволинейных координатах. ион. Обозначения и содержание в основном взяты из Огдена, Симмондса, Грина и Зерны, Басара и Вейхерта и Чиарлета.
Пусть положение точки в пространстве характеризуется тремя координатными переменными (q 1, q 2, q 3) {\ displaystyle (q ^ {1}, q ^ {2}, q ^ {3})}.
Координатная кривая q представляет кривая, на которой q, q постоянны. Пусть x будет вектором положения точки относительно некоторого начала. Тогда, предполагая, что такое отображение и обратное ему отображение существуют и непрерывны, мы можем записать
Поля ψ (x ) называются функциями криволинейных координат в криволинейной системе координат ψ(x) = φ(x).
Координатные кривые q определяются однопараметрическим семейством функций, задаваемым
с фиксированными q, q.
Касательный вектор к кривой xiв точке xi(α) (или к координатной кривой q я в точке x ) равно
Пусть f (x ) будет скалярным полем в пространстве. Тогда
Градиент поля f определяется как
где c - произвольный постоянный вектор. Если мы определим компоненты c из c таковы, что
, затем
Если мы положим f (x) = ψ i ( x) {\ displaystyle f (\ mathbf {x}) = \ psi ^ {i} (\ mathbf {x})}, то, поскольку qi = ψ i (x) {\ displaystyle q ^ {i} = \ psi ^ {i} (\ mathbf {x})}, имеем
, который обеспечивает средства извлечения контравариантной компоненты вектора c.
Если biявляется ковариантным (или естественным) базисом в точке, и если b является контравариантным (или обратным) базисом в этой точке, то
Краткое обоснование этого выбора базиса дается в следующем разделе..
Аналогичный процесс можно использовать для получения градиента a vector field f(x). The gradient is given by
If we consider the gradient of the position vector field r(x) = x, then we can show that
The vector field biis tangent to the q coordinate curve and forms a natural basisat each point on the curve. This basis, as discussed at the beginning of this article, is also called the covariantcurvilinear basis. We can also define a reciprocal basis, or contravariantcurvilinear basis, b. 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.
Since cis arbitrary, we can write
Note that the contravariant basis vector bis perpendicular to the surface of constant ψ and is given by
The Christoffel symbols of the first kind are defined as
To express Γijkin terms of gijwe note that
Since b i,j= bj,iwe have Γijk= Γjik. Using these to rearrange the above relations gives
The Christoffel symbols of the second kind are defined as
in which
This implies that
Other relations that follow are
Еще одно особенно полезное соотношение, которое показывает, что символ Кристоффеля зависит от l y от метрического тензора и его производных, это
Следующие выражения для градиента векторного поля в криволинейных координатах весьма полезны.
Векторное поле v может быть представлено как
где vi { \ displaystyle v_ {i}}- ковариантные компоненты поля, v ^ i {\ displaystyle {\ hat {v}} _ {i}}- физические компоненты и (без суммирования )
- нормализованный контравариантный базис v эктор.
Градиент тензорного поля второго порядка может аналогичным образом быть выражено как
Если мы рассмотрим выражение для тензора в терминах контравариантного базиса, то
Мы также можем написать
Физические компоненты тензорного поля второго порядка могут быть получены с использованием нормированного контравариантного базиса, т. е.
, где штрихованные базисные векторы нормализованы. Это означает, что (опять же без суммирования)
Дивергенция векторного поля (v {\ displaystyle \ mathbf {v}}) определяется как
В терминах компонентов относительно криволинейного базиса
Альтернативным уравнением для расходимости векторного поля является часто используемый. Чтобы вывести это соотношение, напомним, что
Теперь
Отмечая, что из-за симметрии из g {\ displaystyle {\ boldsymbol {g}}},
Recall that if [gij] is the matrix whose components are gij, then the inverse of the matrix is [ g i j ] − 1 = [ g i j ] {\displaystyle [g_{ij}]^{-1}=[g^{ij}]}. The inverse of the matrix is given by
where A are the Cofactor matrix of the components gij. From matrix algebra we have
Plugging this relation into the expression for the divergence gives
A little manipulation leads to the more compact form
The divergence of a second-order tensor field is defined using
where ais an arbitrary constant vector. In curvilinear coordinates,
Лапласиан скалярного поля φ (x ) определяется как
Использование альтернативного выражения для дивергенции векторного поля дает нам
Следовательно,
The curl of a vector field vin covariant curvilinear coordinates can be written as
where
Assume, for thepurposes of this section, that the curvilinear coordinate system is orthogonal, i.e.,
or equivalently,
where g i i = g i i − 1 {\displaystyle g^{ii}=g_{ii}^{-1}}. As before, b i, b j {\displaystyle \mathbf {b} _{i},\mathbf {b} _{j}}are covariant basis vectors and b, bare contravariant basis vectors. Also, let (e, e, e) be a background, fixed, Cartesian basis. A list of orthogonal curvilinear coordinates is given below.
Let r(x) be the position vector of the point xwith respect to the origin of the coordinate system. The notation can be simplified by noting that x= r(x). At each point we can construct a small line element dx. The square of the length of the line element is the scalar product dx• dxand is called the metric of the space. Recall that the space of interest is assumed to be Euclidean when we talk of curvilinear coordinates. Let us express the position vector in terms of the background, fixed, Cartesian basis, i.e.,
Using the chain rule, we can then express dxin terms of three-dimensional orthogonal curvilinear coordinates (q, q, q) as
Therefore, the metric is given by
The symmetric quantity
is called the fundamental (or metric) tensor of the Euclidean space in curvilinear coordinates.
Note also that
where hijare the Lamé coefficients.
If we define the scale factors, hi, using
we get a relation between the fundamental tensor and the Lamé coefficients.
Если мы рассмотрим полярные координаты для R, обратите внимание, что
(r, θ) - криволинейные координаты, а определитель Якоби преобразования (r, θ) → (r cos θ, r sin θ) равно r.
Ортогональные базисные векторы равны br= (cos θ, sin θ), bθ= (−r sin θ, r cos θ). Нормализованные базисные векторы: er= (cos θ, sin θ), eθ= (−sin θ, cos θ), а масштабные коэффициенты - h r = 1 и h θ. = р. Фундаментальный тензор равен g 11 = 1, g 22 = r, g 12 = g 21 = 0.
Если мы хотим использовать криволинейные координаты для вычислений векторного исчисления, необходимо внести корректировки в вычисление линейных, поверхностных и объемных интегралов. Для простоты мы снова ограничим обсуждение тремя измерениями и ортогональными криволинейными координатами. Однако те же аргументы применимы к n {\ displaystyle n}-мерным задачам, хотя в выражениях есть некоторые дополнительные члены, когда система координат не ортогональна.
Обычно при вычислении линейных интегралов нас интересует вычисление
где x (t) параметризует C в декартовых координатах. В криволинейных координатах член
по правилу цепочки . И из определения коэффициентов Ламе,
and thus
Now, since g i j = 0 {\displaystyle g_{ij}=0}when i ≠ j {\displaystyle i\neq j}, we have
and we can proceed normally.
Likewise, if we are interested in a surface integral, the relevant calculation, with the parameterization of the surface in Cartesian coordinates is:
Again, in curvilinear coordinates, we have
и снова используем определение криволинейных координат, чтобы получить
где E { \ displaystyle {\ mathcal {E}}}- это символ перестановки .
В детерминантной форме перекрестное произведение с точки зрения криволинейных координат будет:
In orthogonal curvilinear coordinates of 3 dimensions, where
one can express the gradient of a scalar or vector field as
For an orthogonal basis
The divergence of a vector field can then be written as
Therefore,
We can get an expression for the Laplacian in a similar manner by noting that
Then we have
The expressions for the gradient, divergence, and Laplacian can be directly extended to n-dimensions.
The curl of a vector field is given by
where εijkis the Levi-Civita symbol.
For cylindrical coordinates we have
где
Тогда ковариантный и контравариантный базисные векторы равны
где er, e θ, ez {\ displaystyle \ mathbf {e} _ {r}, \ mathbf {e} _ {\ theta}, \ mathbf {e} _ {z} }- единичные векторы в направлениях r, θ, z {\ displaystyle r, \ theta, z}.
Обратите внимание, что компоненты метрического тензора таковы, что
, что показывает, что базис ортогонален.
Ненулевые компоненты символа Кристоффеля второго типа равны
Нормализованные контрвариантные базисные векторы в цилиндрических полярных координатах:
и физические компоненты вектора v равны
Градиент скалярного поля, f (x ), в цилиндрических координатах может теперь вычисляется из общего выражения в криволинейных координатах и имеет вид
Аналогично, градиент векторного поля, v(x), в цилиндрических координатах может быть показан как
Используя уравнение для расходимости векторного поля в криволинейных координатах, можно показать, что расхождение в цилиндрических координатах составляет
Лапласиан легче вычислить, если учесть, что ∇ 2 f = ∇ ⋅ ∇ f {\ displaystyle {\ boldsymbol {\ nabla}} ^ {2} f = { \ boldsymbol {\ nabla}} \ cdot {\ boldsymbol {\ nabla}} f}. В цилиндрических полярных координатах
Физические компоненты тензорного поля второго порядка - это те, которые получаются, когда тензор выражается в терминах нормированного контравариантного базиса. В цилиндрических полярных координатах эти компоненты следующие:
Используя приведенные выше определения, мы можем показать, что градиент тензорного поля второго порядка в цилиндрических полярных координатах может быть выражен как
Дивергенция тензорное поле второго порядка в цилиндрических полярных координатах может быть получено из выражения для градиента путем сбора членов, в которых скалярное произведение двух внешних векторов в диадических произведениях отлично от нуля. Следовательно,