Асимметричные модуляции ускорения свободного падения в задаче о свободной конвекции в замкнутой полости

Размещено на http://www.allbest.ru/

4

Пермский государственный университет

Асимметричные модуляции ускорения свободного падения в задаче о свободной конвекции в замкнутой полости

А.Б. Мелентьев, Е.Л. Тарунин

Аннотация

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

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

Annotation

Asymmetric modulation of free flow acceleration in free convection task in closed cavity

A. B. Melentyev, E. L. Tarunin Perm State University

Heat convection with the heating from below is considered with math modeling method. The task is heat convection in a square closed cavity. Modulation of free flow acceleration (symmetric and asymmetric) was maintained by moving the cavity in vertical direction. It was found that deviation from oscillations symmetry leads to decreasing of oscillation amplitude of convection characteristics.

Key words: convection; modulation; asymmetry.

1. Введение

Используемая модуляция была вызвана смещением полости в вертикальном направлении по закону

(1)

Вторая производная от смещения (1) создавала дополнительное модуляционное ускорение:

(2)

Параметрами модуляции являются: а - амплитуда, Т₀ - полный период модуляции и параметр асимметрии . Значение параметра асимметрии, равное единице, соответствует симметричной модуляции. Эффекты, связанные с влиянием параметра асимметрии, были исследованы в различных задачах [1]. В наших исследованиях фиксированными были параметр асимметрии , число Прандтля Pr = 1. Рассматривались случаи как без модуляции (амплитуда смещения и частота модуляций равны нулю), так и с модуляцией: амплитуда смещения а = 0.1, частота модуляции изменялась для поиска резонансных эффектов в интервале от 0 до значения частоты в 2 раза большее, чем частота собственных (без модуляции) затухающих колебаний. Рассматривалась как симметричная модуляция (), так и асимметричная ().

Для решения задачи использовались уравнения конвекции в приближении Буссинеска, решались уравнения двухполевым методом [2, 3]:

(3)

где - функция тока и вихрь скорости, есть сумма постоянной и модуляционной компонент числа Грасхофа. Модуляция ускорения свободного падения входит в модуляционную компоненту числа Грасхофа:

(4)

Уравнения (3) записаны в безразмерных переменных функции тока и вихря скорости. В качестве единиц обезразмеривания расстояния и времени были выбраны размер квадратной полости L и характерное вязкое время н/L². Граничные условия соответствовали твердым непроницаемым границам с заданной температурой .

Без модуляции такая задача решалась в работе [1].

2. Результаты численного решения без модуляций

В качестве начального состояния при получении результатов использовалось стационарное решение с возмущение вихря скорости в центре полости . Этому возмущению соответствовало максимальное значение функции тока ш_m0 ? A3.310-³. В качестве амплитуды возмущений обычно использовалось значение A = 100.

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

. (5)

Уравнение Пуассона для функции тока решалось на каждом шаге по времени методом последовательной верхней релаксации с оптимальным параметром релаксации [3] до выполнения неравенства

. (6)

Вихрь скорости на границах вычислялся по формуле Тома. Шаг квадратной пространственной сетки равнялся 1/40 (проверочные расчеты выполнялись с шагом 1/60). В процессе счета следили за тем, чтобы сеточное число Рейнольдса [3] не превосходило 2 (это условие необходимо для выполнения условия устойчивости схемы с центральными разностями в конвективных слагаемых).

Колебательный характер установления стационарного решения для двух значений числа Грасхофа из интервала 710³ < Gr < Gr*? (61±1)10³ приведен на рис. 1.

а б

Рис. 1 Зависимость максимума функции тока от времени при а) Gr = 10⁴, б) Gr = 310⁴

На рисунках показан способ вычисления периода "собственных" колебаний в процессе установления. Как видно, период этих колебаний уменьшается с ростом числа Грасхофа. Полученные значения частоты щ₀(Gr) представлены на рис. 2.

Рис. 2 Зависимость щ₀(Gr)

Методом наименьших квадратов найдена аппроксимационная формула в виде корневой зависимости (пунктирная линия на рис. 2), справедливой для 710³ < Gr < 5010³ с погрешностью не более 20%:

. (7)

В вычислительных экспериментах отслеживался процесс установления стационарного режима с погрешностью менее 1%:

,. (8)

где n - номер шага по времени, - максимумы по всей полости функции тока и вихря скорости, Nu - число Нуссельта (отношение интенсивности теплового потока за счёт конвекции к интенсивности теплового потока за счет теплопроводности через всю полость) [2].

Для стационарного решения подкорректирована корневая зависимость максимума функции тока (закон Ландау) от числа Грасхофа [4]:

, (9)

справедливая c погрешностью не более 4% при 610³ < Gr < 5810³. Заметим, что число Прандтля в расчетах было равным 1 и потому число Релея равно числу Грасхофа.

Зависимость максимума функции тока установившегося течения ш₀ от числа Грасхофа на интервале до Gr = 6010³ представлена на рис. 3 сплошной линией.

Рис. 3 Зависимость ш(Gr)

Штриховые линии определяют характеристики процесса установления. Напомним, что выход на стационарное решение происходил посредством затухающих колебаний. Верхняя линия соответствует первому всплеску максимума функции тока ш_m1 (рис. 1, б), а нижняя линия соответствует первому минимуму после первого всплеска ш_m2. Как видно, с ростом числа Грасхофа увеличивается размах колебаний установления (ш_m1 _ ш_m2).

Рис. 4 Зависимость ш_m1 _ ш_m2 от Gr

Для этой зависимости справедлива на интервале 2010³ < Gr < 6010³ формула с погрешностью не более 2%:

. (10)

3. Результаты численного решения с учётом модуляций

Перейдем к обсуждению результатов расчета при модуляции согласно формуле (4).

Частота модуляций щ соответствует полному периоду модуляции T₀ из формул (1) и (2):

. (11)

Как при симметричных, так и при асимметричных модуляциях реализуется установившийся режим колебаний. Типичный пример установившегося режима колебаний при значении числа Грасхофа Gr = 310⁴ и частоты модуляций щ = 51.4 (соответствующей частоте собственных затухающих колебаний щ₀) представлен на рис. 5.

а б

Рис. 5 Установившиеся колебания при а) симметричной и б) асимметричной модуляции, Gr = 310⁴, щ = 51.4

В вычислительных экспериментах отслеживался процесс установления модуляционных колебаний с погрешностью менее 1%:

, (12)

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

На рис. 6 и 7 представлены изолинии функции тока и температуры в момент максимума и минимума значений функции тока при симметричных модуляционных колебаниях, изображенных на рис. 5,а.

Рис. 6 Картина течения при максимальном значении функции тока колебательного режима при симметричной модуляции, Gr = 310⁴

Рис. 7 Картина течения при минимальном значении функции тока колебательного режима при симметричной модуляции, Gr = 310⁴

На рис. 8 и 9 представлены зависимости амплитуды функции тока, как разности максимального и минимального значений, от частоты модуляции щ при Gr = 10⁴ и Gr = 310⁴.

На рис. 8 и 9 щ₀ - это частота собственных колебаний без модуляции, полученная на рис. 2.

Рис. 8 Амплитудная кривая для Gr = 10⁴

Рис. 9 Амплитудная кривая для Gr = 310⁴

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

Представленные выше результаты соответствовали интервалу чисел Грасхофа от Gr₁ = (5.5±0.5)10³ до Gr*? (61±1)10³. При значениях числа Грасхофа Gr > Gr*? (61±1)10³ процесс установления собственных затухающих колебаний без модуляций приводил, как и в [3, 4], к регулярным колебаниям. Исследование этой области значений числа Грасхофа планируется в дальнейшем.

Выводы

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

2. Показано, что при модуляции в области рассмотренных чисел Грасхофа реализуется установившийся колебательный режим с частотой модуляции как для симметричной, так и для асимметричной модуляции.

3. Найдена зависимость амплитуды колебаний, вызванных модуляцией для симметричной и асимметричной модуляции. Показано наличие резонансного эффекта модуляции в случае симметричной модуляции.

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

Тарунин Е.Л. Обзор особенностей асимметричных колебаний // Проблемы механики и управления. Нелинейные динамические системы / Перм. ун-т. Пермь, 2005. №37. С.169-187.

Гершуни Г.З., Жуховицкий Е.М. Конвективная устойчивость несжимаемой жидкости. М.: Наука, 1972. С.242-255.

Тарунин Е.Л. Вычислительный эксперимент в задаче свободной конвекции: учеб. пособие / Иркутск. ун-т. Иркутск, 1990. 228 с.

Гершуни Г.З., Жуховицкий Е.М., Тарунин Е.Л. Численное исследование конвекции жидкости, подогреваемой снизу // Известия АН СССР. Механика жидкости и газа. 1966, №6.

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