🗊Презентация Обзор методов расчета турбулентных течений

Обзор методов расчета турбулентных течений Лекция 7

Характеристики турбулентности 1. Турбулентность можно рассматривать как совокупность вихрей различного размера, от колмогоровского масштаба (самые мелкие вихри) до когерентных структур (самые крупные вихри с масштабом, соизмеримым с геометрией, например, толщина пограничного слоя или ширина канала). 2. Каждый вихрь характеризуется размером и временем жизни. 3. Процесс передачи энергии от крупных вихрей к мелким играет фундаментальную роль и называется каскадный перенос энергии. 4. Все переменные в турбулентном потоке можно разделить на осредненную и пульсационные составляющие. 5. Турбулентные пульсации как случайные величины рассматриваются в рамках статистической теории турбулентности.

Турбулентные вихри Турбулентность – это совокупность турбулентных вихрей различного размера. 1. Начиная с некоторого малого размера вихри диссипируют в тепло под воздействием вязкости. Характерный минимальный размер вихрей определяется колмогоровским масштабом – местная скорость диссипации на единицу массы – кинематическая вязкость, м2∙с. – кинетическая энергия турбулентного течения, м2∙с-2 – скорость диссипации турбулентной кинетической энергии (имеет размерность энергии, отнесенной к единице массы и единице времени), м2∙с-3

Турбулентные вихри 2. Максимальный размер турбулентных вихрей Turbulence Length Scale (когерентных структур), существующих в потоке, имеет порядок характерного линейного масштаба течения : – для каналов круглого сечения – для пограничного слоя

Каскадный перенос энергии 1. Энергия поступает от осредненного потока к наиболее крупным вихрям. 2. Последовательно передается всё более мелким вихрям (каскадный перенос) 3. Самые мелкие (колмогоровские) вихри диссипируют и передают энергию тепловому движению.

Области энергетического спектра

Физическая картина турбулентности образно выражена в следующем четверостишии, написанном английским физиком Л. Ричардсоном в 1922 г.: Физическая картина турбулентности образно выражена в следующем четверостишии, написанном английским физиком Л. Ричардсоном в 1922 г.: Big whirls make little whirls Which feed on their velocity, Little whirls have smaller ones And so on into viscosity. В поток бурлящий бросив взгляд, ВихрЕй увидишь там каскад: МеньшОй у большего энергию берет, Пока мельчайших вязкость не сотрет.

Статистическая теория турбулентности В рамках статистической теории турбулентности мгновенное поле течения делится на осредненное и пульсационное Таким образом, турбулентное движение состоит как бы из регулярного течения, описываемого осреднёнными значениями скоростей, и из наложенного на него хаотического пульсационного течения. Определение осредненного движения зависит от выбранного способа осреднения. Чаще всего производят осреднение по времени и пространству (Условия осреднения Рейнольдса).

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

Обзор методов расчета турбулентных течений 1. Прямое численное моделирование (метод DNS) Нестационарные уравнения динамики вязкой жидкости (уравнения Навье-Стокса) описывают движение в турбулентном течении вплоть до минимальных масштабов турбулентности. Однако при численном решении этих уравнений для того, чтобы учесть все масштабы, может потребоваться настолько мелкая сетка, что даже современные компьютерные мощности не позволят решить такую задачу. То же относится и к выбору шага численного интегрирования по времени, так как характерное время мелкомасштабной турбулентности очень мало. С другой стороны, именно мелкомасштабная турбулентность играет важнейшую роль при описании турбулентных течений. Поэтому прямое численное моделирование (Direct Numeric Simulation, DNS) турбулентных течений применяется для инженерных расчетов достаточно редко.

Прямое численное моделирование (метод DNS)

Пример расчета с использованием метода DNS

Обзор методов расчета турбулентных течений 2. Моделирование крупных вихрей (LES) Более простой моделью является моделирование крупных вихрей (Large Eddy Simulation, LES). В этом подходе крупные вихри рассчитываются явно, а мельчайшие вихри подсеточного масштаба (Sub-Grid Scale, SGS) моделируются с использованием правил подсеточного замыкания. Основной предпосылкой такого подхода является то, что наибольшие вихри, которые находятся под прямым воздействием граничных условий, несут максимум энергии и должны быть рассчитаны. Этот подход имеют хорошую перспективу.

Метод моделирования крупных вихрей Large Eddy Simulation (LES)

Процедура фильтрации в модели LES

Практическая реализация метода LES

Свойства метода LES

Обзор методов расчета турбулентных течений 3. Осредненные по Рейнольдсу уравнения Навье-Стокса (RANS) В настоящее время наиболее распространенным способом моделирования турбулентности является использование осреднения Рейнольдса, когда вместо уравнений для мгновенных значений параметров используются уравнения для неких осредненных величин. Эти уравнения называются уравнениями Рейнольдса.

Осредненные по Рейнольдсу уравнения Навье-Стокса (RANS)

Гипотеза Буссинеска

Полуэмпирические модели турбулентности

Достоинства метода RANS

Недостатки метода RANS

Обзор методов расчета турбулентных течений 4. Гибридные методы расчета (RANS-LES) При приближении к стенке размер характерных вихревых структур уменьшается пропорционально расстоянию до стенки. Пропорционально должен уменьшаться шаг сетки. Для точного воспроизведения профиля скорости в пограничном слое при существенных числах Рейнольдса это ведет к гигантским вычислительным затратам. Единственный выход - использовать LES только в части области. В остальной части расчетной области использовать более экономичные методы RANS (URANS).

Гибридные методы расчета

Проблемы гибридных методов расчета Серая область • Какие уравнения решаются на границе RANS и LES областей? В RANS подобласти используются уравнения Рейнольдса полученные осреднением по времени. В LES области используются отфильтрованные уравнения, полученные осреднением по пространству (фильтрация). • Хорошо ли разрешаются крупные вихри около границы областей RANS-LES? В RANS области все турбулентные пульсации моделируются эмпирическими зависимостями. В LES области большая часть пульсаций разрешается точно. В окрестности границы происходит взаимное влияние областей. В LES области разрешенные турбулентные пульсации подавляются за счет соседства со стационарной RANS областью. Конвективный перенос приводит к «сносу» RANS решения в LES область и наоборот. В окрестности границы областей возникает проблемная область, называемая «серая область» (grey area).

Проблемы гибридных методов расчета Обеспечение «хорошего» решения в серой области - это основная проблема гибридных методов. Из-за разнообразия течений не существует универсального решения. Проблема «серой области» может быть решена разными способами, например: 1. Поместить границу в малочувствительную область (DES) – окрестность внешней границы пограничного слоя, в которой градиенты скорости сравнительно малы. 2. Использование естественной неустойчивости низкодиссипативных схем на больших градиентах (WMLES, IDDES) – «серая область» расположена внутри пограничного слоя. 3. Введение искусственных пульсаций на границе RANS-LES областей – необходимо создать такие же пульсации, как в «истинной» турбулентности.

Сравнение различных методов расчета турбулентных течений

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

Выводы

Полуэмпирические модели турбулентности ANSYS CFX

Проблемы моделирования

Проблемы моделирования

Истинная и модельная профили скорости

k-ε модель Модель турбулентности из 2 уравнений очень широко используются, т.к. является компромиссом между численными затратами и точностью вычислений. Скорость турбулентных пульсаций и масштаб длины вихря находятся с использованием отдельных уравнений переноса (отсюда термин ‘модель двух уравнений'). k - кинетическая энергия турбулентности, определяется как дисперсия изменений скорости, м2/с2. ε – скорость изменения кинетической энергии, м2/с3. В модель введены пристенные функции, т.е. для этих моделей в уравнения вводятся дополнительные функции, отвечающие за влияние стенок на турбулентность.

k-ε модель

k-ε модель К достоинствам k-ε -моделей относится высокая точность при расчете свободных сдвиговых течений. Они достаточно универсальны и не требуют задания каких-либо дополнительных параметров. Тем не менее, трудности, связанные с их применением в пристенных областях, заставляют исследователей изобретать все новые и новые модели. Однако до настоящего времени никаких предпочтений среди моделей, базирующихся на концепции осреднения по Рейнольдсу, по существу, не наблюдается, поскольку не существует «универсальной» модели турбулентности. Особенности: Буферный слой не моделируется, для расчета скорости у стенки используются пристеночные функции. Благодаря быстрой сходимости и относительно низким требованиям к объему памяти k-ε модель очень популярна при решении промышленных задач. Она не очень точна при моделировании течений с положительным градиентом давления, струйных течений и течений в области с сильно искривленной геометрией. Модель хорошо подходит для решения задач внешнего обтекания тел сложной геометрической формы. Например, k-ε модель можно использовать для моделирования потока вблизи плохо обтекаемого тела

k-ω модель Модель турбулентности - ‘модель двух уравнений‘. Модель позволила существенно улучшить описание пристенных течений в рамках рассматриваемого класса моделей турбулентности благодаря использованию вместо уравнения для скорости диссипации ε уравнения для параметра ω. к - кинетическая энергия турбулентности, определяется как дисперсия изменений скорости, м2/с2. ω – удельная скорость диссипации турбулентной энергии, 1/с3. Является более точной и надежной, в отличие от пред. модели, менее требовательна к качеству сетки при низкорейнольдсовой турбулентности.

k-ω модель

k-ω модель k-ω модель похожа на k-ε, только здесь решается уравнение для удельной скорости диссипации кинетической энергии ω. В данной модели также используются пристеночные функции, поэтому требования к ресурсам памяти здесь те же, что и при использовании k-ε модели. Сходимость при использовании данной модели чуть медленнее и существенно зависит от начального приближения. Зачастую решение, полученное с помощью k-ε модели, используется в качестве хорошего начального приближения для расчета по k-ω модели. Использование k-ω модели дает хорошие результаты в тех задачах, где k-ε модель недостаточно точна, например, при моделировании внутренних течений, течений по сильно искривленным каналам, отрывных и струйных течений. Хорошим примером применения k-ω модели является задача о течении жидкости через колено трубопровода.

Shear Stress Transport (SST) Модель сдвиговых напряжений, модель Ментера Включает в себя об два типа моделей ( k – ε вдали от стенок, k –ω вблизи к стенкам), а также рассчитывает перенос напряжений сдвига турбулентности для определения турбулентной вязкости.

Shear Stress Transport (SST)

Shear Stress Transport (SST) SST-модель представляет собой комбинацию k-ε, и k-ω моделей турбулентности: для расчета течения в свободном потоке используются уравнения k-ε модели, а в области вблизи стенок — уравнения k-ω модели. Пристеночные функции не используются. Данная модель, пожалуй, дает наиболее точные результаты при расчете течений вблизи твердых стенок. Уравнения SST-модели не всегда быстро сходятся, поэтому часто для получения хорошего начального приближения используются решения, полученные с помощью k-ε- или k-ω модели. В тестовом примере задачи об обтекании крылового профиля используется SST-модель турбулентности, полученные результаты хорошо согласуются с экспериментальными данными.

Reynolds Stress Model (RSM) модель Модель Рейнольдсовых напряжений имеет англоязычную аббревиатуру RSM (Reynolds Stress Model) и является одной из самых сложных моделей турбулентности предлагаемых CFX. Эта модель не использует предположение о изотропности турбулентной вязкости, а для замыкания уравнений Навье-Стокса, осредненных по Рейнольдсу, решает уравнения переноса для Рейнольдсовых напряжений совместно с уравнением для скорости турбулентной диссипации ε. Так как модель RSM описывает эффекты кривизны, закрученности, вращения, резкого изменения напряжений между слоями более строго, чем одно- и двух- параметрические модели турбулентности, то она имеет больший потенциал для более точного расчета сложных потоков. Однако RSM модель все-таки имеет некоторые упрощения, которые были приняты для составления уравнений переноса Рейнольдсовых напряжений, что было необходимо для замыкания системы уравнений Навье-Стокса. Использование этой модели турбулентности рекомендуется в случаях, когда анизотропность турбулентного потока оказывает доминируещее влияние на характер турбулентного течения (циклоны, сильно закрученные потоки в камерах сгорания, вращающиеся области, вторичные течения в каналах, вызванные большими нормальными напряжениями).

Основные проблемы, возникающие при использовании моделей рейнольдсовых напряжений 1. Сложность системы дифференциальных уравнений (не менее 7 уравнений для турбулентности), требующих для решения больших вычислительных ресурсов. 2. Эта система уравнений кроме громоздкости отличается сильной нелинейностью, что приводит к необходимости использовать различные ухищрения для улучшения устойчивости и сходимости. К таковым относятся специальные схемы, векторные прогонки и использование сложных алгоритмов демпфирования. При решении столь нелинейной системы уравнений возникают численные проблемы, для решения которых приходится пользоваться специальными приемами. 3. Полученное решение может быть неприемлемым с физической точки зрения. Необходимо проверять и выбраковывать результаты «дефектных» расчетов. Очевидно, что такой инструмент неприменим для «инженерного» использования. 4. Возникают проблемы с постановкой граничных условий на свободных границах для напряжений. 5. При расчете задач теплопередачи необходимо решать еще 3 дополнительных уравнения для корреляций скорости и температуры, что еще больше усложняет систему.

Reynolds Stress Model (RSM) модель Reynolds Stress Model - LRR-IP SSG Reynolds Stress Model - SSG QI Reynolds Stress Model - LRR-IQ

Reynolds Stress Model (RSM) модель

Reynolds Stress Model (RSM) модель