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

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

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

т е. приращение циркуляции скорости по жидкому контуру при перемещении контура за единицу времени, условно называется

ускорением циркуляции.

Наряду с вектором скорости в каждой точке потока в данный момент времени можно определить вектор ускорения соответствую­ щей частицы и тем самым выделить поле ускорения. При этом цир­ куляция ускорения по контуру в данный момент времени опреде­ лится выражением

0)

Теорема Томпсона. Ускорение циркуляции равно циркуляции ускорения.

Рис. 5.8

Рис. 5.9

Для доказательства произведем дифференцирование в правой части равенства (5.10):

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

элементов контура при его перемещении. Производная — б/ вы- d t

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

70

пройдет в

направлении

касательной расстояние

о/Д/, конец--

расстояние

/ , dv, „ Л ,,

элемента ока-

| г», -|------ о/

Л/, и приращение длины

дю

жется равным — - Ы:\1 . За единицу времени удлинение составит

д!

величину — L. Ы, Таким образом,

д/.

ill dl

<0

W>

 

(/)

поскольку

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

по

замкнутому

контуру в том

случае, когда под интегралом находится полный дифференциал, дает нуль. Поэтому

d

V , U :

dt

dV‘ оl. dt

< l )

(0

что и доказывает теорему. Итак,

dr

Д dv

й / .

(5.11)

dt

dt

 

 

(п

ГЛАВА VI

ФУНКЦИЯ ТОКА И ПОТЕНЦИАЛ СКОРОСТИ

Для определения поля скорости в общем случае нужно найти три скалярные функции координат и времени vx, vu и vz. Однако в некоторых случаях движения существует такая скалярная функ­ ция, производные которой по координатам равны проекциям ско­ рости. В этих случаях для определения поля скорости достаточно найти одну лишь эту функцию, что значительно упрощает задачу. Такой функцией в одних случаях служит функция тока, в других— потенциал скорости.

§ 1. ФУНКЦИЯ ТОКА ДЛЯ ПЛОСКО-ПАРАЛЛЕЛЬНОГО ДВИЖЕНИЯ НЕСЖИМАЕМОЙ ЖИДКОСТИ

1. Понятие и важнейшие свойства функции тока. Рассмотри плоско-параллельное движение. Функцией тока будем называть такую скалярную функцию координат и времени, градиент кото-

71

рой равен по модулю скорости, но повернут относительно вектора скорости па 90°, в одном п том же направлении во всех точках

(рис. 6.1 ).

Располагая осп Ох н Оу так, чтобы поворот от первой осп ко

—►

второй происходил бы в том же направлении, что поворот от о к grad Чг, имеем:

gradx 4J’ = —vlh grady 4r= yv.

Иными словами,

Чг

это такая функция, производные которой

по координатам удовлетворяют соотношениям:

 

д у

( 6. 1)

 

дх

 

 

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

что

если ЧЛ есть функция тока, то и 4fi + C,

где С — произвольная

постоянная, также будет функцией тока.

Иначе говоря, функция тока определяется с точностью до аддитив­ ной произвольной постоянной.

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

аdvx

(—■»*)= дх

72

Рис. К.2
дФ
уп= ——. Поэтому
01

так как лишь в этом случае смешанные вторые производные-——

ах оу

д щ «

и ------ , определяемые из (6.1 ), будут равны. Это значит, что

ду дх

функция существует в том и только в том случае, когда выпол­ няется условие

dvr

ov,,

 

=div Vi— 0.

Иначе говоря, введенная вышеуказанным образом функция тока существует лишь в случае плоскопараллельного движения несжи­ маемой жидкости. Следует указать, что функции тока, опреде­ ляемые, правда, по-иному, могут быть введены также и для пло­ скопараллельного движения сжимаемой жидкости и для осесимме­

тричного потока.

Важнейшим свойством функции тока является следующее: во

всех точках линии тока функция тока имеет одну и ту же вели- —> —>

чину. В самом деле, поскольку grad lFJ_t\ а вектор v направлен по

касательной

к линии тока, то

проекция grad XF

на

направление

линии тока,

т. е. производная

дФ

и,

значит, в на-

----- , равна нулю,

 

 

dl ..

 

 

правлении липни тока 'F не меняется. Это означает, что эквиска-

лярные

линии функции Чг, т. е. линии

XF —const, представляют

собой линии тока.

 

 

grad ф

Докажем другое свойство функ­

ции тока: разность

между

значе­

 

ниями функции тока в двух точках

 

равна объему жидкости, протекаю­

 

щей за

единицу ' времени

между

 

этими точками, рассчитанному на

 

единицу высоты потока.

 

 

Действительно, указанный объем

 

равен, очевидно, потоку Q скорости

 

через цилиндрическую поверхность

 

а, опирающуюся на линию А, В,

 

которая

соединяет

обе точки А и

 

В, причем высота поверхности а равна единице.

Q = ) vnda— j vndl.

(=)(Ав)

Но из рис. 6.2 видно, что n„ = grad,'F, т. е.

Учитывая (6.1), легко выразить Q. через Чг. А именно,

 

(

dry

di\.

I d1Ч'

 

(УГ \

 

 

 

дх

ду

I

дх-

г

ду- /

 

Пели

движение

безвихревое

(Ог =

0),

 

d- Ч

d'1Ч

то ЛЧ1— ----

----—0.

Фуикцня

тока в этом

случае удовлетворяет

дх-

ду-

уравнению

Лапласа,

т. е. является гармонической функцией.

 

 

 

2 .

Определение

функции

тока

по заданному полю скорост

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

Найдем функцию F, удовлетворяющую следующим

условиям:

~ = - xy3= f u

~ = 3 * y + 2 # = / s .

(6.2 )

ах

ом

 

Функция F существует не при любых / 1 п /2, а лишь в том слу­

чае, когда последние функции удовлетворяют условию

 

З/i

_

 

ду

дх

 

d-F

гак как только в этом случае, смешанные вторые производные------

дх ду

и -—— будут равны. Легко убедиться в том, что в нашем случае

ду дх

это условие выполняется и, значит, функция F существует.

Для отыскания ее используем сначала первое из условий (6.2). Из него получим

F = \ 'Ixifclx = х2у3+ С , {у),

где о величине Ci(y) можно сказать лишь то, что она не зависит от х, но, вообще говоря, может зависеть от у. Действительно диф­

ференцируя последнее равенство по х, тотчас

же убеждаемся

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

(6.2 ).

Для определения функции С] (у) используем

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

условий (6.2 ), которое требует, чтобы имело место равенство

[* V + C i(0 )l = 3х2у2+2у.

ду

Отсюда находим С\' (у) =2у, Ct(y) = у 2-\-С, где С — произволь­ ная постоянная, и окончательно

F = x 2y3-\-y2~\-C.

7-1

При отыскании функции тока задача сводится к определению

 

дЧ>

ду Vx.

функции VF, удовлетворяющей условиям —- = —щ, И

d ( - v v)

dvx

является

,\'оеднвшпсь в том, что— -г—— = ——, т. е. жидкость

ду

дх

 

несжимаемой, производим определение Чг в указанном выше порядке. Приведем примеры отыскания ЧТ

Q.v

а. Плоский источник. Так как в этом случае 0* = - 2— ,

v,i=

Qy

, то для

определения

\F имеем

равенства:

Л-2 -(-у 2

 

 

дЧ'_

Qy

г)Ч'

Qx

 

 

дх

л'2 + _у5

 

(6.3)

 

 

ду ~ х 2 + у2

Условие существования 'F при этом выполняется, так как вы­ ражения для vx и vy получены с помощью допущения о несжимае­ мости жидкости. Используя теперь первое из равенств (6.3), находим

iF =

г

Qy clx

—Qarctg — +Ci(y).

J

л2+ у2

 

у

Для определения С\{у) используем второе из равенств (6.3):

_д_

L

-Q arc t g ----- f-С, (у)

Qx

 

 

ду

х* -f V*

 

 

 

 

У

 

 

 

откуда С \(у)= 0,

т.

е.

С\ (у) =const = C2.

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

 

 

 

Ч/ =

—Qarctg-----(-const.

 

(6.4)

Уравнение

линий

тока

плоского

источника

найдем,

приравняв

функцию тока

постоянной

 

 

 

 

 

х

величине. Тогда получим Q arctg— ~г

-J-C, = const,

что

равносильно

равенству

х

 

 

— = С3 или у — С^х,

где Сц — произвольная

постоянная.

Отсюда

следует,

что линии

тока представляют собой

лучи,

исходящие

из

начала

координат,

вчем легко убедиться непосредственно.

Вслучае плоского стока функция тока имеет вид

х

( 6. 0)

Q arc tg — -(-const.

75

б. Безвихревое вращение жидкости. В случае вращения протип

часовой

стрелки

мы

имеем:

vv—

 

Су

Сх

- ---------,

v „ = -------- • llo-

 

 

 

 

 

 

Л-2 + у -

д-2 .i_ y'-i

этому

для определения Чг

служат

равенства:

 

 

o' M"

_

Сх

 

 

(Я1;

_

Су

 

<ix

 

х'1

у й

'

ду

л '2 + у-

Условие существования Ч1' выполняется, так как жидкость является несжимаемой. С помощью рассуждений, аналогичных предыдущим, находим:

Чг= —С

In (.v2+//2) -(-const.

((3.6)

Уравнение .шипи тока 4r

= const сведется,

очевидно, к уравне­

нию семейства концентрических окружностей

,v2+_//2=const. (Ре­

зультат совпадает с полученным ранее, см. § 3, гл. II).

 

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

иметь:

 

 

 

 

Чг =

С Iп (д'2-|-;/2) -(-const.

(6.7)

в.

Источннк-рбезвихревое вращение по часовой

стрелке (м

дель стационарного антициклона).

 

В случае сложного движения, складывающегося из двух двп-

 

—>

- >

 

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

d4‘

гРГ

г,,.4

- г,д .

дх

ду

 

 

Очевидно, что если Мг| и Чг9 -■ функции тока каждого из дви­ жений, то сумма их + будет удовлетворять написанным уравнениям, т. е. будет являться функцией тока сложного движе­ ния.

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

х

4F = —Q arc l g --- рС In (x2-pу2) 4-const. (6 .8 )

Уравнение линий тока в этом случае имеет вид

-Q arc tg х , С 1п (х24-у2) = const.

У

76

Если перейти к полярным

 

.V

= clg<f.

координатам, полож ив y

х2-{-у2 — г-, то будем иметь

 

 

 

Q ---- ®j

-^-С In /'2= const

 

откуда

_ _Q_ .

 

 

 

r= ke

ac \

(6.9)

где к — произвольная постоянная.

Легко убедиться в том,

что ли­

нии, соответствующие этому уравнению, представляют собой лога­ рифмические спирали.

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

ражение для функции тока в этом случае:

 

^''z^Qarclg — ~-С In (.v2+ Уг) + const,

(6.10)

V

 

а линии тока будут определяться уравнением вида (6.9), т. е. будут представлять собой логарифмические спирали.

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

§2. ПОТЕНЦИАЛ СКОРОСТИ

1.Понятие и важнейшие свойства потенциала скорости. Потенциалом скорости называется скалярная функция ср, гра­

диент которой равен скорости:

 

 

 

grad ф =

и.

( 6. 11)

Иными словами,

ср — это

функция,

удовлетворяющая

условиям:

ду

 

>

dtp =

dw

( 6. 12)

 

г\,,

г).х

 

~dv

~dz

 

или условию

 

dср vxdx-\-vvdy-\- vzdz.

 

 

 

(6 .1 3 )

77

Очевидно, что если некоторая

функция ср удовлетворяет равен­

ствам (6.11), (6.12) нлп (6.13),

то и функция <p +С, где С — про­

извольная постоянная, будет

удовлетворять им. Таким образом,

потенциал скорости определяется с точностью до аддитивной про­ извольной постоянной.

Эквнскалярные поверхности функции ср называются эквипотен­ циальными поверхностями.

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

дгх

d v v

dvx

d r г

dvy

dvc

ду ~

д х

'д г =

д х '

d z

<?v

(очевидно, что лишь в этом случае смешанные вторые производ­ ные от ср нс будут зависеть от порядка дифференцирования).

 

Равенства (6.14)

равносильны условиям

 

(Н\,

dvx

го1,?> =

0 .

dr.

chy

дх

ду

дх rotv V

VZ - rolvtf= -0

пли условию

rot 0 = 0.

С учетом этих равенств, пз теоремы Стокса следует, что

Г = ф p,dl — 0 .

Таким образом, для существования потенциала скорости необ­ ходимо и достаточно, чтобы движение было безвихревым. Отсюда вытекает, в частности, что в случае потенциального движения ли­ нии тока не могут быть замкнуты. Действительно циркуляция по замкнутой линии тока не может быть равна нулю, ибо суммирова­ ние проекции скорости идет по всему контуру с одним и тем же знаком. С другой стороны, в нашем случае Г = 0. Полученное про­ тиворечие и доказывает высказанное выше утверждение.

Укажем важнейшие свойства ср.

а. Эквипотенциальные поверхности и линии тока взаимно ортогональны, так как в каждой точке поверхности ср = const век-

тор gradcp — v нормалей к поверхности и в то же время совпадает

78

но направлению с касательной к линии тока. (Очевидно, что в на­ правлении движения гр растет). Это же означает, что изолинии ф

пвзаимно ортогональны.

б. В случае несжимаемой жидкости потенциал скорости являе ся гармонической функцией. В самом деле, условие несжимае­ мости означает, что проекции скорости удовлетворяют равенству

div г

 

дг:

(.г.

U.

дх

—L -I___ 1

 

д\>

dz

 

Вводя в него выражения vx, v,„ v: из равенства (6.12), полу­ чаем

д2ср

d2f

д"-о

(6.15)

'дх2 + ~ W

г 'dz2 ~

 

т. е. потенциал скорости представляет собой решение уравнения Лапласа (является функцией гармонической).

Поэтому задачи отыскания поля скорости безвихревого потока несжимаемой жидкости сводятся к отысканию решения уравнения Лапласа, удовлетворяющего граничным условиям, которыми одна задача собственно и отличается от другой.

2. Примеры определения потенциала скорости. Задача об опр делении потенциала скорости по заданному полю скорости сво­ дится к отысканию функции по ее частным производным. В слу­ чае плоско-параллельного движения задача решается подобно тому, как находится функция тока. Приведем примеры определе­ ния ф.

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

до

Qx

до

 

Qy

дх

х2 4- у 2

ду

д'2 -|- у 2

откуда получаем

 

 

 

 

 

(p = Q In (х2-\-у2) -(-const.

(6.16)

Эквипотенциальные поверхности (p = const определяются урав­ нением x2+ « /2= const, т .е. представляют собой семейство соосных круговых цилиндров, ось которых совпадает с осью Oz.

В случае плоского стока

<р— —Q In (х2+г/2)-(-const.

(6 .1 7 )

79

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