Моделирование конвективных течений в CFD-пакетах

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

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

Институт механики сплошных сред УрО РАН

Моделирование конвективных течений в CFD-пакетах

А.В. Евграфова

Д.Ф. Ибраев

А.Н. Сухановский

Аннотация

Проведено численное исследование формирования конвективных течений в цилиндрическом слое жидкости с локализованным нагревом в центральной части. Для численной реализации задачи использовался конечно-элементный программный комплекс ANSYS. Расчеты проводились двумя различными методами, реализованными в CFD-пакетах CFX и FLUENT. Получены интегральные характеристики, характеризующие временную эволюцию течений. Проведено сравнение решений, полученных различными решателями, с экспериментально полученными данными.

Ключевые слова: конвекция; ANSYS CFX; ANSYS FLUENT.

The results of numerical studies of the formation of convective flows in a cylindrical fluid layer with localized heating in the central part are demonstrated. For the numerical solution of the problem we used finite element software suite ANSYS. Calculations were carried out in two different software packages: CFX and FLUENT. Integral characteristics of convection flow are presented. A comparison of solutions obtained by different solvers, with the experimental data is shown.

Key words: convection; ANSYS CFX; ANSYS FLUENT.

Содержание

Введение

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

2. Метод решения

3. Результаты

4. Сравнение результатов

5. Сравнение численных и экспериментальных данных

Заключение

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

Введение

2. Метод решения

Для численной реализации задачи в описанной постановке использовался конечно-элементный программный комплекс для моделирования течений жидкостей и газов ANSYS CFD. Технология ANSYS CFD открывает доступ к программным продуктам ANSYS FLUENT и ANSYS CFX, с помощью которых мы проведем расчеты.

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

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

3. Результаты

Были проведены расчеты для модели цилиндрического слоя жидкости с локализованным нагревом в центре. Радиус модели составляет 150 мм, толщина слоя - 30 мм, радиус нагреваемой области - 50 мм.

Рис. 3. Схематичное изображение конвективного движения

В случае подогрева жидкости в центре модели устанавливается режим течения, схематически изображенный на рис. 3. Локальный нагрев в центральной части дна создает вертикальный и горизонтальный градиенты температуры. Горизонтальный градиент температуры приводит к образованию обратной меридиональной ячейки [2]. Течение в нижней части направлено к центру, а над центральной частью формируется интенсивное подъемное течение. В верхней части слоя радиальное движение направлено к периферии. Таким образом, происходит формирование адвективного течения.

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

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

Рис. 5. Образование конвективных струй в области нагрева

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

4. Сравнение результатов

Были проведены расчеты исходной задачи в различных программных пакетах CFX и FLUENT. В качестве рабочей жидкости использовались вода (динамическая вязкость8.9*10?4 Па·с) и трансформаторное масло (6.7*10-3 Па·с). Физическое время расчетов - 150 с, этого было достаточно для выхода на периодический режим (имеется в виду формирование конвективных струй). Результаты приводятся в размерном виде, так как выбор безразмерных управляющих параметров для такой системы является нетривиальной задачей и это будет сделано в дальнейших работах.

Рис. 6. График зависимости средней кинетической энергии от времени при разной мощности нагрева для воды и масла (CFX)

На рис. 6 представлен график зависимости средней кинетической энергии от времени при разной мощности нагрева, полученного в CFX. Значение средней кинетической энергии в слое жидкости определяется следующим соотношением:

.

График показывает, что увеличение мощности нагрева приводит к заметно более интенсивному конвективному течению жидкости. Сравнивая графики средней кинетической энергии для воды и масла (рис. 6), можно наблюдать, что течение воды более интенсивно по сравнению с течением масла, это обусловлено более высокими значениями потокового числа Грассгофа [4]. Так же наблюдая за поведением изменения кинетической энергии масла, можно увидеть, что оно происходит с некоторой периодичностью.

Изучение режимов конвективного течения масла и воды во FLUENT дало похожие результаты.

Рис. 7. График зависимости средней кинетической энергии от времени при разной мощности нагрева для воды (FLUENT)

Рис. 8. График зависимости средней кинетической энергии от времени при разной мощности нагрева для масла (FLUENT)

Рис. 7 и 8 иллюстрируют изменение кинетической энергии со временем при разной мощности нагрева для масла и воды. График для масла показывает, что периодическое колебание кинетической энергии происходит на сравнительно высоких частотах - признак влияния мелкомасштабных структур, частота возникновения которых высока.

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

5. Сравнение численных и экспериментальных данных

конвективный течение жидкость нагрев

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

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

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

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

Рис. 9. Изменение частоты генерации мелкомасштабных структур н в точке от мощности нагрева Р

Из рис. 9 видно, что графики изменения частоты генерации мелкомасштабных структур, полученных в CFX и FLUENT, почти совпадают. А экспериментальные данные хорошо согласуются с численными при малых мощностях нагрева и существенно различаются при высоких мощностях. Это объясняется тем, что с ростом нагрева реальное течение, исследуемое в эксперименте, теряет осесимметричность, а в двумерных расчетах мы навязываем условие осесимметричности.

Заключение

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

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

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

Авторы благодарны кафедрам МСС и ВТ и теоретической физики Пермского национального исследовательского государ-ственного университета за предоставленную возможность проведения расчетов в пакете Ansys CFX и FLUENT.

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

1. Сухановский А.Н. Формирование дифференциального вращения в цилиндрическом слое жидкости // Вычислительная механика сплошных сред. 2010. Т.3, №2. С.103-115.

2. Batalov V., Sukhanovsky A. and Frick P. Laboratory study of differential rotation in a convective rotating layer // J. Geophys. Astrophys. Fluid Dynam. 2010. Vol.104, №4.P. 349-368. DOI: 10.1080/03091921003759876.

3. Гершуни Г.З., Жуковицкий Е.М. Конвективная неустойчивость несжимаемой жидкости. М.: Наука, 1986. C.7-12.

4. Sukhanovsky A., Batalov V., Teymurazov A., and Frick P. Horizontal rolls in convective flow above a partially heated surface // The European Physical Journal B-Condensed Matter and Complex Systems. 2012. Vol. 85, № 1.

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