Всероссийская открытая научная конференция «Современные проблемы дистанционного зондирования, радиолокации, распространения и дифракции волн» - Муром 2024
УДК 517.598 БОТ: 10.24412/2304-0297-2024-1-140-144
Математические методы исследования резонансного рассеяния на сферических диэлектрических частицах
Ю.В. Шестопалов
Российский Технологический Университет МИРЭА 119454, г. Москва, пр. Вернадского 78. E-mail: shestopalov@mirea. ru
Представлены результаты разработки математических методов, обеспечивающих точное решение многопараметрических обобщенных дисперсионных уравнений (ОДУ) электромагнетизма в комплексной области, в том числе для задач резонансного рассеяния сферическими диэлектрическими частицами. Разработанный численно-аналитический подход позволяет доказать существование и получать резонансные частоты без вычислений для широкого семейства структур планарной, круговой и сферической симметрий. В методе используются анализ функций нескольких комплексных переменных, применение неявных отображений и решение начально-краевых задач для дифференциальных уравнений с заданной точностью.
Ключевые слова: резонанс, рассеяние, неявная функция, комплексная область, сферические диэлектрические частицы
Mathematical methods for the analysis of resonance scattering by spherical dielectric particles
Y.V. Shestopalov
1 Russian Technological University MIREA.
The results are presented of the development of mathematical methods providing exact solution to multiparameter generalized dispersion equations (GDEs) of electromagnetics in the complex domain associated in particular with the resonance wave scattering by spherical dielectric particles. The elaborated numerical-analytical approach enables one to obtain resonance frequencies without computations for a broad family ofplanar-, circular-, and spherical-symmetric structures and employs analysis offunctions of several complex variables, application of implicit mappings and solution to initial-value problems with prescribed accuracy.
Keywords: resonance, scattering, implicit function, complex domain, spherical dielectric particles
Введение
Определение резонансных частот (РЧ) для широкого класса задач рассеяния сводится [1-3] к решению обобщенных дисперсионных уравнений (ОДУ) вида
= 0, ш = (ш;и), где № - комплексный частотный спектральный параметр;
и = (и)у=1 - вектор N > 1 неспектральных параметров (определяющих материальное заполнение и геометрию рассеивателя).
Математически корректно, с точки зрения теории функций нескольких комплексных переменных и неявных комплекснозначных отображений, задача решения ОДУ сводится к определению неявных функций (гиперповерхностей) № = ш(и). Если варьировать только один параметр £ = иу , как это обычно происходит в прикладных исследованиях, то целью является определение неявных комплексных отображений № = — зависимостей спектра РЧ от комплексной переменной
Алгоритм расчета РЧ должен базироваться на возможности корректного сведения задачи к решению функционального или операторного ОДУ. В частности, этот алгоритм может быть основан на анализе коэффициентов разложения в представлениях рассеянного поля вне круга, содержащего сечение металлодиэлектрического рассеивателя, в виде обобщенного ряда Фурье по уходящим цилиндрическим гармоникам
р(р.в)=1п Тп(т) Hn\k0p)eine, где (р, в) - комплексный частотный спектральный параметр;
к0 - волновое число вакуума;
(1)
Нп - функции Ханкеля первого рода порядка п;
р(р, в) - потенциальная компонента рассеянного поля.
Отметим, что такой подход развивался многими авторами, в том числе в [4] и в работах, цитируемых в этой публикации.
Особенности Tn(w) могут классифицироваться как резонансы рассеивающей структуры. Если рассеиватель обладает радиальной симметрией и кусочно-постоянным диэлектрическим заполнением, то коэффициенты разложения можно представить в явном виде как отношения Tn(w) = Tn(w)/xn(^) обобщенных цилиндрических полиномов [5], которые вычисляются явно [1,2]. Обращение в нуль числителя или знаменателя Тп означает [1] соответственно подавление n-й гармоники или резонанс рассеянного поля.
Резонансное рассеяние на сферических диэлектрических частицах является предметом активных исследований уже более 150 лет, начиная с классических работ лорда Рэлея. Многие результаты были обобщены в монографии [6]. Однако до сих пор, насколько нам известно, не были получены строгие доказательства существования РЧ и не изучено их расположение на комплексной плоскости. В [6, 7] и многих других работах анализ свойств РЧ и резонансного рассеяния проводится только на основе результатов, полученных путём численного решения уравнений в комплексной области с использованием коммерческих кодов без обоснования существования корней этих уравнений и исследования их расположения на комплексной плоскости. Такой подход конечно недостаточен для построения математически завершённой теории резонансного рассеяния на диэлектрических частицах и требует соответствующего развития. В настоящей работе такие доказательства проведены по-видимому впервые с использованием метода, основанного на анализе критическихреперных точек (КРТ).
Расчет РЧ в терминах w(u) и Z(Q проводится путем численного решения соответствующих задач Коши для искомых неявных функций с начальными данными, задаваемыми КРТ, которые определяются явно. Последнее является необходимым условием эффективного применения описанного метода вычисления и реализуется на практике для широкого класса структур, обладающих координатной симметрией.
Для простейшего рассеивателя — однородного кругового диэлектрического
цилиндра — ОДУ имеет вид [1, 2]
(1) (1) G(w) = аН^ )(w)Jn(aw) — Н^ )(aw)Jn(w) = 0, ш = (w; а),
с левой частью в виде обобщённого цилиндрического полинома, а КРТ определяются
явно как
(w;a) = (w*;a*) = [hH\ v?/h£], где w = k0a - комплексный частотный спектральный параметр;
а - радиус диэлектрического цилиндра;
е, а = Je - диэлектрическая проницаемость и коэффициент преломления диэлектрика;
(1) (1) и fo^ - нули функций Бесселя /п и Ханкеля Я^ .
Для многослойных радиально-симметричных структур КРТ также вычисляются явно как отношения произведений нулей цилиндрических функций.
Полученные явно без вычислений значения РЧ кругового диэлектрического цилиндра в терминах КРТ сами по себе имеют физический смысл, так как позволяют интерпретировать КРТ как ''суперпозицию'' резонансов внутренней и внешней областей рассеивателя.
РЧ, полученные в явном виде и рассматриваемые как неявные функции (отображения) одного или нескольких комплексных параметров, обобщают понятие дисперсионных кривых — зависимостей РЧ от одного вещественного параметра — и являются способом описания топологии спектра РЧ на комплексной плоскости и динамики изменения РМ в зависимости от параметров.
Результаты [1 — 3] показывают, что исследование такой динамики может раскрыть ранее не изученные свойства спектра и найти важные приложения в анализе резонансного рассеяния.
Математические методы анализа ОДУ для сферической диэлектрической частицы
Изложенная методика эффективно применяется для анализа и определения РЧ сферической диэлектрической частицы. Эта задача сводится [6, 7] к решению уравнения
5n(w) = [wfon)(w)J/n(aw) — ^ [z^)n(zi)]fojl1)(w) = 0, (1)
где w = fc0P - комплексный частотный спектральный параметр;
R - радиус сферической диэлектрической частицы;
е - диэлектрическая проницаемость, а = V^M (далее в данной работе будем предполагать, что д = 1) и z1 = aw;
и fo(l1) - сферические функции Бесселя первого и третьего рода.
После преобразований уравнение (1) сводится при п = 1 к решению эквивалентного уравнения с левой частью в виде обобщённого тригонометрического полинома [5]
7\(w; е) = 4(1)(w; е) sin(aw) + fl(1)(w; е) cos(aw) = 0, (2)
Л(1) (w; е) = £е #(w) + 2(w — t)(ew2 — 1), 5(1)(w; e) = aw[—ie ^(w) — w + i], ^(w) = w2 — iw — 3, i2 = —1.
Уравнения (1) и (2) определяют w = w(e) или w = w(a) как неявные функции параметра е или а. Можно показать, что функция w(e) (или w = w(a)) существует и ее (то есть искомые комплексные РЧ сферической диэлектрической частицы) можно вычислить с заданной точностью, так как известны КРТ (w*; е*), при которых 7\(w*; е*) = 0. Эти КРТ можно найти явно из (а) вида уравнения (1): w* = t, е* = —(Й^)2, s = 1,2, ... (при п = 1),
£(1), s = 1,2, ..., обозначают нули сферических функции Бесселя j'jXz) первого рода; и (б) из вида уравнения (2),
п
aw = -(2fc + 1), fc = 0,±1,±2,... ,
Pfc(4)(w) = 2w3(afc — 1) + w2t(2 — afc) + wafc — 3iafc = 0,
или
>(B)
аш = л"п,п = ±1, ±2,..., Рп( = ш3 + ш2/(Ьп — 1) + шЬп — 3/ЬП = 0. Можно показать, что полиномы РА(л)(ш) и Р^5"*^) имеют для каждого к = 0, ±1, ±2,... и п = ±1, ±2,... один чисто мнимый нуль (р = 1 в Таблице 1 для Р,.(л)(ш))
на отрицательной мнимой полуоси и два комплексных нуля (р = 2.3 в Таблице 1), расположенных симметрично в первом и втором квадрантах комплексной плоскости w.
Таким образом, можно явно определить счётные совокупности РЧ как КРТ. Вычисленные значения этих нулей — искомых комплексных РЧ индекса п = 1 сферической диэлектрической частицы для случая (а) и к = 0 и соответствующие согласно представлениям КРТ значения проницаемости е — представлены в Таблице 1. Отметим, что для р = 1 РЧ в Таблице 1 определяют колебания, потенциальная компонента рассеянного поля которых экспоненциально растёт на бесконечности, а для р = 2,3 убывает и удовлетворяет условию излучения Зоммерфельда.
Материалы с проницаемостями из Таблицы 1, вычисленными согласно значениям КРТ, могут быть классифицированы как специальный класс диамагнетиков. Однако в данном случае важны не сами значения КРТ как таковые, а КРТ в качестве точных начальных значений для последующего вычисления РЧ в физически оправданном диапазоне изменения проницаемости, то есть определение РЧ на соответствующем интервале изменения е.
Действительно далее, РЧ — искомая неявная функция w = w(e) параметра е — вычисляются как решение задачи Коши в комплексной области; например, при п = 1 эта задача Коши имеет вид
де
I I \л/ ■ (- I zz — ■
rdw
f(w; е) =-
дТ\ dw
s = 1.2. 3.
(3)
(^(1)) . так как выполняется условие её однозначной разрешимости — (^
и решается на некотором интервале (— (fs1*) . Е0) изменения параметра е c начальными условиями в КРТ, полученных из условий (а), как задача (3), или (б).
Решения задач Коши (3) существуют в окрестности начальной точки е =
0. и определяют искомые неявные функции, которые являются решениями уравнения Ti(w;e) = 0 при е е (—(^f. Ео).
Наличие КРТ и однозначная разрешимость задач Коши для определения неявных функций w = w(e) или w = w(a) образует в совокупности строгое доказательство существование РЧ и позволяет найти их расположение на комплексной плоскости.
Приближённое вычисление неявной функции w = w(e) как решения (3) проводится с использованием метода Рунге-Кутта при параметризации вещественной переменной е,
(1) 7 "
исходящей из реперной точки е = —(Q )2. расположенной на вещественной оси комплексной плоскости е. Задавая значение Е0 > 1 можно в результате расчетов вычислить РЧ в конечной точке интервала изменения проницаемости для диэлектриков без потерь.
Таблица 1. Комплексные РЧ индекса п = 1 сферической диэлектрической частицы для случая (а) и к = 0 и соответствующие КРТ-значения проницаемости е
p РЧ Значения проницаемости е
1 —1.50775i —1.085378
2 —0.988990 + 0.833506i 0.249854 + 1.453667/
3 0.988990 + 0.833506i 0 .249854 — 1.453667/
Выводы
Разработан численно-аналитический метод расчета РЧ диэлектрических структур планарной, круговой и сферической симметрий, основанный на точном решении многопараметрических ОДУ в комплексной области. Специальные исследования проведены для задач резонансного рассеяния сферическими диэлектрическими частицами. Получены строгие доказательства существования РЧ и изучено их расположение на комплексной плоскости.
Развитый численно-аналитический подход позволяет находить КРТ и РЧ в явном виде для широкого семейства структур, в том числе путём численного определения РЧ как неявных функций — решений соответствующих начально-краевых задач для дифференциальных уравнений с заданной точностью. В этом качестве разработанные методы исследования ОДУ существенно расширяют математический и вычислительный аппарат известных коммерческих продуктов и программных кодов и могут быть использованы для обоснования достоверности результатов и повышения эффективности расчетов.
Для проведения вычислений РЧ как неявных функций с заданной точностью предлагается применять численное решение задач Коши методами Рунге-Кутта высокого порядка. Такие методы после их внедрения могут также существенно дополнить имеющиеся программные коды исследования задач электромагнетизма.
Работа выполнена при поддержке гранта РНФ № 20-11-20087.
Литература
1. Y. Shestopalov, Cloaking: analytical theory for benchmark structures (2021), J. Electr. Waves Appl.. 35, 485-510.
2. Y. Shestopalov, Resonance Scattering by a Circular Dielectric Cylinder (2021), Radio Science, 56, e2020RS007095.
3. Y. Shestopalov, Application of Parametric Fourier Series to the Analysis of Partial Invisibility and Resonance Scattering by Canonical Structures, Proc. 2022 IEEE Int. Conf. on Electromagnetics in Advanced Applications (ICEAA 2022), 5-9 September 2022, Cape Town, South Africa, 47-50.
4. Анютин А.П. Плазмонные резонансы в выпукло-вогнутом наноцилиндре из серебра // Радиотехника и Электроника. 2021. T. 66. № 6. C. 559-564.
5. Y.Shestopalov, Trigonometric and cylindrical polynomials and their applications in electromagnetics (2019), Applicable Analysis, 99, 2807-2822.
6. Климов В.В. Наноплазмоника - М.: Физматлит, 2009. 235 c.
7. V. V. Klimov, A. R. Bekirov, B. S. Lukyanchuk, Trapped modes in particles with a negative refractive index (2023), Optics Letters, 48, 5795-5798.