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

Беляев Ю.Н. Введение в векторный анализ

.pdf
Скачиваний:
401
Добавлен:
29.03.2016
Размер:
1.13 Mб
Скачать

Ÿ 4. Расхождение векторного поля

171

; 2) через поверхность куба, ограниченного плоскостями x = 0; x = 1; y = 0; y = 1; z = 0; z = 1:

96. ~

Вычислить поток векторного поля P (M ) = 4xz~ex y2~ey + yz~ez через поверхность куба, ограниченного плоскостями x = 0; x = 1; y = 0; y = 1; z = 0; z = 1:

Ÿ 4. Расхождение векторного поля

~

Величина элементарного потока P ~end положительна, если

~

угол между векторами P и ~en острый, и отрицательна, если эти векторы образуют тупой угол.

Следовательно, поток (3.48) представляет избыток ½жидкости\, протекающей в сторону положительной нормали ~en, а не абсолютное количество ½жидкости\, прошедшей через поверхность независимо от направления течения.

Поток N , вычисленный для замкнутой поверхности, даст разность между ½жидкостью\, вытекающей из объ¼ма V , ограниченного этой поверхностью, и ½жидкостью\, поступающей в этот объ¼м. Если поток через замкнутую поверхность равен нулю N = 0, то это означает равенство втекающей в объ¼м и вытекающей из объ¼ма ½жидкости\. Если N > 0, то в объ¼ме имеются источники поля, т.е. такие места где поле как-то созда¼тся. Если поток отрицательный, то в объ¼ме есть стоки, где поле как-то уничтожается.

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

~

Рассмотрим отношение потока векторного поля P через замкнутую поверхность к величине объ¼ма V , ограниченного поверхностью ,

1

ZZ

~

 

 

~

(3.55)

V

P d :

Если существует конечный предел отношения (3.55), когда объ¼м V стягивается по произвольному закону к некоторой точке M 2 V , так что площадь поверхности и

172

 

 

 

 

 

Глава 3.

ФУНКЦИИ ТОЧКИ

объ¼м V стремятся

ê íóëþ,

òî

ýòîò

предел

называется

ä è â å ð ã å í ö è å é

4

 

 

 

 

 

 

 

 

 

~

 

, или р а с х о ж д е н и е м, векторного поля P

 

 

 

 

 

~

:

 

 

 

 

 

в точке M и обозначается div P

 

 

 

 

 

 

~

lim

1

ZZ

~

~

dN

(3.56)

 

 

 

 

 

 

div P = V !M V

P d = dV :

 

Здесь dN это поток векторного поля через элементарную замкнутую поверхность, ограничивающую элементарный объ¼м dV .

~

Величина div P есть локальная, точечная характеристика

~

распределения источников и стоков в данном поле P . Согласно определению (3.56), дивергенция векторного поля,

вычисленная в точке M , приближ¼нно равна потоку, выходящему из единицы объ¼ма, окружающего эту точку.

Отметим два очевидных свойства расхождения для суммы полей и произведения векторного поля на произвольную постоянную k:

~

~

~

~

;

(3.57)

div (P1

+ P2) = div P1

+ div P2

 

~

~

 

(3.58)

div kP = k div P :

 

4.1. Расхождение в ортогональной системе координат Построим на координатных линиях q1, q2, q3 криволинейной ортогональной системы координат, проходящих через точку M , элементарный параллелепипед M N KLM 0N 0K0L0 (ñì. ðèñ. 91).

Если координаты точки M обозначить через q1M , q2M , q3M , то гранями рассматриваемого параллелепипеда являются участки координатных поверхностей:

q1 = q1M ; для грани M N KL; q1 = q1M + q1; для грани M 0N 0K0L0;

q2 = q2M ; для грани M M 0L0L; q2 = q2M + q2; для грани N N 0K0K;

4От латинского слова divergence расхождение, этот термин был пред-

ложен Клиффордом в 1878 г.

Ÿ 4. Расхождение векторного поля

173

q3 = q3M ; для грани M M 0N 0N ; q3 = q3M + q3; для грани LL0K0K:

q3

K

K0

Рассмотрим поток векторного поля

 

 

 

 

 

 

~

 

 

q2

P через поверхность элементарного

 

 

параллелепипеда M N KLM 0N 0K0L0.

L

 

L0

 

 

 

Для грани M 0N 0K0L0, поверхность

 

N

N 0

которой обозначим 1, единичным

 

 

 

вектором внешней нормали к объ-

M

M 0

q1

¼му параллелепипеда является пер-

 

 

вый базисный вектор ~e1 системы ко-

 

Ðèñ. 91

 

 

 

ординат.

Площадь элементарной площадки, перпендикулярной первой координатной оси, выражается через коэффициенты Ламэ по формуле d 1 = H2H3dq2dq3. Следовательно, поток векторного поля через эту грань, используя теорему о среднем, можно

представить в виде

 

 

N1 = ZZ

P~ ~e1d = ZZ

P1H2H3dq2dq3 =

1

1

 

= P1H2H3(q1M + q1; q2; q3) q2 q3;

ãäå P1H2H3(q1M + q1; q2; q3) обозначает значение функций P1, H2, H3 (их произведения) в некоторой ½средней\ точке поверхности 1 с координатами q1M + q1, q2 è q3.

Аналогично, для поверхности 2 с вершинами M , N , K, L единичным вектором внешней нормали является ~e1, и поток равен

N2 = ZZ

P~ ~e1d = ZZ

P1H2H3dq2dq3 =

2

2

 

= P1H2H3(q1M ; q2; q3) q2 q3:

Строго говоря, усредн¼нные значения q2 è q3 в формуле для потока N2 могут быть не теми же самыми, что для потока N1.

174

Глава 3. ФУНКЦИИ ТОЧКИ

Но чем меньше будут площади поверхностей 1 è 2, тем меньше возможные различия для второй и третьей координат усредн¼нных точек на этих поверхностях.

По формуле Лагранжа

P1H2H3(q1M + q1; q2; q3) P1H2H3(q1M ; q2; q3) =

= @(P1 H2H3)

 

q1;

 

 

 

 

 

 

 

 

 

 

 

@q1

q1=q1

ãäå q1 2 (q1M ; q1M + q1), поэтому суммарный поток векторного поля через указанные две грани

N1 + N2 = @(P1H2H3)

 

q1 q2 q3:

(3.59)

 

 

 

 

 

 

 

 

 

 

 

@q1

q1=q1

 

Точно так же можно показать, что поток векторного поля

~

P через две грани параллелепипеда, перпендикулярные линии q2, равен

@(P2H3H1)

 

 

 

 

 

 

q1

q2 q3;

ãäå q (q2M ; q2M + q2);(3.60)

 

 

 

 

2 2

@q2

q2=q2

 

~

а поток P через две грани параллелепипеда, перпендикулярные оси q3, åñòü

@(P3H1H2)

 

 

 

 

 

 

 

q q q ;

ãäå q (q ; q

+ q ):(3.61)

 

 

 

 

3 2 3M 3M

 

@q3

q3=q3

1 2 3

3

~

Вычислим теперь расхождение векторного поля P в точке M . Для этого сложим потоки (3.59), (3.60) и (3.61), найд¼м предельное значение полученного выражения, когда объ¼м параллелепипеда стягивается в точку M , и, наконец, разделим его на элементарный объ¼м

 

 

 

dV = H1H2H3dq1dq2dq3:

 

 

В результате получаем

 

 

 

 

 

1

 

@(P1H2H3)

 

@(P2H3H1)

 

@(P3H1H2)

div P~ =

 

 

 

+

 

+

 

: (3.62)

H1H2H3

@q1

@q2

@q3

Ÿ 4. Расхождение векторного поля

175

П р и м е р 75. Р а с х о ж д е н и е в д е к а р т о в о й с и с т е м е к о о р д и н а т. Коэффициенты Ламе: Hx = 1; Hy = 1; Hz = 1:

~

@Px

 

@Py

 

@Pz

 

(3.63)

div P =

@x

+

@y

+

@z

:

П р и м е р 76. Вычислить div ~r.

div ~r = @x@x + @y@y + @z@z = 3 :

4.2.Соленоидальное векторное поле. Векторный

потенциал.

 

~

~

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

, для которого div P

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

П р и м е р 77. Показать, что кулоновское поле

 

~

~r

 

 

 

E = k

r3

 

=0, íà-

(3.64)

является соленоидальным. Постоянная k в формуле (3.64) зависит от системы единиц измерения.

Для вычисления дивергенции подставим компоненты поля (3.64) в формулу (3.63). Находим

div E~ =k

@(xr 3)

+

@(

yr 3)

+

@(zr 3)

! = 0:

(3.65)

@x

 

 

@y

 

@z

Любая вектор-функция

~

 

 

 

которой расхождение

P (M ), äëÿ

равно нулю, может рассматриваться как вихрь некоторого дру-

~

~

называется векторным потенциа-

гого вектора W . Вектор W

~

ëîì ïîëÿ P (M ).

4.3.Лапласово векторное поле. Рассмотрим вектор-

~

ное поле P , которое одновременно является потенциальным и соленоидальным, т.е. вихрь и расходимость этого поля равны нулю:

~

~

~

(3.66)

rot P

= 0;

div P = 0:

5От греческого слова o трубка.

176

Глава 3. ФУНКЦИИ ТОЧКИ

~

Если поле потенциально, то в односвязной области P = grad U . Тогда

~

div P = div grad U = 0:

Последнее равенство есть уравнение Лапласа, которое в декартовой системе координат имеет вид

@2U @2U

 

@2U

 

(3.67)

 

+

 

+

 

= 0:

@x2

@y2

@z2

Поэтому векторное поле, удовлетворяющее условиям (3.66), называется лапласовым и в односвязной области полностью определяется скалярным потенциалом, удовлетворяющим уравнению Лапласа.

Функция U , непрерывная вместе с производными первого и второго порядка и удовлетворяющая уравнению Лапласа, называется гармонической.

П р и м е р 78. Функция

r ; ãäå r = qx2

+ y2 + z2;

1

 

 

 

 

является гармонической всюду, кроме начала координат. Эта функция фундаментальное решение уравнения Лапласа в тр¼хмерном пространстве.

Фундаментальным решением уравнения Лапласа на плоскости называется функция

q

ln x2 + y2:

Она также является гармонической на плоскости всюду, кроме начала координат.

Задачи

 

97.

 

 

 

 

 

~

 

 

Вычислить дивергенцию векторного поля: 1)P (M ) =

3

 

3

3

~

 

 

 

 

x

yz~ex + xy

z~ey + xyz

~ez ; 2)P (M ) = (3x y)~ex +(x 2y)~ey +(2y +

z)~ez ; 3)

~

~

~

2

~

2

~r.

P (M ) = ~r; 4)P (M ) = r ~r; 5)P (M ) = ~r=r

 

; 6)P (M ) = r

Ÿ 5. Символические обозначения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

177

98. Показать, что поле P = y

 

z

 

~ex

+ 2x

 

z

 

~ey 3x

 

y

 

~ez ÿâ-

~

4

 

2

 

 

3

 

2

 

 

2

 

 

2

 

ляется соленоидальным.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

xyz~ez â

99. Для векторного поля P (M ) = xyz~ex + 2xyz~ey

точке M (1; 1; 1) определить: 1) уравнение векторной линии; 2)

расхождение поля; 3) вихрь поля.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100. Проверить, являются ли гармоническими следующие

функции:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'1 = c; '2 = x2 z2; '3 = xy;

'4 = ax + by + cz;

'5 = xyz;

'6 = (x2 + y2)k=2 a sin k arctan x

+ b cos k arctan x :

 

 

 

y

 

 

 

 

 

 

 

 

y

 

Ÿ5. Символические обозначения основных дифференциальных операций

5.1.Символический вектор набла6 это векторный

оператор

@

@

@

 

r = ~ex

 

+ ~ey

 

+ ~ez

 

;

@x

@y

@z

который вв¼л в рассмотрение Гамильтон в 1853 году, поэтому r еще называют оператором Гамильтона. ½Координатами\ этого вектора в декартовой системе координат являются символы @x@ ,

@y@ , @z@ .

При выполнении действий по правилам, установленным ранее для обычных скаляров и векторов, под ½произведением\ символов @x@ , @y@ , @z@ на скаляр ' и вектор ~a будем понимать соответственно

скаляры

@'

;

@'

;

@'

и векторы

@~a

;

@~a

;

@~a

:

@x

@y

@z

@x

@y

@z

 

 

 

 

 

 

 

Применяя оператор r к скалярному полю U (x; y; z), полу- чаем

rU = ~ex

@U

+ ~ey

@U

+ ~ez

@U

= grad U:

(3.68)

 

 

 

@x

@y

@z

6Термин ½набла\ оператору r дал в 1892 Хэвисайд из-за сходства знака

с остовом древнего струнного инструмента с таким названием.

178

Глава 3. ФУНКЦИИ ТОЧКИ

Скалярное произведение векторов r и ~a есть

 

 

r ~a =

 

@ax

+

@ay

+

@az

= div ~a:

 

 

 

@x

@y

 

@z

 

Векторное произведение векторов r и ~a равно:

r ~a = ~ex

@y

@z

+ ~ey @z

@x

+ ~ez @x

 

 

@az

 

@ay

 

 

 

 

@ax

 

 

@az

 

@ay

= rot ~a:

(3.69)

@a@yx =

(3.70)

Поэтому для функций grad U , div ~a, rot ~a часто используются обозначения соответственно7

 

rU;

 

r ~a;

 

 

r ~a:

 

5.2. Оператор Лапласа,

 

или лапласиан, это

 

@2

 

@2

 

 

@2

 

 

 

r2

 

 

+

 

 

 

+

 

 

 

:

 

@x2

@y2

@z2

 

Если применить его к скалярной функции U (x; y; z), то

 

r2U =

@2U

@2U @2U

(3.71)

 

 

 

+

 

 

 

+

 

 

:

@x2

 

 

@y2

 

@z2

Название оператора происходит от уравнения Лапласа

 

@2U

@2U

 

 

 

 

@2U

 

 

 

 

 

 

 

+

 

 

 

+

 

 

= 0;

 

 

@x2

 

@y2

 

 

@z2

 

которое с помощью оператора r2 записывается в виде r2U= 08. Применение оператора Лапласа к вектору ~a дает

r2~a = @2~a + @2~a + @2~a : @x2 @y2 @z2

7Тот факт, что через r можно выразить дивергенцию, ротор и градиент,

открыл в 1862 г. Тэт (Peter Guthrie Tait, 1831-1901).

8Для оператора Лапласа вместо обозначения r2 часто используют символ , и, соответственно, уравнение Лапласа выглядит так U = 0. Ïî-

следнее обозначение в данной работе не используется, поскольку ранее подU мы условились понимать приращение функции U , ñì. ñòð.134.

Ÿ 5. Символические обозначения

 

 

 

179

Òàê êàê

 

 

 

 

 

 

 

 

 

 

 

~a = ~exax + ~ey ay + ~ez az ;

 

 

то для вектора r2~a имеем представление

 

 

r2~a = ~exr2ax + ~ey r2ay + ~ez r2az :

 

5.3. Производная вектора по

другому

вектору

 

 

 

 

 

 

 

~

называется вектор

Производной вектора ~a по другому вектору b

 

@~a

 

@~a

@~a

@~a

 

(3.72)

 

 

=

 

bx +

 

by +

 

bz :

 

~

@x

@y

@z

 

 

@b

 

 

 

 

 

 

 

 

 

Проекции этого вектора на оси Ox, Oy, Oz соответственно равны

@a@xx bx + @a@yx by + @a@zx bz ; @a@xy bx + @a@yy by + @a@zy bz ; @a@xz bx + @a@yz by + @a@zz bz :

Используя определение оператора Гамильтона, выражение (3.72) можно переписать в виде

@~a

~

(3.73)

 

= (b r)~a:

~

@b

 

 

П р и м е р 79. Вычислить производную радиус-вектора ~r по постоянному вектору ~c.

Очевидно,

@~r

 

à

@~c

~

@~c

= ~c;

@~r

= 0:

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

180

Глава 3. ФУНКЦИИ ТОЧКИ

1) Оператор Гамильтона кроме векторных обладает дифференциальными свойствами, поэтому при применении r к произведению двух функций оператор должен подействовать на каждый из сомножителей, то есть r от произведения двух функций, есть сумма двух произведений, в одном из которых r действует на первый множитель, а в другом на второй. Будем отмечать функцию (скалярную или векторную), на которую действует оператор r, вертикальной стрелкой сверху, например

 

 

#

 

#

 

~

~

~

~

 

 

(3.74)

rot (~a b ) = r (~a b ) = r (~a b ) + r (~a b ):

2) Оперируя с r, как с вектором согласно алгебре векторов, нужно преобразовать полученные выражения таким образом, чтобы функция, отмеченная символом #, оказалась на последнем месте в произведении, а оператор r на предпоследнем.

Применим эти правила для вывода некоторых наиболее употребительных формул.

П р и м е р 80. В и х р ь в е к т о р н о г о п р о и з в е д е н и я. Продолжим преобразование выражения (3.74). Раскроем

двойные векторные произведения в правой части (3.74) по фор-

ìóëå:

~

~

 

~

Получим

 

 

~a (b ~c) = b (~a ~c) ~c (~a b):

 

 

#

 

#

#

 

#

#

#

~

~

~ ~

~

~

r(~a b ) + r(~a b ) =~a (rb ) b (r ~a ) +~a (r b ) b (r~a ):

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

 

 

#

 

#

#

#

~

~

~

~

~

rot (~a b ) = (b r) ~a b (r ~a ) + ~a (r b ) (~a r) b :

Используя соотношения (3.69) и (3.73), окончательно перепишем найденную формулу в виде

~

~ ~

@~a

 

~

 

 

 

@b

 

(3.75)

rot (~a b ) = ~a div b b div ~a +

 

 

 

:

~

@~a

 

 

@b