CC BY
26
Математическое моделирование рассеянного поля в рабочей зоне
зеркального коллиматора
Ф. Б. Хлебников1,0, А. Н. Боголюбов1, В. С. Солосин2, Н.Е. Шапкина1
1 Московский государственный университет имени М. В. Ломоносова, физический факультет,
кафедра математики. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2. 2 ФГБУ науки Институт теоретической и прикладной электродинамики РАН (ИТПЭ РАН). Россия, 125412, Москва, ул. Ижорская, д. 13. E-mail: а mnfkhl@gmail.com
Статья поступила 14.07.2015, подписана в печать 01.09.2015.
Построена математическая модель зеркального коллиматора со скругленными краями. Методом интегральных уравнений решена задача расчета компонент электромагнитного поля в рабочей зоне. Представлены результаты численных расчетов.
Ключевые слова: математическое моделирование, зеркальные коллиматоры. УДК: 51-73. PACS: 02.30.Rz.
Введение
Одной из важных задач практической электродинамики является проблема получения в заданном объеме пространства поля, близкого по структуре к полю плоской электромагнитной волны. В частности, эта проблема возникает в задачах связанных с радиолокацией, в которых нужно исследовать процесс облучения тела удаленным источником. Сходные проблемы возникают при измерении диаграмм направленности антенн. Чтобы обеспечить поле, близкое по структуре к плоской волне, необходимо поместить исследуемое тело на достаточно большое расстояние от источника, в зону Фраунгофера, или, как принято говорить в радиолокации, в дальнюю зону. Для этого могут быть использованы открытые полигоны, однако влияние погодных условий, особенностей окружающей среды и других факторов существенно осложняет интерпретацию полученных результатов. Другим способом исследования является формирование поля, близкого к полю плоской волны, в закрытом пространстве измерительной лаборатории. В таких случаях применяют коллиматоры — устройства, преобразующие сферическую волну в плоскую. Измерительные стенды на базе безэховой камеры, в которых плоское поле формируется коллиматором, называются компактными полигонами [1].
Как и в оптике, в электродинамике можно выделить два типа коллиматоров: коллиматоры основанные на линзах и коллиматоры основанные на зеркалах. Существуют также и другие схемы формирования плоского поля, в частности, использующие решетки из простых антенн и голографические — основанные на использовании теневых масок типа зонной пластинки Френеля.
Коллиматоры в виде диэлектрических линз пока не получили широкого распространения из-за большого количества технических сложностей при изготовлении линзы, в частности требованиям к однородности материала линзы и точности обработки
ее поверхности, поэтому большая часть существующих коллиматоров для компактных полигонов — зеркальные коллиматоры.
Рабочей зоной коллиматора называют область, в которой отраженное поле можно с достаточной точностью считать полем плоской волны. Она представляет собой ограниченную цилиндрическую область, ось которой параллельна оси параболоида. Размеры рабочей зоны зависят от частоты падающей волны и размеров самого зеркала, при этом зеркало коллиматора должно быть достаточно большим, чтобы уменьшить влияние краевых дифракционных эффектов. Обычно поперечные размеры коллиматора приблизительно вдвое больше размеров рабочей зоны. Основными причинами неравномерности поля в рабочей зоне компактного полигона являются конструкция коллиматора, свойства радиопоглоща-ющего материала и его расположение в безэховой камере, а также переотражения волны со вспомогательными устройствами.
1. Краевые эффекты на зеркальных коллиматорах
Одной из важнейших задач при использовании коллиматоров является вопрос об уменьшении влияния краевых дифракционных эффектов на поле в рабочей зоне коллиматора. Попытки уменьшения такого влияния начали предприниматься уже с момента появления первых компактных полигонов [2]. Широко известны и повсеместно применяются компактные полигоны с рефлекторами с зазубренными краями. Технология производства таких коллиматоров постоянно совершенствовалась, и на сегодняшний день она позволяет снизить неравномерность поля в рабочей зоне до величин порядка 1 Дб.
Ввиду необходимости в проведении измерений с большей точностью, уже в 1980-е гг. стали активно разрабатываться коллиматоры, позволяющие проводить более точные измерения. Одним из примеров подобных коллиматоров может служить зеркало
22 ВМУ. Физика. Астрономия. № 6
со скругленным краем [3, 4]. В двумерном случае зеркало строится как совокупность параболической части зеркала и скруглений на границах, отклоняющих падающие на края лучи таким образом, чтобы они не попадали в рабочую зону. Примером подобных скруглений могут служить эллиптические скругления.
В работе [4] была предложена модель трехмерного варианта симметричного зеркала со скругленными краями и сложным контуром сопряжения параболической и скругленной частей. В результате исследований методами физической оптики авторами был предложен нетривиальный невыпуклый криволинейный контур сопряжения в виде четырехконечной звезды [5].
Для оценки поведения поля коллиматора с загнутыми краями в рабочей зоне был сделан целый ряд расчетов. В работе [6] можно найти результаты численного эксперимента, проведенного методами геометрической теории дифракции. Эти результаты достаточно хорошо согласуются с результатами, полученными методом физической оптики.
В работе [7] приведены результаты эксперимента, основанного на установке дополнительных рассеивающих пластин на уже функционирующий коллиматор, расположенный в университете штата Огайо, США. Сглаживающая «юбка» была специально изготовлена одним из исследовательских центров под эгидой NASA. Результатом подобной модификации, как и предполагалось, стало существенное сглаживание поля в рабочей зоне. Этот эксперимент наглядно показывает возможность модифицирования коллиматора с зазубренным краем при помощи подобных устройств, тем не менее на практике к такой возможности прибегают редко ввиду большой сложности в производстве «юбки» и переоборудовании камеры.
В качестве примера более подробного теоретического исследования на эту тему можно привести работы [8, 9]. В первой из них рассмотрено построение коллиматора с зазубренным краем специальной формы, угол наклона зубцов которого, в соответствии с методами геометрической оптики, выбран таким образом, чтобы перенаправлять отраженные от ребер лучи мимо рабочей зоны коллиматора. Несмотря на это, сравнение полученного коллиматора с коллиматором, подобным построенному в работе [4], показало, что скругленные края позволяют получить намного более гладкое поле в рабочей зоне. В настоящее время коллиматоры со скругленным краем широко применяются на практике [10].
Отметим, что эффективным способом улучшения поля в рабочей зоне также является использование высококачественных радиопоглащающих материалов. В настоящее время такие материалы применяются повсеместно для покрытия стенок, пола, потолка и приборов внутри безэховой камеры [11, 12]. На сегодняшний день существует достаточно боль-
шое количество радиопоглащающих материалов как узкополосного действия, основанных на явлениях резонанса, так и более универсальных и широкополосных [13, 14, 15, 16].
2. Математическое моделирование поля коллиматора
В работе методом математического моделирования исследуется коллиматор с отогнутой кромкой. Подобные задачи достаточно эффективно решаются с помощью метода конечных элементов [17]. В нашей работе мы используем также весьма эффективный метод интегральных уравнений. Рассмотрим модельную задачу. Пусть электромагнитная волна от линейного источника падает на идеально проводящую бесконечную цилиндрическую поверхность с сечением сложной формы (рис. 1), ось которой параллельна источнику. Здесь поверхность АВ, расположенная ближе к источнику, представляет собой поверхность пароболического цилиндра, фокус которого совпадает с источником, поверхность СО — кругового, а поверхности АС и ВО, соединящие первые две, — цилиндрические поверхности эллиптического типа. Потребуем, чтобы обращенная к источнику часть исследуемой поверхности представляла собой как минимум поверхность типа С1.
А В
Рис. 1. Сечение коллиматора с отогнутым краем
Система координат выбрана так, что ось г направлена вдоль оси цилиндра, поэтому задачу можно рассматривать как двумерную в плоскости ху. Здесь в качестве начала координат выберем фокус параболы, образованной пересечением параболического цилиндра и плоскости ху, ось у направим вдоль оси параболы, а источник поместим в начало координат. Уравнение эллипса (к примеру, АС) запишем как
х2 + 2Ь^ху + Ь22У2 + 2Ьпх + 2Ь23У + Ьзз = 0. (1)
Требуя, чтобы полученный в сечении контур был достаточно гладким в точках сопряжения кривых, запишем условия гладкого сопряжения эллипса АС и параболы АВ:
4й13 Ь12 + 14Ь22 + 8d21Ь13 +
+ 4Л2Ь23 + 4^2 Ь33 + 4^12 = 0, 4й12Ь12 + 13Ь22 + 2й1Ь23 + 4d2Ь13 + 4й21 = 0, (2) (1 + 21 )йЬ12 + 12Ь22 + 2Ь32 + 2d = 0.
Здесь I — х-координата точки сопряжения, d — фокусное расстояние параболы. С учетом этих условий
мы получим семейство эллипсов с двумя свободными параметрами.
Предположим, что камера достаточно велика по сравнению с исследуемым объектом. При такой постановке мы пренебрегаем влиянием стенок камеры и рассматриваем только поле, рассеянное зеркалом коллиматора. Дифрагированное поле можно представить в виде суперпозиции полей электрического и магнитного типа и записать с помощью функций Боргниса.
Для поля электрического типа рассмотрим уравнение Гельмгольца для функции Боргниса u [18] (здесь и далее предполагается зависимость от времени в виде etíjt):
Ди + k2u = f (M) с граничными условиями
и\с = 0
и условием излучения на бесконечности
(3)
(4)
\JU + tku]
tku | = o ( -1 |. ) r^tt VvV
(5)
с граничными контура
условиями Неймана на границе
Гды\
( дП ) +tku)
=o
Ш-
(7)
(8)
Для компонент электромагнитного поля в этом слу чае имеем
ди(х, у)
i du(x, у)
Ex Й !
Ш£ ду
Ez = Hx = Ну = 0,
Ey = —
Hz = u(x, у).
дх
(9) (10)
Введем функцию Грина g(M, Р) свободного пространства как решение уравнения
Ag(M, Р) + к^(М, Р) = -2ж6(гмр) (11)
с условиями излучения на бесконечности с решением в виде функции Ханкеля второго рода нулевого порядка
д(М, Р) = - Н(2)(кгмР). (12)
Для полей электрического и магнитного типа, воспользовавшись свойствами потенциалов простого и двойного слоя, получаем следующие соотношения:
J_
2п
Ou(P) дп
g(M, P) ds = -u0(M),
(13)
2 u(M)+¿
u(P)
дg(M, P) дп
ds = u0(M), (14)
с
Здесь С — контур, полученный путем гладкого сопряжения кривых, описанного выше (рис. 1), £ — волновое число. Задача для поля магнитного типа имеет тот же вид, с точностью до граничного условия. Поле магнитного типа удовлетворяет уравнению
Аи + £2и = / (М), (6)
где М — точка наблюдения на поверхности проводника, Р — точка истока, а функция и0(М) характеризует падающую волну. Для случая возбуждения поверхности нитью электрического тока уравнение для тока на поверхности проводника имеет вид
И(2\кгм) = -¿ф/(Р)щ*'(£\гм-гр|) йз. (15)
С
Зная распределение токов, можно вычислить поле в рабочей зоне коллиматора по формуле 1
С
Задача решается численно методом коллокации.
Перейдем к дискретному аналогу уравнения (15), разделив контур интегрирования на N дуг Сп. Представляя интеграл в формуле (15) в виде суммы интегралов по участкам Сп, С = ^ Сп, получим систему линейных уравнений относительно токов ]п
(2),
Hz(M) = ^ ф/(P)H(2(krP,M) dlP.
H(2) H0
N
(krm) = - 2nYlJ
n= 1
H(2)(k\rM -rP \) ds,
(16)
в которой коэффициентами матрицы будут интегралы от функции Ханкеля
M = - —
Mmn 2п
H(2)(k\rM-rP\) ds.
(17)
Так как дуги Сп выбираются малыми, функцию Хан-келя можно аппроксимировать кусочно-постоянной функцией
Мт,п = --1 АН0(2) (V(хт - Хп)2 + (ут - Уп)2) (18) при т = п. Для т = п
Mn,n = --1 Д (1 - i - ln 2п \ п
jkA
4e
,
(19)
где Ап — длины дуг контуров Сп, 7 — постоянная Эйлера, е — основание натурального логарифма. Таким образом, распределение токов на цилиндре произвольной формы мы можем получить как решение системы линейных уравнений
МтЛ = Ьт (20)
с симметричной матрицей. Тогда искомое поле имеет вид
и = — V /H(2)(kr ) .Д
^app — 2п > 0 ^ ln,mJ Дп-
(21)
На основании приведенной выше модели была составлена программа в среде МаШсаё и произведен ряд расчетов, в частности для нескольких характерных предельных случаев [19].
Для проверки корректности модели были произведены расчеты для простейшего случая цилиндра с круговым сечением, которое получается из контура, изображенного на рис. 1 при стремлении
Рис. 2. Распределение токов на круговом цилиндре и его поле в рабочей зоне
длины дуги AB к нулю. Расположим источник излучения на большом расстоянии. Токи, наведенные на параболоиде в этом случае, как и рассеянное поле, вычисляются аналитически [20]. Как показал численный эксперимент, результаты, полученные с помощью данной математической модели, хорошо согласуются с аналитическим решением (рис. 2).
Более важный с практической точки зрения предельный случай — случай, в котором фокусное расстояние параболы устремляется к бесконечности. В этом случае передняя и задняя стенка сечения г = const вырождаются в отрезки прямых. Так как источник излучения должен располагаться в фокусе параболической дуги, он также становится бесконечно удаленным, а падающая на тело волна — плоской. Был проведен ряд численных экспериментов с различными размерами скругления по отношению к длине волны. В качестве величины, характеризующей размер скругления, была выбрана длина большой полуоси эллипса a.
На графике, изображенном на рис. 3, разные кривые соответствуют разным размерам главной оси скругления a по отношению к длине волны Л. Хорошо видно, что увеличение размеров скругления положительно сказывается на поле в рабочей области, существенно подавляя неоднородности.
Одновременно с этим растет и поле вне рабочей зоны. Несмотря на то, что в рамках данной задачи это поле не рассматривается подробно, при реальных измерениях этот эффект может повлиять на точность эксперимента в безэховой камере из-за переотражений от стенок.
Наибольший интерес представляет исследование общего случая, в котором дуга AB не является вырожденной. В отличие от предыдущего случая, семейство возможных скруглений является многопараметрическим: подходящие эллипсы могут иметь различные углы наклона осей и координаты точки сопряжения с параболической частью контура. Тем не менее, общая тенденция остается неизменной: увеличение размеров закругления при прочих равных
Рис. 3. Поле, отраженное плоским зеркалом с разными размерами скругления
параметрах улучшает поле в рабочей зоне, но одновременно увеличивает и амплитуду поля вне ее.
Использование скругления достаточно больших размеров эффективно подавляет большую часть дифракционных эффектов, как хорошо видно из рис. 4, где представлены распределения поля для разных профилей зеркала, отличающихся размерами скруглений. Если на первом графике, соответствующем наименьшему размеру скругления, присутствуют дополнительные максимумы в рабочей зоне, то уже на втором графике эти экстремумы пропадают и в дальнейшем гладкость поля только улучшается.
Заключение
Построенная математическая модель для двумерного случая, несмотря на свою простоту, позволяет качественно оценить поле в рабочей зоне зеркального коллиматора с отогнутыми краями. Проведенные численные эксперименты показывают как преимущества, так и недостатки подобных коллиматоров: к преимуществам следует отнести существенное улучшение отраженного поля, но достигается оно при достаточно крупных размерах скругленных кромок.
Список литературы
1. Балабуха Н.П., Зубов А.С., Солосин В.С. Компактные полигоны для измерения характеристик рассеяния объектов. М., 2007.
2. Emerson W. // IEEE Trans. Antennas and Propagation. 1973. 21, N 4. P. 484.
3. Pistorius C., Burnside D. // IEEE Trans. Antennas and Propagation. 1987. 35, N 3. P. 342.
4. Gupta I., Ericsen K., Burnside W. // IEEE Trans. Antennas and Propagation. 1990. 38, N 6. P. 853.
5. Lee T., Burnside W., Gupta I. // Antenna Measurement Techniques Association 26th Ann. Mtg. Stone Mountain Park, Ga. 2004, 17-22 Oct.
6. Ellington S., Gupta I., Burnside W. // IEEE Trans. Antennas and Propagation. 1990. 38, N 12. P. 1969.
7. Burnside W., Gilreath M., Kent B., Clerici G. // IEEE Trans. Antennas and Propagation. 1987. 35, N 2. P. 176.
8. Munoz-Acevedo A., Burgos S., Sierra-Castaner M. // Proc. of the 5th European Conf. IEEE Antennas and Propagation (EUCAP). 2011. P. 3586.
9. Lee T., Burnside W. // IEEE Trans. Antennas and Propagation. 1996. 44, N 1. P. 87.
10. Shields M., Fenn A. // Lincoln Laboratory J. 2007. 16, N 2.
11. Appel-Hansen J. // IEEE Trans. Antennas and Propagation. 1973. 21, N 4. P. 490.
12. Никитенко А.В, Зубов А.С., Шапкина Н.Е. // Ма-тем. моделирование. 2014. 26, № 9. С. 18.
13. Severin H. // IRE Trans. Antennas and Propagation. 1956. 4, N 3. P. 385.
14. Sun W., Liu K., Balanis C. // IEEE Trans. Antennas and Propagation. 1996. 44, N 6. P. 798.
15. Amano M., Kotsuka Y. // IEEE Trans. Microwave Theory and Techniques. 2003. 51, N 1. P. 238.
16. Johansson M., Holloway C., Kuester E. // IEEE Trans. Antennas and Propagation. 2005. 53, N 2. P. 728.
17. Боголюбов А.Н., Боголюбов Н.А., Свешников А.Г. // Физические основы приборостроения. 2013. 2, № 1. С. 10.
18. Свешников А.Г., Могилевский И.Е. Математические задачи теории дифракции. М., 2010.
19. Хлебников Ф.Б., Шапкина Н.Е., Солосин В.С. // 12-я Ежегодная научная конференция ИТПЭ РАН. Москва, 4-7 апреля 2011.
20. Хлебников Ф.Б. // 5-я Международная конференция «Акустооптические и радиолокационные методы измерений и обработки информации». Суздаль, 18-19 сентября 2012.
The mathematical simulation of a diffracted field in the target area of a mirror collimator
F.B. Khlebnikov1a, A.N. Bogolyubov1, V.S. Solosin2, N.E. Shapkina1
1 Department of Mathematics, Faculty of Physics, Lomonosov Moscow State University, Moscow 119991, Russia.
2 Federal State Research Institute for Theoretical and Applied Electromagnetics, Russian Academy of Sciences (ITAE RAS), Moscow 125412, Russia.
E-mail: a mnfkhl@gmail.com.
A mathematical model of a mirror collimator with rounded edges is developed. The problem of calculating electromagnetic field components in the target area is solved by the method of integral equations. The results of numerical computations are presented.
Keywords: mathematical simulation, mirror collimator. PACS: 02.30.Rz. Received 14 July 2015.
English version: Moscow University Physics Bulletin 6(2015). Сведения об авторах
1. Хлебников Федор Борисович — аспирант; тел.: (499) 147-63-62, e-mail: mnfkhl@gmail.com.
2. Богобюбов Александр Николаевич — доктор физ.-мат. наук, профессор; тел.: (495) 939-10-33, e-mail: bogan7@yandex.ru.
3. Солосин Владимир Сергеевич — канд. физ.-мат. наук, вед. науч. сотрудник; тел.: (495) 485-95-44, e-mail: svs15105@yandex.ru.
4. Шапкина Наталья Евгеньевна — канд. физ.-мат. наук, доцент; e-mail: neshapkina@mail.ru.
