Н.А. Слёзкин - Динамика вязкой несжимаемой жидкости (1159534), страница 65
Текст из файла (страница 65)
наложенного на прямолинейно-параллельное течение в круглой цилиндрической трубе: д Оф' дар», 1 дф»доы дзф» дпы дф»дтиы д! ы дл ! г дл дг дгда дг дл дгт Для случая основного кругового течения вязкой несжимаемой жидкости в уравнениях (2.10) необходимо положить: де,, диет о жО, о ямО, — вяб, — — О, дф * дл до, — '+ дг до дп', — 1+ дг о "дг;г дт ' ( т ге+ге дт)' де„.доы 1 д)г „ ыдл ' "дг р да+ д (го„) д (го,') ди, + — — '+ — ' =- О. дг де дт >стойчивость ллминленых течении — г+ -г+ — -"-"-+ --.- = О. дг г г дт дг Гели в рассматриваемом случае дополнительно предположить, что проекции вектора скорости поля возмущений и давление не будут зависеть от поляркото угла о до,.
до до, др дч ' дт ' дт ' де — '=О, — — ==О, — -'==О, — =О, то из (2.!5) получим следу>ошие дифференциальные уравнения симметричного поля возмущений, наложенного на основное круговое движение вязкой несжимаемой жидкости: д> г В д>' (, > ' гт! до, 1др дт рдд (го,.) д (го„) дг д" + — =".= О, ('2. 16) тдс д'-' ! д дз дгз ' гдг ' д'' Приближш>ные дифференциальные уравнении (2.9), (2.14) и (2.16) соответственных полей возмущений использовались отдельными автора>>и для исследования устойчивости основных ламннарных течений вязкой несжимаемой жидкости.
Теперь перейдзм к установлению некоторых энергетических соотношений для поля возмущений. Обратимся к полным уравнениям поля возму>ценив (2.5). Умножая первые три уравнения на и', о' н ш' соответственно и складывая, При этих предположениях дифференциальные уравнения поля возмущений представятся в ниле > > дт„о, до, до,„огт>,ч —,, ~бо' т ( '.) (2.15) де+ г дт яда а 21 овщие телвнзиня для возмюценного движения 395 получим: Используя уравнение несжимаемости и выражения компонент вектора вихря, будем иметгн г г и' ди'+ о' ао'+ щ' Лщ' = г .— 2 — (о в' —. юа',) + — (ю а' -- и'в') + + — (и'а' — о'а' )1 — 4 в', =-Й" (';: — ~- ")1+АИ'-а+ ~ ")1+Й-'6+1 "М Первые четыре слагаемых в левой части (2.17) в своеп совокупности представляют собой ииднвидуальнухз произ~годную по времени от кинетической энергии единицы массы в поле возмущений при условии, что переход частиц из одного положения в другое происходит со скоростью основного движения.
Для этой производиои введем отдельное обозначение д д д д де дь ' 'дх 'ду 'дх ' — '= — + и — +о — +щ —. (2.18) Для всех последующих слагаемых левой части (2.17) также введем отдельное обозначение ,а диг, ду,,а да,, удиг де 'г г)4 = — 1и' — +о' — — + щ' — '+ и о ~ — + — )+ дх ду д- )ду дх 7 ' и'.-"'(й'-.'+'д — ')+"'о'(' — "+'— ")) (2 ") Г/г (э д- — д —,— д — д —— 2 2 2 2, мдиг,едег,гдвг — + и, = + о — + тп — + и — + ед" — .+ щ' — + дг дх г ду г дх дх ду дх + (а'Ьи'+о'Ло'+щ'Ьщ'), (2,17) уСТоичньоь~~ вви Используя обозначения (2.18) и (2.19) и предшествующие формулы преобразования, получим из (2,17) следующее энергетическое соотношение для единицы объЕма частиц жидкости в поле возмущений: — '(Р—,) = рМ вЂ” 4рв'+ — (и'(р'+ —, р(>' )+29(о'в' — -ш'в,')~+ + — (о'(р'+-,— р)г' )+2р(ш'в' — и'в'„))+ + — ~ш'(р'+ — р)г' )+ 29(и'в,' — и'в',)~.
(2.20) Для случая плоского поля возму>ценил будем иметь: ь„э — '( — ) = р Ч вЂ” 4йв' -à — ~и'(р' -1--,; э(г' )-4 21>п'в'~+. г> ~ при этом из (2.19) получим: ,э ди> с ди... Гди> ди,т) М = — (и'э — + о' — -(- иги ( — + — )(. дх ду )ду дх ) (2.21) (2.22) Если >ке основное течение будет плоским прямолинейно-параллельным, т. е.
д»> ю,—= . О, д — '=О, и, =- и,(у), то вь>ражение для М примет вид , ии> М =- — и'о —. ду (2.23) Допустим, что мо>хно выбрать такой конечный прямоугольник с площалью 5, на границах которого у = с, у = с( проекции вектора скорости поля возмущений обращаются в нуль, а на других границах х = а, х = д распределение скоростей, давлений и вихрей будет олинаковым. При этих условиях при интегрировании по >лощади рассматриваемого прямоугольника будем нметгн ~ — (и'(р'+ — р)>' )+ 2ро'в'~ г(хг(у = ~ ~и (Р + -2 р"" ) г-2ро™~ >(у = О, з=е ш'=— О, в'= — О, в' =О, в'=.в' х ' и и, следовательно, энергетическое соотношение (2,20) будет предсгавляться в виде й 2) озщиа уиляиання для аозмущйнного даижаиия 397 ~ ~ д [о' (р'+ — 71' ) — 2 и' '~дхду = и:.ь ~о'(р'+ —, р)l' ) — 2ри'м'~ дх =- О, и=а — — рР' с(х ду = р ~ ) М Нх игу — 4)ь ~ ) м' с(х игу.
(2.24) Полученное интегральное соотношение (2.24) представляет собой энергетическое соотношение для частиц жидкости внутри указанного прямоугольника в поле возмущений. Это соотношение показывает, что возрастание кинетической энергии поля возмущений может происходить только тогда, когда величина М будет заведомо положительной и прн этом такой, чтобы значение первого интеграла в правой части (2.24) превосходило значение второго интеграла. г4ля случая плоского прямолинейного течения, для которого Щ)О, дй величина М из (2.23) может быть полозкительной, если проекции вектора скорости поля возмущений будуг иметь разные знаки, т.
е и'о'( О. (2.25) В указанных в предшествующем параграфе статьях, в которых исследование устойчивости ламинарныл течений проводилось с помощью энергетического метода, в качестве допущения принималось, что возрастание со временем кинетической энергии поля возмущений может служить вполне достаточным признаком возникновения неустойчивости исследуемого ламинарного течения. Если принять это допущение, то дальнейшая задача исследования устойчивости прямолинейно-параллельного течения между параллельными стенками будет сводиться к подбору соответственного поля возмущений, удовлетворяющего неравенству (2,25), при котором правая часть равенства (2.24) обращалась бы в нуль и при этом число Рейнольдса исследуемого ламинарного течения принимало бы наименьшее значение. Приравнивая правую часть (2.24) к нулю и подставляя значение М из (2.23) н значение угловой скорости вихря, получим: ~и'и' — 'дхг(у+ ~ ~ ~ — — — ) улду=О. (2.26) .!'.
'' — ' диг ( Р Гди' ди'1г 6'у .),( (дх ду ) Если теперь перейти в равенстве (2.25) к безразмерн~м величинам, полагая х =- '— , у =- — "', и„= — (/((т), и' =. Ыр(Е г), о' = (/ф(Е ч)), Й = — „ устойчивость ллминлтных > знаний (гл. х> то для числа Рейнольлса получим следующее равенство: ) (ф — — ) А'дч К— (2.27) >" (ч] и (; —, ч) ф (-, ч) д( дч Таким образом, задача определения минимально~о значения числа Рейнольдса будет сводиться к опрелелению минимума отношения(2.27) лвух двойных интегралов.
ф 3. Исследование устойчивости ламииарного течения с прямолинейным профилем распределения скоростей Пусть мы имеем дяе параллельные стенки на расстоянии Ь лруг от друга. Если нижняя стенка будет неподвижной, а верхняя будет перемещаться параллельно самой себе со скоростью (У и если перепада давлений в направлении течения не будет, то для основного поля ламинарпого течения межлу параллельными стенками булем иметь прямолинейный профиль распределения скоростей по сечению, т. е. (г и,= ау, о>=0. (34) Вводя безразмерные независимые величины х у (У 0й — — — (х =- —, (3.3) представим дифференциальное уравнение (3.2) в зиле да(у дар — — + т> —.— = — Ь дб'. дт д! и (3.4) Далее, как уже указано в Я 1, функцию тока представим в виде >У = (Уйеын "'>у(г>).
(3.5) Тогла Ь)' = (>Лет>) -"П( — я-'У+ — 1= Еде'>1 -"Пр(т>), дт>л ) (3.6) Для исследования устойчивости данного ламинарного течения по л>етоду чалых колебаний мы должны обратиться к приближенному дифференциальному уравнению (2.9) для функции тока поля возмущений. Подставляя в зто уравнение выражение (3.1) для продольной скорости, получим с.чедующее дифференциальное уравнение с частными производными четззртого порядка: д а>)' (У д лр' де Л- дл ь 31 тзчяниз с пгямолинейным пгоеилам глспгадалания скогоствй 399 и дифференциальное уравнение (3.4) запишется: — я+ф Яг(ат> — >3) — Я") =. О.
ввт дг;-' (3.7) Таким образом, задача свелась к решению однородного дифференциального уравнения с обыкновенными производными второго порядка (3.7) и последующего решсння неоднородного обыкновенного дифференциального уравнения второго порядка — — — яву=у ггзэ' лф (3.8) Если ввести новое независимое комплексное переменное чз+ я(з — ач] (ч>1) Ь то дифференциальное уравнение (3.7) преобразуется: члч — +ау = О.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
