Акустический журнал, 2022, T. 68, № 4, стр. 351-360
К дифракции точечного источника звука на бесконечном клине
М. А. Сумбатян a, *, **, Т. С. Мартынова a, Н. К. Мусатова a
a Южный федеральный университет, Институт математики, механики и компьютерных наук им. Воровича
344090 Ростов-на-Дону, ул. Мильчакова 8а, Россия
* E-mail: masumbatyan@sfedu.ru
** E-mail: sumbatma@mail.ru
Поступила в редакцию 09.09.2021
После доработки 13.01.2022
Принята к публикации 25.01.2022
- EDN: KRXVIZ
- DOI: 10.31857/S0320791922030145
Аннотация
Рассматривается двумерная задача дифракции гармонической звуковой волны, исходящей из точечного источника звука, который расположен вблизи острого угла бесконечного клина несимметрично относительно его граней. Граница считается акустически жесткой. В рамках метода граничных интегральных уравнений при несимметричном расположении источника звука задача сводится к системе из двух интегральных уравнений Фредгольма второго рода. Поведение решения при приближении к окрестности угловой точки определяется условием Мейкснера. В явном виде находится значение давления на конце клина в вершине угла. Производится асимптотическая оценка поведения функции давления на бесконечности. Дискретизация сводит систему основных граничных интегральных уравнений к системе линейных алгебраических уравнений. Предлагается “улучшенная” схема дискретизации с тремя интервалами различной плотности на каждой грани. Построено давление в рассеянном поле.
ВВЕДЕНИЕ
Задача дифракции акустических волн на твердом бесконечном клине является классической [1]. По-видимому, Зоммерфельд, создавая свою теорию дифракции, впервые применил к этой задаче строгие математические методы [2]. Геометрическая теория дифракции Келлера [3] также эффективна в применении к рассматриваемой задаче. Что касается численных методов, то их приложение к дифракции на клине изложено в работах [4–6], где также приводится много других полезных ссылок. Общий подробный обзор известных методов дан в работе [7]. В работах [8–10] также рассматриваются задачи дифракции, однако в качестве областей взяты замкнутые контуры в виде тел вращения. В недавней работе [11] авторами рассмотрен частный случай задачи, когда точечный источник звука расположен симметрично относительно граней клина. Главной целью настоящей работы является обобщение результатов работы [11] на случай произвольного несимметричного расположения источника.
Заметим, что эффективные численные методы на основе граничных элементов для дифракции на телах с углами с успехом развивались ранее [12, 13], однако они применимы лишь к выпуклым препятствиям.
МАТЕМАТИЧЕСКАЯ ФОРМУЛИРОВКА ЗАДАЧИ
Рассмотрим двумерную задачу дифракции гармонической звуковой волны, исходящей из точечного источника звука с координатами $\left( {{{x}_{0}},{{y}_{0}}} \right)$, который расположен вблизи острого угла бесконечного клина с углом раствора 2${{\theta }}$, как показано на рис. 1. Считаем границу акустически жесткой (однородное граничное условие Неймана) и состоящей из двух частей – верхней $l{}_{1}$ и нижней $l{}_{2}$.
В рамках метода граничных интегральных уравнений (ГИУ) [14] при несимметричном расположении источника звука задача дифракции сводится к системе из двух интегральных уравнений Фредгольма второго рода:
(1)
$\left\{ \begin{gathered} {{p}_{1}}\left( {y_{0}^{{(1)}}} \right) - 2\int\limits_{{{l}_{2}}} {\frac{{\partial G\left( {{{y}^{{(2)}}},y_{0}^{{(1)}}} \right)}}{{\partial n_{y}^{{(2)}}}}{{p}_{2}}\left( {{{y}^{{(2)}}}} \right)dl_{y}^{{(2)}}} = \hfill \\ = 2p_{1}^{{{\text{inc}}}}\left( {y_{0}^{{(1)}}} \right),\,\,\,\,y_{0}^{{(1)}} \in {{l}_{1}}, \hfill \\ - 2\int\limits_{{{l}_{1}}} {\frac{{\partial G\left( {{{y}^{{(1)}}},y_{0}^{{(2)}}} \right)}}{{\partial n_{y}^{{(1)}}}}{{p}_{1}}\left( {{{y}^{{(1)}}}} \right)dl_{y}^{{(1)}}} + {{p}_{2}}\left( {y_{0}^{{(2)}}} \right) = \hfill \\ = 2p_{2}^{{{\text{inc}}}}\left( {y_{0}^{{(2)}}} \right),\,\,\,y_{0}^{{(2)}} \in {{l}_{2}}, \hfill \\ \end{gathered} \right.$В системе (1) учтено, что в случае, когда точки ${{y}^{{(1)}}}$ и $y_{0}^{{(1)}}$ (или ${{y}^{{(2)}}}$ и $y_{0}^{{(2)}}$) лежат на одной стороне клина, вклад соответствующей части интегрального оператора равен нулю, так как в этом случае
Поведение решения при приближении к окрестности угловой точки определяется условием Мейкснера (см., например, [7]). Оно гласит, что для того, чтобы физическая энергия сохранялась конечной в любой окрестности угла клина, решение должно удовлетворять следующему условию при $r \to 0$:
Выпишем декартовы координаты точек и нормалей, входящих в первое уравнение системы (1). Для этого в первом уравнении направим ось $x$ по верхней грани клина от угла влево, а ось ξ аналогично по нижней грани. Соответственно, во втором уравнении ось $x$ направим по нижней грани клина от угла влево, а ось ξ аналогично по верхней грани. Тогда в системе координат, связанной естественным образом с горизонтальной и вертикальной осью, имеем
Вектор расстояния между точками ${{y}^{{(2)}}}$ и $y_{0}^{{(1)}}$ равен
Для второго уравнения имеем
Кроме того,
Заметим, что модуль расстояния между “внутренней” и “внешней” точками для двух уравнений имеет один и тот же вид:
Помимо этого легко видеть, что один и тот же вид имеют скалярные произведения вектора расстояния на нормаль в обоих уравнениях. Следовательно, производная функции Грина по нормали вычисляется одинаковым образом для обоих уравнений:
(2)
$\begin{gathered} G = \frac{i}{4}H_{0}^{{(1)}}\left( {kr} \right),\,\,\,\frac{{\partial G}}{{\partial {{n}_{y}}}} = \frac{{\partial G}}{{\partial r}}\frac{{\partial r}}{{\partial {{n}_{y}}}}, \\ \frac{{\partial G}}{{\partial r}} = - \frac{{ik}}{4}H_{1}^{{(1)}}\left( {kr} \right), \\ \end{gathered} $(3)
$\begin{gathered} \frac{{\partial r}}{{\partial {{n}_{y}}}} = \frac{{\left( {r{{n}_{y}}} \right)}}{r} = \\ = \frac{{\left( {x - \xi } \right)\cos \theta \sin \theta + \left( {x + \xi } \right)\cos \theta \sin \theta }}{{\sqrt {{{x}^{2}} + {{\xi }^{2}} - 2x\xi \cos \theta } }} = \\ = \frac{{x\sin 2\theta }}{{\sqrt {{{x}^{2}} + {{\xi }^{2}} - 2x\xi \cos \theta } }}, \\ \end{gathered} $(4)
$\frac{{\partial G}}{{\partial {{n}_{y}}}} = - \frac{{ikx\sin 2\theta }}{{4\sqrt {{{x}^{2}} + {{\xi }^{2}} - 2x\xi \cos \theta } }}H_{1}^{{(1)}}\left( {kr} \right).$Собирая вместе формулы (2)–(4), преобразуем систему ГИУ (1) к следующему виду
(5)
$\left\{ \begin{gathered} {{p}_{1}}\left( x \right) + \frac{{ikx\sin 2\theta }}{2}\int\limits_0^\infty {\frac{{H_{1}^{{(1)}}\left( {kr} \right)}}{r}{{p}_{2}}\left( \xi \right)d\xi } = \hfill \\ = 2p_{1}^{{{\text{inc}}}}\left( x \right),\,\,\,x \in \left( {0,\infty } \right), \hfill \\ \frac{{ikx\sin 2\theta }}{2}\int\limits_0^\infty {\frac{{H_{1}^{{(1)}}\left( {kr} \right)}}{r}} {{p}_{1}}\left( \xi \right)d\xi + {{p}_{2}}\left( x \right) = \hfill \\ = 2p_{2}^{{{\text{inc}}}}\left( x \right),\,\,\,\,x \in \left( {0,\infty } \right). \hfill \\ \end{gathered} \right.$При этом давление в падающем поле в первом уравнении равно $p_{{\text{1}}}^{{{\text{inc}}}}\left( x \right) = \frac{i}{4}H_{0}^{{(1)}}\left( {kr_{0}^{{(1)}}} \right)$, а во втором $p_{{\text{2}}}^{{{\text{inc}}}}\left( x \right) = \frac{i}{4}H_{0}^{{(1)}}\left( {kr_{0}^{{(2)}}} \right)$.
ТОЧНОЕ РЕШЕНИЕ В УГЛЕ
В статье авторов [11] найдено точное значение решения в угле в случае источника, расположенного симметрично относительно граней клина. Найдем значение давления на конце клина в явном виде в рассматриваемом несимметричном случае. Для этого к системе (5) применим замену переменной $\xi = tx$. Тогда интегральный член в первом уравнении примет вид
(6)
$\int\limits_0^\infty {\frac{{H_{1}^{{(1)}}\left( {kx\sqrt {1 - 2t\cos 2{{\theta }} + {{t}^{2}}} } \right)}}{{\sqrt {1 - 2t\cos 2{{\theta }} + {{t}^{2}}} }}} {{p}_{2}}\left( {tx} \right)dt.$Устремим $x \to 0$ в (6) и вынесем значение ${{p}_{2}}(0)$ за знак интеграла. Принимаем во внимание асимптотическое разложение функции Ханкеля для малого аргумента [15]
В итоге в (6) возникает интеграл, который вычисляется в явном виде [16]
После подстановки всех полученных выражений система ГИУ (5) приводит к следующему соотношению
Поскольку, как легко видеть, $p_{{\text{1}}}^{{{\text{inc}}}}\left( 0 \right) = p_{{\text{2}}}^{{{\text{inc}}}}\left( 0 \right)$, то точное решение в угле имеет простой вид:
(7)
${{p}_{1}}\left( 0 \right) = {{p}_{2}}\left( 0 \right) = \frac{\pi }{{\pi - \theta }}{{p}^{{{\text{inc}}}}}(0),$ОЦЕНКА РЕШЕНИЯ НА БЕСКОНЕЧНОСТИ
Перейдем к безразмерным переменным $z = kx,\,\,t = k\xi $, в результате система (5) преобразуется к следующему виду
(8)
$\begin{gathered} \left\{ {\begin{array}{*{20}{c}} {{{P}_{1}}\left( z \right) + \frac{{iz\sin 2{{\theta }}}}{2}\int\limits_0^\infty {\frac{{H_{1}^{{(1)}}\left( {r\left( {t,z} \right)} \right)}}{{r\left( {t,z} \right)}}{{P}_{2}}\left( t \right)dt} = 2p_{1}^{{{\text{inc}}}}\left( z \right),} \\ {\frac{{iz\sin 2{{\theta }}}}{2}\int\limits_0^\infty {\frac{{H_{1}^{{(1)}}\left( {r\left( {t,z} \right)} \right)}}{{r\left( {t,z} \right)}}} {{P}_{1}}\left( t \right)dt + {{P}_{2}}\left( z \right) = 2p_{2}^{{{\text{inc}}}}\left( z \right),} \end{array}} \right. \\ r\left( {t,z} \right) = \sqrt {{{z}^{2}} + {{t}^{2}} - 2zt\cos 2{{\theta }}} ,\,\,\,\,P\left( z \right) = p\left( {\frac{z}{k}} \right), \\ p_{{{\text{1,2}}}}^{{{\text{inc}}}}\left( z \right) = \frac{i}{4}H_{0}^{{(1)}}\left( {\sqrt {{{{\left( {k{{x}_{0}}} \right)}}^{2}} + {{{\left( {k{{y}_{0}}} \right)}}^{2}} + {{z}^{2}} + 2z\left( {k{{x}_{0}}\cos {{\theta }} \mp k{{y}_{0}}\sin {{\theta }}} \right)} } \right). \\ \end{gathered} $Важно отметить, что в (8) ядро интегрального оператора в левых частях обоих уравнений системы не зависит от волнового числа k, а единственная зависимость от данного параметра находится в правых частях. При переходе к дискретной форме матрица соответствующей системы линейных алгебраических уравнений (СЛАУ) также не содержит волновое число k.
В дальнейшем все величины, связанные с линейными размерами, будем выражать в безразмерном виде. В частности, далее величины $k{{x}_{0}}$ и $k{{y}_{0}}$ будем обозначать просто ${{x}_{0}}$ и ${{y}_{0}}$ соответственно.
Ключевым моментом является оценка того, на какой интервал конечной длины $0 < z < L$ следует заменить полубесконечный интервал интегрирования, чтобы обеспечить необходимую точность вычислений при численных расчетах. Для прояснения этого момента проведем асимптотическую оценку функции $P\left( z \right)$ при $z \to \infty $. Применим стандартное асимптотическое разложение функций $r$ и ${{r}_{0}}$:
Примем во внимание асимптотическое поведение функции Ханкеля [15]
В силу того, что различия двух уравнений системы локализуются только в правых частях уравнений и не влияют на общий ход рассуждений, объединим эти два случая вместе:
(9)
$\begin{gathered} {{P}_{{1,2}}}\left( z \right) = \frac{i}{2}{{\left( {\frac{2}{{{{\pi }}z}}} \right)}^{{1/2}}}\exp \left[ {i\left( {z + \left( {{{x}_{0}}\cos {{\theta }} \mp {{y}_{0}}\sin {{\theta }}} \right) - {{{\pi }} \mathord{\left/ {\vphantom {{{\pi }} 4}} \right. \kern-0em} 4}} \right)} \right] - \\ - \,\,\frac{{iz\sin 2{{\theta }}}}{2}\int\limits_0^\infty {{{{\left( {\frac{2}{{{{\pi }}z}}} \right)}}^{{1/2}}}\frac{{\exp \left[ {i\left( {z - t\cos 2{{\theta }} - {{3{{\pi }}} \mathord{\left/ {\vphantom {{3{{\pi }}} 4}} \right. \kern-0em} 4}} \right)} \right]}}{z}} {{P}_{{1,2}}}\left( t \right)dt = \\ = \frac{{{{e}^{{iz}}}}}{{\sqrt z }}\left( {\frac{{\exp \left( {{{{{\pi }}i} \mathord{\left/ {\vphantom {{{{\pi }}i} 4}} \right. \kern-0em} 4}} \right)\exp \left[ {i\left( {{{x}_{0}}\cos {{\theta }} \mp {{y}_{0}}\sin {{\theta }}} \right)} \right]}}{{\sqrt {2{{\pi }}} }} - \frac{{\exp \left( {{{ - 3{{\pi }}} \mathord{\left/ {\vphantom {{ - 3{{\pi }}} 4}} \right. \kern-0em} 4}} \right)\sin 2{{\theta }}}}{{\sqrt {2{{\pi }}} }}\int\limits_0^\infty {\exp \left( { - it\cos 2{{\theta }}} \right)} {{P}_{{1,2}}}\left( t \right)dt} \right) \sim {{A}_{ \pm }}\frac{{\exp \left( {iz} \right)}}{{\sqrt z }}, \\ \end{gathered} $где ${{A}_{ \pm }}$ – некоторая константа. Далее можно показать, что с учетом оценки (9) вклад интегралов по отрезку $(L,\infty )$ имеет порядок $O({1 \mathord{\left/ {\vphantom {1 {{{L}^{2}}}}} \right. \kern-0em} {{{L}^{2}}}})$. Таким образом, при применении численных методов следует брать достаточно большую величину $L$. Эта оценка может быть заметно улучшена введением так называемого “нейтрализатора” Ван дер Корпута [11], который дает экспоненциальное убывание отбрасываемого полубесконечного интеграла.
ЧИСЛЕННОЕ РЕШЕНИЕ ОСНОВНОЙ СИСТЕМЫ ГИУ
Для удобства дальнейшей дискретизации перепишем (8) в виде
(10)
$\begin{gathered} \left\{ {\begin{array}{*{20}{c}} {{{P}_{1}}\left( z \right) + \int\limits_0^L {K\left( {t,z} \right)} {{P}_{2}}\left( t \right)dt = 2p_{1}^{{{\text{inc}}}}\left( z \right),\,\,\,\,0 < z < L,} \\ {\int\limits_0^L {K\left( {t,z} \right)} {{P}_{1}}\left( t \right)dt + {{P}_{2}}\left( z \right) = 2p_{2}^{{{\text{inc}}}}\left( z \right),\,\,\,\,0 < z < L,} \end{array}} \right. \\ K\left( {t,z} \right) = \frac{{iz\sin 2{{\theta }}H_{1}^{{(1)}}\left( {r\left( {t,z} \right)} \right)}}{{2r\left( {t,z} \right)}}. \\ \end{gathered} $Разобьем интервал $0 < z < L$ на $N$ малых интервалов равной длины, и выберем в центре каждого из них соответствующий узел сетки:
Теперь интегрирование по каждому малому интервалу проведем, используя квадратурную формулу Симпсона [17]. С этой целью введем на интервале $t \in \left( {{{z}_{j}} - {h \mathord{\left/ {\vphantom {h 2}} \right. \kern-0em} 2},{{z}_{j}} + {h \mathord{\left/ {\vphantom {h 2}} \right. \kern-0em} 2}} \right)$ три точки начала, конца и середины: $t_{1}^{j} = {{z}_{j}} - {h \mathord{\left/ {\vphantom {h 2}} \right. \kern-0em} 2}$, $t_{2}^{j} = {{z}_{j}} + {h \mathord{\left/ {\vphantom {h 2}} \right. \kern-0em} 2}$, $t_{0}^{j} = {{z}_{j}}$. Формула Симпсона дает [17]:
Применяемая дискретизация сводит систему (10) к СЛАУ размером $2N \times 2N$ относительно вектора $P = \left\{ {{{P}_{j}}} \right\}$:
В матричном виде система (11) записывается так:
Однако наилучшего соотношения точности к количеству узлов удалось достичь при следующей схеме дискретизации. Каждая грань клина разбивается на три интервала $(0,{{L}_{1}}),\;({{L}_{1}},{{L}_{2}}),({{L}_{2}},{{L}_{3}})$ с различными, но равномерными шагами ${{N}_{1}},{{N}_{2}},\;{{N}_{3}}$ соответственно. При этом, чем ближе интервал к вершине угла, тем он меньше по длине и тем большая у него плотность узлов, например
ПРИМЕРЫ РАСЧЕТОВ
Заметим, что для четырех частных углов раствора клина задача дифракции имеет точное решение методом мнимых источников [11]. Это дает хорошую возможность для проверки точности используемого численного метода. Одним из таких углов является случай ${{\theta }} = {{3{{\pi }}} \mathord{\left/ {\vphantom {{3{{\pi }}} 4}} \right. \kern-0em} 4}$, для которого полный угол раствора является тупым: $2{{\theta }} = 270^\circ $. В этом случае в дополнение к одному реальному источнику звука ${{X}_{0}}$ в результате множественных переотражений возникает еще 3 мнимых ${{X}_{1}}$,${{X}_{2}}$,${{X}_{3}}$, расположенных симметрично относительно граней клина. В естественной декартовой системе координаты этих точек в безразмерном виде равны ${{X}_{0}} = \left( {{{x}_{0}},{{y}_{0}}} \right)$, ${{X}_{1}} = \left( {{{y}_{0}},{{x}_{0}}} \right)$, ${{X}_{2}} = \left( { - {{x}_{0}}, - {{y}_{0}}} \right)$, ${{X}_{3}} = \left( { - {{y}_{0}}, - {{x}_{0}}} \right)$. Текущая точка на верхней грани имеет координаты $\bar {X} = {{\left( {z,z} \right)} \mathord{\left/ {\vphantom {{\left( {z,z} \right)} {\sqrt 2 }}} \right. \kern-0em} {\sqrt 2 }}$, если считать, что $z$ – это переменная на положительной полупрямой. Тогда расстояния записываются в виде
Если текущая точка располагается на нижней грани $\bar {X} = {{\left( {z, - z} \right)} \mathord{\left/ {\vphantom {{\left( {z, - z} \right)} {\sqrt 2 }}} \right. \kern-0em} {\sqrt 2 }}$, расстояния равны
Таким образом, явное аналитическое решение в этом частном случае имеет вид
Сравнивая численное решение, построенное методом ГИУ, и аналитическое, построенное методом мнимых источников, можно сказать, что как мнимая, так и вещественная части давления визуально совпадают вдоль всего интервала $0 < z < L$ при общем количестве узлов сетки $N = 2100$, что соответствует трем верным значащим цифрам в числовых значениях. Пример для угла ${{\theta }} = {{3{{\pi }}} \mathord{\left/ {\vphantom {{3{{\pi }}} 4}} \right. \kern-0em} 4}$ (рис. 2, 3) показывает, что разбиение граней на 3 равномерных интервала с более густой сеткой в окрестности угла сводит влияние этой окрестности на численное решение к минимуму.
ДАВЛЕНИЕ В РАССЕЯННОМ ВОЛНОВОМ ПОЛЕ
После решения системы ГИУ (5) можно рассчитать давление в рассеянном поле в любой точке акустической среды. Пусть точка $x\left( {{{x}_{1}},{{x}_{2}}} \right)$ – произвольная точка в естественной декартовой системе координат, связанной с горизонтальным и вертикальным направлениями. Рассеянное давление в ней описывается формулой
(12)
$\begin{gathered} {{p}^{{{\text{sc}}}}}\left( x \right) = \int\limits_l {p\left( y \right)\frac{{\partial G\left( {y,x} \right)}}{{\partial {{n}_{y}}}}d{{l}_{y}}} = \\ = \int\limits_{{{l}_{1}}} {{{p}_{1}}({{y}^{{(1)}}})\frac{{\partial G({{y}^{{(1)}}},x)}}{{\partial n_{y}^{{(1)}}}}dl_{y}^{1}} + \int\limits_{{{l}_{2}}} {{{p}_{2}}({{y}^{{(2)}}})\frac{{\partial G({{y}^{{(2)}}},x)}}{{\partial n_{y}^{{(2)}}}}dl_{y}^{2}.} \\ \end{gathered} $Здесь ${{p}_{1}}({{y}^{{(1)}}})$, ${{p}_{2}}({{y}^{{(2)}}})$ – полное акустическое давление на границах ${{l}_{1}}$ и ${{l}_{2}}$, полученное при решении системы ГИУ. Производная функции Грина вычисляется аналогично вычисленной для ядра интегральных операторов выше. Случай ${{l}_{1}}$:
(13)
$\frac{{\partial G\left( {y,x} \right)}}{{\partial n_{y}^{{(1)}}}} = \frac{{\partial G}}{{\partial r_{{xy}}^{{(1)}}}}\frac{{\partial r_{{xy}}^{{(1)}}}}{{\partial n_{y}^{{(1)}}}} = - \frac{i}{4}H_{1}^{{(1)}}\left( {r_{{xy}}^{{(1)}}} \right)\frac{{\left( {{\mathbf{r}}_{{xy}}^{{(1)}} \cdot {\mathbf{n}}_{y}^{{(1)}}} \right)}}{{\left| {{\mathbf{r}}_{{xy}}^{{(1)}}} \right|}}.$Для кривой ${{l}_{2}}$ в формуле (13) верхний индекс $\left( 1 \right)$ заменится на $\left( 2 \right)$.
Координаты точек и вектора нормали для кривой ${{l}_{1}}$:
Аналогичные величины для кривой ${{l}_{2}}$:
В итоге выражение (12) принимает вид:
(14)
$\begin{gathered} {{p}^{{{\text{sc}}}}}\left( x \right) = \frac{i}{4}\int\limits_0^\infty {\left[ {{{p}_{1}}\left( \xi \right)H_{1}^{{(1)}}\left( {r_{{xy}}^{{(1)}}} \right)\frac{{{{x}_{1}}\sin \theta + {{x}_{2}}\cos \theta }}{{r_{{xy}}^{{(1)}}}}} \right.} + \\ \left. { + \,\,{{p}_{2}}\left( \xi \right)H_{1}^{{(1)}}\left( {r_{{xy}}^{{(2)}}} \right)\frac{{{{x}_{1}}\sin \theta - {{x}_{2}}\cos \theta }}{{r_{{xy}}^{{(2)}}}}} \right]d\xi {\text{.}} \\ \end{gathered} $Таким образом, имея численные значения акустического давления на граничной кривой, по формуле (14) находится структура дифрагированной звуковой волны в произвольной точке акустической среды.
ЗАКЛЮЧЕНИЕ
1. Предложенный метод позволяет распространить результаты предыдущей работы авторов [11] на случай произвольного несимметричного расположения источника звука относительно граней клина. При этом как частный случай получается решение симметричной задачи. Оказалось, что такие результаты как точное значение давления в угле, асимптотическая оценка на бесконечности (9) и аналитическое решение задачи для некоторого ограниченного числа углов могут быть также перенесены с симметричного случая на несимметричный.
2. В качестве метода решения выбран полуаналитический метод ГИУ, который позволил получить важные аналитические результаты и, кроме того, построить устойчивую численную схему. Для дискретизации была выбрана сетка, состоящая из трех интервалов, равномерно разбитых на узлы. Самая плотная сетка располагается в небольшой окрестности угла. Такой способ задания сетки доказывает свою эффективность в многочисленных тестах. Квадратурная формула Симпсона при дискретизации свела ГИУ к СЛАУ, для решения которой использовался метод наименьших квадратов LSQR Пейджа и Сондерса. Тестирование численных результатов по методу ГИУ на задаче, имеющей точное решение методом мнимых источников, показало хорошую точность на всем интервале, в том числе в небольшой окрестности угла.
Заметим, что детальный обзор, посвященный численным методам с адаптивными сетками, представлен в [18].
3. После того, как система ГИУ решена предложенным численным методом, волновое поле в произвольной точке акустической среды определяется в виде некоторых интегралов, выражающихся через найденные граничные значения на гранях клина.
Пример расчета амплитуды полного давления для акустического поля приведен на рис. 4, где источник находится сверху вне изображенной области. На рисунке видно, что волна от источника, распространяясь вправо и вниз, где нет препятствия, практически сохраняет свою круговую форму. В то же время, распространяясь влево и вниз, она взаимодействует с верхней гранью клина, существенно искажаясь. При этом хорошо видна зона тени снизу от клина, причем тень становится все заметнее, если двигаться вблизи нижней грани с удалением от угла.
На рис. 5 показана амплитуда рассеянного полного поля от источника, расположенного вблизи верхней грани. Любопытна структура рассеянного волнового поля сверху от клина. В первом приближении волновое поле здесь формируется путем взаимодействия падающих звуковых лучей и лучей, отразившихся от верхней грани. Структура звукового поля при этом довольно сложная. Волновое поле совсем не похоже на то круговое, которое было бы в волне от источника.
Работа выполнена при поддержке Российского фонда фундаментальных исследований (РФФИ), грант № 19-29-06013.
Список литературы
Скучик Е. Основы акустики. Т. 2. М.: Мир, 1976. 542 с.
Sommerfeld A. Mathematical Theory of Diffraction. Boston: Birkhauser, 2004. https://doi.org/10.1007/978-0-8176-8196-8_2
Keller J.B. Geometrical Theory of Diffraction // J. Opt. Soc. Amer. 1962. P. 116–130. https://doi.org/10.1364/JOSA.52.000116
Hacivelioglu F., Sevgi L., Ufimtsev P.Ya. Numerical evaluations of diffraction formulas for the canonical wedge-scattering problem.// IEEE Antennas Propag. Mag. 2013. V. 55. P. 257–272. https://doi.org/10.1109/MAP.2013.6735531
Apaydin G., Sevgi L.A novel wedge diffraction modeling using method of moments (MoM) // Appl. Comput. Electromagn. Soc. J. 2015. V. 30. P. 1053–1059. https://doi.org/10.1109/TAP.2014.2323414
Lu B., Darmon M., Fradkin L., Potel C. Numerical comparison of acoustic wedge models, with application to ultrasonic telemetry // Ultrasonics. 2016. V. 65. P. 5–9. https://doi.org/10.1016/j.ultras.2015.10.009
Nethercote M.A., Assier R.C., Abrahams I.D. Analytical methods for perfect wedge diffraction: A review // Wave Motion. 2020. 102479. https://doi.org/10.1016/j.wavemoti.2019.102479
Корольков А.И., Шанин А.В., Белоус А.А. Дифракция на вытянутом теле вращения с импедансными границами. Метод граничного параболического уравнения // Акуст. журн. 2019. Т. 65. № 4. С. 440–447. https://doi.org/10.1134/S0320791919040063
Кюркчан А.Г., Маненков С.А. Дифракция поля точечного источника на компактном препятствии в непрерывно-слоистом волноводе // Акуст. журн. 2018. Т. 64. № 5. С. 526–533. https://doi.org/10.1134/S0320791918050052
Андронов И.В. Высокочастотная дифракция на узком гиперболоиде вращения // Акуст. журн. 2017. Т. 63. № 2. С. 129–136. https://doi.org/10.7868/S0320791917010014
Sumbatyan M.A., Martynova T.S., Musatova N.K. Boundary element methods in diffraction of a point-source acoustic wave by a rigid infinite wedge // Enging Anal. Boundary Elem. 2021. V. 125. P. 157–167. https://doi.org/10.1016/j.enganabound.2021.01.017
Dominguez V., Graham I.G., Smyshlyaev V.P. A hybrid numerical-asymptotic boundary integral method for high-frequency acoustic scattering // Numer. Math. 2007. V. 106. P. 471–510. https://doi.org/10.1007/s00211-007-0071-4
Arden S., Chandler-Wilde S.N., Longdon S.A. A collocation method for high-frequency scattering by convex polygons // J. Comp. Appl. Math. 2007. V. 204. P. 334–343. https://doi.org/10.1016/j.cam.2006.03.028
Бреббия К., Теллес Ж., Вроубел Л. Методы граничных элементов. М.: Мир, 1987.
Абрамовиц М., Стиган И. Справочник по специальным функциям. М.: Наука, 1979.
Прудников А.П., Брычков Ю.А., Маричев О.И. Интегралы и ряды. Т. 1. М.: ФИЗМАТЛИТ, 2002.
Форсайт Дж., Малькольм М., Моулер К. Машинные методы математических вычислений. М.: Мир, 1989.
Feischl M., Fuhrer Th., Heuer N., Karkulik M., Praetorius D. Adaptive boundary element methods // Archives Comp. Meth. Enging. 2015. V. 22. P. 309–389. https://doi.org/10.1007/s11831-014-9114-z
Дополнительные материалы отсутствуют.
Инструменты
Акустический журнал