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

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

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

Г Л А В Л X U

ИСПОЛЬЗОВАНИЕ ТЕОРИИ ПОДОБИЯ В ЕИДРОМЕХАНИКЕ

§ 1. УПРОЩЕНИЕ УРАВНЕНИИ

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

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

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

Выпишем исходную систему:

d v

1

 

gk

X- vjyrad

dive

н vA г ;

(12.1)

-5 - + (v‘ V) — ------ grad р

ot

О

 

 

 

 

 

 

 

Ф.

f '/ di v v = 0:

 

 

( 12. 2)

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(12.3)

Из массовых сил

в (12.1)

 

мы учитывали только силу тяжести,

направление которой определяем ортом к.

 

 

 

Далее действуем

аналогично

тому,

как и

в рассмотренном

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

 

H'l'acl,-, pc, ; k

I-

( v, ща<1й divzv, j-

Aof>) .

(12.4)

130

Теперь

легко перейти к той или

иной модели.

Так, процесс

можно считать стационарным,

если,

Sh > 1

. Последнее

не просто

означает /0 -> оо, а

соответствует неравенству

L

 

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

Аналогично,

силой

тяжести

можно

пренебречь, если-=!— 1,

 

 

 

 

гг

что означает

т.

е- опять

же это

результат сопоставления

трех величин v0, L и g. Мы видим, что влияние силы тяжести па­ дает с ростом скорости п уменьшением характерного масштаба.

Вклад вязких членов определяется численным значением мно­

жителя

Влияние вязкости убывает при увеличении скорости и

масштаба, ибо — р 1 эквивалентно неравенству \<£v0L.

Критерий Эйлера всегда 0(1), ибо через силу градиента давле­ ния уравновешиваются все прочие силовые факторы, т. е. это по своей сути реакция среды на приложенные к ней усилия. Она,

естественно, существует всегда. Таким образом

v- = 0 ( 1 ) , или

Ро ~ РсТ'о2-

(12.5)

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

Прологарифмировав и продифференцировав по времени (12.3), получим

 

 

 

и

dp

 

 

dt

/.р

dt

Величина — — а2

как

мы

убедимся далее, есть квадрат ско-

I

(12.2) можно переписать в виде

ростн звука а. Теперь

 

1

dp

 

 

я*

131

Обезразмеривая это уравнение, получим

fv, div r0 = 0.

t<i?n ro dt6

Дли большинства апериодических задач для масштаба времени можно считать справедливой оценку i0~ L l v 0. Подставляя это равенство в последнее уравнение п заменяя />0, согласно (12,5), окончательно получаем

dpа

■Vs - £ r ?a i l i v i ’o = 0 . ( 1 2 . 6 )

Здесь введено обозначение А1 =

а

Величина /VI носит название критерия Маха п, как видим, пред­ ставляет собой отношение характерной скорости к скорости звука. Жидкость можно считать несжимаемой, если /VI . 1, т. е. v Д<7- Оценивая численные значения критериев, мы с учетом заданной точности, можем прийти к той или иной математической модели. При этом малые члены можно из рассмотрения исключить, что всегда ведет к существенным упрощениям задачи.

Следует иметь в виду, что иногда простое отбрасывание сла­ гаемыхможет повести к не вполне корректным результатам. Дело в том, что, производя оценки членов дифференциальных уравне­ ний. мы должны быть уверены, что их порядок во всей исследуе­ мой области сохраняется одинаковым ( < / > = 1 ), т. е. тем самым мы полагаем, что и решение является гладким во всех точках. Не­ приятностей следует ожидать при наличии в уравнении сингуляр­

ностей, или, если имеются значительные

градиенты величин

 

-чс

 

в какой-либо локализованной зоне (тогда

_

Х;

-У--1). В ча-

 

0

 

стностп, это имеет место при наличии резко пульсирующего ха­ рактера /.

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

/= Г 1) + а / 1Ч « / й т . . ■

Подставляя (12.7) в исходные уравнения, начальные п гранич­ ные условия, мы будем получать соответствующие математические формулировки задачи, отбирая последовательно члены при а п, а 1. а,2 и т. д. Это дает нам нулевое, первое, второе, и т. д. приближе-

132

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

Ниже мы проиллюстрируем метод малого параметра на при­ мере вывода уравнений пограничного слоя.

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

§ 2. ДИНАМИЧЕСКИЙ ПОГРАНИЧНЫЙ СЛОЙ

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

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

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

ди

\ II

ди

+

г

ди

 

1

др

 

1

д2и. ^

д1

Т-

ду

 

у дх '

 

 

ах

 

 

 

\ дх2 d i / l ’

ди

 

UV

 

 

ди

1 др

1 д2г’ , д2и \

~д!

"

ж + г

 

dil =

'/

 

ду

 

 

ду2 )'

 

 

 

 

 

ди

ди

 

 

 

 

 

 

 

 

 

 

дх

ду

~

 

 

 

 

где уд= н,

=

и,

v --- продольная

и

поперечная

компоненты

скорости.

характерные

 

масштабы;

,v~.vn,

г/— //0,

v ~ i 'n,

Введем

 

,х 0

; А. ~ р«г, •

«о

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

til

ИпI’o

 

 

ип

^

 

на —ф-, второго— на

д0

, а 3-го — на —- , будем иметь:

л"0

 

 

 

-v0

 

 

_Ai_

ди&

х °

- диб — ___Рп.

?£о_,_

иJ o д/й

0 дхй

'

ИуУ0

б ду0 ~

?«о

<5*о '

133

 

 

 

 

 

V

гУ«6

1 7 “'•О О'2и с

 

 

 

 

 

 

 

 

 

 

л"о У,

Ох?,

' У1 пп У/ б ’

 

 

 

 

 

У,

0 v б

 

dv

Л'о V,,

_

dv с

 

 

 

Ро х о

 

 

 

 

iintn die,

; "й

<3*л +

Уо«..

' 0

дуг,

 

 

 

? J W

’o

 

 

 

 

 

 

 

,

 

d2v 6

 

*л*0 0 -Vб

 

 

 

 

 

 

 

 

 

1 х 0ил Ох 1

 

У1 «о

ду g ’

 

 

 

 

 

 

 

 

 

 

du б , УУо оус,

 

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

дхй ' JV'o дус,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из

уравнения

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

ясно,

что :1Н—’ =

1,

а

из

условия

 

 

 

 

 

 

 

 

 

 

 

Уо Ни

 

 

 

 

равенства

порядков сил трения н

инерции

вытекает,

что

Уп«о

 

 

 

 

 

 

 

 

 

лу, и

 

 

 

 

 

 

Отсюда, вводя число Рейнольдса Re

о

 

получил:

 

 

о “

 

 

 

 

 

 

 

 

_

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У Re

 

 

У Re

 

 

 

 

 

 

Теперь после установления соответствующих масштабов, урав-

нения перепишутся в виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дий

ди5

 

ди б

 

дрй

 

 

1

УУ б

,

дЧи

 

 

ж й У «г,

 

 

 

 

 

 

 

 

 

1

 

 

 

Охс,

'

'(' д у ~

 

дх б ' г

Re

а д

+

dyi

1

(dv с

dv о

 

 

 

 

дРс> ,

 

 

1 д2ой ,

1

д2и й

Re

\дт6

+ Кб Охс,

:

й % б )~

 

ду 6

1

Re2

fog

 

 

 

 

 

 

 

_0£б

dvt

: 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

дхб

ду 6

 

 

 

 

 

 

 

 

Так как О(у0) = 5 1. где бд — толщина

пограничного

слоя,

являю­

щаяся

по

условию

малой

величиной,

то

 

 

можно

 

считать, что

-г-=. Отсюда следует,

что наличие пограничного слоя харак-

V Re

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

теризуется большими значениями числа Re. Принимая величину

—= - за малый параметр,

можно

формально разложить по нему

V Re

 

 

 

искомые величины, т. е.:

 

 

 

и- = и^

-4- -—-

и

°

6 +VRe

('

' ' ’ • ’

Ч — vp

1

г-’о 1У

У

 

У Re

 

 

Ро: -рТ

1

 

 

К Re Р? +

134

Подставляем эти разложения в нашу систему и приравниваем члены с нулевыми степенями при малом параметре. Это дает нам нулевые приближения:

ди6 ,

дий .

он &_

 

"с' дГй ^ ' 1'6

ду(,~"

 

П . _

°Ра

 

 

d (h '

<)рй

;

д*и6

дх'с,

:

ду1

дН(,

dv6

.

дхй

ду6

 

Аналогичный результат получается при Re -О. Возвращаясь

кразмерным переменным, получаем уравнения пограничного слоя

ввиде:

 

ди

,

Он

Он

_

1

Ор

д2и

(12.8)

 

~dt

г и

~г-

 

();/

 

о

дх

 

 

 

дх

 

 

 

 

 

 

 

 

ди

 

0 V

--=0.

 

 

(12.9)

 

 

 

 

дх

 

ду

 

 

 

 

В таком

виде

они

были

впервые

опубликованы в

1904 г.

.'I. Прапдтлем.

 

 

 

 

 

 

 

 

 

,,

On

 

 

 

 

 

 

 

 

 

Условие

ду

= 0

означает,

что

в

пределах

h o i раничного

слоя

давление не меняется т.

е., точнее,

в виду малой толщины его изме­

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

редь,

определяется

по уравнению Бернулли (см. ниже) как:

/ ; = —

rt iv2

— —

I) считается известным из решения уравнения

 

(и — п(х,

для внешней идеальной жидкости, где силами трения можно пре­ небречь) .

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

 

 

 

 

 

 

ди

ии

ди

 

ди

д2и

( 12. 10)

dt-

дх

ду

 

дх

ду2

 

 

 

ди_

да

— 0.

 

( 12.11)

 

дх

+

ои

 

В случае стационарности:

 

 

 

 

 

ди

,

он

-

du

д2и

( 12. 12)

дх

ду

 

ду ­

 

 

 

ди

 

д~о

о.

 

(1 2 .1 3 )

 

~дх

+

 

 

135

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

§Д И НАМИ Ч ЕСКИ И ПОГРАНСЛОЙ НА ПЛАСТИНЕ

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

и— const.

В

этом

случае

уравнение

(12.12)

упростится,

ибо

- da

= 0,

а

(12.13)

останется без

изменений,

так что имеем

си-

и ~Г

dx

 

 

 

 

 

 

 

 

 

 

стему:

 

 

 

ди

,

ди

 

dhi

 

 

 

 

 

 

 

(12.14)

 

 

 

 

дх

----- \- V -----------

= V

_________ *

 

 

 

 

1

ду

 

ду- '

 

 

ди

(12.15)

~~дх

 

Располагая начало координат па передней кромке пластины, направим ось .у по потоку, а ось у перпендикулярно к пластине. Тогда граничные условия примут вид:

и — 0 при у — 0 и 0 <. х

х ,

|

и —>и при у -* со н п

с о ^ х с о

(12.16)

.1

d'lr dxlr

Вводя функцию тока 4х и имея в виду, что и= — — , V— - — ,

дх г)у

приходим к одному уравнению относительно 4х, которое получается заменой и и о через Ч' в уравнении (12.14). Оно имеет вид

(MF

д2Чг

54х д2[[г

д:пУ

(19 17)

ду

дх ду

дх ду2

ду3 •

 

ранпчпые условия:

гР1г

О

при у — 0 п 0

-V .

да

 

 

(12.18)

РЧГ

 

 

и

при у — п

.V /

ду

 

 

 

Поскольку 4х определяется с точностью до аддитивной постоян­ ной. будем считать, что Чх = 0 при у ~ 0 п — сс< х < тс.

Введя масштабы

< lF > = vFo, < х '> = х 0, < у > = у о,

обезраз-

мерим (12.17) и (12.18). Получаем:

 

 

 

 

д'Ус, .

 

а

_

дхс

ду%

... Д*о_

 

дуб

дхбду6

 

упЧ>'0 ду1

 

6 ==0 при у 6

 

О н О

Л”б

 

 

Оуо

 

 

 

 

 

 

 

I 12.18')

а у б

и у а

 

 

 

 

 

 

or-

II

Л-,-

 

 

дус,

-п— при

 

 

 

 

 

 

 

 

 

 

В (12.17')

все члены поделены па общий .множитель

перед ле-

uoii частью

I]/3

,

 

 

цг,

 

 

 

— ^

а в (12.18') на-

1 .

 

 

 

 

хоУ;,

(12.17')

и

у о

безразмерные

комплексы

Полученные в

(12.18')

идентичны приведенным комплексным множителям

в которых

масштабные

коэффициенты

заменены

отношением

характерных

размеров величин. Поэтому из требования подобия вытекает их равенство единице. Таким образом,

v А ' н _И ЗВ __ .

(12.19)

У о ^ ~ П

 

Имеем два соотношения,связывающие масштабы трех перемен­ ных величин д'о, уa, 4J'(). Предположив первоначально, что масштаб л',, нами задан и равен L, определим г/о’н TV

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

 

vL

Ч'„ _

_Д_

 

 

Уо Чг„

it

У

У« п '

 

Теперь вместо

(12.19) имеем

 

 

 

 

 

=

Л Ь .

| .

(12.19')

 

\'l

и

*1'о

 

 

В (12.19') величина у0

входит лишь в один комплекс. Откуда

находим //о= 1

. а затем из второго равенства 4^0=

I vuZ..

\

и

 

 

 

 

Безразмерная функция тока должна быть функцией безразмерных координат:

У'а=/{Х(,, у6),

млн, с учетом полученных выше выражений для масштабов;

ч' = , / ы ! / ( - f , v I /

-

137

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

следует, что все безразмерные

комплексы

должны

объединиться

в такой комбинации,

чтобы нсключилось /.. Например,

ij'

: д-

' ,

-----— —

)

 

 

 

 

 

 

 

 

Y 4 iL \ I.

1

 

)

. Легко видеть,

что а = |1 = —

 

 

 

Тогда получим

 

 

 

 

 

 

 

 

 

 

[

v ИХ

■II

 

 

и

 

(12.20)

 

 

 

'IX

 

 

 

 

 

 

 

 

 

 

 

Мы

видим, что вместо

двух

в решении

задачи осталась

лини.

одна

комплексная

перемена

ц — у

1

"

. I акпм образом, рас­

|

 

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

Легко видеть, что

 

 

_dt]

_

у

1

и

 

 

dp

и

 

 

 

 

 

д\- “ ‘

2М'

I

v

 

ду

 

 

 

 

П оэтому:

 

 

 

 

 

 

 

 

 

 

 

 

 

d'Y

 

 

их

 

■01 —

 

их

^ф(л)

'Щх

do

(7)|

ду

 

 

 

 

 

 

ду

dt ,

dx

 

 

 

 

 

 

 

 

 

do

 

 

 

 

 

 

 

 

 

 

 

 

— и —

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d%

 

 

 

 

 

 

d2xY _

a

dlF J _—

()

j

do

— и

d2o

dn -

/

и

af2cp

dy2 ~

r)y

ду

I

~ “

dy

\

d-t\

~ W ~ d y ^ U

'ix

dif

 

 

 

 

_

d_

/

 

\ _

 

 

d?cp__

 

 

 

 

dip

 

dy

\ dip

J

dW

___

_d_

 

 

dx

dx

I

V U X

 

 

 

 

+ |/

>ux

do

 

£li

о

 

dr{

 

dx

V

'ix

dx?

 

 

 

 

'iu

C( r,) -r-

2

I/

- -

 

 

 

III

Г

 

 

dq

x

j =pMl) — T,

J ’

 

 

 

138

dx ch-

д

1

1 . /

Г

 

.

do

i

1

I /

!!L \ d

Ц,

Г

[?

Ы ~ •'

 

1 =

2

Г

-v 1dr{

 

 

 

 

rf<p"]

дг|

I

1

и

 

 

 

 

 

 

 

dr,\

ду

I

о

X

'

dr,®

Подставляя полученные значения производных в уравнение (12.17) п производя соответствующие сокращения, убеждаемся что оно переходит в обыкновенное дифференциальное уравнение вида

о

rf:,o

i ^

dr?

1 ^

d'-v

_

drf

( 12. 21)

 

Граничные условия (18) перерабатываются следующим обра-

том

 

 

drs

■-=0 при i(

0,

~(R

 

( 12.22)

d<s

 

-> 1 при 1]

cr.

~cht

 

 

Результаты численного решения этого уравнения приводятся в специальной литературе*.

§ <1. РАСПРОСТРАНЕНИЕ СТРУИ

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

Плоская ламинарная струя. Из бесконечно тонкой щели, рас­ положенной по оси z, в направлении оси х истекает струя с беско­ нечно большой скоростью, но конечным импульсом /о (рис. 12.1).

* См., например, Л. Г. Лойцянский. Ламинарный пограничный слой. М..

ГИФМЛ, 1962.

139

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