книги из ГПНТБ / Палагин, Э. Г. Основы гидромеханики учебное пособие для метеорологов
.pdfГ Л А В Л 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) |
дх |
ду |
|
dх |
ду |
||
|
|
|||||
|
ди |
|
д~о |
о. |
|
(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 — |
Л-,- |
|
|
||||
дус, |
-п— при |
|
|
|
|||||
|
1о |
|
|
|
|
|
|
||
В (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т |
Г |
— [? |
Ы ~ •' |
|
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