Получение точных аналитических решений с помощью эвристического метода базовых компонентов

Язык труда и переводы:
УДК:
537.874.6
Дата публикации:
23 ноября 2022, 16:37
Категория:
Перспективные направления исследования необратимых физических процессов
Авторы
Весник Михаил Владимирович
ИРЭ им. В.А. Котельникова РАН
Аннотация:
Впервые с помощью недавно разработанного эвристического метода базовых компонентов получено точное (по сравнению со строгим) аналитическое решение задачи дифракции электромагнитной волны на полуплоскости с двухсторонними импедансными граничными условиями. Под термином «точные» мы подразумеваем формулы, математическая строгость которых не доказана. Поскольку верификацию проводили для всей области входных параметров, нет оснований полагать, что их точность отличается от точности строгого решения.
Ключевые слова:
эвристические подходы, двухсторонние импедансные граничные условия, метод Винера – Хопфа, физическая теория дифракции, метод базовых компонентов, дифракция электромагнитных волн, аналитические решения
Основной текст труда

Введение

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

С целью увеличения эффективности (точности и быстродействия) солверов применяют эвристические решения [1–4] и гибридные методы, сочетающие в себе численные и эвристические подходы [5].

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

К традиционным эвристическим подходам относят: метод геометрической оптики (ГО) [6], метод физической оптики (ФО) [6], геометрическую теорию дифракции (ГТД) [7–9], метод краевых волн (МКВ) [10–12]. Известно, что подход МКВ оказался критически важным при реализации проекта «Стелс». Термин физическая теория дифракции (ФТД) можно применять как для МКВ, так и в расширительном толковании — для всех эвристических методов. Краткое обсуждение особенностей этих методов приведено в [13, 14]. Традиционные эвристические подходы (ГО, ФО, ГТД, МКВ) предлагают эвристические формулы, но не предполагают их дальнейшего уточнения.

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

Недавно разработанный автором метод базовых компонентов (МБК), в отличие от традиционных эвристических подходов, предполагает процесс уточнения («настройки») первичных эвристических формул. При построении первичных эвристических формул применяют решения простейших математически строгих задач дифракции, любые другие проверенные решения (в том числе — приближения традиционных эвристических подходов), а также уже существующие наработки МБК по уточнению аналогичных эвристических формул. Настройку проводят с помощью сравнения первичных эвристических формул с верификационным решением (как правило, численным).

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

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

Двухсторонние импедансные граничные условия давно известны в теории дифракции [15–18]. Непрозрачная полуплоскость описана в [19]. В данной работе рассмотрена полупрозрачная полуплоскость, как и в [20, 21].

Первичные эвристические формулы

Применим МБК для получения эвристического решения задачи дифракции на полуплоскости с двухсторонними импедансными граничными условиями. В качестве первичных эвристических формул выберем выражения [3, 4, 20, 21].

Для TH — поляризации:

fg(R_{TH},T_{TH},\phi ,\phi _{0})={\frac {1}{2}}\left({\frac {1-T_{TH}}{-\cos {\frac {\phi -\phi _{0}}{2}}}}+{\frac {R_{TH}}{-\cos {\frac {\phi +\phi _{0}}{2}}}}\right).                                          (1)

для TE — поляризации:

gf(R_{TE},T_{TE},\phi ,\phi _{0})={\frac {1}{2}}\left({\frac {1-T_{TE}}{-\cos {\frac {\phi -\phi _{0}}{2}}}}+{\frac {R_{TE}}{-\cos {\frac {\phi +\phi _{0}}{2}}}}\right).                                          (2)

Здесь \phi _{0} и \phi — углы точки наблюдения и источника;  R_{TH} и R_{TE} — коэффициенты отражения от безграничной плоскости.

Назовем (1) и (2) обобщенными дифракционными коэффициентами (ОДК) [20].

Верификационное решение

В качестве верификационного решения задачи дифракции на полуплоскости с двухсторонними импедансными граничными условиями используем строгое решение по методу Винера — Хопфа (МВХ) [18, 21]. Будем называть его «решение по МВХ».

Для TH — поляризации решение по МВХ:

fr(X,\phi ,\phi _{0}).                                                               (3)

Для TE — поляризации решение по МВХ:

gr(X,\phi ,\phi _{0}).                                                              (4)

Настроечные функции и эвристические формулы

Способы поиска настроечных функций могут отличаться. В работах [20, 21] эвристические формулы для функций полупрозрачности были выбраны с целью получить наиболее физически ясные и компактные выражения.

Принцип взаимности

Двумерные структуры с неидеальными граничными условиями были рассмотрены в ряде статей [22–25]. Эвристическое решение задачи дифракции для импедансного клина строили с помощью подстановки френелевских коэффициентов отражения в формулы работы [26], затем разными способами пытались добиться наилучшего соответствия решения принципу взаимности и, наконец, верифицировали формулы с помощью решения Малюжинца [27].

В работе [28] были получены новые настроечные функции, обеспечивающие очень высокую точность:

cx_{TH}(m_{TH},X,\phi )=R_{TH}(2X\cdot m_{TH},\phi /2);                                   (5)

cx_{TE}(m_{TE},X,\phi )=R_{TE}(0,5X\cdot m_{TE},\phi /2).                                   (6)

Здесь m_{TH} , m_{TE}   и X — параметры, описывающие граничные условия. Внося изменения в аргументы функций R_{TH}  и R_{TE} , можно получить выражения (5) и (6) для настроечных функций, расчет по которым дает результаты, близкие к результатам верификационного расчета (3) и (4). В результате настройки получаем эвристические формулы:

fh(m_{TH},R_{TH},T_{TH},\phi ,\phi _{0})=fg(R_{TH},T_{TH},\phi ,\phi _{0}){\frac {cx_{TH}(m_{TH},X,\phi )}{cx_{TH}(m_{TH},X,\phi _{s})}}   (TH);           (7)

gh(m_{TE},R_{TE},T_{TE},\phi ,\phi _{0})=gf(R_{TE},T_{TE},\phi ,\phi _{0}){\frac {cx_{TE}(m_{TE},X,\phi )}{cx_{TE}(m_{TE},X,\phi _{s})}}     (TE).           (8)

Здесь углы \phi _{s}=\pi \pm \phi _{0} соответствуют границам «свет — тень» падающего или отраженного поля. Для сравнения (по терминологии МБК — верификации) решения по МВХ с эвристическими формулами (7) и (8) построим отношения

{\frac {fr(X,\phi ,\phi _{0})}{fh(m_{TH},R_{TH},T_{TH},\phi ,\phi _{0})}};                                                      (9)

{\frac {gr(X,\phi ,\phi _{0})}{gh(m_{TE},R_{TE},T_{TE},\phi ,\phi _{0})}}.                                                      (10)

Чем ближе (9) и (10) к единице, тем выше точность.

Количественная оценка точности

Для того чтобы оценить точность эвристических выражений (7) и (8), применяют формулы количественной оценки точности [28]:

Lfh(m_{TH},R_{TH},T_{TH},\phi ,\phi _{0})=\lg {\left({\frac {1}{2\pi }}\int \limits _{\phi =0}^{2\pi }\left|\left|{\frac {fr(X,\phi ,\phi _{0})}{fh(m_{TH},R_{TH},T_{TH},\phi ,\phi _{0})}}\right|-1\right|d\phi \right)};     (11)

Lgh(m_{TE},R_{TE},T_{TE},\phi ,\phi _{0})=\lg {\left({\frac {1}{2\pi }}\int \limits _{\phi =0}^{2\pi }\left|\left|{\frac {gr(X,\phi ,\phi _{0})}{gh(m_{TE},R_{TE},T_{TE},\phi ,\phi _{0})}}\right|-1\right|d\phi \right)}.     (12)

 

Настроечные параметры mTH и mTE

Приведем таблицу значений настроечных параметров m_{TH}  и m_{TE} , найденных феноменологическим способом путем минимизации интегралов (11) и (12).

 

X

mTH

mTE

φ0 = 160º

φ0 = 90º

φ0 = 20º

ПВ

φ0 = 160º

φ0 = 90º

φ0 = 20º

ПВ

10

3,1560

3,1528

3,1555

3,1525

0,9977

0,9997

0,9999

0,9996

150

1,1450

1,1412

1,1423

1,1415

0,9622

0,9371

0,9377

0,9370

300

1,0484

1,0437

1,0443

1,0443

0,8632

0,8334

0,8348

0,8334

500

1,0271

1,0163

1,0169

1,0170

0,7581

0,7221

0,7244

0,7222

1000

1,0436

1,0036

1,0041

1,0044

0,6092

0,5585

0,5601

0,5591

 

Применяя принцип взаимности к формулам (7) и (8), получим аналитические выражения для m_{TH} и m_{TE} [28]:

m_{TH}={\sqrt {\frac {1+{\sqrt {1+{\frac {1}{4}}(W_{0}/X)^{2}}}}{2}}};                                                   (13)

m_{TE}={\sqrt {\frac {-1+{\sqrt {1+4(X/W_{0})^{2}}}}{2(X/W_{0})^{2}}}}.                                                (14)

Здесь W_{0}=120\pi . Значения для m_{TH} и m_{TE} (13) и (14) приведены в таблице в столбцах «ПВ». 

Обсуждение

В данной работе эвристическое решение задаче дифракции на полуплоскости с двусторонними импедансными граничными условиями в форме (7) и (8) получено двумя способами. Отличие состоит в том, как были найдены настроечные параметры    m_{TH} и m_{TE} .

Первый способ — численное определение минимумов параметров количественной оценки (11) и (12). Второй способ — аналитическое определение   m_{TH} и m_{TE}  в виде (13) и (14).

Для заданного вида граничных условий мы получили точные аналитические выражения (7) и (8) с подстановкой (13) и (14), которые можно использовать в практических задачах вместо численного решения. Под термином «точные» мы подразумеваем формулы, математическая строгость которых не доказана. Поскольку верификацию проводили для всей области входных параметров, нет оснований полагать, что их точность отличается от точности строгого решения.

При наличии соответствующих верификационных решений, развитие МБК может происходить в направлениях получения точных эвристических решений, описывающих влияние профиля кромки [3, 4, 29–33], особенностей контура рассеивателя [3, 4, 34–39] и граничных условий [3, 4, 14, 19–21, 28, 40].

Выводы

Для некоторых задач применение МБК позволяет получить аналитические решения, точность которых соответствует точности строгих аналитических решений.

Благодарность

Автор выражает благодарность С.Е. Банкову за предоставление программ для расчета решения по методу Винера — Хопфа и аналитических формул для коэффициентов отражения для безграничной полуплоскости.

Литература
  1. Хенл Х., Мауэ А., Вестпфаль К. Теория дифракции. Москва, Мир, 1964, 428 с.
  2. Kravtsov Y.A., Zhu N.Y. Theory of diffraction: heuristic approaches. Oxford, Alpha Science Int., 2010.
  3. Vesnik M.V. The method of the generalized Eikonal: new approaches in the diffraction theory. Berlin, Walter de Gruyter, 2015.
  4. Весник М.В. Физическая интерпретация математически строгого решения задачи дифракции при помощи эвристических формул. Современная математика. Фундаментальные направления, 2016, т. 62, с. 32–52. URL: http://www.mathnet.ru/links/19673f78aa1f8ed046fe133eda63b15c/cmfd308.pdf
  5. Кюркчан А.Г., Смирнова Н.И., Математическое моделирование в теории дифракции с использованием априорной информации об аналитических свойствах решения. Москва, Медиа Паблишер, 2014, 226 с.
  6. Born M., Wolf E. Principles of optics. Oxford, London, Pergamon Press, 1968.
  7. Keller J.B. Geometrical theory of diffraction. J Opt Soc Am, 1962, vol. 52 (2), pp. 116–130.
  8. James G.L. Geometrical Theory of Diffraction for Electromagnetic Waves. London, Peter Peregrinus Ltd., 1976.
  9. Borovikov V.A., Kinber B.E. Geometrical Theory of Diffraction. London, UK, IEE, 1994.
  10. Уфимцев П.Я. Метод краевых волн в физической теории дифракции. Москва, Советское радио, 1962, 243 с.
  11. Ufimtsev P.Y. Theory of Edge Diffraction in Electromagnetics. Encino, California, Tech Science Press, 2003.
  12. Ufimtsev P.Y. Fundamentals of the Physical Theory of Diffraction. Hoboken, N.J., John Wiley & Sons, Inc., 2007.
  13. Pelosi G., Rahmat-Samii Y., Volakis J.L. High-frequency techniques in diffraction theory: 50 years of achievements in GTD, PTD, and related approaches. IEEE Antennas Propag Magazine, 2013, vol. 32. p. 16.
  14. Vesnik M.V. New possibilities for constructing heuristic solutions to problems of electromagnetic diffraction. MDPI, Eng, 2022, vol. 3, pp. 27–41. DOI: https://doi.org/10.3390/eng3010004
  15. Senior T.B.A. Half plane edge diffraction. Radio Sci, 1975, vol. 10 (6), p. 645.
  16. Senior T.B.A. Diffraction tensors for imperfectly conducting edges. Radio Sci, 1975, vol. 10 (10), p. 911.
  17. Войтович Н.Н., Каценеленбаум Б.З., Коршунова Е.Н., Пангонис Л.И., Переяславец М.Л., Сивов А.Н., Шатров А.Д. Электродинамика антенн с полупрозрачными поверхностями: методы конструктивного синтеза. Москва, Наука. Гл. ред. физ.-мат. лит., 1989, 176 с.
  18. Банков С.Е. Интегральная СВЧ-оптика. Москва, ФИЗМАТЛИТ, 2018.
  19. Банков С.Е., Весник М.В., Кравченко В.Ф. Эвристическое решение задачи дифракции на сверхпроводящей полуплоскости. Радиотехника и электроника, 2020, т. 65, № 4, с. 363–371. DOI: https://doi.org/10.31857/S0033849420040014
  20. Весник М.В. Применение метода базовых компонентов для получения эвристического решения задачи дифракции на полуплоскости с неидеальными граничными условиями. Радиотехника и электроника, 2019, т. 64, № 11, с. 1103–1109. DOI: https://doi.org/10.1134/S003384941911024X
  21. Vesnik M.V., Bankov S.E. Heuristic solution to the problem of diffraction of a TE-polarized electromagnetic wave on a semitransparent half-plane. Waves in Random and Complex Media, 2021. DOI: https://doi.org/10.1080/17455030.2021.1951888
  22. Luebbers R.J. Finite conductivity uniform GTD versus knife edge diffraction in prediction of propagation path loss. IEEE Trans Antennas Propag. 1984, vol. 32, pp. 70–76.
  23. Holm P.D. A new heuristic UTD diffraction coefficient for nonperfectly conducting wedges. IEEE Trans Antennas Propag, 2000, vol. 488.
  24. El-Sallabi H.M., Rekanos I.T., Vainikainen P. A new heuristic diffraction coefficient for lossy dielectric wedges at normal incidence. IEEE Antennas Wirel Propag Lett, 2002, vol. 1, pp. 165–168.
  25. Soni S., Bhattacharya A. New heuristic diffraction coefficient for modeling of wireless channel. Progr Electromagn Res C, 2010, vol. 12, pp. 125–137.
  26. Kouyoumjian R.G., Pathak P.H. A uniform geometrical theory of diffraction for an edge in a perfectly conducting surface. Proc IEEE, 1974, vol. 62, pp. 1448–1461.
  27. Малюжинец Г.Д. Возбуждение, отражение и излучение поверхностных волн от клина с заданными поверхностными импедансами, ДАН СССР, 1958, № 3, с. 752–755.
  28. Vesnik M.V. A technique for obtaining analytical heuristic solutions in problems of diffraction on two-dimensional semi-infinite objects with non-ideal boundary conditionsю Waves in Random and Complex Media, Aug 2022. DOI: https://doi.org/10.1080/17455030.2022.2108160
  29. Vesnik M.V. The Analytical Solution for the Electromagnetic Diffraction on 2-D Perfectly Conducting Scatterers of Arbitrary Shape. IEEE Transactions on Antennas and Propagation, 2001, vol. AP-49, no. 12, pp. 1638–1644.
  30. Весник М.В. Аналитическое решение краевых задач теории дифракции методом обобщенного эйконала. Радиотехника и электроника, 2003, т. 48, № 9, с. 1078–1084.
  31. Весник М.В. Аналитическое решение задачи дифракции электромагнитной волны на двумерной идеально проводящей полупластине при помощи метода обобщенного эй-конала. Радиотехника и электроника, 2008, т. 53, № 2, с. 144–156.
  32. Весник М.В. Аналитическое решение двумерной задачи дифракции электромагнитной волны на усеченном клине. Радиотехника и электроника, 2012, т. 57, № 10, с. 1053–1065.
  33. Vesnik M.V., Kravtsov Yu.A. Section 5.1.7. Diffraction by Bodies with Wedges: Method of Generalised Eikonal (MGE). Kravtsov Yu.A., Ning Yan Zhu. Theory of Diffraction: Heuristic Approaches. Alpha Science International Ltd.Oxford, U.K., 2010.
  34. Vesnik M.V., Ufimtsev P.Y. An Asymptotic Feature of Corner Waves Scattered by Polygonal Plates. Electromagnetics, 1992, vol. 12, no. 3–4, pp. 265–272.
  35. Весник М.В. О возможности построения уточненного эвристического решения в задаче дифракции на плоском угловом секторе. Радиотехника и электроника, 2011, т. 56, № 5, с. 573–586.
  36. Vesnik M.V. Efficiency of Different Heuristic Approaches to Calculation of Electromagnetic Diffraction by Polyhedrons and other Scatterers. Radio Science, 2014, vol. 49, iss. 10, pp. 945–953. DOI: https://doi.org/10.1002/2014RS005520
  37. Весник М.В., Калиничев В.И. Особенности электромагнитной дифракции на идеально проводящем плоском прямоугольнике. Сб. тр. V Всерос. микроволн. конф. Москва, ИРЭ им. В.А.Котельникова РАН, 2017, с. 114–118.
  38. Весник М.В. Свойства дифракционных коэффициентов в задаче дифракции на трехмерном плоском многоугольнике. Ч. 1: Основы метода базовых компонентов. Радиотехника и электроника, 2020, т. 65, № 11, с. 1052–1060. DOI: https://doi.org/10.31857/S0033849420110194
  39. Весник М.В. Свойства дифракционных коэффициентов в задаче дифракции на трехмерном плоском многоугольнике. Ч. 2: Применение метода базовых компонентов для решения практических задач. Радиотехника и электроника, 2020, т. 65, № 12, с. 1147–1159. DOI: https://doi.org/10.31857/S0033849420120189
  40. Vesnik M.V. Physical interpretation of the solution to the problem of diffraction on a half-plane with non-Ideal boundary conditions. Journal of Engineering Research and Sciences, 2022, vol. 1, no. 3, pp. 52–58. DOI: https://doi.org/10.55708/js0103006
Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.