Линейная устойчивость адвективного течения

Размещено на Allbest.ru

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

линейная устойчивость

Ключевые слова: адвективные течения; устойчивость; метод дифференциальной прогонки.

Введение

Одной из важнейших задач гидродинамики является исследование устойчивости конвективных и адвективных течений жидкости, которые могут быть описаны аналитически [1, 2]. Широкое распространение в задачах устойчивости конвективных и адвективных течений получили метод Галёркина, метод пошагового интегрирования с ортогонализацией, метод дифференциальной прогонки и метод сеток. В первых трех случаях задача устойчивости сводится к многократному решению краевой задачи для обыкновенных дифференциальных уравнений. В последнем случае решается краевая задача для одномерной системы уравнений в частных производных.

Адвективное течение, возникающее в плоском горизонтальном слое жидкости при отсутствии вращения для случая, когда температура на обеих границах линейно изменяется с продольной координатой, было впервые описано Р. В. Бирихом аналитически [3]. Особенностью такого течения является отсутствие вертикальной компоненты скорости, вектор скорости в потоке здесь ориентирован перпендикулярно силе плавучести, которая является основной причиной движения. Устойчивость таких течений, возникающих в горизонтальном слое жидкости при отсутствии вращения, исследовалась в основном методом дифференциальной прогонки и изучена достаточно подробно в работах В. М. Мыз-никова совместно с другими авторами [1].

Адвективное течение, возникающее во вращающемся плоском горизонтальном слое жидкости, для случая, когда температура на границах линейно изменяется с продольной координатой, было впервые выведено аналитически С. Н. Аристовым [4]. Как и в случае без вращения, у адвективного течения отсутствует вертикальная компонента скорости, вектор скорости в потоке ориентирован перпендикулярно силе плавучести, однако имеются обе горизонтальные компоненты вектора скорости. Устойчивость адвективного течения, возникающего во вращающемся горизонтальном слое жидкости, исследовалась К. Г. Шварцем с использованием метода сеток в работах [5, 6].

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

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

Рассмотрим плоский горизонтальный слой несжимаемой жидкости, ограниченный твёрдыми плоскостями и вращающийся с постоянной угловой скоростью , где - орт-вектор вертикальной оси. Направление оси вращения совпадает с вертикальной осью координат. Изучение адвективных течений будем производить на основе уравнений конвекции в приближении Буссинеска во вращающейся системе отсчета в декартовой системе координат [1]. Выбрав в качестве единиц измерения длины, времени, скорости, температуры и давления соответственно , , , , (здесь - коэффициент вязкости, - ускорение силы тяжести, - коэффициент теплового расширения, - постоянный горизонтальный температурный градиент на границе слоя, - средняя плотность), мы получим исходные уравнения в безразмерном виде [5]:

(1)

где - вектор скорости, - температура, - давление, - число Грасгофа, - число Тейлора, - число Прандтля. Здесь - коэффициент температуропроводности. На обеих плоскостях задана температура, линейно меняющаяся с продольной координатой , условия прилипания и замкнутости потока

(2)

Профили скорости и температуры возникающего адвективного течения описываются в виде точного решения задачи (1) - (2) [5]

(3)

где , и - компоненты вектора скорости, - действительная часть комплексного значения, - мнимая часть комплексного значения.

Для исследования устойчивости стационарного плоскопараллельного адвективного течения применим метод малых возмущений [1]. Рассмотрим возмущенное течение , , , где , , - малые нестационарные возмущения. Пренебрегаем квадратичными по и слагаемыми и получаем

(4)

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

(5)

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

Учитывая дивергентность возмущений скорости, введем функцию тока возмущений : и вихря возмущения скорости . В качестве малых возмущений рассматриваем нормальные возмущения вида

(6)

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

Подставляя (6) в (4) и избавляясь от давления, получаем систему обыкновенных дифференциальных уравнений вида

(8)

с граничными условиями (5) в виде

(9)

Здесь , , , , , , , .

Краевая задача (8) - (9) решалась методом дифференциальной прогонки с помощью пакета прикладных программ «Гидродинамическая устойчивость» [8].

2. Анализ результатов

Проводились тестовые расчёты для случая отсутствия вращения, которые полностью совпали с результатами работы Мызникова [1, 9].

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

Для значений числа Тейлора, близких к нулю, минимальное критическое значение числа Грасгофа () достигается при . На рис. 1 представлена нейтральная кривая при , .

Рис. 1. Нейтральная кривая при очень слабом вращении (Ta=10-12)

При слабом вращении основные расчеты проводились при . Результаты показали, что минимальное критическое значение числа Грасгофа достигается при . Критические значения числа Грасгофа, вычисленные с различной точностью интегрирования, представлены в табл. 1.

Таблица 1. для различной точности

интегрирования при слабом вращении

Оптимальная точность интегрирования порядка , так как разница с результатами более высокой точности невелика, а экономия времени, по сравнению с решением, полученным при большей точности, очень существенна. Разница между результатами работы [5] составляет менее 0, 03%. Время работы программы при вычислении критического числа Грасгофа с различной точностью интегрирования при представлено в табл. 2.

Таблица 2. Время работы программы при различной точности интегрирования

С ростом числа Тейлора критическое число Грасгофа растет, вращение стабилизирует адвективные течения.

При умеренном вращении устойчивость течения снижается. Основные расчеты проводились при . Критическое значение числа Грасгофа достигается при . Критические значения числа Грасгофа, вычисленные с различной точностью интегрирования, представлены в табл. 3.

Таблица 3. для различной точности интегрирования при умеренном вращении

Разница между результатами работы [5] составляет менее 0, 006%. В этом случае нейтральная кривая имеет специфическую форму «носика» (см. рис. 2).

При быстром вращении основные расчёты проводились при . Критическое значение числа Грасгофа () достигается при . Разница между расчетами, полученными методом дифференциальной прогонки и методом сеток, также составляла менее 0, 006%. С ростом числа Тейлора вращение стабилизирует адвективное течение.

Рис. 2. Нейтральная кривая при умеренном вращении (Ta=500)

Заключение

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

Список литературы

Гершуни Г. З., Жуховицкий Е. М., Непомнящий А. А. Устойчивость конвективных течений. М. : Наука, 1989. 320 с.

Тарунин Е. Л., Шварц К. Г. Исследование линейной устойчивости адвективного течения методом сеток // Вычислительные технологии. 2001. Т. 6, №6. С. 108-117.

Размещено на Allbest.ru