Численное исследование конвективных течений в пакете ANSYS
Численное
исследование конвективных течений в пакете ansys
Содержание
Введение
. Теория
.1 Уравнения тепловой
конвекции
.2 Уравнения движения
.3 Элементы теории
завихренности
. Постановка задачи и метод
решения
.1 Свободная конвекция
.3 Вынужденная конвекция
.3 Рабочая среда ANSYS CFX
.3.1 ANSYS CFX и ANSYS
Workbench
.3.2 Предобработка в CFX-Pre
.3.3 Решатель ANSYS
CFX-Solver
.3.4 Постобработка в ANSYS
CFD-Post
.4 Рабочая среда ANSYS FLUENT
.4.1 Сеточный генератор
GAMBIT
. ANSYS CFX и ANSYS FLUENT
.1 Этапы решения задачи в
ANSYS CFX
.1.1 Создание геометрической
модели в ANSYS DesignModeler
.1.2 Создание сеточной модели
в ANSYS Meshing
.1.3 Предобработка в CFX-Pre
.1.4 Запуск решения в ANSYS
CFX-Solver Manager
.1.5 Постобработка в ANSYS
CFD-Post
.2 Этапы решения задачи в
ANSYS FLUENT
.2.1 Построение сеточной
модели в GAMBIT
.2.2 Построение физической
модели в ANSYS FLUENT
.3 Анализ численной схемы
. Результаты решения задачи
свободной конвекции
4.1 Качественное описание
формирования конвективного течения
.2 Количественное описание
формирования свободной конвекции
.2.1 Результаты CFX
.2.2 Результаты FLUENT
4.2.3 Сравнение результатов
CFX и FLUENT с результатами эксперимента
.3 Результаты решения задачи
вынужденной конвекции
Заключение
Литература
Приложения
Введение
Конвекция - это процесс переноса энергии потоками жидкости или газа.
Конвективные движения являются неотъемлемыми элементами многих природных
процессов, наблюдаемых в атмосфере и океанах Земли, а также течений,
реализуемых в различных технологических устройствах. Это определяет большой
интерес к экспериментальному и численному изучению конвективных процессов [4].
Однако изучение формирования подобных течений в эксперименте требуют
больших временных и материальных затрат. В отличие от эксперимента, численный
подход дает возможность варьировать ряд важных параметров задачи, таких как
вязкость, угловая скорость вращения модели, существенно влияющих на
формирование и поведение конвективных течений.
Цели данной работы:
· провести численное исследование конвективных течений в
программном комплексе ANSYS,
формирующихся вследствие локализованного нагрева в цилиндрическом слое
жидкости. Сравнить результаты расчетов в CFX и FLUENT
для различных режимов течения, сравнить эти результаты с данными, полученными в
эксперименте.
· провести аналогичные расчеты для конвективных течений,
формирующихся вследствие стока через центральное отверстие в неподвижном
цилиндрическом слое жидкости. Интерес к расчетам в неподвижном слое обусловлен
появлением экспериментальных работ [5-7], в которых в близкой постановке
пороговым образом происходило формирование вихря в области стока. Проведение
трехмерных расчетов требует значительных вычислительных ресурсов, поэтому
расчеты провели в осесимметричной постановке. Провести анализ численной схемы,
построить численную модель в программных пакетах ANSYS CFX и ANSYS FLUENT, выявить
наличие вихревых движений в неподвижной численных моделях.
1.
Теория
1.1
Уравнения тепловой конвекции
Жидкость может находиться в механическом равновесии (т.е. в ней может
отсутствовать макроскопическое движение), не находясь при этом в тепловом
равновесии. Уравнение Эйлера для покоящейся жидкости, находящейся в однородном
поле тяжести
- может быть удовлетворено и при непостоянной температуре в
жидкости. При этом, однако, возникает вопрос о том, будет ли такое равновесие
устойчивым. Оказывается, что равновесие будет устойчивым лишь при выполнении
определенного условия. Если это условие не выполняется, то равновесие
неустойчиво, что в результате развития возмущений приводит к появлению в
жидкости беспорядочных течений, стремящихся перемешать жидкость так, чтобы в
ней установилась постоянная температура. Такое движение носит название тепловой
конвекции. Условие устойчивости механического равновесия является,
другими словами, условием отсутствия конвекции. Конвекция наступает при
нарушении следующего условия:
,
где 
- удельная теплоемкость при постоянном давлении,
- температурный коэффициент расширения;
т.е. если температура падает по направлению снизу вверх, причем ее
градиент превышает по абсолютной величине указанное в условии значение.
Выведем уравнения, описывающие конвекцию [1]. Мы будем рассматривать
жидкость как несжимаемую. Это значит, что давление предполагается достаточно
мало меняющимся вдоль жидкости, так что изменение плотности под влиянием
изменения давления можно пренебречь. Это существенно упрощает систему
уравнений, описывающих конвекцию. Что же касается изменения плотности благодаря
неравномерной нагретости жидкости, то этим изменением, конечно, нельзя
пренебречь. Именно оно приводит к появлению сил, вызывающих конвекционное
движение.
Модель, которая описывает конвекцию жидкости, как несжимаемую, называется
моделью Буссинеска [2]. Соответствующие приближенные уравнения обычно
называют уравнениями конвекции в приближении Буссинеска.
В таком случае переменная температура записывается в виде 
, где 
есть некоторое постоянное среднее
значение, от которого отсчитывается неравномерность температуры 
. Будем предполагать, что 
мало по сравнению с 
. Плотность жидкости напишем в виде 
с постоянным
.
Ввиду малости изменения температуры мало также и вызываемое им изменение
плотности 
, причем можно написать:
, (1.1)
где 
- температурный коэффициент расширения.
В давлении же 
величина 
не будет постоянной.
Это давление, соответствующее механическому равновесию при постоянных
(равных и 
) температуре и плотности.
Оно меняется с высотой согласно гидростатическому уравнению
,
где координата z
отсчитывается вертикально вверх.
Учитывая все эти соображения, получаем полную систему уравнений,
описывающих термогравитационную конвекцию несжимаемой жидкости в приближении
Буссинеска:
Уравнение Навье-Стокса:
. (1.2)
Уравнение теплопроводности:
.
(1.3)
Уравнение непрерывности:
, (1.4)
где 
- вектор скорости, 
- изменение давления, 
- изменение температуры, 
- коэффициент кинематической
вязкости, 
- гидростатическое значение плотности, 
- коэффициент
температуропроводности, 
- коэффициент теплопроводности, 
- удельная теплоемкость при
постоянном давлении.
Возвращаясь к допущениям, сделанным при выводе уравнений (1.2) - (1.4),
отметим, что основным моментом в приближении Буссинеска является предположение
о том, что рассматривается в некотором смысле «слабая» конвекция: вызванные
неоднородностью температуры отклонения плотности от среднего значения
предполагаются настолько малыми, что ими можно пренебречь во всех уравнениях,
кроме уравнения движения (1.2), где отклонение учитывается лишь в члене с
подъемной силой. Разумеется, учет неоднородности плотности лишь в уравнении движения
означает некоторую непоследовательность приближения Буссинеска. Однако
сравнение результатов решения уравнений конвекции (1.2) - (1.4) с обширным
экспериментальным материалом с определенностью свидетельствует о том, что эти
уравнения достаточно хорошо отражают все важнейшие особенности тепловой
конвекции в лабораторных масштабах.
В частных случаях температура или тепловой поток могут быть заданы
непосредственно на границах полости, которую заполняет жидкость. При этом
уравнения и граничные условия будут содержать следующие параметры: характерную
длину полости L, характерную
разность температур Θ, время τ, характеризующее нестационарность
внешних условий, и параметры жидкости ν, χ, gβ.
Из этих величин можно
построить 2 независимые безразмерные комбинации:
так называемые числа Рэлея (R), Прандтля (P). Число Прандтля
зависит только от свойств самого вещества жидкости; основной же характеристикой
конвекции как таковой является число Рэлея.
Два течения подобны, если их числа R и P
одинаковы. Теплопередачу при конвекции в поле тяжести характеризуют числом
Нуссельта:
,
где 
- коэффициент теплопередачи между твердыми телами и
жидкостью, 
- коэффициент теплопроводности среды, l - характерный размер.
Конвективное движение может быть как ламинарным, так и турбулентным.
Наступление турбулентности определяется числом Рэлея - конвекция становится
турбулентной при очень больших значениях R.
Ламинарное течение - течение, при котором жидкость или газ перемещается слоями
без перемешивания и пульсаций (то есть беспорядочных и быстрых изменений
скорости и давления).
Турбулентное течение - явление, заключающееся в том, что при увеличении скорости
течения жидкости или газа в среде самопроизвольно образуются многочисленные
нелинейные фрактальные волны и обычные, линейные различных размеров, без
наличия внешних, случайных, возмущающих среду сил и/или при их присутствии.
1.2
Уравнения движения
Рассмотрим движения несжимаемой вязкой жидкости с неизменными физическими
свойствами. Уравнения, описывающие движение несжимаемой вязкой жидкости,
выражают сохранение массы и количества движения:
(2.1)
(2.2)
Здесь
- скорость частицы, измеряемая в системе координат,
вращающейся с постоянной угловой скоростью 
; 
представляют соответственно радиус-вектор частицы,
время, давление, плотность, кинематическую вязкость и массовую силу, отнесенную
к единице массы. Массовая сила предполагается консервативной 
, так что ее вместе с центробежной силой и Р можно
записать в форме редуцированного давления
(2.3)
Это
упрощает уравнение (2.2)
(2.4)
Полная
форма конвективного ускорения 
используется
чаще, чем инвариантное векторное представление.
Уравнения
движения в инерциальной системе координат получаются из предыдущих уравнений,
если положить в них 
. Формула
связывает
между собой скорости частицы в инерциальной и вращающейся системах.
На
твердых непроницаемых поверхностях вязкая жидкость должна двигаться вместе с
этими поверхностями, так как скольжение вдоль них или пересечение их
невозможно. Если поверхность проницаемая, можно задать нормальную компоненту
скорости, но требование отсутствия скольжения или относительной тангенциальной
скорости остается в силе. В общем случае граничная поверхность имеет части,
равномерно вращающиеся с угловой скоростью 
. По
отношению к системе, вращающейся со скоростью 
,
граничное условие на такой поверхности есть
(2.5)
Постановка
задачи завершается описанием начального поля скорости
(2.6)
Задача,
таким образом, состоит в том, чтобы решить уравнения (2.1) и (2.4) в фиксированной
области с граничными условиями (2.5) и (2.6).
Пусть
характеризуют типичную длину, время и относительную
скорость движения частицы. Замена переменных 
их
нормированными значениями 
позволяет привести уравнения к безразмерному виду:
(2.7)
(2.8)
с соответствующими граничными условиями. (Значок Λ
отмечает единичный
вектор.) При этом появляются два важных безразмерных параметра: число Экмана
(2.9)
и
число Россби
(2.10)
Первый является грубой мерой отношения типичной силы вязкости к силе
Кориолиса и есть, по существу, обратное число Рейнольдса. Подобно этом число
Россби - отношение конвективного ускорения к ускорению Кориолиса - дает общую
оценку относительного значения нелинейных членов. Число Экмана очень мало в
большинстве тех случаев, когда преимущественно проявляются эффекты вращения.
Практически величина 10-5 является для него обычной и в последующем
предположение Е << 1 используется без дальнейших оговорок. Число Россби
имеет порядок единицы или меньше; в линейной теории его значение принимается
бесконечно малым.
Безразмерное уравнение вихря во вращающейся системе есть
(2.11)
Где
(2.12)
Тривиальное
решение этих уравнений, 
= 0, отвечает, конечно, нетривиальному состоянию
твердого вращения, наблюдаемому во вращающейся системе. В инерциальной системе
соответствующая размерная скорость есть просто 
.
Очевидно, что вязкая жидкость, заключенная в закрытый равномерно вращающийся
сосуд стремится с течением времени к этому естественному состоянию жесткого
вращения.
1.3
Элементы теории завихренности
Нужно
будет различать абсолютный вихрь 
,
измеренный в инерциальной системе, и относительный вихрь 
, измеренный в равномерно вращающейся системе. Они
связаны соотношением
(3.1)
Такое
же различие устанавливается для циркуляции по замкнутому контуру
Некоторое
преобразование последнего интеграла приводит к форме
(3.2)
Здесь
- проекция поверхности 
,
ограниченной контуром L, на плоскость, перпендикулярную вектору 
; 
-
единичный вектор, нормальный к .
Применение
теоремы Стокса дает равносильное определение циркуляции
(3.3)
Завихренность
в некоторой точке пропорциональна мгновенному угловому моменту сферического
элемента жидкости (жесткого) в этой точке. Разумеется, мгновенная угловая
скорость частицы есть как раз 
. Линия в
жидкости, повсюду касательная к ,
называется вихревой линией; вихревые линии, проходящие через каждую точку малой
замкнутой кривой, образуют вихревую трубку. Если площадь поперечного сечения
трубки 
мала, величина 
имеет
одно и то же значение повсюду вдоль трубки и называется ее интенсивностью. Это
следует из теоремы о дивергенции
(3.4)
примененной
к объему вихревой трубки, заключенному между сечениями 
и 
.
Согласно (3.3), этот результат равносилен утверждению о том, что циркуляция по
некоторому контуру, опоясывающему боковую поверхность трубки, постоянна. Другое
заключение, вытекающее из (3.4), состоит в том, что вихревые линии не могут
начинаться или оканчиваться в жидкости, они либо замкнуты, либо оканчиваются на
твердых границах.
Структуру
поля завихренности можно изобразить с помощью вихревых трубок равной
интенсивности, пронизывающих среду. Положение каждой трубки задается ее
центральной вихревой линией. Плотность трубок в некотором месте пропорциональна
завихренности, возрастающей по мере сближения трубок. При растяжении вихревой
трубки постоянной интенсивности завихренность возрастает, компенсируя
уменьшение площади поперечного сечения.
Аккуратная
аппроксимация уравнения завихренности в инерциальном пространстве
(3.5)
дает
соотношение между конвективной скоростью изменения абсолютной завихренности
(левая часть) и тремя процессами, вызывающими это изменение (правая часть). Эти
процессы в том порядке, как они записаны, есть растяжение и поворот вихревых
линий, возникновение завихренности за счет изменения плотности и диффузионный
перенос завихренности от соседних элементов.
Ниже
кратко сформулированы несколько классических теорем о завихренности невязких
жидкостей.
Теорема
Лагранжа утверждает, что если в
невязкой среде с постоянной плотностью поле завихренности равно нулю в
начальный момент, то оно остается нулевым всегда, 
. (Условия, разумеется, исключают все способы
возникновения завихренности.)
Теорема
Кельвина вытекает из формулы для
полного изменения циркуляции по замкнутому контуру, движущемуся вместе с жидкостью:
(3.6)
Следовательно,
если жидкость однородна или баротропна Р = Р (ρ) и диссипация отсутствует, то циркуляция по контуру,
движущемуся вместе с жидкостью, остается все время постоянной. Теорема
Гельмгольца следует прямо отсюда и утверждает, что при тех же условиях
вихревая трубка движется вместе с жидкостью и ее интенсивность остается
постоянной. Эти важные теоремы, будучи записаны для вращающейся системы
координат, показывают, как изменяются относительная завихренность и циркуляция
под действием основного вращения.
2.
Постановка задачи и метод решения
2.1
Свободная конвекция
Рассматривается цилиндрический слой жидкости. Полагается, что на верхней
границе выполняется условие проскальзывания, на боковой стенке и на дне -
условие прилипания (схематично исследуемая модель и циркуляция жидкости
представлены на рисунке 1)
Рис. 1. Схематичное изображение расчетной области и меридиональной циркуляции
В центральной части цилиндр нагревается за счет постоянного теплового
потока; область нагрева показана на схеме черным прямоугольником. Отвод тепла
осуществляется через свободную верхнюю границу, на которой, кроме
проскальзывания, также задается постоянный поток тепла (отрицательный).
Количество вводимого в слой тепла задается равным количеству тепла, отводимому
с верхней границы.
Граничные условия:
Полагается, что на боковых стенках и дне выполняется условие прилипания
, (4)
На верхней границе - условие проскальзывания:
, 
. (5)
На верхней границе и в области нагрева также задаются потоки тепла:
, 
. (6)
2.2
Вынужденная конвекция
Рассматривается цилиндрический слой жидкости. Полагается, что на боковой
стенке и на дне выполняется условие прилипания (схематично исследуемая модель
представлена на рисунке 3)
Рис. 2. Схематичное изображение расчетной области
Граничные условия:
Полагается, что на боковых стенках и дне без стока выполняется условие
прилипания
, (7)
На верхней границе и в области стока задаются значения скорости:
, 
(8)
Исследуемые
задачи решаются в осесимметричной постановке (не учитывается зависимость от азимутальной
координаты φ),
течение жидкости предполагается
ламинарным и описывается системой управляющих уравнений в размерной
формулировке, включающей:
- уравнение сохранения полной энергии
(9)
уравнение
количества движения
; (10)
- уравнение непрерывности
; (11)
В
уравнениях (9-11) используются обозначения: 
- полная
энтальпия при наличии вращения; ω - угловая скорость вращения жидкости; r -
радиус-вектор; λ
- коэффициент теплопроводности; ρ - плотность жидкости; p -
давление; U - вектор скорости; T - температура; τ - тензор вязких напряжений; 
- член, учитывающий силу Кориолиса и центробежную
силу; 
- источник импульса, образованный силами плавучести; β - коэффициент объемного расширения; g - вектор
гравитационного ускорения; 
-
источник энергии.
Для численной реализации задачи в описанной постановке использовались
конечно-элементные программные комплексы ANSYS CFX и ANSYS FLUENT. Принятые допущения об осевой
симметрии задачи и ламинарном характере течения позволили значительно сократить
количество расчетных узлов и время вычислений.
2.3
Рабочая среда ANSYS CFX
ANSYS CFX - мощный
инструмент для оптимизации процесса разработки и технологической подготовки в
области вычислительной динамики жидкостей и газов. ANSYS CFX совмещает в себе
передовую технологию решателя с современным пользовательским интерфейсом и
адаптивной архитектурой, что делает этот инструмент доступным как для
разработчиков, владеющих общими инженерными знаниями, так и для специалистов в
области гидродинамики, работающими с моделью и ее свойствами на глубоком
уровне. CFX позволяет детально изучить оборудование и процессы изнутри,
повысить эффективность, увеличить срок службы и оптимизировать процессы.
2.3.1
ANSYS CFX и ANSYS Workbench
Модуль ANSYS CFX полностью интегрирован в расчетную среду ANSYS Workbench
- платформу, объединяющую все инструменты инженерного моделирования компании
ANSYS. Адаптивная архитектура позволяет пользователю выполнять любые действия
от стандартного анализа течения жидкости или газа до обработки сложных
взаимодействующих систем. Пользователи могут легко оценить производительность
во множественных расчетных точках или сравнить несколько альтернативных
конструкций. Для решения задач из различных расчетных дисциплин в рамках
платформы ANSYS Workbench можно получить доступ к общим для всех расчетов
инструментам, таким как инструменты для работы с геометрией ANSYS DesignModeler
и с сеткой ANSYS Meshing.
Геометрия: модуль ANSYS DesignModeler специально разработан для создания
геометрии, ее подготовки к моделированию. Эта полностью параметрическая среда
является простой в использовании, обладает прямыми и двусторонними связями со
всеми основными CAD-пакетами и используется как средство согласования и ввода
геометрии во все программные продукты ANSYS для последующих инженерных
расчетов.
Создание сетки: для получения точных результатов расчетов в области
вычислительной динамики жидкостей и газов необходима самая совершенная
технология создания сетки. ANSYS Meshing предоставляет множество технологий
построения сетки в одном приложении. Это дает возможность выбрать наиболее
подходящий способ построения сетки на основе применения подхода
последовательного построения сетки на каждой из деталей сборки. конвективный течение нагрев жидкость
2.3.2
Предобработка в CFX-Pre
Физический препроцессор ANSYS CFX - это современный, интуитивно понятный
интерфейс для подготовки к моделированию динамики жидкостей или газов. В
дополнение к основному режиму работы встроенный мастер помогает пользователю
пройти процесс подготовки общих расчетов течения жидкостей или газов.
Мощный язык программирования дает возможность задать описание проблемы в
числовом виде, как в случае со сложными граничными условиями, авторскими
моделями материалов или дополнительными уравнениями переноса. Адаптивная
архитектура CFX-Pre также позволяет пользователю создавать собственные панели
графического пользовательского интерфейса, чтобы стандартизировать ввод для
выбранных приложений.
2.3.3
Решатель ANSYS CFX-Solver
Сердцем модуля ANSYS CFX является передовой алгебраический многосеточный
сопряженный решатель, использующий технологию Coupled Algebraic Multigrid,
являющуюся ключом к получению точных результатов в короткие сроки. Параметры
решателя, граничные условия могут быть скорректированы во время выполнения
расчета, при этом нет необходимости останавливать решатель. Решатель ANSYS CFX
использует схему дискретизации второго порядка по умолчанию, обеспечивая
получение максимально точных результатов.
Использование технологии сопряженных решателей ANSYS CFX дает
значительные преимущества при проведении любого расчета и позволяет получить
устойчивые и масштабируемые решения для задач динамики жидкостей и газов.
2.3.4
Постобработка в ANSYS CFD-Post
Постпроцессинг результатов графического и количественного анализа,
полученных в ANSYS CFX, выполняется в ANSYS CFD-Post. Возможность написания
сценариев, полная автоматизированность, генерация отчетов дают пользователю
максимальное количество информации, полученной в результате расчетов.
2.4
Рабочая среда ANSYS FLUENT
Программный модуль ANSYS FLUENT имеет широкий спектр возможностей
моделирования течений жидкостей и газов для промышленных задач с учетом
турбулентности, теплообмена, химических реакций. К примерам применения FLUENT
можно отнести задачи обтекания крыла, горение в печах, внешнее обтекание
нефтедобывающих платформ, течение в кровеносной системе, конвективное
охлаждение сборки полупроводника, вентиляция в помещениях, моделирование
промышленных стоков. Специализированные модели горения, аэроакустики,
вращающихся/неподвижных расчетных областей, многофазных течений серьезно
расширяют области применения базового продукта.FLUENT - это удобный,
отказоустойчивый инструмент, позволяющий даже новичкам достигать высокой
производительности труда. Интеграция модуля ANSYS FLUENT в рабочую среду ANSYS
Workbench, а также возможность использования модуля ANSYS CFD-Post для
обработки результатов создает комплексное решение для выполнения инженерного
анализа в области моделирования течений жидкостей и газов.
Технология ANSYS FLUENT является лидером по числу сложных физических
моделей, предлагаемых для расчетов на неструктурированных сетках. Представлены
наборы элементов различных форм: четырехугольники и треугольники для двумерных
расчетов, гексаэдры, тетраэдры, полиэдры, призмы, пирамиды для трехмерных
расчетов.
Сетки можно строить при помощи инструментов компании ANSYS или при помощи
инструментов сторонних производителей.FLUENT включает в себя сложные числовые
отказоустойчивые решатели, в том числе сопряженный решатель с алгоритмом на
основе давления, расщепленный решатель с алгоритмом на основе давления и два
решателя с алгоритмами на основе плотности, что позволяет получать точные
результаты для практически любого течения.
2.4.1 Сеточный генератор GAMBIT
FLUENT предлагает универсальное решение для создания расчетных сеток
любых типов - набор программных продуктов, состоящий из GAMBIT, Tgrid и
G/Turbo. Последний является своеобразным аналогом TurboGrid, но значительно
уступает ему по функциональным возможностям. Отметим, что в отличие от ANSYS
CFX, FLUENT поддерживает некомфорные сетки, сгенерированные в GAMBIT.
Качество получаемых на основе проведения вычислительного эксперимента
результатов напрямую зависит от качества построенной расчетной сетки.
Предпроцессор GAMBIT позволяет быстро создавать и обрабатывать геометрии
исследуемых процессов. GAMBIT имеет единый интерфейс для создания
геометрических моделей и построения сетки. Кроме того, в процессе его
использования производится запись всех операций построения в текстовый файл
(журнал событий), что дает возможность пользователю легко отслеживать путь
создания геометрической модели, исправлять ошибки, перестраивать модели и
сетки, а также интегрировать результаты в другие программы.может импортировать
геометрию из различных CAD программ. Часто во время импорта геометрии возникают
проблемы нестыковки объектов. GAMBIT имеет встроенный модуль поиска и
исправления подобных проблем, что значительно упрощает процесс создания модели.обладает
мощными возможностями для создания двухмерных и трехмерных расчетных областей
непосредственно внутри программы, начиная от построения линий и сплайнов,
кончая созданием твердотельных объектов. Встроенный специализированный модуль
G/Turbo позволяет автоматизировать создание геометрических моделей лопаточных
машин.обладает мощным генератором сеток, позволяющим создавать разнообразные
типы сеток: структурированную гексаэдальную сетку, автоматическую
(неструктурированную) гексаэдальную и тетраэдальную сетки. Кроме того, в нем
имеется возможность создания пограничных слоев с комбинированными сетками.
После построения сетки пользователь имеет возможность проверить ее качество по
разнообразным параметрам (скошенность элементов, соотношение сторон).
3. ANSYS CFX и ANSYS FLUENT
CFX и ANSYS FLUENT - это основные продукты для задач гидрогазодинамики
общего назначения, предлагаемые компанией ANSYS, Inc.
Оба решателя разрабатывались в течение десятилетий независимо друг от
друга и обладают несколькими существенными отличиями, несмотря на некоторые
схожие черты. Оба модуля основаны на методе контрольных объемов, дающем высокую
точность, и используют решатель по давлению, что позволяет применять эти
продукты для решения широкого круга инженерных задач. Основные отличия состоят
в способе интегрирования уравнений течения жидкостей и в стратегиях решения
уравнений.
Решатель ANSYS CFX использует сетку конечных элементов (числовые значения
в узлах сетки), схожую с теми, что используется в анализе прочности, для
дискретизации области. В отличие от ANSYS CFX, решатель ANSYS FLUENT использует
сетку конечных объемов (числовые значения в центрах ячеек). В итоге оба подхода
формируют уравнения для конечных объемов, которые обеспечивают сохранение
значений потока, что является необходимым условием для точных решений задач
гидрогазодинамики. В ANSYS CFX особый упор сделан на решение основных уравнений
движения (сопряженная алгебраическая сетка), а ANSYS FLUENT предлагает
несколько подходов к решению (метод на основе плотности, расщепленный метод на
основе давления, сопряженный метод на основе давления). Оба решателя содержат в
себе самые ценные возможности физического моделирования для получения
максимально точных результатов.
3.1
Этапы решения задачи в ANSYS CFX
Запускаем программную платформу ANSYS Workbench. В левой
панели инструментов «Системы анализа» (Analysis Systems) выбираем анализ Fluid Flow (CFX).
В рабочей области «Project
Schematic» появится модуль в виде структурной
схемы, в которой каждому этапу соответствует раздел, содержащий объекты
расчетной модели:
Для каждого объекта возможен ввод и редактирование свойств.
Рис. 4.Схематический вид готового проекта в среде ANSYS Workbench
3.1.1
Создание геометрической модели в ANSYS DesignModeler
Для создания новой геометрической модели нужно нажать правой кнопкой мыши
«Geometry» в структурной схеме модуля Fluid Flow (CFX)
(Рис. 3). Появится контекстное меню, которое позволяет создать новую геометрию
«New Geometry», импортировать данные геометрии,
передать существующую геометрии в другие системы анализа для работы в исходном
проекте, стереть данные, просмотреть свойства геометрической модели. Выбрав «New Geometry» в контекстном меню откроется
программа для работы с геометрией ANSYS DesignModeler.
После открытия DesignModeler сразу появится диалоговое окно с выбором
единицы измерения длины. Расчет будем проводить в миллиметрах.
Графический интерфейс программы состоит из окон с графиком
модели 3D «Graphics», со схемой модели «Tree Outline», с детальным обзором функций «Details View», и с разнообразными панелями инструментов.
Аналогична организация графического интерфейса других программ, которых мы
будем использовать в дальнейшем.
Рис. 5.Схематический вид начального проекта в ANSYS DesignModeler
Создание геометрической модели начнем с построения точек. Для этого в
главной панели нажимаем Create -> Point 
. В нижнем левом окне «Details View» появятся функции для детального задания точки (Рис.
6). Выбираем тип точки «Construction Point», задание
точки через координаты «Manual Input». Вводим
значения координат в соответствующие поля. Для добавления новой точки следует
нажать правой кнопкой мыши и выбрать «Add New Point Group». После того, как ввели все нужные точки нужно нажать
кнопку «Generate» 
. Чтобы сгенерировались все
введенные точки. Это следует делать всякий раз при добавлении нового
геометрического элемента.
Далее нужно создать прямые линии для соединения точек. Для этого в
главной панели нажимаем Concept ->Lines From Point 
.
Удерживая мышкой клавишу «Ctrl»,
соединяем нужные точки в графической области. В окне «Details View» в поле «Operation» выбрать «Frozen».
После этого строим поверхность через созданные линии. Для этого в главной
панели нажимаем Concept->Surfaces From Edges 
. Удерживая мышкой клавишу «Ctrl»,
выбираем линии в графической области.
Теперь нужно сделать копию этой поверхности, повернутую относительно
вертикальной оси Y. Для этого в
главной панели нажимаем Create ->Body Operation 
. В окне «Details View» выбираем поверхность, которую хотим повернуть, тип
операции «Rotate», сохранение исходной поверхности,
ось, относительно которой поворот, и угол поворота.
Так как наша задача осесимметричная, угол берется минимальным для
уменьшения времени расчета, путем уменьшения количества узлов в сеточной
модели.
Чтобы связать поверхности, аналогичным образом используем операции и .
Для создания объема воспользуемся функцией . В окне «Details View» выбираем все поверхности для «сшивки», тип операции
«Sew», слияние всех поверхностей «Merge Bodies», «Create Solids».
В итоге получили геометрическую модель задачи (Рис. 8).
Следует заметить, что наши поверхности не проходят через центр координат,
мы берем отступ в 1 мм. Причина этого в том, что решатель выдает плохие
результаты в центре координат.
Закрываем DesignModeler, все данные по геометрии автоматически
сохраняются в проекте Workbench. Можно переходить к следующему этапу - создание сеточной модели.
.1.2
Создание сеточной модели в ANSYS Meshing
Для создания новой сеточной модели следует выбрать «Mesh» в структурной схеме модуля Fluid Flow (CFX)
(Рис. 3). В открывшейся программе Meshing в «Outline» выбрать «Mesh» - в нижней части программы «Details of Mesh» появятся опции для работы с сеткой.
Инструменты для создания сетки позволяют генерировать сеточные модели для
разных типов анализа. Во вкладке «Defaults» выбираем метод вычислительной гидрогазодинамики «CFD», решатель «CFX».
Качество сеточной модели влияет на точность, сходимость и скорость
получения решения. Пространственное разрешение возьмем 1 мм. Тестирование
показало, что дальнейшее сгущение расчетной сетки не влияет на результаты. ANSYS Meshing позволяет задавать особые настройки сетки в требуемой
зоне модели используя команду Mesh Control -> Sizing 
. Для нас же представляет наибольший интерес поверхности
симметричные относительно вертикальной оси. Во вкладке «Sizing» выбираем минимальное и максимальное
разрешение сетки «Min Size», «Max Size». Для применения настроек сетки нужно нажать кнопку «Update» 
. Во вкладке «Statistics» можно посмотреть полученные числа
узлов и элементов.
Закрываем ANSYS Meshing, переходим к следующему этапу
CFX-Pre.
3.1.3
Предобработка в CFX-Pre
ANSYS CFX-Pre реализует процесс определения физики задачи. Физический
препроцессор импортирует сетку, созданную на предыдущем шаге. Это следующий шаг
постановки задачи, на котором определяются физические модели, на основе которых
будет происходить симуляция процесса, а также их основные параметры и
характеристики. CFX-Pre позволяет определить начальные и граничные условия
процесса (входные, выходные параметры), модели теплообмена.
Чтобы открыть программу, следует выбрать «Setup» в структурной схеме модуля Fluid Flow (CFX)
(Рис. 3)
Параметры расчетной области
Чтобы задать параметры для всей расчетной области нужно создать домен,
выбрав Insert -> Domain
в панели главного меню. Откроется
окно с деталями.
Во вкладке «Basic Settings» выбираем нашу расчетную область в
поле «Location» (Рис. 11). Задаем тип расчетной
области - жидкость «Fluid Domain». Вещество - вода («Material» - «Water»). Давление всей области «Preference Pressure» - 1 atm. Указываем, что конвекция есть «Buoyancy Model» - «Buoyant»,
т.е. задаем ускорение свободного падения против оси Y (Gravity Y Dirn. =
-9.8). Buoy.Ref. Temp=20 [C].
Если требуется задать вращение модели в поле «Domain Motion» следует выбрать «Rotating», ввести угловую скорость вращения и
ось, относительно которой вращается модель.
Переходим во вкладку «Fluid Models». Во
вкладке «Heat Transfer» выбираем изотермический процесс «Isothermal». Рассматриваем ламинарное течение «Laminar», без окисления и радиации.
Для определения начальных условий переходим во вкладку «Initialization». Определяем начальные значения для
компонентов вектора скорости, относительное давление «Relative Pressure» выбираем равным нулю, начальную
температуру жидкости - 20 [C].
Нажимаем кнопку ОК.
Для использования другого вещества (масла) в качестве жидкости следует
выбрать Insert -> Material в панели главного меню 
.
В появившемся окне деталей во вкладке «Basic Settings» выбираем тип вещества - чистое
вещество «Pure Substance», группу жидкостей с постоянными
свойствами «Material Group» - «Constant Property
Liquids», термодинамическое состояние «Thermodynamic State» выбираем жидким «Liquid». Переходим во вкладку «Material Properties», чтобы задать свойства материала
(Рис. 12). В поле «Equation
of State» задаем значение молярной массы и плотности вещества.
В поле «Specific Heat Capacity» - удельную теплоемкость. В поле «Reference State» - значения для температуры и давления. В поле «Transport Properties» - значение динамической вязкости «Dynamic Viscosity» и коэффициент теплопроводности «Thermal Conductivity». В «Buoyancy Properties» - коэффициент теплового расширения
«Thermal Expansivity».
Граничные условия
Для создания нового граничного условия нужно выбрать Insert -> Boundary 
. Появится окно для задания
параметров граничного условия. Во вкладке «Basic Settings» выбираем тип «Boundary Type» и область граничного условия. Во вкладке «Boundary Details» выбираем параметры для массы и моментов потока
жидкости «Mass And Momentum» и тип теплопередачи «Heat Transfer».
Тепловая конвекция:
В области с нагревателем задается условие прилипания и положительный поток
тепла «Heat Flux» [W/m2] - мощность нагрева 30 Вт, поделенная на площадь
нагреваемой области.
Начнем с того, что зададим граничное условие на внешней боковой стенке.
Выбираем область и задаем тип граничного условия - стенка «Wall». Условие прилипания - «No Slip Wall», нулевой поток тепла - «Adiabatic».
На внутренней боковой стенке и на части дна без нагревателя такие же
граничные условия.
На верхней области задается условие проскальзывания «Free Slip Wall» и отрицательный поток тепла «Heat Flux» [W/m2] - мощность нагрева 30 Вт, поделенная на площадь
охлаждаемой области.
Вынужденная конвекция:
Зададим граничное условие на внешней боковой стенке. Выбираем область и
задаем тип граничного условия - стенка «Wall». Вращающиеся стенки «Frame Type» - «Rotating». Условие прилипания - «No Slip Wall».
На внутренней боковой стенке и на части дна вне отверстия такие же
граничные условия.
В области стока жидкости в центральной части модели задается тип
граничных условий - «Inlet».
Задаем вращение - «Rotating». И
указываем величину скорости, нормальной к области.
В области истока жидкости в верхней части модели задается тип граничных
условий - «Inlet». Задаем вращение - «Rotating». И указываем величину скорости,
нормальной к области.
Для реализации осесимметричной постановки задачи надо задать граничные
условия на поверхностях, симметричных относительно вертикальной оси. Нужно
выбрать Insert -> Domain Interface 
.
В открывшемся окне выбираем тип условия «Interface Type» - «Fluid Fluid». Выбираем
соответствующие поверхности. Для задания симметрии относительно оси, выбираем
опцию «Rotational Periodicity» и ось симметрии Y.
В итоге получили физическую модель нашей задачи (Рис. 13).
CFX-Pre также позволяет настроить параметры
для вывода данных. Для этого следует перейти в окне «Outline» в «Output Control» 
. В случае нестационарной задачи
панель управления выводом позволяет управлять, какие данные будут записываться
в файлы вывода и с каким промежутком времени. Результаты могут быть написаны на
особых стадиях решения при записи резервных файлов после заданного числа
итераций. Данные мониторинга позволяют также показывать процесс решения задачи
в реальном времени.
Для записи результатов в определенные моменты времени следует перейти во
вкладку «Trn Results», добавить новый элемент (item), и выбрать нужные параметры для записи «Output Variables List». Также следует указать интервал времени между
записями «Time Interval». Нас интересуют компоненты
скорости, и полная температура с интервалом в 1 секунду. Закрываем CFX-pre.
3.1.4
Запуск решения в ANSYS CFX-Solver Manager
CFX-Solver Manager представляет собой графический интерфейс
пользователя, который позволяет задавать параметры для вычислений: управлять
процессом решения CFX-Solver в интерактивном режиме, определять
входные данные файла решателя, запускать или приостанавливать CFX-Solver, контролировать процесс решения задачи, устанавливать
решатель для проведения параллельных вычислений.
Чтобы открыть его, в окне Workbench следует выбрать «Solution» в структурной схеме модуля Fluid Flow (CFX) (Рис. 3). В появившемся диалоговом
окне можно настроить некоторые параметры решения задачи и запустить CFX-Solver Manger.
После завершения вычислений будет записан общий файл результатов и файлы
интересующих нас параметров (компоненты скорости и полная температура) для
заданных моментов времени. Это основные данные, которые будут использованы и
обработаны в CFD-Post.
CFX-Solver Manger может интерполировать решение данных текущей сеточной
модели на новую сеточную модель. Также позволяет экспортировать файл
результатов в другие форматы данных.
Можно закрывать CFX-Solver Manger и переходить к обработке данных в CFD-Post .
3.1.5
Постобработка в ANSYS CFD-Post
ANSYS CFD-Post - это программа,
предназначенная для анализа, визуализации и представления результатов,
полученных в ходе решения задачи посредством ANSYS CFX-Solver. Для этого
используются следующие средства:
· визуализация геометрии и исследуемых областей;
· векторные графики для визуализации направления и величины
потоков;
· визуализация изменения скалярных величин (такие как
температура, давление) внутри исследуемой области.
Графики, изображения и видео, полученные в результате анализа решения
задачи можно сохранить в виде отдельных файлов.
Чтобы открыть CFD-Post, в окне Workbench следует выбрать «Results» в структурной схеме модуля Fluid Flow (CFX)
(Рис. 3).
Для начала введем график линий уровня с использованием цветной
градиентной заливки «Contour Plot» на
интересующую нас поверхность, выбрав Insert -> Contour 
. С помощью «Contour Plot» можно будет посмотреть изменение какой-либо
физической величины в данной области.
В окне «Details» следует выбрать область и
переменную для построения графика. Так же можно выбрать число контуров
прорисовки для точности.
Для того чтобы посмотреть поле температуры в конкретный момент времени,
следует перейти в Tools -> Timestep Selector 
и выбрать момент времени из списка.
Для наблюдения за изменением процесса во времени есть функция,
позволяющая делать анимации. Надо перейти в Tools -> Animation 
. В открывшемся окне нужно выбрать
объекты для анимации «Timesteps», можно регулировать скорость анимации, количество просмотров, и
сохранить анимацию в различных форматах видео.
Также можно поменять настройки цветовой шкалы «Legend»: изменить расположение шкалы в окне просмотра «3D Viewer», выбрать переменную, формат текста и прочее. Post
позволяет построить график зависимости какой-либо величины от времени в
некоторой точке.
Чтобы добавить точку, нужно выбрать в панели главного меню Insert -> Location -> Point 
. В окне «Details», выбрав метод «Method» определения точки (через координаты «XYZ» или по номеру узла), ввести
соответствующие значения. Так же можно определить форму, цвет и размеры точки.
Чтобы добавить график, следует выбрать Insert -> Chart 
. В окне «Details» во вкладке «General» надо выбрать тип графика «XY - Transient» (т.е. зависимость от времени), добавить заголовок графика. В следующей
вкладке «Data Series» в поле «Location» выбрать точку наблюдения температуры. Во вкладке «Y Axis» в качестве переменной «Variable» выбрать интересующую нас величину.
Нажав кнопку «Apply», построится график зависимости от
времени в точке. Нажав кнопку «Export»,
можно экспортировать данные в точке в качестве таблицы значений.
Функция экспорта данных CFD-Post позволяет экспортировать значения
необходимых переменных в каждом узле в интересующей области в один определенный
момент времени. Формируется текстовый документ с колонками значений в каждом
узле. Для этого надо выбрать File -> Export
-> Export… 
. Откроется окно с параметрами экспорта (Рис. 17), где во
вкладке «Option» можно указать название и формат
файла вывода, интересующую область модели и требуемые переменные. Во вкладке «Formatting» можно выбрать точность значений, а
также определить тип разделителя для колонок в файле вывода.
Эта функция экспорта позволяет экспортировать файл данных только для
одного конкретного момента времени. Чтобы экспортировать файл данных для всех
моментов времени, следует использовать функцию сессий, иначе потребуется
вручную делать экспорт файлов данных для каждого момента времени.
Принцип работы сессий заключается в том, чтобы записать один сеанс работы
в проекте ANSYS CFD-Post в виде программного кода на
языке программирования CCL (ANSYS CFX Command Language) и использовать этот сеанс в цикле
данного проекта или в других подобных проектах в ANSYS.
Для начала следует создать файл сессии. Для этого нужно перейти в Session -> New Session… 
и сохранить сессию в формате CFD-Post Session «*.cse». Чтобы начать запись сеанса,
следует перейти Session
-> Start Recording
. С этого момента все действия в
CFD-Post, в том числе и экспорт данных, будут записываться в созданный файл
сессии. После проведения всех нужных действий для записи в файл сессии, нужно
остановить запись, т.е. перейти Session -> Stop Recording 
. Теперь можно просмотреть программный код записанной сессии,
открыв файл сессии «*.cse»
обычным текстовым редактором. Для изменения программного кода (добавления
цикла) следует учитывать особенности синтаксиса языка программирования CCL. Добавив команду цикла и сохранив
все изменения в программном коде файла сессии, следует «проиграть» эту сессию,
выбрав Session -> Play Session…
. В результате выполнятся все
команды, записанные в файл сессии, т.е. в цикле будут экспортированы файлы
данных для каждого момента времени.
Далее следует обработать полученные файлы данных, извлекая из них нужные
значения и высчитывая соответствующие выражения. Для этого можно
воспользоваться любым языком программирования, поддерживающим функцию чтения
текстового документа. Мы же в этой работе будем использовать систему Wolfram Mathematica, которая позволяет импортировать
полученные данные в свою систему в удобном формате, обрабатывать их и
экспортировать результаты в виде таблиц и графиков.
Программный код сессии с циклом и код программы, выполненный в системе Wolfram Mathematica, находятся в Приложениях.
3.2
Этапы решения задачи в ANSYS FLUENT
3.2.1
Построение сеточной модели в GAMBIT
Качество получаемых на основе проведения вычислительного эксперимента
результатов напрямую зависит от качества построенной расчетной сетки. GAMBIT
имеет единый интерфейс для создания геометрических моделей и построения сетки.
Запускаем программу GAMBIT.
Начнем с построения геометрической модели.
. Выберем кнопку «GEOMETRY COMMAND BUTTON»
.
1.1. Далее, ниже, для создания точек выберем кнопки: 
→ 
.
Надо последовательно ввести соответствующие координаты точек в поля x, y, z и нажимать
кнопку «Apply».
.2. Для создания прямых линий перейти - 
→ 
.
Нажать кнопку 
. В появившемся окне выбрать точки для их соединения и
перенести в правую часть. Нажать «Apply».
.3. Для создания плоскости следует перейти - 
→ 
.
Нажать кнопку . В появившемся окне выбрать кнопку «All->». Нажать «Apply».
. Выберем кнопку «MESH COMMAND BUTTON» 
.
2.1. Выбираем опции для создания сетки на плоскости 
→ 
.
Нажав кнопку , появляется окно, где надо выбрать нашу плоскость.
В поле «Interval Size» задается размер элемента. Нажать «Apply».
Задаем граничные условия на сторонах модели и тип среды.
. Выберем кнопку «ZONES COMMAND BUTTON» 
.
3.1. Для задания граничных условий следует выбрать кнопку 
.
Нажимая кнопку , в появившемся окне выберем соответствующую сторону. В поле
«Type» выбираем тип граничных условий.
Присваиваем имя для стороны в поле «Name». Нажимаем «Apply». Так
следует проделать для всех сторон. В случае тепловой конвекции четыре стороны
будут типа «WALL», одна сторона - ось симметрии -
типа «AXIS».
В случае вынужденной конвекции две стороны будут типа «WALL», одна сторона - ось симметрии -
типа «AXIS», две стороны типа «VEL-OUTLET».
.2. Для задания типа среды следует выбрать кнопку 
.
Нажимая кнопку , в появившемся окне выберем соответствующую плоскость «Faces». Выбрать тип - жидкость «FLUID». Выбираем имя для нашей среды в «Name». Нажать «Apply».
Так мы построили сеточную модель нашей задачи.
В ходе построения можно просмотреть нашу модель в разных плоскостях,
нажав кнопку 
в нижней части панели инструментов.
Получившуюся сеточную модель нужно экспортировать. Для этого перейдем в
главное меню File → Export → Mesh…
. В появившемся окне выбрать название файла, поставить метку «Export 2D Mesh». Нажать «Apply».
В процессе построения модели производится запись всех операций в текстовый
файл (журнал событий). Это позволяет использовать его как программный код. Для
сохранения журнала событий надо перейти File → View File. В
появившемся окне скопировать код в текстовый документ.
Программный код операций Gambit,
находится в Приложениях.
3.2.2
Построение физической модели в ANSYS FLUENT
Запускаем ANSYS FLUENT.
Для того чтобы открыть файл сетки, следует выбрать в главном меню File → Read → Mesh.
Слева расположено основное меню для построения физической модели.
) В разделе «General» основного
меню задаются параметры решателя. В колонке «Solver» выбираем тип решателя «Pressure-Based», в «Time»
нестационарный процесс «Transient», «Velocity Formulation» - «Absolute», осесимметричную задачу в «2D Space» - «Axisymmetric». Ставим метку «Gravity»,
задаем ускорение свободного падения против оси x.
) В разделе «Models» задаем
основные уравнения: «Energy»
- On, «Viscous» - Laminar.
) В «Materials» задаем свойство жидкости. Нажимаем
кнопку «Create». В появившемся окне задаем свойства
жидкости. Нажимаем кнопку «Create».
) В разделе «Cell Zone Conditions» для нашей среды выбираем тип «fluid». В «Edit» выбираем нашу жидкость.
Для того чтобы, задать вращение модели ставим галочку напротив «Frame motion». Задаем угловую скорость в поле «Rotational Velocity».
) В «Boundary Conditions» задаем граничные условия на наши
стороны.
Вынужденная конвекция: для этого нужно выбрать соответствующую сторону и
перейти по кнопке «Edit». Для зоны «Inlet» в поле «Momentum» следует выбрать метод задания скорости «Velocity Specification Method» -
«Magnitude, Normal to Boundary». Нам нужно задать относительную скорость стока, поэтому в поле
«Reference Frame» выбираем опцию «Relative to Adjacent Cell Zone». Теперь нужно задать величину скорости в «Velocity Magnitude».
Для зоны «Outlet» все также,
только нужно выбрать абсолютную скорость в «Reference Frame» - «Absolute».
Свободная конвекция: перейдя по кнопке «Edit», откроется окно, где можно выбрать условие
прилипания «No Slip», условие проскальзывания «Specified Shear», задать тепловой поток «Heat Flux» на соответствующих сторонах.
6) В разделе «Reference Value» задаются исходные величины
некоторых параметров: плотность, вязкость, температура. Так же следует выбрать
нашу зону в «Reference Zone».
) В «Solution Methods» выбираем опции решателя. Выбираем в «Pressure-Velocity Coupling» схему «Scheme» - «PISO». В «Spatial Discretization»
выбираем для давления «Pressure» - «Presto!». Это наиболее подходящие опции для
нашего течения.
) В «Solution Initialization» выбрать начальное значение
температуры. Следует нажать кнопку «Initialize», чтобы установить начальное состояние задачи.
) В «Calculation Activities» можно настроить экспорт данных по
ходу решения задачи, выбрать интересующие параметры для вывода, указать
директорию для записи.
) В «Run Calculation» можно выбрать адаптивный или
фиксированный типы шага, задать размер шага, число шагов. Для того чтобы
запустить решение нужно нажать кнопку «Calculate».
По окончанию решения появится окно с уведомлением «Calculation completed».
Полученные данные можно просмотреть и обработать в постпроцессоре «CFD-Post», перейдя File →
Export to CFD-Post.
.3
Анализ численной схемы
В решение нашей задачи используется сеточный численный метод, основанный
на дискретной аппроксимации уравнений. Основным требование, предъявляемым к
таким сеточным методам, является, прежде всего, обеспечение высокой точности
(малой численной ошибки) получаемых результатов при минимально необходимом
вычислительном ресурсе.
Поэтому очень важно правильно выбрать параметры решения задачи, которые
влияют на точность и длительность расчета данной численной модели. Основным
параметром является шаг сетки. Нашей целью будет выбор такого шага сетки, чтобы
относительная погрешность интересующих нас результатов не превосходила 1-2%.
Рассмотрим
модель вынужденной конвекции. Радиус модели составляет 20 см, толщина слоя -
20см, радиус области стекания жидкости - 2см. На верхней границе втекает
жидкость со скорость 
. В нижней центральной части жидкость вытекает со
скоростью 
. Схематично область моделирования изображена на
рисунке 19.
Были
получены результаты решения задачи для сеточных моделей с разными шагами сетки:
0.5 мм, 1 мм, 2 мм, 3 мм. График изменения средней кинетической энергии со
временем, представленный на рисунке 20, показывает насколько отличаются
результаты решения в зависимости от шага сетки.
Рис. 20. Зависимость средней кинетической энергии от времени
Приняв за точное значение кинетической энергии при шаге сетки 0.5 мм,
относительная погрешность для кинетической энергии при разных шагах сетки будет
соответственно: 1 мм - 1.2 %, 2 мм - 3.2 %, 3 мм - 6.6 %.
Таким образом, результат, полученный для модели с шагом сетки 1 мм,
обладает высокой точностью, меньшим потреблением вычислительных ресурсов по
сравнению с шагом сетки 0.5 мм и является оптимальным выбором для нашей задачи.
4.
Результаты решения задачи свободной конвекции
4.1
Качественное описание формирования конвективного течения
Были проведены расчеты для модели цилиндрического слоя жидкости с
локализованным нагревом в центре. Радиус модели составляет 150 мм, толщина слоя
- 30 мм, радиус нагреваемой области - 50 мм.
В случае подогрева жидкости в центре модели устанавливается режим
течения, схематически изображенный на рис. 21. Локальный нагрев в центральной
части дна создает вертикальный и горизонтальный градиенты температуры.
Горизонтальный градиент температуры приводит к образованию обратной
меридиональной ячейки [5]. Течение в нижней части направлено к центру, а над
центральной частью формируется интенсивное подъемное течение. В верхней части
слоя радиальное движение направлено к периферии. Таким образом, происходит
формирование адвективного течения.
На фоне основного течения, который занимает весь слой, возникают
вторичные движения в виде конвективных струй, которые уносятся основным
потоком (Рис. 25). Они образуются около нагретой горизонтальной поверхности
вследствие неустойчивого распрделения температуры в пограничном слое.
Существенное влияние на формирование вторичных конвективных течений оказывает
распределение температуры - горизонтальный градиент температуры.
Набегающий вдоль дна поток холодной жидкости уносит с собой восходящие
конвективные струи. Конвективные струи формируются на некотором расстоянии от
края нагреваемой области, если тепловой поток поддерживается постоянным, то это
расстояние практически не изменяется. Они возникают периодически по времени,
частота их образования возрастает, а расстояние между точками возникновения
уменьшается с ростом подводимого теплового потока.
Вторичные течения существенно влияют на процессы тепло-массобмена в
пограничном слое, поэтому их исследование представляет большой интерес, как для
фундаментальных, так и для прикладных задач.
4.2
Количественное описание формирования свободной конвекции
.2.1
Результаты CFX
Были рассмотрены конвективные течения с постоянным локализованным
нагревом для воды и масла. Физическое время расчетов 200 секунд, этого было
достаточно для выхода на периодический режим. Начальная температура для обеих
жидкостей 20 оС. Для каждой жидкости расчеты проводились при режимах
с вариацией мощности нагрева P:
5 Вт 10 Вт, 20 ВТ, 30 ВТ, 40 ВТ, 50 ВТ.
Были построены графики зависимости средней кинетической энергии, времени
выхода в квазистационарное движение от мощности нагрева и график зависимости
средней кинетической энергии от времени для всех режимов. Так же построен
график, показывающий изменение периода появления конвективной струи в точке А
(37.5 мм, 3.9 мм, 0 мм) в зависимости от мощности.
Сравнение результатов для режимов воды и масла
На рисунке 27 представлен график зависимости средней кинетической энергии
от времени при разной мощности нагрева. Значение средней кинетической энергии в
слое жидкости определяется следующим соотношением: 
. График показывает, что увеличение
мощности нагрева приводит к заметно более интенсивному конвективному течению
жидкости.
Так же видно, что с некоторого момента времени интенсивность течения
начинает падать. Сравнивая графики средней кинетической энергии для воды и
масла (Рис. 27), можно наблюдать, что течение воды более интенсивно по
сравнению с течением масла. Так же наблюдая за поведением изменения
кинетической энергии масла, можно увидеть, что оно происходит с некоторой
периодичной пульсацией.
Действительно, наблюдая за течением основного потока жидкости, можно
увидеть, что со временем происходит установление периодического режима.
Из графика, представленного на рисунке 28, видно, что время установления
периодического режима у воды больше чем у масла.
Время установления периодического режима t* было подсчитано вручную - были исследованы графики
зависимости полной температуры от времени в точке А (Рис. 29), и найдены
моменты времени, после которого колебания температуры периодически повторялись.
Из этого графика (Рис. 28) видно, что t* уменьшается с увеличением мощности нагрева P.
Изменение средней кинетической энергии в зависимости от P происходит схожим линейным
образом, единственное отличие заключается в том, что величина средней
кинетической энергии в случае воды больше той же величины в случае масла.
Как было замечено ранее, конвективные струи возникают в фиксированных,
регулярно расположенных точках поверхности, если тепловой поток поддерживается
постоянным. Конвективные струи возникают периодически по времени.
Интересно рассмотреть поведение таких мелкомасштабных структур при разной
мощности нагрева. Добавив точку наблюдения А в непосредственной близости к
очагу образования конвективной струи, можно наблюдать изменение полной
температуры в этой точке при разной мощности нагрева.
Рисунок 33 иллюстрирует изменение частоты генерации мелкомасштабных
структур от мощности нагрева в фиксированной точке А для воды и масла. Частота
генерации мелкомасштабных структур в масле больше чем в воде.
4.2.2 Результаты FLUENT
Для режимов конвективного течения масла и воды во FLUENT получаются похожие результаты.
Рисунки 32 и 33 иллюстрируют изменение кинетической энергии со временем
при разной мощности нагрева для масла и воды. График для масла показывает, что
периодическое колебание кинетической энергии происходит на высоких частотах -
признак влияния мелкомасштабных структур, частота возникновения которых высока.
График для воды показывает, что колебания происходят с низкочастотными
пульсациями, связанное с течением основного потока конвективной ячейки.
4.2.3 Сравнение результатов CFX и FLUENT с результатами эксперимента
В эксперименте были получены следующие значения частот генерации
мелкомасштабных структур для масла при разной мощности нагрева.
|
P, Вт
|
5,7
|
10
|
15
|
20,3
|
26,3
|
30
|
35
|
40
|
45
|
|
ν,
1/с
|
0,068
|
0,113
|
0,146
|
0,179
|
0,211
|
0,246
|
0,282
|
0,288
|
0,314
|
Можно сравнить результаты эксперимента с результатами, полученными в CFX и FLUENT.
Из графика, представленного на
рисунке, видно, что значения частоты генерации мелкомасштабных структур,
полученных в CFX и FLUENT, являются схожими. А экспериментальные данные хорошо
согласуются с численными данными при малых мощностях нагрева и существенно
различаются при высоких мощностях. Это объясняется тем, что с ростом нагрева
реальное течение, исследуемое в эксперименте, теряет осесимметричность, а в
двумерных расчетах мы навязываем условие осесимметричности. Таким образом,
нужно переходить к трехмерным расчетам и искать механизм который обуславливает
отличие осесимметричного случая и трехмерного, при этом в эксперименте
необходимо проверить существование низкочастоных колебаний кинетической энергии
основного течения.
4.3
Результаты решения задачи вынужденной конвекции
Были проведены расчеты для модели цилиндрического слоя жидкости со стоком
в центре. Радиус модели составляет 20 см, толщина слоя - 20см, радиус области
стока жидкости - 2см.
В случае стока жидкости в центре модели, устанавливается режим течения,
схематически изображенный на рис. 36.
На рисунках 37 и 38 представлены мгновенные поля скорости течения в
радиальном сечении в разных программных пакетах.
Были проведены расчеты для разных режимов течения с вариацией скорости
стока (м/с): 0,01, 0,02, 0,04, 0,08, 0,16.
Нашей целью в данной постановке является определение наличия вихревых
движений в слое жидкости. Для того чтобы выявить возникновение вихревых
движений в зоне стока, были изучены поля азимутальной компоненты скорости в
радиальном сечении. Результаты для всех режимов течения показали, что
азимутальная компонента скорости во всем сечении равна нулю. Таким образом, в
осесимметричной постановке задачи в неподвижном слое жидкости вихревые движения
не возникают. Это означает, что это явление имеет трехмерную природу и нужно
переходить к трехмерным расчетам для изучения вихревых движений в неподвижном
слое жидкости.
Заключение
В данной работе были исследованы конвективные течения в цилиндрическом
слое жидкости с локализованным нагревом в центральной области в пакете ANSYS. Результаты проведенного
исследования показали, что варьирование мощности нагрева в рамках используемой
численной модели не приводит к существенным изменениям в структуре течений. От
жидкости и мощности нагревателя зависят только количественные характеристики:
изменение интенсивности движения и частоты появления конвективной струи. С
увеличением мощности нагрева увеличивается средняя кинетическая энергия
течения, и частота появления вторичных структур в виде конвективных струй.
В ходе решения задачи было проведено сравнение решений, полученных
различными решателями CFX и FLUENT и экспериментом: значения частоты
генерации мелкомасштабных структур, полученных в CFX и FLUENT,
являются схожими. На масле и CFX и FLUENT хорошо согласуются с экспериментом, но только для малых мощностей
нагрева. Таким образом, нужно переходить к трехмерным расчетам и искать
механизм который обуславливает отличие осесимметричного случая и трехмерного.
Были исследованы конвективные течения в неподвижном цилиндрическом слое
жидкости со стоком в центре. В осесимметричной постановке задачи в неподвижном
слое жидкости вихревые движения не возникают. Это означает, что это явление
имеет трехмерную природу и нужно переходить к трехмерным расчетам для изучения
вихревых движений в неподвижном слое жидкости.
Литература
1. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Том VI. Гидродинамика. - М.: Наука. - 1986.
- C. 306-308.
2. Гершуни Г.З., Жуковицкий Е.М. Конвективная неустойчивость
несжимаемой жидкости. - М.: Наука. - 1986. - C.7-12.
3. Педлоски Дж. Геофизическая гидродинамика. Том I. - М.: Мир. - 1984. - C. 26-34.
. Сухановский А.Н. Формирование дифференциального вращения в
цилиндрическом слое жидкости// Вычислительная механика сплошных сред. - 2010. -
Т. 3, № 2. - С. 103-115
5. Shinji Yukimoto, Hiroshi Niino, Takashi
Noguchi, Ryuji Kimura, Frederic Y. Moulin. Structure of a bathtub vortex:
importance of the bottom boundary layer // Theor. Comput. Fluid Dyn. (2010)
. A. Andersen, T. Bohr, B. S Tenum, J. Juul
Rasmussen. The bathtub vortex in a rotating container // J. Fluid Mech. (2006),
vol. 556, pp. 121-146.
. R. Fernandez-Feria and E. Sanmiguel-Rojas.
On the appearance of swirl in a confined sink flow // PHYSICS OF
FLUIDS. 2000
Приложения
Программный код сессии CFD-Post Session «*.cse»
! for ($i=0; $i <= 200; $i++) {
> load timestep=$i*10:Export Data = Element Heat FluxFile
Format = ANSYSReference Temperature = 0.0 [K]Specify Reference Temperature =
OffSupplemental HTC = 0.0 [W m^-2 K^-1]Profile Type = Inlet
VelocityConnectivity = OffCoord Frame = GlobalFile =
C:/Users/c3po/Desktop/Report/OIL_REPORT/30WOIL/30WOIL_test/test$i.txtGeometry =
OnNode Numbers = OnNull Data = OnType = GenericUnits System = CurrentVariable
Type = CurrentFile Information = OffHeader = OffList = Domain Interface 1 Side
1Token = null= On= 8= " "Variables = X,Y,ZList = Total Temperature,
Velocity u, Velocity v, Velocity wBrackets = ()Display = Scalar
>export
!}
Код программы, выполненный в системе Wolfram Mathematica
n = 200;[i = 0, i <= n, i++,
file[i] =
StringJoin[{"E:\\Report\\WATER_REPORT\\30W\\30W_test\\",
"test", ToString[i], ".txt"}];[i] = Import[file[i],
"Table"]];= Length[table[0]];[i = 0, i <= n, i++, at[i] =
Sum[table[i][[j]][[5]], {j, l}]/l]= Table[{i, at[i]}, {i, 0, 200, 1}];=
ListPlot[temptab, Joined -> True, PlotRange -> Full, Ticks ->
{Table[i, {i, 0, 200, 20}]}, AxesLabel -> {"s", K}, AxesStyle
-> Directive[FontSize -> 22], AspectRatio -> 1/3, PlotStyle ->
Directive[Thick], PlotLabel -> Style["Средняя температура",
28]];["C:\\Users\\c3po\\Desktop\\Report\\OIL_REPORT\\50WOIL\\50WOIL_\temtab.gif",
Lt, ImageSize -> {936,
396.5}];["C:\\Users\\c3po\\Desktop\\Report\\OIL_REPORT\\50WOIL\\50WOIL_\temtab.txt",
temptab, "Table"];[i = 0, i <= n, i++,
ae[i] = 10^7 Sum[\[Pi]*table[i][[j]][[2]]
(table[i][[j]][[6]]^2 + table[i][[j]][[7]]^2 + table[i][[j]][[8]]^2)/l, {j,
l}]];= Table[{i, ae[i]}, {i, 0, 200, 1}];= ListPlot[entab, Joined -> True,
PlotRange -> Full, Ticks -> {Table[i, {i, 0, 200, 20}]}, AxesLabel ->
{"s","10^7m^2/s^2"}, AxesStyle -> Directive[FontSize
-> 22], AspectRatio -> 1/3, PlotStyle -> Directive[Thick], PlotLabel ->
Style["Средняя кинетическая энергия",
28]];["E:\\Report\\WATER_REPORT\\30W\30W_entab.gif", Le,
ImageSize -> {936,
396.5}];["E:\\Report\\WATER_REPORT\\30W\\30W_entab.txt", entab,
"Table"];
Программный код Gambit
vertex: vertex.2> vertex create coordinates 0 0.15
0vertex: vertex.3> vertex create coordinates 0.03 0.15 0vertex: vertex.4>
vertex create coordinates 0.03 0 0vertex: vertex.5> edge create straight
"vertex.1" "vertex.2"edge: edge.1> edge create straight
"vertex.2" "vertex.3"edge: edge.2> edge create straight
"vertex.3" "vertex.4"edge: edge.3> edge create straight
"vertex.4" "vertex.5"edge: edge.4> edge create straight
"vertex.1" "vertex.5"edge: edge.5> face create wireframe
"edge.1" "edge.2" "edge.3" "edge.4"
"edge.5" realface: face.1> window modify shade> face mesh
"face.1" map size 0.001generated for face face.1: mesh faces =
4500.> physics create "heat" btype "WALL" edge
"edge.1"Boundary entity: heat> physics create "bottom"
btype "WALL" edge "edge.2"Boundary entity: bottom>
physics create "side" btype "WALL" edge
"edge.3"Boundary entity: side> physics create "top"
btype "WALL" edge "edge.4"Boundary entity: top> physics
create "axis" btype "AXIS" edge "edge.5"Boundary
entity: axis> physics create "fluid" ctype "FLUID" face
"face.1"Continuum entity: fluid