Вычислительная аэроакустика

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

История

Зарождение вычислительной аэроакустики, скорее всего, можно отнести только к середине 1980-х годов, когда появилась публикация Хардина и Ламкина. который утверждал, что

« … область вычислительной механики жидкости быстро развивалась в последние несколько лет и теперь дает надежду, что «вычислительная аэроакустика», где шум вычисляется непосредственно на основе первых принципов определения непрерывных полей скорости и завихренности, могла бы быть возможным, [...] »

Позже в публикации 1986 г. те же авторы ввели аббревиатуру CAA. Первоначально этот термин использовался для подхода с низким числом Маха (расширение поля акустических возмущений вокруг несжимаемого потока), как это описано в разделе EIF . Позже, в начале 1990-х годов, растущее сообщество CAA подхватило этот термин и широко использовало его для любого численного метода, описывающего шумовое излучение аэроакустического источника или распространение звуковых волн в неоднородном поле потока. Такими численными методами могут быть методы интегрирования в дальней зоне (например, FW-H ), а также прямые численные методы, оптимизированные для решения (например, ) математической модели, описывающей возникновение и/или распространение аэродинамического шума. Благодаря быстрому развитию вычислительных ресурсов в этой области за последние три десятилетия произошел впечатляющий прогресс.

Методы

Подход прямого численного моделирования (DNS) к CAA

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

Гибридный подход

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

1. Интегральные методы
   1. Аналогия Лайтхилла
   2. Интеграл Кирхгофа
   3. FW-H
2. ЛИ
3. Псевдоспектральный
4. ЕИФ
5. ОБЕЗЬЯНА

Интегральные методы

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

Аналогия Лайтхилла

Также называется « Акустическая аналогия ». Чтобы получить аэроакустическую аналогию Лайтхилла, основные уравнения Навье-Стокса перестраиваются. Левая часть представляет собой волновой оператор, который применяется к возмущению плотности или возмущению давления соответственно. Тогда правая часть идентифицируется как источники звука в потоке жидкости. Поскольку аналогия Лайтхилла следует непосредственно из уравнений Навье-Стокса без упрощений, все источники присутствуют. Некоторые из источников затем идентифицируются как турбулентный или ламинарный шум. Звуковое давление в дальней зоне затем выражается как объемный интеграл по области, содержащей источник звука. Источниковый термин всегда включает в себя физические источники и такие источники, которые описывают распространение в неоднородной среде.

Волновой оператор аналогии Лайтхилла ограничен условиями постоянного потока за пределами зоны источника. Никакое изменение плотности, скорости звука и числа Маха не допускается. Различные условия среднего потока по аналогии идентифицируются как сильные источники с противоположным знаком, как только акустическая волна проходит через него. Часть акустической волны удаляется одним источником и излучается новая волна, чтобы зафиксировать другую скорость волны. Это часто приводит к очень большим объемам с сильными источниками. Было предложено несколько модификаций оригинальной теории Лайтхилла для объяснения взаимодействия звука и потока или других эффектов. Чтобы улучшить аналогию Лайтхилла, различные величины внутри волнового оператора, а также разные волновые операторы рассматриваются с помощью следующих аналогий. Все они получают измененные исходные термины, которые иногда позволяют более четко понять «настоящие» источники. Акустические аналогии Лиллея, Пирс, Хау и Мёринг Это лишь некоторые примеры аэроакустических аналогий, основанных на идеях Лайтхилла. Все акустические аналогии требуют объемного интегрирования по исходному термину.

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

Интеграл Кирхгофа

Кирхгоф и Гельмгольц показали, что излучение звука из ограниченной области источника можно описать, окружив эту область источника управляющей поверхностью — так называемой поверхностью Кирхгофа. Тогда звуковое поле внутри или снаружи поверхности, где не допускаются никакие источники и применяется волновой оператор слева, может быть создано как суперпозиция монополей и диполей на поверхности. Теория следует непосредственно из волнового уравнения. Силу источника монополей и диполей на поверхности можно рассчитать, если известны нормальная скорость (для монополей) и давление (для диполей) на поверхности соответственно. Модификация метода позволяет даже рассчитывать давление на поверхность, исходя только из нормальной скорости. Нормальную скорость можно определить, например, с помощью FE-моделирования движущейся конструкции. Однако модификация, позволяющая избежать известного акустического давления на поверхности, приводит к проблемам при рассмотрении замкнутого объема на его резонансных частотах, что является основной проблемой реализации их метода. Интегральный метод Кирхгофа находит применение, например, в Методы граничных элементов (МГЭ). Ненулевая скорость потока учитывается путем рассмотрения движущейся системы отсчета со скоростью внешнего потока, в которой происходит распространение акустической волны. Повторяющиеся применения метода могут привести к возникновению препятствий. Сначала рассчитывается звуковое поле на поверхности препятствия, а затем вводится препятствие путем добавления источников на его поверхности, чтобы исключить нормальную скорость на поверхности препятствия. Изменения среднего поля потока (скорости звука, плотности и скорости) можно учесть аналогичным методом (например, двойной взаимностью BEM).

FW-H

Метод интеграции Фоукса Уильямса и Хокингса основан на акустической аналогии Лайтхилла. Однако путем некоторых математических модификаций в предположении об ограниченной области источника, которая окружена управляющей поверхностью (поверхность FW-H), объемного интеграла можно избежать. Поверхностные интегралы по монопольным и дипольным источникам остаются. В отличие от метода Кирхгофа, эти источники следуют непосредственно из уравнений Навье-Стокса по аналогии Лайтхилла. Источники вне поверхности FW-H можно учесть с помощью дополнительного объемного интеграла по квадрупольным источникам, следующего из тензора Лайтхилла. Однако при рассмотрении тех же предположений, что и в линейной теории Кирхгофа, метод FW-H равен методу Кирхгофа.

Линеаризованные уравнения Эйлера

Учитывая небольшие возмущения, наложенные на однородный поток средней плотности $\rho _{0}$, давление $p_{0}$ и скорость по оси X $u_{0}$, уравнения Эйлера для двумерной модели представлены в виде:

${\frac {\partial \mathbf {U} }{\partial t}}+{\frac {\partial \mathbf {F} }{\partial x}}+{\frac {\partial \mathbf {G} }{\partial y}}=\mathbf {S}$,

где

$\mathbf {U} ={\begin{bmatrix}\rho \\u\\v\\p\\\end{bmatrix}}\ ,\ \mathbf {F} ={\begin{bmatrix}\rho _{0}u+\rho u_{0}\\u_{0}u+p/\rho _{0}\\u_{0}v\\u_{0}p+\gamma p_{0}u\\\end{bmatrix}}\ ,\ \mathbf {G} ={\begin{bmatrix}\rho _{0}v\\0\\p/\rho _{0}\\\gamma p_{0}v\\\end{bmatrix}},$

где $\rho$, $u$, $v$ и $p$ – переменные акустического поля, $\gamma$ соотношение удельных теплоемкостей $c_{p}/c_{v}$, для воздуха при 20 °C $c_{p}/c_{v}=1.4$, и исходный термин $\mathbf {S}$ справа представлены распределенные нестационарные источники. Применение LEE можно найти в исследованиях шума двигателей.

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

Псевдоспектральный

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

ЕИФ

Расширение о несжимаемом потоке

ОБЕЗЬЯНА

Уравнения акустических возмущений

См. статью Р. Эверта и В. Шредера «Уравнения акустических возмущений, основанные на разложении потока посредством фильтрации источника».

См. также

• Аэроакустика
• Акустическая теория

Источники

• Лайтхилл, М.Дж., «Общее введение в аэроакустику и атмосферные звуки», отчет ICASE 92-52, Исследовательский центр НАСА в Лэнгли, Хэмптон, Вирджиния , 1992 г.

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

• Примеры в аэроакустике от НАСА
• Вычислительная аэроакустика в Центральной школе Лиона
• Вычислительная аэроакустика в Левенском университете
• Вычислительная аэроакустика в Берлинском техническом университете
• Сценарий лекции CAA из Технического университета Берлина.