Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Введение в тензорный анализ (120

..pdf
Скачиваний:
8
Добавлен:
15.11.2022
Размер:
525.77 Кб
Скачать

Copyright ОАО «ЦКБ «БИБКОМ» & ООО «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 = gTαj ,

где (gij) — матрица, обратная к матрице (gij):

ggαj = δji ,

получим из тензора T поднятием индекса. Если сначала опустить индекс, а затем поднять (или наоборот), то получится первоначально заданный тензор. В самом деле,

gβiSβj = gβigβαTjα = ggβα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 ))x1, . . . , ξ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