Вычислительная гидродинамика - Computational fluid dynamics

Вычислительная гидродинамика ( CFD ) - это раздел механики жидкости, который использует численный анализ и структуры данных для анализа и решения проблем, связанных с потоками жидкости . Компьютеры используются для выполнения расчетов, необходимых для моделирования набегающего потока жидкости и взаимодействия жидкости ( жидкостей и газов ) с поверхностями, определяемыми граничными условиями . С помощью высокоскоростных суперкомпьютеров могут быть достигнуты лучшие решения, которые часто требуются для решения самых крупных и сложных проблем. Постоянные исследования позволяют получить программное обеспечение, которое повышает точность и скорость сложных сценариев моделирования, таких как трансзвуковые или турбулентные потоки. Первоначальная проверка такого программного обеспечения обычно выполняется с использованием экспериментального оборудования, такого как аэродинамические трубы . Кроме того, для сравнения можно использовать ранее выполненный аналитический или эмпирический анализ конкретной проблемы. Окончательная проверка часто выполняется с использованием полномасштабных испытаний, таких как летные испытания .

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

Предпосылки и история

Компьютерное моделирование высокоскоростного воздушного потока вокруг космического корабля "Шаттл" при входе в атмосферу .
Моделирование ГПВРД Hyper-X на скорости -7 Маха

Фундаментальной основой почти всех задач CFD являются уравнения Навье – Стокса , которые определяют множество однофазных (газовых или жидких, но не обоих) потоков жидкости. Эти уравнения можно упростить, удалив члены, описывающие вязкие действия, чтобы получить уравнения Эйлера . Дальнейшее упрощение путем удаления членов, описывающих завихренность, дает полные уравнения потенциала . Наконец, для малых возмущений в дозвуковых и сверхзвуковых потоках (не трансзвуковых или гиперзвуковых ) эти уравнения могут быть линеаризованы для получения линеаризованных потенциальных уравнений.

Исторически впервые были разработаны методы для решения линеаризованных потенциальных уравнений. Двумерные (2D) методы, использующие конформные преобразования обтекания цилиндра в обтекание профиля, были разработаны в 30-е годы прошлого века.

Одним из самых ранних типов вычислений, напоминающих современные CFD, являются вычисления Льюиса Фрая Ричардсона в том смысле, что в этих вычислениях использовались конечные разности и физическое пространство делилось на ячейки. Несмотря на то, что они потерпели неудачу, эти расчеты вместе с книгой Ричардсона « Прогноз погоды с помощью числового процесса» заложили основу для современной CFD и численной метеорологии. Фактически, в ранних расчетах CFD в 1940-х годах с использованием ENIAC использовались методы, близкие к методам из книги Ричардсона 1922 года.

Доступные компьютерные возможности ускорили развитие трехмерных методов. Вероятно, первая работа с использованием компьютеров для моделирования потока жидкости в соответствии с уравнениями Навье – Стокса была выполнена в Лос-Аламосской национальной лаборатории в группе T3. Эту группу возглавил Фрэнсис Х. Харлоу , которого многие считают одним из пионеров CFD. С 1957 и до конца 1960 - х лет, эта группа разработала множество численных методов для имитации переходных двумерный потоков жидкости, таких как частицы в ячейке метод, жидкости в ячейке метод, завихренность функции тока метод, и маркер-и-клетка метод . Метод завихренности-функции потока Фромма для двумерного переходного несжимаемого потока был первым в мире исследованием сильно искривленных течений несжимаемой жидкости.

Первая статья с трехмерной моделью была опубликована Джоном Хессом и АМО Смитом из Douglas Aircraft в 1967 году. Этот метод дискретизировал поверхность геометрии с помощью панелей, что привело к появлению этого класса программ, получивших название панельных методов. Сам их метод был упрощен, поскольку он не включал подъемные потоки и, следовательно, в основном применялся к корпусам кораблей и фюзеляжам самолетов. Первый код панели подъема (A230) был описан в статье, написанной Полом Руббертом и Гэри Саарисом из Boeing Aircraft в 1968 году. Со временем более совершенные трехмерные коды панели были разработаны в Boeing (PANAIR, A502), Lockheed (Quadpan) , Дуглас (HESS), McDonnell Aircraft (MACAERO), NASA (PMARC) и аналитические методы (WBAERO, USAERO и VSAERO). Некоторые (PANAIR, HESS и MACAERO) были кодами более высокого порядка, использующими распределения сингулярностей более высокого порядка, в то время как другие (Quadpan, PMARC, USAERO и VSAERO) использовали отдельные особенности на каждой панели поверхности. Преимущество кодов более низкого порядка состояло в том, что они работали намного быстрее на компьютерах того времени. Сегодня VSAERO превратилась в многопользовательский код и является наиболее широко используемой программой этого класса. Он использовался при разработке многих подводных лодок , надводных кораблей , автомобилей , вертолетов , самолетов и, в последнее время, ветряных турбин . Его родственный код, USAERO, представляет собой метод неустойчивой панели, который также использовался для моделирования таких вещей, как высокоскоростные поезда и гоночные яхты . Код NASA PMARC из ранней версии VSAERO и производной от PMARC, названной CMARC, также коммерчески доступен.

В двумерной области был разработан ряд кодов панелей для анализа и проектирования аэродинамического профиля. В коды обычно включен анализ пограничного слоя , чтобы можно было моделировать эффекты вязкости. Ричард Эпплер  [ де ] разработал код ПРОФИЛЯ, частично при финансировании НАСА, который стал доступен в начале 1980-х годов. Это было вскоре после чего Марк дрела «s XFOIL код. И PROFILE, и XFOIL включают двумерные коды панелей со связанными кодами пограничного слоя для работы по анализу профиля. PROFILE использует метод конформного преобразования для инверсного дизайна крылового профиля, в то время как XFOIL имеет как конформное преобразование, так и метод обратной панели для конструирования крылового профиля.

Промежуточным шагом между панельными кодами и кодами полного потенциала были коды, в которых использовались уравнения малых трансзвуковых возмущений. В частности, широкое распространение получил трехмерный код WIBCO, разработанный Чарли Боппе из Grumman Aircraft в начале 1980-х годов.

Разработчики обратились к кодам полного потенциала, поскольку методы панели не могли рассчитать нелинейный поток, присутствующий на околозвуковых скоростях. Первое описание способов использования уравнений полного потенциала было опубликовано Эрлом Мурманом и Джулианом Коулом из Boeing в 1970 году. Фрэнсис Бауэр, Пол Гарабедян и Дэвид Корн из Института Куранта при Нью-Йоркском университете ( Нью-Йоркский университет ) написали серию из двух статей. Широко использовавшиеся пространственные коды полного потенциала аэродинамического профиля, наиболее важная из которых была названа программой H. Дальнейшее развитие программы H было разработано Бобом Мельником и его группой в Grumman Aerospace под названием Grumfoil. Энтони Джеймсон , первоначально работавший в Grumman Aircraft и Курантском институте Нью-Йоркского университета, работал с Дэвидом Коги над разработкой важного трехмерного кода полного потенциала FLO22 в 1975 году. После этого появилось много кодов полного потенциала, кульминацией которых стал код Boeing Tranair (A633), который по-прежнему активно используется.

Следующим шагом были уравнения Эйлера, которые обещали дать более точные решения трансзвуковых потоков. Методология, использованная Джеймсоном в его трехмерном коде FLO57 (1981), использовалась другими для создания таких программ, как программа Lockheed TEAM и программа MGAERO IAI / Analytical Methods. MGAERO уникален тем, что является кодом структурированной декартовой сетки, в то время как большинство других таких кодов используют структурированные сетки, приспособленные к телу (за исключением очень успешного кода CART3D НАСА, кода SPLITFLOW от Lockheed и NASCART-GT Технологического института Джорджии). Энтони Джеймсон также разработал трехмерный код AIRPLANE, в котором использовались неструктурированные тетраэдрические сетки.

Что касается двумерной области, то Марк Дрела и Майкл Джайлз, тогда аспиранты Массачусетского технологического института, разработали программу ISES Euler (фактически набор программ) для проектирования и анализа аэродинамического профиля. Этот код впервые стал доступен в 1986 году и получил дальнейшее развитие для проектирования, анализа и оптимизации одно- или многоэлементных аэродинамических поверхностей в виде программы MSES. MSES находит широкое применение во всем мире. Производным от MSES для проектирования и анализа аэродинамических профилей в каскаде является MISES, разработанный Гарольдом Янгреном, когда он был аспирантом Массачусетского технологического института.

Уравнения Навье – Стокса были конечной целью разработки. Двумерные коды, такие как код ARC2D НАСА Эймса, впервые появились. Был разработан ряд трехмерных кодов (ARC3D, OVERFLOW , CFL3D три успешных взносы НАСА), что приводит к многочисленным коммерческим пакетам.

Иерархия уравнений потока жидкости

CFD можно рассматривать как группу вычислительных методологий (обсуждаемых ниже), используемых для решения уравнений, управляющих потоком жидкости. При применении CFD критическим шагом является решение, какой набор физических допущений и связанных уравнений необходимо использовать для решения данной проблемы. Чтобы проиллюстрировать этот шаг, ниже резюмируются физические допущения / упрощения, сделанные в уравнениях потока, который является однофазным (см. Многофазный поток и двухфазный поток ), однокомпонентным (т. Е. Он состоит из одного химического вещества), не -реагирующий и (если не указано иное) сжимаемый. Тепловым излучением пренебрегают и учитывают массовые силы, обусловленные гравитацией (если не указано иное). Кроме того, для этого типа потока следующее обсуждение подчеркивает иерархию уравнений потока, решаемых с помощью CFD. Обратите внимание, что некоторые из следующих уравнений можно вывести более чем одним способом.

  • Законы сохранения (CL): это самые фундаментальные уравнения, рассматриваемые с помощью CFD, в том смысле, что, например, все следующие уравнения могут быть выведены из них. Для однофазного сжимаемого потока, состоящего из одного вида, учитывают сохранение массы , сохранение количества движения и сохранение энергии .
  • Законы сохранения континуума (CCL): начните с CL. Предположим, что масса, импульс и энергия сохраняются локально : эти величины сохраняются и не могут «телепортироваться» из одного места в другое, а могут перемещаться только непрерывным потоком (см. Уравнение неразрывности ). Другая интерпретация состоит в том, что каждый начинает с CL и предполагает сплошную среду (см. Механику сплошной среды ). Полученная система уравнений является незамкнутой, поскольку для ее решения требуются дополнительные соотношения / уравнения: (а) определяющие соотношения для тензора вязких напряжений ; б) определяющие соотношения для диффузионного теплового потока ; (c) уравнение состояния (EOS), такое как закон идеального газа ; и (d) калорическое уравнение состояния, связывающее температуру с такими величинами, как энтальпия или внутренняя энергия .
  • Сжимаемые уравнения Навье-Стокса (C-NS): начните с CCL. Предположим, что ньютоновский тензор вязких напряжений (см. Ньютоновская жидкость ) и тепловой поток Фурье (см. Тепловой поток ). C-NS необходимо дополнить УС и калорическим УС, чтобы получить замкнутую систему уравнений.
  • Несжимаемые уравнения Навье-Стокса (I-NS): начните с C-NS. Предположим, что плотность всегда и везде постоянна. Другой способ получить I-NS - это предположить, что число Маха очень мало и что разница температур в жидкости также очень мала. В результате уравнения сохранения массы и импульса отделяются от уравнения сохранения энергии, поэтому нужно решать только первые два уравнения.
  • Сжимаемые уравнения Эйлера (EE): начните с C-NS. Предположим, что поток без трения и диффузного теплового потока отсутствует.
  • Слабо сжимаемые уравнения Навье-Стокса (WC-NS): начните с C-NS. Предположим, что изменение плотности зависит только от температуры, а не от давления. Например, для идеального газа используйте , где - это удобно определенное эталонное давление, которое всегда и везде постоянно, - плотность, - удельная газовая постоянная и - температура. В результате WK-NS не улавливает акустические волны. В WK-NS также часто пренебрегают членами, работающими под давлением, и вязким нагревом в уравнении сохранения энергии. WK-NS также называют C-NS с приближением малого числа Маха.
  • Уравнения Буссинеска: начните с C-NS. Предположим, что вариациями плотности всегда и везде можно пренебречь, за исключением гравитационного члена уравнения сохранения импульса (где плотность умножается на ускорение свободного падения). Также предположим, что различные свойства жидкости, такие как вязкость , теплопроводность и теплоемкость , всегда и везде постоянны. Уравнения Буссинеска широко используются в микромасштабной метеорологии .
  • Сжимаемые усредненные по Рейнольдсу уравнения Навье – Стокса и сжимаемые усредненные по Фавру уравнения Навье-Стокса (C-RANS и C-FANS): начните с C-NS. Предположим, что любая переменная потока , такая как плотность, скорость и давление, может быть представлена ​​как , где - среднее по ансамблю любой переменной потока, а - возмущение или отклонение от этого среднего значения. не обязательно маленький. Если является классическим усредненным по ансамблю (см. Разложение Рейнольдса ), можно получить усредненные по Рейнольдсу уравнения Навье – Стокса. И если это среднее по ансамблю, взвешенное по плотности, то получаются усредненные по Фавру уравнения Навье-Стокса. В результате и в зависимости от числа Рейнольдса диапазон масштабов движения значительно уменьшается, что приводит к гораздо более быстрым решениям по сравнению с решением C-NS. Однако информация теряется, и результирующая система уравнений требует закрытия различных незамкнутых членов, особенно напряжения Рейнольдса .
  • Уравнения идеального или потенциального потока : начните с EE. Предположим нулевое вращение жидких частиц (нулевая завихренность) и нулевое расширение потока (нулевая дивергенция). Результирующее поле потока полностью определяется геометрическими границами. Идеальные потоки могут быть полезны в современном CFD для инициализации моделирования.
  • Линеаризованные сжимаемые уравнения Эйлера (LEE): начните с EE. Предположим, что любая переменная потока , такая как плотность, скорость и давление, может быть представлена ​​как , где - значение переменной потока в некотором эталонном или базовом состоянии, а является возмущением или отклонением от этого состояния. Кроме того, предположим, что это возмущение очень мало по сравнению с некоторым эталонным значением. Наконец, предположим, что он удовлетворяет «собственному» уравнению, например EE. LEE и многие его разновидности широко используются в вычислительной аэроакустике .
  • Уравнение звуковой или акустической волны : начните с LEE. Пренебрегайте всеми градиентами и и предполагайте, что число Маха в исходном или базовом состоянии очень мало. Полученные в результате уравнения для плотности, количества движения и энергии можно преобразовать в уравнение давления, которое дает хорошо известное уравнение звуковой волны.
  • Уравнения мелкой воды (SW): рассмотрим поток у стены, где интересующий масштаб длины, параллельный стене, намного больше, чем интересующий масштаб длины нормали к стене. Начнем с EE. Предположим, что плотность всегда и везде постоянна, пренебрегаем составляющей скорости, перпендикулярной стенке, и считаем скорость, параллельную стенке, пространственно постоянной.
  • Уравнения пограничного слоя (BL): начните с C-NS (I-NS) для сжимаемых (несжимаемых) пограничных слоев. Предположим, что рядом со стенами есть тонкие области, где пространственные градиенты, перпендикулярные стене, намного больше, чем градиенты, параллельные стене.
  • Уравнение Бернулли: начните с EE. Предположим, что изменение плотности зависит только от изменения давления. См . Принцип Бернулли .
  • Устойчивое уравнение Бернулли: начните с уравнения Бернулли и предположите устойчивый поток. Или начните с EE и предположите, что поток устойчивый, и проинтегрируйте полученное уравнение вдоль линии тока.
  • Уравнения Стокса или ползущего потока: начните с C-NS или I-NS. Пренебрегайте инерцией потока. Такое предположение может быть оправдано, когда число Рейнольдса очень мало. В результате результирующая система уравнений является линейной, что значительно упрощает их решение.
  • Двумерное уравнение потока в канале: рассмотрим поток между двумя бесконечными параллельными пластинами. Начнем с C-NS. Предположим, что течение стационарное, двумерное и полностью развитое (т. Е. Профиль скорости не изменяется в продольном направлении). Обратите внимание, что это широко используемое полностью разработанное предположение может быть неадекватным в некоторых случаях, таких как некоторые сжимаемые микроканальные потоки, и в этом случае его можно заменить полностью разработанным на местном уровне предположением.
  • Одномерные уравнения Эйлера или одномерные уравнения газовой динамики (1D-EE): начните с EE. Предположим, что все величины потока зависят только от одного пространственного измерения.
  • Уравнение потока Фанно : рассмотрим поток внутри воздуховода с постоянной площадью и адиабатическими стенками. Начнем с 1D-EE. Предположите установившийся поток, отсутствие гравитационных эффектов, и введите в уравнение сохранения импульса эмпирический член для восстановления эффекта трения стенки (которым пренебрегают в EE). Чтобы замкнуть уравнение потока Фанно, необходима модель для этого члена трения. Такое закрытие включает в себя предположения, зависящие от проблемы.
  • Уравнение потока Рэлея . Рассмотрим поток внутри воздуховода с постоянной площадью и либо неадиабатическими стенками без объемных источников тепла, либо адиабатическими стенками с объемными источниками тепла. Начнем с 1D-EE. Предположите установившийся поток, отсутствие гравитационных эффектов, и введите в уравнение сохранения энергии эмпирический член для восстановления эффекта теплопередачи стенок или влияния источников тепла (которым пренебрегают в EE).

Методология

Во всех этих подходах соблюдается одна и та же основная процедура.

  • Во время предварительной обработки
    • В геометрии и физические границы проблемы могут быть определены с помощью системы автоматизированного проектирования (CAD). Оттуда данные можно соответствующим образом обработать (очистить) и извлечь объем текучей среды (или область текучей среды).
    • Объем , занимаемый жидкость разделяется на отдельные клетки (сетки). Сетка может быть однородной или неоднородной, структурированной или неструктурированной, состоящей из комбинации гексаэдрических, тетраэдрических, призматических, пирамидальных или многогранных элементов.
    • Определено физическое моделирование - например, уравнения движения жидкости + энтальпия + излучение + сохранение частиц.
    • Определены граничные условия. Это включает в себя определение поведения и свойств жидкости на всех ограничивающих поверхностях области жидкости. Для переходных задач также определены начальные условия.
  • Моделирование запускается и уравнения решаются итеративно как стационарные или переходный процесс .
  • Наконец, для анализа и визуализации полученного решения используется постпроцессор.

Методы дискретизации

Стабильность выбранной дискретизации обычно устанавливается численно, а не аналитически, как в простых линейных задачах. Особое внимание следует уделять тому, чтобы при дискретизации изящно обрабатывались прерывистые решения. Уравнения Эйлера и уравнения Навье – Стокса допускают удары и контактные поверхности.

Некоторые из используемых методов дискретизации:

Метод конечных объемов

Метод конечных объемов (FVM) - это общий подход, используемый в кодах CFD, поскольку он имеет преимущество в использовании памяти и скорости решения, особенно для больших задач, турбулентных потоков с высоким числом Рейнольдса и потоков с преобладанием источников (например, сгорание).

В методе конечных объемов основные дифференциальные уравнения в частных производных (обычно уравнения Навье-Стокса, уравнения сохранения массы и энергии и уравнения турбулентности) преобразовываются в консервативную форму, а затем решаются для дискретных контрольных объемов. Такая дискретизация гарантирует сохранение потоков через определенный контрольный объем. Уравнение конечного объема приводит к основным уравнениям в виде

где - вектор сохраняющихся переменных, - вектор потоков (см. уравнения Эйлера или уравнения Навье – Стокса ), - объем элемента контрольного объема и - площадь поверхности элемента контрольного объема.

Метод конечных элементов

Метод конечных элементов (МКЭ) используется в структурном анализе твердых тел, но также применим к жидкостям. Однако рецептура МКЭ требует особого внимания, чтобы гарантировать консервативное решение. Формулировка МКЭ адаптирована для использования с управляющими уравнениями гидродинамики. Хотя МКЭ должен быть тщательно сформулирован, чтобы быть консервативным, он намного более устойчив, чем подход конечных объемов. Однако FEM может потребовать больше памяти и иметь более медленное время решения, чем FVM.

В этом методе формируется взвешенное остаточное уравнение:

где - невязка уравнения в вершине элемента , - уравнение сохранения, выраженное на элементной основе, - весовой коэффициент и - объем элемента.

Метод конечных разностей

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

где есть вектор консервативных переменных, и , и являются потоки в , и направления соответственно.

Метод спектральных элементов

Метод спектральных элементов - это метод конечных элементов. Это требует, чтобы математическая задача (уравнение в частных производных) была сформулирована в слабой формулировке. Обычно это делается путем умножения дифференциального уравнения на произвольную тестовую функцию и интегрирования по всей области. Чисто математически тестовые функции совершенно произвольны - они принадлежат бесконечномерному функциональному пространству. Ясно, что бесконечномерное функциональное пространство не может быть представлено на дискретной сетке спектральных элементов; здесь начинается дискретизация спектрального элемента. Самое главное - это выбор функций интерполяции и тестирования. В стандартном МКЭ низкого порядка в 2D для четырехугольных элементов наиболее типичным выбором является билинейный тест или интерполирующая функция формы . Однако в методе спектральных элементов интерполяционные и тестовые функции выбираются как полиномы очень высокого порядка (обычно, например, 10-го порядка в приложениях CFD). Это гарантирует быструю сходимость метода. Кроме того, необходимо использовать очень эффективные процедуры интегрирования, поскольку количество интеграций, которые необходимо выполнить в числовых кодах, велико. Таким образом, используются квадратуры интегрирования Гаусса высокого порядка, поскольку они обеспечивают наивысшую точность при наименьшем количестве выполняемых вычислений. В настоящее время существует несколько академических кодов CFD, основанных на методе спектральных элементов, и еще несколько в настоящее время находятся в стадии разработки, поскольку в научном мире появляются новые схемы с временным шагом.

Решеточный метод Больцмана

Решеточный метод Больцмана (LBM) с его упрощенной кинетической картиной на решетке обеспечивает эффективное с вычислительной точки зрения описание гидродинамики. В отличие от традиционных методов CFD, которые решают уравнения сохранения макроскопических свойств (т. Е. Массы, импульса и энергии) численно, LBM моделирует жидкость, состоящую из фиктивных частиц, и такие частицы выполняют последовательные процессы распространения и столкновения по дискретной решетке. В этом методе используется дискретная в пространстве и времени версия кинетического эволюционного уравнения в форме Больцмана Бхатнагара-Гросса-Крука (BGK) .

Метод граничных элементов

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

Схемы дискретизации высокого разрешения

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

Модели турбулентности

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

Модели турбулентности можно классифицировать на основе вычислительных затрат, которые соответствуют диапазону масштабов, которые моделируются по сравнению с разрешенными (чем больше масштабов турбулентности будет разрешено, тем выше разрешение моделирования и, следовательно, тем выше стоимость вычислений). Если большая часть или все турбулентные масштабы не моделируются, вычислительные затраты будут очень низкими, но компромисс заключается в снижении точности.

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

Усредненное по Рейнольдсу Навье – Стокса

Внешняя аэродинамика модели DrivAer , рассчитанная с использованием URANS (вверху) и DDES (внизу)
Моделирование аэродинамического пакета Porsche Cayman (987.2) .

Усредненные по Рейнольдсу уравнения Навье – Стокса (RANS) являются старейшим подходом к моделированию турбулентности. Решается ансамблевая версия основных уравнений, которая вводит новые кажущиеся напряжения, известные как напряжения Рейнольдса . Это добавляет тензор неизвестных второго порядка, для которого различные модели могут обеспечивать разные уровни закрытия. Распространенное заблуждение состоит в том, что уравнения RANS неприменимы к потокам с изменяющимся во времени средним потоком, потому что эти уравнения являются «усредненными по времени». Фактически, одинаково можно рассматривать и статистически нестационарные (или нестационарные) потоки. Иногда это называют УРАН. В усреднении по Рейнольдсу нет ничего, что могло бы предотвратить это, но модели турбулентности, используемые для замыкания уравнений, действительны только до тех пор, пока время, в течение которого происходят эти изменения среднего значения, велико по сравнению с временными масштабами турбулентного движения, содержащего большую часть энергия.

Модели RANS можно разделить на два основных подхода:

Гипотеза Буссинеска
Этот метод включает использование алгебраического уравнения для напряжений Рейнольдса, которое включает определение турбулентной вязкости и, в зависимости от уровня сложности модели, решение уравнений переноса для определения турбулентной кинетической энергии и диссипации. Модели включают k-ε ( Лаундера и Сполдинга ), модель смешанной длины ( Прандтль ) и модель нулевого уравнения (Себечи и Смит ). Модели, доступные в этом подходе, часто называют количеством уравнений переноса, связанных с методом. Например, модель длины смешения является моделью «нулевого уравнения», потому что уравнения переноса не решаются; представляет собой модель «Два уравнения» , потому что два уравнения переноса (один для и один для ) решено.
Модель напряжения Рейнольдса (RSM)
Этот подход пытается фактически решить уравнения переноса для напряжений Рейнольдса. Это означает введение нескольких уравнений переноса для всех напряжений Рейнольдса, и, следовательно, такой подход требует гораздо больших затрат ресурсов центрального процессора.

Моделирование больших вихрей

Объемный рендеринг вихревого пламени без предварительного смешивания, моделируемый LES.

Моделирование крупных вихрей (LES) - это метод, в котором самые мелкие масштабы потока удаляются с помощью операции фильтрации, а их влияние моделируется с использованием моделей подсеточного масштаба. Это позволяет разрешить самые большие и наиболее важные масштабы турбулентности, при этом значительно снижая вычислительные затраты, связанные с наименьшими масштабами. Этот метод требует больших вычислительных ресурсов, чем методы RANS, но намного дешевле, чем DNS.

Моделирование отдельных вихрей

Моделирование отдельных вихрей (DES) - это модификация модели RANS, в которой модель переключается на формулировку подсеточного масштаба в областях, достаточно мелких для расчетов LES. Области вблизи твердых границ, где турбулентный масштаб длины меньше максимального размера сетки, получают режим решения RANS. Поскольку турбулентный масштаб длины превышает размер сетки, области решаются с использованием режима LES. Следовательно, разрешение сетки для DES не так требовательно, как для чистого LES, что значительно снижает стоимость вычислений. Хотя DES изначально был сформулирован для модели Спаларта-Аллмараса (Spalart et al., 1997), он может быть реализован с другими моделями RANS (Стрелец, 2001) путем соответствующего изменения масштаба длины, который явно или неявно задействован в модели RANS. . Таким образом, в то время как DES на основе модели Спаларта – Аллмараса действует как LES с моделью стены, DES, основанный на других моделях (например, модели с двумя уравнениями), ведет себя как гибридная модель RANS-LES. Генерация сети более сложна, чем в случае простого RANS или LES, из-за переключателя RANS-LES. DES - это незональный подход, который обеспечивает единое плавное поле скорости в областях RANS и LES решений.

Прямое численное моделирование

Прямое численное моделирование (DNS) разрешает весь диапазон турбулентных масштабов длины. Это минимизирует влияние моделей, но стоит очень дорого. Вычислительные затраты пропорциональны . DNS не поддается обработке для потоков со сложной геометрией или конфигурациями потоков.

Когерентное моделирование вихря

Подход когерентного вихревого моделирования разбивает поле турбулентного потока на когерентную часть, состоящую из организованного вихревого движения, и некогерентную часть, которая представляет собой случайный фоновый поток. Это разложение выполняется с использованием вейвлет- фильтрации. Этот подход имеет много общего с LES, поскольку он использует разложение и разрешает только отфильтрованную часть, но отличается тем, что не использует линейный фильтр нижних частот. Вместо этого операция фильтрации основана на вейвлетах, и фильтр можно адаптировать по мере развития поля потока. Фарж и Шнайдер протестировали метод CVS с двумя конфигурациями потока и показали, что когерентная часть потока демонстрирует энергетический спектр, демонстрируемый полным потоком, и соответствует когерентным структурам ( вихревым трубкам ), в то время как некогерентные части потока составляют однородный фон. шум, который не имеет организованных структур. Гольдштейн и Васильев применили модель FDV к моделированию больших вихрей, но не предполагали, что вейвлет-фильтр полностью устраняет все когерентные движения из масштабов подфильтра. Используя фильтрацию как LES, так и CVS, они показали, что в диссипации SFS преобладает когерентная часть поля потока SFS.

PDF методы

Методы функции плотности вероятности (PDF) для турбулентности, впервые представленные Лундгреном , основаны на отслеживании одноточечной PDF скорости , которая дает вероятность того, что скорость в точке находится между и . Этот подход аналогичен кинетической теории газов , в которой макроскопические свойства газа описываются большим числом частиц. Методы PDF уникальны тем, что их можно применять в рамках ряда различных моделей турбулентности; основные отличия заключаются в форме уравнения переноса PDF. Например, в контексте моделирования больших вихрей PDF-файл становится фильтрованным PDF-файлом. Методы PDF также могут использоваться для описания химических реакций и особенно полезны для моделирования химически реагирующих потоков, поскольку термин химический источник является замкнутым и не требует модели. PDF обычно отслеживается с помощью методов лагранжевых частиц; в сочетании с моделированием крупных вихрей это приводит к уравнению Ланжевена для эволюции частиц подфильтра.

Вихревой метод

Вихревой метод представляет собой бессеточный метод моделирования турбулентных течений. Он использует вихри в качестве вычислительных элементов, имитируя физические структуры в турбулентности. Вихревые методы были разработаны как безсеточная методология, которая не будет ограничена фундаментальными эффектами сглаживания, связанными с сеточными методами. Однако, чтобы быть практичными, вихревые методы требуют средств для быстрого вычисления скоростей от вихревых элементов - другими словами, они требуют решения конкретной формы задачи N тел (в которой движение N объектов связано с их взаимными влияниями. ). Прорыв произошел в конце 1980-х годов с разработкой алгоритма быстрого мультиполя (FMM), разработанного В. Рохлиным (Йельский университет) и Л. Грингардом (Институт Куранта). Этот прорыв проложил путь к практическому вычислению скоростей от вихревых элементов и является основой успешных алгоритмов.

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

  • Он практически не имеет сетки, что исключает многочисленные итерации, связанные с RANS и LES.
  • Все проблемы рассматриваются одинаково. Никаких вводных данных для моделирования или калибровки не требуется.
  • Возможно моделирование временных рядов, которые имеют решающее значение для правильного анализа акустики.
  • Малый и крупный масштабы точно моделируются одновременно.

Метод удержания завихренности

Метод удержания завихренности (ВК) - это метод Эйлера, используемый при моделировании турбулентных следов. Он использует подход, подобный уединенной волне, для получения стабильного решения без численного расширения. VC может захватывать мелкомасштабные объекты с точностью до 2 ячеек сетки. В рамках этих функций решается нелинейное разностное уравнение в отличие от конечно-разностного уравнения . ВК аналогичен методам захвата ударных волн , где выполняются законы сохранения, так что основные интегральные величины вычисляются точно.

Линейная модель вихря

Модель линейных вихрей - это метод, используемый для моделирования конвективного перемешивания, происходящего в турбулентном потоке. В частности, он предоставляет математический способ описания взаимодействий скалярной переменной в векторном поле потока. Он в основном используется в одномерном представлении турбулентного потока, поскольку может применяться в широком диапазоне масштабов длины и чисел Рейнольдса. Эта модель обычно используется в качестве строительного блока для более сложных представлений потока, поскольку она обеспечивает прогнозы с высоким разрешением, которые справедливы для большого диапазона условий потока.

Двухфазный поток

Моделирование орды пузырей методом объемной жидкости

Моделирование двухфазного потока все еще находится в стадии разработки. Были предложены различные методы, включая метод объема жидкости, метод установки уровня и отслеживание фронта . Эти методы часто предполагают компромисс между поддержанием четкой границы раздела или сохранением массы. Это очень важно, поскольку оценка плотности, вязкости и поверхностного натяжения основана на значениях, усредненных по границе раздела. Лагранжевые многофазные модели, которые используются для дисперсных сред, основаны на решении лагранжевого уравнения движения для дисперсной фазы.

Алгоритмы решения

Дискретизация в пространстве дает систему обыкновенных дифференциальных уравнений для нестационарных задач и алгебраических уравнений для стационарных задач. Неявные или полунеявные методы обычно используются для интегрирования обыкновенных дифференциальных уравнений, создавая систему (обычно) нелинейных алгебраических уравнений. Применение итераций Ньютона или Пикара дает систему линейных уравнений, несимметричную при наличии адвекции и неопределенную при наличии несжимаемости. Такие системы, особенно в 3D, часто слишком велики для прямых решателей, поэтому используются итерационные методы, либо стационарные методы, такие как последовательная сверхрелаксация, либо методы подпространства Крылова . Методы Крылова, такие как GMRES , обычно используемые с предварительным условием , работают, минимизируя невязку по последовательным подпространствам, генерируемым оператором предобусловливания.

Преимущество Multigrid - асимптотически оптимальная производительность при решении многих задач. Традиционные решатели и предобуславливатели эффективны для уменьшения высокочастотных компонентов остатка, но для уменьшения низкочастотных компонентов обычно требуется много итераций. Работая в нескольких масштабах, multigrid уменьшает все компоненты остатка на аналогичные коэффициенты, что приводит к независимому от сетки количеству итераций.

Для неопределенных систем предварительные кондиционеры, такие как неполная факторизация LU , аддитивная Шварца и многосеточные, работают плохо или полностью выходят из строя, поэтому для эффективного предварительного кондиционирования необходимо использовать структуру задачи. Методами, обычно используемыми в CFD, являются алгоритмы SIMPLE и Uzawa, которые демонстрируют скорость сходимости, зависящую от сетки, но недавние достижения, основанные на блочной факторизации LU в сочетании с многосеточной факторизацией для результирующих определенных систем, привели к предварительным кондиционерам, которые обеспечивают скорость сходимости, не зависящую от сетки.

Неустойчивая аэродинамика

CFD сделала большой прорыв в конце 70-х годов с введением LTRAN2, двумерного кода для моделирования колеблющихся профилей, основанного на трансзвуковой теории малых возмущений, разработанной Ballhaus и соавторами. Он использует алгоритм переключения Мурмана-Коула для моделирования движущихся ударных волн. Позже AFWAL / Boeing расширили его до 3-D с использованием развернутой разностной схемы, что привело к LTRAN3.

Биомедицинская инженерия

Моделирование кровотока в аорте человека

CFD-исследования используются для уточнения характеристик аортального кровотока в деталях, которые выходят за рамки возможностей экспериментальных измерений. Для анализа этих условий извлекаются CAD-модели сосудистой системы человека с использованием современных методов визуализации, таких как МРТ или компьютерная томография . На основе этих данных реконструируется трехмерная модель, и может быть рассчитан поток жидкости. Необходимо учитывать такие свойства крови, как плотность и вязкость, а также реалистичные граничные условия (например, системное давление). Таким образом, можно анализировать и оптимизировать поток в сердечно-сосудистой системе для различных приложений.

CPU против GPU

Традиционно моделирование CFD выполняется на процессорах. В последнее время моделирование также выполняется на графических процессорах. Обычно они содержат более медленные, но больше процессоров. Для алгоритмов CFD, которые обладают хорошей производительностью параллелизма (т.е. хорошим ускорением за счет добавления большего количества ядер), это может значительно сократить время моделирования. Неявные жидкие частицы и методы решетки-Больцмана являются типичными примерами кодов, которые хорошо масштабируются на графических процессорах.

Смотрите также

использованная литература

Примечания

  • Андерсон, Джон Д. (1995). Вычислительная гидродинамика: основы с приложениями . Наука / Инженерия / Математика. McGraw-Hill Science. ISBN 978-0-07-001685-9.
  • Патанкар, Сухас (1980). Числовой теплообмен и поток жидкости . Серия Hemisphere по вычислительным методам в механике и теплотехнике. Тейлор и Фрэнсис. ISBN 978-0-89116-522-4.

внешние ссылки