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

Мат Анал / Kratn_2

.pdf
Скачиваний:
21
Добавлен:
16.10.2021
Размер:
495.14 Кб
Скачать

Теорема 4.3 [1] (теорема Стокса). Пусть Ω – ограниченная область в R3, внутри которой лежит гладкая ориентированная поверхность Σ+ с кусочно-гладким краем l+ и направление на l+ согласовано с

~

ориентацией Σ+. Пусть также функция f(P ) непрерывно дифференцируема в Ω. Тогда

IZZ

~ ~

~

(f, dl) =

(rot f, d~σ).

l+

Σ+

Теорема Стокса означает, что циркуляция дифференцируемого век-

~

торного поля f по произвольной кусочно-гладкой замкнутой кривой l+ рав-

~

на потоку вектора rot f через поверхность Σ+, ограниченную этой кривой l+.

Отметим, что теорема Грина [1] является частным случаем теоремы Стокса и соответствует тому случаю, когда Σ является частью плоскости

0xy, ~n+ = ~k и f~ = fx~i + fy~j. При этом rot f~ =

∂xy

∂yx

~k.

 

∂f

 

∂f

 

 

 

~

~

~ ~

Пример 4.12. Вычислить циркуляцию вектора f = zi + xj + yk, по

окружности x2 +y2 +z2 = R2, x+y +z = R, в положительном направлении

(обход против часовой стрелки).

 

 

 

 

~

~

~

~

~

~

~

~

Так как f

= zi

+ xj

+ yk, то rot f

= i

+ j

+ k. За поверхность Σ,

ограниченную данной окружностью, примем сам круг, образованный сечением шара x2 + y2 + z2 = Rz плоскостью x + y + z = R. Центр этого

r

круга 00

R

;

R

;

R

, его радиус R1

= R

 

2

. Единичный вектор нормали

3

3

 

 

1

 

 

3

 

 

 

3

3

 

 

 

 

 

 

 

 

 

 

~

~

 

~

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~n =

 

(i + j + k ). Так как (rot f, ~n) =

 

=

 

 

3, то находим

 

 

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2πR2

I (f~, dl~) =ZZ (rot f~, d~σ) =ZZ (rot f~, ~n )dσ = 3ZZ

dσ = 3πR12 =

 

.

3

l+

 

Σ+

 

 

 

 

 

Σ+

 

 

 

 

 

 

 

Σ+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

~

~

 

 

~

 

 

 

Пример 4.13. Найти циркуляцию вектора2f =2yi−22zj +2xk вдоль эл-

липса, образованного сечением гиперболоида 2x −y +z

= R

плоскостью

y = x, в положительном направлении.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

~

~

~

~

 

~

~

~

 

 

 

 

 

 

 

Так как f = yi

− 2zj + xk, то rot f = 2i

− j − k. За поверхность Σ,

ограниченную кривой l, примем часть секущей плоскости, лежащей внутри эллипса. Единичный вектор нормали, направленный в нужную сторону,

1

 

~

~

3

 

 

 

 

 

имеет вид ~n =

 

 

(i

− j ). Поэтому (rot~a, ~n ) =

 

. Значит,

 

2

2

dσ =

I (f~, dl~) =

ZZ (rot f~, d~σ ) =

ZZ (rot f~, ~n )dσ = 2

ZZ

 

 

 

 

 

 

 

3

 

 

 

l+

Σ+

Σ+

 

 

 

Σ+

 

60

3

 

 

 

 

2

 

 

 

 

=

 

 

πR

2R = 3πR

,

2

так как эллипс имеет полуоси R

 

и R и его площадь, как известно, равна

2

произведению полуосей на число π.

 

 

 

Упражнения. Используя формулу Стокса решить задачи:

 

 

 

 

 

 

~ 2~ 2~ 2~

4.33. Найти циркуляцию вектора f = z i+x j +y k по сечению сферы

x2 + y2 + z2 = R2 плоскостью x + y + z = R в положительном направлении.

~

3~

3~

3~

4.34. Найти циркуляцию вектора f = z

i

+ x j + y

k по сечению ги-

перболоида 2x2 − y2 + z2 = R2 плоскостью x + y = 0 в положительном

направлении.

2~

 

 

2

 

2 ~

 

~

~

 

+ y

 

4.35. Найти циркуляцию вектора f = y

i + xyj + (x

 

)k по кривой,

вырезаемой в первом октанте из параболоида x2 + y2

= Rz плоскостями

x = 0, y = 0, z = R, в положительном направлении.

 

 

 

 

 

~

 

~

 

 

~

~

~

4.36. Найти поток вихря поля (rot f) вектора f = yi + zj + xk через

часть поверхности z = 2(1 − x2 − y2), лежащей над плоскостью 0xy, в

направлении внешней нормали.

 

 

 

~

~

~

~

4.37. Найти циркуляцию вектора f = (y + z)i

+ (z + x)j + (x + y)k

вдоль эллипса l = {x = a sin2 t; y = 2a sin t cos t; z = a cos2 t, 0 ≤ t ≤ π},

пробегаемого в направлении возрастания t.

 

 

 

~

~

~

~

4.38. Найти циркуляцию вектора f = (y+z)i+(x+z)j+(x+y)k вдоль

окружности, являющейся сечением сферы x2 + y2 + z2 = R2 плоскостью

x + y + z = 0, в положительном направлении.

R3

 

4

πR3.

 

3

πR4. 4.35.

Ответы: 4.33.

 

4.34.

 

 

 

. 4.36. −π. 4.37. 0.

3

2

3

4.38.0.

5.ЭЛЕМЕНТЫ ТЕОРИИ ПОЛЯ

5.1. Оператор Гамильтона и его применение

Все операции векторного анализа можно выразить при помощи оператора Гамильтона – символического вектора r, определяемого равенством

~

~

~

 

r = i

∂x

+ j

∂y

+ k∂z .

Применяя известные операции умножения вектора на скаляр, скалярного и векторного произведений двух векторов, находим:

 

 

~

∂ϕ

~

∂ϕ

~

∂ϕ

= rϕ;

grad ϕ = i

∂x

+ j

∂y

+ k

∂z

∂ϕ

= (~e, rϕ) = (~e, r)ϕ;

 

 

 

∂~e

61

 

 

 

 

 

 

 

~

 

∂fx

 

∂fy

 

 

 

∂fz

 

 

~

 

 

 

 

 

 

 

 

 

div f =

 

 

 

+

 

 

 

 

+

 

 

 

 

= (r, f );

 

 

 

 

 

rot f~ = ∂yz

∂x

∂y

 

∂z

 

∂yx ~k =

 

 

 

− ∂zy ~i +

 

∂zx

∂xz

~j +

∂xy

 

 

 

 

 

 

∂f

 

 

∂f

 

 

 

 

 

∂f

 

 

 

 

∂f

 

 

 

 

∂f

∂f

 

 

 

 

 

 

 

 

 

 

 

~i

 

~j

 

~k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= det

 

 

 

= r × f~ = [r, f~].

 

 

 

 

 

 

 

 

∂x

∂y

 

∂z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fx

 

fy

 

fz

направлению от скалярной функции

 

∂ϕ

По аналогии с производной

по

 

 

вводится понятие производной по направлению единичного вектора ~e

 

∂~e

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

от векторной функции f(x, y, z):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂f

~

~

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

~

 

 

 

 

 

 

∂~e

= (~e, r)f

= i(~e, grad fx) + j(~e, grad fy) + k(~e, grad fz) =

 

 

 

 

 

 

 

 

 

 

 

∂fx ~

 

 

∂fy ~

 

 

∂fz ~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

i +

 

 

j +

 

 

k.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂~e

∂~e

 

∂~e

 

 

 

 

 

 

 

С помощью оператора Гамильтона удобно выполнять дифференциаль-

ные операции векторного анализа над сложными выражениями.

 

 

 

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

 

 

Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

grad(ϕg) = r(ϕg) = r(ϕ g) + r(ϕ g)

 

 

 

( указывает функцию, на которую “действует” оператор). Но

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r(ϕ g) = grϕ = g grad u,

 

 

r(ϕ g) = ϕrg = ϕ grad g.

 

Таким образом, grad(ϕg) = g grad ϕ + ϕ grad g.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 5.2. Найти rot[f,~c ], где ~c – постоянный вектор.

 

 

 

По известной формуле векторной алгебры h~a, [~b,~c ]i = (~a,~c )~b −(~a,~b )~c,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~0, имеем rot[f~,~c ]

 

hr, [f~,~c ]i

 

учитывая соотношение r, [f~,~c ]

=

 

=

=

"#

 

 

 

 

 

 

= r, [f~,~c ]

+ r, [f~,~c ] = (r,~c ) f~ ) − (r, f~ )~c. Но (r,~c ) f~ = (~c, r) f~.

~

~

~

 

~

~

Далее, (r, f)~c = ~c(r, f) = ~c div f. Таким образом, rot[f,~c ] = (~c, r)f−

~

 

 

 

 

 

−~c div f.

 

 

 

 

 

Можно образовать 5 дифференциальных операций второго порядка:

1) div grad ϕ = (r, r)ϕ = r2ϕ =

ϕ (лапласиан функции ϕ);

 

62

2)

rot

~

 

 

r r ~; 3)

 

 

~

r r

~

; 4)

r [r

~

] ;

 

 

grad ϕ = [ , ]ϕ grad div f =

(

, f )

 

div rot f = ,

, f

 

5)

rot rot f =

 

r, [r, f ] .

r

2

 

 

 

 

 

 

 

 

 

 

того, операцию

2 ~

 

 

 

 

 

 

 

Кроме

 

h

i

 

 

 

 

 

 

 

 

 

т. е. рассматривать операцию r f. Вторая и четвертая операции приводят к нулевому вектору и нулю соответственно:

r r ≡ ~ ~ r r~ ≡ rot grad ϕ = [ , ]ϕ 0, div rot f = , [ f ] 0.

Это следует из векторного смысла оператора r: в первом случае формально имеем векторное произведение двух коллинеарных векторов, а во втором

– смешанное произведение компланарных векторов. Упражнения. Доказать справедливость следующих формул:

 

~

~

~

 

~

~

~

5.1. div(ϕf ) = ϕ div f + f grad ϕ.

5.2. div[f, ~g ] = ~g rot f

− f rot ~g.

~

~

 

~

 

~

~

~

5.3. rot(ϕf ) = ϕ rot f + [grad ϕ, f ].

5.4. grad([f, ~g ]) = (f, r)~g + (~g, r)f+

~

~

 

~

~

~

~

~

+[f, rot ~g ] + [~g, rot f ].

5.5. rot[f, ~g ] = (f, r)~g + (~g, r)f + f div ~g −~g div f.

 

~

~

~

 

 

 

 

5.6. rot rot f = grad div f −

f.

 

 

 

 

5.2. Потенциальное поле

В физике принято называть функции, заданные в R3, термином “поле”. Будем также использовать эту терминологию: если рассматривается функция f : R3 → R, будем говорить о скалярном поле, если функция

~

R

3

→ R

3

– о векторном поле.

 

 

 

f :

 

 

 

 

 

 

 

 

 

 

~

~

~

~

 

Определение 5.1 [1]. Векторное поле f

= fxi

+ fyj

+ fzk называ-

ется потенциальным (в области Ω), если существует такая непрерывно дифференцируемая на Ω скалярная функция ϕ: R3 → R, что в каждой точке Ω выполнено равенство

~

(5.1)

f = grad ϕ.

~

Скалярное поле ϕ называется при этом потенциалом поля f.

~

Теорема 5.1 [1]. Если поле f потенциально в области Ω и потенциал ϕ является дважды непрерывно дифференцируемой в Ω функцией,

то

 

 

~

~

(5.2)

rot f = 0.

Теорема 5.1 означает, что равенство (5.2) является необходимым условием потенциальности гладкого поля. Это же условие оказывается и до-

~

статочным для того, чтобы непрерывно дифференцируемое поле f имело потенциал.

63

~

Теорема 5.2 [1]. Если векторное поле f непрерывно дифференциру-

~

емо в односвязной области Ω и выполнено условие (5.2), то f является потенциальным полем, т. е. существует потенциал ϕ, удовлетворяющий равенству (5.1). Если P0 – некоторая фиксированная точка области Ω, потенциал ϕ определяется единственным образом и может быть вычислен по формуле

P

 

 

ϕ(P ) = ϕ(P0) + Z (f~, dl~).

 

(5.3)

P0

 

 

~

 

 

В физических задачах чаще всего вектор f является силой, при этом

говорят о силовом поле. В этом случае интеграл

Z

(f~, dl~) называется

l+(P1,P2)

~

работой A силового поля f вдоль пути l+(P1, P2) с началом в точке P1 и концом в точке P2. Для потенциального поля работа поля зависит от точек P1 и P2 и не зависит от пути, соединяющего эти точки, т. е. в этом случае A = A(P1, P2). Из равенства (5.3) следует, что для потенциального поля A(P, P0) = ϕ(P ) − ϕ(P0) и, значит (из-за аддитивности криволинейного интеграла), A(P1, P2) = ϕ(P2) − ϕ(P1) для любых точек P1 и P2, лежащих в области Ω. Это известный физический закон: работа поля равна разности

потенциалов конечной и начальной точек.

 

 

 

 

2

 

 

 

 

2~

 

 

 

 

 

 

 

 

 

 

 

~

~

 

 

 

 

~

Пример 5.3. Найти потенциал поля f = 2xyi + (x

 

2yz)j

− y k.

Убедимся сначала, что поле потенциально:

 

 

 

 

 

 

 

 

 

 

 

∂fz

 

∂fy

= −2y,

∂fx

∂fz

∂fy

 

 

∂fx

 

 

 

 

 

=

 

 

 

 

=

 

= 0,

 

=

 

 

= 2x.

 

 

∂y

 

∂z

∂z

∂x

∂x

∂y

 

 

 

 

~

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, rot f

≡ 0. За путь интегрирования примем ломаную 0ABP ,

где 0(0, 0, 0), A(X, 0, 0), B(X, Y, 0), P (X, Y, Z). Находим:

 

 

 

 

 

 

 

 

 

Z

 

 

A

 

 

 

B

 

 

 

 

 

P

 

 

 

ϕ(X, Y, Z) =

 

(f~, dl~) + c = Z (f~, dl~) + Z (f~, dl~) + Z (f~, dl~) + C,

 

 

 

0ABP

 

0

 

 

 

A

 

 

 

 

 

B

 

 

 

 

 

 

~

~

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

(f, dl ) = 2xy dx + (x

 

− 2yz)dy − y

dz.

 

 

A

 

 

Так как на 0A имеем y = z = 0, dy = dz = 0, следовательно, Z0

(f~, dl~) = 0.

Аналогично на AB имеем x = X, dx = 0, z = 0, dz = 0, поэтому

BY

ZZ

~ ~

2

2

Y.

(f, dl ) = X

dy = X

A

0

 

 

64

На BP имеем x = X, y = Y , dx = dy = 0, значит,

P Z

ZZ

~ ~

Y

2

dz = −Y

2

Z.

(f, dl ) = −

 

 

B

0

 

 

 

 

Таким образом, ϕ(X, Y, Z) = X2Y − Y 2Z + C.

Упражнения. Найти потенциалы следующих плоских и трехмерных полей:

 

~

 

 

 

2

y − y

3 ~

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

3

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.7. f = (3x

 

)i + (x

 

 

 

− 3xy )j.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

sin(2x) cos(2y)i + cos(2x) sin(2y)j

 

 

 

 

 

 

 

 

 

 

 

 

 

5.8. f =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

pcos ~

sin

2 y

+

 

 

x2

 

 

 

2

~

 

 

 

 

 

 

 

 

2 ~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 x

 

 

 

 

 

 

 

 

 

 

 

sin2 x cos2 y

 

 

 

 

 

 

 

 

 

 

 

 

 

5.9. f = (yz − xy)i + xz −

 

 

 

 

 

 

+ yz j + (xy + y z)k.

 

 

 

 

 

 

2

 

 

 

 

 

 

 

5.10. f~ =

z

x2

~i + x

y2

~j +

y z2

 

~k.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

z

 

 

 

 

 

 

1

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

5.11. f~ =

y2

z2

 

2x3

~i+ x2

z2

2y3

 

~j+ x2 + y2

+

2z3

~k.

 

 

 

 

 

 

z

 

 

 

y

 

 

 

 

 

 

 

yz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

x

 

 

 

 

 

xz

 

 

 

 

y

 

x

 

xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

m

 

 

 

 

 

 

5.12. Найти потенциал гравитационного поля f

= −

r

3

~r, создавае-

мого массой m, помещенной в начале координат.

 

 

 

 

k k

 

 

 

 

 

 

Ответы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.7. ϕ(x, y) = x3y − xy3 + c.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.8.

(

 

 

 

) = 2p

cos

 

 

 

x2y

 

 

 

 

 

 

y2z2

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

ϕ

x, y

 

 

 

 

 

 

 

 

 

 

 

 

 

2 x sin2 y + sin2 x cos2 y + c

 

 

 

 

 

 

 

 

 

 

 

 

5.9. ϕ(x, y, z) = xyz −

 

 

 

 

 

+

 

 

 

 

 

+ c.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.10. ϕ(x, y, z) =

y

+

z

 

+

x

 

+ c.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

y

 

 

 

 

 

 

xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yz

 

 

 

 

 

 

xz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.11. ϕ(x, y, z) =

 

 

 

+

 

 

 

 

+ c.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

y2

 

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

5.12. ϕ(x, y, z) = p . x2 + y2 + z2

5.3. Соленоидальное поле Определение 5.2 [1]. Непрерывно дифференцируемое векторное по-

~

~

~

 

ле ~v = vxi + vyj + vzk называется соленоидальным в области Ω, если

 

 

div ~v = 0

(5.4)

в каждой точке Ω.

65

Определение 5.3 [1]. Если для поля ~v

существует такая непре-

~

 

рывно дифференцируемая функция A, что

 

~

(5.5)

rot A = ~v

~

в области Ω, то A называется векторным потенциалом поля ~v.

Важным является следующее утверждение.

Теорема 5.3 [1]. Если поле ~v имеет в области Ω дважды непрерыв-

~

но дифференцируемый векторный потенциал A, то ~v соленоидально в Ω.

Аналогично случаю потенциального поля не только из равенства (5.5) следует соотношение (5.4), но и наоборот, из условия (5.4) следует суще-

~

ствование вектора A, удовлетворяющего (5.5). Это означает, что любое соленоидальное поле имеет векторный потенциал.

Пример 5.4. Доказать, что для любого дважды дифференцируемого

~

трехмерного векторного поля f поле вихрей соленоидально. Имеем:

rot f~ =

∂yz

∂zy ~i +

∂zx

∂xz

~j +

∂xy

∂yx

~k.

 

 

∂f

 

∂f

∂f

 

∂f

 

∂f

 

∂f

 

Учитывая равенство смешанных производных второго порядка, получаем

 

 

 

 

 

 

 

∂f

 

∂f

 

 

 

 

 

∂f

∂f

 

 

∂f

 

 

∂f

 

 

 

 

div rot f~

=

 

 

 

 

z

 

 

 

y

+

 

 

 

x

z

 

+

 

 

 

 

y

 

 

x

≡ 0.

∂x

∂y

 

∂z

∂y

∂z

∂x

∂z

∂x

 

∂y

Упражнения. Проверить соленоидальность следующих полей:

 

 

 

 

~

 

 

 

2

 

 

 

3

~

3

 

 

2

~

 

~

2~

2 ~

 

 

 

 

2

 

 

2

~

5.13. f = (x

y + y )i + (x

−xy )j. 5.14. f = xy i + x yj

−(x

+ y

)zk.

 

 

 

x

 

 

 

 

y

 

 

 

 

 

(x + y) ln z

 

 

 

 

 

 

~

~

 

 

 

(x

2

 

 

2

 

~

 

f~ =

 

~i +

 

 

~j

 

 

 

~k

 

 

 

f~ =

xi − yj

+

 

− y

 

)zk

5.15.

 

 

 

 

 

 

 

. 5.16.

 

 

(x2

+ y2)3/2 .

 

yz

 

xz

 

 

 

xy

 

 

 

 

 

 

px2 + y2

 

5.4.Применение криволинейных координат в векторном анализе

Напомним, что для

дифференцируемого скалярного поля f : R3

R

 

 

 

 

 

3

 

 

 

 

в том случае, когда положение точки в R

 

описывается декартовыми ко-

ординатами, вектор grad f определяется равенством

 

 

 

 

∂f

~

∂f

~

 

∂f

~

 

 

grad f = ∂x i +

∂y j + ∂z k.

 

 

Пусть теперь в R3 заданы криволинейные ортогональные координаты (ξ, η, ζ). Как показано в [1], в этом случае в каждой точке P определен ко-

 

1

 

∂x

 

∂y

 

∂z

 

T

ординатный ортогональный базис {~eξ, ~eη, ~eζ }, где ~eξ =

 

 

,

 

,

 

;

Hξ

∂ξ

∂ξ

∂ξ

66

 

1

 

∂x

 

 

∂y

 

∂z

 

T

1

 

∂x

 

∂y

 

∂z

 

T

~eη =

 

 

,

 

 

,

 

; ~eζ =

 

 

 

,

 

,

 

– представления базисных

Hη

∂η

∂η

∂η

 

Hζ

∂ζ

∂ζ

∂ζ

 

 

 

 

 

 

 

 

 

 

 

~ ~ ~

 

 

 

 

 

 

векторов в исходном базисе {i, j, k}, а Hξ, Hη, Hζ – коэффициенты Ламе,

определенные в [1]. В [1] было получено представление градиента в криволинейных ортогональных координатах:

1

 

∂f

1

 

 

 

∂f

1

 

 

 

∂f

 

grad f =

 

 

 

 

 

 

~eξ +

 

 

 

 

 

 

 

 

 

~eη +

 

 

 

 

 

 

 

 

 

~eζ .

Hξ

 

∂ξ

Hη

∂η

Hζ ∂ζ

В частности, в цилиндрической системе координат

 

 

 

 

 

 

∂f

 

1 ∂f

∂f

 

 

 

 

 

 

 

grad f =

 

 

 

~eρ +

 

 

 

 

 

 

~eϕ +

 

 

 

~ez,

 

 

 

∂ρ

ρ

∂ϕ

∂z

 

а для сферической системы координат

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂f

1

 

 

 

 

 

∂f

1

 

 

∂f

 

grad f =

 

 

~eρ +

 

 

 

~eϕ +

 

 

 

~eθ.

∂ρ

 

ρ sin θ

∂ϕ

ρ

∂θ

В [1] было получено представление дивергенции в криволинейных ор-

тогональных координатах:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div f~ =

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(HξHηfζ ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(HηHζ fξ) +

 

 

(HξHζ fη) +

 

 

 

 

HξHηHζ

∂ξ

∂η

∂ζ

 

В частности, для цилиндрической системы координат

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 ∂

 

 

 

 

 

 

 

 

 

 

 

∂fϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 ∂

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 ∂fϕ

 

 

∂fz

 

div f~ =

 

 

 

 

 

 

(ρfρ) +

 

 

 

+

 

 

(ρfz) =

 

 

 

 

 

 

 

(ρfρ) +

 

 

 

 

 

 

+

 

 

,

ρ

∂ρ

∂ϕ

∂z

ρ

∂ρ

ρ ∂ϕ

∂z

а для сферической системы координат

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div f~ =

 

 

 

 

1

 

 

 

 

 

 

 

 

2 sin θfρ) +

 

 

(ρfϕ) +

 

(ρ sin θfθ) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ2 sin θ

∂ρ

 

 

∂ϕ

∂θ

 

=

 

1

 

 

2fρ) +

1

 

 

 

∂fϕ

+

1

 

 

 

 

(sin θfθ) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ2

∂ρ

ρ sin θ ∂ϕ

 

ρ sin θ ∂θ

 

 

 

 

 

так как Hρ = 1, Hϕ = ρ sin θ, Hθ = ρ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В [1] было получено представление rot f в цилиндрических и сфери-

ческих координатах соответственно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

∂f

 

 

 

 

∂f

 

 

 

 

 

 

∂f

 

 

 

 

 

∂f

 

 

 

 

 

 

 

 

1

 

 

 

∂(ρf

)

 

 

 

∂f

 

rot f~ =

 

 

 

 

z

 

 

ρ

 

ϕ

 

~eρ

 

+

ρ

 

 

z

~eϕ +

 

 

 

 

 

 

 

 

ϕ

 

 

 

ρ

~ez,

ρ

∂ϕ

 

 

∂z

 

∂z

 

 

∂ρ

ρ

 

 

 

 

∂ρ

 

 

∂ϕ

 

 

 

 

 

 

1

 

 

 

 

 

 

∂f

 

 

 

∂(sin θf

)

 

 

 

 

 

1

 

 

∂f

 

 

 

 

 

∂(ρf )

 

 

 

 

 

rot f~ =

 

 

 

 

 

θ

 

 

 

 

 

 

ϕ

 

~eρ

+

 

 

 

ρ

 

 

 

θ

 

~eϕ +

 

 

 

ρ sin θ

 

∂ϕ

 

 

 

∂θ

 

 

 

 

 

ρ

 

∂θ

 

 

∂ρ

 

 

67

 

 

 

 

 

 

 

+ρ sin θ

sin θ

 

∂ρϕ

)

∂ϕρ ~eθ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

∂(ρf

 

 

∂f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 5.5. Перейти к цилиндрическим координатам в выражении

 

 

 

 

 

 

 

~

 

~

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

векторного поля f~ =

xi + yj

− zk

 

и найти div f~ и rot f~.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 + y2 + z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

случае xi + yj = ~r, то

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как в данномp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f~ =

ρ~eρ − z~ez

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ2 + z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По полученным формулам

находим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

+ ρ ∂zz

 

 

 

 

 

 

 

 

 

 

 

 

 

div f~ = ρ

 

∂ρρ

+

∂ϕ

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

∂(ρf

)

 

 

∂fϕ

 

 

∂f

 

 

 

 

 

 

 

 

 

 

 

=

1

 

 

2ρ(ρ2 + z2) − ρ3

ρ

2 + z2) − z2

 

=

 

2z2

 

;

 

 

 

 

 

 

 

 

2 + z2)3/2

 

 

ρ

2 + z2)3/2

 

 

 

 

2 + z2)3/2

 

 

 

~ez =

rot f~ = ρ ∂ϕz

∂zϕ

~eρ +

∂zρ

∂ρz

~eϕ + ρ

∂ρϕ

)

∂ϕρ

 

1 ∂f

 

∂f

 

 

 

 

 

 

∂f

 

∂f

 

 

 

 

 

1

 

 

∂(ρf

 

∂f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ρz

= −2 + z2)3/2 ~eϕ.

Упражнения.

5.17. Перейти к цилиндрическим координатам

~ ~ ~ − p 2 2~ f = xzi + yzj z x + y k

~~

инайти div f и rot f.

5.18.Дано скалярное поле в цилиндрических координатах

f(ρ, ϕ, z) = ρ2ϕ + z2ϕ3 − rϕz + 3.

Найти grad f.

5.19. Дано векторное поле в цилиндрических координатах

~

 

2

+ ρ)~eρ + zρϕ~eϕ + (z

2

+ ρϕ)~ez.

f(ρ, ϕ, z) = (z

 

 

~

~

 

 

 

 

Найти div f и rot f.

5.20. Дано скалярное поле в сферических координатах f(ρ, ϕ, θ) = ρ2ϕθ + ρ3ϕ2 + ϕ + θ2 + 1.

Найти grad f.

5.21. Дано векторное поле в сферических координатах

~

2

2

3

)~eϕ + (ρ

4

+ 1)~eθ.

f(ρ, ϕ, θ) = ρ

ϕ~eρ + (ρϕ

θ + ρ

 

68

 

 

 

~

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найти div f и rot f.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.17.

~

 

ρz ~e

~e

 

 

 

 

 

~

 

 

 

 

 

 

 

~

 

ρ

 

z ~e

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

=

 

( ρ

 

z), div f = 2z − ρ, rot

2

= (

 

 

+

) ϕ.

 

 

 

 

 

 

5.18. grad f = (2ρϕ − ϕz)~eρ + ρ +

3z2ϕ

 

− ρz

~eϕ + (2zϕ3 − ρϕ)~ez.

 

 

 

ρ

 

 

 

 

 

 

 

 

 

~

 

z2

 

 

 

 

 

2

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.19. div f =

ρ

 

+ z + 3z

 

+ 2,

rot f = (1 − ρϕ)~eρ + (2z − ϕ)~eϕ + 2zϕ~ez.

 

5.20. grad f = (2ρϕθ+3ρ2ϕ2)~eρ+

 

1

 

 

2θ+2ρ3ϕ+1)~eϕ+

1

2

ϕ+2θ)~eθ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ cos θ

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

~

2

 

 

 

 

2ϕθ

 

 

4 + 1) tg θ

 

 

~

 

 

 

2

 

2

 

2

 

 

5.21. div f = ρ

ϕ +

cos θ

 

 

 

ρ

 

 

 

, rot f = [tg θ(ϕ

θ + ρ

)

−ϕ

]~eρ

 

4 + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~eϕ

+ h2θ + 4ρ2

 

i~eθ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ cos θ

cos θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соседние файлы в папке Мат Анал