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

книги из ГПНТБ / Палагин, Э. Г. Основы гидромеханики учебное пособие для метеорологов

.pdf
Скачиваний:
13
Добавлен:
20.10.2023
Размер:
8.24 Mб
Скачать

При этих условиях углы будут удовлетворять следующим соот­ ношениям:

 

 

Л6| =

Д а+ Д р,

 

 

Откуда

 

Дбг = А а — Ар,-

 

 

 

 

 

 

 

 

 

 

 

 

Аа = 7j- (А б| -f-Абг )',•

 

 

 

 

1

 

 

 

 

 

 

 

Д р = — (А0 ,-Л(-Ь),

 

Разделив эти равенства

на промежуток временй At и перейдя

к пределу при А/ -» 0,- получаем:

 

 

 

 

 

d z

 

1

,( d & t

,

d

\

 

 

 

 

 

 

 

 

d t

 

2

*[ d t

' r I t

 

 

 

d 3

1 ( d \

 

d . b . .

 

 

d t

у

1

 

d

t

 

 

 

 

\ d t

 

 

В этих

выражениях

dz

представляет

собой

угловую скорость

 

вращения

отрезков против

часовой

 

 

d?

стрелки,

— угловую ско-

рость перекоса прямого угла,

У

I

I

I

,4

J С

о

Рис. 2.23

d'jj

и---- угловую скорость по-

ворота отрезков против часо­ вой стрелки.

Две последние производ­ ные можно выразить через производные по координатам от проекций скорости. Для этого рассмотрим движения отрезка АС, параллельного оси Ох, в направлении оси Оу, от­ влекшись пока от его переме­

щения

в направлении

осей Ох

и Оу (рис. 2.23).

Точка

А

центр

частицы—перемещается

X в направлении

Оу

со

ско­

ростью

vv, точка

С — со

ско-

 

dv

1

Ах-

 

ростыо У„+

Y

 

в

результате этого за малый промежуток времени

АI

точка

А проходит расстояние АА" = vyA(, точка С —

рас-

ю

/

dv,,

1

\

 

стояние CC — (

I'yH—

^-Д.т)ДЛ

Поэтому разность этих рас­

стояний имеет величину

 

 

 

 

С'С" = СС" -

С С =

dv

1

- Л -

• -=- Лл-Д/,.

 

 

 

 

ох

2

угол поворота отрезка в направлений от Ох к Оу равен

 

 

 

С С"

dv,,

Т Лл

а угловая скорость поворота

Д§, dv,,

Дt дх

Таким образом угловая скорость поворота отрезка, параллель­ ного Ох, в направлении от Ох к Оу определяется производной

dv,,

. Движение отрезка АС в направлении оси Ох может при­

вести к растяжению пли сжатию отрезка, но вклада в поворот отрезка около точки А не внесет. Точно так же не вызовет доба­ вочного поворота в плоскости хОу и движение отрезка в направле­ нии оси Oz (хотя оно и может вызвать поворот в плоскости xOz).

Подобным же образом нетрудно показать, что жидкий отрезок, параллельный оси Оу, поворачивается в плоскости хОу в направ-

d v x

лении от Оу к Ох с угловой скоростью ^ ■•

В результате этих рассуждений мы имеем:

5

dv,,

do2

dv,

dt

дх

ИГ

~dv

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

d$

1

I dv,. (

dvx

dt

2

l дх

ду I

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

2 \ дх 1 dy I

3’

41

Аналогичным образом легко показать, что fli и 02 отражают пер.екосы частиц соответственно в плоскостях yOz и xOz.

ГЛАВА III

ПОТОК и ДИВЕРГЕНЦИЯ СКОРОСТИ

§ I. ФЙЗИЧЕСКЙЙ СМЫСЛ ПО'ГОкА ВЕКТОРА СКОРОСТИ

•"V

1.Поток вектора. Представим' себе, что'в поле вектора а (х, у,

задана некоторая поверхность а (рис. 3.1). Будем полагать, что эта поверхность является ориентированной, т. е. в каждой точке ее задано определенное положительное направление нормали (opt

нормали н), причем это направление при перемещении вдоль поверхности меняется непрерыв­ но. Таким образом, мы пола­ гаем, что задана. не только геометрическая форма поверх­ ности но и показаны ее «внеш­

няя» (та,

куда

направлены

нормали)

п «внутренняя» сто­

роны.

Если

поверхность

является замкнутой, то приня­ то направлять нормали изну­ три вовне объема, ограничен­

 

ного

поверхностью.

Разобьем

 

поверхность

на

элементарные

Рис. 3.1

площадки

Да,-

п

составим

произведения

вида

А,-,Да,-,

 

где

о,-„ — проекция

вектора

а,-, взятого в произвольной точке площадки Да,-, на нормаль. Про­ суммируем все подобные произведения и найдем предел этой сум­ мы при условии, что все Да,- стремятся к нулю, а число их неогра­ ниченно возрастает. Этот предел

называется потоком вектора■а через поверхность а.

По своему смыслу это скалярная величина, численно равная количеству а, проходящему через поверхность а. В частном слу­ чае а может быть замкнута. Тогда <3 = фа,Да дает разность между

входящим и выходящим потоками, т. е., в конечном счете, накопле-

>

и не (или наоборот) а в объеме ■, стягиваемом поверхностью а.

2. Расход жидкости. Поставим себе задачей установить физ ческий смысл потока вектора скорости через заданную поверх­ ность о:

Q— | vnda,

Ь)

Для этого представим себе сначала, что поверхность о является неподвижной «контрольной» поверхностью, через которую жидкость способна свободно протекать. Легко видеть, (рис. 3.2), что через площадку da за единицу времени' проходит объем

жидкости, заключенный в цилиндре с основанием da и образую­

щей V. Величина объема равна vnda. В силу этого vnda соответ­ ствует объему жидкости, протекающей через da, взятому со зна­ ком «+», если жидкость течет в сторону нормали, н со знаком «—», если жидкость течет в противоположном направлении. Интеграл

| vnda равен разности между объемом жидкости, протекающей

(з)

в сторону, куда направлены нормали, и объемом жидкости, проте­ кающей в противоположном направлении. Таким образом, поток скорости через поверхность представляет собой объемный расход, т. е. объем жидкости, протекающей через эту поверхность (счита­ ющуюся неподвижной) за единицу временц. Равенство потока нулю при наличии движения возможно лишь в том случае, когда объем жидкости, протекающей в сторону, куда направлены нор­ мали, равен объему жидкости, протекающему «внутрь» поверх­ ности.

Поток скорости через замкнутую поверхность Q— С\)unda (нор-

С

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

43

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

Если вместо v„ рассматривать ро„, то величина /VI = |

?v„da

а

 

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

/VI = ф p v ndcst и м ы п о л у ч а е м м а с с у , в ы х о д я щ у ю ( в х о д я щ у ю ) и з

Объема, также за единичный интервал времени. Ясно, что в слу­ чае несжимаемой жидкости накопление массы в объеме т не имеет места.

§ 2. Д И В Е Р Г Е Н Ц И Я С К О РО С Т И

1. Физический смысл дивергенции скорости. В курсах матема-

тпческого анализа* дивергенция вектора а— {ах, а„, ог) опреде­ ляется как

г)ал.

da v

d a .

illv а ■ ~7ьГ

dv

~дг

или, через посредство оператора Гамильтона, как скалярное про­ изведение:

(I iVci— (V -а).

Там же выводится формула связи между потоком вектора и его дивергенцией (формула Гаусса-Остроградского):

ф a „ d a = \ div a dr,

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

div а=

dv _

dv,,

dv,

(3.1)

—_

j___

L

j------ i

 

dx

dy

 

' dz

 

* См., например, Смирнов В. И. «Курс высшей математики», т. II. ГИФМЛ,

М., 1958.

44

или

div v = (V - v )

\ vnda— | divn^T.

з

(З.Г)

(3.2)

Последнее выражение расшифровывается следующим образом:

объем жидкости,

выходящий за единицу

времени за

пределы т

(или, наоборот,

остающийся

в нем) равен объемному

интегралу

от дивергенции

скорости. (В случае вытекания

дивергенция отри­

цательна) .

 

 

 

интерпретации

->

Отсюда легко перейти к физической

div и. Для

этого, воспользовавшись теоремой о

среднем, вынесем в (3.2)

{Iiv v из-под знака

интеграла.

Тогда

 

 

 

 

 

 

ф vnda— div и |

d-..

 

 

 

Поскольку ([) vndo равен разности входящих Q2 и выходящих Q\

о

 

 

 

 

 

величиной т=Дт,

потоков, т. с. Q2—Q, = AQ, т о , считая т малой

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

 

 

 

 

d

 

 

 

 

 

Переходя к пределу, получим

 

 

 

 

 

 

div v =

ЛQ

dQ

 

 

 

 

lim -г— =

- j -

 

 

 

 

Лт>0

dt

 

 

 

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

Все сказанное выше может быть перефразировано и по отно- ->

шенпю к секундной массе pv, для которой выражение для дивер­ генции должно иметь вид

 

dpv v

dpvy

dpv,

 

div pv — dx

r dy

dz

 

Физический смысл

•>

 

>

Просто вместо

div pv тот же, что и

div о.

объема в приведенных

выше

определениях будет

фигурировать

масса.

 

 

 

 

15

2. Запись плоской дивергенции в натуральных координатах.

Дивергенцию двухмерного движения ciivui будем называть пло­ ской. Для горизонтального движения (параллельно хОу) она за­ пишется как

d v d v .

div сг

дх dv

Запишем ее в натуральных координатах, направив, как и ранее, ось Ох (/) вдоль горизонтальной линии тока, а ось Оу (п) по нормали к ней влево от направ­ ления движения. Тогда имеем

(рис. 3.3).

 

dv

d't

 

(liv :l|~

~Ж ~*r V ~dn '

(3-3)

I I •

 

CjV

при­

Производная ——, равная

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

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

убывает. Производная

представляет собой угол поворота

дп

(в радианах) вектора скорости при смещении по нормали к линии

п

гока па единицу длины. Величина v —— представляет сооои часть

дп

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

дЬ

—> 0, т. е. лнинн тока в направлении движения расходятся, и

дп

величиной отрицательном, если д'? < 0, т. е. линии тока сходятся.

дп

В виде иллюстрации приведем примеры дивергенции модуля и дивергенции направления различных знаков (рис. 3 .4).

46

В общем случае негоризонтального движения дивергенция скорости может быть представлена в виде

 

 

 

д vy

О V.

д v.

 

div

дх

dv

Hz =

div t'i + H z ’

д v v

д г\.

представляет

собой дивергенцию гори­

где div«,=

- 3

 

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

 

OV

 

dv.

div V-

НГ "

с/п

(3.4)

H z '

где с — модуль скорости

горизонтального

движения.

о)^ > 0 -

и &

= 0

dl

on

dl

on

 

в)Ж < 0 и Ш >0

■ > ж > 0 v i l < 0

Рис. 3.‘!

3. Определение вертикальной скорости при движении несж маемой жидкости. Из физического смысла дивергенции скорости ясно, что для несжимаемой жидкости она равна нулю, ибо в этом случае не может быть накопления или уменьшения жидкости в единичном объеме. (Входит столько же сколько и выходит). Таким образом,

div V— дх

d v v

d v z

0 .

 

+ HV

Hz

 

(В случае плоского движения

- ►

д г \

d v v

 

div t>i=

—з----1— г-^-=0).

 

 

дх

dv

При горизонтальном движении, пользуясь натуральными коор­ динатами, будем иметь

dv

+ v

Й

= 0.

lil

 

д/г

 

Из последнего соотношения следует, что в случае плоско-па­ раллельного движения несжимаемой жидкости положительная дивергенция направления (расхождение линий тока в направле­ нии движения) сопровождается отрицательной дивергенцией мо­ дуля, т. е. убыванием модуля скорости вдоль линий тока (рис. 3.5, а). Напротив, при сближении линий тока в направлении движения имеет место возрастание скорости (рис. 3.5,6).

6)

aj

Ь)

в)

 

 

e iiie s iiu

 

 

h i

h i -

Рис. З.й

 

Рис. 3.0

 

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

пользоваться равенством divo = ().

Если направить ось Oz по вертикали вверх, го поток можно представить в форме

,d v 2

 

 

die с'Н

0 .

 

где div i>i— дивергенция

горизонтальной

составляющей

скорости,

d v x

дъ\,

dv

-{-г1

дв

 

равная либо

л - ^ ~ , либо

. Таким образом,

 

 

dv .

Cl

 

 

 

____ * --- __

 

 

Это значит, что положительная

дивергенция горизонтальной

составляющей скорости

(т. е. вектора скорости ветра на

синопти­

ческих картах) сопровождается уменьшением вертикальной со­

ставляющей

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

снизу

вверх (рис.

3.6),

Напротив,

отрицательным значениям

div щ

соответствует

рост

vz с высотой.

 

 

 

48

Полагая начало координат находящимся на земной поверх­ ности и интегрируя последнее равенство от нуля до некоторой вы­

соты /г, имеем

Л

 

vz(h)

— j

cliv

dz,

 

 

 

U

 

 

 

ИЛИ

 

 

 

 

 

 

vz(h) = —h (div щ) cp ,

 

где (divo)Cp

есть среднее значение

дивергенции

горизонтальной

составляющей

скорости. Так, например, в области повышенного

давления в слое трения

имеет

 

■"V

(рис. 3.7,а), что

место divu,> 0

обусловлено действием сил трения. Этому соответствует v,(h) <сО,

т. е. в области повышенного

давления имеют место нисходящие

движения воздуха. Напротив,

в области пониженного давления,

гдебш щ сО (рис. 3.7,6), vz( h ) > 0.

а)

6)

и 1 1 | L J !

t I I t t t t f

Рис. 3.7

Приведем несколько примеров грубой оценки порядка вели­ чины иг.

Пр и м е р 1. В области антициклона с круговыми изобарами радиальная составляющая скорости на расстоянии г=500 км от центра в среднем равна 1 м/сек до высоты h — 1 км. Оценим сред­

нее значение скорости нисходящих движений на высоте 1 км.

—>

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

div t'i =

2т. г h v r

2vr

■ кг2 h

г

 

4 Зак. 112

49

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