Отражение акустических волн от пористой преграды с «открытой» границей

№53-3,

физико-математические науки

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

Похожие материалы

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

Звукопоглощающая способность пористых материалов обусловлена их пористой структурой и наличием большого числа открытых сообщающихся между собой пор, максимальный диаметр которых обычно не превышает 2 мм. Большая удельная поверхность материалов, создаваемая стенками открытых пор, способствует активному преобразованию энергии звуковых колебаний в тепловую энергию вследствие потерь на трение. Упругие свойства скелета материала и наличие воздуха, заключённого в его порах, обусловливают гашение энергии удара и вибрации, что способствует снижению структурного и ударного шума [1, 6, 7].

В данной работе рассматривается отражение акустических волн от пористой преграды с «открытой» границей (рис. 1). Когда граница «открытая», то газ, содержащийся в порах, расположенных непосредственно на границе, соприкасается с флюидом, окружающим пористую преграду.

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

Изображение пористой преграды с «открытыми» порами на границе
Рисунок 1. Изображение пористой преграды с «открытыми» порами на границе

Найдем коэффициенты отражения и прохождения волны. Пусть гармоническая плоская волна с круговой частотой \omega, распространяясь по однородному флюиду, падает перпендикулярно первой границе «флюид – пористая среда». Тогда движение слева от границы (x<0) является суперпозицией двух волн – падающей и отраженной [8, 9], давление p и скорость \upsilon которых определяются по формулам

p_{e}^{(o)} =A_{p,e}^{(o)} \exp \left[i\cdot \left(K_{e} x-\omega \cdot t\right)\right], \upsilon _{e}^{(o)} =A_{\upsilon ,e}^{(o)} \exp \left[i\cdot \left(K_{e} x-\omega \cdot t\right)\right] [1]

p_{e}^{(r)} =A_{p,e}^{(r)} \exp \left[i\cdot \left(-K_{e} x-\omega \cdot t\right)\right], \upsilon _{e}^{(r)} =A_{\upsilon ,e}^{(r)} \exp \left[i\cdot \left(-K_{e} x-\omega \cdot t\right)\right] [2]

Волновые числа «быстрой» Ka и «медленной» Kb волн вычислялись для трех случаев с учетом:

  1. межфазных сил;
  2. межфазного теплообмена;
  3. межфазных сил и межфазного теплообмена [6, 7].

Нижний индекс e соответствует параметрам волны, распространяющейся по флюиду, а верхние индексы o, r – параметрам падающей и отраженной волн. Амплитуды давления A_{p,e}^{(m)} и скорости A_{\upsilon ,e}^{(m)} (m=o, r) падающей и отраженной гармонических волн связаны соотношениями

A_{p,e}^{(o)} =\rho _{e0}^{0} C_{e} A_{\upsilon ,e}^{(o)}, A_{p,e}^{(r)} =-\rho _{e0}^{0} C_{e} A_{\upsilon ,e}^{(r)} [3]

Здесь \rho _{e0}^{0} – плотность флюида, окружающего пористую среду, Ce – скорость распространения волны во флюиде, Ke – волновое число.

При падении волны из флюида на пористую среду в однородной среде возникает отраженная волна [2 - 5], а в пористой среде – прошедшая волна, представляющая собой суперпозицию двух волн – «быстрой» и «медленной». Для определения коэффициентов отражения и прохождения падающей волны запишем граничные условия, количество которых должно соответствовать количеству возникающих волн. На границе $x=0$ должны выполняться следующие условия:

1. непрерывность давления на границе между окружающим преграду флюидом и газом, находящимся в открытых порах на границе:

p_{e} =p_{g} [4]

2. непрерывность нормальной компоненты усредненной по объему скорости на поверхности раздела:

\upsilon _{e} =\alpha _{g0} \upsilon _{g}^{} +\alpha _{s0} \upsilon _{s}^{} [5]

3. равенство сил, действующих на единицу площади поверхности (т.е. суммарные напряжения), по обе стороны от этой поверхности:

p_{e} =-\sigma _{s}^{*} +p_{g} [6]

Здесь p_{e} =p_{e}^{\left(0\right)} +p_{e}^{\left(r\right)}, \upsilon _{e} =\upsilon _{e}^{\left(0\right)} +\upsilon _{e}^{\left(r\right)} – результирующие возмущение давления и скорость частиц на границе со стороны однородной среды соответственно.

Результирующие возмущения эффективного напряжения в скелете \sigma _{s}^{*} и давления в газовой фазе определяются по формулам

\sigma _{s}^{*} =\left[A_{\sigma ,a} \exp \left(iK_{a} x\right)+A_{\sigma ,b} \exp \left(iK_{b} x\right)\right]\exp \left(-i\omega \, t\right) [7]

p_{g} =\left[A_{p,a} \exp \left(iK_{a} x\right)+A_{p,b} \exp \left(iK_{b} x\right)\right]\exp \left(-i\omega \, t\right) [8]

Здесь A_{\sigma a}A_{p,a} и A_{\sigma b}A_{p,b} – соответственно амплитуды напряжения и давления «быстрой» и «медленной» волн, распространяющихся по скелету и газовой фазе. Нижние индексы a и b соответствуют параметрам «быстрой» и «медленной» волн.

Результирующие скорости движения скелета и газа в порах находятся из выражений

\upsilon _{s} =\left[A_{s,a} \exp \left(iK_{a} x\right)+A_{s,b} \exp \left(iK_{b} x\right)\right]\exp \left(-i\omega t\right) [9]

\upsilon _{g} =\left[A_{g,a} \exp \left(iK_{a} x\right)+A_{g,b} \exp \left(iK_{b} x\right)\right]\exp \left(-i\omega t\right) [10]

Зависимости амплитуд возникающих волн от амплитуды падающей волны определены из граничных условий и уравнений движения. После этого найдем коэффициенты отражения N_{1}^{(op)} и прохождения «медленной» Mb и «быстрой» Ma волн на первой границе «флюид – пористая преграда»:

M_{b} =2\left[1-z_{b}/z_{a}+\psi (C_{b}-C_{a}z_{b}/z_{a})\right] [11]

M_{a}=2/\left[1-z_{a}/z_{b}+\psi (C_{a}-C_{b}z_{a}/z_{b})\right] [12]

N_{1}^{\left(op\right)} =M_{b} +M_{a} -1 [13]

Здесь параметры \psi , Cj, zj, B1, \beta и \xi определяются из выражений

\psi =\alpha_{g0}(1+\chi_{T})\rho _{e0}^{0} C_{e}/\rho _{g0}^{0}C_{g}^{2}, C_{j} =\omega /K_{j},

z_{j} =-\chi _{\mu} C_{s}^{2} \left(B_{1}/C_{g}^{2}-1/C_{j}^{2})/\beta\xi\left(1+i\chi_{V}\right), j=a,b,

B_{1} =\left(1+\chi _{T} \right)\left(1+i\chi _{V} \alpha _{s0} \right), \beta =\rho _{g0}^{0} /\rho _{s0}^{0} , \xi =\alpha _{s0}^{0} /\alpha _{g0}^{0}.

Здесь \alpha _{g0} – начальная объемная доля газовой фазы, \alpha _{s0} – начальная объемная доля твердой фазы. Коэффициент \chi _{\mu } учитывает влияние вязкости материала скелета пористой среды на эволюцию «быстрой» и «медленной» волн. Коэффициент \chi _{T} учитывает влияние теплообмена между скелетом и газом, насыщающим поры среды, на динамику «быстрой» и «медленной» волн, распространяющихся в пористой среде. Влияние нестационарных сил межфазного взаимодействия на динамику «быстрой» и «медленной» волн учитывается величиной \chi _{V} .

На рис. 2 представлены результаты расчетов, полученные по формулам [11] – [13]. Параметры фаз взяты при температуре среды 300 {\rm K} и давлении p_{0} =10^{5} Па. Преграда находится в воздухе. Штриховые линии на рис. 2 соответствуют случаю, когда в расчетах учитываются лишь межфазные силы, т.е. получены без учета межфазного теплообмена, а сплошные линии построены с учетом межфазных сил и межфазного теплообмена. Линии 1, 2 и 3 соответствуют следующим случаям: линии 1 – радиус пор a_{0} =10^{-3} м, содержание объемной доли газа \alpha _{g0} =0.9, линии 2 – a_{0} =10^{-3}м, \alpha _{g0} =0.7, линии 3 – a_{0} =10^{-4} м, \alpha _{g0} =0.9.

Расчеты показали, что во всем частотном диапазоне действительная часть коэффициента отражения величина положительная. В области низких частот (\omega <1000 с-1) с уменьшением радиуса пор среды модуль коэффициента отражения увеличивается. В области высоких частот модуль коэффициента отражения слабо зависит от радиуса пор среды. Уменьшение содержания объемной доли газа приводит к увеличению модуля коэффициента отражения во всем частотном диапазоне. Межфазный теплообмен лишь незначительно уменьшает модуль коэффициента отражения.

Зависимости модуля коэффициента отражения на первой границе преграды от круговой частоты. Внешняя однородная среда воздух.
Рисунок 2. Зависимости модуля коэффициента отражения на первой границе преграды от круговой частоты. Внешняя однородная среда воздух.

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

Список литературы

  1. Володин С.В., Дмитриев В.Л., Хусаинов И.Г. Распространение линейных волн во влажных насыщенных газом пористых средах // Теплофизика высоких температур. — 2009. — Т. 47, № 5. — С. 734-740.
  2. Хусаинов И.Г. Воздействие акустическим полем на насыщенную жидкостью пористую среду // Современные проблемы науки и образования. – 2014. – № 6; URL: http://www.science-education.ru/120-15160 (дата обращения: 31.10.2014).
  3. Хусаинов И.Г. Динамика релаксации давления в полости с плоско-параллельными стенками после ее опрессовки // Современные проблемы науки и образования. — 2014. — № 5; URL: http://www.science-education.ru/119-15159 (дата обращения: 31.10.2014).
  4. Хусаинов И.Г. Отражение акустических волн в цилиндрическом канале от перфорированного участка // Прикладная математика и механика. – 2013. – Т. 77, № 3. – С. 441-451.
  5. Хусаинов И.Г. Оценка качества перфорации скважины акустическим методом // Современные проблемы науки и образования. – 2014. – № 5; URL: http://www.science-education.ru/119-14505 (дата обращения: 09.09.2014).
  6. Хусаинов И.Г. Эволюция импульса давления при прохождении через пористую преграду, расположенную в воде // Фундаментальные исследования. – 2014. – № 11–12. – С. 2645-2649.
  7. Хусаинов И.Г., Дмитриев В.Л. Исследование эволюции волнового импульса при прохождении через пористую преграду // Прикладная механика и техническая физика. — 2011. — Т. 52, № 5 (309). — С. 136-145.
  8. Хусаинов И.Г., Хусаинова Г.Я. Исследование параметров пласта методом опрессовки // Современные проблемы науки и образования. – 2014. – № 3; URL: http://www.science-education.ru/117-13813 (дата обращения: 04.07.2014).
  9. Хусаинова Г.Я. Исследование температурных полей при стационарном течении аномальных жидкостей // Автоматизация. Современные технологии. – 2016. – № 7. – С. 13-16.