Акустический журнал, 2021, T. 67, № 2, стр. 126-137
Два подхода к решению задачи дифракции на сфере Януса
Д. В. Крысанов a, А. Г. Кюркчан a, b, c, С. А. Маненков a, *
a Московский технический университет связи и информатики
111024 Москва, ул. Авиамоторная 8а, Россия
b Фрязинский филиал Института радиотехники и электроники им. В.А. Котельникова РАН
141190 Фрязино, Московской обл., пл. Введенского 1, Россия
c ФГУП Центральный научно-исследовательский институт связи
111141 Москва, 1-й проезд Перова поля 8, Россия
* E-mail: mail44471@mail.ru
Поступила в редакцию 14.11.2020
После доработки 12.12.2020
Принята к публикации 22.12.2020
Аннотация
На основе метода продолженных граничных условий разработаны два алгоритма численного решения задачи дифракции плоской волны на сфере Януса в виде проницаемого шара, частично покрытого абсолютно мягким или абсолютно жестким сферическим экраном. Выполнено сравнение результатов расчета интенсивности рассеянного поля, полученных с использованием предложенных методов, с результатами, найденными методом T-матриц. Проведена проверка точности выполнения оптической теоремы и невязки краевого условия в случае условий Дирихле и Неймана на экране. Построены угловые зависимости интенсивности рассеянного поля для различных углов раскрыва отражающего экрана.
ВВЕДЕНИЕ
Частицы, покрытые тонкими экранами (в частности, частицы Януса), представляют большой интерес в антенной технике, медицине и биологии. Например, частицы Януса используются в качестве нано- или микропловцов, управляемых бегущими акустическими волнами, для точного позиционирования лекарств. Их также можно использовать как микрозонды или сенсоры в силу того, что разные части частиц Януса дают различные отклики на внешнее воздействие. Асимметрия частиц Януса позволяет частицам легко перемещаться, что открывает возможности для создания микромоторов [1–3]. Несмотря на прикладную значимость частиц Януса, рассеяние волн на таких структурах исследовано достаточно слабо. В литературе имеется ряд работ, посвященных как акустической, так и электромагнитной задаче дифракции на сфере Януса [3–7]. Для решения данной задачи дифракции использовались, например, метод моментов, основанный на решении переопределенных алгебраических систем [3], метод интегральных уравнений [4], методы, основанные на теории сумматорных уравнений [5] и метод T-матриц [6, 7].
В настоящей работе рассмотрена акустическая задача дифракции плоской волны на проницаемом шаре, частично покрытом бесконечно тонким абсолютно мягким или абсолютно жестким сферическим экраном. Задача рассеяния на покрытом шаре решалась методом продолженных граничных условий (МПГУ), который успешно применялся ранее к решению широкого круга двумерных и трехмерных задач теории дифракции [8–11]. Основная идея метода состоит в том, что граничное условие краевой задачи “переносится” на вспомогательную поверхность, расположенную на небольшом расстоянии от исходной поверхности экрана в области, где ищется поле. Возможность такого подхода обусловлена вещественной аналитичностью дифракционного поля. Таким образом, например, если граничное условие выполняется на поверхности экрана, то оно приближенно выполняется и в некоторой “окрестности” этой поверхности.
Решение рассматриваемой задачи дифракции при помощи МПГУ основано на использовании функции Грина (ФГ) проницаемого шара. При этом исходная краевая задача сводится к интегральному уравнению первого рода относительно поля или его нормальной производной на поверхности экрана с ядром, выражающимся через ФГ шара. Для решения интегрального уравнения использовались две методики. Первая основана на кусочно-постоянной аппроксимации неизвестной функции по двум угловым координатам (в сферической системе координат). С целью алгебраизации задачи использовался метод Крылова–Боголюбова. При этом осевая симметрия задачи не учитывалась. Данный подход применялся ранее авторами для решения близкой задачи дифракции электромагнитного поля на плоском экране, расположенном на границе раздела двух сред [11]. Второй метод основан на использовании осевой симметрии рассматриваемой задачи дифракции. При помощи проектирования двумерного интегрального уравнения на стандартный базис Фурье по азимутальному углу, задача сводилась к решению бесконечного набора одномерных интегральных уравнений относительно зенитного угла ${\theta }$. Для решения указанных одномерных интегральных уравнений применялся базис из полиномов Якоби с соответствующим весом, учитывающим особенность на краю сферического экрана. При этом применялся метод коллокации, аналогичный методу, описанному в работе [12]. При таком подходе удается добиться большей точности, чем при использовании первого метода. Отметим, что данный подход является обобщением двумерной задачи дифракции на цилиндрическом экране, рассмотренной в работе [13].
Основной трудностью при решении рассматриваемой задачи дифракции, как при использовании первой, так и второй методики, является расчет ФГ проницаемого шара, так как ФГ выражается в виде медленно сходящегося ряда, непригодного для численных расчетов. С целью ускорения сходимости указанного ряда выделялась асимптотика n-го члена ряда при $n \to \infty $ (эту величину можно назвать сингулярной частью ФГ). Сингулярная часть ФГ суммировалась аналитически с использованием производящей функции полиномов Лежандра. Оставшаяся (регулярная) часть ФГ представляла собой достаточно быстро сходящийся ряд.
1. ПОСТАНОВКА ЗАДАЧИ
Рассмотрим математическую постановку задачи. Пусть однородный шар радиуса a покрыт бесконечно тонким сферическим экраном $S$ с углом раскрыва ${\text{2}}{{{\theta }}_{{\text{J}}}}$. Введем сферическую систему координат, причем ось z направим вдоль оси рассматриваемого тела вращения (сферы Януса). Геометрия задачи изображена на рис. 1. Предполагаем, что волновые числа и плотности сред вне и внутри шара равны ${{k}_{1}},\;{{{\mu }}_{1}}$ и ${{k}_{2}},\;{{{\mu }}_{2}}$ соответственно. Таким образом, волновое поле вне и внутри шара удовлетворяет однородным уравнениям Гельмгольца:
(1)
$\begin{gathered} \Delta U + k_{1}^{2}U = 0,\,\,\,\,r > a, \\ \Delta U + k_{2}^{2}U = 0,\,\,\,\,{\text{0}} < r < a, \\ \end{gathered} $(3)
${{\left. {\frac{{\partial U}}{{\partial n}}} \right|}_{{r = a,{ \theta } < {{{\theta }}_{{\text{J}}}}}}} = 0.$При ${\theta } \in \left( {{{{\theta }}_{{\text{J}}}},{\pi }} \right)$ выполнены условия сопряжения
(4)
$\left\{ \begin{gathered} \left[ U \right] = 0, \hfill \\ \left[ {\frac{1}{{\mu }}\frac{{\partial U}}{{\partial r}}} \right] = 0, \hfill \\ \end{gathered} \right.$(5)
$\begin{gathered} {{U}_{{{\text{пад}}}}} = \\ = \exp \left( { - i{{k}_{1}}r\left( {\sin {\theta }\sin {{{\theta }}_{{\text{0}}}}\cos {\varphi } + \cos {\theta }\cos {{{\theta }}_{{\text{0}}}}} \right)} \right). \\ \end{gathered} $Здесь ${{{\theta }}_{{\text{0}}}}$ – угол падения плоской волны. Рассеянное поле ${{U}^{1}}$ удовлетворяет условию излучения на бесконечности
(6)
$\mathop {\lim }\limits_{r \to \infty } r\left( {\frac{{\partial {{U}^{1}}}}{{\partial r}} + i{{k}_{1}}{{U}^{1}}} \right) = 0.$Полное поле удовлетворяет также условию Мейкснера на краю сферического экрана.
2. РЕШЕНИЕ ЗАДАЧИ ПРИ ПОМОЩИ МПГУ
Применим МПГУ для решения поставленной задачи дифракции. С этой целью представим поле вне проницаемого шара в виде
(7)
$U({\mathbf{r}}) = {{U}^{0}}({\mathbf{r}}) - \int\limits_S {J({\mathbf{r}}{\kern 1pt} ')G({\mathbf{r}},{\mathbf{r}}{\kern 1pt} ')ds{\kern 1pt} '} $(8)
$U({\mathbf{r}}) = {{U}^{0}}({\mathbf{r}}) + \int\limits_S {J({\mathbf{r}}{\kern 1pt} ')\frac{{\partial G({\mathbf{r}},{\mathbf{r}}{\kern 1pt} ')}}{{\partial r{\kern 1pt} '}}ds{\kern 1pt} '} $(10)
${{G}_{0}} = \frac{{{{e}^{{ - i{{k}_{1}}R}}}}}{{4{\pi }R}},\,\,\,\,R = \left| {{\mathbf{r}} - {\mathbf{r}}{\kern 1pt} '} \right|,$(11)
${{G}_{1}} = \frac{{{{k}_{1}}}}{{4{\pi }i}}\sum\limits_{n = 0}^\infty {(2n + 1){{V}_{n}}h_{n}^{{(2)}}({{k}_{1}}r{\kern 1pt} ')h_{n}^{{(2)}}({{k}_{1}}r){{P}_{n}}(\cos {\gamma })} ,$(12)
$\begin{gathered} \cos \gamma = \sin \theta \sin \theta {\kern 1pt} '\cos \psi + \cos \theta \cos \theta {\kern 1pt} ', \\ \psi = \varphi - \varphi {\kern 1pt} ', \\ \end{gathered} $(13)
${{V}_{n}} = \frac{{{{{\mu }}_{{12}}}{{j}_{n}}({{k}_{1}}a){{{\eta }}_{n}}({{k}_{2}}a) - {{{\eta }}_{n}}({{k}_{1}}a){{j}_{n}}({{k}_{2}}a)}}{{{\xi }_{n}^{{}}({{k}_{1}}a){{j}_{n}}({{k}_{2}}a) - {{{\mu }}_{{12}}}h_{n}^{{(2)}}({{k}_{1}}a){{{\eta }}_{n}}({{k}_{2}}a)}},$(14)
${{{\mu }}_{{12}}} = \frac{{{{{\mu }}_{1}}}}{{{{{\mu }}_{2}}}},\,\,\,\,{{{\eta }}_{n}}(x) = xj_{n}^{'}(x),\,\,\,\,{{{\xi }}_{n}}(x) = xh_{n}^{{(2){\kern 1pt} '}}(x),$Заметим, что первичное поле вне шара имеет вид
(15)
$\begin{gathered} {{U}^{0}}({\mathbf{r}}) = {{U}_{{{\text{пад}}}}}({\mathbf{r}}) + \\ + \,\,\sum\limits_{n = - \infty }^\infty {{{i}^{{ - n}}}(2n + 1){{V}_{n}}h_{n}^{{(2)}}({{k}_{1}}r){{P}_{n}}(\cos {{{\gamma }}_{0}})} , \\ \end{gathered} $В соответствии со стандартной схемой МПГУ подставим далее формулу (7) или (8) в граничное условие (2) или (3), поставленное на вспомогательной поверхности ${{S}_{{\delta }}}$, смещенной на малое расстояние ${\delta }$ от поверхности $S$ [8–11]. В результате задача сведется к решению интегрального уравнения первого рода относительно функции $J({\mathbf{r}}{\kern 1pt} ')$:
(16)
$\int\limits_S {G({\mathbf{r}},{\mathbf{r}}{\kern 1pt} ')J({\mathbf{r}}{\kern 1pt} ')dS{\kern 1pt} '} = {{U}^{0}}({\mathbf{r}}),\,\,\,\,{\mathbf{r}} \in {{S}_{{\delta }}}$(17)
$\int\limits_S {\frac{{{{\partial }^{2}}G\left( {{\mathbf{r}},{\mathbf{r}}{\kern 1pt} '} \right)}}{{\partial r\partial r{\kern 1pt} '}}J({\mathbf{r}}{\kern 1pt} ')dS{\kern 1pt} '} = - \frac{{\partial {{U}^{0}}}}{{\partial r}},\,\,\,\,{\mathbf{r}} \in {{S}_{{\delta }}}$С использованием сферических координат, интегральные уравнения (16) и (17) могут быть записаны в виде
(18)
$\begin{gathered} \int\limits_0^{2{\pi }} {\int\limits_0^{{{{\theta }}_{{\text{J}}}}} {K({\theta },{\varphi },{\theta }{\kern 1pt} {\text{'}},{\varphi }{\kern 1pt} {\text{'}})J({\theta }{\kern 1pt} {\text{'}},{\varphi }{\kern 1pt} {\text{'}})\sin {\theta }{\kern 1pt} {\text{'}}d{\theta }{\kern 1pt} {\text{'}}d{\varphi }{\kern 1pt} {\text{'}}} = B({\theta },{\varphi })} , \\ {\theta } \in \left[ {0,{{{\theta }}_{{\text{J}}}}} \right],\,\,\,\,{\varphi } \in \left[ {0,2{\pi }} \right], \\ \end{gathered} $(19)
$\begin{gathered} K({\theta },{\varphi },{\theta }{\kern 1pt} {\text{'}},{\varphi }{\kern 1pt} {\text{'}}) = {{a}^{2}}{{\left. G \right|}_{{r{\kern 1pt} ' = a,r = a + {\delta }}}}, \\ B({\theta },{\varphi }) = {{\left. {{{U}^{0}}} \right|}_{{r = a + {\delta }}}}, \\ \end{gathered} $(20)
$\begin{gathered} K({\theta },{\varphi },{\theta }{\kern 1pt} {\text{'}},{\varphi }{\kern 1pt} {\text{'}}) = {{a}^{2}}{{\left. {\frac{{{{\partial }^{2}}G}}{{\partial r\partial r{\kern 1pt} '}}} \right|}_{{r{\kern 1pt} ' = a,r = a + {\delta }}}}, \\ B({\theta },{\varphi }) = - {{\left. {\frac{{\partial {{U}^{0}}}}{{\partial r}}} \right|}_{{r = a + {\delta }}}}. \\ \end{gathered} $Как было указано во Введении, ряд ФГ (11) сходится очень медленно в силу того, что отношение ${a \mathord{\left/ {\vphantom {a {(a + {\delta }}}} \right. \kern-0em} {(a + {\delta }}}) \approx 1$ и, следовательно, непригоден для численных расчетов. Поэтому воспользуемся асимптотикой членов ряда при $n \to \infty $. Во-первых, заметим, что
(21)
${{V}_{n}}h_{n}^{{(2)}}({{k}_{1}}a) \sim - \frac{{{{{\mu }}_{{12}}} - 1}}{{{{{\mu }}_{{12}}} + 1}}{{j}_{n}}({{k}_{1}}a),$(22)
${{V}_{n}}{{{\xi }}_{n}}({{k}_{1}}a) \sim \frac{{{{{\mu }}_{{12}}} - 1}}{{{{{\mu }}_{{12}}} + 1}}{{{\eta }}_{n}}({{k}_{1}}a)$(23)
$\begin{gathered} G = \frac{2}{{{{{\mu }}_{{12}}} + 1}}{{G}_{0}} + \frac{{{{k}_{1}}}}{{4{\pi }i}}\sum\limits_{n = 0}^\infty {(2n + 1)h_{n}^{{(2)}}({{k}_{1}}r){{P}_{n}}(\cos {\gamma }) \times } \\ \times \,\,\left( {{{V}_{n}}h_{n}^{{(2)}}({{k}_{1}}a) + \frac{{{{{\mu }}_{{12}}} - 1}}{{{{{\mu }}_{{12}}} + 1}}{{j}_{n}}({{k}_{1}}a)} \right) \equiv \\ \equiv \frac{2}{{{{{\mu }}_{{12}}} + 1}}{{G}_{0}} + \frac{{{{k}_{1}}}}{{4{\pi }i}}\sum\limits_{n = 0}^\infty {{{u}_{n}}{{P}_{n}}(\cos {\gamma }),} \\ \end{gathered} $(24)
$\begin{gathered} \frac{{{{\partial }^{2}}G}}{{\partial r\partial r{\kern 1pt} '}} = \frac{{2{{{\mu }}_{{12}}}}}{{{{{\mu }}_{{12}}} + 1}}\frac{{{{\partial }^{2}}{{G}_{0}}}}{{\partial r\partial r{\kern 1pt} '}} + \frac{1}{{ar}}\frac{{{{k}_{1}}}}{{4{\pi }i}} \times \\ \times \,\,\sum\limits_{n = 0}^\infty {(2n + 1){{{\xi }}_{n}}({{k}_{1}}r){{P}_{n}}(\cos {\gamma }) \times } \\ \times \,\,\left( {{{V}_{n}}{{{\xi }}_{n}}({{k}_{1}}a) - \frac{{{{{\mu }}_{{12}}} - 1}}{{{{{\mu }}_{{12}}} + 1}}{{{\eta }}_{n}}({{k}_{1}}a)} \right) \equiv \\ \equiv \frac{{2{{{\mu }}_{{12}}}}}{{{{{\mu }}_{{12}}} + 1}}\frac{{{{\partial }^{2}}{{G}_{0}}}}{{\partial r\partial r{\kern 1pt} '}} + \frac{1}{{ar}}\frac{{{{k}_{1}}}}{{4{\pi }i}}\sum\limits_{n = 0}^\infty {{{v}_{n}}{{P}_{n}}(\cos {\gamma })} , \\ \end{gathered} $(25)
$\begin{gathered} {{j}_{n}}(x) \sim \frac{{{{x}^{n}}}}{{(2n + 1)!!}}{{{\tilde {j}}}_{n}}(x),\,\,\,\,{{{\eta }}_{n}}(x) \sim \frac{{n{{x}^{n}}}}{{(2n + 1)!!}}{{{{\tilde {\eta }}}}_{n}}(x), \\ h_{n}^{{(2)}}(x) \sim i\frac{{(2n - 1)!!}}{{{{x}^{{n + 1}}}}}\tilde {h}_{n}^{{(2)}}(x), \\ {{{\xi }}_{n}}(x) \sim - i\frac{{n(2n - 1)!!}}{{{{x}^{{n + 1}}}}}{{{{\tilde {\xi }}}}_{n}}(x), \\ \end{gathered} $(26)
$\begin{gathered} {{{\tilde {j}}}_{n}}(x) = 1 - \frac{{{{x}^{2}}}}{{4n}} + \left( {\frac{3}{8}{{x}^{2}} + \frac{{{{x}^{4}}}}{{32}}} \right)\frac{1}{{{{n}^{2}}}} - \\ - \,\,\left( {\frac{9}{{16}}{{x}^{2}} + \frac{{{{x}^{4}}}}{8} + \frac{{{{x}^{6}}}}{{384}}} \right)\frac{1}{{{{n}^{3}}}}, \\ \end{gathered} $(27)
$\begin{gathered} \tilde {h}_{n}^{{(2)}}(x) = 1 + \frac{{{{x}^{2}}}}{{4n}} + \left( {\frac{{{{x}^{2}}}}{8} + \frac{{{{x}^{4}}}}{{32}}} \right)\frac{1}{{{{n}^{2}}}} + \\ + \,\,\left( {\frac{{{{x}^{2}}}}{{16}} + \frac{{{{x}^{4}}}}{{16}} + \frac{{{{x}^{6}}}}{{384}}} \right)\frac{1}{{{{n}^{3}}}}, \\ \end{gathered} $(28)
${{{\tilde {\eta }}}_{n}}(x) = 1 - \frac{{{{x}^{2}}}}{{4n}} + \left( {\frac{{{{x}^{4}}}}{{32}} - \frac{{{{x}^{2}}}}{8}} \right)\frac{1}{{{{n}^{2}}}} + \left( {\frac{3}{{16}}{{x}^{2}} - \frac{{{{x}^{6}}}}{{384}}} \right)\frac{1}{{{{n}^{3}}}},$(29)
$\begin{gathered} {{{{\tilde {\xi }}}}_{n}}(x) = 1 + \left( {1 + \frac{{{{x}^{2}}}}{4}} \right)\frac{1}{n} + \left( {\frac{{{{x}^{4}}}}{{32}} - \frac{{{{x}^{2}}}}{8}} \right)\frac{1}{{{{n}^{2}}}} + \\ + \,\,\left( {\frac{{{{x}^{6}}}}{{384}} - \frac{{{{x}^{4}}}}{{32}} - \frac{{{{x}^{2}}}}{{16}}} \right)\frac{1}{{{{n}^{3}}}}. \\ \end{gathered} $Поясним ускорение сходимости ряда, через который выражается ядро интегрального уравнения (18), на примере условия Дирихле на поверхности экрана. Запишем ФГ (23) в виде
(30)
$\begin{gathered} G = \frac{2}{{{{{\mu }}_{{12}}} + 1}}{{G}_{0}} + \frac{{{{k}_{1}}}}{{4{\pi }i}}\sum\limits_{n = 0}^\infty {{{u}_{n}}{{P}_{n}}(\cos {\gamma })} = \\ = \frac{2}{{{{{\mu }}_{{12}}} + 1}}{{G}_{0}} + \frac{{{{k}_{1}}}}{{4{\pi }i}}\sum\limits_{n = 0}^\infty {\left( {{{u}_{n}} - u_{n}^{\infty }} \right){{P}_{n}}(\cos {\gamma })} + \\ + \,\,\frac{{{{k}_{1}}}}{{4{\pi }i}}\sum\limits_{n = 0}^\infty {u_{n}^{\infty }{{P}_{n}}(\cos {\gamma })} , \\ \end{gathered} $(32)
${\chi }_{n}^{1} = \left( {{{{\tilde {V}}}_{n}}\tilde {h}_{n}^{{(2)}}({{k}_{1}}a) - \frac{{{{{\mu }}_{{12}}} - 1}}{{{{{\mu }}_{{12}}} + 1}}{{{\tilde {j}}}_{n}}({{k}_{1}}a)} \right)\tilde {h}_{n}^{{(2)}}({{k}_{1}}r),$(33)
${{\tilde {V}}_{n}} = \frac{{{{{\mu }}_{{12}}}{{{\tilde {j}}}_{n}}({{k}_{1}}a){{{{\tilde {\eta }}}}_{n}}({{k}_{2}}a) - {{{{\tilde {\eta }}}}_{n}}({{k}_{1}}a){{{\tilde {j}}}_{n}}({{k}_{2}}a)}}{{{{{{\tilde {\xi }}}}_{n}}({{k}_{1}}a){{{\tilde {j}}}_{n}}({{k}_{2}}a) + {{{\mu }}_{{12}}}\tilde {h}_{n}^{{(2)}}({{k}_{1}}a){{{{\tilde {\eta }}}}_{n}}({{k}_{2}}a)}}.$При этом ограничиваемся только одним членом разложения. Теперь с использованием соотношения
(34)
$\sum\limits_{n = 0}^\infty {\frac{{{{\kappa }^{{n + 1}}}{{P}_{n}}(\cos {\gamma })}}{{n + 1}}} = \int\limits_0^\kappa {\frac{{dt}}{{\sqrt {1 - 2\cos {\gamma }t + {{t}^{2}}} }}} ,$(36)
$\tilde {K} = \left\{ {\frac{{2{{a}^{2}}}}{{{{{\mu }}_{{12}}} + 1}}{{G}_{0}} - \frac{{{{a}^{2}}}}{{4{\pi }r}}\frac{{{{A}_{0}}}}{\kappa }\int\limits_0^\kappa {\frac{{dt}}{{\sqrt {1 - 2\cos {\gamma }t + {{t}^{2}}} }}} } \right\},$(37)
$\tilde {\tilde {K}} = - i\frac{{{{k}_{1}}{{a}^{2}}}}{{4{\pi }}}\sum\limits_{n = 0}^\infty {\left( {{{u}_{n}} - u_{n}^{\infty }} \right){{P}_{n}}(\cos {\gamma })} ,$(38)
$\begin{gathered} \int\limits_0^\kappa {\frac{{dt}}{{\sqrt {1 - 2\cos \gamma t + {{t}^{2}}} }}} = \\ = \left[ \begin{gathered} - \ln \left( {1 - \kappa } \right),{\text{ }}\cos \gamma = 1, \hfill \\ \ln \left( {\frac{{\kappa - \cos \gamma + \sqrt {{{\kappa }^{2}} - 2\cos \gamma \kappa + 1} }}{{1 - \cos \gamma }}} \right),{\text{ }}\cos \gamma \ne 1. \hfill \\ \end{gathered} \right. \\ \end{gathered} $Аналогично, в случае абсолютно жесткого экрана c использованием равенства
(39)
$\sum\limits_{n = 0}^\infty {n{{\kappa }^{{n - 1}}}{{P}_{n}}(\cos \gamma )} = \frac{{\cos \gamma - \kappa }}{{{{{\left( {1 - 2\cos \gamma \kappa + {{\kappa }^{2}}} \right)}}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}}}},$(40)
$\begin{gathered} \tilde {K} = \left\{ {\frac{{2{{\mu }_{{12}}}{{a}^{2}}}}{{{{\mu }_{{12}}} + 1}}\frac{{{{\partial }^{2}}{{G}_{0}}}}{{\partial r\partial r{\kern 1pt} '}} - \frac{a}{{4\pi {{r}^{2}}}} \times } \right. \\ \times \,\,\left[ {\frac{{{{A}_{1}}\kappa (\cos \gamma - \kappa )}}{{{{{\left( {1 - 2\cos \gamma \kappa + {{\kappa }^{2}}} \right)}}^{{{3 \mathord{\left/ {\vphantom {3 2}} \right. \kern-0em} 2}}}}}} + \frac{{{{A}_{2}}}}{{{{{\left( {1 - 2\cos \gamma \kappa + {{\kappa }^{2}}} \right)}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}}} + } \right. \\ \left. {\left. { + \,\,\frac{{{{A}_{3}}}}{\kappa }\int\limits_0^\kappa {\frac{{dt}}{{\sqrt {1 - 2\cos \gamma t + {{t}^{2}}} }}} } \right]} \right\}, \\ \end{gathered} $(41)
$\tilde {\tilde {K}} = - i\frac{{{{k}_{1}}a}}{{4{\pi }r}}\sum\limits_{n = 0}^\infty {\left( {{{v}_{n}} - v_{n}^{\infty }} \right){{P}_{n}}(\cos {\gamma })} .$Коэффициенты ${{A}_{1}},{\text{ }}{{A}_{2}},{\text{ }}{{A}_{3}}$ находим из асимптотики следующей величины (при $n \to \infty $):
(42)
${\chi }_{n}^{2} = {{n}^{2}}{{{\tilde {\xi }}}_{n}}({{k}_{1}}r)\left[ {{{{\tilde {V}}}_{n}}{{{{\tilde {\xi }}}}_{n}}({{k}_{1}}a) - \frac{{{{{\mu }}_{{12}}} - 1}}{{{{{\mu }}_{{12}}} + 1}}{{{{\tilde {\eta }}}}_{n}}({{k}_{1}}a)} \right].$При этом сохраняем только три члена разложения. В силу громоздкости соответствующих выкладок, для нахождения коэффициентов ${{A}_{0}},{{A}_{1}},{\text{ }}{{A}_{2}},{\text{ }}{{A}_{3}}$ использовались символьные вычисления в среде Maple. Отметим, что члены рядов в формулах (37) и (41) имеют порядок не больше $\frac{1}{{{{n}^{{{5 \mathord{\left/ {\vphantom {5 2}} \right. \kern-0em} 2}}}}}}$, то есть сходятся достаточно быстро.
3. ЧИСЛЕННОЕ РЕШЕНИЕ ДВУМЕРНОГО ИНТЕГРАЛЬНОГО УРАВНЕНИЯ НА ОСНОВЕ ДВУХ ПОДХОДОВ
Рассмотрим вопрос о численном решении уравнений (18) при помощи двух методов. Первый подход основан на использовании кусочно-постоянной аппроксимации неизвестной функции с последующим применением метода коллокации. А именно, вводим сетку:
(43)
$\begin{gathered} {{{\theta }}_{q}} = \frac{{{{{\theta }}_{{\text{J}}}}}}{{{{N}_{1}}}}\left( {q - \frac{1}{2}} \right),\,\,\,\,q = 1,2,...,{{N}_{1}}, \\ {{{\varphi }}_{l}} = \frac{{2{\pi }}}{{{{N}_{2}}}}\left( {l - \frac{1}{2}} \right),\,\,\,\,l = 1,2,...,{{N}_{2}}. \\ \end{gathered} $Неизвестную функцию $J(\theta {\kern 1pt} ',\varphi {\kern 1pt} ')$ ищем в виде:
где ${{f}_{q}}({\theta }{\kern 1pt} {\text{'}})$ и ${{w}_{l}}({\varphi }{\kern 1pt} {\text{'}})$ – импульсные функции:(45)
$\begin{gathered} {{f}_{q}}({\theta }{\kern 1pt} {\text{'}}) = \left[ \begin{gathered} 1,{\text{ }}{{{\theta }}_{q}} - {{{{\Delta }_{1}}} \mathord{\left/ {\vphantom {{{{\Delta }_{1}}} 2}} \right. \kern-0em} 2} < {\theta }{\kern 1pt} {\text{'}} < {{{\theta }}_{q}} + {{{{\Delta }_{1}}} \mathord{\left/ {\vphantom {{{{\Delta }_{1}}} 2}} \right. \kern-0em} 2}, \hfill \\ 0,{\text{ }}else, \hfill \\ \end{gathered} \right. \hfill \\ {{w}_{l}}({\varphi }{\kern 1pt} {\text{'}}) = \left[ \begin{gathered} 1,{\text{ }}{{{\varphi }}_{l}} - {{{{\Delta }_{2}}} \mathord{\left/ {\vphantom {{{{\Delta }_{2}}} 2}} \right. \kern-0em} 2} < {\varphi }{\kern 1pt} {\text{'}} < {{{\varphi }}_{l}} + {{{{\Delta }_{2}}} \mathord{\left/ {\vphantom {{{{\Delta }_{2}}} 2}} \right. \kern-0em} 2}, \hfill \\ 0,{\text{ }}else, \hfill \\ \end{gathered} \right. \hfill \\ \end{gathered} $(46)
$\begin{gathered} \sum\limits_{q = 1}^{{{N}_{1}}} {\sum\limits_{l = 1}^{{{N}_{2}}} {\left( {{{{\tilde {K}}}_{{pj,ql}}} + {{{\tilde {\tilde {K}}}}_{{pj,ql}}}} \right){{c}_{{ql}}}} } = {{B}_{{pj}}}, \\ p = \overline {1,{{N}_{1}}} ,\,\,\,\,j = \overline {1,{{N}_{2}}} , \\ \end{gathered} $(47)
${{\tilde {K}}_{{pj,ql}}} = \int\limits_{{{{\varphi }}_{l}} - {{{{\Delta }_{2}}} \mathord{\left/ {\vphantom {{{{\Delta }_{2}}} 2}} \right. \kern-0em} 2}}^{{{{\varphi }}_{l}} + {{{{\Delta }_{2}}} \mathord{\left/ {\vphantom {{{{\Delta }_{2}}} 2}} \right. \kern-0em} 2}} {\int\limits_{{{{\theta }}_{q}} - {{{{\Delta }_{1}}} \mathord{\left/ {\vphantom {{{{\Delta }_{1}}} 2}} \right. \kern-0em} 2}}^{{{{\theta }}_{q}} + {{{{\Delta }_{1}}} \mathord{\left/ {\vphantom {{{{\Delta }_{1}}} 2}} \right. \kern-0em} 2}} {\tilde {K}({{{\theta }}_{p}},{{{\varphi }}_{j}},{\theta }{\kern 1pt} {\text{'}},{\varphi }{\kern 1pt} {\text{'}})d{\theta }{\kern 1pt} {\text{'}}d{\varphi }{\kern 1pt} {\text{'}}} } ,$В последнем интеграле применяем приближенную формулу, являющуюся двумерным аналогом формулы прямоугольников:
(49)
${{\tilde {\tilde {K}}}_{{pj,ql}}} = \tilde {\tilde {K}}({{\theta }_{p}},{{\varphi }_{j}},{{\theta }_{q}},{{\varphi }_{l}}){{\Delta }_{1}}{{\Delta }_{2}},$(50)
$\begin{gathered} {{B}_{{pj}}} = \sum\limits_{n = 0}^\infty {{{i}^{{ - n}}}(2n + 1){{P}_{n}}\left( {\cos \left( {{{{\gamma }}_{0}}({{{\theta }}_{p}},{{{\varphi }}_{j}})} \right)} \right)} \times \\ \times \,\,\left( {{{j}_{n}}({{k}_{1}}r) + {{V}_{n}}h_{n}^{{(2)}}({{k}_{1}}r)} \right) \\ \end{gathered} $(51)
$\begin{gathered} {{B}_{{pj}}} = - {{k}_{1}}\sum\limits_{n = 0}^\infty {{{i}^{{ - n}}}(2n + 1){{P}_{n}}\left( {\cos \left( {{{{\gamma }}_{0}}({{{\theta }}_{p}},{{{\varphi }}_{j}})} \right)} \right)} \times \\ \times \,\,\left( {j_{n}^{'}({{k}_{1}}r) + {{V}_{n}}h_{n}^{{(2)'}}({{k}_{1}}r)} \right) \\ \end{gathered} $В качестве результата решения задачи дифракции мы рассматривали диаграмму рассеяния частицы Януса, которая в данном случае определяется по формулам:
(52)
${{U}^{1}} \approx g(\theta ,\varphi )\frac{{\exp ( - i{{k}_{1}}r)}}{{{{k}_{1}}r}},\,\,\,\,r \to \infty .$Здесь
(54)
${{g}^{0}}(\theta ,\varphi ) = i\sum\limits_{n = 0}^\infty {(2n + 1){{V}_{n}}} {{P}_{n}}\left( {\cos \left( {{{\gamma }_{0}}(\theta ,\varphi )} \right)} \right)$(55)
$\begin{gathered} {{g}^{1}}(\theta ,\varphi ) = - \frac{1}{{4\pi }}\sum\limits_{q = 1}^{{{N}_{1}}} {\sum\limits_{l = 1}^{{{N}_{2}}} {{{c}_{{ql}}}\left[ {\exp (i{{k}_{1}}a\cos {{\gamma }_{{ql}}})} \right.} } + \\ + \,\,\left. {\sum\limits_{n = 0}^\infty {{{i}^{n}}(2n + 1){{V}_{n}}h_{n}^{{(2)}}({{k}_{1}}a){{P}_{n}}(\cos {{\gamma }_{{ql}}})} } \right]{{a}^{2}}\sin {{\theta }_{q}}{{\Delta }_{1}}{{\Delta }_{2}} \\ \end{gathered} $(56)
$\begin{gathered} {{g}^{1}}(\theta ,\varphi ) = \frac{1}{{4\pi }}\sum\limits_{q = 1}^{{{N}_{1}}} {\sum\limits_{l = 1}^{{{N}_{2}}} {{{c}_{{ql}}}\left[ {i\cos {{\gamma }_{q}}\exp (i{{k}_{1}}a\cos {{\gamma }_{{lq}}})} \right.} } + \\ + \,\,\left. {\sum\limits_{n = 0}^\infty {{{i}^{n}}(2n + 1){{V}_{n}}h_{n}^{{(2)'}}({{k}_{1}}a){{P}_{n}}(\cos {{\gamma }_{{lq}}})} } \right] \times \\ \times \,\,{{\left( {{{k}_{1}}a} \right)}^{2}}\sin {{\theta }_{q}}{{\Delta }_{1}}{{\Delta }_{2}} \\ \end{gathered} $(57)
$\cos {{\gamma }_{{ql}}} = \sin {{\theta }_{q}}\sin \theta \cos \left( {\varphi - {{\varphi }_{l}}} \right) + \cos {{\theta }_{q}}\cos \theta .$Опишем основные моменты решения интегрального уравнения (18) при помощи второй методики. Учтем осевую симметрию частицы Януса. Разложим неизвестную функцию в ряд Фурье
(58)
$J(\theta {\kern 1pt} ',\varphi {\kern 1pt} ') = \sum\limits_{m = - \infty }^\infty {{{I}_{m}}(\theta {\kern 1pt} ')\exp (im\varphi {\kern 1pt} ').} $В результате подстановки формулы (58) в интегральное уравнение (18) и проектирования обеих частей равенства на базис Фурье, получим следующую бесконечную систему одномерных интегральных уравнений
(59)
$\begin{gathered} \int\limits_0^{{{{\theta }}_{{\text{J}}}}} {{{K}_{m}}} ({\theta },{\theta }{\kern 1pt} {\text{'}}){{I}_{m}}({\theta }{\kern 1pt} {\text{'}})\sin {\theta }{\kern 1pt} {\text{'}}d{\theta }{\kern 1pt} {\text{'}} = {{B}_{m}}({\theta }), \\ m = 0, \pm 1, \pm 2,...,\,\,\,\,{\theta } \in [0,{{{\theta }}_{{\text{J}}}}], \\ \end{gathered} $(60)
${{\tilde {K}}_{m}} = \frac{{{{k}_{1}}{{a}^{2}}}}{{{{\mu }_{{12}}} + 1}}{{S}_{m}} - \frac{a}{2}{{A}_{0}}{{F}_{m}},$(61)
$\begin{gathered} {{{\tilde {\tilde {K}}}}_{m}} = - i\frac{{{{k}_{1}}{{a}^{2}}}}{2}\sum\limits_{n = |m|}^\infty {\left( {{{u}_{n}} - u_{n}^{\infty }} \right)\frac{{(n - m)!}}{{(n + m)!}}} \times \\ \times \,\,P_{n}^{m}(\cos \theta )P_{n}^{m}(\cos \theta {\kern 1pt} '), \\ \end{gathered} $(62)
$\begin{gathered} {{B}_{m}}({\theta }) = \sum\limits_{n = |m|}^\infty {{{i}^{{ - n}}}(2n + 1)\frac{{(n - m)!}}{{(n + m)!}}P_{n}^{m}(\cos {\theta })} \times \\ \times \,\,P_{n}^{m}(\cos {{{\theta }}_{0}})\left( {{{j}_{n}}({{k}_{1}}r) + {{V}_{n}}h_{n}^{{(2)}}({{k}_{1}}r)} \right) \\ \end{gathered} $в случае абсолютно мягкого экрана и
(63)
$\begin{gathered} {{{\tilde {K}}}_{m}} = \frac{{{{k}_{1}}\kappa {{\mu }_{{12}}}}}{{{{\mu }_{{12}}} + 1}} \times \\ \times \,\,\left( {\frac{{{{\partial }^{2}}{{S}_{m}}}}{{\partial \rho \partial \rho {\kern 1pt} '}}\rho \rho {\kern 1pt} ' + \frac{{{{\partial }^{2}}{{S}_{m}}}}{{\partial z\partial \rho {\kern 1pt} '}}z\rho {\kern 1pt} ' + \frac{{{{\partial }^{2}}{{S}_{m}}}}{{\partial \rho \partial z{\kern 1pt} '}}\rho z{\kern 1pt} '\,\, + \frac{{{{\partial }^{2}}{{S}_{m}}}}{{\partial z\partial z{\kern 1pt} '}}zz{\kern 1pt} '} \right) - \\ - \,\,\frac{\kappa }{2}\left( {{{A}_{1}}} \right.\left( {\frac{{\tilde {S}_{{m + 1}}^{I} + \tilde {S}_{{m - 1}}^{I}}}{2}\rho \rho {\kern 1pt} ' + \tilde {S}_{m}^{I}(zz{\kern 1pt} ' - {{a}^{2}})} \right) + \\ \left. { + \,\,{{A}_{2}}{{{\tilde {S}}}_{m}} + \frac{1}{a}{{A}_{3}}{{F}_{m}}} \right), \\ \end{gathered} $(64)
$\begin{gathered} {{{\tilde {\tilde {K}}}}_{m}} = - i\frac{{{{k}_{1}}\kappa }}{2}\sum\limits_{n = |m|}^\infty {\left( {{{v}_{n}} - v_{n}^{\infty }} \right)\frac{{(n - m)!}}{{(n + m)!}}} \times \\ \times \,\,P_{n}^{m}(\cos \theta )P_{n}^{m}(\cos \theta {\kern 1pt} '), \\ \end{gathered} $(65)
$\begin{gathered} {{B}_{m}}({\theta }) = - {{k}_{1}}\sum\limits_{n = |m|}^\infty {{{i}^{{ - n}}}(2n + 1)\frac{{(n - m)!}}{{(n + m)!}}} \times \\ \times \,\,P_{n}^{m}(\cos {\theta })P_{n}^{m}(\cos {{{\theta }}_{0}})\left( {j_{n}^{'}({{k}_{1}}r) + {{V}_{n}}h_{n}^{{(2)'}}({{k}_{1}}r)} \right) \\ \end{gathered} $в случае абсолютно жесткого. Здесь
(66)
$\begin{gathered} {{S}_{m}} = \frac{1}{{2\pi }}\int\limits_0^{2{\pi }} {\frac{{\exp ( - i{{k}_{1}}R - im{\psi })}}{{{{k}_{1}}R}}} d{\psi }, \\ {{{\tilde {S}}}_{m}} = \frac{1}{{2\pi }}\int\limits_0^{2{\pi }} {\frac{{\exp ( - im{\psi })}}{{{{k}_{1}}R}}} d{\psi }, \\ \end{gathered} $(67)
$\begin{gathered} \tilde {S}_{m}^{\operatorname{I} } = \frac{1}{{2\pi }}\int\limits_0^{2{\pi }} {\frac{{\exp ( - im{\psi })}}{{{{{({{k}_{1}}R)}}^{3}}}}} d{\psi }, \\ {{F}_{m}} = \frac{1}{{2\pi }}\int\limits_0^{2{\pi }} {F({\psi })\exp ( - im{\psi })} d{\psi }, \\ \end{gathered} $Для решения одномерных интегральных уравнений (59) воспользуемся разложением неизвестных функций ${{I}_{m}}({\theta }{\kern 1pt} {\text{'}})$ в ряд по ортогональным полиномам Якоби с весом, учитывающим особенность на краю экрана. Предварительно сделаем замену переменной $x = \cos {\theta }$. В результате
(68)
${{I}_{m}}(x) = \frac{1}{{\sqrt {x - {{x}_{0}}} }}\sum\limits_{q = 0}^{L - 1} {a_{q}^{m}P_{q}^{{(0,{{ - 1} \mathord{\left/ {\vphantom {{ - 1} 2}} \right. \kern-0em} 2})}}\left( {2\frac{{x - {{x}_{0}}}}{{1 - {{x}_{0}}}} - 1} \right)} $(69)
${{I}_{m}}(x) = \sqrt {x - {{x}_{0}}} \sum\limits_{q = 0}^{L - 1} {a_{q}^{m}P_{q}^{{(0,{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2})}}\left( {2\frac{{x - {{x}_{0}}}}{{1 - {{x}_{0}}}} - 1} \right)} $(70)
${{x}_{p}} = {{x}_{0}} + \frac{1}{2}({{t}_{p}} + 1)(1 - {{x}_{0}}),\,\,\,\,p = \overline {1,L} ,$(71)
$P_{L}^{{(0,{{ - 1} \mathord{\left/ {\vphantom {{ - 1} 2}} \right. \kern-0em} 2})}}({{t}_{p}}) = 0,\,\,{\text{либо}}\,\,P_{L}^{{(0,{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2})}}({{t}_{p}}) = 0,$(72)
$\sum\limits_{q = 0}^{L - 1} {\left( {\tilde {G}_{{pq}}^{m} + \tilde {\tilde {G}}_{{pq}}^{m}} \right)a_{q}^{m}} = B_{p}^{m},\,\,\,\,p = \overline {1,L} ,$(73)
$\tilde {G}_{{pq}}^{m} = 2\sqrt {1 - {{x}_{0}}} \int\limits_0^1 {{{{\tilde {K}}}_{m}}({{x}_{p}},{{x}_{0}} + {{{\tau }}^{2}}(1 - {{x}_{0}})){{P}_{{2q}}}({\tau })d{\tau }} ,$(74)
$\begin{gathered} \tilde {\tilde {G}}_{{pq}}^{m} = - i\frac{{{{k}_{1}}{{a}^{2}}}}{2}\sum\limits_{n = |m|}^\infty {\left( {{{u}_{n}} - u_{n}^{\infty }} \right)\frac{{(n - m)!}}{{(n + m)!}}P_{n}^{m}({{x}_{p}}} ) \times \\ \times \,\,\int\limits_{{{x}_{0}}}^1 {P_{n}^{m}(x)} {{P}_{{2q}}}\left( {\tau } \right)\frac{{dx}}{{\sqrt {x - {{x}_{0}}} }}, \\ \end{gathered} $(75)
$\begin{gathered} B_{p}^{m} = \sum\limits_{n = |m|}^\infty {{{i}^{{ - n}}}(2n + 1)\frac{{(n - m)!}}{{(n + m)!}}P_{n}^{m}({{x}_{p}})P_{n}^{m}} (\cos {{{\theta }}_{0}}) \times \\ \times \,\,\left( {{{j}_{n}}({{k}_{1}}r) + {{V}_{n}}h_{n}^{{(2)}}({{k}_{1}}r)} \right) \\ \end{gathered} $(76)
$\begin{gathered} \tilde {G}_{{pq}}^{m} = 2\sqrt {{{{(1 - {{x}_{0}})}}^{3}}} \times \\ \times \,\,\int\limits_0^1 {{{{\tilde {K}}}_{m}}({{x}_{p}},{{x}_{0}} + {{{\tau }}^{2}}(1 - {{x}_{0}})){{P}_{{2q + 1}}}({\tau }){\tau }d{\tau }} , \\ \end{gathered} $(77)
$\begin{gathered} \tilde {\tilde {G}}_{{pq}}^{m} = - i\frac{{{{k}_{1}}\kappa }}{2}\sum\limits_{n = |m|}^\infty {\left( {{{v}_{n}} - v_{n}^{\infty }} \right)\frac{{(n - m)!}}{{(n + m)!}}P_{n}^{m}({{x}_{p}}} ) \times \\ \times \,\,\int\limits_{{{x}_{0}}}^1 {P_{n}^{m}(x)} \frac{{{{P}_{{2q + 1}}}\left( {\tau } \right)}}{{\tau }}\sqrt {x - {{x}_{0}}} dx, \\ \end{gathered} $(78)
$\begin{gathered} B_{p}^{m} = - {{k}_{1}}\sum\limits_{n = |m|}^\infty {{{i}^{{ - n}}}(2n + 1)\frac{{(n - m)!}}{{(n + m)!}}P_{n}^{m}({{x}_{p}})P_{n}^{m}} (\cos {{{\theta }}_{0}}) \times \\ \times \,\,\left( {j_{n}^{'}({{k}_{1}}r) + {{V}_{n}}h_{n}^{{(2)'}}({{k}_{1}}r)} \right) \\ \end{gathered} $Приведем формулы для диаграммы рассеяния при использовании второй методики. Диаграмма вновь имеет вид (53). При этом
(79)
${{g}^{1}}({\theta },{\varphi }) = \sum\limits_{m = - \infty }^\infty {(g_{m}^{{10}}({\theta }) + g_{m}^{{11}}({\theta }))\exp (im{\varphi })} ,$(80)
$\begin{gathered} g_{m}^{{10}}({\theta }) = - \frac{{{{{\left( {{{k}_{1}}a} \right)}}^{2}}}}{2}{{i}^{m}}\sum\limits_{q = 0}^{L - 1} {a_{q}^{m}} \times \\ \times \,\,\int\limits_{{{x}_{0}}}^1 {{{J}_{m}}\left( {{{k}_{1}}a\sin {\theta }\sqrt {1 - {{x}^{2}}} } \right)\exp (i{{k}_{1}}a\cos {\theta }) \times } \\ \times \,\,{{P}_{{2q}}}\left( {\tau } \right)\frac{{dx}}{{\sqrt {x - {{x}_{0}}} }}, \\ \end{gathered} $(81)
$\begin{gathered} g_{m}^{{11}}({\theta }) = - \frac{{{{k}_{1}}{{a}^{2}}}}{2}\sum\limits_{q = 0}^{L - 1} {a_{q}^{m}} \times \\ \times \,\,\sum\limits_{n = |m|}^\infty {{{i}^{n}}(2n + 1)\frac{{(n - m)!}}{{(n + m)!}}P_{n}^{m}(\cos {\theta }){{V}_{n}}h_{n}^{{(2)}}({{k}_{1}}a) \times } \\ \times \,\,\int\limits_{{{x}_{0}}}^1 {P_{n}^{m}(x){{P}_{{2q}}}\left( {\tau } \right)\frac{{dx}}{{\sqrt {x - {{x}_{0}}} }}} \\ \end{gathered} $(82)
$\begin{gathered} g_{m}^{{10}}({\theta }) = \frac{{{{{\left( {{{k}_{1}}a} \right)}}^{2}}}}{2}{{i}^{m}}\sum\limits_{q = 0}^{L - 1} {a_{q}^{m}} \times \\ \times \,\,\int\limits_{{{x}_{0}}}^1 {\left[ {\sin {\theta }\sqrt {1 - {{x}^{2}}} J_{m}^{'}\left( {{{k}_{1}}a\sin {\theta }\sqrt {1 - {{x}^{2}}} } \right) + } \right.} \\ \left. { + \,\,ix\cos {\theta }{{J}_{m}}\left( {{{k}_{1}}a\sin {\theta }\sqrt {1 - {{x}^{2}}} } \right)} \right] \times \\ \times \,\,\exp (i{{k}_{1}}a\cos {\theta })\frac{{{{P}_{{2q + 1}}}\left( {\tau } \right)}}{{\tau }}\sqrt {x - {{x}_{0}}} dx, \\ \end{gathered} $(83)
$\begin{gathered} g_{m}^{{11}}({\theta }) = \frac{{{{k}_{1}}a}}{2}\sum\limits_{q = 0}^{L - 1} {a_{q}^{m}} \times \\ \times \,\,\sum\limits_{n = |m|}^\infty {{{i}^{n}}(2n + 1)\frac{{(n - m)!}}{{(n + m)!}}P_{n}^{m}(\cos {\theta }){{V}_{n}}{{{\xi }}_{n}}({{k}_{1}}a) \times } \\ \times \,\,\int\limits_{{{x}_{0}}}^1 {P_{n}^{m}(x)\frac{{{{P}_{{2q + 1}}}\left( {\tau } \right)}}{{\tau }}\sqrt {x - {{x}_{0}}} dx} \\ \end{gathered} $4. ЧИСЛЕННЫЕ РЕЗУЛЬТАТЫ
Рассмотрим результаты численного моделирования. Для тестирования разработанных методов было проведено сравнение с результатами, полученными при помощи метода Т-матриц в работе [7]. Волновые числа и плотности сред вне и внутри шара были равны ${{k}_{1}} = 1,\;{{{\mu }}_{1}} = 1$ и ${{k}_{2}} = 1.5,\;{{{\mu }}_{2}} = 1.5$, соответственно. Параметр ${\delta }$ при использовании МПГУ во всех случаях брался равным ${{10}^{{ - 3}}}$. Угол полураскрыва экрана был равен ${{{\theta }}_{{\text{J}}}} = 90^\circ $. Угол падения первичной волны ${{{\theta }}_{{\text{0}}}} = 0^\circ $. Рассматривались два различных радиуса шара ${{k}_{1}}a = 3$ и ${{k}_{1}}a = 6$. Число точек коллокации для первого метода ${{N}_{1}} = 40,{{N}_{2}} = 160$, для второго метода $L = 32$. В качестве рассматриваемой характеристики бралась интенсивность рассеянного поля, определяемая следующим образом:
На рис. 2 и 3 приведены угловые зависимости интенсивности рассеяния для метода Т-матриц (кривая 1) и предлагаемых подходов на основе МПГУ (кривая 2 соответствует первому, а кривая 3 – второму методу). На рис. 2 видно, что для условия Дирихле на экране результаты, полученные всеми тремя методами, совпадают с графической точностью. Для условия Неймана (рис. 3) результаты расчета интенсивности рассеянного поля, которые получены при помощи подхода, основанного на использовании кусочно-постоянной аппроксимации (кривая 2), отличаются меньше от соответствующих результатов, полученных методом Т-матриц, чем при использовании второго метода. Данный факт можно объяснить тем, что базис (69) учитывает особенность неизвестной функции на краю экрана, в то время как в методе Т-матриц и в первом методе на основе МПГУ эта особенность не учитывается. Отметим также, что в силу наличия второй нормальной производной ФГ в случае условия Неймана на экране, точность расчета диаграммы при помощи МПГУ несколько ниже, чем в случае условия Дирихле.
В качестве критерия проверки правильности полученных результатов рассматривалась оптическая теорема, которая записывается в виде [16]
где(86)
$\sigma = \frac{1}{{4\pi }}\int\limits_0^{2\pi } {\int\limits_0^\pi {{{{\left| {g(\theta ,\varphi )} \right|}}^{2}}\sin \theta d\theta d\varphi } } .$Для оценки точности выполнения оптической теоремы рассчитывалась относительная разность левой и правой частей в формуле (85):
(87)
${{\Delta }_{{{\text{отн}}}}} \equiv \frac{{\left| {\sigma + \operatorname{Im} \left( {g({{\theta }_{0}},{\text{0}})} \right)} \right|}}{\sigma }.$Результаты этой проверки для первого и второго из предлагаемых подходов при рассмотренных выше параметрах приведены в табл. 1 и 2. Как было отмечено ранее, для условия Неймана на экране точность выполнения оптической теоремы ниже, чем для условия Дирихле, поэтому рассматривался случай абсолютно жесткого экрана. Из табл. 1 видно, что при использовании первого подхода относительная разность правой и левой частей равенства (87) – величина ${{\Delta }_{{{\text{отн}}}}}$ – уменьшается с ростом числа точек коллокации ${{N}_{1}}$ и не превышает 3 × 10–2. При одинаковых значениях ${{N}_{1}}$ для тела меньших размеров точность выполнения оптической теоремы выше. Для второго подхода оптическая теорема выполняется с более высокой точностью в случае большего размера тела, а именно не превосходит 3 × 10–3 для ${{k}_{1}}a = 3$ и 2 × 10–4 для ${{k}_{1}}a = 6$.
Таблица 1.
${{N}_{1}},\;{{N}_{2}} = 4{{N}_{1}}$ | $\sigma $ | $ - \operatorname{Im} \left( {g({{\theta }_{0}},{\text{0}})} \right)$ | ${{\Delta }_{{{\text{отн}}}}}$ |
---|---|---|---|
${{k}_{1}}a = 3$ | |||
15 | 2.489595 | 2.515150 | 1.03 ×10–2 |
25 | 2.544564 | 2.564611 | 7.88 ×10–3 |
35 | 2.570533 | 2.587121 | 6.45 ×10–3 |
40 | 2.578923 | 2.594292 | 5.96 ×10–3 |
${{k}_{1}}a = 6$ | |||
15 | 14.917236 | 15.227829 | 2.08 ×10–2 |
25 | 15.002267 | 15.247020 | 1.63 ×10–2 |
35 | 15.060940 | 15.255580 | 1.29 ×10–2 |
40 | 15.081490 | 15.257716 | 1.17 ×10–2 |
Таблица 2.
L | $\sigma $ | $ - \operatorname{Im} \left( {g({{{\theta }}_{0}},{\text{0}})} \right)$ | ${{\Delta }_{{{\text{отн}}}}}$ |
---|---|---|---|
${{k}_{1}}a = 3$ | |||
16 | 2.634956 | 2.640329 | 2.04 ×10–3 |
32 | 2.636773 | 2.642852 | 2.30 ×10–3 |
48 | 2.633772 | 2.638545 | 1.81 ×10–3 |
64 | 2.632997 | 2.637633 | 1.76 ×10–3 |
${{k}_{1}}a = 6$ | |||
16 | 15.24457 | 15.25526 | 7.01×10–4 |
32 | 15.22389 | 15.21330 | 6.96 ×10–4 |
48 | 15.25218 | 15.26120 | 5.91×10–4 |
64 | 15.26292 | 15.26505 | 1.39 ×10–4 |
В качестве еще одной проверки разработанных алгоритмов, для второго из рассматриваемых подходов была вычислена невязка краевого условия, определяемая формулой:
(88)
$\begin{gathered} \Delta (\theta ,\varphi ) = \\ = \sum\limits_{m = - \infty }^\infty {\left[ {\int\limits_0^{{{\theta }_{J}}} {{{K}_{m}}} (\theta ,\theta {\kern 1pt} '){{I}_{m}}(\theta {\kern 1pt} ')\sin \theta {\kern 1pt} 'd\theta {\kern 1pt} '\,\, - {{B}_{m}}(\theta )} \right] \times } \\ \times \,\,\exp (im\varphi ), \\ \end{gathered} $На рис. 6 и 7 изображены угловые зависимости интенсивности рассеянного поля от угла раскрыва сферического экрана. Кривая 1 соответствует дифракции на проницаемом шаре, который не покрыт экраном. Соответствующая кривая получена при помощи подстановки в формулу (84) диаграммы первичного поля ${{g}^{0}}(\theta ,\varphi )$ вместо $g(\theta ,\varphi )$. Кривой 6 на рисунках показаны зависимости интенсивности при дифракции плоской волны на абсолютно мягкой (рис. 6) или абсолютно жесткой (рис. 7) сфере соответствующего волнового размера, полученные при помощи модифицированного метода дискретных источников (ММДИ) [17]. Кривые 2–5 соответствуют углам раскрыва экрана, равным 45°, 90°, 135° и 179°. Волновой размер сферы Януса и угол падения волны – ${{k}_{1}}a = 6,{\text{ }}{{{\theta }}_{0}} = 0$. Материальные параметры сред и волновые числа те же, что и для предыдущих рисунков. Результаты, приведенные на рисунках, получены при помощи первой из рассматриваемых методик при ${{N}_{1}} = 40,{\text{ }}{{N}_{2}} = 160$. На рисунках видно, что в случае, когда экран почти полностью покрывает шар (кривая 6), график интенсивности рассеянного поля совпадает с результатами для идеально отражающей сферы, что соответствует физической картине рассматриваемого явления. Видно также, что в случае абсолютно мягкого и абсолютно жесткого экранов имеется резкий максимум интенсивности в направлении угла падения плоской волны. В случае условия Дирихле величина максимума имеет наибольшее значение для ${{{\theta }}_{{\text{J}}}} = 179^\circ $ (когда экран вырождается в сферу). При отсутствии экрана (то есть при дифракции на проницаемом шаре) максимум в направлении падения волны имеет еще большее значение (чем для покрытого шара). В направлении обратного рассеяния (при ${\theta } = 180^\circ $) также имеется максимум интенсивности, который принимает наибольшие значения при ${{{\theta }}_{{\text{J}}}} = 45^\circ $ и ${{{\theta }}_{{\text{J}}}} = 90^\circ $. В случае абсолютно жесткого экрана величина максимума интенсивности в направлении падения плоской волны имеет существенно большее значение для ${{{\theta }}_{{\text{J}}}} = 45^\circ $ по сравнению с другими углами раскрыва экрана. Уровень обратного рассеяния также максимален при ${{{\theta }}_{{\text{J}}}} = 45^\circ $.
ЗАКЛЮЧЕНИЕ
На основе МПГУ разработаны два алгоритма численного решения задачи дифракции плоской волны на сфере Януса в виде проницаемого шара, частично покрытого абсолютно мягким или абсолютно жестким сферическим экраном. Выполнено сравнение результатов расчета интенсивности рассеянного поля, полученной с использованием предложенных методов с результатами, найденными методом T-матриц. Показано хорошее совпадение результатов как при применении обеих методик на основе МПГУ, так и при использовании метода Т-матриц. Проведена проверка точности выполнения оптической теоремы и невязки краевого условия в случае условий Дирихле и Неймана на экране. Показано, что погрешность выполнения оптической теоремы не превосходит $3 \times {{10}^{{ - 2}}}$ при использовании первой методики и $3 \times {{10}^{{ - 3}}}$ в случае применения второго метода. Показано, что невязка краевого условия при применении второй методики имеет малые значения на большей части осевого сечения экрана. Таким образом, можно сделать вывод о том, что вторая методика дает более точные результаты. Построены и исследованы угловые зависимости интенсивности рассеянного поля для различных углов раскрыва отражающего экрана. Показано существенное различие между поведениями угловых зависимостей интенсивности в случае абсолютно мягкого и абсолютно жесткого экранов. Установлено, что для выбранных параметров задачи максимальный уровень обратного рассеяния достигается при ${{{\theta }}_{{\text{J}}}} = 45^\circ $.
Работа выполнена при частичной финансовой поддержке Российского фонда фундаментальных исследований (проекты № 18-02-00961, 19-02-00654).
Список литературы
Zhang J., Grzybowski B.A., Granick S. Janus particle synthesis, assembly, and application // Langmuir. 2017. V. 33. № 28. P. 6964–6977.
Lattuada M., Hatton T.A. Synthesis, properties and applications of Janus nanoparticles // Nano Today. 2011. V. 6. № 3. P. 286–308.
Kim D., Avital E.J., Miloh T. Sound scattering and its reduction by a Janus sphere type // Advances in Acoustics and Vibration. 2014. V. 2014. Article ID 392138.
Gillman A. An integral equation technique for scattering problems with mixed boundary conditions // Advances Comput. Math. 2017. V. 43. P. 351–364.
Шестопалов В.П. Сумматорные уравнения в теории дифракции. Киев: Наук. думка, 1983.
Hawkins S.C., Rother T., Wauer J. A numerical study of acoustic scattering by Janus spheres // J. Acoust. Soc. Am. 2020. V. 147. № 6. P. 4097–4105.
Rother T. Sound Scattering on Spherical Objects. Heidelberg: Springer, 2020.
Кюркчан А.Г., Анютин А.П. Метод продолженных граничных условий и вейвлеты // Докл. Росс. Акад. наук. 2002. Т. 385. № 3. С. 309–313.
Кюркчан А.Г., Смирнова Н.И. Математическое моделирование в теории дифракции с использованием априорной информации об аналитических свойствах решения. М.: Медиа Паблишер, 2014.
Кюркчан А.Г., Маненков С.А. Гибридный подход к решению задачи дифракции на плоских экранах // Акуст. журн. 2015. Т. 61. № 3. С. 302–310.
Кюркчан А.Г., Маненков С.А. Решение задачи дифракции на плоском экране, расположенном в плоскослоистой среде, с помощью метода продолженных граничных условий // Радиотехника и электроника. 2020. Т. 65. № 7. С. 644–652.
Gerolymatos P.G., Manenkov A.B., Tigelis I.G., Amditis A.J. Metal iris influence on guided-mode diffraction // J. Opt. Soc. Am. 2006. V. 23. № 6. P. 1333–1339.
Маненков С.А. Применение различных базисов при решении задачи дифракции на незамкнутых экранах // Радиотехника и электроника. 2014. Т. 59. № 3. С. 246–252.
Справочник по специальным функциям с формулами, графиками и математическими таблицами // Под ред. Абрамовица М. и Стиган И.М. М.: Наука, 1979. 832 с.
Маненков С.А. Задача дифракции электромагнитного поля на неоднородном теле с осевой симметрией // Радиотехника и электроника. 2018. Т. 63. № 1. С. 3–13.
Шендеров Е.Л. Излучение и рассеяние звука. Л.: Судостроение, 1989.
Маненков С.А. Новая версия модифицированного метода дискретных источников применительно к задаче дифракции на теле вращения // Акуст. журн. 2014. Т. 60. № 2. С. 129–136.
Дополнительные материалы отсутствуют.
Инструменты
Акустический журнал