Дифференциальные уравнения в акустической инженерии: применение звукового распространения и шумоконтроля

∂²u/∂t² = c² ⁇ ²u

Здесь u обозначает звуковое давление, ttcc ⁇ ² — лаплакийский оператор, представляющий сумму вторых пространственных производных. Это простое, но мощное выражение предполагает, что среда однородна, изотропна и не имеет потерь. На практике реальные среды вводят сложности, такие как вязкость, теплопроводность и турбулентность, которые изменяют волновое уравнение в более общую форму — обобщенное волновое уравнение или уравнение лужайкой волны.

Происхождение из физических принципов

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

• Уравнение непрерывности: ∂ρ/∂t + ⁇ ·(ρv) = 0, где ρ — плотность, а v — скорость частиц.
• Уравнение Эйлера: ρ ∂v/∂t + ⁇ p = 0 для потока невязкости, при этом p представляет давление.
• Уравнение состояния: p = c² ρ для малых возмущений.

Устранение скорости и плотности приводит к волновому уравнению выше. Это вывод подчеркивает прямую связь между фундаментальной физикой и математическими моделями, используемыми в акустической инженерии.

Варианты волнового уравнения

Для различных применений требуются изменения стандартного волнового уравнения.

• Уравнение Гельмгольца: Для устойчиво-устойчивых (временно-гармонических) звуковых полей временная часть отделяется, давая ⁇ ²p + k²p = 0, где k = ω/c — волновое число, а ω — угловая частота.
• Связанное волновое уравнение: Учитывает средний поток (например, ветер или вентиляцию) путем включения термина, включающего градиент скорости.
• Модифицированное пористостью волновое уравнение: Для распространения звука через пористые материалы включены дополнительные демпфирующие и инерционные термины.

Каждый вариант позволяет инженерам адаптировать аналитическое или численное решение к конкретной физической ситуации.

Решение волнового уравнения: аналитические и численные методы

В то время как волновое уравнение имеет решения замкнутой формы для простых геометрий (например, бесконечные плоские волны, сферические волны в свободном пространстве, волны в прямоугольных комнатах), приложения реального мира часто требуют численных подходов.Два широко используемых метода - метод конечных элементов (FEM) и метод граничных элементов (BEM) [FLT: 3].

Метод конечных элементов (FEM) в акустике

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

• Моделирование сложных внутренних пространств (комнат, транспортных средств, кабин самолетов).
• Включая импедансные граничные условия для поглощающих стен.
• Связанные структурно-акустические проблемы (например, передача звука через панели).

Коммерческое программное обеспечение, такое как COMSOL Multiphysics ^® и Ansys ^® , обеспечивает надежные FEM-решатели для акустики, позволяя инженерам визуализировать распределение давления и оптимизировать дизайн итеративно.

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

BEM сводит проблему только к границе, используя функции Грина для связи давления и нормальной скорости на поверхности. Этот метод идеально подходит для внешних проблем (например, шума, излучаемого вентилятором или двигателем), потому что он по своей сути удовлетворяет радиационному состоянию Зоммерфельда (распад волн на бесконечности). BEM требует решения плотной системы уравнений, но он избегает сетки полного объема, что делает его эффективным для больших внешних доменов.

Методы трассировки лучей и энергетические методы

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

Приложения для контроля шума

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

Звуковая абсорбция и уравнение Гельмгольца

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

∂p/ ∂n = -j k β p, где β — специфическое акустическое допускание.

Итерационно изменяя различные толщины и расположения материала, инженеры могут минимизировать уровни звукового давления на целевых частотах. Например, в промышленном канале глушитель, спроектированный с использованием этих уравнений, может ослаблять шум вентилятора на 10-20 дБ в широком диапазоне частот.

Активный шумоконтроль (ANC)