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

Chast_8_Teoria_polya

.pdf
Скачиваний:
20
Добавлен:
23.02.2015
Размер:
2.2 Mб
Скачать

Таким образом, по данной формуле поток сводится к

 

интегралу 1-го рода по поверхности от скалярного

0

произведения вектора a P на нормаль n0 (P ) к этой

 

поверхности (иначе: от проекции поля a P

на

0

 

 

нормаль к поверхности n0

(P )).

xy

 

 

 

3.3.1. Проектирование на одну координатную плоскость

Пусть поверхность задана явно уравнением

z f (x, y)

и однозначно

проектируется в область Dxy на координатной плоскости Oxy .

 

 

 

 

 

 

 

 

 

 

 

fx , fy , 1 ,

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

n

n0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fx

 

 

 

 

 

 

 

 

 

 

 

 

 

fy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

fx

 

fy

 

 

fx

 

fy

 

 

 

fx

 

 

fy

 

1

 

 

и поток вектора a a P a x, y, z

через эту поверхность равен

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dxdy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

fy

2

 

 

 

 

a,n0 d

a,n0

 

 

 

 

 

 

 

a,n0

 

 

fx

 

1

dxdy

 

 

 

cos

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ax x, y, f x, y fx ay x, y, f x, y fy az x, y, f x, y dxdy ,

Dxy

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

Аналогичные формулы получаются при проектировании на другие координат-

ные плоскости для поверхностей вида x f ( y, z)

и y f (x, z) .

 

 

 

 

 

 

 

3.3.2. Проектирование на три координатные плоскости

 

 

 

 

 

 

Пусть поверхность задана (неявно) уравнением F x, y, z 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n

 

 

 

 

 

 

 

 

 

, n0

 

 

 

 

=

 

 

 

 

 

 

 

 

Fx

,

Fy

, Fz ,

 

 

Fx

 

Fy

Fz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

2

 

 

 

 

2

 

 

 

2

 

 

2

 

 

2

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fx

 

Fy

Fz

 

 

 

 

 

 

Fx

 

Fy

 

Fz

 

 

 

 

 

 

Fx

 

Fy

 

Fz

 

 

 

20

Пусть , , - углы, которые образует нормаль с осями координат. Тогда

орт n0 имеет координаты: n0 (cos ,cos ,cos ) . Так как a ax ,ay ,az , то

a,n0 ax cos ay cos az cos

и

a,n0 d = ax cos d ay cos d az cos d .

Рассмотрим отдельные слагаемые: az cos d . Если поверхность описыва-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ется уравнением z z(x, y) ,

 

 

а поле a P

в поверхностном интеграле берётся в

точке P , для любой его компоненты координата z

выражается через x и

y ,

az x, y, z az x, y, z x, y ,

d

 

 

dxdy

 

 

,

и

az cos d =

 

 

cos

 

 

 

 

 

 

 

 

dxdy

 

 

 

 

 

 

 

 

 

 

 

az

(x, y, z(x, y))cos

 

 

 

 

az (x, y, z(x, y))dxdy .

 

 

 

 

cos

 

 

 

 

 

 

 

 

Dxy

 

 

 

 

 

Dxy

 

 

 

 

 

 

 

 

 

Знак (+) соответствует острому углу между нормалью и осью z (cosγ > 0), знак

(–) – тупому углу между нормалью и осью z (cosγ < 0). Аналогично,

ax cos d ax (x,(y, z), y, z)dydz ,

 

 

 

Dyz

 

 

 

ay

cos d ay (x, y(x, z), z)dxdz ,

 

 

 

Dxz

 

и окончательно имеем:

 

 

 

 

(x,(y, z), y, z)dydz ay (x, y(x, z), z)dxdz az (x, y, z(x, y))dxdy .

(a,d ) ax

 

Dyz

 

Dxz

Dxy

1). Знаки перед слагаемыми соответствуют знакам направляющих косинусов нормали cos ,cos ,cos .

2). Вычисление потока векторного поля сводится к вычислению трёх двойных интегралов при условии, что поверхность взаимно однозначно проектируется на все три координатные плоскости. Если это не имеет места, поверхность нужно разбить на однозначно проектирующиеся участки.

3). Указанная формула устанавливает связь между потоком и поверхностным интегралом 2-го рода

(a,d ) (ax (x, y, z)dydz ay (x, y, z)dxdz + az (x, y, z)dxdy) .

 

 

21

3.4. Физический смысл потока

Пусть a P - поле скоростей некоторой жидко-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сти a V , а - произвольная поверхность в поле, то-

 

 

 

 

 

 

гда: a,d ,n d =

 

V

 

cos d = Ï ð V d - объём

 

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

n

 

 

 

Pi

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

столба жидкости с основанием d

и высотой Ï ðn0V ,

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е. объем жидкости, протекающей через площадку d

 

 

 

 

 

 

в единицу времени в направлении n0 . Суммируя по по-

 

 

 

 

 

 

 

 

 

 

 

 

-

 

 

 

 

 

 

 

 

 

верхности , получаем, что a d

поток жидкости, протекающей через

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

поверхность в единицу времени.

 

 

 

 

 

 

 

 

 

 

 

 

ПРИМЕР. Вычислить поток векторного

 

поля радиус-

 

 

 

 

 

 

вектора a r(x, y, z) через внешнюю сторону цилиндра

 

 

Z

0

 

3

 

 

 

 

(H– высота, R- радиус).

 

 

 

 

 

 

 

 

 

 

 

0

2

 

 

 

 

 

 

 

 

 

 

 

Решение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

a(r) r ;, 1 2 3

 

 

 

 

 

 

 

X

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2 3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П1 a,d = r,n0 d =…

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

на нормаль к 1

равна R}

 

 

 

{ r,n0 R , из рисунка ясно, что проекция r

 

 

 

 

 

…= Rd R d 2 R2 H .

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П2 r

,d = r,n0 d = …

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

{из рисунка ясно, что проекция r

на n0

по 2

равна H,

 

 

 

 

 

 

т.е. a, n0 2 H }

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

…= Hd H d R2 H .

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П3 a, d = r, n0 d =0.

 

 

 

 

 

 

 

 

 

 

2

 

 

 

3

 

 

3 R2H.

 

 

 

 

 

 

 

 

 

 

П П П

2

П

3

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

ПРИМЕР. Вычислить поток векторного поля a(r) y2 j zk через всю поверх-

z x2 y2

ность (нормаль внешняя): : .

z 2.

22

Решение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2

 

и пред-

 

 

 

Z 0 2

П2

Разобьем поверхность на две части

 

 

 

Z=2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ставим поток в виде П П1 П2 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1 П1

П1

 

 

 

 

 

 

 

 

 

 

 

 

 

a(r) (0, y

2

, z) ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a,d

=

a,n0 d ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n (2x,2y, 1) ,

 

n

 

 

 

 

 

 

4x

2

4y

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,

n0

 

 

 

 

(знак вы-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

бирается «+», так как cos( ) 0 ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2y

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

n0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4y

2 1

 

 

 

4x2 4y2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4x2

 

 

 

 

 

 

4x2 4y2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П1

 

 

2y3 z

 

 

 

 

 

 

 

 

 

 

 

 

 

2y3 z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

4x

2

4y

2

1dxdy =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4x

2

 

4y

2

 

 

 

 

 

 

4x

2

4y

2

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

Dxy

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2y3 (x2 y2 ))dxdy =…

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dxy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{перейдем в полярную систему координат}

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

... dφ d (2 3 sin3 φ 2 ) ... 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a,d =

a,n0 d ...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{ n0 (0;0;1) (a,n0 ) z }

 

 

 

 

 

 

 

 

 

 

 

 

 

... zdxdy 2 dxdy =2 (

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)2 4 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dxy

 

 

 

 

 

 

 

 

 

 

Dxy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П П1 П2

2 4 2 .

 

 

 

 

 

 

ПРИМЕР. Найдите поток вектора a xyi yz j xzk через часть сферы

 

x2 y2 z2 1, расположенную в первом октанте (нормаль внешняя).

 

Решение:

П ax dydz ay dxdz az dxdy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{компоненты поля и области интегрирования обладают

 

 

 

 

 

 

симметрией относительно замены x y z и

 

 

 

 

 

 

 

 

 

 

 

 

 

Dxy

 

D yz Dxy Dxz }

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П 3 azdxdy 3 x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x2

y2 dxdy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dxy

 

 

 

Dxy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ 2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 d cos

1

2

d

=

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23

Важно отметить, что cosα, cosβ, cosγ положительны, перед всеми интегралами берется знак (+), так как сторона поверхности - внешняя.

3.5. Дивергенция векторного поля Дивергенция - это дифференциальная и локальная (зависит от точки) ко-

личественная характеристика векторного поля. Пусть вектор-функция

a(P) ax i ay j az k имеет непрерывные частные производные первого поряд-

ка по всем переменным.

ОПРЕДЕЛЕНИЕ. Дивергенцией векторного поля a a(Ð)в точке Р(x,y,z)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ax

(x, y, z)

 

 

ay (x, y, z)

az (x, y, z)

называется число

 

 

diva

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, или, опуская ар-

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ay

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ax

 

 

 

 

 

 

az

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

гументы:

 

diva

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

Используя оператор Гамильтона (набла):

 

 

 

x

 

y

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

j

 

k

 

 

,

 

дивергенцию можно записать в виде скалярного произве-

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

y

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дения diva ( ,a) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.5.1. Свойства дивергенции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1. Линейность div( a b) diva divb , где и - произвольные постоян-

 

 

ные.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Пусть u u(x, y, z)

- скалярное поле, тогда div(u a) u diva (a gradu).

Доказательство:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div(u a)

(u a

 

)

 

 

(u ay )

 

 

(u a

)

 

=

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

u

 

ay u

 

 

 

a

 

u

 

 

 

 

a

x

 

 

 

 

ay

 

 

 

 

a

 

)= agradu u diva .

=

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

+u (

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

x

 

 

 

 

 

 

y

 

 

 

 

z

 

 

 

 

 

 

 

x

 

 

 

 

 

y

 

 

 

 

 

z

 

 

ПРИМЕР. 1). a

 

 

 

 

 

xi yj zk .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z 1 1 1 3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

divr

 

 

 

 

 

 

x

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2). a (c ,c

 

,c

) ,

 

diva

 

 

c

 

 

c

 

 

 

 

c

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 1

 

 

 

y

2

 

 

 

z

3

 

 

 

 

 

 

 

 

24

 

 

3.5.2. Физический смысл дивергенции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(a,d )

 

 

 

 

 

Поскольку величина

 

 

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

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a,d

 

 

 

пространственной области G, то

lim

 

 

diva

есть плотность потока в

V

 

точке М.

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Точки

поля,

в

которых

дивергенция

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

т.е.

diva(M ) 0 0 , называют источниками векторного поля, а точки, в которых дивергенция отрицательна, div a(M ) 0 0 - стоками векторного поля.

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

ОПРЕДЕЛЕНИЕ. Величину diva(M ) называют мощностью источника или

стока.

3.6. Физический смысл потока через замкнутую поверхность

Рассмотрим замкнутую поверхность , ограничи-

a

вающую объем G в векторном поле a a(P) скоростей

0

течения несжимаемой жидкости.

Поток вектора a a(P) через поверхность

(a,) равен количеству жидкости, протекаю-

0

P

2

a

P 0

1

щему через поверхность в единицу времени. Обозна-

0

чим единичный вектор внешней нормали n0 . Вектор-

 

ные линии входят и выходят из замкнутой поверхности . В точке P1 угол

 

 

 

 

 

 

 

 

 

 

 

(a,n0 )

 

 

; это означает, что жидкость втекает внутрь поверхности.

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В точке выхода P2 (a,n0 )

 

 

, жидкость вытекает. Поток векторного поля

2

 

 

 

 

 

 

 

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

Если П>0, жидкости вытекает больше, чем втекает, следовательно, в области G есть источники поля.

Если П<0, втекает жидкости больше, чем вытекает, значит в G есть стоки.

25

Если П=0, это означает, что источников и стоков нет или они компенсируют друг друга.

3.7. Теорема Остроградского - Гаусса

Если в некоторой области G трёхмерного пространства, ограниченной замкнутой кусочно-гладкой поверхностью , задано непрерывно дифференци-

руемое векторное поле a axi ay j azk , то поток

векторного поля

a

 

через внешнюю сторону замкну-

 

той поверхности

 

равен

тройному интегралу от

X

a

x

 

ay

 

a

z

 

по области G, ограничен-

 

функции

 

 

 

 

 

 

 

 

y

z

 

 

x

 

 

 

 

 

 

 

 

 

 

 

ной поверхностью :

 

 

 

 

 

 

 

ay

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ax

 

 

az

 

 

 

 

 

 

(a d )

 

 

 

 

 

dxdydz ,

 

 

 

 

 

y

z

 

 

 

 

 

 

 

 

 

 

V

x

 

 

 

 

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

Доказательство:

Z

 

2

 

3 V

 

1

 

Y

xy

Г

 

а) Рассмотрим область G, правильную в направлении оси Oz, которую будем называть элементарной Hz областью. Это означает, что снизу и сверху она ограничена поверхностями: 1 : z1 z(x, y) и 2 : z2 z(x, y) соответственно, а сбоку цилиндрической поверхностью 3 с образующими, параллельными оси Oz и направляющей Г.

Рассмотрим одно слагаемое:

 

a

z21

a

 

z2 (x, y)

 

z

dxdydz dxdy

z

dz

dxdyaz (x, y, z)

z (x, y)

z

z

G

 

D

z

 

 

D

1

 

 

xy

1

 

 

xy

 

= dxdy az (x, y, z2 (x, y)) az (x, y, z1(x, y) =

Dxy

= az (x, y, z2 (x, y))dxdy az (x, y, z1(x, y))dxdy =…

Dxy

Dxy

 

{на 2

cos 0 , а на 1

cos 0.

Учитывая, что dxdy cos d , получаем:

на 2 : dxdy cos d ,

на 1 : dxdy cos d } …= az (x, y, z)cos d + az (x, y, z)cos d =…

2 1

26

Добавим интеграл по 3

az (x, y,z)cos d в полученную сумму,

так как на 3

 

 

3

 

 

cos всюду равен нулю,

а следовательно, и az (x, y, z)cos d 0.

Тогда

 

 

3

 

 

 

 

 

…= az (x, y, z)cos d + az (x, y, z)cos d + az (x, y, z)cos d =

2

 

1

3

 

 

 

az (x, y, z)cos d = az (x, y, z)cos d .

 

 

1 2 *

 

 

 

б) Рассмотрим пространственную область G, которую можно разбить на n

 

 

 

 

 

 

 

 

n

 

 

 

элементарных областей Hz типа, т.е.

G Gk . Докажем, что и в этом случае

 

 

 

 

 

 

 

 

k 1

 

 

 

справедлива теорема Остроградского-Гаусса.

 

 

 

Пусть k

,

k ,

k

- нижняя, верхняя и боковая части поверхности k ,

1

 

2

3

 

 

 

 

 

 

 

 

ограничивающей область Gk ,

 

 

 

 

 

 

тогда

 

 

 

az

 

 

n

 

az

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dxdydz

 

dxdydz

 

 

 

 

 

 

 

 

G

z

 

k 1

G

z

 

 

 

 

 

 

 

 

 

k

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= az cos d ,

 

az

cos d + az cos d + az cos d

k 1

 

k

 

 

 

k

 

k

 

 

 

 

 

1

 

 

 

2

 

3

 

 

 

так как интегралы по k

равны нулю, а по поверхности k

и k

составляют в

3

1

2

 

сумме интеграл по поверхности k .

в) Аналогично для Hx и Hy областей справедливо:

 

ay

dxdydz ay cos d ;

a

x

dxdydz ax cos d .

y

x

G

 

G

 

Складывая почленно, получаем утверждение теоремы.

1). Векторная форма записи теоремы Остроградского-Гаусса имеет вид:

(ax cos ay cos az cos )d (

a

x

 

ay

 

a

z

)dxdydz ,

 

 

y

 

 

 

G

x

 

z

27

где cos ,cos ,cos - координаты единичного вектора внешней нормали.

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

 

 

 

(a

n0 )d diva dxdydz .

 

 

G

Поток векторного поля (вектора) через внешнюю сторону замкнутой поверхности равен тройному интегралу от diva по области G, ограниченной .

Применение теоремы Остроградского – Гаусса

Вычисление объемов.

ПРИМЕР. Пусть a(x, y, z) r ; r (x, y, z) ; div(a) div(r ) 1 1 1 3.

3dxdydz x cos d y cos d z cos d ,

Z

G

 

 

V

1

 

(xdydz ydxdz zdxdy) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

Вычисление потоков.

 

 

 

 

 

 

 

ПРИМЕР. Вычислите поток поля

2

i y

2

j z

2

k

a

x

 

 

через замкнутую поверхность .

 

 

 

 

 

 

 

Решение:

 

 

 

 

 

 

X

x2 y2 z2 R2 ,

Σ:

z 0 z 0 .

 

 

 

 

 

 

(2x 2y 2z)dxdydz =…

 

 

(a,d ) = divadxdydz =

 

 

 

 

 

 

G

G

 

 

 

 

 

 

 

 

{перейдём в сферическую систему координат}

 

 

 

 

2

 

R

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

R4

…= 2 d

d r

dr(2r cos sin 2r sin sin 2r cos ) =

 

 

.

 

2

0

0

0

 

 

 

 

 

 

 

 

 

 

 

ПРИМЕР. Найдите поток поля

 

2

i y

2

j z

2

k через

 

 

a

x

 

 

 

 

x2 y2 z2 R2 ,

внешнюю сторону полусферы: :

z 0.

Решение:

Воспользуемся результатами предыдущей задачи. Замкнем поверхность 1 поверхностью 2 , которая представляет собой часть плоскости XOY.

0

Y

xy

Z

1

 

0

1

 

 

xy

 

 

Y

X

2

 

2

28

 

 

 

 

R4

 

 

 

 

 

 

 

 

 

П П1 П2

 

 

, П1 П П2 ,

2 a,n0

d =

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x

2

y

2

2

 

2

 

2

d 0 ,

=…{ a

n0

 

 

z

) (0,0,1) z

}… = z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.к. на 2

z 0 и П1 П П2

=

R4

.

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

3.8. Инвариантное определение дивергенции

Пусть a a(P) - векторное поле, удовлетворяющее условию теоремы Остроградского – Гаусса. Пусть точка M - произвольная точка области G. Выберем поверхность , охватывающую область G. Из теоремы Остроградского – Гаус-

 

 

 

са следует, что (a,n0 )d divadxdydz .

 

 

G

Воспользуемся теоремой о среднем, согласно которой существует такая точка

 

 

 

 

 

 

 

 

 

 

 

ad

 

 

 

 

 

М1, принадлежащая G, что diva |M1

V ad ;

diva(M1)

 

 

, где

V

V

 

 

 

 

 

 

 

 

объем G. Пусть стягивается в точку М, тогда М1М, а diva(M1) → diva(M ) ,

 

 

 

 

(a,d )

 

 

 

diva(M ) lim

 

 

.

V

 

M

 

 

Поскольку правая часть выражения не зависит от системы координат (инвариантна), то инвариантно и данное ОПРЕДЕЛЕНИЕ дивергенции.

4. ЛИНЕЙНЫЙ ИНТЕГРАЛ В ВЕКТОРНОМ ПОЛЕ

4.1. Понятие линейного интеграла

Рассмотрим кусочно-гладкую кривую L и дугу AB (обозначение AB ) и

векторное поле a (ax ,ay ,az ), непрерывное на L. Разо-

бьем дугу AB произвольным образом точками A0, A1, …An на n частей. Обозначим ri - вектор, стягивающий

концы дуги Ai-1Ai . Выберем точку Pi Ai-1Ai . Найдём скалярное произведение (a(Pi ), ri ) и просуммируем

n

по всем участкам дуг Sn (a(Pi ), ri ). Вычислим пре-

i 1

 

Z

Ai

Ai 1

i 1

i

 

i 1

 

Y

X

дел lim S

 

 

lim

n

 

 

n

 

(a(P), r ) .

 

 

 

 

i

i

n

 

 

max( ri

0)

i 1

 

 

ОПРЕДЕЛЕНИЕ. Если этот предел существует и не зависит от способа разбиения дуги AB на отдельные участки и от выбора точки Pi, то он называ-

29

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]