Abstract

Abstract

The governing equations of fluid and gas mechanics are considered such as continuity equation and Navier – Stokes equations. The application of Reynolds decomposition for closure problem of turbulence is described. Boussinesq hypothesis which allows introducing the conception of turbulent viscosity is considered. The most widely-known k–ε models are described as well as their advantages and disadvantages.

Full Text

В наши дни пакеты численного моделирования открывают огромные возможности для инженеров и исследователей из самых разных областей. Путем простого нажатия ряда кнопок и ввода необходимых входных данных исследователь может получать решения для комплексных междисциплинарных задач. Между тем все эти пакеты основаны на фундаментальных законах механики, и каждый входной параметр является критически важным для получения достоверного и соответствующего действительности решения. В нашей статье рассматриваются основные понятия механики жидкости и газа. Рассматривается уравнение Навье – Стокса и дается возможное физическое объяснение основных его компонент. Также рассматриваются основные модели, принятые для описания турбулентности и широко применяемые в пакетах численного моделирования. Ставится задача рассмотреть теоретический базис, необходимый для дальнейшего численного решения проблем из области турбулентности. Основные понятия механики жидкости. Согласно макроскопической модели вещества [1] жидкость и газ представляют собой сплошную текучую изотропную ньютоновскую среду с непрерывным распределением массы и других физических величин. Вспомним несколько основных понятий, применяемых в механике жидкости и газа (далее будем говорить о механике жидкости, подразумевая, что модели, описывающие поведение жидкости, также пригодны для описания поведения газа). Текучесть среды – свойство неограниченной деформируемости среды, т.е. способность изменять свою форму под действием сколь угодно малых сил, если жидкость не сдерживается какими-либо стенками. Сплошность или неразрывность среды – способность заполнять весь объем, занимаемый материалом тела, без всяких пустот, общность свойств любой части среды и среды в целом. Изотропность среды – независимость всех физических величин и свойств среды от направления. Ньютоновская среда – среда, в которой касательные напряжения прямо пропорциональны градиенту скорости (или скорости угловых деформаций). Кроме поля скоростей также рассматриваются скалярные величины: плотность ; температура ; тензоры напряжений и скоростей деформаций. Каждая из данных величин является функцией координат и времени. Одним из основных уравнений в механике жидкости и газа является уравнение Навье – Стокса . (1) Данное уравнение подробно описывает изменение скорости жидкости по времени с помощью четырех компонент. Первая из них, , показывает, как дивергенция влияет на скорость. Ее физический смысл можно наглядно объяснить на примере течения реки [2]. Так, когда река сужается, скорость воды в ней возрастает, и наоборот, когда река расширяется, скорость воды уменьшается. Вторая компонента, , показывает, как влияет на движение изменение давления, особенно на направленность движения от областей с более высоким давлением. Также ясно, что чем больше плотность жидкости, тем труднее ей осуществлять перемещение. Следующая компонента , где – кинематическая вязкость, показывает влияние, оказываемое на частицу со стороны соседних частиц. Чем больше вязкость, тем, соответственно, больше величина данного влияния. И четвертая компонента, , характеризует влияние, оказываемое на данную жидкость со стороны любой другой силы. Другим фундаментальным уравнением является уравнение неразрывности (2) Для несжимаемой жидкости и уравнение приобретает вид (3) Уравнения (2), (3) справедливы как для идеальной, так и для реальной жидкости [3]. Таким образом, для ламинарной жидкости мы получаем систему из четырех уравнений: три уравнения Навье – Стокса в проекциях на оси и уравнение неразрывности для четырех неизвестных: три компонента вектора скорости и гидродинамическое давление. Декомпозиция Рейнольдса. Турбулентная жидкость характеризуется колебаниями скорости во всех направлениях и имеет бесконечное число степеней свободы [4]. Решение уравнений Навье – Стокса для турбулентной жидкости затруднено, так как в данном случае уравнения эллиптические, нелинейные и содержат по две неизвестных величины. Жидкость в данном случае хаотическая, диффузионная, диссипативная и прерывистая. Существует несколько путей решения данной проблемы. Одним из них является декомпозиция Рейнольдса, согласно которой произвольную величину можно записать как сумму ее среднего значения и отклонения [5]: . Такая декомпозиция будет давать набор уравнений, описывающих некоторое среднее поле жидкости. В результате мы получим усредненные по Рейнольдсу уравнения Навье – Стокса, которые также называются уравнениями Рейнольдса, а также усредненное уравнение неразрывности. Уравнение неразрывности в компонентах для несжимаемой жидкости имеет вид (4) Тогда для усредненной скорости (5) Вычитая (5) из (4), получаем уравнение неразрывности для отклонения Используя (4), можно записать уравнение (1) в компонентах следующим образом: , (6) где – напряжения в жидкости, определяются по формуле (7) Соотношения (7) являются определяющими соотношениями для ньютоновской жидкости. называется дельтой Кронекера и определяется как Используя декомпозицию Рейнольдса, уравнение (6) можно записать в следующем виде [5]: (8) Это уравнение известно как уравнение Рейнольдса. Данное уравнение достаточно похоже на уравнение (6) и отличается лишь дополнительным слагаемым в правой части . Это слагаемое называется напряжениями Рейнольдса и представляет собой симметричный тензор второго порядка, состоящий из шести независимых компонент. Таким образом, для турбулентной жидкости имеются все те же четыре уравнения и уже десять неизвестных: три компонента скорости, гидродинамическое давление и шесть напряжений Рейнольдса. Данная проблема носит название проблемы замыкания турбулентности. Стандартная модель турбулентности . Чтобы замкнуть турбулентность, необходимо определить связь между напряжениями по Рейнольдсу и параметрами осредненного течения. Эту связь определяют с помощью различных моделей турбулентности [6]. В этих моделях принимаются определенные допущения, на основе которых вводится недостающее число уравнений, что позволяет найти все неизвестные. Одним из допущений является введение турбулентной вязкости, которое впервые осуществил Буссинеск. Турбулентную динамическую вязкость он ввел по аналогии с динамической вязкостью (9) Далее перейдем непосредственно к получению стандартной модели из двух уравнений, которая сегодня рассматривается как стандартная модель для описания турбулентности и решения инженерных задач. В данной модели вводятся два важных понятия – генерация Р и диссипация ε. Физический смысл генерации турбулентности Р заключается в порождении новых вихрей и пульсаций, которые и образуют турбулентность [7]. Диссипация ε, напротив, представляет собой рассеивание больших вихрей на более малые, приводит к усреднению течения и уменьшению турбулентности. Два уравнения переноса позволяют рассматривать турбулентность в пространстве и времени. Данная модель является полуэмпирической и опирается на феноменологический поход и результаты, полученные опытным путем. Выполнив некоторые алгебраические преобразования и умножив на (8) можно привести к следующему виду [8]: (10) Определим кинетическую энергию турбулентности как и подставим ее в (10), принимая : (11) Второе слагаемое правой части (11), по определению [8], является диссипацией (12) тогда как четвертое слагаемое правой части выражения (1.1), включая минус, по определению [8], является генерацией Р: . (13) Далее делаем допущение [8], что . (14) Учитывая (12)–(14), уравнение (11) можно записать в следующем виде: . (15) Уравнение (15) является уравнением для кинетической энергии – параметр, обеспечивающий нужную размерность для слагаемого с . Обычно принимается Уравнение для диссипации аналитически не выводится и просто записывается по аналогии с (15): , (16) где обеспечивает нужную размерность, а константы вводятся, поскольку форма уравнения (16) лишь предполагается, но не выводится аналитически. В пакете Ansys Fluent уравнения cтандартной модели применяются в несколько ином, модернизированном виде. Его можно получить из (15) и (16) путем алгебраических преобразований, и он также описывается создателями модели [9], www.ansys.com: (17) В данной системе уравнений представляет турбулентную кинетическую энергию, образованную от средних градиентов скорости. Принимая гипотезу Буссинеска, ее можно выразить по формуле , где , (18) – плотность газа, – инвариант тензора деформаций, ; – кинетическая энергия выталкивающей силы, , где – турбулентная постоянная Прандтля для энергии, – компонента вектора гравитации в i-м направлении; – коэффициент температурного расширения, , где Т – температура. – константа, определяющая степень воздействия выталкивающей силы на определяется по формуле где – компонента скорости жидкости, параллельная скорости гравитации и – компонента скорости жидкости, перпендикулярная скорости гравитации. для слоев жидкости, для которых направление скорости жидкости параллельно вектору гравитации, для слоев жидкости, для которых направление скорости жидкости перпендикулярно вектору гравитации. – вклад переменного расширения при турбулентности сжатия в общую скорость диссипации. Данную величину следует учитывать при большом числе Маха. Ее обязательно учитывать, когда моделируется сжимаемый идеальный газ. где – число Маха для турбулентной жидкости, , где а – скорость звука, Остальные константы определены из экспериментов для фундаментальных турбулентных жидкостей и имеют следующие значения: RNG модель. RNG модель была получена при помощи теории ренормализованных групп [10]. Она имеет схожую форму со стандартной моделью, но включает следующие улучшения: – имеет дополнительный член в уравнении для который улучшает точность вычислений для жидкостей с высокими скоростями деформаций; – в модели учтено влияние завихренности на турбулентность, что увеличивает точность для высокозавихренных жидкостей; – данная теория предлагает аналитические формулы для турбулентных чисел Прандтля, тогда как стандартная модель использует заданные пользователем постоянные значения; – RNG модель предлагает аналитически полученные формулы для эффективной вязкости, которая предназначена для жидкостей с низкими числами Рейнольдса. Тем не менее эффективное использование этой опции зависит от правильного рассмотрения пристеночной области. Данные улучшения делают RNG модель более точной и надежной, позволяя эффективно применять ее для более широкого класса жидкостей по сравнению со стандартной моделью. Уравнения RNG модели имеют следующий вид: (19) Далее рассмотрим значение величин, новых по сравнению с (17). – обратные эффективные числа Прандтля для и соответственно. В (19) означает эффективную вязкость. Данная вязкость приблизительно равна из стандартной модели для высоких чисел Рейнольдса. Для низких чисел Рейнольдса создателями модели [10] предлагается дополнительное дифференциальное уравнение, позволяющее более точно вычислить Главное отличие RNG модели от стандартной заключается в дополнительном члене в уравнении для вычисляется по формуле , (20) где Значение может стать более очевидным, если записать второе уравнение (20) в следующем виде: , где . В случае когда вносит положительный вклад, становится больше, чем Таким образом, для жидкостей со слабыми или умеренными скоростями деформаций RNG модель дает результаты, схожие с полученными при помощи стандартной модели. При больших скоростях деформаций, вносит отрицательный вклад, становится меньше, чем , снижается и следовательно, эффективная вязкость. В результате, в жидкостях с большими скоростями деформаций RNG модель дает меньшую турбулентную вязкость, чем стандартная модель. Константы и имеют следующие значения: . Реальная модель. По сравнению со стандартной моделью данная модель имеет два существенно важных отличия [11]: – реальная модель содержит альтернативную формулировку для турбулентной вязкости; – модифицированное уравнение переноса для скорости диссипации было получено из точного уравнения для переноса среднеквадратичных колебаний завихренности. Данная модель удовлетворяет точным математическим ограничениям по напряжениям Рейнольдса, вытекающим из физики турбулентной жидкости. Неравенство Шварца и положительный знак напряжений по Рейнольдсу накладывают некоторые ограничения на Значения данной константы меняются в зависимости от свойств жидкости и местоположения, они достаточно точно были получены экспериментальным путем для разных условий. Другая проблема заключается в том, что уравнение для скорости диссипации не всегда работает достаточно хорошо. Например, хорошо известна аномалия круглого поперечного сечения струи. Скорость распространения для круглого сечения описывается достаточно хорошо, тогда как для ассиметричного результаты получаются неудовлетворительными. Реальная модель учитывает эти недостатки с помощью следующих улучшений: во-первых, предлагается новая формула для определения турбулентной вязкости, первоначально предложенная еще Рейнольдсом; во-вторых, используется новое уравнение для диссипации основанное на динамическом уравнении среднеквадратичных колебаний завихренности. Ограничением является то, что можно получить нефизичные турбулентные вязкости в ситуациях, когда вычислительная область содержит как зоны с турбулентностью, так и зоны со стационарной жидкостью. Уравнения реальной модели имеют следующий вид: При этом В реальной модели уравнение для такое же, как и в стандартной модели. В то же время уравнение для отличается существенно. Одной положительной чертой является то, что правая сторона уравнения для не содержит . Считается, что это обеспечивает лучший перенос спектральной энергии. Другое преимущество заключается в том, что при не возникает деления на ноль. Хорошие результаты для этой модели были получены для вихревых однородных жидкостей со сдвиговыми напряжениями, для свободных жидкостей и разделенных жидкостей. Для всех этих случаев были получены существенно лучшие результаты, чем при использовании стандартной модели. Также были получены хорошие результаты для ассиметричного поперечного сечения струи. В реальной модели нахождение турбулентной вязкости существенно отличается от других моделей. Хотя сама вязкость также находится по формуле (18), уже не является константой и находится по формуле , где где – тензор средних скоростей вращения; – угловая скорость. Константы где . В данном случае является функцией средней скорости вращения и средней скорости деформаций, угловой скорости и переменных турбулентности и . Остальные константы имеют следующие значения: Выводы и рекомендации. В теории турбулентности еще не разработана универсальная теория, позволяющая одинаково успешно находить решения для всех классов задач. модели являются наиболее широко распространенными и применимыми в численных пакетах. Результаты, полученные с их помощью, получили практическое подтверждение для широкого класса функций. Между тем не следует забывать, что данные модели имеют как преимущества, так и недостатки, и существует огромное количество альтернативных моделей. Выбор модели должен осуществляться в зависимости от условий конкретной задачи, а результаты решения должны быть тщательно проверены.