Н.А. Слёзкин - Динамика вязкой несжимаемой жидкости (1161662), страница 71
Текст из файла (страница 71)
.Рд ()д,„а) ддГд (й,„а) (5.1 1) представляет собой цилиндрическую функцию общего вида, Постоянные Ад и Аз подобраны так, чтобы были выполнены граничные условия для и,, т. е. Лд/д()д„,Ь)+ А Ддд (Й,„Ь) = О, Л,/,(й а) +Л.М,(й а) =О. (гл. кг гстойчивость ллминлгных течений где Н,„= ~ г3,*(ймг) тг. ь Подставляя (5.9) в первое уравнение (5.6), получим: ч ~ — -+ — — — ~ — + Л' ~ из~ — 2СЛ ~ а,2д(й,„г). (5.13) Гнзиз ! лаз / 1 лб м=з Решение уравнения (5.13) без правой части представляет собой цилиндрическую функцию Л„(Ыг) = А„),(Гл'г)+А Ид(ЬГг), где Аз и А,— произвольные постоянные. Решение уравнения (5.!3) с правой частью можно представить в виде ряда Х К,.А(д.,г) ю=а коэффициенты которого могут быть определены после подстановки этого ряда в уравнение (5.!2) в виде 2С,а,„ '(Л + Лг) (5.
14) Таким образом, для функции а. мы получим: из(г) =Азу (1Л'г)-(-А М„(1Л'г)+~ Ь„,2~(й,„г). (5,15) Так как функция (5.10) обращается в нуль на границах, то для удовлетворения граничных условий (5.8) для из постоянные Аз и А, необходимо положить равными нулю. Обратимся теперь к уравнению (5.6) для и . На основаник рекуррентных формул имеем: — „— ~е(йг) = — ~~ (йг). 1 а Уравнению для из без правой части К Гнал 1 низ г А+ ',' 1 — 0 йг~ дгз ' г Мг На основании теории рядов фурье — Бесселя коэффициенты а,„будут представляться в виде а ат = ~ гиа (г) Еа Ямг) Ыг, 1 (5.12) и„,, ь й 5) кгкговов движанив мвждт ввсконвчными цилиндглми 425 можно удовлетворить, полагая иа = Аз+ Ар/о(!1,'г)+ Агро(!Л'г).
Следовательно, общее решение второго уравнения(5.6) можно искать в виле "а(г) = Аз+Аз!о (!Лг)+ МЧо (!Л г)+ Х ау~~о (Льаг) (5.16) Подставляя (5.16), (5.!5) и (5.9) во второе уравнение (5.6), полу- чим следующее уравнение для определения постоянных а(„;. +(Л' +й')Е,((г «) = ~У «(Л' + й~)а Е,(й г)+ ааа 1 ж=а +2(Са+ — ~) ~~» ., а 2 (й,„г). (5.17) аа 1 Подставляя (5.!6) н (5.9) в третье уравнение (5.6) и учитывая рекуррентное соотношение ха(х) хо(х) а Еа(х) получим: ~ (л а +ЛИ )Лр(й г)+Л(Аь+АеУо(!Лг)+А М((Лг)) =О.
(518) на=а К уравнениям (5.!7) и (5.18) необходимо присоединить два уравнения, которые мы можем получить, удовлетворяя выражением (5.16) граничным условиям (5.8) для ие(г). Дальнейшие вычисления, проведенные в работе Тэйлора '), приводят к бесконечной однородной системе уравнений для постоянных Аа, А„ и а . Приравнивая нулю определитель этой системы уравнений, получим характеристичаское или вековое уравнение, связывающее величины р и Л'=! ха+ — с заданными параметрами аадачи вы на, (а и а. Подробный анализ этого уравнения проводится в цитированной работе Тэйлора в предположении, что разность а+а радиусов цилиндров а — (а мала по сравнению с полусуммой — .
2 На основании этого анализа получается, что если цилиндры вращаются в одну сторону, то круговое движение жидкостн будет всегда устойчивым при выполнении следующего неравенства; еаааааа ( наля. (5.19) а) Т а у!ос О., Ргос. йоу. 8ос. (А), т. 223, !923, 426 тстойчнвОсть ллминлвных твчвний (гл. х~ В работе Сайнджа ') покавано, что критерий устойчивости (5.19) можно доказать и не прибегая к предположению о малости разности радиусов цилиндров по сравнению с их полусуммой.
Если же неравенство (5.19) не выполняется, т, е, если аг — > —. ее Ьв (5. 20) или если концентрические цилиндры вращаются в разные стороны, то круговое движение частиц вязкой жидкости теряет свою устойчивость, как только число Рейнольдса (относящееся, например, к внешлз,„,1 нему цилиндру, т. е. К =. — 1 превысит своа критическое значение, достаточно близкое к значению, установленному из экспериментов.
-глгг -л;гг -lлл -гйт -лгг и ьр тыл гьп лпт Рис. 102. г) 3 у и й е 1., Ргос. йоу, Бес. (А), 1.опбоп, г. 167, 19Ж На рис. 102 представлена кривая, отделаюшая область устойчивости кругового движения от области неустойчивости, подсчитанная а работе Тэйлора для случая, когда Ь = 3,55 сж н а = 4,035 сж, так что ()= а 1е — 11 = 1,292.
По горизонтальной оси отложены значения — "', а по вер- Ь) тикальной — —. Есля внешний цилиндр находится в покое, а внутренний вращается или оба цилиндра аращаются з одну сторону, то появление неустойчивости характеризуется образованием ряда вихрей в плоскости меридионального сечения, заполняющих вса пространство между поверхностями цилиндров, прн этом направления вращений этих вихрей чередуются.
Такое образование вихрей хорошо подтверждается экспериментально. Окрашенная жидкость, первоначально распределенная 6 61 ов истойчивости лвижания взвешенной частицы ч27 в виде тонкого слоя по поверхности внутреннего цилиндра, впоследствии свзртывается в виде чередующихся колец, охватывающих центры вихрей. Если цилиндры вращаются в разные стороны, то неустойчивост>, кругового движения жидкости проявляется в образовании двух рядов вихрей, из которых один.
имеющий большую интенсивность, располагается вблизи внутреннего цилиндра, а второй — с меньшей интенсивностью — вблизи внешнего. В опытах окрашенная жидкость собиралась только вокруг вихрей с большей интенсивностью, а в области, примыкающей к внешнему цилиндру, вода оставалась прозрачноя. Места расположения центров вихрей, установленные на основании вычислений, очень хорошо подтверждены опытами; отклонения вычисленных значений критического числа Рейнольдса от соответственных опытных значений не превышают 2",' . В 6. Об устойчивости движения взвешенной частицы в ламинарном потоке Первые экспериментальные исследования Людвига и Гагена, послужившие оснояаниеи для постановки вопроса об устойчиности и неустойчивости ламинарного движения в цилиндрической трубке, были проведены с помощью наблюдений над движением примешанных к воде видимых частиц.
Такими частицами в опытах Гагена были опилки темного янтаря. В экспериментальных исследованиях Репнольдса наблюдения проводились за поведением тонкой окрашенной струйки, вводимой в поток прозрачной жидкости. Таким образом, при экспериментальных исследованиях производились лишь местные возмущения исследуемого ламинарного течении.
а при теоретических исследованиях, рассмотренных в предылущих параграфах, возмущения накладывались на все течение в целом. В этом и заключаются расхождения между подходами а теории и в экспериментах. Поэтому всякая попытка приблизить подход в теории по вопросу об устойчивости ламинарного течения жидкости к тому подходу, который использовался в ряде зкспериментоа, может представлять известный интерес. Рассмотрим движение частицы, авелэнноя каким-либо образом в поток вязкой жидкости. Чтобы составить дифференциальные уравнения движения такая частицы, необходимо учесть, по возможности, все основные силы воздействия на частицу со стороны окружающая жидкости, находящейся в движении. Будем эти силы относить к единице объема оассматриваемой частицы.
К основным силам следует отнести, во-первых, силу веса эа вычетом статической силь> Архимеда Р = — К1р — р)у (6Л>) тле р,— плотность частицы, а у — единичный вектор вертикальной оси у. Во-вторых, силу сопротивления, пропорциональную в первом гстойчизость ллмиилгных твчвний (гл. хг приближении разности скоростей частицы и местной скорости потока ), = — д,()г — и), (6.2) где )г — вектор скорости частицы. е/ — вектор местной скорости потока и Й вЂ” коэффициент сопротивлении, который для шаровой частицы ралиуса а будет равен й = — —. Я и 2 лз' (6.3) Помимо этих двух сил, необходимо учесть и боковую силу.
Зело в том, что благодаря непрерывному распределению вихрей в ламмнарном потоке относительное обтекание потоком частицы будет всегда циркуляционным и, следовательно, всегда будет возникать боковая сила, пропорциональная по теореме Н. Е. Жуковского относительной скорости набегающего потока н циркуляции. При обтекании плоско- параллельным потоком круглого цилиндра длины ! полная подъемная сила представляется в виде рГ) ), а подъемная сила, приходящаяся на единицу объема цилиндра, будет равна р)г —, Г Таким образом, боковую силу, приходящуюся на единицу объема частицы, можно представить в виде Р,=д.,р(() — ~)Х ((), (6.4) где й,— безразмерный коэффициент, с помощью которого можно приближенно учесть необходимые поправки на те допущения, которые связаны с распространением теоремы Жуковского, применимой к плоско-параллельному обтеканию бесконечного цилиндра, на рассматриваемый случай пространственного обтекания частицы произвольной формы.
Если к тому же учесть, что формула Жуковского о подъемной силе оправдывается экспериментально при малых углах где à — циркуляция, У в скорость потока на достаточном удалении и 5 в площадь сечения цилиндра плоскостью, перпендикулярной к его образующей. При применении этой теоремы Жуковского к обтеканию частицы мы принимаем за скорость набегающего потока разность скоростей Гà — У, за вектор присоединенного вихря — вектор вихря рассматриваемого потока го) сГ, за Ья — сечение частицы той плоскостью, которая перпендикулярна к плоскости векторов () — гг и го)К и, наконец, за величину циркуляпии по контуру сечения Ь5 мы можем принять согласно теореме Стокса поток вектора вихря через ЬЗ, т.