Математический анализ II Учебное Пособие
.pdfQ = ∫∫ax (x ,y ,z )dy dz +ay (x ,y ,z )dzdx +az (x ,y ,z )dxdy .
S
Вспоминая, что под знаком тройного интеграла в правой части стоит дивергенция векторного поля a , тогда окончательно можно написать
∫∫andS = ∫∫∫div adxdydz .
S T
Выясним теперь механический смысл дивергенции векторного поля. Допустим, что векторное поле есть поле скоростей текущей жидкости
v(x ,y ,z ,t ) и пусть плотность этой жидкости ρ =const . Возьмём в поле
вектора v замкнутую поверхность S , ограничивающую малый объём. По теореме Остроградского
QT = ∫∫∫div vdxdydz .
T
Применяя теорему о среднем, получим QT = (div v)M vT , где M - некоторая "средняя точка", лежащая в телеT . Отсюда
(div v) |
|
=Q . |
|
M |
T |
|
|
|
Сжимая тело в точку, в пределе проучим |
||
(div v) |
= limQT . |
|
M |
|
λ→0 vT |
|
|
Вывод: дивергенция векторного поля v даёт нам расход жидкости из точечного источника в единицу времени, т.е. удельную силу источника. Заметим, что если div v в данной точке положительна, то это значит, что в этой
точке происходит втекание жидкости (т.е. в этой точке сток).
Принимая во внимание это рассуждение, можно дать механическое истолкование теоремы Остроградского:
расход жидкости из тела T , ограниченного поверхностью S , в единицу времени t равен сумме попарных произведений сит источников.
Заметим, что попутно мы доказали независимость дивергенции векторного поля от выбора системы координат. Действительно, с одной стороны
div a = ∂∂axx + ∂∂ayy + ∂∂azz ,
а с другой - это удельная сила источника, а независимость последней от системы координат очевидна.
114
3. Q 3 = ∫∫[x cos λ +(y +1)cos μ +z cosν]dS , n3 (0,0, −1) , т.е.
S 3
|
|
|
|
|
|
|
cos λ = 0 , cos μ , cosν = −1, |
|
|
|
|||||||||||||||||
Q = −∫∫zdS , т.к. на поверхности S 3 |
|
z = 0 . |
|
|
|
|
|
|
|
||||||||||||||||||
S 3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4. Q 4 = ∫∫[x cos λ +(y +1)cos μ +z cosν]dS . |
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
S 4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Поверхность S 4 имеет уравнение z =1 −x −y , следовательно, |
|
||||||||||||||||||||||||||
|
|
|
p (x ,y ) = |
∂z (x ,y ) |
|
= −1, q(x ,y ) = |
∂z (x ,y ) |
= −1, |
|
|
|||||||||||||||||
|
|
|
|
|
∂x |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|||||
тогда dS = |
p 2 (x ,y ) +q 2 (x ,y ) +1 = |
3 . Поверхностный интеграл здесь вы- |
|||||||||||||||||||||||||
числяется по верхней стороне поверхности S 4 , значит направляющие коси- |
|||||||||||||||||||||||||||
нусы нормали n4 будут равны: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
cos λ = |
1 |
, cos μ = |
1 |
, cosν = |
1 |
. |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
3 |
3 |
|
|
|
|||||||||||||||
Тогда получим |
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
1 |
|
1 |
|
|
1 |
|
|
1 |
|
|
|
|
1 |
|
|
|
|
1 |
|
|||||||
Q4 =∫∫ x |
|
|
+(y +1) |
|
+z |
|
|
|
dS = |
∫∫ x |
|
|
|
|
+(y +1) |
|
|
+(1−x −y ) |
|
3dxdy = |
|||||||
|
3 |
3 |
|
|
|
|
3 |
3 |
3 |
||||||||||||||||||
S4 |
|
|
|
3 |
S3 |
|
|
|
|
|
|
|
|
||||||||||||||
= 2∫∫dxdy =1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
S 3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Окончательно: Q =Q1 +Q 2 +Q 3 +Q 4 = 0 − |
1 |
+ 0 +1 = |
1 . |
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
2 |
|
|
|
|
3-й метод решения. Вычислим тот же самый поток с помощью поверхностного интеграла второго рода Q = ∫∫axdydz +aydzdx +a zdxdy .
|
|
|
S |
|
|
В нашем случае |
Q = ∫∫xdydz +(y +1)dzdx +zdxdy , как и в предыду- |
||||
|
|
S |
|
|
|
щем случае, поток Q представим в виде суммы четырёх потоков соответст- |
|||||
венно, через поверхности S1 , S 2 , |
S 3 , S 4 : |
|
|
||
1. На S1 |
x = 0 , dx = 0 , а значит Q1 = 0 . |
|
|
||
2. На S 2 |
y = 0 , dy = 0 , а сторона поверхности, по которой вычисляется |
||||
интеграл, нижняя, значит |
|
1−∫x dz = − |
|
||
|
|
Q 2 = −∫∫dzdx = −∫1 dx |
1 |
||
|
|
S 2 |
0 |
0 |
2 |
|
|
|
|||
3. На S 3 |
z = 0 , dz = 0 , сторона поверхности нижняя значит Q 3 = 0 . |
||||
4. Q4 = |
∫∫ |
xdydz +(y +1)dzdx +zdxdy = |
|
S4 (верхн. стор)
116
′ |
на противоположное, тогда получим |
Изменим направление нормали n |
|
∫∫andS = ∫∫andS |
|
S1 |
S 2 |
Следовательно, мы можем сделать вывод, что поток соленоидального векторного поля через поперечное сечение векторной трубки есть величина постоянная и не зависит от площади сечения. Эта величина называется напряжением векторной трубки.
Если соленоидальное поле является полем скоростей текущей несжимаемой жидкости v(vx ,vy ,vz ) , то мы имеем
∂v |
x |
(x ,y ,z ) |
+ |
∂vy (x ,y ,z ) |
+ |
∂v |
z |
(x ,y ,z ) |
= 0 |
(1) |
|
|
∂x |
∂y |
|
|
∂z |
||||
|
|
|
|
|
|
|
|
Значит, если в пространстве, где течёт жидкость, нет ни источников, ни стоков, то проекции скоростей связаны соотношением (1), которое называется уравнением неразрывности.
Рассмотрим теперь некоторое скалярное поле U (x ,y ,z ) и найдём его градиент:
a =gr adU (x ,y |
,z ) = |
∂U i + |
∂U |
j + |
∂U k . |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
∂x |
|
∂y |
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|||
Найдём далее дивергенцию векторного поля a : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
div a =divgr adU |
= U |
= ∂2U + |
∂2U |
|
+ ∂2U . |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
∂x 2 |
|
∂y 2 |
|
|
|
∂z 2 |
|
||||||||||
Определение 2. |
|
|
|
|
|
|
r |
|
r |
r |
|
|
∂2 |
|
|
|
|
∂2 |
|
|
|
|
∂2 |
|
||||
Оператор |
|
2 |
= = |
|
|
|
+ |
|
|
|
|
|
|
+ |
|
|
|
|
называется |
|||||||||
|
∂x 2 |
|
∂y |
2 |
|
|
∂z 2 |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
оператором Лапласа и обозначается символом |
, т.е. |
|
|
|
|
|
|
|
||||||||||||||||||||
r2 |
|
∂2 |
|
∂2 |
|
|
∂2 |
, а U = |
∂2U |
|
+ |
∂2U |
+ |
∂2U |
||||||||||||||
= |
= |
|
+ |
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
∂x 2 |
∂y 2 |
|
∂z 2 |
∂x |
2 |
|
∂y 2 |
∂z |
2 |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
называется Лапласианом функцииU (x ,y ,z ) . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Определение 3. Если функцияU (x ,y ,z ) такова, что |
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
∂2U |
+ ∂2U |
+ ∂2U |
= 0 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
∂x 2 |
|
|
∂y 2 |
|
∂z 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
т.е. U = 0 , то функция называется гармонической.
Теорема. Для того, чтобы поле градиента какой-нибудь скалярной функции было соленоидальным, необходимо и достаточно, чтобы эта функция была бы гармонической.
Доказательство. Необходимость. Поле градиента соленоидально, |
|||
ur uur |
r uur |
= U . Значит, U |
= 0 , а это означает, что функ- |
т.е. U |
= 0 , но U |
цияU (x ,y ,z ) гармоническая.
Достаточность. Пусть функцияU (x ,y ,z ) - гармоническая, т.е.
118
|
|
Докажем, однако, эту теорему подробнее. |
Обозначим |
|
|
r |
|
тогда |
||||||||||||||||||||||||||||||||||||||||||||||
|
|
c = ×a , |
||||||||||||||||||||||||||||||||||||||||||||||||||||
r |
|
|
|
∂cx |
|
∂cy |
|
|
|
|
|
∂cz |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
|
∂a z |
|
∂ay |
|
||||||||
c |
= |
+ |
+ |
, |
|
|
|
|
|
но |
|
|
|
|
координаты |
|
|
|
вихря |
x |
= |
− |
, |
|||||||||||||||||||||||||||||||
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
∂z |
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
∂x |
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|||||||||||
c |
y |
= |
∂ax |
− |
∂a z , c |
y |
= |
|
∂ay |
− |
∂ax . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
∂z |
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ |
|
|
|
∂ ∂a |
|
∂ |
|
|
|
|
|
|
|
||||||||||||||||||||||
r |
|
|
|
∂ |
∂ |
|
|
|
∂a ∂ ∂ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
a z − |
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
ax |
|
|
|
|
a z |
+ |
|
|
|
|
|
|
|
|
y |
− ax = |
|
|
|
|
|
|
|
|||||||||
c |
= |
|
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
∂z |
|
|
|
∂y |
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
∂x |
∂y |
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
∂x |
|
|
∂z |
|
|
∂y |
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
= |
|
∂2a |
z |
|
− |
∂2ay |
|
|
+ |
|
∂2a |
x |
− |
|
∂2a |
z |
|
+ |
|
∂2ay |
− |
|
∂2a |
x |
|
= 0 . |
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
∂z ∂x |
|
|
∂z ∂y |
|
∂x |
|
|
|
|
∂x |
∂y |
|
∂y |
∂z |
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
∂y ∂z |
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Теорема доказана.
Заметим, что при доказательстве теоремы мы предположили существование вторых смешанных производных от проекций вектора a и их непрерывность.
Допустим, что в поле вектора a лежит замкнутый контур l , имеющий в каждой своей точке касательную. Перекинем через контур l двухстороннюю поверхность S , имеющую в каждой своей точке нормаль n к поверхности. При таких предположениях о контуре l и поверхности S справедлива доказанная ранее формула Стокса.
Докажем теперь векторную форму теоремы Стокса.
Теорема Стокса. (векторная форма).
Циркуляция векторного поля по замкнутому контуру равна потоку его вихря через поверхность, натянутую на этот контур.
Доказательство. Итак, пусть через замкнутый контур, обладающий указанными выше свойствами, перекинута (натянута) двухсторонняя поверхность S . Рассмотрим циркуляцию векторного поля по замкнутому кон-
туру l : Ц= ∫aτdS .
l
Пусть касательный вектор τ образует углы α , β и γ с координатными
осями O x , O y и O z , тогда aτ =прτra = a τ0 =ax cosα +ay cos β +a z cosγ . Подставим aτ в выражение циркуляции, тогда получим
Ц= ∫ ax cosα +ay cos β +a z cosγ dS или, принимая во внимание связь меж-
l
ду криволинейными интегралами первого и второго рода, получим
Ц= ∫axdx +aydy +a zdz .
l
Воспользуемся теперь формулой Стокса:
|
|
∂a |
|
∂ay |
|
∂a |
|
∂a |
|
∂ay |
|
∂a |
|
||
∫axdx +aydy +azdz =∫∫ |
|
z − |
|
dydz + |
|
x − |
|
z dzdx + |
|
− |
|
x dxdy . |
|||
|
∂z |
|
|
∂x |
|
||||||||||
l |
S |
∂y |
|
|
∂z |
∂x |
|
|
∂y |
||||||
|
|
|
|
|
|
120 |
|
|
|
|
|
|
|
|
|
При этом предполагается, что наблюдатель, у которого нормаль к поверхности проходит от ног к голове, обходит контур в таком направлении, что поверхность S остаётся слева.
Нетрудно заметить, что в правой части последней формулы находится поток вихря векторного поля a через поверхность S , т.е. имеем
∫aτdS = ∫∫(r ot a)n dS .
l S
Теорема доказана.
Теорема 2. Вихрь векторного поля не зависит от выбора координатных осей.
Доказательство. Рассмотрим в поле вектора a произвольное направ-
ление, заданное некоторым вектором n . |
|
|
|
|
|
|||||||
z |
|
|
|
|
|
|
Проведём окружность l |
радиуса r , |
||||
C |
1 |
|
|
плоскость круга S , ограниченного кон- |
||||||||
|
l 2 : y +z =1 |
|
туром |
l , перпендикулярна |
вектору n |
|||||||
l 3 |
|
(рис. 9). По теореме Стокса |
|
|||||||||
|
|
|
|
|||||||||
0 |
B |
y |
|
|
|
|
∫aτdS = ∫∫cndS , |
|
||||
|
|
|
|
l |
|
|
S |
|
||||
x 1A |
1 |
|
|
где cn =прn rot a , т.е. c =rot a . |
||||||||
l1 :x 2 +y 2 =1 |
|
Тогда по теореме о среднем имеем |
||||||||||
|
рис. 9 |
|
|
∫aτdS =[cn ]M |
πr 2 , |
|
||||||
|
|
|
|
l |
|
|
|
|
|
|
|
|
где M - некоторая "средняя точка", лежащая в круге S . Тогда |
|
|||||||||||
|
|
[cn ]M |
= |
1 |
∫l aτdS . |
|
||||||
|
|
πr 2 |
|
|||||||||
Отсюда следует, что c |
= lim |
1 |
|
∫l |
a |
dS , где M - центр круга S . Выра- |
||||||
|
|
|
||||||||||
|
[ n ]M |
r →0 |
πr 2 |
τ |
|
|
|
|
|
жение в правой части равенства не зависит от выбора координатных осей. Значит cn тоже от них зависит, а так как n - любое направление, то вектор
r ot a не зависит от выбора координатных осей.
Пример. Вычислить циркуляцию векторного поля a =xz i +x j +y k по
|
|
линии пересечения конуса |
|
n |
|
x 2 +y 2 −(z −1)2 = 0 |
|
M |
τ |
с координатными плоскостями, лежа- |
|
l |
щей в первом октанте, непосредственно |
||
|
|||
|
и по теореме Стокса (рис. 10). |
||
рис. 10 |
|
||
|
|
||
|
|
121 |
Решение.
1) Непосредственное вычисление циркуляции. Контур l можно разбить на три части: l 1 , l 2 и l 3 , лежащие в координатных плоскостях xO y , yO z и
zOx соответственно, таким образом циркуляция Ц=Ц1+Ц2+Ц3 , где
Ц1= ∫xzdx +xdy +ydz . На кривой l1 : z = 0 , dz = 0 x = 1 −y 2 , y [0,1].
l1 |
|
|
|
Следовательно, Ц1= ∫1 |
1 −y 2dy = |
π |
. Далее |
0 |
|
4 |
|
Ц2= ∫xzdx +xdy +ydz . На кривой l2 : y +z =1, x = 0, dx = 0, z [0,1], т.е.
|
|
l2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ц2= ∫1 |
(1 −z )dz = |
1 . И, наконец, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
0 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ц3= ∫xzdx +xdy +ydz . |
На кривой l3 : |
x +z =1, |
y = 0 , |
dy = 0 , x [0,1], |
||||||||||||||||||
|
|
l3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
следовательно, Ц3= ∫1 x (1 −x )dx = |
1 . |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
0 |
|
|
|
|
6 |
+ 1 |
+ 1 |
= π + |
2 . |
|
|
|
|
|||
Окончательно Ц=Ц1+Ц2+Ц3 = |
π |
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
4 |
2 |
6 |
6 |
3 |
|
|
|
|
||
2) Вычисление циркуляции по теореме Стокса: |
|
|
|
|
||||||||||||||||||
Ц= |
∫∫S |
∂a |
|
∂a y |
∂a |
|
|
∂a |
|
|
|
∂a y |
|
∂a |
|
|
||||||
|
∂y |
∂z |
|
|
∂z |
|
∂x |
|
|
|
∂x |
|
∂y |
|
||||||||
|
|
|
z − |
|
dydz + |
|
|
x |
− |
|
z dzdx + |
|
|
− |
|
x |
dxdy . |
|||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||
Подставим сюда az =xz , ay |
=x , az =y , получим |
|
|
|
|
Ц1= ∫∫dy dz +xdzdx +dzdy . Перейдём в правой части к поверхностному
S |
|
|
|
∫∫[ |
|
|
|
|
|
|
|
|
|
|
] |
|
|
|
|
|
|
|||||
интегралу первого рода Ц= |
|
|
|
|
|
|
|
|
|
|
|
, |
где |
интеграл |
||||||||||||
|
1 cos |
λ +x cos μ +1 cosν |
dS |
|||||||||||||||||||||||
|
|
|
|
|
|
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
вычисляется по верхней стороне поверхности S . |
|
|
|
|
|
|
|
|
|
|||||||||||||||||
Уравнение поверхности S : z =1 − |
|
x 2 +y 2 , |
|
следовательно, |
|
|
|
|||||||||||||||||||
|
p (x ,y ) = |
∂z |
= |
|
−x |
|
|
, q(x ,y ) |
= |
∂z |
= |
−y |
|
|
, |
|
|
|||||||||
∂x |
x 2 +y 2 |
|
∂y |
x 2 +y 2 |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
dS = 1 + p 2 (x ,y ) +q 2 (x ,y ) = 2 ; |
|
|
|
|
|
|
|
||||||||||||||||||
cos λ = |
−p (x ,y ) |
|
= |
|
|
x |
|
; |
cos μ = |
|
|
−q(x ,y ) |
|
= |
|
|
y |
|
; |
|||||||
1 + p 2 +q 2 |
|
2 x 2 +y 2 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
1 + p 2 +q 2 |
|
|
2 x 2 +y 2 |
||||||||||||||
|
|
|
cosν = |
|
1 |
|
= |
|
1 |
. |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
1 |
+ p 2 +q 2 |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
122 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|