Электронный журнал «Труды МАИ». Выпуск № 41
www.mai.ru/science/trudv/
УДК 538.574.6
Распространение электромагнитной волны в плоскослоистом диэлектрике вблизи нуля диэлектрической проницаемости
И. С. Васильев, И. П. Козлов
Аннотация.
Рассматривается строгое решение задачи нормального падения плоской электромагнитной волны на плоскослоистый диэлектрик конечной толщины вблизи нуля
диэлектрической проницаемости . Приводятся результаты расчетов.
Ключевые слова: строгое решение; нормальное падение; плоская электромагнитная волна; плоскослоистый диэлектрик; ноль диэлектрической проницаемости.
Введение
Исследования задачи нормального падения электромагнитной волны на плоскослоистый диэлектрик представляют интерес из-за выявленной в нуле в в случае среды без поглощения критической точки, вблизи которой решение качественно зависит от малых
изменений параметров физической задачи [1]. При линейной зависимости в среде без
поглощения средний (по времени) поток энергии ~ чб7Г' " ' ^ при £ где
и напряженности электрического и магнитного полей двух ортогональных поляризаций. Вблизи нуля в: имеет место аномалия коэффициента отражения, который почти мгновенно
£ = , £ —'■> Оч
возрастает до единицы; поверхность (при ) сворачивается в точку; решение
в р. неустойчиво по направлению распространения волны , так возмущение ™ при любом
£ приводит К 15 ~~' " ~ при ^ Н \ Е° £п51П (Д#„) СОПЗI _ инвариант на сдое.
Введение поглощения делает математическую задачу устойчивой. Этот случай как раз и исследуется в настоящей статье, задача решается с использованием функций Эйри. В связи
с выявленной критической точкой в нуле £ задача представляет как теоретический, так и практический интерес. Полученные результаты расчетов находятся в соответствии с
классическими представлениями о поведении полей вблизи нуля диэлектрической
проницаемости [2]. Рассматривается длинноволновое приближение, которое может широко применяться на практике.
Решение задачи
Рассматривается нормальное падение плоской электромагнитной волны на изотропный неоднородный плоский слой диэлектрика с поглощением произвольной
толщины от 11 до ■ при соответствующих значениях диэлектрической
проницаемости с и
(рис. 1).
До слоя """ и за слоем среда однородная, причем за слоем имеется только проходящая волна. Слой состоит из двух частей. На первом
- -
участке от с до
1 (Ч
)
Рис. 1. Линейный слой сопряженный в окрестности нуля е с нелинейным. Кривые 1-2 соответствуют различным аппроксимациям линейного слоя.
диэлектрическая проницаемость зависит от координаты линейно
¿1
а =
¿ОЫ
= < О
а на втором
гч г.,- £.:Ч ~ £г(кг) = В2/[4ГЫ- кх)2~\ ~ п
участке от 1 до ( , ) нелинейно 1 , слои гэлея [3], где
дг(ё') = —— const
. На слое функция
непрерывна вместе с
к = 2тт /Я Л г т-г
производной, где , - длина волны в свободном пространстве. Полученное
решение позволяет при £° 1 производить расчеты для малых а (при 1 расчет
функций Эйри ограничен по параметру
Будем искать решение стационарного волнового уравнения
а* я,
[¿кг1
(к¿)ЕХ = О,
(1)
на неоднородном слое от а до ® (рис. 1). Поле на однородном полубесконечном слое при
£ = £„
представим суперпозицией
Аа
падающей (с амплитудой " ) и отраженной (с амплитудой ) волн
а поле справа от точки
г = г.
■ проходящей волной
Пусть
- = -,-
определяется пересечением продолжения линеиного слоя с осью
т Ь' = t- ¿¿Г (кг.")* ,г = (кгГ.>г(кг)) кг.-, = 1 /а
Тогда, производя замену переменных: к , ( ^ ^ ^ и ,
г<га = (Ь^ - ЮДк^3
принимает вид
(3)
(4)
волновое уравнение (1) при
йг:1 л
Его решение будет
«СО г'С0
где 4 функции Эйри в обозначениях работы [4, доп. 2], а Б и й постоянные
коэффициенты, причем
Удовлетворение граничным условиям непрерывности тангенциальных составляющих
полей (2), (4) и их производных справа (знак +) и слева (знак -) от точки
Е\, _ = Е
йЕ _ ¿Е
хаГ ¿Е
приводит к неполной системе уравнений относительно
Д
сС-
С Г
коэффициента отражения и постоянных коэффициентов и
(5)
(6)
Д
а С- — пс _ 2 ¿V, г'пкг,
с С-
где принято ~ , а тогда
Д
Исключая из уравнений (5), (6), получим уравнение относительно Б и О
где
Для сшивания линейной и нелинейной частей слоя поле около точки
- = -1
при
£ = £-,
представим в виде
Е = .4 (1 — й
ю
справа
где ' ^ определяют поле в локальной системе координат с центром 01 , при
г = О А1 = А™е
и
(в системе координат с центром 0, где
Аналогичным образом как в точке ~ ~а, удовлетворение граничным условиям в
точке ~ после исключения из полученной системы уравнений определяется
ниже) и ряда преобразований приводит к следующим выражениям для F и G
С - -Ву1 Г = Ви1
(8)
Б =
где
— ч - ц ^ ^ ^ _ ь ! ■ 11 .
На этом заканчивается математическое решение задачи. Поле ^ на участке от 3 до
г ДсС'
1 определяется подстановкой (8) в (4), а коэффициент отражения находится из (5)
В свою очередь поле определяется выражением (7), где
Решение для нелинейной части слоя, до слоя (при £ _ £ 1) и за слоем (при 5 ~ 5 лг)
среда однородная, в системе координат с центром 01, при -1 = 0, через элементарные
¡1С
функции для коэффициента отражения ^ и поля ^ Е 1 на слое представляется в виде:
где
Отметим, что линейный слой при сдвиге фаз волны ^ <:' ^ на слое в качестве длинноволнового приближения можно заменить нелинейным слоем.
Результаты расчетов
Результаты расчетов полей при 1 и в области в > 0 показываются на
рис.2 в зависимости от толщины линейного слоя, определяемого величиной
1^(01
а = ёв/с^кг). Первый максимум поля ■ - имеет резонансный характер. Резкое увеличение
£ = 0 л
напряженности электрического поля в точке, где , может привести к необходимости решения нелинейной задачи о распространении волн. Постепенное возрастание максимумов поля при увеличении толщины слоя объясняется свойством асимптотики поля [2]. Расчет
Е, = 10
-30
показывает, что при соответствует известной теории.
filc| = 1
(совпадение до 15 - го знака после запятой), что
п , £ = £(кЕ)
На рис. 3 дается зависимость * на нелинейном участке слоя при разных
значениях поглощения, определяемого мнимой частью
. Точка 8=0
й Е > О £ < 0 т-г
обходится как сверху при , так и снизу при . Проведено сравнение характеристик
линейного и нелинейного слоев вблизи нуля 8 при совпадении £ в начале слоя (£ а и
1£ £ £ £ *) и Яе( ) в конце ( 1 и соответственно Л). Расчеты показывают, что при
в случае линеиного и нелинейного слоев -
а 0.82- 2ремя ВЫЧИсления поля для линейного слоя на 1-2 порядка больше, чем для
нелинейного слоя. Это объясняется тем, что решение для нелинейного слоя представляется через элементарные функции. На рис. 4 представлены результаты расчетов полей при
изменении поглощения в линейном слое для
е. = -3 X 10
—з
и менее. В области
отрицательных значений 8 поле затухает по экспоненте, как это и следует из теории [4].
Рис. 2. Зависимость поля Е(-) стоячей волны, е" = 0, sa= 1, s1 = 1.2-10"5 , от толщины линейного слоя при изменении
а ее de/d(kz), кривые: 1 - а = -2. МО"5, 2 - а = -2.0-10"4, 3-а = -1.9-10"2, 4 - а = -1.П-10"1, 5 - а = -2.0-10"3.
Рис. 3. Поведение функции к(к/) при различных значениях 1т(81), кривые: 1 - 1т(Б!) = 2.4-10"5, 2 - 1т(8!) = 1.14-10"5,
3 - 1т(8!) = 0.
1® 1 .-у ..... 80} 1 4 \ > ■ ч Д v \ U __________¿q
\ / Т. /> \ / V» |
-4-2 0 1 4
t
Рис. 4. ПоляД-) при а = -2.МО"5, sa=0.12, к, = -3-10"3 (t, = -3.9) и изменении поглощения в линейном слое: 1 - б" = 0, 2 - s"= -8Н0"6, 3 - s"= -8i-10~\
Заключение
Решена задача о распространении волн в слое конечной толщины вблизи точки f
Расчетным путем показано, что для нормального падения волны на плоскослоистый
~ , ~ ds/dfkz) Е
диэлектрик с линеинои функцией 4 * при малых величина v у напряженность поля
вблизи нуля в может существенно возрастать. Задача о распространении волн в общем
случае является нелинейной. Решение этой задачи возможно методом прогонок при
использовании метода самосогласованных конечных разностей [1]. Полученное
длинноволновое приближение для линейного слоя может иметь широкое применение.
Библиографический список
1. Козлов И.П. Исследование прохождения электромагнитной волной плоского слоя диэлектрика вблизи критической точки // Письма в ЖТФ. - 2000. - Т.26, Вып. 14. - С.28-35.
2. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред // М.:Наука - 1973.
3. Стретт Дж. В. (Лорд Рэлей) Теория звука (1, 148 б) // М.:Гостехиздат - 1955.
4. В.А.Фок, Проблемы дифракции и распространения электромагнитных волн// М.: "Сов. радио". - 1970. - 560с.
Сведения об авторах.
Васильев Иван Сергеевич, аспирант Московского Государственного Университета Леса (государственного технического университета).Ул. Горького, 14Б, кв. 357, Королев, 141080; тел.: 762-17-31; e-mail: boy2k@list.ru
Козлов Игорь Петрович, д.т.н. Ул. Первомайская, 1, кв. 44, Мытищи-5, 141005, 516-5169.