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

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

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

Так как J0= v-Q -too или 0 — скорость, Q — секундный рас­ ход), то при Q-- 0, о ■■ . Это означает, что в случае весьма узкой щели, скорость велика, а расход мал. Истекающую струю можно'

трактовать

как пограничный

 

слон

в

окружающей

ее

большой

массе жидкости. Последнюю считаем неподвижной,

откуда слс-

дуст

dp

п _

 

 

 

 

 

 

 

 

 

 

 

 

1огда в случае стационарности уравнения движения

и неразрывности, как и ранее,

примут вид:

 

 

 

 

 

 

 

 

Он

 

ди

'

 

О2и

 

 

(12.2:5)

 

 

 

 

Ох

 

оу

 

дц- '

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Он

 

 

O v

 

 

 

(12.24)

 

 

 

 

77~

 

~ду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(жидкость несжимаема).

 

 

 

 

 

 

 

 

 

Опять вводим функцию тока:

 

 

 

 

 

 

 

 

 

 

и

сМг

 

Г

 

гЛ'

 

 

 

 

 

 

 

Оу ’

 

 

Ох '

 

 

 

За пулевую линию тока г!' ==()

примем ось Ох.

Тогда уравнение

н граничные условия запишутся как:

 

 

 

 

 

 

 

 

tFF

<2ЛГ

onF

ЛЧГ

 

(РЧГ

 

(12.25)

 

 

fly

Ох fly

Ох

ду­

 

ду3

 

 

 

 

 

 

 

 

 

 

0 Ч !

 

 

— о

 

 

 

 

 

 

 

 

 

 

ду-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(12.2(5)

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~оГ

 

>

"X

 

 

 

 

 

 

Первое

из условий

(12.26)

следует

из симметрии

струн, что

 

Ov

второе

вытекает

из

того,

что v -» 0 , когда

означает

^ - = 0 , а

V—►±

. Поставленной

задаче

 

отвечает

тривиальное

решение

Чг = 0.

Однако это противоречит заданию

конечного импульса, ко­

торый

при

отсутствии внешних

сил

 

= 0 должен

быть конеч­

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

(12.23), которое с учетом

(12.24)

может быть записано в виде

 

0_

 

 

uv)

0 4

 

Ох

 

 

ду2

 

 

 

 

Интегрируя его

от

- ■ до

+ со , получим

d

)

ytddy+lyuv]

г

ди

с/х

!

ду

 

00

 

 

 

 

■НО

Так как при у

зо u = v =

ди

 

О, то

 

ду

 

 

 

 

 

 

 

 

 

d_

 

tddy

0.

 

 

 

dx

 

 

 

 

 

 

 

 

 

Откуда

 

 

 

 

 

 

 

uw2r/t/ — c o n s t =

/ о,

 

пли

 

 

 

 

 

 

 

\ ■'

j дЧг у

,

,

(1-2.27)

 

( d 7 )

d !,= u

 

 

 

Cv_

 

 

 

 

 

Вводя характерные

масштабы -v~.v0, г/ —//о, Мг~ Ч гц п перейдя

к безразмерным комплексам,

можно

(12.25) —(12.27) записать как

ачп

ду 6

Тогда

cPWi,

d'V6 d2W6 _

va

(12.28)

 

1 1

дхй

ду %

Ч'оу0 ду 1

 

 

 

 

^'мо,

 

 

0.

 

 

 

дуЪ

Уй-°

 

 

 

 

 

(12 29)

 

дЧ<6

у-

>

->0;

 

 

 

dy<,

 

 

 

х

 

 

( 0 W

 

 

 

 

(12.30)

J

\ Оуй5V

liy- —

 

 

 

 

Л‘

 

у

•'*о

Л

 

 

 

 

Уо ' Ч'и Уо' р 'Г

 

Поскольку масштабы ,v(), /у(), МП не заданы, то комбинируя без­ размерные переменные и комплексы, придем к зависимости вида

v 70 х

•Ч'о

/ (’О,

I де

.V

•V

141

Перейдя к повой переменной, вместо (12.28) —(12.30) получим

 

rfV

1

<d ± X

, d4_

 

(12.31)

 

 

 

LV d-ц I

 

: / ' d r

0 ;

 

 

 

 

 

 

 

 

 

 

cPf

 

=

0.

 

 

 

 

 

 

 

drf

 

 

 

 

(12.32)

 

 

 

 

 

 

 

 

 

 

 

 

 

df

 

 

 

 

 

 

 

 

 

 

,,>±

 

0 ;

 

 

 

 

 

 

 

(1 t\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i f f ) " - ’ -

 

 

 

(12.33)

 

 

 

 

 

 

 

Опуская

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

уравнения *,

выпишем лишь

окон чателы1ый результат:

 

 

 

 

 

 

 

 

'Г =1.651 у Л

Л X

tli

0.275 \

'

Л -

. JL

 

Расход в струе

 

 

 

 

 

 

 

-v=3

 

 

 

 

 

 

 

 

 

 

 

Q -= >Г (.о) -

«Г ( -

оо) =

3,302

Ay f —^ —

 

Скорости:

 

 

 

 

 

 

 

 

 

 

 

,/ = 0,454

j /

 

 

 

 

 

1

 

 

 

 

f/’ v.v

 

cli2 (0,275 v| Г

 

°,55°

j / ф - ■

 

0,550 l]

 

th (0.275 p)

 

 

cli2 (0,275 p)

 

 

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

пропорционально Л'ь.

 

 

 

 

 

 

 

 

 

.Максимальная скорость при р = 0

 

 

 

 

 

 

 

„,,, = 0,454

1

/ - Д

 

 

 

 

 

 

 

 

 

 

 

ГУ- V.V

 

 

 

Приняв

за

полуширине

 

струп

b значение

//, при

котором

// = 0,0161 //,„, получим

 

 

 

 

 

 

 

 

 

Ых) ~ 10,91

Тогда

и

Решение можно найти и литературе (см. сноску на стр. 139J.

М 2

Р а з д е л III

ДИНАМИКА ИДЕАЛЬНОЙ ЖИДКОСТИ

ГЛАВА ХШ

УРАВНЕНИЯ ИДЕАЛЬНОЙ ЖИДКОСТИ

§ I. З А П И С Ь У Р А В Н Е Н И Й В Ф О Р М Е Э Й Л Е Р А

Выше мы определили идеальную жидкость как жидкость, ли­ шенную трения. По, во-первых, в реальных средах вязкость всегда есть, во-вторых, если для какого-либо случая нами получено реше­ ние уравнении с вязкими членами, а затем осуществлен переход v —>(), то полученный в итоге результат нс всегда эквивалентен ре­ шению уравнении, в которых члены с вязкостью опущены сразу. Поэтому более правильно говорить о модели идеальной жидкости, справедливой в тех случаях, когда вязкие эффекты малы в срав­ нении с прочими силовыми факторами.

Уравнения движения идеальной жидкости можно получить из

uL

(12.4) при R e = —— >0 0 . Это означает, что либо велики харак­

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

стенки.

Итак, устремляя Re-»-'4, мы получим систему уравнений ди­ намики идеальной жидкости. Перейдя после этого к размерному виду, будем иметь:

d v ,.

г

dv ,.

,

,

<>vx

 

 

~

дх

'

■'

dy

 

 

 

 

 

0 v y

+ Tv

dvy

 

 

dvv

+

Vz

dt

~dx

^

 

*y

dv,

+ И-

dv,

1

 

 

+

ri

~df

dx

'

-v

dy

dv v dt

dz

1!

VQ1^

Ъо

 

1 dp

0dx

1dp

9

dy

1

dp

0

dz

1 4 3

В в е к т о р н о й ф о р м е :

 

 

 

 

uv

{г-х) г =

др

■г F.

(13. Г)

dt

дп

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

Уравнение неразрывности имеет вид такой же, как н ранее.

 

- ^ - + p i i i v w =

0.

(13.2)

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

равна пулю, так что оно имеет вид

 

 

С

dT -\-р div и — div

ОТ

(1о.З)

v

dt

г)и

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

замкнута:

 

 

 

 

Р = / ( / Л Л .

 

(13.-1)

В отношенпн замыкания системы можно сказать

то же, что н

в случае вязкой жидкости, т. е. для баротроппой жидкости урав­ нения (13.1), (13.2), (13.4) могут быть разрешены без использова­ ния (13.3), а в бароклпнноп среде следует решать всю систему

( 13.1)-(13.4).

§ 2. НАЧАЛЬНЫЕ И ГРАНИЧНЫЕ УСЛОВИЯ

Начальные условия, как и ранее, означают, что для исходного

момента времени t = l0 заданы

поля всех

искомых величин как

функции координат, т. е. vx(x, //,

т, /0) = / , (.у, //,

г), р(х, //,

г, /п) =

= / 2(.v, у, z) н т. д.

 

 

 

 

Различают кинематические н динамические граничные условия.

I. Кинематические граничные

условия

-это

условия,

опреде­

ляющие значения скорости на границах потока.

 

 

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

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

144

ты между жидкостью и стейком, а это противоречит представлению о сплошном заполнении жидкостью пространства. Таким образом, должно иметь место ц„ = 0, т. е. скорость должна быть направлена по касательном к стенке.

Эти условия можно выразить в аналитической форме. Если )(х, у, z ) —- 0 представляет собой уравнение поверхности стенки, то

орт п нормали к поверхности имеет проекции:

Of

Of

fix

ду

Условие v„ — 0 равносильно равенству (v-n)— 0, т. е.

dj ,

di '

г

°i - О

 

~дй т

 

Тг

Если границей потока является подвижная непроницаемая стенка или свободная поверхность, то, поскольку жидкость не может протекать сквозь стенку или свободную поверхность п сплошным образом заполняет пространство, нормальная состав­ ляющая скорости v„ любой частицы жидкости, соприкасающейся со стенкой пли свободной поверхностью, и нормальная составляю­ щая скорости v„ ст соответствующей точки стенки пли поверх­ ности должны быть равны:

V,, V ,1 ст .

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

lim г (х, у, г) = v 0.

X>сс

2. Динамические граничные условия— это условия, определяю­ щие значения давления на границах потока. Исходя из закона о равенстве действия и противодействия, получаем, что давление, оказываемое на жидкость внешней средой или стенкой, и давле-

10 Зак. 112

145

нне, испытываемое частицами жидкости, расположенными на стенке пли на свободной поверхности, должны быть равны по ве­ личине.

Отсюда следует, что на свободной поверхности давление жид­ кости равно атмосферному давлению.

§3. РАЗЛИЧНЫЕ ФОРМЫ ЗАПИСИ УРАВНЕНИИ ДВИЖЕНИЯ

I.Уравнения движения в натуральных координатах. При рас­

смотрении

плоского движения

удобно

пользоваться

за­

писью уравнений движения в натуральных

координатах. Исполь­

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

 

 

о:•

 

dv

 

 

I

op

 

 

 

 

 

 

 

di

 

'

 

о д ! '

 

13.5)

 

 

ten

 

 

 

 

1

dr.

 

 

 

R

 

F

 

 

 

 

 

di

 

 

 

 

dn

 

 

где l:i

и Fn

касательная

и

нормальная

проекции напряжения

массовой силы.

 

 

 

 

 

 

 

 

 

 

 

2. Уравнения движения в форме И. С. Громеко. Преобразуем

левую

часть

первого

из

уравнений

движения в форме Эйлера

(13.1):

d

v

rk х

 

 

d v ,

 

 

dv л

 

 

 

 

 

 

~дГ

■' Vx' ~дхТ~

 

Г"v ~ dдуtf

Г■' ~ д Т

Для этого возьмем тождество

 

 

 

 

 

 

 

 

 

 

г-

 

е‘л-J—j—zjt/2 L*r2

 

 

 

дифференцирование которого дает

 

 

 

 

 

 

(

- -

)

 

 

I

ГУ

дс

+

r.

dv.

 

 

~д.\

dx

~ d T ’

 

дх

' 2

>

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d vx

^

о

I

o-

,

dvy

 

 

dv,

 

v ~dx

Hie

\ 1 Г

J

ГуИ7

V*

 

~d~x~'

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

uv х

д

г)д

dvz \i •

n v

X \

~дГ

Их -9 /

dz

дх ' h " \ дх

oy

/ °y

146

Учитывая, что

[aXb] = (aubz—azby) /+ (azbx—axbz) /-f (axb,,—avhx) k.

мы можем предыдущее выражение записать как

d v ,.

д ( v~

д(

дх

— l + [rol t'Xi'l.v,

 

и первое из уравнений движения в форме Эйлера приобретает вид

dv v

д

|1 °2 \

- -

1

dp

(13.6)

~дГ +

Hie

\~2

)+ [rotvX v]x = X

fj

dx

1

 

 

 

Аналогичные преобразования, примененные ко второму и третьему уравнениям, дают:

dv.,

д

1f

v~

\

1

dp

(13.6')

Ж + ~W

(

j + [ r o l t ' X t ' ] //= y —

0

dy '

 

 

 

dv.

д

/

V1

\

]

dp

(13.6")

И Г +

~дг

 

 

J + [rot o X » ]r= - Z —

0

П Т ’

 

 

 

 

Эти три уравнения носят название уравнений движения в форме Громеко. Их можно представить в виде одного векторного уравне­ ния:

ди

grad — + [ r o t o X f ] = ^ -

grad р

(13.7)

ПТ

 

 

 

 

3. Уравнение А. А. Фридмана. Возьмем ротор от обеих частей векторного уравнения движения в форме Громеко:

rot

dv

-rot grad

) -frol [rot oX^l =rot F — rot

ПТ

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

..

. да

д

-*

 

 

 

1)

rot -дг =

ж;

го* я ;

 

 

 

 

д(

dt

 

 

 

 

2)

rot grad tp = 0;

 

 

 

3)

rot (aX~b) =

(XV) a — (a- V)

b+a div b — b div a;

4)

rot (cp a) =<p rot a+[grad cpXa]-

 

 

Учитывая эти соотношения,

получаем

следующие выражения:

..

—>

д

—>

дЯ

-■*,

 

. dv

*

;

•)

го1 ~ д Г ~

ж

r o i v = n f

(y -

rotcx

10*

147

2) ro t g r a d | Д^— j - - 0:

3) rot (rot v'Xv | — (c1• Л) о -- (O . c-f-У • div v\

4) ml

| Д-

grad p j — — -Д |gradp

. grad p | = Дг |grad/;X

1 ,ц;гad p ] .

 

эти выражения в наше уравнепне п учитывая, что

Подставляя

 

 

()U

("Д V jii =

dil

 

 

 

dt

' d t

получим

 

 

 

 

 

 

 

dil

(£2 ■V | t*-J-£2 div v = rot F-\-

1

grad дХ.grad p |. (13.8)

dt

 

 

 

 

 

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

Уравнение (13.8) показывает, что изменение вихря в точке может происходить за счет: конвективного переноса (член

(у-\7)У), деформаций жидкого элемента ((О • V) и), изменения

его объема (Odix'c), а также ввиду действия внешних сил (если

ro lF -^О) п наличия фактора бароклпнпостп | 4 r ['grad />Xgi‘ad р] j.

Очевидно,

что в случае баротропной

среды, поскольку поверхности

/)= consl

и р= const

совпадают, векторное произведение равно

нулю.

 

 

 

 

Если

массовые

силы имеют

потенциал

( / • '= - - grad

rot F— —rot grad U= Q) п жидкость баротроппа, то левая часть уравнения (13.8), именуемая гсльмгольцианом, равна пулю. Эго обстоятельство обуславливает удобство исследования движений идеальной жидкости в указанной форме записи.

ГЛАВА XIV

ИНТЕГРАЛЫ УРАВНЕНИЙ ДВИЖЕНИЯ ИДЕАЛЬНОЙ ЖИДКОСТИ

§I. УРАВНЕНИЕ БЕРНУЛЛИ В ОБЩЕЙ ФОРМЕ

1.Функция давления. Функцией давления называется величина

(14.1)

148

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

Установим физический смысл функции давления. Результиру­ ющая поверхностных сил, приложенных к единице объема, равна— —grad/?, а результирующая поверхностных сил, действующих на

°тас1 р единицу массы, равна - ------—. Работа сил градиента давления

при перемещении единицы массы на элементарный отрезок dr равна

KradP

' {др

1 др

rfv

dz

dp

 

% dX+

f

t

 

 

 

OX

 

Oil

'

dx

 

 

а при перемещении

единицы массы

из

точки

А с

давлением р

в точку В с давлением, равным

нулю, определяется

выражением

 

 

О

 

 

Р

 

 

 

Г/grad/;

л и>)

V

!'

I)

 

Р

Таким образом, функция давления представляет собой работу, которую должны совершить силы градиента давления, чтобы пе­ реместить единицу массы жидкости из данной точки в точку с ну­ левым давлением, пли, что то же самое, работу, которую нужно приложить к единице массы, чтобы, преодолевая силы градиента давления, переместить ее нз точки, где /; = 0, в данную точку.

Далее, элементарная работа сил градиента давления, прило­ женных к единице массы

i>

ciP(p),

о

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

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

2

др_

dx

- 2

dv + i f i

dz

о

дх

 

dij

'

dz

 

 

ад-

dx-

dP

dv

dP

dz

 

 

dy

 

dz

 

149

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