
книги из ГПНТБ / Картвелишвили Н.А. Потоки в недеформируемых руслах
.pdfскорость подъема пузырьков, заполняющих с одинаковой концен трацией s пространство между двумя горизонтальными плоскостями, определяется по формуле
|
|
ш = Ш о ( 1 - 5 ) , |
|
|
(18.15) |
|
где «о — скорость всплывания |
(гидравлическая |
крупность) |
уеди |
|||
ненного пузырька. |
|
|
|
|
|
|
В качестве |
гипотезы |
V I I |
примем |
вместо |
выражения |
(18.15) |
формулу |
|
|
|
|
|
|
|
|
ш=ш 0 (1 — S). |
|
|
||
Теперь можно написать уравнение работ |
|
|
||||
им |
р j ( l — s) |
и dF-\-ps j |
stCdF\=t0yJU |
|
||
— |
\ F |
|
|
|
|
|
или
1— aS |
R 7c 1 (1 — aS) U ' |
v |
' |
где R — гидравлический радиус, a
?(l-S) |
+ |
9sS |
|
|
g |
|
|
объемный вес водовоздушной |
смеси. Для однофазного |
потока |
|
г 0 |
f/2 |
|
(18.17) |
|
С2 |
|
|
где С — коэффициент Шези. Примем, |
что формула (18.17) |
сохра |
няет силу и для аэрированного потока. В этом и будет заключаться гипотеза V I I I . Основанием для этой гипотезы служит то, что в край
них случаях 5 = 0 (поток воды) |
и 5 = 1 |
(поток |
воздуха) |
формула |
|||
(18.17) верна, а значит |
можно |
ожидать, |
что |
она будет верна и |
|||
в промежуточных случаях. Из выражений |
(18.12), |
(18.16) |
и (18.17) |
||||
имеем |
|
|
|
|
|
|
|
; * = - Й | |
{ 0 - a S ) W |
+ |
И |
- 5 |
) } • |
08.18) |
Гипотеза IX заключается в том, что формулу (18.18) можно при менять и для неустановившегося движения. Это есть обобщение обычной для гидравлики гипотезы — сопротивления в неустановив шемся движении могут определяться по тем же формулам, что и в установившемся. В результате всех этих преобразований членов
уравнения (18.11) ему можно придать такой вид:
^ - в 5 > - ^ ~ т ^ и о - |
- 5 ) ] + |
+ ( l _ » 2 5 ) 4 - ^ + « - ^ - J + A |
(18.19) |
|
где |
|
|
a = a' + a" - & 3 S, p ^ - ^ S , |
Л = 1 , 2, 3, 4. |
|
Складывая уравнения неразрывности (18.6) и (18.8) и учитывая, что Q + q = UF, получаем (при W — Ws = 0)
|
|
|
d U F |
• d |
F |
-0. |
|
|
(18.20) |
|
|
|
дх |
' |
dt |
|
|
|
|
Используя |
выражение |
(18.20), |
можно |
придать |
уравнению |
||||
(18.19) окончательный вид |
|
|
|
|
|
|
|||
' 0 - e 5 |
) - ^ - = T ^ [ 0 - ^ ) w |
+ ^ 5 |
( 1 - 5 ) ] + |
||||||
, 1 — hS |
dU |
, |
£/ |
dt/ |
a — 1 + |
ft2S |
t7 |
, |
|
1 |
a^ |
^ ~ |
g |
dx |
|
g* |
|
F |
dt ' |
+ |
(18.2!) |
В частном случае неаэрированного потока, когда 5 = 0, это урав нение переходит в уравнение (13.25), если в последнем отбросить член, учитывающий поверхностное натяжение.
Вернемся теперь к допущениям, которые были приняты при вы воде гидравлических уравнений неустановившегося потока.
В случае однофазного (не аэрированного) потока вывод уравне ния (13.25) не требует, в сущности говоря, никаких предположений, кроме тех, которые вытекают из использования одномерной идеа лизации течения, т. е. являются совершенно обычными для всей гидравлики вообще. Останавливаться на этих предположениях нет необходимости. Нужно лишь подчеркнуть еще раз, что принятый способ вывода гидравлических уравнений из общих уравнений гид родинамики позволяет свести к минимуму число предположений та кого рода, в частности, он делает ненужными неоправданное пред положение о гидростатическом распределении давлений в потоке
9* |
131 |
и весьма искусственные построения, обычно применяемые при учете влияния кривизны струй в уравнениях движения.
Значительно сложнее обстоит дело в случае двухфазного |
(аэри |
||||
рованного) |
потока, ибо для замыкания уравнений Франкля, |
из ко |
|||
торых мы |
в этом |
случае исходили, |
требуется такая |
информация |
|
о взаимодействии |
фаз между собой и с граничными |
поверхностями |
|||
потока, которой |
наука в данное |
время не располагает. В |
связи |
с этим полезно дать сводку гипотез, принятых при выводе гидрав лических уравнений аэрированного потока, с некоторыми дополни тельными ком ментар 11 ям и.
Гипотеза I — поток имеет четко выраженную свободную поверх ность. В рамках одномерной идеализации эта гипотеза совершенно необходима, так как случай, когда она не выполняется, т. е. когда движение воды с воздушными пузырьками в нижних слоях потока плавно переходит в его верхних слоях в движение водяных капель в воздухе, принципиально нельзя исследовать, отвлекаясь от стра тификации фаз по вертикали. Иными словами, этот последний слу чай требует по меньшей мере двухмерной идеализации.
Гипотеза I I — вязкость внутри жидкой и газообразной фаз и на поверхностях их контакта во внимание не принимается. В гидрав лической идеализации от этой гипотезы можно отказаться за счет соответствующего усложнения выкладок. Но делать это нецелесо образно, ибо те тонкие рассуждения, к которым тогда придется прибегнуть, не дадут ничего нового: дифференцированный анализ вязких н турбулентных напряжений все равно бесполезен, так как он будет сведен на нет необходимостью принять достаточно грубую
гипотезу о гидравлических сопротивлениях |
(гипотеза V I I I ) , |
о кото |
рой подробнее будет сказано ниже. |
А в уравнении |
|
Гипотеза III — пренебрежение членом |
(18.12). |
Эта гипотеза основана на том, что А есть квадратичный член по от ношению к вертикальным составляющим скоростей течения. А так как эти составляющие в десятки раз меньше горизонтальных, то их квадратами вполне можно пренебречь при изучении явлений, опре
деляющихся в своих главных чертах |
горизонтальными составляю |
||
щими. |
|
|
|
Гипотеза IV — отсутствует обмен |
воздухом |
и водяной пылью |
|
между потоком и атмосферой. Эта гипотеза вынужденная, так |
как |
||
о процессе этого обмена сейчас ничего неизвестно, по крайней |
мере |
||
о его количественной стороне. Сфера |
действия |
гипотезы IV очень |
ограничена. В задаче об устойчивости равномерного потока, к кото рой и применяются (в данной монографии) полученные с ее помо щью уравнения, она не приводит к противоречиям. Но, например, при исследовании замедленного (в результате подпора) движения аэрированного потока, из которого обязательно выделяется воздух, эта гипотеза неверна даже в первом приближении.
Гипотеза V — средняя концентрация воздуха в потоке постоянна. Эта гипотеза может быть (но не обязательно будет) справедливой только в тех условиях, в которых справедлива гипотеза IV. Для того чтобы отказаться от нее хотя бы в рамках действия гипо-
тезы IV, нужен очень подробный гидродинамический анализ двух фазного потока, а этот анализ еще не выполнен: он невозможен, пока не будет установлена достоверная связь между скоростями движения фаз.
Гипотеза V I — количество движения, |
обусловленное |
пульсаци |
|
ями скоростей, связано с количеством |
движения, |
обусловленным |
|
осредненными скоростями, такой же зависимостью, |
что |
и в одно |
фазном потоке. Эта гипотеза так же, как и соответствующее соот
ношение для однофазных потоков, |
по-видимому, |
не вызывает |
осо |
бых возражений. |
|
|
|
Гипотеза VII—гидравлическая |
крупность системы пузырьков |
||
в потоке связана с гидравлической |
крупностью |
уединенного |
пу |
зырька через среднее, а не через локальное воздухосодержание. Не обходимость в такой гипотезе вытекает из одномерной идеализации, т. е. из того, что мы отвлекаемся от стратификации воздушных пу зырьков по вертикали.
Гипотеза V I I I — связь между средними касательными напряже ниями на дне потока, удельным весом водовоздушной смеси, сред ней скоростью ее течения и коэффициентом Шези такая же, как и в случае однофазного потока. Опытов, которые прямо подтверж дали бы эту гипотезу, нет. Основание для ее принятия, указанное выше и состоящее в том, что эта гипотеза безусловно верна в двух крайних случаях (поток воды и поток воздуха), носит необходимый, но не достаточный характер: если бы рассматриваемая гипотеза не была верна в этих случаях, ее следовало бы отвергнуть, но из того, что она верна в этих случаях, еще не следует, что она верна вообще. Заменять ее каким-либо более сложным соотношением сейчас нет никаких оснований, ибо предложить что-либо более обоснованное при нынешнем состоянии изученности сопротивлений в двухфазных потоках просто нельзя. Вопрос о сопротивлениях до сих пор далеко не решен даже для однофазного потока, ибо рекомендации класси ческой гидравлики в этой области внутренне противоречивы, о чем уже говорилось в § 16. Если такое положение справедливо для од нофазных потоков, то ждать каких-либо теоретически или экспери ментально обоснованных результатов для много более сложного случая двухфазного потока, разумеется, сейчас не приходится.
Гипотеза IX — законы сопротивлений в установившемся и неус тановившемся движениях одинаковы. Эта гипотеза заимствована из однофазной гидравлики, где она принимается повсеместно, при отсутствии каких-либо экспериментальных данных о сопротивле ниях в неустановившемся потоке. Но эта гипотеза стоит под сомне нием, о чем также говорилось в предыдущей главе.
Как видим, попытка анализа гипотез IV, V и V I I , а также в осо бенности гипотез V I I I и IX поднимает две проблемы принципиаль ного значения: проблему гидродинамики аэрированного потока и проблему гидравлических сопротивлений, каждая из которых тре бует самостоятельных специальных исследований.
Глава VI
УС Т О Й Ч И В О С Т Ь
ТЕ Ч Е Н И Й
19. Общие определения и критерии устойчивости
Задачи, связанные с устойчивостью различных движений, и в ча стности течений жидкости, играют огромную роль в современной науке и технике. Сущность этих задач сводится к следующему. Пусть некоторая система (машина, текущая жидкость или вообще любое естественное образование или искусственное устройство, описываемое некоторыми функциональными уравнениями) нахо дится в состоянии движения А І . Дадим этому движению то или иное изменение, например (если система описывается дифференциаль ными уравнениями) изменим начальные условия. Такое изменение называют возмущением, возникшее в результате возмущения, дру гое движение Лг называется возмущенным движением, а прежнее движение Аі — невозмущенным движением. Если различие между движениями А І п An с течением времени становится и дальше оста ется малым в некотором (точно определяемом в конкретных слу чаях) смысле, то движение А І называется устойчивым по отноше нию к данному возмущению. Если же различие между движениями А І и Ао С течением времени не только становится малым, но и стре мится к нулю, то движение А І называется асимптотически устойчи вым по отношению к данному возмущению.
Если бы уравнения движения тех систем, с которыми приходится сталкиваться на практике, допускали точные решения в квадрату рах, то теория устойчивости не могла бы возникнуть: все свойства движения, включая устойчивость или неустойчивость, отражались бы этими решениями и в какой-то специальной теории попросту не было бы необходимости. Потребность в теории, позволяющей оце нивать устойчивость движений, не решая уравнений, описывающих эти движения, возникла именно потому, что точные решения прин ципиально возможны лишь в простейших и практически очень мало интересных случаях, а обычные приближенные решения для этой цели не годятся, ибо речь идет о поведении системы в течении не ограниченного времени.
Обратимся к точным определениям понятия устойчивости. Пусть движения некоторой динамической системы описываются функцио
нальным уравнением |
|
U7(u, г, f)=0, |
(19.1) |
в котором |
г — вектор |
пространственных координат; и — вектор не |
известных |
функций; |
W — функциональный векторный оператор. |
Когда речь идет о движении жидкости, последнее может пони маться как актуальное или как усредненное турбулентное, в пер вом случае и — вектор действительных компонентов скорости и да вления, а во втором — усредненных. Поэтому дальнейшие опреде ления будут относиться либо к актуальному, либо к усредненному движению.
Определим каким-нибудь способом норму ||u|j вектора и, напри мер,
Здесь N — число компонент Ui вектора и. Способ введения нормы в данном случае не играет роли, важно только, чтобы были выпол нены требования, предъявляемые к норме: норма должна быть не
отрицательна, инвариантна по отношению к нумерации |
компонент, |
|||||||
обращаться |
в нуль |
тогда |
и только тогда, |
когда все |
компоненты |
|||
равны нулю; |
норма |
суммы двух векторов |
должна |
быть |
не более |
|||
суммы их норм. |
|
|
|
|
|
|
|
|
Пусть u* (г, t) —решение системы |
(19.1), устойчивость |
которого |
||||||
исследуется, называемое |
невозмущенным движением, |
a u** (г, t) — |
||||||
некоторое другое решение системы (19.1), |
называемое |
возмущен |
||||||
ным движением. |
|
|
|
|
|
|
|
|
Если существуют такие два положительных числа |
є и б = 6 (є), |
|||||||
что из неравенства |
|
|
|
|
|
|
|
|
|
|
||ц**(г, *„)-«*(г, |
^ 0 )||< 8 |
|
|
(19.2) |
||
вытекает для всех t>ton |
для всех г в области определения |
решений |
||||||
системы (19.1), неравенство |
|
|
|
|
|
|||
|
|
||и**(г, 0 - и * ( г , |
О Н О , |
|
|
(19.3) |
то система устойчива в смысле Ляпунова. Если же при выполнении неравенства (19.2)
Ит 1 u**(г, 0 - и * ( г , t)\\=Q,
то система асимптотически устойчива в смысле Ляпунова. В даль нейшем устойчивость в смысле Ляпунова будем называть устойчи востью в L смысле, или L-устойчивостью. Этот вид устойчивости
имеет колоссальное практическое и научное |
значение. Но вместе |
с тем L-устойчивость не всегда определяет |
действительную устой |
чивость системы. Можно указать случай, очень важный для гидро логии и гидротехники, когда L-устойчивая система оказывается не устойчивой в ином, очень существенном смысле. Речь идет о само произвольном возникновении бегущих волн в открытых руслах большого уклона. Из многочисленных наблюдений известно, что та
кие волны, возникающие от случайных |
малых возмущений, имеют |
в каждом сечении потока ограниченную |
амплитуду, т. е. для них |
выполняется условие (19.3). Обозначим через єо(г) точную верх нюю границу є при данном г. Эта граница, вообще говоря, зависит от г и, в частности, она может увеличиваться вдоль потока, т. е. при перемещении бегущих волн вниз по течению их амплитуды могут возрастать. Практика обычно выдвигает требование, чтобы такого возрастания не было, т. е. чтобы устойчивость была равномерной по длине потока. В общем случае (потоки в двумерной и трехмерной идеализациях) требование невозрастания волн при их перемещении не сводится к требованию равномерной /.-устойчивости. Поэтому наряду с L-устойчивостыо рассмотрим другой вид устойчивости. Условие
|u**(r, t)~u*(r, |
0 1 = max |
|
(19.4) |
|
г |
|
|
определяет функцию ^ = 0(г) или семейство таких |
функций. Если |
||
для каждой функции этого семейства |
можно указать |
ограниченную |
|
положительную функцию / (г) такую, что из неравенств |
|
||
||u**(r, g - u * ( r , |
* 0 ) I < / ( r ) , |
|
(19.5) |
^ < Є ( г ' ) < Є > ( г " ) |
|
(19.6) |
|
следует, что в некоторой подобласти области определения |
решений |
||
системы (19.1) имеет место неравенство |
|
|
|
«и** [г", Є ( г " ) ) - и * ( г " , G(r")}||<||u** (г', 0 ( r ' ) } - u * { r - ' , |
Є(г'))||, |
||
|
|
|
(19.7) |
то система устойчива в этой подобласти. Этот вид устойчивости бу дем называть устойчивостью в 5 смысле, или 5-устойчивостью.
Если рассматривать возникновение |
бегущих воли |
в призмати |
||||||
ческом |
русле, |
невозмущениое |
движение |
в котором равномерно, |
то |
|||
5-устойчивость |
(для |
этого случая) будет не |
только |
равномерной |
||||
по длине L-устойчивостью, но и более того, |
5-устойчивость в си |
|||||||
стеме |
координат, движущейся |
вдоль потока |
со скоростью волны, |
|||||
будет |
тождественна |
обычной |
L-устойчивости |
в неподвижной |
си |
|||
стеме координат. Однако из всего этого нельзя делать |
вывода, |
что |
||||||
L - и 5-устойчивости |
совпадают. |
|
|
|
|
Исследование устойчивости течений сейчас реально только по отношению к возмущениям, дающим бесконечно малые изменения давлений, скоростей и ускорений. В этом случае задача очень упро щается тем, что уравнения возмущенного движения можно линеари зовать. Устойчивость по отношению к таким возмущениям назы вают устойчивостью в линейном приближении или устойчивостью' в малом. Для систем, описываемых обыкновенными дифференци альными уравнениями, доказано Ляпуновым [38], а для систем,, описываемых дифференциальными уравнениями в частных произ водных,— Ю. И. Неймарком [46], что из L-устойчивости в малом вытекает /.-устойчивость в большом, т. е. если обеспечена /--устой чивость в малом, то существуют такие конечные возмущения, по от ношению к которым /--устойчивость сохраняется. Доказательство'
аналогичного утверждения для б'-устончивости до настоящего вре мени не получено, но это не умаляет целесообразности исследова ний S-устойчивости в малом: пока не известны факты, свидетельст вующие о том, чтобы выполнение условий S-устойчивости в малом было недостаточно для обеспечения S-устойчивости течения по от ношению к любым наблюдаемым в действительности возмущениям.
Рассмотрим более подробно устойчивость в малом призматиче
ского потока. Определение |
призматического потока было дано |
|
в § 15. Здесь мы несколько |
обобщим это |
определение, включив |
в класс призматических также некоторые |
потоки в одномерной и |
двумерной идеализации. Мы будем называть призматическим лю
бой поток, описываемый |
одномерными, |
двумерными |
(например, |
|||||||||
в осесимметричной задаче) или трехмерными |
уравнениями, |
если |
||||||||||
коэффициенты уравнений |
возмущенного |
движения |
не зависят |
от / |
||||||||
и от х. Уравнение |
(19.1), например для трехмерного случая, |
может |
||||||||||
быть записано в координатной форме |
|
|
|
|
|
|
|
|
||||
|
Wt(u, |
v, w, р, |
х, Ь, г, ()=0, |
i=l, |
2, |
3, |
4. |
|
(19.8) |
|||
Здесь х, f>, г — цилиндрические координаты; |
и, v, w — проекции |
|||||||||||
вектора скорости на координатные направления х, |
г; первые три |
|||||||||||
уравнения (i = l , 2, 3) |
есть проекции |
векторного |
уравнения |
движе |
||||||||
ния на координатные |
направления; последнее уравнение — уравне |
|||||||||||
ние неразрывности. Положим в (19.8) |
|
|
|
|
|
|
|
|
||||
|
и=и0-{-Аи, |
v=vQ-\~Av, w=w0-\-hw, |
р=р0-{-Ар, |
|
(19.9) |
|||||||
где «о, "Ой,дао,ро — значения и, v, w, |
р |
в невозмущенном движении; |
||||||||||
Аи, Av, Aw, Ар — возмущения. Подставив значения |
(19.9) |
в урав |
||||||||||
нения |
(19.8), считая |
невозмущенное |
движение |
установившимся, |
||||||||
раскладывая левые части |
уравнений |
(19.8) |
в ряд Тейлора |
по Аи, |
||||||||
Av, Aw, Ар и удерживая только низшие степени возмущений, |
полу |
|||||||||||
чаем |
(учитывая, |
что |
Wi — некоторые |
дифференциальные |
опера |
|||||||
торы) |
систему однородных дифференциальных |
уравнений |
в част |
|||||||||
ных производных |
по х, •б-, г, t относительно |
Ди, Av, |
Aw, Ар |
с коэф |
фициентами, которые могут зависеть от т> и г, но не зависят ни от х,
ни от t. |
Эта система может быть сведена к |
одному однородному |
уравнению, частные решения которого суть |
|
|
|
ep' + f*<]>(&)W(r). |
(19.10) |
Здесь |
р и q — комплексные постоянные; |
Ф (•&) и ^¥ (г) — функ |
ции, зависящие соответственно только от •& и только от г. Эти функ ции ограничены, ибо в противном случае решение (19.10) не имеет физического смысла. Если ось х помещается внутри потока, то Ф(Ф), очевидно, периодическая функция периода 2я: в противном случае не будет однозначного соответствия значений Ф(ф) точкам сечения потока, нормального к оси х.
Связь между р, q, Ф, W получается при подстановке частного решения (19.10) в уравнение, из которого оно получено. Это приво дит к одному соотношению, связывающему р и q, следовательно, одно из этих чисел выражается через другое.
Решение (19.10) можно записать так: |
|
|
|
ехр (л- Re q-\-t Reр) |
[cos (л' Im |
q-\-t\mp)-\- |
|
+Y~^\s\n(x\mq |
+ t\mp)] |
Ф(&)іР(г). |
(19.11) |
Экстремум выражения (19.11) имеет место при |
|
||
x\mq + t\mp=k-%-, |
k=0, |
1, 2, . . . |
(19.12) |
Иными словами, уравнение (19.12) есть в данном случае урав нение семейства функций t = © (г), получаемых из условия (19.4). Скорость с перемещения экстремума (19.11), т. е. гребня волны воз мущения есть
с= — Im р\\т q.
Подставляя t из (19.12) в (19.11), получаем зависимость экс тремума волны возмущения от координаты х
е х Р [kJT • Щ ехр {(Re q - Rep Jgft) х\ • Ф (в) W (г).
При о 0 волна перемещается по оси х (прямая |
волна) |
и будет |
затухать, если |
|
|
R e < 7 - R e / > - ] ^ - < 0 . |
|
(19.13) |
Если рассматривать движение не в нашей неподвижной |
системе |
|
координат, а в системе, движущейся со скоростью с |
в направлении |
нашей первоначальной оси х, то 5-устойчивость сведется к L-устой- чивости, о чем уже упоминалось выше. Но с переходом к устойчиво сти в большом, а в рамках устойчивости в малом с переходом к не призматическим потокам возможность такого сведения исключа ется, ибо во втором случае скорость волны заранее не известна, нее определение требует решения соответствующих уравнений, а в пер
вом случае, сверх того, эта скорость зависит |
еще |
от |
неизвестных |
||
функций (например, от глубины |
потока). Поэтому |
указанная |
воз |
||
можность имеет ограниченное значение. |
оси х |
|
|
|
|
При с < 0 волна перемещается |
навстречу |
(обратная |
вол |
||
на) и будет затухать, если |
|
|
|
|
|
R e < 7 - R e / > - ] ^ - > 0 . |
|
|
(19.14) |
||
Представляя q и р в тригонометрической форме |
|
|
|
||
q — p (cos «р + У — 1 sin ср), |
р=а (cos ф + У — 1 |
sin <!>), |
|
||
напишем критерии затухания (19.13) и (19.14) так: |
|
|
|
||
coscp - sinc? - ^| - < 0 |
|
|
(19.15) |
при sin ф/sin ф < 0 , т. е. для прямой волны, и |
|
|
|
c o s c ? - s i n < p - ^ > 0 |
|
(19.16) |
|
при sin <p/sin чр>0, т. е. для обратной |
волны. Критерии |
(19.15) и |
|
(19.16) можно объединить и записать в виде |
|
|
|
sin (4> — <?)= — s sin |
|
(19.17) |
|
где є — произвольное положительное |
число для |
прямой |
волны и |
произвольное отрицательное число для обратной |
волны. |
|
-IX. |
О |
ті |
2п |
Зп to |
|
I |
|
II |
|
Рис. 21. Области устойчивости и неустойчивости по 5-
крптерига |
для призматических потоков. |
/ — прямые волны; |
/ / —обратные волны. / — затухающие; |
|
2 — нарастающие. |
Области существования, затухания и возрастания прямых и об
ратных волн, согласно критерию |
(19.17), представлены на рис. 21. |
На этом рисунке видно, что если |
некоторым ф и г|з соответствует |
прямая волна, то, сохранив гр неизменным и заменив ф на ф + я , по лучаем обратную волну, затухающую, если затухающей была пря мая волна, и нарастающую в противоположном случае. Условие