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

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

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

Предыстория и история

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

Моделирование ГПВРД Hyper-X в работе на Маха . скорости -7

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

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

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

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

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

В двумерной области был разработан ряд панельных кодов для анализа и проектирования профиля крыла. Коды обычно включают анализ пограничного слоя , что позволяет моделировать вязкие эффекты. Ричард Эпплер разработал код PROFILE, частично при финансовой поддержке НАСА, который стал доступен в начале 1980-х годов. Вскоре за этим последовал Марка Дрела код XFOIL . И PROFILE, и XFOIL включают двумерные панельные коды со связанными кодами пограничного слоя для анализа профиля профиля. PROFILE использует метод конформного преобразования для проектирования обратного профиля крыла, в то время как XFOIL имеет как конформное преобразование, так и метод обратной панели для проектирования профиля крыла.

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

Моделирование космического корабля SpaceX во время входа в атмосферу

Разработчики обратились к кодам Full Potential, поскольку панельные методы не могли рассчитать нелинейный поток, присутствующий на околозвуковых скоростях. Первое описание способа использования уравнений полного потенциала было опубликовано Эрлом Мурманом и Джулианом Коулом из Boeing в 1970 году. Фрэнсис Бауэр, Пол Гарабедян и Дэвид Корн из Института Куранта при Нью-Йоркском университете (NYU) написали серию двумерных кодов профиля профиля полного потенциала, которые получили широкое распространение, наиболее важным из которых была названа программа 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 Энтони Джеймсон также разработал трехмерный код САМОЛЕТА. в котором использовались неструктурированные тетраэдральные сетки.

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

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

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

Иерархия

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

• Законы сохранения (CL). Это наиболее фундаментальные уравнения, рассматриваемые с помощью CFD, в том смысле, что, например, из них можно вывести все следующие уравнения. Для однофазного одновидового сжимаемого потока учитываются сохранение массы , сохранение погонного момента и сохранение энергии .
• Законы сохранения континуума (CCL): начните с CL. Предположим, что масса, импульс и энергия локально сохраняются: эти величины сохраняются и не могут «телепортироваться» из одного места в другое, а могут перемещаться только непрерывным потоком (см. уравнение неразрывности ). Другая интерпретация состоит в том, что каждый начинает с CL и предполагает наличие сплошной среды (см. Механику сплошной среды ). Полученная система уравнений является незамкнутой, так как для ее решения необходимы дополнительные соотношения/уравнения: а) определяющие соотношения для тензора вязких напряжений ; (б) определяющие соотношения для диффузионного теплового потока ; (c) уравнение состояния (EOS), такое как закон идеального газа ; и (d) уравнение состояния калорий, связывающее температуру с такими величинами, как энтальпия или внутренняя энергия .
• Сжимаемые уравнения Навье-Стокса (C-NS): начните с CCL. Предположим, что существует ньютоновский тензор вязких напряжений (см. Ньютоновская жидкость ) и тепловой поток Фурье (см. Тепловой поток ). C-NS необходимо дополнить EOS и калорическим EOS, чтобы получить замкнутую систему уравнений.
• Несжимаемые уравнения Навье-Стокса (I-NS): начните с C-NS. Предположим, что плотность всегда и везде постоянна. Другой способ получить I-NS — предположить, что число Маха очень мало. и что разница температур в жидкости также очень мала. В результате уравнения сохранения массы и сохранения импульса отделяются от уравнения сохранения энергии, поэтому нужно решить только первые два уравнения.
• Сжимаемые уравнения Эйлера (EE): начните с C-NS. Предположим, что существует течение без трения и без диффузионного теплового потока.
• Слабо сжимаемые уравнения Навье-Стокса (WC-NS): начните с C-NS. Предположим, что изменения плотности зависят только от температуры, а не от давления. Например, для идеального газа используйте $\rho =p_{0}/(RT)$, где $p_{0}$ представляет собой удобно определяемое эталонное давление, всегда и везде постоянное, $\rho$ плотность, $R$ – удельная газовая постоянная , $T$ это температура. В результате WC-NS не улавливает акустические волны. В WC-NS также часто пренебрегают условиями работы давления и вязкого нагрева в уравнении сохранения энергии. WC-NS также называют C-NS с приближением малого числа Маха.
• Уравнения Буссинеска: начнем с C-NS. Предположим, что изменения плотности всегда и везде пренебрежимо малы, за исключением гравитационного члена уравнения сохранения импульса (где плотность умножает гравитационное ускорение). Также предположим, что различные свойства жидкости, такие как вязкость , теплопроводность и теплоемкость , всегда и везде постоянны. Уравнения Буссинеска широко используются в микромасштабной метеорологии .
• Сжимаемые уравнения Навье-Стокса, усредненные по Рейнольдсу , и сжимаемые уравнения Навье-Стокса, усредненные по Фавре (C-RANS и C-FANS): начните с C-NS. Предположим, что любая переменная потока $f$, такие как плотность, скорость и давление, можно представить как $f=F+f''$, где $F$ среднее по ансамблю любой переменной потока, и $f''$ является отклонением или отклонением от этого среднего значения. $f''$ не обязательно мал. Если $F$ является классическим средним по ансамблю (см. разложение Рейнольдса ), получают усредненные по Рейнольдсу уравнения Навье – Стокса. И если $F$ является взвешенным по плотности средним по ансамблю, получаются усредненные по Фавре уравнения Навье-Стокса. В результате, в зависимости от числа Рейнольдса, диапазон масштабов движения значительно сокращается, что приводит к гораздо более быстрому решению по сравнению с решением C-NS. Однако информация теряется, и полученная система уравнений требует замыкания различных незамкнутых членов, в частности, напряжения Рейнольдса .
• Уравнения идеального потока или потенциального потока : начните с EE. Предположим, что вращение жидкости и частиц равно нулю (нулевая завихренность) и нулевое расширение потока (нулевая дивергенция). Результирующее поле течения полностью определяется геометрическими границами. Идеальные потоки могут быть полезны в современных CFD для инициализации моделирования.
• Линеаризованные уравнения Эйлера сжимаемой жидкости (LEE): Начните с ЭЭ. Предположим, что любая переменная потока $f$, такие как плотность, скорость и давление, можно представить как $f=f_{0}+f'$, где $f_{0}$ - значение переменной потока в некотором эталонном или базовом состоянии, и $f'$ является возмущением или отклонением от этого состояния. Далее, предположим, что это возмущение $f'$ очень мала по сравнению с некоторым эталонным значением. Наконец, предположим, что $f_{0}$ удовлетворяет «своему» уравнению, такому как EE. LEE и его многочисленные вариации широко используются в вычислительной аэроакустике .
• звуковой волны или Уравнение акустической волны : начните с LEE. Пренебрегайте всеми градиентами $f_{0}$ и $f'$и предположим, что число Маха в исходном или базовом состоянии очень мало. Полученные уравнения для плотности, импульса и энергии можно преобразовать в уравнение давления, что дает хорошо известное уравнение звуковой волны.
• Уравнения мелкой воды (SW). Рассмотрим течение вблизи стенки, где интересующий масштаб длины, параллельной стене, намного больше, чем интересующий масштаб длины, перпендикулярный стене. Начните с ЭЭ. Предположим, что плотность всегда и везде постоянна, пренебрегаем составляющей скорости, перпендикулярной стенке, и считаем скорость, параллельную стене, пространственно-постоянной.
• Уравнения пограничного слоя (BL): начните с C-NS (I-NS) для сжимаемых (несжимаемых) пограничных слоев. Предположим, что рядом со стенами есть тонкие области, где пространственные градиенты, перпендикулярные стене, намного больше, чем параллельные стене.
• Уравнение Бернулли: начнем с ЭЭ. Предположим, что изменения плотности зависят только от изменений давления. См. Принцип Бернулли .
• Устойчивое уравнение Бернулли. Начните с уравнения Бернулли и предположите устойчивый поток. Или начните с EE и предположите, что поток устойчив, и проинтегрируйте полученное уравнение вдоль линии тока.
• Уравнения потока Стокса или ползущего потока: начните с C-NS или I-NS. Пренебрегаем инерцией потока. Такое предположение может быть оправдано, когда число Рейнольдса очень мало. В результате результирующая система уравнений является линейной, что значительно упрощает их решение.
• Двумерное уравнение потока в канале: рассмотрим поток между двумя бесконечными параллельными пластинами. Начните с ЦНС. Предположим, что течение стационарное, двумерное и полностью развитое (т. е. профиль скорости не меняется в продольном направлении). Обратите внимание, что это широко используемое полностью разработанное предположение может быть неадекватным в некоторых случаях, например, в некоторых сжимаемых микроканальных течениях, и в этом случае оно может быть заменено локально полностью разработанным предположением.
• Одномерные уравнения Эйлера или одномерные уравнения газовой динамики (1D-EE): начните с EE. Предположим, что все величины потока зависят только от одного пространственного измерения.
• Уравнение потока Фанно : рассмотрим поток внутри воздуховода с постоянной площадью и адиабатическими стенками. Начните с 1D-EE. Предположим, что поток устойчивый, нет эффектов гравитации, и введите в уравнение сохранения импульса эмпирический член, чтобы восстановить эффект трения о стенки (пренебрегаемый в EE). Чтобы закрыть уравнение потока Фанно, необходима модель для этого члена трения. Такое замыкание включает в себя предположения, зависящие от проблемы.
• Уравнение потока Рэлея . Рассмотрим течение внутри канала постоянной площади и либо с неадиабатическими стенками без объемных источников тепла, либо с адиабатическими стенками с объемными источниками тепла. Начните с 1D-EE. Предположим, что поток устойчивый, нет эффектов гравитации, и введите в уравнение сохранения энергии эмпирический член, чтобы восстановить эффект теплопередачи стенок или эффект источников тепла (пренебрегаемый в EE).

Методология

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

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

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

Дополнительная информация: Дискретизация уравнений Навье – Стокса.

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

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

Метод конечного объема

Основная статья: Метод конечного объема

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

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

${\frac {\partial }{\partial t}}\iiint Q\,dV+\iint F\,d\mathbf {A} =0,$

где $Q$ вектор сохраняющихся переменных, $F$ – вектор потоков (см. уравнения Эйлера или уравнения Навье – Стокса ), $V$ - объем элемента контрольного объема, а $\mathbf {A}$ – площадь поверхности элемента контрольного объема.

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

Основная статья: Метод конечных элементов

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

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

$R_{i}=\iiint W_{i}Q\,dV^{e}$

где $R_{i}$ - невязка уравнения в вершине элемента $i$, $Q$ — уравнение сохранения, выраженное на элементной основе, $W_{i}$ - весовой коэффициент, а $V^{e}$ — объем элемента.

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

Основная статья: Метод конечных разностей

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

${\frac {\partial Q}{\partial t}}+{\frac {\partial F}{\partial x}}+{\frac {\partial G}{\partial y}}+{\frac {\partial H}{\partial z}}=0$

где $Q$ — вектор сохраняющихся переменных, а $F$, $G$, и $H$ являются потоками в $x$, $y$, и $z$ направления соответственно.

Метод спектрального элемента

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

Метод спектральных элементов представляет собой метод типа конечных элементов. Для этого требуется, чтобы математическая задача (уравнение в частных производных) была представлена ​​в слабой формулировке. Обычно это делается путем умножения дифференциального уравнения на произвольную тестовую функцию и интегрирования по всей области. Чисто математически тестовые функции совершенно произвольны — они принадлежат бесконечномерному функциональному пространству. Очевидно, что бесконечномерное функциональное пространство не может быть представлено на сетке дискретных спектральных элементов; здесь начинается дискретизация спектрального элемента. Самым важным является выбор интерполяционных и проверочных функций. В стандартном МКЭ низкого порядка в 2D для четырехугольных элементов наиболее типичным выбором является билинейный тест или интерполирующая функция вида $v(x,y)=ax+by+cxy+d$. Однако в методе спектральных элементов интерполирующие и тестовые функции выбираются в виде полиномов очень высокого порядка (обычно, например, 10-го порядка в приложениях CFD). Это гарантирует быструю сходимость метода. Кроме того, необходимо использовать очень эффективные процедуры интегрирования, поскольку число интегрирований, выполняемых в числовых кодах, велико. Таким образом, используются квадратуры интегрирования Гаусса высокого порядка, поскольку они достигают наибольшей точности при наименьшем количестве выполняемых вычислений.В настоящее время существует несколько академических CFD-кодов, основанных на методе спектральных элементов, а еще несколько находятся в стадии разработки, поскольку в научном мире возникают новые схемы временного шага.

метод

Основная статья: Решеточные методы Больцмана

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

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

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

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

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

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

Основная статья: Метод граничных элементов

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

высокого

Основная статья: Схема высокого разрешения

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

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

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

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

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

Рейнольдсу Навье –

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

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

Моделирование аэродинамического пакета Porsche Cayman (987.2)

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

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

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

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

Модель стресса Рейнольдса (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 решений.

IDDES Моделирование BMW от Karel Motorsports. Это тип моделирования DES, выполненный в OpenFOAM. График представляет собой коэффициент давления.

численное

Основная статья: Прямое численное моделирование

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

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

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

PDF-методы

Методы функции плотности вероятности (PDF) для турбулентности, впервые представленные Лундгреном , основаны на отслеживании одноточечной PDF скорости, $f_{V}({\boldsymbol {v}};{\boldsymbol {x}},t)d{\boldsymbol {v}}$, что дает вероятность скорости в точке ${\boldsymbol {x}}$ находясь между ${\boldsymbol {v}}$ и ${\boldsymbol {v}}+d{\boldsymbol {v}}$. Этот подход аналогичен кинетической теории газов , в которой макроскопические свойства газа описываются большим количеством частиц. Методы PDF уникальны тем, что их можно применять в рамках множества различных моделей турбулентности; основные различия заключаются в форме уравнения переноса PDF. Например, в контексте моделирования больших вихрей PDF-файл становится отфильтрованным PDF-файлом. Методы PDF также можно использовать для описания химических реакций, и особенно полезны для моделирования химически реагирующих потоков, поскольку член химического источника является закрытым и не требует модели. PDF обычно отслеживается с использованием методов лагранжевых частиц; в сочетании с моделированием больших вихрей это приводит к уравнению Ланжевена для эволюции частиц подфильтра.

Метод ограничения завихренности

Основная статья: Удержание завихренности

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

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

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

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

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

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

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

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

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

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

аэродинамика

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

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

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

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

Процессор против графического процессора

Традиционно моделирование CFD выполняется на центральных процессорах.

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

См. также

• Теория элемента лезвия
• Граничные условия в гидродинамике
• Моделирование кавитации
• Центральная разностная схема
• Вычислительная магнитогидродинамика
• Метод дискретных элементов
• Метод конечных элементов
• Метод конечных объемов для нестационарного течения
• Плавная анимация
• Метод погруженных границ
• Решеточные методы Больцмана
• Список пакетов программного обеспечения для конечных элементов
• Бессеточные методы
• Полунеявный метод движущихся частиц
• Динамика многочастичных столкновений
• Многопрофильная оптимизация проектирования
• Численные методы в механике жидкости
• Оптимизация формы
• Гидродинамика сглаженных частиц
• Стохастический метод Эйлера Лагранжа
• Моделирование турбулентности
• Визуализация (графика)
• Аэродинамическая труба

Примечания

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

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

• Курс: Вычислительная гидродинамика – Суман Чакраборти ( Индийский технологический институт Харагпур )
• Курс: Численные методы PDE для ученых и инженеров , Лекции в открытом доступе и коды для числовых PDE, включая современный взгляд на сжимаемую CFD