![](/user_photo/_userpic.png)
книги из ГПНТБ / Палагин, Э. Г. Основы гидромеханики учебное пособие для метеорологов
.pdfт е. приращение циркуляции скорости по жидкому контуру при перемещении контура за единицу времени, условно называется
ускорением циркуляции.
Наряду с вектором скорости в каждой точке потока в данный момент времени можно определить вектор ускорения соответствую щей частицы и тем самым выделить поле ускорения. При этом цир куляция ускорения по контуру в данный момент времени опреде лится выражением
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
так как лишь в этом случае смешанные вторые производные-——
ах оу
д щ «
и ------ , определяемые из (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) |
Иными словами, |
ср — это |
функция, |
удовлетворяющая |
условиям: |
|
ду |
|
VЛ> |
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