Беляев Ю.Н. Введение в векторный анализ
.pdfŸ 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 Хэвисайд из-за сходства знака
с остовом древнего струнного инструмента с таким названием.
Ÿ 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 |
|
|
|
|