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

Лосик М.В. - Лекции по векторному и тензорному анализу

.pdf
Скачиваний:
152
Добавлен:
09.06.2015
Размер:
420.3 Кб
Скачать

51

3.10.Дифференциальные операции второго порядка

Мы будем называть дифференциальной операцией первого поряд-

êà взятие градиента скалярного поля или взятие либо дивергенции, либо ротора векторного поля.

Дифференциальной операцией второго порядка называется

последовательное применение двух дифференциальных операций первого порядка.

По определению у нас есть следующие дифференциальные операции

второго порядка: div grad ϕ, rot grad ϕ, grad div˜a, div rot~a è rot rot~a.

~

Из них две операции тождественно равны нулю: rot grad ϕ = 0 согласно теореме 22 и div rot~a = 0 согласно теореме 25.

Дифференциальная операция второго порядка div grad ϕ называет-

ся оператором Лапласа (или лапласианом) и обозначается ϕ. Выразим оператор Лапласа в декартовых координатах и через оператор Гамиль-

òîíà r. Применяя формулы (2.4) и (3.9), мы получаем

 

2ϕ

 

2ϕ

 

2ϕ

ϕ =

 

=

 

+

 

.

∂x2

∂y2

∂z2

Применяя формулы (3.14) и (3.17), мы получаем

ϕ = (r, rϕ) = (r, r)ϕ = r2ϕ.

Выразим оставшиеся дифференциальные операции второго порядка че- рез оператор r, используя формулы (3.17) и (3.21),

grad div˜a = r(r,~a),

rot rot~a = [r, [r,~a]].

3.11.Выражение дивергенции и ротора векторного поля в ортогональных криволинейных координатах

Пусть задана ортогональная криволинейная система координат с коор-

динатами u, v, w. Напомним (см. раздел 1.6), что в этом случае в каждой точке области задания координат имеется ортонормированный ба-

çèñ ~eu, ~ev, ~ew и коэффициенты Ламе Hu, Hv, Hw. Мы будем считать базис

~eu, ~ev, ~ew

 

 

íèÿ

положительным, тогда по определению векторного произведе-

 

[~eu, ~ev] = ~ew, [~ev, ~ew] = ~eu, [~ew, ~eu] = ~ev.

(3.28)

52

Заметим, что, взяв одно из этих равенств за начальное, мы получим остальные два циклическими перестановками индексов u, v, w.

Пусть векторное поле ~a задано в этой системе координат, т.е.

~a = au~eu + av~ev + aw~ew.

Используя теоремы 19 и 22, мы получаем

div ~a = (grad au, ~eu) + au div ~eu + (grad av, ~ev) + av div ~ev+

 

 

 

 

 

 

(grad aw, ~ew) + aw div ~ew,

(3.29)

rot~a = [grad au, ~eu] + aurot~eu + [grad av, ~ev] + avrot~ev+

 

 

 

 

 

 

[grad aw, ~ew] + awrot~ew.

(3.30)

Так как слагаемые, содержащие градиенты компонент

au, av, aw

 

найти с помощью формулы (2.5), нам нужно найти

 

 

 

 

 

 

 

 

 

 

дивергенции и роторы,можно

базисных векторных полей

~eu, ~ev, ~ew. Начнем с отыскания роторов этих

полей.

 

 

Возьмем, например, ~eu. Согласно формуле (2.5), мы имеем

~eu =

Hugrad u. Тогда, применяя теорему 22, мы получаем

 

 

rot~eu = rot (Hugrad u) = [grad Hu, grad u]+

 

 

 

 

 

 

 

 

Hurot grad u = [grad Hu, grad u].

Применяя снова формулу (2.5), мы окончательно получаем

 

 

 

rot~eu =

1

[grad Hu, ~eu].

 

 

(3.31)

 

 

 

 

 

 

Hu

 

 

Аналогично мы получаем следующие формулы:

 

 

 

 

1

 

 

1

 

 

 

 

rot~ev =

 

[grad Hv, ~ev],

rot~ew =

 

[grad Hw, ~ew].

(3.32)

Hv

Hw

Согласно формуле (2.5), последней из формул (3.28) и формулам (3.31) и (3.32), мы получаем

div ~eu = div [~ev, ~ew] = ~ewrot ev − ~evrot~ew =

1

(grad Hv~ev~ew)−

Hv

 

1

(grad Hv~ew~ev) =

 

 

1

 

∂Hv

+

1

 

∂Hw

.

 

 

 

 

 

 

 

 

 

 

Hw

 

 

 

HvHu ∂u HwHu ∂u

Легко видеть, что полученную формулу можно записать в виде

div ~eu =

1

 

∂(HvHw)

.

 

 

 

 

HuHvHw

 

 

 

 

 

 

 

 

 

∂u

 

 

 

 

53

Это равенство и формула (2.5) позволяет получить следующее равенство:

(~eu, grad au) + au div ~eu =

1 ∂au

+

 

au

 

∂(HvHw)

=

 

 

 

 

 

 

 

 

Hu ∂u

HuHvHw

 

∂u

 

 

 

 

 

 

 

 

1

 

 

∂(auHvHw)

.

(3.33)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

HuHvHw

∂u

 

Переставляя в последней формуле циклически индексы

 

u, v è w, ìû

получаем следующие равенства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(~ev, grad av) + av div ~ev =

1

 

 

∂(avHuHw)

,

(3.34)

HuHvHw

 

 

 

 

 

 

 

 

∂v

 

 

 

 

 

(~ew, grad aw) + aw div ~ew =

1

 

 

 

∂(awHuHv)

.

(3.35)

 

 

 

 

 

 

 

HuHvHw

∂w

 

 

 

 

 

Заменяя с помощью формул (3.33), (3.34) и (3.35) пары слагаемых в формуле (3.29), мы, наконец, получаем выражение дивергенции векторного поля в ортогональных криволинейных координатах

div ~a = HuHvHw

u∂u

+

∂v

+

w∂wu

v

 

. (3.36)

1

 

∂(a HvHw)

 

∂(avHuHw)

 

∂(a H

H

)

 

Теперь рассмотрим следующую сумму:

[grad au, ~eu] + aurot~eu.

Применяя формулу (3.31), мы получаем

[grad au, ~eu] + aurot~eu = [grad au, ~eu] +

 

au

[grad Hu, ~eu] =

 

 

 

 

 

Hu

[grad au +

au

~eu] =

1

[grad (auHu), ~eu].

 

 

 

 

Hu

 

 

Hu

Циклически переставляя индексы u, v è w, получаем

1

 

 

 

 

 

[grad av, ~ev] + avrot~ev =

 

[grad (avHu), ~ev],

Hv

(3.37)

(3.38)

 

 

[grad aw, ~ew] + awrot~ew =

1

[grad (awHw), ~ew].

(3.39)

 

 

Hw

Заменяя с помощью формул (3.37), (3.38) и (3.39) пары слагаемых в

формуле (3.30), мы получаем

 

 

 

 

 

 

rot~a =

1

[grad (auHu), ~eu]+

 

 

 

 

 

 

 

 

 

 

 

 

 

Hu

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

[grad (avHu), ~ev] +

 

[grad (awHw), ~ew].

 

 

 

Hv

Hw

 

 

 

 

 

 

 

 

 

 

 

 

 

 

54

Учитывая (3.28) и ортонормированность базиса

~eu, ~ev, ~ew, мы можем вы-

числить первую координату

 

 

 

 

 

 

образом:

 

 

(rot~a)u векторного поля rot~a следующим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(rot~a)u = (rot~a, ~eu) =

1

(grad (avHu), ~ev, ~eu)+

 

 

 

 

 

 

 

 

 

 

 

Hv

 

1

 

 

 

 

1

 

 

 

 

(3.40)

 

 

 

(grad (awHw), ~ew, ~eu) =

 

(~ew, grad (avHv))+

 

Hw

Hv

 

Hw (~ev, grad (awHw)) =HvHw

 

∂v

∂w

.

1

 

 

 

1

 

 

 

∂(awHw)

 

∂(avHv)

 

Циклически переставляя индексы u, v и w, получаем остальные координаты векторного поля rot~a

(rot~a)v = HwHu

 

∂w

∂u

 

,

(3.41)

1

 

∂(auHu)

 

∂(awHw)

 

 

(rot~a)w = HuHv

 

∂u

∂v

 

,

(3.42)

1

 

 

∂(avHv)

 

∂(auHu)

 

 

 

Используя формулы (3.40), (3.41) и (3.42), окончательно получаем

следующую формулу для ротора в криволинейных ортогональных коор-

динатах:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

∂(awHw)

 

 

∂(avHv)

 

 

 

 

 

 

 

 

rot~a =

 

 

 

 

 

 

 

 

 

 

~eu+

 

 

 

 

 

HvHw

 

∂v

∂w

 

 

 

 

 

 

 

 

 

1

 

∂(auHu)

∂(awHw)

~ev+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

HwHu

 

∂w

 

 

∂u

~ew. (3.43)

 

 

 

 

 

 

 

 

 

 

1

 

 

 

∂(avHv)

∂(auHu)

 

 

 

 

 

 

 

 

 

 

HuHv

 

∂u

 

∂v

 

В частности, напишем эти формулы для цилиндрических и сфери- ческих координат. Напомним (см. раздел 1.6), что для цилиндрических

координат ρ, ϕ, z, Hρ = 1, Hϕ = ρ è Hz = 1, а для сферических координат

r, θ, ϕ Hr = 1, Hθ = r è Hϕ = r sin θ. Поэтому формулы (3.36) и (3.43) в этом случае переписываются следующим образом:

 

div ~a = ρ

 

∂ρρ

+

∂ϕ

+

 

∂z

 

,

 

 

 

1

 

 

 

∂(ρa

)

 

 

∂(aϕ)

 

∂(ρaz)

 

 

1

r2 sin θa

)

 

∂(r sin θa

 

∂(ra

)

div ~a =

 

 

 

(

 

 

 

r

 

+

 

 

θ)

+

 

ϕ

 

r2 sin θ

 

 

 

∂r

 

 

 

 

∂θ

 

 

 

∂ϕ

 

rot~a =

1

 

∂(az)

∂(ρaϕ)

~eρ+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

∂ϕ

∂z

∂ρz

 

 

~eϕ +

ρ

∂ρϕ

 

 

 

 

 

 

 

 

∂zρ

)

 

)

 

 

 

 

 

 

 

∂(a

 

∂(a

)

 

1

 

∂(ρa

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

(aρ) ~ez, ∂ϕ

rot~a = r sin θ

 

∂θ ϕ

∂ϕθ

 

 

~er+

 

 

 

 

 

 

 

 

 

 

55

 

)

 

 

 

 

 

 

 

 

 

 

 

1

 

 

∂(sin θa )

 

∂(ρa

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~eθ + r

 

 

 

 

 

 

 

 

 

r sin θ

∂ϕr

)

∂r

ϕ

 

∂rθ

)

∂θr

 

~eϕ.

 

1

 

 

∂(a

 

∂(r sin

θa

)

1

 

∂(ra

 

∂(a

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Глава 4.

Элементы тензорного анализа

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

4.1.Определение тензора и примеры тензоров

Пусть L n-мерное действительное линейное пространство. Мы далее

воспользуемся обозначениями раздела 3.1. Рассмотрим два базиса (ei)

(i = 1, . . . , n) è (ei0 ) (i0 = 10, . . . , n0) пространства L. Напишем разложения первого базиса по второму и второго базиса по первому:

n0

n

 

ei =

Aii0 ei0 , ei0 = Aii0 ei.

(4.1)

iX0 0

X

 

=1

i=1

 

В дальнейшем мы будем использовать следующее соглашение: если в некотором выражении имеются два одинаковых индекса, один из которых стоит вверху, а другой внизу, то по этим индексам производится суммирование по всем значениям этого индекса, даже если это суммирование и не указано специально. Поэтому, чтобы не было путаницы, в выражении не может быть более двух одинаковых индексов, а если их точно два, то они не могут стоять одновременно вверху или внизу. Поэтому равенства (4.1) могут быть записаны следующим образом

ei = Aii0 ei0 , ei0 = Aii0 ei.

(4.2)

57

скаля-

 

 

 

 

 

 

 

 

58

Напомним, что матрицы

i0

è

i

взаимно обратны, т.е.

i0

i

i0

Обратим внимание на(Aòî,)÷òî(Aобъекты,)

которые были рассмотреныA A = δ .

 

i

 

i0

 

 

i

j0

j0

ранее, а именно, векторы, линейные операторы и p-линейные формы ре-

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

ра базиса. Посмотрим, как преобразуются эти числа при замене базиса.

Рассмотрим несколько примеров.

 

(xi) относительно

1. Вектор x L задается набором своих компонентi

базиса

 

, которые определяются равенством

x = x ei

. Используя (4.2),

найдем его компоненты

 

 

(ei)

 

 

 

 

Отсюда следует, что

(xi0 ) относительно базиса (ei0 ): x = xiei = xiAii0 ei0 .

 

 

 

 

xi0 = Aii0 xi.

 

(4.3)

2. Согласно разделу 3.1, iлинейный оператор f : L → L определяется

заданием своей матрицы (aj) относительно базиса (ei), которая опреде-

ляется равенством

 

i

 

 

этой матрицы, согласно (3.5), преобразуются.Кроме того, приследующимзамене базисаобразом:элементы

 

 

 

f(ej) = ajei

 

 

 

 

 

 

aji00 = Aii0 Ajj0 aji .

 

(4.4)

3. Согласно теореме 20, p-линейная форма f(x1, . . . , xp) однозначно

определяется заданием набора чисел fi1...ip = f(ei1 , . . . , eip ), зависящих от

выбора базиса (ei). Выясним, как эти числа преобразуются при замене базиса.

fi

 

...i

 

 

= f(ei0 , . . . , ei0

) = f(A

i1

ei1 , . . . , A

ip

ei0

) =

 

 

0

p

0

i0

i0

 

 

1

 

1

p

 

 

 

 

p

 

 

 

 

 

 

 

 

1

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

i1

ip

f(ei1

i1

ip

fi1...ip ,

 

 

 

 

 

 

 

Ai0

. . . Aip0

, . . . , eip ) = Ai0

. . . Aip0

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

где мы использовали последовательно линейность формы

f ïî âñåì åå

аргументам

x1

, . . . , xp

. Таким образом, набор чисел

fi1

...ip

 

 

ïî

 

 

 

 

 

 

 

 

 

 

преобразуется

 

следующему закону:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

...ip0

= A

i1

. . . A

ip

f

i1...ip

.

 

 

 

(4.5)

 

 

 

 

i10

ip0

 

 

 

 

 

 

 

i10

 

 

 

 

 

 

 

 

Каждое действительное число c можно рассматривать как постоянную функцию на множестве базисов пространства L c законом преобразования: c0 = c. Мы будем называть такую постоянную функцию

ðîì.

Заметим, что все рассмотренные выше объекты являются функциями

на множестве базисов пространства L, удовлетворяющими некоторому закону преобразования при замене базиса. Теперь мы дадим общее определение, которому удовлетворяют все рассмотренные выше объекты.

59

Определение 23. Говорят, что задан тензор T типа (p, q) (p, q ≥ 0) на n-мерном линейном пространстве L, если каждому базису (ei) ïðî-

странства L поставлен в соответствие набор np+q чисел T i1...ip

j1...jq , заданных для любых i1, . . . , ip, j1, . . . , jq = 1, . . . , n, которые при замене базиса

(ei) базисом (ei0 ) преобразуются следующим образом:

i10 ...ip0

i10

ip0

j1

jq

i1...ip

 

Tj10 ...jq0

= Ai1

· · · Aip Aj10

· · · Ajq0

Tj1...jq .

(4.6)

Индексированные числа T i1...ip T j1...jq называются компонентами тензора

относительно базиса (ei).

Заметим, что уравнения преобразования компонент тензора зависят от типа тензора.

Сравним это общее определение с приведенными выше примерами. Все уравнения преобразования (4.3), (4.4), (4.5) и уравнение преобразования скаляра являются частными случаями общего преобразования

компонент тензора из определения 23 для различных выборов p è q. Пример 1 показывает, что вектор является тензором типа (1, 0). Пример 2 показывает, что линейный оператор можно считать тензором типа (1, 1). Пример 3 показывает, что p-линейная форма может рассматриваться как

тензор типа (0, p). Наконец, скаляр является тензором типа (0, 0).

Введем теперь несколько операций, производимых над тензорами.

1. Сложение тензоров.

i1...ip

 

i1...ip

 

Пусть P = (Pj1...jq ) è Q = (Qj1...jq ) два тензора одного типа (p, q).

Определим их сумму P + Q равенством

 

i1...ip

i1...ip

i1...ip

(P + Q)j1...jq

= Pj1...jq

+ Qj1...jq ,

 

i1...ip

которое определяет компоненты (P +Q)j1...jq величины P +Q относитель-

но базиса (ei). Покажем, что P + Q является тензором типа (p, q), т.е. эти компоненты удовлетворяют уравнениям (4.6). Действительно, используя

уравнения преобразования компонент тензоров P è Q, получаем

i0 ...i0

i0

...i0

i0 ...i0

 

 

 

 

 

 

 

(P + Q) 1

p

= P 1 p + Q 1

p =

 

 

 

 

 

 

 

j10 ...jq0

j10 ...jq0

j10 ...jq0

 

 

 

 

 

 

 

i10

 

ip0

j1

jq

 

i1...ip

 

i10

ip0 j1

 

jq

i1...ip

 

Ai1

· · · Aip Aj10 · · ·

Ajq0

Pj1...jq

+ Ai1 · · · Aip Aj10

· · · Ajq0

Qj1...jq

=

 

 

i10

 

ip0

j1

 

jq

i1...ip

i1

...ip

=

 

 

 

Ai1

· · · Aip Aj10

· · · Ajq0

Pj1...jq

+ Qj1...jq

 

 

 

 

 

 

 

 

 

i10

ip0

j1

 

jq

i1...ip

 

 

 

 

 

 

 

 

Ai1 · · · Aip Aj10

· · · Ajq0 (P + Q)j1...jq .

Отметим, что операция сложения тензоров содержит в качестве частных случаев операции сложения скаляров, векторов, линейных операто-

ðîâ è p-линейных форм.

относитель-
(P · Q)i1...ip+r j1...jq+s
...ip+r ,
...jq+s
= P i1...ip Qip+1 j1...jq jq+1

 

 

 

 

 

 

60

2. Умножение тензоров.

i1

...ir

 

Пусть P = (P

i1

...ip

) тензор типа (p, q) è Q = (Q

) тензор типа

j1...jq

j1...js

 

 

 

(r, s). Определим их произведение (тензорное произведение) P ·Q равенством

которое определяет компоненты (P ·Q)i1...ipr P ·Q j1...jq+s величины

но базиса (ei). Покажем, что P · Q является тензором типа (p + r, q + s), т.е. эти компоненты удовлетворяют уравнениям (4.6). Действительно, ис-

пользуя уравнения преобразования компонент тензоров P è Q, получаем

i0 ...i0

(P · Q) 1 p+r

j10 ...jq0 +s

i10

ip0

Ai1

· · · Aip

= P

 

i0 ...i0

 

i0

...i0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 p Q p+1

 

p+r =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j10 ...jq0

 

jq0 +1...jq0

+s

 

 

 

 

 

 

 

 

 

 

 

 

 

j1

 

 

j

q

i ...i

p

i0

 

 

i0

j

q+1

 

j

q+s

i

p+1

...i

p+r

 

 

 

1

 

p+1

· · ·

 

p+r

 

 

 

 

 

Aj10

· · · Ajq0

Pj1...jq Aip+1

Aip+r Ajq0 +1

· · · Ajq0

+s Qjq+1...jq+s =

 

i10

 

 

i0

 

j1

 

jq+s

i1...ip

 

ip+1...ip+r

 

 

 

 

 

 

 

 

p+r

 

 

 

 

 

 

 

 

Ai1

· · · Aip+r Aj10

· · · Ajq0 +s Pj1...jq

Qjq+1...jq+s

=

 

 

 

 

 

 

 

 

 

 

 

 

 

i10

 

i0

 

 

j1

 

 

jq+s

 

 

 

i1...ip+r

 

 

 

 

 

 

 

 

 

 

p+r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ai1

· · · Aip+r Aj10

· · · Ajq0 +s

(P · Q)j1...jq+s .

Отметим, что операция умножения тензоров содержит в качестве частных случаев операции умножения чисел, умножения числа на век-

тор, умножения числа на линейный оператор и умножения числа на p- форму.

3. Свертывание (свертка) тензора.

Пусть P = (P i1...ip ) тензор типа (p, q) è p, q > 0. Выберем целые

j1...jq

числа r, s такие, что 1 ≤ r ≤ p è 1 ≤ s ≤ q. Определим свертку SsrP тензора P равенством

(SrP )i1...ip−1

= P i1...ir−1kir...ip−1

 

s j1...jq−1

 

j1...js−1kjs...jq−1

 

которое определяет компоненты (SrP )i1...ip−1

SrP относитель-

 

s

j1...jq−1 величины

s

но базиса (ei). Покажем, что SsrP является тензором типа (p − 1, q − 1). Для простоты обозначений мы рассмотрим частный случай когда s = p

è r = q. Используя уравнения преобразования компонент тензора P , получаем

p

i0

...i0

 

i0 ...i0

 

k0

1

i0

 

k0

j1

jq

 

l

i1...ip

1k

 

1

p

 

1

1 p

 

1

 

i0

p

 

1

 

Aj10

 

1

 

 

(Sq P )j10 ...jq0

−1

= Pj10 ...jq0

−1k0

= Ai1

· · · Aip−1 Ak

· · · Ajq0

−1 Ak0 Pj1...jq−1l .

Учитывая равенство Al

Ak0

= δl

 

 

 

 

 

 

 

l, мы получаем

 

 

 

 

 

 

 

 

 

k0

 

k

 

k и суммируя по

 

 

 

 

 

 

 

 

 

p

i0

...i0

1

i0

 

i0

1

j

jq

1

i1...ip

1k

 

 

 

 

 

 

 

 

 

1

p

 

p

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

· · · Ajq0

 

 

 

 

 

 

 

 

 

 

 

 

 

(Sq

P )j10 ...jq0

−1

= Ai1

· · · Aip−1 Aj10

−1 Pj1...jq−1k

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i0

 

i0

 

1

j

jq

1

p

i1...ip

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

· · · Ajq0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ai1

· · · Aip−1 Aj10

−1

(Sq

P )j1...jq−1

,