Теоретическое и экспериментальное исследование эффекта волнового обтекания цилиндрических объектов в СВЧ диапазоне

Теоретическое и экспериментальное исследование эффекта волнового обтекания цилиндрических объектов в СВЧ диапазоне

МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ

Учреждение образования

«Гомельский государственный университет

имени Франциска Скорины»

Дипломная работа

Теоретическое и экспериментальное исследование эффекта волнового обтекания цилиндрических объектов в СВЧ диапазоне

Реферат

Ключевые слова: маскировка, метод волнового обтекания, двухвитковая спираль, искусственная структура, металлический цилиндр, прямолинейный проводник.

Объект исследования: теоретическое и экспериментальное исследование искусственных структур, созданных на основе двухвитковых спиралей и прямолинейных проводниках, окружающие металлический цилиндр.

Предмет исследования: распределение электрических и магнитных полей в исследуемой структуре. Частотные зависимости интенсивности прошедшей волны в заданном диапазоне частот.

Методы исследования: компьютерное моделирование рассеяния плоской электромагнитной волны на искусственных структурах в программах ANSYS HFSS, Comsol Multiphysics и MathCAD. Методика изготовления искусственных образцов и проведение эксперимента.

Цель практической работы: теоретическое исследование эффекта волнового обтекания различных объектов СВЧ волнами. Исследование и моделирование задачи рассеяния плоской электромагнитной волны на любой искусственной структуре. Изготовление искусственных образцов. Ознакомление с принципом работы аппаратуры, используемой в эксперименте и методикой проведения эксперимента. Анализ полученных результатов.

Введение

Во все времена людей интересовали способы маскировки. Мечта о полной невидимости нашла воплощение в сказках многих народов, а также в фантастических романах и фильмах. Всевозможные плащи и шапки-невидимки, а также другие сказочные “устройства” придавали героям произведений невидимость и множество связанных с ней преимуществ. Мечты человечества о невидимости материальных объектов и попытки их воплотить существуют и в наши дни. В XX в. такие попытки только усилились с развитием радиолокации и звуколокации, а одним из примеров, показывающим практическую значимость невидимости, явилось создание технологии “Стелс”, которая основана на придании укрываемым объектам специальной формы или покрытии их специальной твёрдой оболочкой.

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

Технологии обнаружения в последние десятилетия не отставали от технологий маскировки: для каждого нового способа маскировки очень быстро находился свой способ обнаружения объекта по косвенным признакам. В связи с этим имеется острая потребность в такой маскировке, которая минимизировала бы все эти косвенные признаки. В последние годы было предложено несколько новых типов маскировки. Наибольшее развитие благодаря своей универсальности получила идея маскировки по методу волнового обтекания. Популярный обзор метода волнового обтекания приведён в [1].

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

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

Основной задачей данной работы является теоретическое и экспериментальное исследование эффекта волнового обтекания цилиндрических объектов СВЧ волнами. Моделирование рассеяния плоской электромагнитной волны на любой искусственной структуре. Снятие полученных результатов. Изготовление искусственных образцов. Проведение эксперимента и анализ полученных данных.

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

Подробно рассматривается первая экспериментальная демонстрация маскировки по методу волнового обтекания, так как именно этот первый успех и стал отправной точкой нового направления в технике маскировки.

Рассматривается строгое решение рассеяния плоской электромагнитной волны на идеально проводящем цилиндре. Решение данной задачи позволяет проверить достоверность решение этой же задачи с помощью программ ANSYS HFSS и Comsol Multiphysics.

Выполнен расчет параметров искусственной структуры для реализации возможности огибания цилиндрических объектов СВЧ волнами. На основании полученных результатов в моделировании, был изготовлен искусственный образец, на основе оптимизированных двухвитковых спиралей, для реализации метода волнового обтекания. Был проведен эксперимент для подтверждения полученных данных в моделировании.

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

Так же исследуется искусственная структура на основе не резонансных прямолинейных проводников. Данная структура смоделирована с помощью программы ANSYS HFSS и на основании полученных данных был проведен эксперимент.

1 Теоретическое исследование маскировки тел методом волнового обтекания

.1 Различные методы маскировки тел

Способы маскировки чрезвычайно интересовали людей во все времена. Мечта о полной невидимости нашла воплощение в сказках многих народов, а также в фантастических романах и фильмах. Всевозможные плащи и шапки-невидимки, а также другие сказочные “устройства” придавали героям произведений невидимость и множество связанных с ней преимуществ. Мечты человечества о невидимости материальных объектов и попытки их воплотить существуют и в наши дни.

Мечты человечества о невидимости материальных объектов и попытки их воплотить существуют и в наши дни. В XX в. такие попытки только усилились с развитием радиолокации и звуколокации, а одним из примеров, показывающим практическую значимость невидимости, явилось создание технологии “Стелс”, которая основана на придании укрываемым объектам специальной формы и/или покрытии их специальной твёрдой оболочкой.

Для маскировки, сокрытия или уменьшения заметности материальных тел придумано много различных методов. Популярный обзор этих методов приведён в [2].

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

Всем известный пример такой маскировки под цвет - камуфляжная одежда военных, охотников и др. У представителей животного мира различают сразу несколько видов маскировки под цвет: криптическая окраска под цвет окружающей среды (примеры: зелёные ящерицы, кузнечики, гусеницы в травянистой местности, жёлто-бурые саранча); дизруптивная окраска - маскировка на основе контрастной окраски, уменьшающая заметность в местности со сложным рельефом (зебры, бурундуки); скрадывающая окраска, основанная на эффекте светотени, при котором более освещенные участки тела животного (спина) окрашены в более тёмный цвет (встречается у многих животных). Ещё один способ достижения невидимости, заимствованный у природы, - маскировка не под цвет, а под свойства окружающей среды. Например, у медуз и некоторых рыб показатель преломления тела почти равен показателю преломления морской воды. Объект, обладающий таким свойством, практически невозможно заметить даже при самом внимательном разглядывании. К сожалению, этим свойством нельзя наделить произвольный объект, поэтому такой способ маскировки имеет очень ограниченную область применения.

Появление локационной техники значительно расширило понятие невидимости. Раньше под словом “невидимый” подразумевали нечто недоступное зрению, однако сейчас это понятие иногда применяют к объектам, которые не способен обнаружить радар. У противолокационной маскировки есть свои особенности, по сравнению с оптической. Эти особенности связаны, прежде всего, с тем, что локатор “видит” объекты в собственном отражённом излучении. Если излучение локатора будет полностью поглощено каким-то предметом, то это будет воспринято так же, как если бы оно ушло в пустоту, поскольку в обоих случаях отражённая волна отсутствует. Получается, что этот предмет невидим. Однако глаз, в отличие от радара, видит объекты в свете внешнего источника, поэтому объект с поглощающим покрытием хорошо видим почти на любом фоне, так как имеет чёрный цвет.

Для противолокационной маскировки применяют и рассеивающие покрытия. Например, технология “Стелс” основана на придании самолёту специфической угловатой формы, благодаря которой он рассеивает падающее радиоизлучение не назад (в направлении излучателя), а в каком-то совсем другом направлении (кстати, такая форма значительно ухудшает аэродинамику самолёта), тогда как в оптической области рассеивающее покрытие не придаёт объекту незаметность, а лишь размывает его очертания. Однако технологии маскировки, подобные “Стелс”, теряют эффективность, если локационная система состоит из нескольких радаров: отражённое от объекта излучение одного радара попадает в поле зрения другого. С помощью системы локаторов можно обнаружить и поглощающие объекты: каждый радар "видит" уже не только в собственном, но и во внешнем излучении, поэтому поглощающий излучение объект будет заметен радар у так же, как глазу заметен чёрный предмет на светлом фоне.

Технологии обнаружения в последние десятилетия не отставали от технологий маскировки: для каждого нового способа маскировки очень быстро находился свой способ обнаружения объекта по косвенным признакам (тени, рассеянному излучению или по путному следу). В связи с этим имеется острая потребность в такой маскировке, которая минимизировала бы все эти косвенные признаки. В последние годы было предложено несколько новых типов маскировки. Наибольшее развитие благодаря своей универсальности получила идея маскировки по методу волнового обтекания.

.2 Идея метода волнового обтекания

Считается, что кардинально новый метод маскировки предложили в 2006 г. сэром Джоном Пендри (Империал Колледж, Великобритания) с коллегами [3,4], известным специалистом в области электродинамики левых сред и метаматериалов, а его изобретение было основано на математических работах другого английского учёного - Ульфа Леонардта [5,6]. Предложенный ими метод маскировки был назван H.H. Розановым методом волнового обтекания [2], а в англоязычной литературе за этим методом закрепился термин cloaking.

Их основная идея состоит в том, что маскирующее покрытие должно искривлять фронт падающего электромагнитного излучения, заставляя лучи огибать объект, а на выходе принимать прежнее направление (рисунок 1). Тогда внешне всё выглядит так, как будто на пути волны ничего не было, т.е. наблюдатель вместо объекта, защищенного маскирующей оболочкой, увидит пустое пространство позади неё. Для того чтобы наблюдатель не заметил никаких признаков присутствия неоднородности, нужно соблюсти ещё одно условие: оптическая длина пути каждого луча в оболочке должна быть такой же, как если бы он распространялся прямолинейно в свободном пространстве. Иначе лучи, прошедшие через оболочку, просто не “впишутся” в общую картину: они будут интерферировать с лучами, не взаимодействовавшими с оболочкой, искажая картину поля.

Рисунок 1 - Траектории лучей в маскирующей оболочке

Осуществлять идею предлагается за счёт неоднородности вещества маскирующего покрытия. Действительно, явление преломления луча на границе двух сред всем хорошо знакомо. Но если показатель преломления среды изменяется непрерывно, то луч непрерывно преломляется и его траектория представляет собой гладкую кривую. Возникает вопрос: как же рассчитать параметры маскирующего покрытия, а именно тензоры диэлектрической и магнитной проницаемостей, чтобы получить нужные лучевые траектории?

В [7] предлагается не менее интересная, чем сама идея искривления лучей, процедура. Эта процедура заключается в том, что внутри оболочки мысленно создаётся пространство с искривлённой метрикой, для чего находится необходимое преобразование координат. В основе процедуры лежит тот факт, что уравнения Максвелла инвариантны по отношению к пространственным преобразованиям. Эту процедуру легче всего пояснить на следующем простом примере. Пусть нужно создать маскирующее покрытие в виде шарового слоя , внешней средой которого является вакуум (или воздух). Сделаем преобразование, переводящее шар радиусом  с центром в начале координат в этот шаровой слой

                                            (1)

Уравнения Максвелла

                                                      (2)

после таких координатных преобразований формально сохраняют свой вид

                                                    (3)

где  и  - компоненты поля в новых координатах,  и  - тензоры с диагональными компонентами

(4)

Поскольку вид уравнений одинаков, то и решения будут иметь одинаковый вид по отношению к своим системам координат. Это означает, что среда, в которой создано распределение параметров (4), будет искривлять прямой луч так же, как преобразование (1) искривляет прямую линию, пересекающую шар с радиусом  (рисунок 2).

Рисунок 2 - а) Прямая линия и вектор  в исходной (декартовой) системе координат; б) Декартова сетка, те же линия и вектор в новых координатах

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

Итак, маскирующее покрытие является анизотропной средой. Возникает вопрос: почему каждый луч преломляется в ней единственным образом? Где же обыкновенный и необыкновенный лучи? Ведь анизотропные среды, как правило, снимают поляризационное вырождение для волны. Авторы метода [7] объясняют такую сохранность поляризационного вырождения тем, что у маскирующего покрытия тензоры диэлектрической и магнитной проницаемостей равны между собой. Действительно, среда маскирующего покрытия - это преобразованный вакуум, который обладает тем же свойством (). Но вакуум преломляет луч единственным образом, так почему же маскирующая оболочка должна расщеплять его надвое? Помимо этих общих рассуждений, авторы предлагают вполне строгое доказательство отсутствия двулучепреломления в маскирующем покрытии. Определяя оператор  выражением

                                                                        (5)

и обозначая , они выводят дисперсионное соотношение для маскирующей оболочки в виде

                                                  (6)

Уравнение (6) - это уравнение четвёртого порядка относительно , однако оно имеет только два различных решения (с противоположными знаками). В средах с двулучепреломлением, напротив, все четыре решения дисперсионного уравнения различны.

Если рассчитывать параметры среды не для шара, а, например, для куба, то некоторые компоненты проницаемостей окажутся отрицательными. Это можно наглядно показать на рисунке 3.

Рисунок 3 - Схема падения лучей на маскирующую оболочку (светлая область) прямоугольного сечения (лучи падают слева направо)

Лучи, проходящие выше штрихпунктирной линии (которая делит прямоугольник пополам), преломляются на входе в оболочку: как обычно, падающий и преломлённый лучи находятся по разные стороны нормали к границе раздела. Лучи, проходящие ниже штрихпунктирной линии, преломляются отрицательно: падающий и преломлённый лучи оказываются по одну сторону от нормали, т.е. здесь уже необходима левая среда (понятие левой среды будет рассмотрено в следующем параграфе). В то же время нетрудно понять, что любой падающий на оболочку кругового сечения луч будет преломляться положительно. Отрицательное преломление также необходимо присутствует в маскирующих оболочках многих других форм. Прежде всего, это относится к невыпуклым оболочкам, например, в форме сердца (такие формы будут рассмотрены ниже). Но вместе с тем понятно, что если необходимо замаскировать любое, даже невыпуклое тело, то это можно сделать с помощью сферической оболочки достаточно большого размера. Но что действительно необходимо для маскировки, так это наличие анизотропных градиентных материалов с компонентами проницаемостей меньшими единицы, или в некоторых случаях - отрицательными. Таких сред в природе не существует, но можно создать искусственно материалы или среду, у которых диэлектрическая и магнитная проницаемости будут отрицательные. Такие материалы (среды) называются метаматериалами.

Понятно, что невозможно создать оболочку с параметрами (4). Обсуждение о достижимости идеальной невидимости было рассмотрено в статье [8]. Например, приведём один из аргументов против возможности “универсального” маскирующего покрытия, позволяющего скрывать объекты с произвольными электродинамическими характеристиками. Для электромагнитного поля это возможно, если объект идеально экранирован от внешнего излучения произвольного вида, т.е. если напряженности поля тождественно равны нулю в объёме объекта. Но тогда решением уравнений Максвелла, например, для монохроматического излучения и сред с конечными значениями и кусочно-гладким пространственным профилем элементов тензоров ε и μ служит тождественный нуль во всем пространстве, что противоречит наличию излучателя. Можно прийти к выводу, что для идеальной невидимости или физически несодержательной неаналитичности пространственной зависимости ε и/или μ, или также физически недопустимых их бесконечных значений. Конечно, и в реальных средах в рамках полных (нестационарных) уравнений Максвелла поле может быть равным нулю за границей раздела (в области, в которую излечение не успело проникнуть за время действия источника). Но эта граница не может быть стационарной, она движется со световой скоростью [9], так что излучение непрерывно действующего источника со временем заполняет всё пространство, затекая и в объём маскируемого объекта. Т.е. маскировка с сингулярностью невозможна!

В последнее время разработаны два варианта образования оболочек без сингулярностей [10], для которых вполне реализуем метод волнового обтекания на практике, поэтому проблема с сингулярностью оказывается решаемой! Так же стоит отметить, что возникали ситуации, в которых, казалось бы, неразрешимые проблемы, но всегда находилось нетривиальное решение, позволяющее выйти из тупика. Так будем же оптимистами!

1.3 Отрицательный показатель преломления среды на основе метаматериалов

Метаматериал - композиционный материал, свойства которого обусловлены не столько свойствами составляющих его элементов, сколько искусственно созданной периодической структурой из макроскопических элементов, обладающих произвольными размерами и формой [11]. В очень грубом приближении такие макроскопические элементы можно рассматривать как искусственно внесенные в исходный материал атомы чрезвычайно больших размеров. Искусственная периодическая структура модифицирует диэлектрическую и магнитную проницаемости исходного материала. Разработчик метаматериалов имеет возможность выбора различных свободных параметров (размеры структур, форма, постоянный и переменный период решетки элементов, образующих структуру). Одно из возможных свойств метаматериалов - отрицательный коэффициент преломления, который проявляется при одновременной отрицательности диэлектрической и магнитной проницаемостей материала.

Интерес к материалам с отрицательным коэффициентом преломления начинается с работы советского физика В.Г. Веселаго, которая была опубликована в журнале “Успехи физических наук” за 1967 г. [12]. В работе было указано на возможность существования материала с отрицательным коэффициентом преломления, который был назван “левосторонним”. Материал с положительным коэффициентом преломления по аналогии следует назвать “правосторонним”. Автор пришел к заключению, что в левостороннем материале почти все известные оптические явления распространения волн существенно изменяются, хотя в то время материалы с отрицательным коэффициентом преломления еще не были известны.

В своей основополагающей работе В.Г. Веселаго показал, что рефракция - отклонение электромагнитной волны при прохождении границы раздела двух сред - изменяется в материалах с отрицательным коэффициентом преломления. В условиях, когда оба материала имеют одинаковый знак коэффициента преломления, волна, пересекая границу раздела, появляется на противоположной стороне линии, проходящей перпендикулярно к этой границе (нормаль к поверхности). Однако если один материал имеет положительный коэффициент преломления, а другой - отрицательный, волна будет появляться на той же стороне нормали, с которой она подходила к границе раздела.

Также особым свойством метаматериалов является характерная дисперсия волн, распространяющихся в такой среде. Уравнение распространения электромагнитных волн в изотропной среде имеет вид

                                                                            (7)

                                                                                        (8)

Решающим фактором, определяющим дисперсию электромагнитной волны, является коэффициент преломления. Так называемые ”левосторонние“ среды имеют отрицательный коэффициент преломления.

Уравнение (7) получено на основе теории Максвелла. Для сред, у которых диэлектрическая ε и магнитная µ проницаемости одновременно положительные, три вектора электромагнитного поля - электрический , магнитный  и волновой  - образуют так называемую правую систему векторов

                                                               (9)

Такие среды соответственно называют “правыми”. Среды, у которых ε и µ одновременно отрицательные, электрический , магнитный  и волновой вектор  образуют левую систему векторов. Такие среды называют “левосторонними”. В англоязычной литературе описанные материалы называют right- и left-handed materials, или сокращенно RHM (правые) и LHM (левые) материалы соответственно. Для характеристики распространения волны и ее дисперсии необходимо дать определение фазовой и групповой скоростей волны. В случае плоской гармонической волны фазовая скорость вдоль волнового вектора есть скорость движения поверхности равных фаз и ее можно выразить следующим образом

                                                                                        (10)

Подставляя сюда выражение для ω из (8), получаем на первый взгляд очевидное соотношение

                                                                                         (11)

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

Групповая скорость - это параметр, характеризующий скорость распространения “группы волн”, т. е. распространения более или менее хорошо локализованной квазимонохроматической волны Групповая скорость обычно интерпретируется как скорость перемещения максимума амплитудой огибающей квазимонохроматического волнового пакета. В одномерном случае групповая скорость вычисляется из закона дисперсии

                                                                                      (12)

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

На рисунке 4 показаны две искусственные линии передачи, собранные из сосредоточенных конденсаторов и индуктивных элементов. Для двух этих искусственных линий передачи получаем разные законы дисперсии.

Рисунок 4 - Искусственные линии передачи, собранные из сосредоточенных конденсаторов и индуктивных элементов: а) правосторонняя линия; б) левосторонняя линия

Для правосторонней линии имеем

                                                                                   (13)

где  - погонная индуктивность проводников линии,  - погонная емкость между проводниками линии.

Для левосторонней линии имеем

                                                                       (14)

где  - погонная обратная индуктивность, включенная между проводниками линии передачи,  - погонная обратная емкость, включенная последовательно в проводники линии передачи.

Погонная обратная индуктивность определяется следующим образом:

                                                       (15)

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

Более 30 лет прошло после публикации статьи В.Г. Веселаго, пока левосторонние материалы (ЛСМ) не были изобретены и экспериментально продемонстрированы. Эти ЛСМ не были природными веществами, как предполагал В.Г. Веселаго, а были искусственными структурами, эффективно однородными. Эти структуры были продемонстрированы Д.Р. Смитом с сотрудниками из калифорнийского университета в Сан-Диего [13]. В основу была положена работа Дж.В. Пендри [14-16] из лондонского Империал колледжа. Дж.В. Пендри предложил структуру плазмонного типа, имеющую отрицательную диэлектрическую проницаемость на частотах, соответствующих СВЧ-диапазону. Структура представляла собой сетку тонких металлических проволочек. Если электрическое поле параллельно осям проволочек, то в проволочках возбуждается ток, создающий эквивалентные электрические дипольные моменты. Усреднение дипольных моментов образует диэлектрическую проницаемость структуры с частотной зависимостью плазмонного типа (рисунок 5(a))

Рисунок 5 - Структура (a) - сетка тонких металлических проволочек. Если электрическое поле параллельно осям проволочек, то в проволочках возбуждаются электрические дипольные моменты. Структура (б) - кольцевые разомкнутые резонаторы. На частотах выше резонансной частоты эффективная магнитная проницаемость такой структуры отрицательна

Эффективная диэлектрическая проницаемость плазмонного типа равна

                                                                      (16)

где  - частота электронной плазмы, ζ - параметр затухания. Частота электронной плазмы и параметр затухания рассчитываются как

                                           (17)

где p - расстояние между проволочками, a - радиус проводника проволочки, c - скорость света, σ - удельная электропроводность материала проводника. Проводники расположены на расстоянии , где λg - длина волны, распространяющейся в структуре, что позволяет считать структуру электродинамически однородной. Структура кольцевых разомкнутых резонаторов (КРР) обладает резонансом (рисунок 5(б)), причем на частотах выше резонансной эффективная магнитная проницаемость такой структуры отрицательна

                                                              (18)

где  - резонансная частота кольцевого резонатора, ζ - параметр затухания. Резонансная частота кольцевого резонатора равна

                                (19)

где a - внутренний радиус меньшего кольца, δ - радиальное расстояние между кольцами, R1 - погонное сопротивление металлического проводника.

Д.Р. Смит с сотрудниками сформировали комбинацию структур [17] Дж.В. Пендри с отрицательной диэлектрической и магнитной проницаемостями (рисунок 6).

Рисунок 6 - Комбинация структур с отрицательной диэлектрической и отрицательной магнитной проницаемостями (экспериментальная реализация)

Образованная структура обладает отрицательным показателем преломления и СВЧ диапазоне.

.4 Первая экспериментальная демонстрация маскировки по методу волнового обтекания

О первом экспериментальном осуществлении маскировки методом волнового обтекания объектов оболочками из метаматериалов следует рассказать отдельно. Успех этих экспериментов во многом обусловил ту волну интереса, которую вызвал в научном мире этот метод маскировки. Он же выявил и основные проблемы, над решением которых специалисты работают в последние годы. Экспериментам предшествовало исследование [18], которое включало в себя численное моделирование обтекания в бесконечно длинной цилиндрической оболочке кругового сечения. Взаимодействие оболочки с линейно поляризованной электромагнитной волной моделировалось методом FDTD (Finite Difference Time Domain), не использующим приближение геометрической оптики и позволяющим рассчитать волновую картину полностью (full-wave simulation). Параметры цилиндрического покрытия были рассчитаны по методу координатных преобразований:

       (20)

в котором использовано линейное преобразование

                                                  (21)

где а и b - внутренний и внешний радиусы цилиндрической оболочки соответственно. Картина взаимодействия линейно поляризованной волны, вектор  которой параллелен оси цилиндра, с таким маскирующим покрытием показана на рисунке 7(а). Эффект маскировки налицо: излучение не проникает внутрь полости и практически не рассеивается на поверхности оболочки, фронт волны после прохождения через оболочку снова становится плоским.

Расчёты в работе [18] являлись подготовкой к экспериментам, поэтому её авторы не остановились на идеальной маскирующей оболочке. Понимая, что в действительности невозможно создать оболочку с непрерывным распределением параметров и без диссипативных потерь, они смоделировали также взаимодействие плоской волны со слоистым покрытием (восемь слоев, каждый из которых имеет свой показатель преломления) и покрытием с потерями.

Рисунок 7 - Распределение нормированной амплитуды электрического поля вблизи цилиндрической маскирующей оболочки. Тонкие горизонтальные линии показывают траектории лучей, направление потока энергии поля во всех случаях слева направо: а) покрытие с идеальными параметрами (20); б) покрытие с идеальными параметрами (20) и тангенсом угла потерь, равным 0,1; в) восьмислойное приближение идеального покрытия; г) покрытие с упрощёнными параметрами (23)

Из рисунка 7(в), в видно, что дискретность оболочки несколько ухудшает картину (траектории лучей вне оболочки перестают быть прямолинейными из-за рассеяния), а ненулевое поглощение приводит к появлению тени. К сожалению, пока очень трудно воплотить в жизнь покрытие с параметрами (20), так как все компоненты проницаемостей пространственно неоднородны. Однако если, например, вектор  падающей волны параллелен оси цилиндра оz, то картина становится двумерной и z-компоненты проницаемостей можно положить.

Тогда уравнения Максвелла внутри оболочки примут следующий вид:

(22)

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

                                      (23)

Несомненно, это сильно упрощает дело, поскольку здесь только одна компонента пространственно неоднородна. На рисунке 7(г) показана маскирующая оболочка с параметрами (23). По искривлённым траекториям лучей вне оболочки видно, что рассеяние во всех направлениях весьма значительно. Тем не менее эффект маскировки имеет место, так как излучение все же не проникает в область .

Для экспериментов была создана оболочка, обладающая всеми вышеперечисленными отклонениями от идеала: упрощёнными материальными параметрами (23), дискретностью, поглощением. Кроме того, цилиндр имел конечную (и очень небольшую) высоту. Покрытие содержало десять слоев, каждый из которых, в свою очередь, состоял из множества кольцевых резонаторов (наиболее распространённый тип структурных элементов метаматериалов). Резонаторы каждого слоя имели свои геометрические размеры, которые были подобраны так, чтобы достичь зависимости , максимально приближенной к (23). Радиусы слоев выбирались таким образом, чтобы на их окружностях укладывалось по целому числу резонаторов. Оболочка и параметры искусственной среды показаны на рисунке 8.

Рисунок 8 - Цилиндрическая маскирующая оболочка (фоновое изображение) с графиками  (наклонная линия),  (нижняя горизонтальная прямая) и  (верхняя горизонтальная прямая)

В качестве маскируемого объекта взят полый проводящий цилиндр радиусом м. Из рисунка 9, на котором представлены результаты экспериментов, видно, что, несмотря на все свои несовершенства, цилиндрическая оболочка всё же даёт частичный маскирующий эффект, приводящий к уменьшению заметности маскируемого объекта: проводящий цилиндр, покрытый этой оболочкой, значительно меньше рассеивает излучение, чем тот же цилиндр без оболочки (рисунок 9).

Рисунок 9 - Распределение нормированной амплитуды электрического поля: A) Покрытие с идеальными параметрами (20), (тонкие тёмные кривые показывают направление потока энергии поля); B) Покрытие с упрощёнными параметрами (23); C) Экспериментальные измерения с цилиндром без покрытия; D) Экспериментальные измерения для замаскированного цилиндра

Именно этот первый успех и стал отправной точкой нового направления в технике маскировки.

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

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

.5 Разнообразие форм маскирующих покрытий

Эксперименты с цилиндрической оболочкой вызвали появление целой лавины публикаций, посвященных расчётам маскирующих покрытий других форм. Надо сказать, что во всех появившихся работах приведены результаты лишь численных результатов (моделирование). Действительно, технологии на тот момент не позволяли создавать достаточно качественные оболочки. Из-за несовершенства технологий эксперименты (кроме первого, который доказывал реальность эффекта маскировки) теряли смысл, поскольку картина получалась бы слишком размытой для того, чтобы делать выводы о свойствах, обусловленных формой покрытия.

Прежде всего, заметим, что и цилиндрическая оболочка может не хуже сферической скрывать любой объект, который можно в неё поместить. Однако для большинства практических приложений цилиндр неудобен, например, сферическая оболочка в этом смысле более универсальна. В самом деле, при прямом падении волны на цилиндрическую оболочку (фронт волны параллелен образующей цилиндра) поток энергии поля вдоль образующей отсутствует и в любом поперечном сечении картина поля одинакова, т.е. она является плоской (поэтому такие покрытия называют 2В-цилиндрическими). Тем не менее, такие условно плоские задачи значительно менее громоздки, чем полноценные трёхмерные. Собственно, благодаря этой относительной простоте цилиндрические оболочки и заслужили столько внимания. Итак, перечислим основные исследованные формы оболочек.

Эллиптический цилиндр

Эллиптическая форма удобна своей универсальностью: в зависимости от отношения полуосей эллипс может быть очень вытянутым или почти круглым и тем самым удовлетворять разнообразным практическим задачам. На это обращается внимание в [20], где содержится вывод тензоров  и  для концентрической эллиптической оболочки.

Приведем для справки выражения, полученные в [20]

Здесь к - отношение большой и малой полуосей, одинаковое для внутреннего и внешнего эллипсов, а и b - их малые полуоси,  и . При этом использовалось нелинейное преобразование:

                                            (25)

В результате такая оболочка дает симметричный эффект. На рисунке 10 показано распределение амплитуды электрического поля для различных углов падения излучения.

Рисунок 10 - Распределение нормированной амплитуды электрического поля для различных углов падения излучения (по отношению к горизонтальной оси) на эллиптическую маскирующую оболочку: а) 0°; б) 90°; в) 45°

Произвольная форма

Параметры оболочки любой формы также можно рассчитать методом координатных преобразований. В общем случае не существует преобразований, переводящее односвязную область в подобную ей двусвязную (например, если граница области является ломаной линией). Тогда  и  задают кусочно, делая сразу несколько преобразований.

Метод неравномерных рациональных бисплайнов (Non-Uniform Rational В-Spline - NURBS), используемый в [21], в этом смысле является более точным. В этом методе гладкая граница маскирующей оболочки аппроксимируется несколькими кривыми Безье второго порядка. Эти кривые могут представлять собой любые конические сечения (эллипсы, параболы, гиперболы), в зависимости от параметров. Например, форма сердца описывается всего двумя кривыми Безье, соответствующими его левой и правой половинам (рис. 11).

Рисунок 11 - Распределение нормированной амплитуды электрического поля при падении излучения на маскирующую оболочку в форме сердца; стрелками показаны направления потока энергии

Ещё один интересный метод расчёта цилиндрических маскирующих оболочек с произвольным сечением изложен в работе [22]. Предложенный подход интересен тем, что он позволяет избежать подобия внутренней и внешней поверхностей оболочки. На рисунке 12 показана маскировка одной из оболочек, параметры которой получены этим методом.

Рисунок 12 - Распределение нормированной амплитуды электрического поля для цилиндрической маскирующей оболочки с произвольной формой внешней и внутренней границ

Такая маскировка имеет важные особенности: во-первых, маскируемый объект не является выпуклым, во-вторых, детали его геометрии имеют размер менее длины волны.

2. Исследование рассеяния плоской электромагнитной волны о металлический цилиндр

.1 Теоретический расчет рассеяния плоской электромагнитной волны о бесконечный круглый цилиндр

Предположим, что на бесконечный круглый цилиндр радиуса а падает плоская электромагнитная волна, распространяющаяся в направлении, перпендикулярном оси oz (рисунок 13). Требуется определить вторичное электромагнитное поле, возникающее в пространстве.

Рисунок 13 - Дифракции плоской волны на бесконечном круглом цилиндре в заданной системе координат

При решении поставленной задачи целесообразно рассматривать два случая поляризации падающей волны относительно оси oz: вектор  параллелен Oz и вектор  параллелен oz. Общий случай произвольной поляризации падающей волны может быть получен наложением этих типов полей. В соответствии с исходным условием проекции векторов поля падающей волны будут иметь вид:

      (26)

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

Для решения этой задачи целесообразно использовать цилиндрическую систему координат ρ, φ, z. Так как в рассматриваемом случае векторы  и  падающей волны не зависят от координаты z, а параметры цилиндра одинаковы на всем интервале , вторичное электромагнитное поле также не должно зависеть от z. Система уравнений для комплексных амплитуд проекций векторов вторичного поля поперечно-электрического типа в установившемся режиме имеет вид:

         (27)

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

                              (28)

где  - проекция вектора  падающей волны, касательная к поверхности цилиндра, а  - аналогичная проекция вектора  вторичного поля. Для решения системы уравнений (27) выразим из первого и второго уравнений неизвестные  и  через  и подставим найденные значения в последнее уравнение. В результате получим

                                                 (29)

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

                  (30)

где  - функция Бесселя,  - функция Ханкеля второго рода, m - произвольная постоянная,  - амплитуда падающей волны.

На рисунке 14 изображены нормированные кривые изменения амплитуд напряженности  вторичного поля () в зависимости от угла φ, построены при различных значениях параметра .

Рисунок 14 - Диаграмма направленности вторичного поля при дифракции плоской электромагнитной волны на проводящем цилиндре

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

2.2 Описание программного обеспечения и граничных условий модели

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

Используется пакет программ Ansoft HFSS, предназначенный для расчета параметров электромагнитных полей в сложных СВЧ устройствах. Данные программы включают черчение трехмерных объектов, процесс расчета, включающую несколько методов решения граничных задач и пост-процессор для обработки и детального анализа полученных результатов. Электродинамическое моделирование в Ansoft HFSS основано на использовании метода конечных элементов (Finite Element Method, FEM). Использование метода конечных элементов обеспечивает высокую степень универсальности численных алгоритмов, которые оказываются весьма эффективными для широкого круга задач от анализа волноводных и полосковых структур до моделирования антенн и сложных невзаимных устройств, содержащих гиротропные среды [24].

Рассмотрим идеализированную физическую модель в HFSS (рисунок 15). Она представляет собой идеально проводящий цилиндр, который находится в центре воздушного параллелепипеда (бокса). Бокс служит границами пространства, где налагаются граничных условий (ГУ). На грани бокса ставят такие граничные условия как идеально поглощающие слой (PML) и бесконечные периодические ГУ (Master и Slave).(Perfect Material Layer) - слой идеального поглощения, состоящий из несколько слоев специальных материалов, которые поглощают излучаемые волны. Идеально согласованные слои PML являются виртуальными материальными объектами, которые полностью поглощают электромагнитные поля, падающие на них. Эти материалы анизотропны. Слой PML смоделирован так, что с учетом отсутствия отражения направленных волн, структура продлевается, постепенно увеличивая толщину отдельных слоев к бесконечности. Поверхности, на которые она нагружена, излучают в направлении, в котором распространяется волна. На рисунке 15 слои PML, изображены голубым цветом, окружают боковые грани бокса.

Рисунок 15 - Вид модели в HFSS

Master и Slave - ведущие и ведомые границы, которые дают возможность моделировать периодические плоскости, где E-поле на одной поверхности соответствует E-полю на другой плоскости с заданной разностью фаз. Они фиксировано устанавливают E-поле в каждой точке на границе Slave, соответствуют E-полю с разностью фаз в каждой соответствующей точке на границе Master. Они накладываются на верхнюю и нижнюю грани бокса. В отличие от границ симметрии, на этих границах поле E не является тангенциальным или нормальным. Единственное условие состоит в том, что поля на этих двух границах должны иметь равные величины и направление (или равные величины и встречные направления).

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

                                                                                  (31)

Задержка по фазе рассчитывается по углам сканирования в сферической системе координат (φ и θ).

Граничные условия Master и Slave полезны для моделирования устройств типа бесконечных антенных решеток.

.3 Тестовая задача для проверки работоспособности используемой программы

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

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

Для примера решим задачу, которая была разобрана в параграфе 2.1. Постановка задачи: на идеально проводящий бесконечный цилиндр радиусом 0,0995 метра падает плоская электромагнитная волна, вектор напряженности электрического поля параллелен оси цилиндра, с длиной волны 0,1 метр. Нужно построить диаграмму направленности рассеянного поля, используя теоретический расчет в параграфе 2.1 и сравнить полученный результат с пакетом программ ANSYS HFSS.

Поставленная задача была смоделирована и решена. На рисунке 16 представлены диаграммы направленности (диаграммы рассеяния) электрического поля, рассчитанные теоретически, используя программу MathCAD и HFSS.

Рисунок 16 - Диаграмма рассеяния при дифракции плоской волны на проводящем цилиндре: а) теоретический расчет в MathCAD, б) моделирование в HFSS

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

Рисунок 17 - Распределение напряженности электрического поля (ANSYS HFSS)

Для графического сравнения распределения электрического поля можно использовать программу COMSOL Multiphysics. Эта программа использует решение уравнений в частных производных. В данной программе смоделирован бесконечный цилиндр, на который падает плоская электромагнитная волна. Физическая модель является идеализированной и полностью соответствует поставленной задаче, решенной в параграфе 2.1, и численно рассчитанной с помощью MathCAD. Распределение электрического поля показано на рисунке 18.

Рисунок 18 - Распределение напряженности электрического поля (COMSOL Multiphysics)

Следовательно, можно сделать вывод о том, что программа ANSYS HFSS корректно рассчитывает поля с высокой точность.

. Разработка искусственной структуры на основе двухвитковых спиралей для реализации возможности огибания цилиндрических объектов СВЧ волнами

3.1 Исследование возможности использования анизотропных параметров среды для практического применения метода волнового обтекания

В главе 1 были рассмотрены различные виды и формы скрывающих покрытий. Создать среду с параметрами (20) в настоящее время невозможно, поэтому рассматриваются более простые преобразования, которые возможно реализовать практически. Параметры среды задаются тензорами диэлектрической и магнитной проницаемости, поэтому необходимо найти компоненты тензоров для получения явления огибания электромагнитной волной цилиндрических объектов. На данный момент, наиболее простыми для практического использования являются компоненты тензоров, представленные в статье [25], которые даны в цилиндрической системе координат и имеют вид:

,                (32)

где а и b - внутренний и внешний радиусы цилиндрической оболочки, r расстояние между внутренним и внешним радиусами оболочки .

Для наглядного представления явления огибания в COMSOL Multiphysics смоделирован бесконечный цилиндр, окруженный оболочкой, в которой заданы компоненты тензоров (32). На цилиндр падает плоская электромагнитная волна с длиной волны 0,1 м., ось цилиндра направлена вдоль оси z.

Рисунок 19 - Цилиндр, окруженный искусственной оболочкой при а=0,1 м; b=0,2 м.

На рисунке 19 показаны, для сравнения, цилиндр без оболочки и цилиндр с маскирующей оболочкой, в которой заданы преобразования (32), при падении плоской электромагнитной волны.

На практике достаточно трудно осуществить, чтобы показатель преломления в оболочке изменялся непрерывно, поэтому его заменяют дискретной структурой. Подробный анализ дискретной оболочки рассмотрен в [26,27]. На рисунке 20 показано распределение напряженности электрического поля для дискретной оболочки с различным количеством слоев.

Рисунок 20 - Распределение напряженности электрического поля для дискретной оболочки: а) 2 слоя; б) 10 слоев

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

.2 Расчет параметров искусственной структуры на основе двухвитковых спиралей

Для реализации параметров среды (32) можно использовать различные метаматериалы. Первая реализация искусственной среды с параметрами (32) была рассмотрена в статье [28]. В качестве метаметериала в статье использовали каноническую спираль (рисунок 21(а)). Предполагалось, что полученная искусственная среда на канонических спиралях будет скрывать объект сразу для двух поляризаций, так как она имеет одинаковый отклик для горизонтальной и вертикальной поляризаций падающей волны.

Вертикальной поляризацией будем называть поляризацию падающей волны (Е-волна), у которой вектор напряженности электрического поля  параллелен оси исследуемого образца (ось z), соответственно, горизонтальная поляризация - поляризация падающей волны (H-волна), перпендикулярная оси образца. Однако стоит заметить, что в статье [28] показан результат только для одной поляризации падающей волны.

Для создания искусственной среды, которая работала бы сразу при двух поляризациях падающей волны, можно использовать обычную гладкую спираль (рисунок 21(б)). Преимуществом гладкой спирали по сравнению с канонической спиралью является простота изготовления и реализация в любом диапазоне частот. Для создания искусственной среды с параметрами (32), на основе гладких спиралей, необходимо использовать правосторонние и левосторонние спирали, для компенсирования (уничтожения) киральности в среде.

Рисунок 21 - Геометрия спиралей: а) каноническая спираль; б) гладкая правосторонняя двухвитковая спираль

Для получения искусственной среды с параметрами (32), необходимо знать диэлектрическую и магнитную проницаемости, которые могут быть выражены через диэлектрическую и магнитную восприимчивости среды.

Как известно, каждая спираль характеризуется одновременно диэлектрической, магнитной и киральной восприимчивостью. Метод расчета всех компонент тензоров восприимчивостей для любых спиралей с заданными параметрами был рассмотрен в [29,30]. Расчет компонент тензоров производился в программе ANSYS HFSS. Как указывалось выше, программа предназначена для расчета параметров электромагнитных полей в сложных СВЧ устройствах, поэтому идеально подходит для поставленной задачи.

Формулы для расчета необходимых компонент тензоров восприимчивостей спирали приведены ниже:

) Ex≠0, Ну≠0 - эти компоненты векторов определяют поляризацию падающей плоской электромагнитной волны.

, , ,

, ,                                     (33)

) Ey≠0, Hx≠0:

, , ,   (34)

) Ez≠0, Hx≠0:

, , ,

, ,                          (35)

) Ex≠0, Hz≠0:

, , , (36)

Здесь ,  - напряженность электрического поля падающей электромагнитной волны, f - частота, на которой проводятся измерения,  - напряженность электрического поля волны, излучаемой спиралью (нижний индекс обозначает направление, вдоль которого ведется измерение, а верхний - проекцию вектора на данное направление).

Для создания искусственной среды будем использовать двухвитковые спирали, так как они обладают преимуществом, по сравнению с одновитковыми спиралями. В статье [31] показан теоретический расчет параметров идеально проводящей двухвитковой спирали, для которой выполняется равенство диэлектрической, магнитной и киральной восприимчивостей. Такая спираль обладает одинаковым откликом для любой поляризации падающей волны. В статье [27] предоставлены данные осевых компонент тензоров двухвитковой спирали с параметрами как в статье [31], у которых расхождение в значениях компонент составляет порядка 10%. Этот результат можно улучшить, используя метод расчета компонент тензоров восприимчивостей, можно оптимизировать параметры двухвитковой спирали. Оптимизированные параметры двухвитковой спирали являются следующими:

Двухвитковая спираль: NВ=2, α=6,35°; L=0,05 м; r=3,3∙10-3 м; h=2,3∙10-3 м; d=0,5∙10-3 м,

где NВ - число витков спирали, α - угол подъема спирали относительно плоскости, перпендикулярной оси спирали, L - длина проволоки, из которой изготовлена спираль, r - радиус витка, h - шаг спирали, d - диаметр проволоки. Эти оптимальные параметры были использованы при компьютерном моделировании электромагнитных свойств спиралей.

На рисунке 22 показаны осевые компоненты тензоров восприимчивости оптимизированной двухвитковой спирали.

Рисунок 22 - График частотной зависимости действительных и мнимых частей осевых компонент тензоров восприимчивостей двухвитковой спирали

Полученные данные рассматриваются в дальней зоне в диапазоне частот от 3 до 3,25 ГГц, где наблюдается главный частотный резонанс токов. Из графика видно, что различия в значениях осевых компонент тензоров составляет не более 1%.

Как известно, чтобы использовать формулу (32) необходимо рассчитать компоненты тензоров диэлектрической и магнитной восприимчивостей, а через них диэлектрическую и магнитную проницаемость среды. Относительная диэлектрическая и магнитная проницаемость среды может быть рассчитана, используя формулу Клаузиуса-Моссоти [32]:

                                 (33)

где  и  - диэлектрическая и магнитная восприимчивости среды, n объемная концентрация частиц в среде.

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

Исходя из данных на рисунке 22 и используя формулы (32) и (33) можно рассчитать необходимую объемную концентрацию частиц (спиралей) в среде.

                                                    (34)

Количество частиц в среде для цилиндрического слоя рассчитывается по формуле как

                   (35)

где V - объем слоя, S - площадь слоя, h - высота слоя.

Для расчета необходимого количества спиралей можно использовать формулу (35), при этом задавая произвольно радиус цилиндра и радиус маскирующей оболочки, но параметр r должен изменяться непрерывно. Этого очень трудно добиться, поэтому делают дискретную оболочку, где параметр r разбивается (дискретизируется), в пределе между радиусом оболочки и радиусом цилиндра (), на конечное число, которое соответствует числу слоев в оболочке.

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

                                        (36)

где  - толщина первого слоя, r1 - расстояние от центра до середины первого слоя. Аналогично рассчитывается количество спиралей в остальных слоях. На рисунке 23 показана схема дискретной оболочки, состоящая из четырех слоев, окружаемая скрываемый цилиндр.

Рисунок 23 - Схема дискретной оболочки

Параметр r берется равным от центра всей структуры до середины каждого слоя, так как именно на этом расстоянии будут располагаться спирали, которые и сформируют слои оболочки.

Расчет искусственной структуры проводится при следующих параметрах:

радиус цилиндра a=7,75 см;

радиус оболочки b=15,5 см;

толщина структуры h=1 см;

восприимчивость двухвитковой спирали на частоте 3,14 ГГц равна минус 13,8 см3.

Значение действительной части восприимчивостей двухвитковой спирали выбирается максимально отрицательной, для того, чтобы диэлектрическая восприимчивость среды изменялась в пределах от единицы до нуля (см. формулу 32).

Основываясь на теоретических расчетах и параметрах структуры, были рассчитаны диэлектрическая проницаемость среды и количество двухвитковых спиралей в каждом слое. Данные приведены в таблице 1.

Таблица 1 - Параметры дискретной оболочки

Схема расстановки спиралей в образце выбрана с учетом симметрии относительно цилиндра и изображена на рисунке 24.

Рисунок 24 − Схема расстановки спиралей в образце

Основываясь на рассчитанных параметрах, было проведено моделирование бесконечной многослойной искусственной структуры с использованием программы ANSYS HFSS. Результаты моделирования представлены на рисунке 25, в виде распределения напряженности электрического поля для вертикальной поляризации падающей волны.

Рисунок 25 − Распределение напряженности электрического поля на частоте 3,01 ГГц для вертикальной поляризации: а) металлический цилиндр без оболочки; б) металлический цилиндр, окруженный искусственной оболочкой

Из рисунка 25 видно, что фронт волны частично восстанавливается для вертикальной поляризации на 3,01 ГГц, но при этом присутствует достаточно большое рассеяние на структуре.

На рисунке 26 показано распределения напряженности магнитного поля для горизонтальной поляризации падающей волны.

Рисунок 26 − Распределение напряженности магнитного поля на частоте 2,97 ГГц для горизонтальной поляризации: а) металлический цилиндр без оболочки; б) металлический цилиндр, окруженный искусственной оболочкой

Из рисунка 26 следует, что напряженность магнитного поля за цилиндром с искусственной оболочкой лишь немного превышает амплитуду магнитного поля за цилиндром без оболочки.

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

                                                                               (37)

где  - комплексная величина напряженности электрического поля;  комплексно-сопряженная величина напряженности электрического поля; S - площадь плоскости, созданной за искусственной структурой.

На рисунке 27 показана частотная зависимость интенсивности прошедшей волны относительно падающей при различной поляризации падающей волны.

Рисунок 27 − Частотная зависимость интенсивности прошедшей волны относительно падающей при горизонтальной и вертикальной поляризации

Из рисунка 27 следует, что интенсивность прошедшей волны при вертикальной поляризации на частоте 3,01 ГГц за цилиндром с искусственной оболочкой, состоящей из двухвитковых спиралей больше на 18%, чем интенсивность волны для цилиндра без оболочки. Интенсивность прошедшей волны при горизонтальной поляризации на частоте 2,97 ГГц за цилиндром с искусственной оболочкой больше на 13%, чем для цилиндра без оболочки.

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

.3 Методика изготовления экспериментального образца на основе двухвитковых спиралей. Проведение эксперимента

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

Для получения более точного соответствия спиралей рассчитанным параметрам была разработана технология изготовления спиральных элементов: медная проволока наматывается на заранее изготовленный шаблон (рисунок 28), это обеспечивает взаимную идентичность двухвитковых спиралей и точное соблюдение заданных размеров. Для снятия проволоки с шаблона с сохранением параметров, вдоль продольной бороздки, нанесённой на шаблоне, на проволоке ставятся метки, соответствующие полному витку спирали. Затем проволока скручивается с шаблона, и по меткам нарезаются спиральные элементы, полностью соответствующие расчетным параметрам. На рисунке 29 показана фотография изготовленных правой и левой двухвитковых спиралей.

Рисунок 28 - Фотография шаблонов для изготовления правосторонних и левосторонних двухвитковых спиралей

Рисунок 29 - Фотография изготовленных правосторонних и левосторонних двухвитковых спиралей с углом подъема 6,35º

На основании схемы (рисунок 24) был изготовлен искусственный образец, состоящий из гладких правосторонних и левосторонних двухвитковых спиралей для реализации возможности огибания цилиндрических объектов СВЧ волнами. В качестве подложки выбран пенопласт, так как этот материал является радиопрозрачным в исследуемом диапазоне частот и не влияет на результаты измерений. На рисунке 30 показаны фотографии одного слоя маскирующего покрытия.

Рисунок 30 − Фотографии искусственной оболочки на основе двухвитковых спиралей

Следовательно, была изготовлена искусственная структура, состоящая из 12 слоев, в совокупности представляющую собой многослойную объемную структуру высотой 16 см (рисунок 31), а также металлический цилиндр радиусом 7,75 см и высотой 17 см. На каждой пенопластовой пластине были закреплены по два слоя двухвитковых спиралей.

Рисунок 31 − Фотографии объемной многослойной искусственной структуры

Для проведения эксперимента используется следующая измерительная аппаратура:

) Генератор сигналов высокочастотный Г4-80 (рисунок 33), рассчитан на диапазон частот 2,56 - 4 Ггц. Предназначен для регулировки, проверки и исследования радиоэлектронной аппаратуры в отапливаемых помещениях в пределах рабочих температур от +5°С до +40°С и относительной влажности до 95% при температуре воздуха +30°С.

Рисунок 33 − Генератор сигналов высокочастотный Г4-80

Генератор сигналов высокочастотный выполняет функции источника СВЧ сигналов, калиброванных по частоте, уровню выходной мощности и по параметрам импульсной модуляции [33].

2) Приемник измерительный П5-5Б (рисунок 34), предназначен для абсолютных и относительных измерений мощности синусоидального сигнала, измерения плотности потока мощности и напряженности электромагнитного поля при наличии измерительных антенн с известными параметрами, а также для индикации пикового значения напряжения сигналов непрерывной генерации, амплитудно-модулированных и импульсно-модулированных сигналов. Приемник измерительный обеспечивает возможность измерений в диапазоне частот 2,35 - 4 ГГц.

а − блок измерительного приемника, б − блок питания

Рисунок 34 − Приемник измерительный П5-5Б

3) Для эксперимента, в качестве излучающей и приемной антенны, используется рупорная антенна П6-23А (рисунок 35) в диапазоне сверхвысоких частот.

− кабель; 2 − волновод; 3 − рупор

Рисунок 35 − Рупорная антенна П6-23А

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

Для измерения и испытания радиотехнической аппаратуры, излучающей в свободное пространство, применяются безэховые камеры. Безэховая камера представляет собой помещение, облицованное изнутри радиопоглощающим материалом (PML) с целью уменьшения отражения от стен и обеспечения в некотором объёме камеры - безэховой зоне - заданного малого уровня отражений, то есть условий, приближающихся к условиям «свободного пространства». Стенки современных безэховых камер с высокой эффективностью экранирования от магнитных и электрических полей выполняются из железосодержащих сплавов. Тонкий стальной экран позволяет гасить магнитное поле на 60 дБ.

Используя в эксперименте вышеперечисленную аппаратуру, исследования были проведены в безэховой камере в диапазоне частот от 2,85 до 3,55 ГГц. Схема проведения эксперимента показана на рисунке 32.

− приемник измерительный П5-5Б (2,35-4,0 ГГц); 2 − генератор ГЧ-80 (2,56-4,0 ГГц); 3 − стена из радиопоглощающего материала; 4,5 − антенны измерительные П6-23А; 6 - искусственная оболочка (пенопласт, с включениями из медных двухвитковых спиралей); 7 - металлический цилиндр; 8 − безэховая камера; 9 - камера, в которой фиксируются измерения

Рисунок 32- Схема проведения эксперимента (вид сверху)

Наблюдатель (экспериментатор) находиться в камере (см. рисунок 32), в которой находиться измерительная аппаратура, и фиксирует показания на приемнике. Вручную переключает на генераторе частот нужную частоту и калибрует ее на измерительном приемнике. Для каждой отдельной частоты нужно калибровать приемник для установления точного показания в пределах погрешности прибора, поэтому процесс снятия данных проходит довольно долгое время.

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

На рисунке 33 показаны экспериментальные кривые частотной зависимости интенсивности прошедшей волны относительно падающей при различной поляризации падающей волны. Для сравнения приведены кривые из рисунка 27, которые были рассчитанны в моделировании.

Рисунок 33 − Частотная зависимость интенсивности прошедшей волны относительно падающей

Из рисунка 33 следует, что интенсивность прошедшей волны при вертикальной поляризации для экспериментальной кривой на частоте 3 ГГц за цилиндром с искусственной оболочкой, состоящей из двухвитковых спиралей больше на 12%, чем интенсивность волны за цилиндром без оболочки. Интенсивность прошедшей волны при горизонтальной поляризации на частоте 3 ГГц за цилиндром с искусственной оболочкой, больше на 15%, чем для цилиндра без оболочки. Так же наблюдается пик на частоте 3,09 ГГц, при этом, интенсивность прошедшей волны за искусственной структурой превышает интенсивность за цилиндром на 17%. Экспериментальные кривые приблизительно сходяться с кривыми, полученными в моделировании. Погрешность экспериментальных данных составляет не более 13%.

Различия экспериментальных кривых и кривых, полученных в моделировании, объясняются наличием краевых эффектов на искуственном образце, так как структуры конечной, в то время как в моделировании, физическая модель является идеализированной и бесконечной.

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

цилиндр электромагнитный волна двухвитковый

Глава 4. Оптимизация исследуемой искусственной структуры на основе двухвитковых спиралей

.1 Оптимизация искусственной структуры, состоящей из одного слоя двухвитковых спиралей

Для улучшения результатов, полученных в главе 3, проведена оптимизация исследуемой структуры. В процессе исследования было установлено, что достаточно одного ряда двухвитковых спиралей. Параметры структуры являются следующими:

a=7,75 см; b=11 см; δ=32,5 см; N=36,

где а - радиус металлического цилиндра, b - радиус искусственной оболочки, δ - расстояние от края цилиндра до спиралей, N - количество спиралей в одном слое.

Рисунок 34 − Схема расстановки спиралей в образце

На рисунке 35 показано распределение напряженности электрического поля для искусственной оболочки с одним рядом оптимизированных двухвитковых спиралей при вертикальной поляризации.

Рисунок 35 − Распределение напряженности электрического поля на частоте 2,93 ГГц при вертикальной поляризации: а) металлический цилиндр без искусственной оболочки; б) металлический цилиндр, окруженный искусственной оболочкой

На рисунке 36 показано распределения напряженности магнитного поля для горизонтальной поляризации падающей волны.

Рисунок 36 − Распределение напряженности магнитного поля на частоте 2,85 ГГц при горизонтальной поляризации: а) металлический цилиндр без искусственной оболочки; б) металлический цилиндр, окруженный искусственной оболочкой

Из рисунков 35 и 36 видно, что фронт волны за искусственной структурой частично восстанавливается для вертикальной поляризации на 2,93 ГГц и для горизонтальной поляризации на 2,85 ГГц.

На рисунке 37 показана частотная зависимость интенсивности прошедшей волны относительно падающей. Интенсивность прошедшей волны измеряется на расстоянии равным 30 см от края металлического цилиндра.

Рисунок 37 − Частотная зависимость интенсивности прошедшей волны относительно падающей

Из рисунка 37 видно, что интенсивность волны за искусственной структурой, в диапазоне частот 2,65 - 2,97 ГГц, при двух поляризациях падающей волны больше, чем интенсивность волны за металлическим цилиндром.

.2 Разработка и изготовление искусственного образца

На основании схемы (рисунок 34) был изготовлен искусственный образец, состоящий из правосторонних и левосторонних двухвитковых спиралей для реализации возможности огибания цилиндрических объектов СВЧ волнами. На рисунке 38 показана фотография искусственной оболочки, состоящая из одного ряда оптимизированных двухвитковых спиралей.

Рисунок 38 - Фотография искусственной оболочки на основе двухвитковых спиралей

Искусственная структура состоит из 14 слоев, в совокупности представляющую собой многослойную объемную структуру высотой 15 см (рисунок 39), а также металлический цилиндр радиусом 7,75 см и высотой 17 см. На каждой пенопластовой пластине были закреплены по два слоя спиралей.

Рисунок 39 − Фотографии многослойной искусственной структуры

На рисунке 40 показаны экспериментальные кривые частотной зависимости интенсивности прошедшей волны относительно падающей при различной поляризации падающей волны. Для сравнения приведены кривые из рисунка 37, которые были рассчитанны в моделировании.

Рисунок 40 − Частотная зависимость интенсивности прошедшей волны относительно падающей

Из рисунка 40 видно, что интенсивность прошедшей волны, при вертикальной поляризации, для экспериментальной кривой на частоте 2,84 ГГц за цилиндром с искусственной оболочкой, состоящей из двухвитковых спиралей больше на 31%, чем интенсивность волны за цилиндром без оболочки. Интенсивность прошедшей волны, при горизонтальной поляризации, на частоте 2,95 ГГц за цилиндром с искусственной оболочкой, больше на 30%, чем интенсивность волны за цилиндром без оболочки. В моделировании, интенсивность прошедшей волны, при вертикальной поляризации, на частоте 2,94 ГГц за искусственной структурой больше на 35%, чем за металлическим цилиндром. Интенсивность прошедшей волны, при горизонтальной поляризации, на частоте 2,85 ГГц за цилиндром с искусственной оболочкой, больше на 33%, чем интенсивность волны за металлическим цилиндром.

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

.3 Экспериментальные исследования искусственной структуры на основе прямолинейных проводников

В данном параграфе рассматривается рассеяние плоской электромагнитной волны только с вертикальной поляризацией, на искусственной структуре, состоящей из металлического цилиндра, окруженном медными прямолинейными проводниками. Прямолинейные проводники в исследуемом диапазоне частот являются не резонансными. В процессе исследования были найдены оптимальные параметры структуры в исследуемом диапазоне частот. Параметры искусственной структуры являются следующими:=77,5 мм; R=91 мм; δ=13,5 мм; h=36 мм; l=30 мм; r0=0,5 мм; N= 30,

где r - радиус цилиндра, R - радиус структуры, δ - расстояние от края цилиндра до проводников; h - толщина одного слоя структуры, l - длина проводника, r0 - радиус проводника, N - количество проводников в одном слое структуры. На рисунке 41 приведена схема расстановки прямолинейных проводников относительно скрываемого цилиндра.

Рисунок 41 − Схема расстановки прямолинейных проводников в образце

На рисунке 42(а) показана геометрия структуры, состоящей из 4 слоев прямолинейных проводников, расположенных симметрично относительно скрываемого объекта.

Рисунок 42 - Геометрия структуры: а) в моделировании; б) в эксперименте

Основываясь на рассчитанных параметрах искусственной структуры, было проведено моделирование бесконечной многослойной искусственной структуры на основе прямолинейных проводников.

Результаты моделирования представлены на рисунке 25, в виде распределения напряженности электрического поля при вертикальной поляризации падающей электромагнитной волны.

Рисунок 43 - Распределение напряженности электрического поля на частоте 3,1 ГГц: а) металлический цилиндр без искусственной оболочки; б) металлический цилиндр, окруженный искусственной оболочкой на основе прямолинейных не резонансных проводников

Из рисунка 43(б) следует, что плоский фронт волны за цилиндром с искусственной структурой на частоте 3,1 ГГц восстанавливается, при этом рассеяние волны уменьшается.

На основании полученных результатов был изготовлен искусственный образец, состоящий из 5 слоев прямолинейных не резонансных проводников (рисунок 42(б)) и металлического цилиндра. В качестве крепления проводников в структуре используется пенопласт.

На рисунке 33 показаны экспериментальные кривые частотной зависимости интенсивности прошедшей волны относительно падающей при вертикальной поляризации падающей волны и приведены кривые, рассчитанные в моделировании.

Рисунок 44 - График частотной зависимости интенсивности прошедшей волны относительно падающей

Из рисунка 44 следует, что интенсивность прошедшей волны относительно падающей, при вертикальной поляризации, для экспериментальной кривой на частоте 3,05 ГГц за цилиндром с искусственной оболочкой, состоящей из прямолинейных проводников, больше на 44%, чем интенсивность волны за цилиндром без оболочки. В моделировании, интенсивность прошедшей волны, при вертикальной поляризации, на частоте 3,12 ГГц за искусственной структурой, больше на 50%, чем интенсивность волны за металлическим цилиндром.

Так как прямолинейные проводники являются не резонансными, то восстановление волнового фронта наблюдается во всем исследуемом диапазоне частот 2,85 - 3,55 ГГц. Из экспериментальной кривой следует, что максимум интенсивность волны за цилиндром с искусственной оболочкой наблюдается в диапазоне частот от 3,03 до 3,15 ГГц.

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

Заключение

Технологии обнаружения в последние десятилетия не отставали от технологий маскировки: для каждого нового способа маскировки очень быстро находился свой способ обнаружения объекта по косвенным признакам (тени, рассеянному излучению или по путному следу). В связи с этим имеется острая потребность в такой маскировке, которая минимизировала бы все эти косвенные признаки. В последние годы было предложено несколько новых типов маскировки. Наибольшее развитие благодаря своей универсальности получила идея маскировки по методу волнового обтекания.

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

Теоретически рассмотрена задача рассеяния плоской электромагнитной волны на идеально проводящем цилиндре. Построены диаграммы рассеяния при различном параметре ka. Проверено решение этой задачи с помощью программ ANSYS HFSS и Comsol Multiphysics. Рассчитанные теоретически диаграммы рассеяния сходятся с высокой точностью с диаграммами рассеяния, построенными с помощью компьютерных программ, даже при использовании различных методов расчета.

Рассмотрены свойства дискретной оболочки, позволяющие перейти к реализации данной оболочки на практике.

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

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

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

Так же рассмотрена искусственная структура на основе не резонансных прямолинейных проводников. Данная структура смоделирована с помощью программы ANSYS HFSS и на основании полученных данных проведен эксперимент. Из полученных данных в моделировании следует, что восстановление волнового фронта наблюдается во всем исследуемом диапазоне частот 2,85 - 3,55 ГГц. Из экспериментальной кривой следует, что максимум интенсивность волны за цилиндром с искусственной оболочкой наблюдается в диапазоне частот от 3,03 до 3,15 ГГц. Поэтому, рассмотренная искусственная структура, на основе прямолинейных не резонансных проводников, при вертикальной поляризации падающей волны, дает лучший эффект огибания цилиндрических объектов СВЧ волной по сравнению с искусственными образцами, рассмотренными в главе 3 и параграфах 4.1, 4.2.

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

Список использованных источников

1.     Дубинов, А. Е. Маскировка материальных тел методом волнового обтекания// А. Е. Дубинов, Л. А. Мытарева, УФН. - 2010. - том 180, №5. - С. 475-501.

2.     Розанов Н. Н. / Н. Н. Розанов // Природа №6, 3, 2008 г.

3.      Pendry, J. B. Controlling electromagnetic fields / J. B. Pendry, D. Schurig, D. R. Smith // Science. - 2006. - Vol. 312. − P. 1780−1782.

4.     Schurig, D. Metamaterial electromagnetic cloak at microwave frequencies / D. Schurig [et al.] // Science. - 2006. - Vol. 314. - P. 977−980.

5.     Leonhardt, U. Optical conformal mapping Science. - 2006. - Vol. 312. - pp. 1777-1780.

6.      Leonhardt U, Philbin T. G. General relativity in electrical engineering // Leonhardt U, Philbin T. G. // 2006 New J. Phys. 8 247

.        Schurig, D. Calculation of material properties and ray tracing in transformation media / Schurig D., Pendry J. B., Smith D. R // Optics Express, Vol. 14, Issue 21, pp. 9794-9804 (2006).

8.      Розанов, Н.Н. Достижима ли идеальная невидимость даже при монохроматическом излучении? / Н.Н. Розанов // УФН том 181 №7 2011 г. с. 787-788.

.        Борн, М. В. Основы оптики / М. Вольф Борн // М.: Наука, 1970 г.

10.  Дубинов, А.Е. Возможна ли маскировка без сингулярностей? Новые идеи маскировки / А.Е Дубинов, Л.А. Мытарева // УФН том 182 №3 2012 г. с. 337-341.

11.   Вендик, И.Б. Метаматериалы и их применение в технике сверхвысоких частот / И.Б. Вендик, О.Г. Вендик Журнал технической физики, 2013, том 83, вып. 1, с. 3-28.

12.   Веселаго, В. Г. Электродинамика веществ с одновременно отрицательными значениями e и m / В. Г. Веселаго // УФН. − 1967. − Т. 92, вып. 3. − С. 517−526.

13.   Smith D.R., Padilla W.J., Vier D.C., Nemat-Nasser S.C., Schultz S. // Phys. Rev. Lett. 2000. Vol. 84. N 18. P. 4184-4187.

14.    Pendry J.B., Holden A.J., Stewart W.J., Youngs I. // Phys. Rev. Lett. 1996. Vol. 76. N 25. P. 4773-4776.

.        Pendry J.B., Holden A.J., Robbins D.J., Stewart W.J. // IEEE Trans. Microw. Theory. Tech. 1999. Vol. 47. N 11. P. 1075-1084.

.        Shelby R.A., Smith D.R., Schultz S. // Science. 2001. Vol. 292, P. 77-79

.        Cummer, S. A. Full-wave simulations of electromagnetic cloaking structures / Steven A. Cummer*, Bogdan-Ioan Popa, David Schurig, and David R. Smith // Phys. Rev. E 74, 036621 (2006) [5 pages].

.        Cai W., Chettiar U.K., Kildishev A.V., Shalaev V.M. // Nature Photonics. 2007. V.1. №4. P.224-227.

.        Jiang, W. X. Arbitrarily elliptical-cylindrical invisible cloaking / Wei Xiang Jiang, Tie Jun Cui, Guan Xia Yu, Xian Qi Lin, Qiang Cheng and Jessie Yao Chin 2008 J. Phys. D: Appl. Phys. 41 085504.

.        Jiang, W. X. Analytical design of conformally invisible cloaks for arbitrarily shaped objects / Wei Xiang Jiang, Jessie Yao Chin, Zhuo Li, Qiang Cheng, Ruopeng Liu, and Tie Jun Cui // Phys. Rev. E 77, 066607 (2008) [6 pages]

.        Yao, K. Two-dimensional dissimilar electromagnetic cloak for irregular regions / K. Yao, C. Li, F. Li // Applied Physics B, 2009, V. 96, pp 355-361.

23.  Гольштейн, Л. Д. Электромагнитные поля и волны / Л. Д. Гольштейн, Н. В. Зернов // Изд. 2-ое перераб. и дополненное, М. Изд-во “Советское радио”, 1971 г., с. 359-365.

24.   Банков, С.Е. Проектирование СВЧ устройств и антенн с Ansoft HFSS / С.Е. Банков, А.А. Курушин // М.: «Солон».- 2009.- C. 194-197.

25.    Asghar, M. Electromagnetic cloaking with a mixture of spiral inclusions// M. Asghar et al. - Proc. Metamaterials’ 2007, Rome, Italy. - P. 957.

26.    Фаняев И.А., Безотражательное огибание цилиндрических объектов электромагнитными волнами, Творчество молодых ¢2012 Сборник научных работ студентов и аспирантов, Гомель 2012, с. 229-232.

.        Фаняев И.А., Огибание цилиндрической поверхности методом волнового обтекания, Дни студенческой науки, Материалы XLI студенческой научно-практической конференции, Гомель, 2-3 мая 2012 года, с. 58.

28.    Guven, K. Electromagnetic cloaking with canonical spiral inclusions// K. Guven, E. Saenz, R. Gonzalo, E. Ozbay and S. Tretyakov. - New Journal of Physics. - 2008, 115037. - no. 10.

.        I.A. Faniayeu, V.S. Asadchy, The calculation of all tensors components of the dielectric, magnetic and chiral polarizabilities for helices with the optimal parameters, 11th International Conference on Global Research and Education in Engineers for Better Life (inter-Academia), 27-30 August, 2012, Budapest, Hungary, pp. 497-504.

.        I.A. Faniayeu, V.S. Asadchy, I.V. Semchenko, S.A. Khakhomov, Calculation and analysis of the tensors of electric, magnetic and chiral polarizabilities of the helices with optimal shape, Metamaterials '2012: The Sixth International Congress on Advanced Electromagnetic Materials in Microwaves and Optics, St. Petersburg, September 17 - 22, 2012.

.        Semchenko, I.V. Radiation of Circularly Polarized Electromagnetic Waves by the Artificial Flat Lattice with Two-Turns Helical Elements// I.V. Semchenko, A.L. Samofalov, S.A. Khakhomov. - Bianasotropics' 2004. - 10th International Conference on Complex Media and Metamaterials, Het Pand, Chent, Belgium. - P. 236-239.

32.   Tretyakov, S. A. Analytical Modeling in Applied Electromagnetics / S. A. Tretyakov. − Norwood: Artech House, 2003. − 260 p.

33.    Генератор сигналов высокочастотный Г4−80. Описание, инструкция по эксплуатации. − Горький, 1977. − 102 с.