Московский энергетический институт (Технический университет), кафедра Антенных устройств и распространения радиоволн
СРАВНЕНИЕ ВАРИАНТОВ ИМПЕДАНСНЫХ ГРАНИЧНЫХ УСЛОВИЙ ПРИ ПАДЕНИИ ПЛОСКОЙ ЭЛЕКТРОМАГНИТНОЙ ВОЛНЫ НА ДИЭЛЕКТРИЧЕСКОЕ ПОЛУПРОСТРАНСТВО
И.В. Белькович., М.А. Жексенов, Д.А. Козлов
Аннотация
электромагнитный волна импедансный диэлектрический
Рассмотрены три варианта импедансных граничных условий (ИГУ) на примере падения плоских электромагнитных волн перпендикулярной и параллельной поляризаций на границу раздела воздуха и диэлектрического полупространства. Поверхностный импеданс второй среды в ИГУ берется равным импедансу при падении волны: по нормали к поверхности (ИГУ Леонтовича-Щукина), под углом Брюстера, по касательной к поверхности. Изучены зависимости относительной погрешности модуля и максимальной абсолютной ошибки фазы коэффициента отражения от мнимой части диэлектрической проницаемости.
Ключевые слова: плоская электромагнитная волна, параллельная и перпендикулярная поляризации, диэлектрическое полупространство, отражение, импедансные граничные условия.
Abstract
Three variants of impedance boundary conditions in case of plane wave incidence on interface between air and dielectric are examined. Parallel and perpendicular polarizations are considered. In given boundary conditions the surface impedance of the second medium is taken equal to impedance in case of: normal incidence (Leontovich-Schukin IBC), incidence at an Brewster angle, grazing incidence. Angular dependences of reflection coefficient magnitude error of all IBCs and dependecies of maximum absolute error of reflection coefficient phazor from imaginary part of complex permittivity are studied.
Keywords: plane electromagnetic wave, parallel and perpendicular polarizations, reflection, dielectric halfspace, impedance boundary conditions.
Введение
Приближенные импедансные граничные условия (ИГУ), предложенные в работах А. Н. Щукина [1] и М. А. Леонтовича [2], стали широко применяться в теоретических исследованиях с появлением сборника [3], в котором была опубликована статья [2] и ряд других работ, использовавших ИГУ. М. А. Леонтовичем [2] были в общем виде сформулированы ИГУ на поверхности тел с комплексными диэлектрической е и магнитной м проницаемостями
, (1)
где - внешняя нормаль к поверхности тела, и - компоненты электрического и магнитного полей во внешней области, - импеданс среды. В работе [2], где принята гауссова система единиц, Z0=1. При использовании международной системы единиц Z0=120·р Ом, а е и м - относительные проницаемости среды. Представление (1) основано на допущении, что волна распространяется в среде с потерями по законам геометрической оптики.
Далее в [2] были поставлены два вопроса: «во-первых, дать оценку ошибок, связанных с изображение поля в виде волны, распространяющейся внутрь тела по законам геометрической оптики, и, во-вторых, выяснить, при каких условиях это поле может быть представлено в виде волны такого вида».
Ответ на первый вопрос М. А. Леонтович извлек из работы С. М. Рытова [4]. В [4] с помощью метода малых возмущений было построено разложение электромагнитного поля по малому параметру, в качестве которого была выбрана величина, обратная модулю комплексного показателя преломления среды. Ответ на второй вопрос был сформулирован следующим образом: граничное условие (1) применимо при выполнении следующих требований: «глубина проникновения в тело и длина волны в нем малы по сравнению с длиной волны в окружающем пространстве, по сравнению с расстояниями от источников поля и по сравнению с радиусами кривизны поверхности тела. Изменения е и м тела на длине волны (или на длине, равной глубине проникновения), малы».
Отметим, что из работы С. М. Рытова [4] следуют не только ИГУ (1), но и приближенные граничные условия (ПГУ) с учетом вклада следующего порядка малости.
Дальнейшее развитие ПГУ было связано с ослаблением требований, наложенных в [2], и определением ПГУ из решений более общих модельных задач. В ряде работ (см. [5] и цитированные там работы) было предложено использовать в качестве ИГУ поверхностный импеданс для скользящего падения плоской волны на диэлектрическое полупространство. Поверхностный импеданс, найденный для случая падения плоских волн на кусочно - однородную плоскослоистую среду [6], может рассматриваться как ПГУ для искривленных поверхностей. ИГУ были предложены для сильно неоднородной по нормали к поверхности плазмы в [7,8]. Результаты анализа различных вариантов ПГУ систематизированы в монографиях [9-11], см. также недавние работы [12-14].
Вместе с тем ряд вопросов, связанных с ИГУ, остается до сих пор не изученным. Так, с физической точки зрения часто представляет интерес не поверхностный импеданс, а коэффициент отражения (КО) волны, и оценку применимости ИГУ разумно определять по погрешности расчета КО. В [9] оценена точность ИГУ Леонтовича-Щукина в задаче падения плоской волны на однородное диэлектрическое полупространство Отмечено, что в случае нормального падения формулы для КО, использующие ИГУ, являются точными для любого значения показателя преломления. Утверждается, что при фиксированном показателе преломления их точность уменьшается и является наихудшей при скользящем падении. В [9] приведен также пример использования граничного условия, в котором импеданс второй среды берется при угле Брюстера.
Цель данной работы - более детально изучить влияние выбора различных ИГУ на величину КО на примере простейшей модельной задачи о падении плоских электромагнитных волн на диэлектрическое полупространство с потерями при не слишком больших значениях модуля диэлектрической проницаемости. В прикладном плане такая задача возникает при анализе распространения и дифракции электромагнитных волн в естественных условиях (рассеяние волн на земной и морской поверхности, льдах и т.д.) и в искусственных диэлектрических структурах.
Исходные формулы
Выражения для коэффициентов отражения плоских волн линейной поляризации от диэлектрического полупространства с потерями хорошо известны [6]. Пусть - отсчитываемый от нормали угол падения плоской волны из среды 1 (свободное пространство) в среду 2 (полупроводящая среда). В случае перпендикулярной поляризации КО по электрическому полю равен
, (2)
В формуле (2) - величины поверхностного импеданса сред 1 и 2, определенные через отношения тангенциальных составляющих электрического и магнитного полей, падающих в среде 1 и преломленных в среде 2 волн.
Аналогично, в случае параллельной поляризации КО по магнитному полю имеет вид
, где , . (3)
Результаты расчетов
В этой части проведено исследование влияния величин действительной и мнимой части диэлектрической проницаемости второй среды на модуль и фазу КО. Наряду с точным выражением для КО рассмотрены три приближённых выражения, отличающиеся углом , при котором фиксируется поверхностный импеданс второй среды :
1. нормальное падение волны (), при этом тангенциальные компоненты полного поля во внешней области связаны с помощью граничного условия Леонтовича -Щукина (1);
2. скользящее падение волны на полупроводящее полупространство ();
3. падение волны под углом полного преломления (углом Брюстера ).
Коэффициенты отражения, соответствующие каждому из указанных случаев, будем обозначать , , i=1, 2, 3.
При использовании ИГУ Леонтовича-Щукина приближенное значение импеданса совпадает для обеих поляризаций падающего поля, т.е.
. В двух других случаях приближенные значения импедансов различаются для каждой поляризации, поэтому вторая среда характеризуется диагональным тензором импеданса, и, сответственно, в формуле (1) скалярная функция импеданса должна быть заменена тензором.
Была проведена серия расчетов модулей и фаз КО по строгим и приближенным формулам для следующих значений диэлектрической проницаемости . Такой выбор обусловлен тем, что они включают типичные для сухого грунта и ряда строительных материалов значения проницаемости в диапазоне метровых - сантиметровых волн. Заметим, что ИГУ для полупространства принято вводить только при наличии потерь в среде, чтобы исключить отражения от объектов, которые могут находиться внутри среды. Тем не менее мы рассмотрим случай среды без потерь, как предельный.
Результаты расчета модулей КО по строгим формулам для двух поляризаций приведены на рис. 1 а) и рис. 2 а) для диэлектрической проницаемости второй среды, равной . Угол Брюстера для данной среды равен . На рис. 1 б) и рис. 2 б) приведены зависимости относительной погрешности , КО для трех вариантов ПГУ. Погрешность определяется следующим образом
, (4)
где - точное выражение для КО при определенной поляризации,
- приближенное выражение для КО той же поляризации, i=1, 2, 3.
На рис. 1в), 2в) для сравнения приведены относительные погрешности поверхностных импедансов.
Сплошные кривые на рис. 1б), 1в) и 2б), 2в) соответствуют заданию ИГУ при , пунктирные - при , штриховые - при .
Из зависимостей рис. 1б), 2б) видно, что относительная погрешность всех трех выражений для приближенных КО обеих поляризаций обращается в нуль при угле падения, стремящемся к . В этом пределе проявляется приципиальное отличие от поведения погрешности поверхностного импеданса в вариантах 1 и 3 (рис.1в), 2в)), которая является достаточно большой. Заметим в этой связи, что Леонтовичем было проведено сравнение приближенного коэффициента отражения Френеля со строгим для параллельной поляризации и сделан вывод, что ошибка приближенного коэффициента отражения, имеет порядок [2]. Однако эта оценка несправедлива при скользящих углах падения, когда погрешность приближенного коэффициента отражения стремится к нулю, как .
Обратимся к сравнительному анализу погрешностей приближенных КО.
При использовании ИГУ Леонтовича-Щукина (вариант 1) КО имеют нулевую погрешность при нормальном падении плоской волны обеих поляризаций. В случае перпендикулярной поляризации относительная погрешность модуля КО в первом варианте ИГУ максимальна при промежуточных углах (в данном случае при угле ).
а)
б)
в)
Рис. 1 Зависимости модуля строгогоКО а), относительных погрешностей приближенных выражений для КО б), относительных погрешностей поверхностных импедансов в) от угла падения в случае перпендикулярной поляризации
а)
б)
в)
Рис. 2 Зависимости модуля строгогоКО а), относительных погрешностей его приближенных выражений б), относительных погрешностей поверхностных импедансов в) от угла падения в случае параллельной поляризации
Второй вариант ИГУ дает наибольшую относительную погрешность вычисления КО по сравнению с другими вариантами при нормальном падении волны. При скользящем падении, естественно, этот вариант ИГУ дает наилучшие результаты.
При угле Брюстера относительная погрешность КО параллельной поляризации для первого и второго вариантов ИГУ , т.к. строгий КО равен нулю. Отметим, что при значениях КО, близких к нулю, информативнее становится абсолютная погрешность .
Использование вариантов 1 и 2 ИГУ приводит к тому, что положение угла Брюстера смещается для этих приближений по сравнению с истинным значением угла Брюстера. Соответственно, для углов, в которых приближенные выражения КО равны нулю, относительная погрешность , что непосредственно вытекает из выражения (4). В тех точках, где точное значение КО , погрешность .
Третий вариант приближенного расчета КО имеет нулевую погрешность при падении волны перпендикулярной поляризации под углом Брюстера. Относительная погрешность КО в этом варианте максимальна при и равна .
Сравним результаты расчета приближенных КО между собой. Для среды без потерь в случае перпендикулярной поляризации метод вычисления модуля КО с использованием ИГУ Леонтовича - Щукина (вариант 1) в диапазонах углов падения волны 0°ч48° и 0°ч39° имеет меньшую погрешность, чем второй и третий варианты ПГУ соответственно. Для обеих поляризаций погрешность КО при использовании первого варианта ИГУ меньше погрешностей вариантов 2 и 3 в широком секторе углов, прилегающем к нормальному падению (0о-40о), а погрешность КО при использовании третьего варианта ИГУ меньше погрешностей остальных вариантов в окрестности угла Брюстера (40о -70о ). Если учесть, что при скользящих углах падения малую погрешность КО дают второй и третий варианты задания ИГУ, то можно рекомендовать последний для использования в широком диапазоне углов падения, включающем скользящие углы.
Что касается фаз приближенных выражений для КО, то в случае среды без потерь при перпендикулярной поляризации они полностью совпадают с фазой точного выражения во всем секторе углов падения . В случае параллельной поляризации только фаза третьего варианта КО совпадает с , фазы остальных приближенных выражений для КО сдвинуты относительно точного результата или, что тоже самое, сдвинуты величины углов Брюстера . Для первого и второго вариантов ИГУ они равны и соответствено.
Рассмотрим случай среды с потерями, приняв диэлектрическую проницаемость второй среды равной . Обозначения кривых КО на рис. 3 и 4 аналогичны рис. 1 и 2. На рис. 3 а) и 4 а) изображены угловые зависимости КО при наличии потерь во второй среде. Так как е комплексная величина, угол Брюстера уже не является действительным