Структура сжимаемых вихревых течений Куэтта-Тэйлора Текст научной статьи по специальности «Математика»

УДК 531.1

М. И. Липатов, 3. С. До

Московский физико-технический институт (государственный университет)

Структура сжимаемых вихревых течений Куэтта^Тэйлора

Исследованы течения Куэтта-Тэйлора сжимаемого газа при больших числах Рейнольдса. На основе асимптотического анализа при больших числах Рейнольдса сформулированы математические модели и определены параметры подобия. Построена диаграмма характерных режимов течения при разных длинах волн возмущенного течения.

Ключевые слова: течения Куэтта-Тэйлора, влияние сжимаемости, структура течения.

1. Введение

Исследованию течений несжимаемой жидкости посвящено много работ. В то же время с точки зрения технологических приложений и фундаментальной науки анализ течений сжимаемого газа имеет несомненный интерес. Обычно важным методом исследований являются численные методы анализа. Был опубликован ряд статей [1], в которых представлены результаты решения уравнений Навье-Стокса для сжимаемого газа.

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

На практике могут реализовываться и течения с конечными числами Рейнольдса. Течение в таком случае описывается обычными уравнениями Навье-Стокса, нелинейными или линейными, в зависимости от рассматриваемых задач (описания среднего течения или исследовании устойчивости). Понятно, что для решения таких задач единственным способом является применение численных методов, основанных на использовании конечных разностей или представлениями в виде разложений в ряды.

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

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

При анализе течений с большими числами Рейнольдса возникает вопрос о сосуществовании инерционного и вязкого течений. Под вязким понимается течение, описывающееся только диссипативной частью уравнений Навье-Стокса. На самом деле даже без предположения об отсутствии инерционности, только принимая условие одномерности течения, можно получить описание с диссипативными членами (как в случае несжимаемой жидкости [3-4].

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

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

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

Более интересными оказываются течения, в которых действительно сказывается влияние вязкости, а именно в таких, для которых размер вихря в направлении у сравним с величиной зазора.

Естественно начать анализ с таких длин волн возмущений или вихрей, размер которых в направлении г также совпадает с величиной зазора.

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

2. Постановка задачи

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

Для координат и функций течения вводятся следующие обозначения: ¿^1/^1, гК\, гЩ,ш, иШ\, ■шШх, р0Ш2, р3р, НШ2, аШ\ соответственно для време-

ни, осевой, радиальной и угловой координат и для соответствующих компонентов вектора скорости, для статического давления, плотности, полной энтальпии, для коэффициента вязкости и для скорости звука. Индекс «1» соответствует параметрам внутреннего цилиндра (в случае если внутренний цилиндр покоится, используется скорость внешнего цилиндра). Индекс «в» соответствует значениям термодинамических функций газа в условиях отсутствия вращения цилиндров.

3. Система оценок

Оценим вначале перепад давления между стенками цилиндров, возникающий под действием центробежных сил

Ар — (1)

предполагая, что скорость в окружном направлении порядка единицы, й - ширина зазора между цилиндрами, а % = Я-1 - кривизна.

Индуцируемая этим перепадом скорость в направлении г (предполагая, что перепад давления в радиальном направлении такой же, как и в трансверсальном) оценивается следующим образом:

■ш — х1 /2. (2)

Эта оценка, строго говоря, имеет место, если силы вязкости несущественны. Мы рассматриваем именно режимы, для которых или вязкость несущественна, или граничные, для которых вязкость и инерционные силы сравнимы.

Поскольку мы ограничиваемся режимами, где от окружной координаты ничего не зависит и неустойчивость (в отличие от вихрей Тэйлора-Гертлера на пластине) развивается как нестационарный процесс, то уравнение неразрывности в невырожденном виде должно содержать производные как в радиальном направлении, так и в трансверсальном. Тогда при заданном масштабе по г, характеризующемся величиной Л:

■ш V х1/2Л3/2

— — - ил и V —---. (3)

ли л

Оценим теперь в уравнении радиального импульса влияние вязкости

' - ¿или 4 - ^^, где е = Ее-/2, (4)

а а2 X

то есть имеется такая иерархия. Если зазор и длина волны сравнимы, тогда, если они совпадают с этим вязким размером, получатся параболизованные уравнения Навье-Стокса, в которых вязкость работает в двух направлениях: в радиальном и трансверсальном. Если эти размеры больше вязкого размера, тогда у нас получаются уравнения Эйлера.

Далее, если длина волны больше зазора, тогда есть предельная длина волны, при которой в течении проявится влияние сил вязкости, но уже только в одном направлении -радиальном. Поэтому

др дри дру 1 + -к + £+""=0 (5)

§8? + + I2 + + 1 (»»< = Н{2 [2#■«(Iй - ¡т -1»)] +

от в. ^ в. вт г" Н^ 3 ^ в.. г" \ в. вт г ' J

+ з[ ¡Г М — у? + 2 ¡г — Г] + I? + ¡г]}, "От + + ~ёг~ + + Ру2 — = Не {[ Угг 1? + §1] +

+ 2 . г ¡8 I 2¡У _ ГП _ 2 г ¡8 _ I 2„] + 3§т + 2§т г]} ЗНе [ §т + 2и],

в™ + ¡(у.7) + ¡м + 21 = 1 ¡_ г^] + 1 ¡_ г а (—7 + ¡г)] + от 1 1 вт гг Не а. ^ а. J Не ат ат ^ г ат ' J

I [_ад | ]

+ Не [ г + !г ],

(6)

(7)

(8)

а ++ Рдг +1 =Не{¡¡г 1 ^ + ¡г+ ¡¡г] + + н I? + ¡Г ] +1 И—И + 2 —«] +1И [—И + 2 ¡Г —«] + [ I? + ¡Г ] + + §Г ¡? + )] +1 ¡¡Г ¡1 +2 ¡Г — Г )] +а?[ I? +¡1] +1 !г[—+2 ¡Г — +

+__1 ^(„УТч ..§Т + (. §Т)] +

+ РГ(7-1)М2 Не [§2 ч г Ут + вт )] +

+Не У. [^ ¡7 ] + Не ¡Г 7 + )] + Н1е 7 + §7)].

где Ые = , Рг - число Прандтля, М - число Маха.

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

Линии АС соответствуют вихри, на которые влияет вязкость, а характерный размер зазора много меньше, чем длина волны (рис. 1).

Итак, введем следующие представления координат и функций течения:

г = \х\, г = ёг1, £ = дГ1/211\, ад = -шо + ..., и = д1/2и1 + ..., V = в3/2\-1у1 + ..., р = ро + dpl + ..., р = ро + ... Подстановка разложений в исходную систему уравнений дает

^о + + ^ = 0,

т1 ог1 ог1

¡¡^ + ¡¡о?1 + ¡21 + й^™ = { 2 [2 у. П1П2^( I? — ¡Г )] +

+ 1 [П1^ ¡Г 1—Ё + 2 §г ]}, 1 ;

дрото + д(рои1то) + д(роУлто)

ди

дх1

дг1

д

-Пл Щ —

дг

дw дг

- П1П-л I

д ( ^ дг \дг

[Хт + т- + ]П2 + Ш - ро< = П1П2{[дм% + П2%] + + I дтг^-дх + 2 тИг]},

(12) (13)

дрЕ + дриН + дpvH

дг

дх

+

дг

рг(7-11)м 2 [П1 П2 ж (Р> ш) + П1П21 дг (Р жт)] +

+ П1П2 д [ПИ дт ] +П1П-1 д [Ы дЦ)]

-1 д

' дт ^

где Пл = -ЦЗ/2, П| =

(14)

Рис. 1. Диаграмма различных возможных режимов точения Куэтта Тойлора

П1

дс.ляст отношение диффузионных эффектов к инерционным эффектам. По существу, это

П2

рины зазора к длине волны пространственных возмущений (вихрей). Напомним, что во всех случаях мы имеем дело с безразмерными параметрами.

Общая система уравнений соответствует точке А, для которой оба введенных параметра подобия сохраняют конечное значение.

Для дальнейших) анализа отметим, что возможны различные предельные переходы. Если мы будем двигаться вдоль линии АВ, параметр П1 = 0, а параметр П2 сохраняет конечное значение. Система уравнений для этого предельного перехода имеет вид:

дро + дроУл + дро Ул = о

дЪ1 дх1 дг1 ' дрои1 1 дрои2 1 дрл 1 дрои1у1

+

дгл ' дгл + дгл + дгл дро-шо , д(рощ-шо) , д(роУл'Шо)

ди

+

дгл

+

дгл

=0,

0,

дроУл + дроЩУл + дроу"2

дП

дгл

дгл

тт2 дрл 2

П2 + д^ - РоЩ = 0,

дpоWо + д (ро Щ'Шо) + д(роУ^о) = о

дгл

дгл

дрЕ + дриН + друН о

дг

дг

дг

(15)

(16)

(17)

(18)

(19)

(20)

Соответственно, параметры, находящиеся вдоль линии АС, приводят к режиму течения, в котором вязкость существенна (в одном направлении):

дро + дррщ + дро vi dt\ dzi дг\

0,

dp0u\ dp0u2

dU

+

1 + dpi + dpouivi

dz\ dz\

дгл

д

дг1

JL ()

^ дг1 дг1

0,

др1 2 п

---poW2 = 0,

or1

дрЕ dpuH dpvH

1

dt

dz

dr Pr(7 - 1)M2

д

1 dr \dУ

^ д i du ^ ^ ^v \ dr \ dz dr J

+ П dr

i dw\ \dr J

4. Выводы

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

Литература

1. Kao К.Н., Chow C.Y. Linear stability of compressible Tavlor-Couette flow // Physics of Fluids A: Fluid Dynamics. - 1992. - V. 4, I. 5. - P. 984-996.

2. Боголепов В.В., Липатов И.И. Влияние сжимаемости на структуру вихрей Гертлера // Механика жидкости и газа. — 1997. — Т. 1. — С. 36-47.

3. Кочип Н.Е., Кибель И.А., Розе Н.В. Теоретическая гидромеханика. — М. : Физматлит, 1963.

4. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. — Т. 6. Гидромеханика. — М. : Наука, 1986.

Поступила в редакцию 17.07.2013.