Введение в тензорный анализ (120
..pdfCopyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
В отличие от векторных полей отображение ковекторных полей всегда можно определить по формуле
(f (ω))x(ξ) = ωf(x)((f )x ξ), |
|
|||||||||||||
где ξ — вектор из Tx(R2). |
|
|
|
|
|
|
|
|
|
|
|
|
||
Если обозначить ϕ = f ω, можно получить |
|
|||||||||||||
ω2 |
|
|
∂f1 |
|
∂f1 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
ϕ2 |
|||||
∂f2 |
|
∂f2 |
||||||||||||
ω1 |
= |
|
∂x1 |
|
∂x2 |
|
|
ϕ1 |
, |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
∂x |
|
|
|
|
|||
или |
|
|
|
|
∂f α |
|
|
|
|
|
||||
|
|
ωi = |
|
|
|
|
ϕα. |
|
|
|
|
|||
|
|
|
∂x |
i |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Если отображение f обратимо и g = f−1, то последнюю формулу можно переписать в виде
ϕi = |
∂gα |
ωα. |
(1.14) |
|
∂xi |
||||
|
|
|
2. ТЕНЗОРНЫЕ ПОЛЯ
2.1. Определение тензора
Введем определение тензорного поля. Для упрощения обозначений ограничимся случаем поля типа (2, 1), общий случай ничем не отличается. Пусть B(U) — векторное пространство векторных полей в области U, а B (U) — векторное пространство ковекторных полей в той же области. Зададим функцию, линейную по каждому аргументу:
T : B(U) × B(U) × B (U) → F(U). |
(2.1) |
Определим линейность функции T :
T (X + Y, W ; ϕ) = T (X, W ; ϕ) + T (Y, W ; ϕ);
T (fX, W ; ϕ) = fT (X, W ; ϕ),
11
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
где X, Y, W — векторные поля; ϕ — ковекторное поле; f — функция; то же верно для остальных аргументов функции T .
Пусть в области пространства, в которой рассматривается поле T , задана система координат x1, x2, x3. Положим
T |
∂ |
, |
∂ |
, dxk = Tijk . |
|
|
|||
∂xi |
∂xj |
Функции Tijk называют координатами тензорного поля T в данной системе координат.
Пусть z1, z2, z3 — другая система координат в той же области
пространства. Обозначим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
T |
|
∂ |
|
∂ |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
, dzk = Tij0k. |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
∂zi |
|
∂zj |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
Согласно формулам |
(1.9) и (1.10) имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
∂ |
|
∂ |
|
|
|
|
|
|
∂xα ∂ |
∂xβ |
|
|
∂ |
∂zk |
|
|
|||||||||||||||
Tij0k = T |
|
|
|
, |
|
|
|
, dzk = T |
|
|
|
|
, |
|
|
|
|
, |
|
|
dx |
γ |
= |
|||||||||||
∂zi |
∂zj |
∂zi |
∂xα |
∂zj |
∂xβ |
∂xγ |
||||||||||||||||||||||||||||
|
∂xα ∂xβ ∂zk |
∂ |
|
|
∂ |
∂xα ∂xβ ∂zk |
|
|
||||||||||||||||||||||||||
= |
|
|
|
|
|
T |
|
, |
|
|
, dxγ = |
|
|
|
|
|
|
Tαβγ , |
|
|||||||||||||||
∂zi |
∂zj |
∂xγ |
∂xα |
∂xβ |
∂zi |
|
∂zj |
∂xγ |
|
т. е. формулу для преобразования координат тензорного поля типа (2, 1):
Tij0k = |
∂xα ∂xβ ∂zk |
T |
γ . |
(2.2) |
|||||
|
|
|
|
|
|||||
∂zi ∂zj ∂xγ |
|||||||||
|
|
αβ |
|
Таким образом, каждому тензорному полю типа (2, 1) и системе координат x1, x2, x3 ставится в соответствие набор из 32+1 = 27 функций Tijk (x), преобразующихся при переходе к системе коор-
динат z1, z2, z3 по формуле (2.2). Можно показать, что задание такого соответствия определяет тензорное поле в смысле данного выше определения.
Итак, имеем два определения тензорного поля: 1) как линейного отображения (см. формулу (2.1)), 2) с помощью сопоставления каждой системе координат функций Tijk (x) координат тензорного поля, которые преобразуются при замене системы координат по формуле (2.2).
12
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Изложенное выше верно для тензоров любого типа (p, q). Тензорное поле в этом случае называют p раз ковариантным и q раз контравариантным. Число p + q называют рангом тензорного поля.
Упражнение 2.1. Дайте определение тензоров типа (2, 0), (1, 1) и (0, 2). Напишите для них аналоги формулы (2.2).
Тензорные поля возникают при решении многих физических задач.
Пример 2.1. Пусть в точке x проводника приложен вектор E(x) напряженности электрического поля. Тогда в этой же точке возникает вектор j(x) плотности тока. Зависимость j(x) от E(x) линейна (этот факт доказан экспериментально). Если система координат фиксирована, имеем
jl(x) = cls(x)Es(x), l = 1, 2, 3.
Поскольку зависимость j(x) от E(x) линейна, числа cls(x) при изменении системы координат изменяются по формуле (2.2), т. е. являются координатами тензора — тензора теплопроводности. В нашей терминологии имеем тензорное поле типа (1, 1). Отметим, что если среда однородна и изотропна, то матрица (cls) скалярная, т. е. j = cE; в общем случае это не так.
2.2. Операции над тензорными полями
Рассмотрим Tqp(U) — множество тензорных полей типа (p, q). Тензорные поля одного типа можно складывать, умножать на числа, используя формулы, аналогичные формулам (1.3) для векторных полей. Легко проверить, что множество Tqp(U) с этими операциями является векторным пространством (бесконечномерным). Пусть, например, T, S T12(U), λ R. Рассмотрим тензорные поля R = T + S и V = λT . Ясно, что для координат этих тензорных полей справедливы следующие равенства:
Ri |
(x) = T i |
(x) + Si |
(x); V i |
(x) = λT i |
(x). |
jk |
jk |
jk |
jk |
jk |
|
Отметим, что тензорные поля можно умножать не только на числа, но и на функции. Если, например, T T12(U), а f — функция и W = fT , (W T12(U)), то
Wjki (x) = f(x)Tjki (x).
13
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Положим |
M |
|
|
|
T(U) = Tpq(U) |
|
p,q |
(прямая сумма пространств). В множестве T(U) также есть структура линейного пространства, порожденная структурами векторных пространств в каждом множестве из множеств Tqp(U).
Введем в T(U) одну важную операцию — тензорное произве-
q |
s |
|
|
q+s |
|
|
|
дение. Пусть T Tp(U), а S Tr(U). Определим |
|
|
|
|
|||
|
|
тензорное про- |
|||||
изведение полей T и S как тензорное поле R = T |
S Tp+r(U). |
||||||
|
q |
|
, r |
|
s |
|
, |
Объясним эту конструкцию на примере: p = 2, N= 0 |
|
= |
|
= 1 |
|
общий случай будет очевиден. Положим
R(X, Y, Z; ϕ) = T (X, Y ) ∙ S(Z; ϕ).
Можно легко проверить, что построенный объект есть линейная функция своих аргументов и, следовательно, является тензорным полем.
Упражнение 2.2. Докажите следующие свойства тензорного умножения:
1) |
|
T |
|
|
|
q |
|
|
|
( |
1 |
+ T2) S = T1 |
S + T2 S, где |
T1, T2 Tp(U), S |
|||||
s |
|
||||||||
Tr(U); |
|
|
|
q |
, S2 |
|
|||
2) |
T (S1 + S2) = T S1 + T S2, где |
T Tp(U), S1 |
Tsr(U);
3)(λT ) S = T (λS) = λ(T S), где T Tqp(U), S Tsr(U),
λ R;
4) (T S) R = T (S R), где T Tqp(U), S Tsr(U), R Tkm(U).
Свойства 1— 4 показывают, что множество T(U) является кольцом относительно операций сложения и тензорного умножения. Линейные операции и тензорное умножение удовлетворяют обычным аксиомам кольца и векторного пространства, а также обладают свойством
(λT ) S = T (λS) = λ(T S),
где T, S T(U), λ R. Такие алгебраические структуры называют
алгебрами. Таким образом, построили алгебру тензорных полей в
области U, т. е. тензорную алгебру.
14
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Пусть U R3 — область и x1, x2, x3 — произвольная система координат в этой области. Как объяснялось выше (см. разд. 1),
определены векторные поля |
|
∂ |
, |
∂ |
, |
∂ |
и ковекторные поля |
||||||||
∂x1 |
∂x2 |
∂x3 |
|||||||||||||
|
|
|
|
|
|
|
|
|
∂ |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
dx1, dx2, dx3, причем для каждой точки x U векторы |
|
x, |
|||||||||||||
∂x1 |
|||||||||||||||
|
∂ |
|
∂ |
|
|
|
|
|
|
|
|
||||
|
|
x, |
|
|
x образуют |
базис касательного пространства |
|||||||||
∂x2 |
∂x3 |
Tx(R3), линейные функции (dx1)x, (dx2)x, (dx3)x — базис соответствующего сопряженного пространства.
Теорема 2.1. Любое тензорное поле T T12(U) может быть
записано в виде |
|
|
|
|
|
|
|
|
|
T = Tijk dxi dxj |
∂ |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
||||
∂xk |
|
|
|
|
|||||
Д о к а з а т е л ь с т в о. Обе части данной формулы являются |
|||||||||
|
значения на наборе |
|
|
∂ |
|
∂ |
γ |
||
|
|
|
|
|
|
|
|||
|
∂xα |
|
∂xβ , dx |
||||||
тензорными полями. Их |
γ |
|
, |
рав- |
|||||
ны значениям функции |
Tαβ, т. е. координаты тензоров совпадают, |
||||||||
тогда совпадают и сами тензорные поля. |
|
|
|
|
|
|
|
Теорема справедлива и для тензоров произвольного типа (p, q). Упражнение 2.3. Напишите формулы для тензоров второго
ранга, аналогичные формуле, приведенной в теореме 2.1.
Еще одна важная операция в тензорной алгебре — операция
свертки.
Пусть задано тензорное поле T типа (2, 1). Его сверткой по первому верхнему и первому нижнему индексам называют тензорное поле R = C11T (C11 — операция свертки), координаты которого вычисляют по формуле
Ri(x) = Tααi(x).
Докажем, что мы действительно получили координаты тензорного поля R типа (1, 0). С помощью формулы (2.2) проверим справедливость формулы
Ri0(x) = ∂xβ Rβ(x). ∂zi
15
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Используя формулу (1.8), получим |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
R0 |
(x) = T |
0α |
(x) = |
∂xs |
|
∂xt ∂zα |
T |
r |
(x) = δ |
s |
∂xt |
T |
r |
(x) = |
|||||||||
αi |
|
|
α |
|
|
i |
|
|
r |
st |
r |
|
i |
st |
|||||||||
i |
|
|
|
∂z |
|
∂z |
|
∂x |
|
|
∂z |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
∂xβ |
|
|
|
|
|
|
||||||||
|
|
|
= |
∂xt |
|
s |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
Tst(x) = |
|
|
|
|
Rβ(x), |
|
|
|
|
|
|
|||||||
|
|
|
∂zi |
∂zi |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
что и требовалось доказать.
Можно брать свертку тензорного поля T типа (2, 1) и по другой паре индексов, т. е. можно построить тензор R0 = C21T . В этом случае проводят свертку по первому верхнему и второму нижнему индексам.
Упражнение 2.4. Найдите координаты тензора R0.
Пусть задано тензорное поле T типа (p, q). Его сверткой по первому верхнему и первому нижнему индексам будет тензорное поле R = C11T типа (p − 1, q − 1).
Пример 2.2. Пусть ω = ωidxi — ковекторное поле, X = ξjXj —
векторное поле. Возьмем |
тензор |
ω |
X и применим операцию |
||
1 |
|
|
|||
свертки. Получим тензор |
Ci1 |
(ω X) |
типа (0, 0), т. е. функцию, |
||
и эта функция равна ωi ξ . |
С учетом формулы (1.6) получим |
ω(X) = C11(ω X).
Геометрическое определение тензора второго ранга приведено в работе [2].
2.3. Тензоры в пространствах с метрикой
Пусть в каждом касательном пространстве Tx(R3) задано скалярное произведение. Положим
gij(x) = ((Xi)x, ( Xj)x).
Тогда функции gij(x) представляют собой координаты так называемого метрического тензора. Его можно записать в виде
g = gijdxi dxj.
Очевидно, что матрица (gij) симметрическая и квадратичная
форма
ds2 = gijdxidxj
положительно определена.
16
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Пример 2.3. Евклидова метрика в пространстве R3 есть ds2 = dx2 + dy2 + dz2.
Если в пространстве задана метрика (метрический тензор), то определены операции поднимания и опускания индексов.
Пусть, например, T — тензорное поле типа (1, 1). Возьмем тензорное произведение g T , а затем применим операцию свертки. В результате получим тензор типа (2, 0): S = C11(g T ). Его координаты вычисляют по формуле
Sij = gαiTjα.
Говорят, что тензор S получен из тензора T опусканием индекса. Тензор R с координатами
Rij = giαTαj ,
где (gij) — матрица, обратная к матрице (gij):
giαgαj = δji ,
получим из тензора T поднятием индекса. Если сначала опустить индекс, а затем поднять (или наоборот), то получится первоначально заданный тензор. В самом деле,
gβiSβj = gβigβαTjα = giβgβαTjα = δiαTjα = Tji.
Пример 2.4. Пусть задана евклидова метрика, т. е. gij = δij. Тогда
Sij = gαiTjα = δαiTjα = Tji.
Таким образом, в случае евклидовой метрики верхние и нижние индексы не различаются.
2.4.Действие отображений на тензорные поля
Вразд. 1 было показано, как действуют отображения на векторные и ковекторные поля, т. е. на тензорные поля типов (0, 1) и (1, 0) соответственно. Рассмотрим общий случай. Пусть f : U → V , где U Rn, V Rm — гладкое отображение. Поскольку невозможно определить отображение векторных полей f : T1(U) → T1(V ),
17
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
это нельзя сделать и для полей типа Tq(U) при q > 1. В то же время отображение f : Tp(V ) → Tp(U) определить можно. Положим
|
(f (T ))x(ξ1, . . . , ξp) = Tf(x)((f )x ξ1, . . . , (f )x ξp), |
(2.3) |
|||||
где ξ1, . . . , |
ξp Tx(Rn). |
|
|
|
|
||
Пример 2.5. Рассмотрим в пространстве R3 евклидову мет- |
|||||||
рику, приведенную в |
примере |
2.3. Пусть задана |
поверхность |
||||
r = r(u, v), (u, v) U R2, где |
r = r(x, y, z), x = x(u, v), |
y = |
|||||
= y(u, v), |
z = z(u, v). |
Имеем, |
таким |
образом, |
отображение |
||
|
→ R |
|
|
|
|
3 |
|
f : U |
|
3 , которое индуцирует тензор |
h = f g типа (2, 0) в |
области U, где g — евклидова метрика в пространстве R . Обозначив x = x1, y = x2, z = x3, u = z1, v = z2, с помощью формулы (2.3) получим
|
∂xα |
∂xβ |
|
∂xα ∂xβ |
|
∂xα ∂xβ |
||||||||
hij = g |
|
Xα, |
|
Xβ |
= |
|
|
|
g Xα, Xβ = |
|
|
|
δαβ. |
|
∂zi |
∂zj |
∂zi |
∂zj |
∂zi |
∂zj |
|||||||||
В классических обозначениях имеем |
|
|
|
|
|
|||||||||
h11 = (ru, ru); h12 = (ru, rv), h22 = (rv, rv), |
||||||||||||||
где ru = (xu, yu, |
zu ), rv = (xv, yv, zv). Мы получили первую |
|||||||||||||
квадратичную форму поверхности: |
|
|
|
|
|
|||||||||
|
|
|
|
ds2 = hijdxidxj. |
|
|
|
|
|
|||||
Пример 2.6. Вычислим первую квадратичную форму сферы |
||||||||||||||
единичного радиуса. Зададим |
|
|
|
|
|
|||||||||
|
x = cos u cos v, |
y = cos u sin v, z = sin u. |
||||||||||||
Значит, |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
ru = (− sin u cos v, − sin u sin v; |
cos u) ; |
||||||||||||
|
|
rv = (− sin v cos u, cos v cos u; |
0) , |
|
|
|
|
откуда
ds2 = du2 + cos2u dv2.
Упражнение 2.5. Пусть поверхность задана как график функции
z = f (x, y) .
18
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Докажите, что
ds2 = 1 + fx2 dx2 + 2fxfydxdy + 1 + fy2 dy2.
Пусть теперь отображение f : U → V — диффеоморфизм. Тогда отображение f˜: Tqp(U) → Tqp(V ) можно построить. Покажем, как это можно выполнить. Следует обратить внимание на то, что отображения f и f действуют в разные стороны. Поэтому векторное поле сносится с помощью отображения f , а ковекторное — с помощью отображения f−1 . Если (x1, x2), (z1, z2) — системы координат соответственно в областях U и V и отображение f задается формулами
z1 = z1 x1, x2 ; |
z2 = z2 x1, x2 , |
то координаты тензора f˜T типа (1,1) вычисляются по формуле, аналогичной формулам (1.13), (1.14):
i |
∂zi ∂xβ |
|
|||
f˜T j = |
|
|
|
Tβα. |
(2.4) |
∂xα |
∂zj |
Такие же формулы имеют место и для тензорного поля произвольного типа (p, q).
3. ПРОИЗВОДНАЯ ЛИ
Рассмотрим дифференцирование тензорных полей. Назовем отображение D: T(U) → T(U) дифференцированием, если оно аддитивно, сохраняет тип тензора и обладает свойствoм
D (T S) = DT S + T DS |
(3.1) |
для всех T, S T(U).
Дифференцирование полностью определяется своим действием на функции, векторные и ковекторные поля. Это следует из формулы (3.1) и теоремы 2.1.
3.1. Фазовый поток
Для упрощения обозначений считаем, что n = 2. Пусть в области U задано векторное поле X = ξiXi. С этим полем связана
19
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
система дифференциальных уравнений |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
dx1 |
= ξ1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dx2 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
dt |
x1(t), x2(t) ; |
|
|||||||
|
|
|
|
|
|
|
|
|
2 |
|
1 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
= ξ |
x (t), x |
|
(t) . |
|
||||
|
|
|
|
|
|
|
|
|
|
||||||||
|
i |
|
i |
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 1, 2, называют фазовыми кривыми векторного |
|||||||||||
Линии x = x (t), i |
|
|
|
|
|
|
|
|
|
|
|
||||||
поля X. Обозначим |
|
|
|
|
|
|
|
|
|
|
|
||||||
фазовые |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
F i |
x1 |
, x2 |
|
= xi = xi |
t, x1, x2 |
, i = 1, 2, |
(3.2) |
|||||||
|
|
|
t |
|
0 |
0 |
|
|
|
|
0 |
|
0 |
|
|
|
|
|
кривые поля X с начальными условиями |
|
|||||||||||||||
|
|
|
|
|
|
|
xi|t=0 = x0i , i = 1, 2. |
|
|
||||||||
Формула (3.2) задает отображение |
|
|
|
|
|
|
|||||||||||
|
|
Ft : x01, x02 → (x1 t, x01, x02 , x2 t, x01, x02) |
|
||||||||||||||
|
|
|
зависящее от параметра t (сдвиг за время t вдоль фазо- |
||||||||||||||
области |
U, |
|
|
|
|
|
|
|
|
|
|
|
|
|
вой кривой). Из теории обыкновенных дифференциальных уравнений известно, что отображение Ft, определенное при малых t в окрестности данной точки (x10, x20), является диффеоморфизмом. При малых t явный вид отображения Ft таков:
xi t, x10, x20 = xi0 + tξi x10, x20 + o(t), i = 1, 2.
С той же точностью o(t) матрица Якоби отображения Ft имеет вид
∂xi = δi + t ∂ ξi + o(t). ∂xj0 j ∂xj0
Матрица Якоби обратного отображения выглядит образом:
∂xi0 = δi − t ∂ ξi + o(t). ∂xj j ∂xj
(3.3)
следующим
(3.4)
3.2. Производная Ли тензорного поля
Приведем основное определение производной Ли тензорного поля. Ограничимся частным случаем тензоров типа (2, 1).
20