Тензоры в криволинейных координатах

редактировать

Криволинейные координаты могут быть сформулированы в тензорном исчислении с важными приложениями в физика и инженерия, особенно для описания перемещения физических величин и деформации материи в механике жидкости и механике сплошных сред.

Содержание
  • 1 Векторная и тензорная алгебра в трехмерных криволинейных координатах
    • 1.1 Преобразования координат
      • 1.1.1 Якобиан
    • 1.2 Векторы в криволинейных координатах
    • 1.3 Тензоры второго порядка в криволинейных координатах
      • 1.3.1 Метрика тензор и отношения между компонентами
    • 1.4 Альтернативный тензор
    • 1.5 Векторные операции
      • 1.5.1 Карта идентичности
      • 1.5.2 Скалярное (точечное) произведение
      • 1.5.3 Векторное (перекрестное) произведение
    • 1.6 Тензорные операции
      • 1.6.1 Идентификационная карта
      • 1.6.2 Действие тензора второго порядка на вектор
      • 1.6.3 Внутреннее произведение двух тензоров второго порядка
      • 1.6.4 Определитель тензора второго порядка
    • 1.7 Связь между криволинейными и декартовыми базисными векторами
  • 2 Векторное и тензорное исчисление в трехмерных криволинейных координатах
    • 2.1 Основные определения
    • 2.2 Касательный вектор к координатные кривые
    • 2.3 Градиент
      • 2.3.1 Скалярное поле
      • 2.3.2 Векторное поле
      • 2.3.3 Символы Кристоффеля первого рода
      • 2.3.4 Символы Кристоффеля второго рода
      • 2.3.5 Явное выражение для градиента векторного поля
      • 2.3.6 Представление физического векторного поля
    • 2.4 Тензорное поле второго порядка
      • 2.4.1 Явное выражение для градиента
      • 2.4.2 Представление физическое тензорное поле второго порядка
    • 2.5 Дивергенция
      • 2.5.1 Векторное поле
      • 2.5.2 Тензорное поле второго порядка
    • 2.6 Лапласиан
      • 2.6.1 Скалярное поле
    • 2.7 Завиток векторное поле
  • 3 Ортогональные криволинейные координаты
    • 3.1 Метрический тензорв ортогональных криволинейных координатах
      • 3.1.1 Пример: полярные координаты
    • 3.2 Линейные и поверхностные интегралы
      • 3.2.1 Линейные интегралы
      • 3.2.2 Поверхностные интегралы
    • 3.3 Град, curl, div, лапласиан
  • 4 Пример: Цилиндрические полярные координаты
    • 4.1 Представление физического векторного поля
    • 4.2 Градиент скалярное поле
    • 4.3 Градиент векторного поля
    • 4.4 Дивергенция векторного поля
    • 4.5 Лапласиан скалярного поля
    • 4.6 Представление физического тензорного поля второго порядка
    • 4.7 Градиент секунды тензорное поле второго порядка
    • 4.8 Дивергенция тензорного поля второго порядка
  • 5 См. также
  • 6 Ссылки
  • 7 Внешние ссылки
Векторная и тензорная алгебра в трехмерных криволинейных координатах
Примечание : соглашение Эйнштейна о суммировании суммирования по повторяющимся индексам используется ниже.

Элементарная векторная и тензорная алгебра в криволинейных координатах используется в некоторой старой научной литературе по механике и физика и может быть незаменим для понимания работ начала и середины 1900-х годов, например текста Грина и Зерны. Некоторые полезные rel В этом разделе приведены элементы алгебры векторов и тензоров второго порядка в криволинейных координатах. Обозначения и содержание в основном взяты из Огдена, Нагди, Симмондса, Грина и Зерны, Басара и Вейхерта и Чиарлета.

Преобразования координат

Рассмотрим две системы координат с координатными переменными (Z 1, Z 2, Z 3) {\ displaystyle (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}})}{\displaystyle (Z^{\acute {1}},Z^{\acute {2}},Z^{\acute {3}})}, который мы будем представлять вкратце, просто Z i {\ displaystyle Z ^ {i}}{\displaystyle Z^{i}}и Z i ´ {\ displaystyle Z ^ {\ строго {i}}}{\displaystyle Z^{\acute {i}}}соответственно и всегда предполагаем, что наш индекс i {\ displaystyle i}iпроходит от 1 до 3. Мы будем предполагать,что эти системы координат встроены в трехмерное евклидово пространство. Могут использоваться координаты Z i {\ displaystyle Z ^ {i}}{\displaystyle Z^{i}}и Z i ´ {\ displaystyle Z ^ {\ строго {i}}}{\displaystyle Z^{\acute {i}}}. чтобы объяснить друг друга, потому что, перемещаясь по линии координат в одной системе координат, мы можем использовать другую для описания нашего положения. Таким образом, координаты Z i {\ displaystyle Z ^ {i}}{\displaystyle Z^{i}}и Z i ´ {\ displaystyle Z ^ {\ строго {i}}}{\displaystyle Z^{\acute {i}}}являются функциями друг друга

Z я = fi (Z 1 ´, Z 2 ´, Z 3 ´) {\ displaystyle Z ^ {i} = f ^ {i} (Z ^ {\ строго {1}}, Z ^ {\ острый {2}}, Z ^ {\ острый {3}})}{\displayst yle Z^{i}=f^{i}(Z^{\acute {1}},Z^{\acute {2}},Z^{\acute {3}})}для i = 1, 2, 3 {\ displaystyle 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}})}{\displaystyle Z^{i}=Z^{i}(Z^{\acute {1}},Z^{\acute {2}},Z^{\acute {3}})=Z^{i}(Z^{\acute {i}})}для i ´, i = 1, 2, 3 {\ displaystyle {\ строго {i}}, i = 1,2,3}{\displaystyle {\acute{i}},i=1,2,3}

Эти три уравнения вместе также называются координатами преобразование из Z i ´ {\ displaystyle Z ^ {\ строго {i}}}{\displaystyle Z^{\acute {i}}}в Z i {\ displaystyle Z ^ {i}}{\displaystyle Z^{i}}. мы обозначим это преобразование как T {\ displaystyle T}T. Поэтому мы представим преобразование из системы координат с координатными переменными Z i ´ {\ displaystyle Z ^ {\ строго {i}}}{\displaystyle Z^{\acute {i}}}в систему координат с координатами Z i { \ displaystyle Z ^ {i}}{\displaystyle Z^{i}}как:

Z = T (z ´) {\ displaystyle Z = T ({\ строго {z}})}{\displaystyle Z=T({\acute {z}})}

Аналогичным образом мы можем представить Z я ´ {\ displaystyle Z ^ {\ строго {i}}}{\displaystyle Z^{\acute {i}}}как функция от Z i {\ displaystyle Z ^ {i}}{\displaystyle Z^{i}}следующим образом :

Z я ´ = gi ´ (Z 1, Z 2, Z 3) {\ displaystyle Z ^ {\ строго {i}} = g ^ {\ строго {i}} (Z ^ {1}, Z ^ {2}, Z ^ {3})}{\displaystyle Z^{\acute {i}}=g^{\acute {i}}(Z^{1},Z^{2},Z^{3})}для i ´ = 1, 2, 3 {\ displaystyle {\ строго {i}} = 1,2,3}{\displaystyle {\acute {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})}{\displaystyle Z^{\acute {i}}=Z^{\acute {i}}(Z^{1},Z^{2},Z^{3})=Z^{\acute {i}}(Z^{i})}для i ´, i = 1, 2, 3 {\ displaystyle {\ строго {i}}, i = 1,2,3}{\displaystyle {\acute{i}},i=1,2,3}

Эти три уравнения вместе также называются преобразованием координат из Z i {\ displaystyle Z ^ {i}}{\displaystyle Z^{i}}до Z i ´ 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

v = v k b k {\displaystyle \mathbf {v} =v^{k}\,\mathbf {b} _{k}} \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 relation

b i ⋅ b j = δ j i {\displaystyle \mathbf {b} ^{i}\cdot \mathbf {b} _{j}=\delta _{j}^{i}} \mathbf{b}^i\cdot\mathbf{b}_j = \delta^i_j

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}} \mathbf{v} = v_k~\mathbf{b}^k

The components vkare the covariantcomponents of the vector v {\displaystyle \mathbf {v} }\mathbf {v} .

Second-order tensors in curvilinear coordinates

A second-order tensor can be expressed as

S = 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}} \boldsymbol{S} = S^{ij}~\mathbf{b}_i\otimes\mathbf{b}_j = S^{i}_{~j}~\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 components

The quantities gij, g are defined as

g 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}} 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 have

v 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}} 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 by

v ⋅ 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}} \mathbf{v}\cdot\mathbf{b}^i = v^k~\math bf{b}_k\cdot\mathbf{b}^i = v^k~\delta^i_k = 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}} \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}} \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}} \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 by

S 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}} S^{ij} = g^{ik}~S_k^{~j} = g^{jk}~S^i_{~k} = g^{ik}~g^{jl}~S_{kl}

The alternating tensor

In an orthonormal right-handed basis, the third-order alternating tensor is defined as

E = ε 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}} \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 as

E = 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}} \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 that

E 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]} \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}} \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}} \mathcal{E}_{ijk} = J~\varepsilon_{ijk} = \sqrt{g}~\varepsilon_{ijk}

Similarly, we can show that

E 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}} \mathcal{E}^{ijk} = \cfrac{1}{J}~\varepsilon^{ijk} = \cfrac{1}{\sqrt{g}}~\varepsilon^{ijk}

Vector operations

Identity map

The identity map Idefined by I ⋅ v = v {\displaystyle \mathbf {I} \cdot \mathbf {v} =\mathbf {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}}{\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) product

The scalar product of two vectors in curvilinear coordinates is

u ⋅ 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}}{\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) product

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}u_{j}v_{k}\mathbf {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}}{\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}}\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}}{\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}}\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}}{\displaystyle \mathbf {e} _{p}\times \mathbf {e} _{q}=\varepsilon _{ipq}\mathbf {e} _{i}}

and

b 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}.}{\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}}{\partialq^{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}}}}{\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},}{\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}}{\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 \m athbf {b} _{s}]u^{m}v^{n}\mathbf {b} ^{s}={\mathcal {E}}_{smn}u^{m}v^{n}\mathbf {b} ^{s}}

Tensor operations

Identity map

The identity map I {\displaystyle {\mathsf {I}}}\mathsf{I}defined by I ⋅ v = v {\displaystyle {\mathsf {I}}\cdot \mathbf {v} =\mathbf {v} }\mathsf{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 {\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}}{\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 vector

The action v = S ⋅ u {\displaystyle \mathbf {v} ={\boldsymbol {S}}\cdot \mathbf {u} }\mathbf{v} = \boldsymbol{S}\cdot\mathbf{u}can be expressed in curvilinear coordinates as

v 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}}{\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 tensors

The inner product of two second-order tEnsors U = S ⋅ T {\ displaystyle {\ boldsymbol {U}} = {\ boldsymbol {S}} \ cdot {\ boldsymbol {T}}}\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} }{\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}}{\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}}}{\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]}{\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}}\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}). }{\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}} \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}} ~.} \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}} \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} = [\boldsymbol{F}^{\rm{T}}\cdot\boldsymbol{F}]_{ij}

Пусть J: = det F {\ displaystyle J: = \ det {\ boldsymbol {F}}}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] ~.} \left[\mathbf{b}_1,\mathbf{b}_2,\mathbf{b}_3\right] = \det\boldsymbol{F}\left[\math bf{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} \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})} 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}]}{\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}} 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}}}} \det[g^{ij}] = \cfrac{1}{J^2}

Следовательно, используя тот факт, что [gij] = [gij] - 1 {\ displaystyle [g ^ {ij}] = [g_ {ij}] ^ {- 1}}[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}} \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}}}} \mathbf{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}} \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}} \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^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}} \math bf{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} \mathbf{x}_i(\alpha) = \boldsymbol{\varphi}(\alpha, q^j, q^k) ~,~~ i\ne j \ne 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}}}} \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(\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}} [\boldsymbol{\nabla}f(\mathbf{x})]\cdot\mathbf{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})} 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}} [\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})}f(\mathbf{x}) = \psi^i(\mathbf{x}), то, поскольку qi = ψ i (x) {\ displaystyle q ^ {i} = \ psi ^ {i} (\ mathbf {x})}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}} [\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}} [\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}} [\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 that

c = ∂ 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}}}} \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}} \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 by

b i = ∇ ψ i {\displaystyle \mathbf {b} ^{i}={\boldsymbol {\nabla }}\psi ^{i}} \mathbf{b}^i = \boldsymbol{\nabla}\psi^i

Christoffel symbols of the first kind

The Christoffel symbols of the first kind are defined as

b 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}} \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 that

g 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}}} \begin{align} g_{ij,k} = (\mathbf{b}_i\cdot\mathbf{b}_j)_{,k} = \mathbf{b}_{i,k}\cdot\mathbf{b}_j + \mathbf{b}_i\cdot\mathbf{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{align}

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}]} \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 kind

The Christoffel symbols of the second kind are defined as

Γ i j k = Γ j i k {\displaystyle \Gamma _{ij}^{k}=\Gamma _{ji}^{k}} \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}}\cfrac{\partial \mathbf{b}_i}{\part ial 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}}}} \Gamma_{ij}^k = \cfrac{\partial \m athbf{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}} \cfrac{\partial \mathbf{b}^i}{\partial q^j} = -\Gamma^i_{jk}~\mathbf{b}^k ~;~~ \boldsymbol{\nabla}\math bf{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)} \Gamma^k_{ij} = \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}}} \begin{align} \boldsymbol{\nabla}\mathbf{v} = \left[\cfrac{\partial v^i}{\partial q^k} + \Gamma^i_{lk}~v^l\right]~\mathbf{b}_i\otimes\mathbf{b}^k \\[8pt] = \left[\cfrac{\partial v_i}{\partial q^k} - \Gamma^l_{ki}~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}} \mathbf{v} = v_i~\mathbf{b}^i = \hat{v}_i~\hat{\mathbf{b}}^i

где vi { \ displaystyle v_ {i}}v_{i}- ковариантные компоненты поля, v ^ i {\ displaystyle {\ hat {v}} _ {i}}\hat{v}_i- физические компоненты и (без суммирования )

b ^ i = bigii {\ displaystyle {\ hat {\ mathbf {b}}} ^ {i} = {\ cfrac {\ m athbf {b} ^ {i }} {\ sqrt {g ^ {ii}}}}} \hat{\mathbf{b}}^i = \cfrac{\mathbf{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}} \boldsymbol{\nabla}\boldsymbol{S} = \cfrac{\partial \boldsymbol{S}}{\partial 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}} \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^l_{ki}~S_{lj} - \Gamma^l_{kj}~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 { выровнено}}} \begin{align} \boldsymbol{\nabla}\boldsymbol{S} = \left[\cfrac{\partial S^{ij}}{\partial q^k} + \Gamma^i_{kl}~S^{lj} + \Gamma^j_{kl}~S^{il}\right]~\mathbf{b}_i\otimes\mathbf{b}_j\otimes\mathbf{b}^k \\[8pt] = \left[\cfrac{\partial S^i_{~j}}{\partial q^k} + \Gamma^i_{kl}~S^l_{~j} - \Gamma^l_{kj}~S^i_{~l}\right]~\mathbf{b}_i\otimes\mathbf{b}^j\otimes\mathbf{b}^k \\[8pt] = \left[\cfrac{\partial S_i^{~j}}{\partial q^k} - \Gamma^l_{ik}~S_l^{~j} + \Gamma^j_{kl}~S_i^{~l}\right]~\mathbf{b}^i\otimes\mathbf{b}_j\otimes\mathbf{b}^k \end{align}

Представление физического тензорного поля второго порядка

Физические компоненты тензорного поля второго порядка могут быть получены с использованием нормированного контравариантного базиса, т. е.

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}} \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}}}} \hat{S}_{ij} = S_{ij}~\sqrt{g^{ii}~g^{jj}}

Дивергенция

Векторное поле

Дивергенция векторного поля (v {\ displaystyle \ mathbf {v}}\mathbf {v} ) определяется как

div ⁡ v = ∇ ⋅ v = tr (∇ v) {\ displaystyle \ operatorname {div} ~ \ mathbf {v} = {\ boldsymbol {\ nabla}} \ cdot \ mathbf {v} = {\ text {tr}} ({\ boldsymbol {\ nabla}} \ mathbf {v})} \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}} \boldsymbol{\nabla}\cdot\mathbf{v} = \cfrac{\partial v^i}{\partial q^i} + \Gamma^i_{\ell i}~v^\ell = \left[\cfrac{\partial v_i}{\partial q^j} - \Gamma^\ell_{ji}~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}} \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]} \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}}}\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}}}} 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 }} \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}]}[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}}}[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 have

g = 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}} 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}}}} [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 }} \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}})} \boldsymbol{\nabla} \cdot \mathbf{v} = \cfrac{1}{\sqrt{g}}~\frac{\partial }{\partial q^i}(v^i~\sqrt{g})

Second-order tensor field

The 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})} (\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}}} \begin{align} \boldsymbol{\nabla}\cdot\boldsymbol{S} = \left[\cfrac{\partial S_{ij}}{\partial q^k} - \Gamma^l_{ki}~S_{lj} - \Gamma^l_{kj}~S_{il}\right]~g^{ik}~\mathbf{b}^j \\[8pt] = \left[\cfrac{\partial S^{ij}}{\partial q^i} + \Gamma^i_{il}~S^{lj} + \Gamma^j_{il}~S^{il}\right]~\mathbf{b}_j\\[8pt] = \left[\cfrac{\partial S^i_{~j}}{\partial q^i} + \Gamma^i_{il}~S^l_{~j} - \Gamma^l_{ij}~S^i_{~l}\right]~\mathbf{b}^j \\[8pt] = \left[\cfrac{\partial S_i^{~j}}{\partial q^k} - \Gamma^l_{ik}~S_l^{~j} + \Gamma^j_{kl}~S_i^{~l}\right]~g^{ik}~\mathbf{b}_j \end{align}

лапласиан

Скалярное поле

Лапласиан скалярного поля φ (x ) определяется как

∇ 2 φ: = ∇ ⋅ (∇ φ) {\ displaystyle \ na bla ^ {2} \ varphi: = {\ boldsymbol {\ nabla}} \ cdot ({\ boldsymbol {\ nabla}} \ varphi)} \nabla^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}})} \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}}}} \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)} \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}} \boldsymbol{\nabla}\times\mathbf{v} =\mathcal{E}^{rst} v_{s|r}~ \mathbf{b}_t

where

v s | r = v s, r − Γ s r i v i {\displaystyle v_{s|r}=v_{s,r}-\Gamma _{sr}^{i}~v_{i}} v_{s|r} = v_{s,r} - \Gamma^i_{sr}~v_i
Orthogonal curvilinear coordinates

Assume, 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}}} \mathbf{b}_i\cdot\mathbf{b}_j = \begin{cases} g_{ii} \text{if } i = j \\ 0 \text{if } i \ne 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}}} \mathbf{b}^i\cdot\mathbf{b}^j = \begin{cases} g^{ii} \text{if } i = j \\ 0 \text{if } i \ne j, \end{cases}

where g i i = g i i − 1 {\displaystyle g^{ii}=g_{ii}^{-1}}g^{ii} = g_{ii}^{-1}. As before, b i, b j {\displaystyle \mathbf {b} _{i},\mathbf {b} _{j}}\mathbf{b}_i, \mathbf{b}_jare 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 coordinates

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.,

x = ∑ i = 1 3 x i e i {\displaystyle \mathbf {x} =\sum _{i=1}^{3}x_{i}~\mathbf {e} _{i}} \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) as

d 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}} \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 by

d 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}} \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 quantity

g 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}} 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 that

g 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}} 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, using

b 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}} \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)}{\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}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} \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 |} \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}} {\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}}} \begin{align} \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{align}

Now, since g i j = 0 {\displaystyle g_{ij}=0}{\displaystyleg_{ij}=0}when i ≠ j {\displaystyle i\neq j} i \ne 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}}}} \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 integrals

Likewise, 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}\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|} \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 }} {\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} = \sum_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}}} \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 \part ial 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}}}{\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}}}\begin{vmatrix} \mathbf{e}_1 \mathbf{e}_2 \math bf{e}_3 \\ \\ \sum_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} \\ \\ \sum_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, Laplacian

In orthogonal curvilinear coordinates of 3 dimensions, where

b 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}}}} \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}} \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 basis

g = 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}} 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})} \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}}}} 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)} \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 that

g 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\}} g^{li}~\frac{\partial \varphi}{\partial q^l} = \left\{ g^{11}~\frac{\partial \varphi}{\partial q^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)} \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}}}} \nabla\times\mathbf{v} = \frac{1}{h_1h_2h_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 coordinates

For 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\}} (x_1, x_2, x_3) = \mathbf{x} = \boldsymbol{\varphi}(q^1, q^2, q^3) = \boldsym bol{\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} \}} \{\psi^1(\mathbf{x}), \psi^2(\mathbf{x}), \psi^3(\mathbf{x})\} = (q^1, q^2, q^3) \equiv (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 0 <r <\infty ~, ~~ 0 <\theta <2\pi ~,~~ -\infty <z <\infty

Тогда ковариантный и контравариантный базисные векторы равны

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 { выровнено}}} \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{align}

где er, e θ, ez {\ displaystyle \ mathbf {e} _ {r}, \ mathbf {e} _ {\ theta}, \ mathbf {e} _ {z} }\mathbf{e}_r, \mathbf{e}_\theta, \mathbf{e}_z- единичные векторы в направлениях r, θ, z {\ displaystyle r, \ theta, z}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} g^{ij} = g_{ij} = 0 (i \ne 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} \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}} \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})} (\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}} \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}}} \begin{align} \boldsymbol{\nabla}\mathbf{v} = \cfrac{\partial 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}_\th eta\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}}} \begin{align} \boldsymbol{\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}\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]} \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}}} } \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^2f}{\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^2f}{\partial \theta^2} + \cfrac{\partial^2 f}{\partialz^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}}}{\displaystyle {\begin{aligned}{\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_{\thetaz}\\[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{aligned}}}

Градиент секунды тензорное поле -порядка

Используя приведенные выше определения, мы можем показать, что градиент тензорного поля второго порядка в цилиндрических полярных координатах может быть выражен как

∇ 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}}} \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{\partialS_{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\otim es\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\math bf{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{\partial 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]~\math bf{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\otimes\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}}}{\displaystyle {\begin{aligned}{\boldsymbol {\nabla }}\cdot {\boldsymbol {S}}={\frac {\partial S_{rr}}{\partial r}}~\mathbf {e}_{r}+{\frac {\partial S_{r\theta }}{\partial r}}~\mathbf {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 {\partial S_{\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 also
References
Notes
Further reading
External links
Последняя правка сделана 2021-06-10 13:37:59
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте