CC BY
39
УДК 537.84
А. В. Проскурин, А. М. Сагалаков
Устойчивость течения Озеена возле кругового цилиндра* A. V. Proskurin, A.M. Sagalakov A Stability of Oseen Flow by Circular Cylinder
Исследуется устойчивость течения Озеена возле кругового цилиндра. Работа носит методический характер. Изучены возможности метода на основе функций Рвачева. Построена структура, удовлетворяющая граничным условиям и приближенно представляющая решение задачи. Неизвестные коэффициенты находились методом коллокаций. Течение устойчиво в области применимости формулы Озеена. Метод прост и экономичен, имеет высокий порядок точности. Ключевые слова: гидродинамическая устойчивость, течение Озеена, уравнения с частными производными, спектральные методы.
A stability of Oseen flow by circular cylinder is investigated. The work has methodical character. Rvachev function is applied for illustrating possibilities of the method. The authors build the structure that satisfies the boundary conditions and represents approximately solution. Unknown factors were found by a collocation method. The flow is stable in the field of applicability of the Oseen formula. The method is simple and exponentially converging.
Key words: hydrodynamic stability, Oseen flow, partial differential equations, spectral methods.
Исследование устойчивости и бифуркаций течений представляет собой сложную и интересную проблему. Классическая постановка задачи устойчивости течений вязкой жидкости приводит к задаче Орра-Зоммерфельда [1, 2] - краевой задаче для системы обыкновенных дифференциальных уравнений. Это возможно благодаря наличию однородных направлений, относительно которых возмущение можно представить в виде элементарных волновых решений с определенным набором волновых чисел, например, для течения в плоском канале, течения в круглой трубе, течений между коаксиальными цилиндрами. С другой стороны, есть много течений, в которых картина развития возмущений неоднородна вдоль или поперек потока: течение в прямоугольной трубе, обтекание крыла, обтекание цилиндра. В таких случаях задача устойчивости сводится к краевой задаче для уравнений с частными производными (см. обзор: [3]).
Направим ось X вдоль направления движения жидкости, а ось у — перпендикулярно, обе перпендикулярно оси цилиндра (см. рис. 1). Уравнение Навье-Стокса для функции тока имеет вид
д к дш дАш дш дАш 1 -
— Аш------—+ ^—*- =—А>, (1)
дї дх ду ду дх Re
Ш
где у - функция тока; Яе=--- - число Рейнольдса;
V
и - скорость набегающего потока; й - радиус цилиндра; V - коэффициент вязкости.
В соответствии с [3], решение (1) ищем в виде
W = Wo = w(x, y)ect, (2)
f(r) - решение Озеена [4],
где у0 = y
л/Х^7
f(r)= ^-М-гіпг +M2r +M 3r 1, M2 = 1 + iMi, M3 = -iMi, Mi = 2ln7 406 , У)єС - малое возмущение, ^(x,y)
He
- амплитуда, С = Х + 11 - определяет декремент затухания возмущения. Подставляя (2) в (1), получим
С (у» + Ууу ) = ^ (^ + 2Уху + Ууууу )-
-2У-(1 - У2 + ¥уух),
где производные обозначены нижними индексами.
Рис. 1. Конфигурация потока
Работа выполнена при финансовой поддержке федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» (контракт 14.740.11.0355).
i15
физика
Зададим условия на границе расчетной области Q (см. рис. 1). На поверхности цилиндра возмущение скорости обращается в нуль. На входе в Q возмущение также задано равным нулю. На выходе из расчетной области Q, следуя [5], можно задать нулевой градиент скорости вдоль оси x . В этом случае возмущение может покидать расчетную область с ненулевой амплитудой. В работах [3, 6], наоборот, использовались нулевые граничные условия. При малых числах Рейнольдса возмущения локализованы возле поверхности цилиндра, так что примем
¥ = ¥х =Уу = 0. (4)
Решение (3), удовлетворяющее граничным условиям, будем искать в виде
V = °W£ а„Т, ( x ) Т, ( У ), (5»
I, j
где T - многочлены Чебышева первого рода; а, -
неизвестные коэффициенты; cr(Q) - уравнение границы. Уравнение границы составим при помощи R -конъюнкции [7]
^(Q) = ^ + °у + <г , (6)
где <Jsq, <Jcyl - уравнения прямоугольника и круга соответственно.
Уравнение круга имеет вид acyl = 1 - x2 - у2, уравнение для прямоугольника получим из уравнений полос <Jx = Ly- у2 и ау = L - x2 при помощи выражения (6). Аналитические преобразования производились в системе Maxima, выражения для cr(Q) и особенно его производные громоздки, и мы их приводить не будем.
Определив множество точек коллокации Гаусса-Лобатто x. = Lcos(7li/n),i = 0,... ,n , у. = Ly cosfnj / n), j=0,...,k, получим задачу на собственные значения
Av = CBv, (7)
где v = {a00, a01,..., an(k-1), ank } .
Первоначально задача на собственные значения (7) решалась с использованием QZ алгоритма из библиотеки LAPACK, которая ранее использовалась нами в расчетах устойчивости плоского течения Пуазейля в магнитном поле (см.: [8, 9]) вплоть до чисел Рейнольдса порядка 107. Однако вычислительные затраты этого метода в данном случае оказались очень велики, поэтому для расчетов были использованы методы Крылова [10, 11], реализованные библиотеками SLEPc и ARPACK. Численные эксперименты показали эффективность обратных итераций с нулевым сдвигом.
Как известно [4], формулы Озеена справедливы при числах Рейнольдса меньше единицы и только вблизи цилиндра. Невозможно построить решение,
удовлетворяющее граничным условиям на бесконечности. Поэтому примем следующие значения параметров: Ly = 2, Ь = 3 . Основываясь на результатах, приведенных в [12], можно предположить, что собственные функции при малых числах Рейнольдса ограничиваются данной областью. Количество точек коллокации по оси х составляло 80, по оси у — 70.
На рисунке 2 представлена зависимость коэффициента затухания X от числа Рейнольса. При увеличении числа Рейнольсда X увеличивается и при Re>1 стремится к нулю. Значения X меняются при изменении области расчетов, но сама зависимость Х^е) остается монотонно возрастающей, а течение устойчивым.
Рис. 2. Зависимости X (Re)
•3-2-10 1 2 3
х
Рис. З. Действительная часть собственной функции с m = 0 при Re = 0.1
х
Рис. 4. Действительная часть собственной функции с m = 2 при Re = 0.1
На рисунках 3-4 представлена действительная часть собственных функций с т = 0 и т = 2. Так как течение Озеена симметрично относительно осей X и у , то и данные собственные функции попарно симметричны относительно оси цилиндра с собственными функциями т = 1 и т = 3 соответственно. На рисунке 3 линии тока имеют вид одиночного вихря за цилиндром, а на рисунке 4 — двух вихрей с противоположным вращением.
Из полученных результатов можно сделать вывод об устойчивости течения возле кругового цилиндра в области применимости формулы Озеена. Вид собственных функций соответствует структуре переходных режимов течения, которые хорошо изучены независимыми методами, численными и экспериментальными. Предлагаемый метод намного проще аналогов [3, 5, 6], более экономичен, имеет высокий порядок точности.
Библиографический список
1. Henningson D. S., Schmid P. J. Stability and transition in shear flows. — New York, 2001.
2. Гольдштик М. А., Штерн В. Н. Гидродинамическая устойчивость и турбулентность. — Новосибирск, 1977.
3. Theofilis V. Advances in global linear instability analysis of non-parallel and three-dimensional flows // Progress in Aerospace Sciences. — 2003. — No 39.
4. Бэтчелор Д. Введение в динамику жидкости. — М., 1973.
5. Barkley D. Confined three-dimensional stability analysis of the cylinder wake. 2004 [Электронный ресурс]. — URL: http://arXiv.org/abs/physics/0405153v2
6. Gonzalez L., Theofilis V, Sherwin S. J. High-order methods for the numerical solution of the BiGlobal linear stability eigenvalue problem in complex geometries // International Journal for Numerical Methods in Fluids. — 2011. — Vol. 65. Issue 8.
7. Кравченко В. Ф., Рвачев В. Л. Алгебра логики, атомарные функции и вейвлеты в физических приложениях. — М., 2006.
8. Проскурин А. В., Сагалаков А. М. Новая ветвь неустойчивости магнитогидродинамического течения Пу-азейля в продольном магнитном поле // Письма в журнал технической физики. — 2008. — Т. 34, № 5.
9. Проскурин А. В., Сагалаков А. М. Устойчивость течения Пуазейля при наличии продольного магнитного поля // Журнал технической физики. — 2012. — Т. 82, вып. 5.
10. Saad J. Numerical methods for large eigenvalue problems. — Manchester, 1992.
11. Krylov-Schur methods in SLEPC. SLEPc technical report STR-7 / V. Hernandez, J. E. Roman, A. Tomas, V. Vidal. 2007 [Электронный ресурс]. — URL: www.grycap.upv.es/slepc/documentation / reports / str7. pdf.
12.Rogers S., Kwak D. An upwind-differencing scheme for theincompressible Navier-Stokes equations // NASA Technical Memorandum 101051. — 1988.
