ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА 2011 Математика. Механика. Информатика
Вып. 1(5)
УДК 532.5
Асимметричные модуляции ускорения свободного падения в задаче о свободной конвекции в замкнутой полости
А. Б. Мелентьев, Е. Л. Тарунин
Пермский государственный университет, Россия, 614990, Пермь, ул. Букирева, 15 a.b.melentyev@mail.ru; 8 902 79-97-681
Методами математического моделирования исследованы эффекты, возникающие при модуляции ускорения свободного падения в случае тепловой конвекция жидкости, подогреваемой снизу. Решалась задача: конвекция в плоской полости квадратного сечения. Модуляция силы тяжести (симметричная и асимметричная) была образована перемещением полости в вертикальном направлении. Выяснено, что отклонение от симметрии колебаний приводит к уменьшению амплитуды колебаний характеристик конвективного течения.
Ключевые слова: конвекции; модуляции; асимметрия.
Введение
Используемая модуляция была вызвана смещением полости в вертикальном направлении по закону
УО ) =
\а С08(ю/), 0 < t < 7^
[-а соб(ю2 (, -Т1)), Т1 < t < Т0. ю1 = п / 7|,
®2 = П/(70 - 71).
(1)
Вторая производная от смещения (1) создавала дополнительное модуляционное ускорение:
1-аю,2 cos а, 0 < t < 71, ...
т==1 У, Т, т '<,<т (2)
[аю2 cosю2(t-71), 71 <, <70.
Параметрами модуляции являются: а -амплитуда, Т0 - полный период модуляции и параметр асимметрии % = а1/ а2 = (Т0 - Т1) / Т1 . Значение параметра асимметрии, равное еди-
нице, соответствует симметричной модуляции. Эффекты, связанные с влиянием параметра асимметрии, были исследованы в различных задачах [1]. В наших исследованиях фиксированными были параметр асимметрии % = 2, число Прандтля Рг = 1. Рассматривались случаи как без модуляции (амплитуда смещения и частота модуляций равны нулю), так и с модуляцией: амплитуда смещения а = 0.1, частота модуляции изменялась для поиска резонансных эффектов в интервале от
0 до значения частоты в 2 раза большее, чем частота собственных (без модуляции) затухающих колебаний. Рассматривалась как симметричная модуляция (% = 1), так и асимметричная (% = 2).
Для решения задачи использовались уравнения конвекции в приближении Бусси-неска, решались уравнения двухполевым методом [2, 3]:
© А. Б. Мелентьев, Е. Л. Тарунин, 2011
dm .dw dm dw dm ^ . ,dT
— + (———----------——) = Am + Gr (t)—,
dt dy dx dx dy dx
dT dw dT
— + (—— dt dy dx
Aw+m = o,
dT) = 1 AT,
dx dy Pr
(3)
где /, а - функция тока и вихрь скорости,
Gr(t) = Gr0 + Grmod(t) есть сумма постоянной
и модуляционной компонент числа Грасхофа. Модуляция ускорения свободного падения входит в модуляционную компоненту числа Грасхофа:
g (t) =
-a(—) cos — ,0 < t < Tj,
T
a(
T -T *0 Jj
T
)2 cos
x(t - Tj)
To - Tj
(4)
Tj < t < T0,
Уравнения (3) записаны в безразмерных переменных функции тока и вихря скорости. В качестве единиц обезразмеривания расстояния и времени были выбраны размер квадратной полости L и характерное вязкое время v/L . Граничные условия соответствовали твердым непроницаемым границам с заданной температурой Т |г = 1 - у .
Без модуляции такая задача решалась в работе [1].
2. Результаты численного решения без модуляций
В качестве начального состояния при получении результатов использовалось стационарное решение с возмущение вихря скорости в центре полости р0 = Ах(1 - х)у (1 - у) . Этому возмущению соответствовало максимальное значение функции тока щт0 ~ А-3.3-10-3. В качестве амплитуды возмущений обычно использовалось значение А = 100.
Решение системы выполнялось по явной двухслойной схеме с аппроксимацией первых производных конвективных слагаемых центральными разностями. Шаг по времени вычислялся по формуле с запасом, обеспечивающим выполнение условия устойчивости:
At =-
h2
12 + h wm
(5)
Уравнение Пуассона для функции тока решалось на каждом шаге по времени методом последовательной верхней релаксации с оптимальным параметром релаксации [3] до выполнения неравенства
maxl W -V"j 1
ij ______________
< o.oj.
(6)
l w l
I r max I
Вихрь скорости на границах вычислялся по формуле Тома. Шаг квадратной пространственной сетки равнялся 1/40 (проверочные расчеты выполнялись с шагом 1/60). В процессе счета следили за тем, чтобы сеточное число Рейнольдса [3] не превосходило 2 (это условие необходимо для выполнения условия устойчивости схемы с центральными разностями в конвективных слагаемых).
Колебательный характер установления стационарного решения для двух значений числа Грасхофа из интервала
7-103 < Gr < Gr*~ (6 j±l)-103 приведен на рис. 1.
¥
-
Ж-4-
:чг 1
.¡_ _ Т :
0 0.2 0.4
12
10
8
6
4
2
О
Vml J
Wo j
kJ
: i - T 4 -i ...24 i
•iTTTTVVl
0.1 0.2
б
t
0.3
Рис. 1. Зависимость максимума функции тока от времени при а) Gr = 104, б) Gr = 3-104
На рисунках показан способ вычисления периода "собственных" колебаний в процессе установления. Как видно, период этих колебаний уменьшается с ростом числа Грас-хофа. Полученные значения частоты ю0^г) представлены на рис. 2.
О 20 000 40 000 60 000
Рис. 2. Зависимость m0(Gr)
а
Методом наименьших квадратов найдена аппроксимационная формула в виде корневой зависимости (пунктирная линия на рис. 2), справедливой для 7-103 < Gr < 50-103 с погрешностью не более 20%:
со0(Ог) * 0.632>/От -5621. (7)
В вычислительных экспериментах отслеживался процесс установления стационарного режима с погрешностью менее 1%:
кт+х-к„і+1^"+1 -ки"і+іс-^і ;001. (8)
і ах І ІКи"І К^і ■ ,
где П - номер шага по времени, Ктах, ®тах -максимумы по всей полости функции тока и вихря скорости, Ыи - число Нуссельта (отношение интенсивности теплового потока за счёт конвекции к интенсивности теплового потока за счет теплопроводности через всю полость) [2].
Для стационарного решения подкорректирована корневая зависимость максимума функции тока (закон Ландау) от числа Грас-хофа [4]:
К (От) * 0,063>/От - 5011, (9)
справедливая с погрешностью не более 4% при 6-103 < От < 58-103. Заметим, что число Прандтля в расчетах было равным 1 и потому число Релея равно числу Грасхофа.
Зависимость максимума функции тока установившегося течения ^0 от числа Грасхофа на интервале до От = 60-103 представлена на рис. 3 сплошной линией.
ф .Л г"'
Г ІІ п
> Ж Р и
і * . . • ^ 1» * * * » 1
ф п 12
о 9
О 20 000 40 000 60 000
Рис. 3. Зависимость щ(Ог)
Штриховые линии определяют характеристики процесса установления. Напомним, что выход на стационарное решение происходил посредством затухающих колебаний. Верхняя линия соответствует первому всплеску максимума функции тока (рис. 1, б), а нижняя линия соответствует первому минимуму после первого всплеска щт2. Как видно, с ростом числа Грасхофа увеличивается размах колебаний установления (щт1 - щт2).
ф у рг
А у /
у Г
Ґ
у Г
У у 6г
О 20000 40000 60 000
Рис. 4. Зависимость ут1 - ут2 от Gr
Для этой зависимости справедлива на интервале 20-103 < Gr < 60-103 формула с погрешностью не более 2%:
(/т1 - /т2) * 0 003 15 ■ ^ - 2 96 . (10)
3. Результаты численного решения с учётом модуляций
Перейдем к обсуждению результатов расчета при модуляции согласно формуле (4).
Частота модуляций ш соответствует полному периоду модуляции Т0 из формул (1) и (2): а = 2л/Т0. (11)
Как при симметричных, так и при асимметричных модуляциях реализуется установившийся режим колебаний. Типичный пример установившегося режима колебаний при значении числа Грасхофа Gr = 3 ■ 104 и частоты модуляций ш = 51.4 (соответствующей частоте собственных затухающих колебаний со о) представлен на рис. 5.
Рис. 5. Установившиеся колебания при а) симметричной и б) асимметричной модуляции, Gr = 3-104, т = 51.4
В вычислительных экспериментах отслеживался процесс установления модуляционных колебаний с погрешностью менее 1%:
1^
( Я+1) тах
г
(п)
| | Ыи - +
(п + 1) тах
-Ыи(п) |
тах I
0.01 (12)
1^1 I Ыи^х |
где п - номер колебательного периода (один период - от минимума до следующего минимума), ^тах, Ыитах - максимум функции тока
по всей полости и числа Нуссельта через полость за один период колебания.
На рис. 6 и 7 представлены изолинии функции тока и температуры в момент максимума и минимума значений функции тока при симметричных модуляционных колебаниях, изображенных на рис. 5,а.
Рис. 6. Картина течения при максимальном значении функции тока колебательного режима при симметричной модуляции, Gr = 3 -104
Рис. 7. Картина течения при минимальном значении функции тока колебательного режима при симметричной модуляции, Gr = 3 -104
На рис. 8 и 9 представлены зависимости амплитуды функции тока, как разности максимального и минимального значений, от частоты модуляции т при Gr = 10 Gr = 3-104
и
На рис. 8 и 9 ю0 - это частота собственных колебаний без модуляции, полученная на рис. 2.
■шп
- сим.модул.
- асим.модул.
Рис. 8. Амплитудная кривая для Gr = 10
1,5
0,5
. фтах 1 \
1 : ч
- 1 \ 1 \
- 1 1 1 1
-шГ1
40
— сим.модул.
80
120 асим.модул.
4
Рис. 9. Амплитудная кривая для Gr = 3-10
Как видно из рис. 8 и 9, симметричные модуляции с частотой, равной частоте собственных колебаний ш0, приводят к резонансу (прерывистая линия), однако при асимметричных модуляциях такого эффекта резонанса не наблюдается (сплошная линия), что аналогично результатам, полученным в механической задаче [1].
Представленные выше результаты соответствовали интервалу чисел Грасхофа от Gr1 = (5.5±0.5)-103 до Gr*~ (61±1)-103. При значениях числа Грасхофа Gr > Gr*« (61±1)-103
процесс установления собственных затухающих колебаний без модуляций приводил, как и в [3, 4], к регулярным колебаниям. Исследование этой области значений числа Грасхофа планируется в дальнейшем.
Выводы
1. Найдены зависимости характеристик колебательного режима установления стационарного решения (частоты и амплитуд) в области чисел Грасхофа, где существует
стационарное решение в отсутствии модуляции.
2. Показано, что при модуляции в области рассмотренных чисел Грасхофа реализуется установившийся колебательный режим с частотой модуляции как для симметричной, так и для асимметричной модуляции.
3. Найдена зависимость амплитуды колебаний, вызванных модуляцией для симметричной и асимметричной модуляции. Показано наличие резонансного эффекта модуляции в случае симметричной модуляции.
Список литературы
1. Тарунин Е.Л. Обзор особенностей асимметричных колебаний // Проблемы механики и управления. Нелинейные динамические системы / Перм. ун-т. Пермь, 2005. №37. С. 169—187.
2. Гершуни Г.З., Жуховицкий Е.М. Конвективная устойчивость несжимаемой жидкости. М.: Наука, 1972. С.242-255.
3. Тарунин Е.Л. Вычислительный эксперимент в задаче свободной конвекции: учеб. пособие / Иркутск. ун-т. Иркутск, 1990. 228 с.
4. Гершуни Г.З., Жуховицкий Е.М., Тарунин Е.Л. Численное исследование конвекции жидкости, подогреваемой снизу // Известия АН СССР. Механика жидкости и газа. 1966, №6.
Asymmetric modulation of free flow acceleration in free convection task in closed cavity
A. B. Melentyev, E. L. Tarunin
Perm State University, Russia, 614990, Perm, Bukireva st., 15 a.b.melentyev@mail.ru; +7 902 79-97-681
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.