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

книги из ГПНТБ / Сахарников Н.А. Высшая математика учебник

.pdf
Скачиваний:
98
Добавлен:
27.10.2023
Размер:
24.22 Mб
Скачать

выше, придем к выводу, что div а (Л/) характеризует мощность источника или стока поля в точке М.

Поэтому точки, в которых div а (М) О 0, называются

источ­

никами векторного поля\ точки, в которых

div а (М) < 0 ,

назы­

ваются стоками векторного поля а(М).

 

 

 

 

Заметим, что иногда в физике под интенсивностью источника

(стока) понимают не div а (М), а —

div а (М).

 

 

 

Ф о р м а л ь н ы е с в о й с т в а

д и в е р г е н ц и и .

Пусть

а (М)

и b (М) — векторные функции точки и cp (М) — скалярная

функция точки. Тогда

 

 

 

 

 

 

 

 

1°. div (а

b) ~

diva-j-div Ь,

 

 

 

 

 

(24)

2°. div (ф а) = ф div а

(grad ф, а).

 

 

 

 

 

(25)

Действительно,

согласно

определению

дивергенции имеем *

1)

div (а ;-Ь):

 

і

 

дах .

 

j

dbx

 

дх

 

дх

 

г

дх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div а -- div b,

 

 

 

 

2)

аіѵ (<р«) =

І |

^ Ы . .

. _

ф ( - ^ і

-

; - .

. . )

+

 

 

 

+

( - 5 - “ ' +

 

Ч>div “

(gr»d'f.

«)■

 

 

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

diva(/H)==0. (26)

Рассмотрим произвольную векторную трубку между двумя ее произвольными сечениями и S 2 (рис. 144); поверхность рас­ сматриваемого участка трубки обозначим S. Согласно условию (26) формула Остроградского (22) дает

И ап da I j апder + 11 а,, da О.

По определению векторной трубки векторы а и п на поверх­ ности S ортогональны, их скалярное произведение равно нулю, и поэтому равен нулю интеграл по S. Если в интеграле по поверх­ ности S 2 изменить направление нормали, то придем к равенству

I l

anda= I l anda,

s,

s2

* О значении многоточий см. примечание в конце п. 184-

выражающему следующее положение: поток соленоидалъного век­ тора через поперечное сечение векторной трубки сохраняет по­ стоянную величину, ее называют интенсивностью векторной трубки.

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

Пусть источник находится в начале координат. Тогда

Е = 75-г,

где г = хі-\- у] + zk.

 

Рис. 144.

 

 

Проекции вектора Е будут Ех =

Е„

ЧУ

£ -

i l

Их

гз

 

r 3 •

производные

соответственно равны:

 

 

 

 

 

дЕг

7Г (7'2~ Зх2),

дЕи

ГЬ (г2- Ш

дЕг

2-

3z2).

 

дх

ду

 

dz

 

 

 

 

 

 

 

 

Поэтому, если г 0, то

дивергенция

вектора

напряженности

равна

 

дЕх

дЕи

дЕг

 

 

 

 

 

 

divE:

 

№ + У2 +*)]■

 

 

дх

ду

dz = -%■№>*S

 

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

divE(M) = 0.

(27)

Рассмотрим величину потока век­ тора напряженности через замкну­ тую поверхность S всторону ее внешней нормали. Возможны два случая. Если поверхность S не со­ держит внутри себя источника, то по формуле Остроградского, прини­ мая во внимание условие (27), полу­ чим

 

 

\ \ E nd a ^ 0.

(28)

 

 

s

 

Пусть S — сфера с центром в начале координат. Тогда

 

^ E nda ^

^ - ^ r . ^ d a = ^-2

d a = - ^ - 4 л й 2 = 4яд.

(29)

s

"s

s

 

Этим доказано, что поток вектора напряженности через по­ верхность сферы с центром в источнике в направлении внешней нормали не зависит от радиуса сферы и равен 4яд.

Докажем, что этот результат верен для любой гладкой замкну­ той поверхности S, содержащей заряд внутри себя. Для этого рассмотрим сферу 6Д с центром в начале координат, имеющую столь малый радиус, что она целиком содержится внутри S. Но формуле Остроградского для области В, ограниченной поверх­ ностями S и S lt имеем

JJ Ændtr + Jj

Enda = JJJ divEdx.

S

S t

B

Здесь правая часть равна нулю согласно формуле (27). Изменим направление нормали к сфере и в соответствии с формулой (29) получим

JJ Endo — J J Endo —4яд.

SS 1

199.Формула Стокса. * Формула Стокса связывает поверх­ ностный и криволинейный интегралы. Пусть поверхность S обла­ дает следующими свойствами: 1) это гладкая (или кусочно-глад­ кая) поверхность, ограниченная гладким (или кусочно-гладким) контуром I, 2) прямые, параллельные координатным осям, пере­ секают S не более чем в одной точке, 3) поверхность двухсторон­ няя; выберем ту сторону S, на которой выполняется условие

cosy = cos (n, z) > 0 .

Обозначим через X проекцию I на пло­

скость Оху, через А — проекцию S на плоскость Оху (рис. НБ).

Пусть поверхность

S задана уравнением z — f (х, у) в об­

ласти А. Направляющие косинусы нормали п найдем так же, как в п. 195:

 

cos а =

 

,

cos ß= —

 

 

I l

+p2 + q2

1

Kl + PH -g*

 

 

cos Y :

 

 

 

 

 

 

 

 

 

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

УI + р2+92

 

 

 

 

 

 

 

 

р cos Y =

— cos a,

 

g cos у — — cosß.

(30)

На поверхности S заданы непрерывно дифференцируемые

функции

Р (М ),

Q (М)

и R (М).

интеграл

в двойной,

а затем

Преобразуем

криволинейный

в поверхностный:

 

 

 

 

 

 

 

J Р (х, у, z) dx = J

P (x, y , / (x,

y)) dx =

 

 

i

 

1 i

 

 

и

 

= [ P{x, g, f(x, y)) dx = — I j (Py + P ’z -Zy)dxdy =

 

 

I

 

III

JAJ

IV

 

=

—■J j (Py + P'zq) cos y d a ^

j" j (P'zcos ßPÿ cos Y) da.

 

s

 

 

v

s

 

 

* Джон Габриель Стокс (1819—1903) — английский фішж и математик.

Здесь равенство I основано на том, что кривая I лежит на поверх­ ности S, и поэтому для всех точек I выполнено соотношение между координатами z = / (х, у). Переход II основан на том, что подын­

тегральная функция

зависит лишь от х и у, и поэтому согласно

правилу

вычисления

криволинейного

интеграла его

величина

не изменится при замене контура I

на К. Переход III основан на

формуле

Грина (18)

п. 190 в случае Q е_-.: 0. Переход IV основан

на формуле

 

z) cos у da = JJ

 

 

 

 

 

\ \ Р { х ,

у ,

Р (х,

у,

f(x, y))dxdy,

 

S

 

 

 

А

 

 

 

 

вытекающей из

соотношений (12)

и (15).

Переход V

выполнен

с помощью второй из формул (30).

 

 

 

 

Аналогично получим

 

 

 

 

 

 

[ Q dy -

J j (Qx cos у —Qz cos a) da,

 

 

i

 

 

s

 

 

 

 

 

 

J R dz = j j

(Ry cos a — R'x cos ß) da.

 

 

i

 

 

s

 

 

 

 

 

Сложив эти равенства,

придем к соотношению

 

 

 

 

j P dz -f Q dy + R dz =

 

 

 

 

 

I

 

 

 

 

 

 

= j j

\Шу — <?г) cos a + (p; — i?;)cosß + (0i — Py) cos Y}

da. (31)

s

 

 

 

 

 

 

 

 

 

Перейдем к поверхностному интегралу второго рода, положив согласно примечанию 1 п. 193 в правой части равенства (31)

cos a da ~ dy dz, cos $da ~dxdz, cos Y da ~ dx dy.

Получим формулу Стокса в координатной форме

j P dx-f Q dy + R dz = i

= J j {Ry — Q'z)dydz-ir (P ’z— R'x)dxdz^{Q 'x— Py)dxdy. (32) s

Пр и м е ч а н и е 1. Формула Грина есть тот частный случай формулы Стокса, когда S — область А плоскости Оху с конту­ ром I и Л ET 0.

Пр и м е ч а н и е 2. Формула (32) верна и в случае, когда прямые, параллельные координатным осям, пересекают S более

чем в одной точке, если выполнено следующее п р а в и л о о б ­ х о д а к о н т у р а : наблюдатель, обходящий I и направленный

по n *, должен иметь от себя слева часть поверхности S, непо­ средственно прилегающую к обходимому участку контура. Для доказательства этого утверждения достаточно разделить S вспо­ могательными контурами lk на части, в которых выполнено условие пересечения прямыми, параллельными координатным осям, не более чем в одной точке. Далее надо написать формулу Стокса для каждой части S и результаты сложить. При этом получится формула (32) для всей поверхности S в целом, потому что криво­ линейные интегралы по lk в сумме будут равны нулю, принимая во

внимание правило обхода контура.

а (М)

--

200.

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

P (М) і + Q (М ) j -f- R (М ) к. Левая часть формулы (32)

дает

 

величину работы этого векторного поля вдоль замкнутого кон­ тура /, т. е. циркуляцию вектора а (М) по контуру I. Согласно равенствам (8) и (11) п. 188 имеем

J P dx -J- Q dy -j- R dz =

J a dr

------J" а%ds.

i

 

i

i

О п р е д е л е н и е .

Ротором,

или

вихрем, вектора а (М )

называется вектор rot а (М), определяемый равенством

rot а (М) = (R'yQ'z) i -b (P'z R'x) \ + {Q'x Ру) k.

Его удобно представить символической формулой

 

і

3

к

rota(M )

д

д

д

дх

ду

dz

 

 

Р

Q

R

которую следует понимать в смысле равенства (33). Проекции вектора rot а по определению суть

(rota )x*=Ry — Q'z, (rota )y ^ P z’ — Rx,

(33)

(34)

 

(rota )z = Qx — Py.

(35)

Пусть

n = cos ai + cos ßj + cos yk есть нормаль

к поверх­

ности S,

соответствующая выбранной стороне поверхности. Тогда

подынтегральная функция правой части равенства (31) предста­ вляет скалярное произведение векторов n и rot а, а вся правая часть равенства (31) равна потоку вектора rot а через поверх­

ность S. Формулу Стокса (31)

можно представить

в следующей

в е к т о р н о й ф о р м е :

 

 

J a dr = j |

ro ta -n d o ,

(36)

is

*То есть n пронизывает наблюдателя от ног к голове.

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

Принимая во внимание формулу (И) п. 188, равенство (35) можно записать в виде

J ax(M)ds = JJ ro t„ a (M)da,

(37)

I

s

 

где ax (il/) — проекция

вектора a (M) на касательную

к кон­

туру I в точке М\ rot„ а — проекция rot а на нормаль п к поверх­

ности S в точке М.

 

 

о п р е д е ­

 

 

 

И н в а р и а н т н о е

 

 

 

л е н и е

вихря вектора

а получим с

 

 

 

помощью формулы Стокса (37). Для

 

 

 

этого фиксируем точку М и какое-либо

 

 

 

направление

и (начало п в точке М).

 

 

 

Пусть S — бесконечно

малая плоская

 

 

 

область

с площадью

о,

содержащая

 

 

 

точку М и перпендикулярная п. Напи­

 

 

 

шем формулу

(37)

и

по

теореме о

среднем для поверхностного

интеграла получим

 

 

 

 

 

 

 

 

 

 

 

J аТ (М ) ds = rot„ а (Ж*) о.

 

 

 

 

I

 

 

 

 

 

 

 

Отсюда следует равенство

 

I ах ds

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot„ а (Ж) =

lim

/-------,

(38)

 

 

 

 

 

 

(О - м

а

 

 

которое

позволяет

сформулировать

инвариантное определение

вихря вектора: вихрем вектора а (Ж) в точке М

называется век­

тор rot а (М), проекция

которого на любое направление п равна

пределу

отношения

циркуляции

вектора

а

вдоль

контура беско­

нечно малой

площадки,

содержащей

М

и

перпендикулярной п,

к величине этой площадки.

 

 

понятия вихря вектора

Ф и з и ч е с к о е

с о д е р ж а н и е

выясним с помощью его инвариантного определения. Для этого будем пока трактовать а как вектор скорости течения однородной несжимаемой жидкости. Поместим в точке М потока жйдкости бесконечно малое колесико с лопастями, расположенными по окружности перпендикулярно плоскости колесика (рис. 146). Под воздействием потока жидкости такое колесико будет, вообще говоря, вращаться со скоростью, зависящей от направления оси колесика п. Каждая точка окружности I вращается со скоростью ѵ, направленной по касательной к I. Величина этой скорости во

всех точках окружности одинакова и равна среднему значению касательной составляющей вектора а к линии I, т. е. среднему значению величины ах (М ) на окружности I, а именно

I

 

 

Теперь из формулы (38) следует, что

 

 

rotrtа (М)

2й>,

(39)

где (о — угловая скорость вращения

колесика. Если

векторы п

и rot а ортогональны, то rot„ а = 0 и

со = 0. Величина со дости­

гает наибольшего значения, когда векторы п и rot а коллинеарны, при этом

IW>

= - | ' l r o ta l'R и I rot а I = -2t?1j|-б- =

2сонаиб. -

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

наибольшая угловая скорость вращения равна

 

®наие=4“ ІГ 0 Іа І'

(40)

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

П р и м е р 1. Если жидкость течет с постоянной по величине и напра­ влению скоростью а, то будут постоянными и проекции этого вектора. По­ этому rot а = 0 и течение жидкости безвихревое.

П р il м с р 2- Если скорость а течения жидкости имеет проекции

 

Р = —ау,

Q сох,

R — 0,

(41)

где

ю — постоянная, то по формуле (33) получим

 

 

гоta = 2cok.

 

(42)

Из

(41) следует, что ѵ = | а | =

со У х 2 +

у2 = cor.

Поэтому вся жидкость

как целое вращается вокруг оси O z с постоянной угловой скоростью со. Дей­ ствительно, частицы жидкости, находящиеся на окружности с данным радиу­ сом г (в плоскости, перпендикулярной оси O z ) , движутся с постоянной скоро-

и

СТЬЮ V = ШГ II постоянной угловой скоростью (Û = — .

Каждая частица жидкости при движении вокруг оси O z участвует в двух движениях: в переносном движении со скоростью а (—соу, сох, 0) и во враща­

тельном движении, в котором мгновенная угловая скорость вращения каждой

1 частицы равна со, что следует из формул (40) и (42), так как ювр= — | 2сок|=

= ш. Следовательно, угловая скорость вращения каждой частицы совпадает с угловой скоростью макроскопического движения жидкости вокруг оси O z H равна со.

Ф о р м а л ь н ы е

с в о й с т в а

р о т о р а .

 

Пусть

а (М)

и b (М) — векторные

поля

и

cp (М) — скалярное

поле.

Тогда

1°. rot (а

b) •= rot а |

rot b,

 

 

 

 

 

 

 

2°. rot (cp а) фrot а -|-grad ф Xа.

 

 

 

 

 

 

Действительно, по формуле (33) имеем *

 

 

 

 

1) ro t(а .

Ь)

Г д (az-i-Ьг)

д(ау 'ГЬи)-\

,

 

 

 

 

L

 

ду

 

 

dz

J

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ daz

дау '

 

 

•••

, (

dbz

dby\ .

,

 

 

 

 

V ду

dz. .г

 

 

1 \

ду

dz

)

 

 

 

 

 

 

 

- д (ф(7г)

д (ФЙД 1 ; .

 

_

 

 

 

 

2) ro t(ф а) = [

 

ду

 

 

 

 

 

 

 

 

-

 

 

dz

J 1 ' ‘ '

 

 

 

 

 

-|- grad фX a.

 

 

 

 

+

 

 

!S )i + ---

''■prota

 

 

 

 

 

 

 

 

 

 

 

 

 

201.

Безвихревые

векторные

ноля.

Векторное поле а (М )

в области В называется безвихревым, если в каждой точке этой

области ротор вектора а (М) равен нулю:

 

 

 

 

 

 

 

 

 

 

rot а (М) ==0.

 

 

 

 

(43)

Следовательно, если поле безвихревое, то согласно (35) - вы­

полнено условие

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dÿ__ _дР_

 

 

 

 

дР _ дЕ

 

 

(44)

 

 

дх

 

 

ду

ду

dz

dz

дх

 

 

 

 

 

 

 

которое совпадает с условием потенциальности поля. Поэтому

всякое

безвихревое

поле

потенциально,

а

всякое

потенциальное

поле безвихревое.

В и. 192 установлено, что для потенциальности векторного поля а (М ) необходимо и достаточно, чтобы оно было полем гра­ диентов некоторой скалярной функции а (М) = grad и (М).

Следовательно, если выполнено это условие, то поле потен­ циально и по доказанному оно безвихревое, т. е. выполнено усло­ вие (43). Поэтому имеем

rotgrad и (М ) = 0,

(45)

если функция и дважды непрерывно дифференцируема.

Таким образом, доказано, что поле градиентов всегда безвих­ ревое.

Докажем, что поле вихрей соленоидально, т. е. имеет место

формула

(46)*

divrot а (М) = 0.

* О значении многоточий см. в и. 184. 22 Заказ Ц4

Д о к а з а т е л ь с т в о . Рассмотрим векторное поле а и поле его вихрен rot а. Составим выражение дивергенции поля вихрен. С помощью (35) получим

divrot а -^ (ro t аД-4~ (rot a)„-f

(rot a)t -■=

если P, Q и R дважды непрерывно дифференцируемы. Следова­ тельно, поле вихрей не имеет источников и стоков.

Поле а (М) называется гармоническим, пли лапласовым, если

оно безвихревое

и соленоидальное,

а

т. е.

rot а (М)

О

и div а (М) — 0.

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

grad

и (М), причем

потенциал этого поля и (М ) удовлетворяет уравнению Лапласа div grad и ~ 0 или Аи = 0. Например, гармоническим полем является электростатическое поле, создаваемое точечным источ­ ником.

§ 34. ДИФФЕРЕНЦИАЛЬНЫЕ ОПЕРАЦИИ ТЕОРИИ ПОЛЯ

ИИХ ВЫРАЖЕНИЕ В РАЗЛИЧНЫХ СИСТЕМАХ КООРДИНАТ

202.Дифференциальные операции второго порядка. Пусть даны дважды непрерывно дифференцируемые в области В функ­ ции ср ( М) и а (TIP). Выберем систему координат и вычислим вели­ чины

grad ф (TIP), div а (Ж), ro ta (TIP).

(1)

Действия, в результате выполнения которых получены вели­ чины (1), можно назвать дифференциальными операциями первого порядка. Они произведены над функциями ф (TIP) и а (TIP) соот­ ветственно.

Если принять величины (1) за исходные (это функции точки TIP), то путем выполнения над ними дифференциальных операций первого порядка (тем самым будут выполнены дифференциальные операции второго порядка над функциями ф (TIP) и а (TIP)) можно составить только следующие величины:

1. graddiva, 2. divgradф, 3. divrota, 4. гоідгагіф, 5. rotrot а

В п. 201 рассмотрены две из этих величин rot grad ф и div rot а, они тождественно равны нулю. Найдем для остальных величин их выражения через исходные данные.

1. По формулам (7) п. 184 и (21) п. 197 получаем

graddiv а = (div a)* i -f (div a)y j

(div a)2k;

grad div a (Pxx -■ Q x y ~ , B Xz)

i ;

(PXy "

Q y y B y z ) 3-f-

~Г (Px2 Q y z

P-'zz) k.

 

Оператор Лапласа* от скалярной функции ф(М) опреде­ ляется равенством

 

 

 

 

 

Аф =

фхж

у фу у -г фгг,

 

 

 

 

(3)

а от векторной

функции

а (М ) — равенством

 

 

 

 

 

 

 

 

 

Д а == аХх л

&уу~г ®z2-

 

 

 

(4)

Из

(4) следует,

что

если

а

P i

(J j

 

 

/Ік. то

 

 

 

 

 

Л а --

АР і

' з -j- Л/î k.

 

 

(5)

2. Из формул (7) п. 184 и (21) п. 197 следует равенство

divgrad ф--

(grad ф)*+-^- (grad ф)„ +

 

(grad ф)г =

 

 

 

дх

\

дх

 

 

ду

f ÜSL4) 4-

(

dq>\

 

 

 

 

 

 

\

ду

 

dz

\

dz

) ‘

 

Поэтому имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

divgrad ф=Аф.

 

 

 

 

 

(6)

3.

Исходя из определения

ротора, последовательно получаем

 

 

і

 

3

 

 

к

 

 

 

 

 

 

 

 

 

rotrot а

д

 

д

 

 

д

 

 

l(Qx-Py)y-(Pz-RxYZ]

дх

 

ду

 

 

dz

 

 

rot* а

 

rot,, а

гоС а

 

 

 

 

 

 

 

 

 

 

=-■КК х

-1- Q "xy -і- R 'x z ) —

( Р х х Л - Р ' у у

+

P "zz)] І + ... .

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

 

rotrot а = graddiv а — А а.

 

 

(7)

 

 

 

 

 

 

С и м в о л и к а

 

Г а м и л ь т о н а

**.

 

Введем

символйче-

ский вектор

(называемый набла)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д

i

.

9

i ,

 

д

 

 

 

/оч

и условимся

в

следующем:

1) будем

понимать

под

«произведе-

нием» символов

д

д

 

д

 

 

 

 

 

 

 

 

 

 

 

 

^

и — на скаляр ф соответственно величины

частных производных

да>

9 ф

и

9 ф

,

 

 

 

 

это

«произведение»

 

 

~

г~

(конечно,

не обладает переместительным свойством),

2) алгебраические дей­

ствия

с вектором набла

будем выполнять

по правилам действий

*Пьер Симон Лаплас (1749—1827) — французский астроном, математик

ифизик.

**Уильям Роуан Гамильтон (1805—1865) — английский математик.

22*

339

Соседние файлы в папке книги из ГПНТБ