УДК 53(075.8)
А. М. Макаров, Л. А. Л у н е в а, К. А. М а к а р о в
О СОПРЯЖЕНИИ ПЛОСКИХ ГАРМОНИЧЕСКИХ ВОЛН НА ПОВЕРХНОСТИ РАЗДЕЛА ДВУХ ОДНОРОДНЫХ ИЗОТРОПНЫХ СРЕД В КЛАССИЧЕСКОЙ ЭЛЕКТРОДИНАМИКЕ
Описание взаимодействия электромагнитной плоской гармонической волны с границей раздела двух материальных сред, если хотя бы одна из них является проводящей, помимо известных четырех условий сопряжения касательных составляющих векторов напряженности электрического и магнитного полей требует выполнения закона сохранения стороннего поверхностного электрического заряда, что эквивалентно условию скачка нормальных составляющих вектора электрического смещения. Проведен анализ физического содержания "комплексного" вектора электрического смещения, являющегося следствием использования понятия комплексной диэлектрической проницаемости среды.
Отражение и преломление плоской гармонической электромагнитной волны на плоской границе раздела двух однородных изотропных прозрачных сред (диэлектриков) в классической электродинамике изучено досконально. Законы отражения и преломления (в частности, закон Снеллиуса) вошли в элементарный курс физики, формулы Френеля — в университетский курс общей и прикладной физики [1-6]. В случае взаимодействия плоской гармонической электромагнитной волны с границей раздела двух непроводящих однородных и изотропных сред (диэлектрик-диэлектрик) обращаются в нуль объемные плотности сторонних электрических зарядов и объемные плотности токов проводимости, а также поверхностная плотность сторонних электрических зарядов и поверхностная плотность токов проводимости на границе раздела сред. При этом оказывается, что непрерывность нормальных к границе раздела компонент векторов магнитной индукции и электрического смещения является следствием непрерывности векторов напряженностей электрического и магнитного полей [2]. Особенности явления отражения-преломления плоской гармонической электромагнитной волны от границы раздела двух сред, если вторая среда является проводящей (т.е. поглощающей), в рамках линейной классической электродинамики обычно рассматривают с помощью введения комплексной диэлектрической проницаемости е. При этом практически общепринятым является утверждение, что формальной замены действительной диэлектрической проницаемости е на комплексную
достаточно для получения исчерпывающей информации о рассматриваемом явлении: "Проводящая среда в математическом отношении отличается от диэлектрика лишь тем, что в уравнении для проводящей среды вместо абсолютной диэлектрической проницаемости е входит комплексная диэлектрическая проницаемость" [5, с. 272]. Следствием этого является использование условий непрерывности касательных компонент векторов напряженностей электрического и магнитного полей (при нулевой поверхностной плотности токов проводимости) на границе раздела двух сред по аналогии со случаем падения электромагнитной волны на границу раздела двух непроводящих сред. Однако, в рассматриваемом случае непрерывность нормальных по отношению к границе раздела составляющих электромагнитного поля должна быть проверена.
Ниже покажем, что в случае падения электромагнитной волны на границу раздела двух однородных изотропных проводящих сред поверхностная плотность сторонних электрических зарядов отлична от нуля. При этом выясним, почему формальный переход от действительной к комплексной диэлектрической проницаемости не приводит к ошибкам при определении коэффициентов отражения и преломления
и каково физическое содержание векторного поля О = е0еЕ, нормальные компоненты которого непрерывны на границе раздела проводящих сред (монохроматическое приближение), в сравнении с вектором
О = ее0 Е.
Рассмотрим систему уравнений классической электродинамики в линейном приближении для случая однородной изотропной среды:
ШуО = р; = 0;
- дБ Гт - дО (1)
к*Е = --; гоШ = , + -; 00
О = еоеЕ; Б = ^Я; = 7 Е.
Все физические величины в системе уравнений (1) — действительные, обозначения общепринятые; 7 — электрическая проводимость; коэффициенты в материальных уравнениях среды считаются независящими от пространственных координат и времени.
Ищем решение системы уравнений (1) в форме бегущих плоских гармонических волн вида
а = а0 ехр (¿(к •г — шЬ)), (2)
где а0 — постоянная амплитуда; к — волновой вектор; г — радиус-вектор точки наблюдения; ш — круговая частота волны; Ь — время. Волновой вектор к в выражении (2) является постоянной векторной
величиной. Система уравнений (1) переходит в алгебраическую систему для комплексных амплитуд (или мгновенных физических величин):
гк • D = р; (3)
гк • В = 0; (4)
гк х E = гшВ; (5)
ik х H = j — iuD. (6)
Система уравнений (3)-(6) может быть записана отдельно как для первой, так и для второй среды, причем в первой среде имеет место и падающая, и отраженная волна. Параметры падающей волны обозначены символами без верхних индексов, параметры отраженной волны выделим единичным, а преломленной волны — двойным штрихом. Вектор единичной нормали к плоской границе раздела двух сред п направим из первой среды во вторую и отметим соответствующими нижними индексами принадлежность величин к первой или второй среде.
На плоской неподвижной границе раздела двух проводящих сред должны выполняться следующие условия сопряжения для результирующего электромагнитного поля. Условие непрерывности касательных компонент векторов напряженности электрического поля имеет вид
п х (п х Ех) = п х (п х Е2); (7)
условие непрерывности нормальных компонент векторов магнитной индукции —
п • Ёх = п • В2; (8)
условие непрерывности касательных компонент векторов напряженности магнитного поля, если на поверхности раздела сред отсутствуют поверхностные токи проводимости (последнее должно быть установлено в процессе решения задачи), может быть записано как
п х (п х Нх) = п х (п х Н2); (9)
условие для скачка нормальных компонент векторов О с учетом возможности наличия на границе раздела поверхностной плотности сторонних электрических зарядов —
п • О - п • 01 = а (10)
и условие сохранения электрического заряда на поверхности раздела —
п • 32 - п • 31 = гша. (11)
Последнее из условий записано с учетом гармонического изменения во времени поверхностной плотности сторонних электрических зарядов
и в предположении, что отсутствует их поверхностное растекание. Эта посылка будет оправдана, если будет выполнено, в частности, условие отсутствия поверхностных токов проводимости на границе раздела двух сред.
Решение в форме плоских гармонических бегущих волн обуславливает одинаковую круговую частоту всех волн и геометрические условия сопряжения
п х (п х к) = п х (п х к') = п х (п х к''). (12)
Начало декартовой системы координат {х, у, г} совместим с границей раздела, положительное направление оси г совпадает с направлением нормали к границе раздела, нижнее полупространство (г < 0) заполнено средой 1, верхнее (г > 0) — средой 2. Развернутая форма соотношений (12) имеет вид
кХ кх; ку ку; кХ К ; ку ку. (13)
Кроме того, нормальные компоненты волнового вектора падающей волны и волны, отраженной от поверхности раздела двух сред, связаны соотношением
К = —К. (14)
Соотношение (14) справедливо, поскольку модули рассматриваемых векторов одинаковы в силу распространения волн в одной и той же среде (дисперсионные уравнения тождественны друг другу), а касательные компоненты одинаковы в силу условий (13).
Условия (7) и (9) в рассматриваемой системе координат приводят к уравнениям
Ех + ЕХ = ЕХ'; Еу + ЕУ = ЕУ'; (15)
Нх + Н'х = Н'Х; Ну + Н'у = Н'У. (16)
В дальнейших рассуждениях соотношения (15) и (16) будем рассматривать в качестве основополагающих, привлекая по мере необходимости уравнения (3)-(6) для каждой из рассматриваемых волн с учетом среды их распространения.
Непрерывность нормальных составляющих векторов магнитной индукции. Из уравнения (5) следует
Б = кхЕ. (17)
ш
Нормальная к границе раздела компонента вектора магнитной индукции может быть записана в форме
Бг = кхЕу — куЕх. (18)
ш
Условие непрерывности нормальных компонент вектора магнитной индукции (8) принимает вид
кхЕу — ку Ех + к'х Е'у — к'у Е'х = к'£Е'У — куЕ'£. (19)
Легко видеть, что при выполнении условий сопряжения касательных компонент волновых векторов (13) и условия непрерывности касательных компонент вектора напряженности электрического поля (15) уравнение (19) удовлетворяется тождественно. Таким образом, уравнение (8) в рамках рассматриваемой задачи может рассматриваться как следствие условий непрерывности касательных компонент вектора напряженности электрического поля на границе раздела сред.
Скачок нормальных составляющих вектора б. Из уравнения (6) следует соотношение
б = к Х Й - 1 ^ 3 (20)
ш
Нормальная к границе раздела компонента вектора б может быть записана в форме
б = кхНУ — куНх — Ъ • (21)
ш
С учетом зависимости (21) потребуем выполнения условия (10): к" Н" — к"Н" — 3 "
кхНу к'У Нх
Ш
kxHy kyНх ijz i кхНу kyНх ijz
+ -^-^ = ^ (22)
ш ш
В полученом выражении для зависимости величины поверхностной плотности сторонних электрических зарядов а используем условия (13) сопряжения касательных компонент волновых векторов и условие (16) непрерывности касательных компонент векторов напряженности магнитного поля:
Я- 3* + 3*) = ъша. (23)
Соотношение (23) совпадает с уравнением (11) сохранения электрического заряда на границе раздела сред в предположении, что на ней отсутствуют поверхностные токи проводимости или их природа такова, что не имеет места поверхностное растекание электрического заряда. Предположение об отсутствии поверхностных токов проводимости согласуется с ранее принятым предположением о непрерывности касательных компонент векторов напряженности магнитного поля.
Таким образом, условие скачка нормальных компонент вектора б и условие сохранения электрического заряда на поверхности раздела двух материальных сред в рамках рассматриваемой проблемы не
являются независимыми друг от друга. Этого можно было ожидать: система (1) уравнений классической электродинамики неявно содержит в себе уравнение сохранения электрического заряда. Переход к интегральной форме системы (1) и использование ее для формулировки условий сопряжения на границе раздела материальных сред только выявляет органическую взаимосвязь составляющих электромагнитного поля.
Скачок нормальных компонент вектора О на границе раздела материальных сред органически связан с законом сохранения электрического заряда на этой поверхности. Этот скачок исчезает в нескольких случаях: если обе материальные среды являются непроводящими (З' = З = 3 = 0); если обращаются в нуль нормальные компоненты вектора объемной плотности электрического тока проводимости (з'г = Уг = Зг = 0), что имеет место при определенной поляризации падающей плоской электромагнитной волны.
Комплексная диэлектрическая проницаемость среды и комплексный вектор О. Рассмотрим правую часть уравнения (6). Введем комплексный вектор О соотношением
—шО + З = —1шБ. (24)
Поскольку в левую часть соотношения (24) вектор напряженности электрического поля Е входит линейно, естественно предположить, что это свойство имеет место и для правой части:
О = ££оЕ, (25)
где е — комплексная диэлектрическая проницаемость материальной среды, которая выражается формулой
е = е + г—. (26)
еош
Удобство введения комплексного вектора О определяется не только упрощением уравнения (6):
гк х Н = —гшО, (27)
но и тем, что нормальные компоненты комплексного вектора О зависят только от касательных компонент вектора напряженности магнитного поля:
Ог = кхНу — куНх. (28)
ш
Рассматривая скачок нормальных компонент комплексного вектора О на границе раздела материальных сред, запишем
b" H"— b" H"
D"z - (Dz + D'z) = bXHy bvH*
2
Ш
кхНу — ку Нх + кх Ну — ку Нх\ (29)
ш ш ) '
Видно, что при выполнении условий (13) сопряжения касательных компонент волновых векторов и условий (16) непрерывности касательных компонент векторов напряженности магнитного поля правая часть соотношения (29) обращается в нуль. Нормальные компоненты
"комплексного" вектора Б не испытывают скачка на границе раздела двух материальных сред вследствие условия непрерывности касательных составляющих вектора напряженности магнитного поля.
Известно, что для определения векторного поля в безграничном пространстве необходимо и достаточно задать его источники — дивергенцию и ротор [7]. Для векторного поля в форме плоской гармонической бегущей волны (2) источниками поля являются правые части соотношений
%к •а = ф, %к х а = Д. (30)
Здесь ф и / — известные функции, структура которых необходимо согласована с выражением (2). Вычислим векторное произведение вектора %к на второе из соотношений (30):
%к х (%к х а) = %к х /. (31)
Раскрывая двойное векторное произведение и используя первое из соотношений (30), получаем
а =-^^-. (32)
к •к
В соответствии с определением (27) имеем
Г Д Т т гк(гк • Н) Н • (к •к)
%к • Б = 0; %к х Б =--Ц---?-. (33)
%ш %ш
Если учесть уравнение (4) и линейность материальных уравнений среды, то окажется, что источник векторного поля Б действительно пропорционален вектору напряженности магнитного поля:
т д Н(к•к)
%к х Б =--^-. (34)
%ш
Аналогичные построения для вектора Б с учетом того, что материальная среда является однородной, линейной и изотропной, с учетом закона сохранения электрического заряда, приводят к зависимостям
• г ,=; ^ -ч* %к х (%к х Н — ])
%к • Б = 0; %к х Б =---^. (35)
Видно, что источники векторных полей б и б различаются между собой, что влечет за собой и различие в физическом содержании и математической форме интегральных соотношений на границе раздела материальных сред.
Выводы. Использование комплексной диэлектрической проницаемости ё для проводящей среды в граничных условиях задачи о падении плоской гармонической электромагнитной волны на границу раздела двух однородных изотропных сред позволяет получить формулы Френеля. Вместе с тем свойство непрерывности нормальных компонент
к границе раздела сред присуще комплексным векторам б = ёе0Е. Это свойство является следствием непрерывности касательных компонент напряженностей электромагнитного поля (т.е. отсутствия поверхностных токов проводимости). Нормальные компоненты вектора б = ££0Е на границе раздела испытывают скачок, равный поверхностной плотности сторонних электрических зарядов а, при этом на поверхности раздела выполняется закон сохранения электрического заряда. Отмеченное выше условие непрерывности нормальных компонент комплексных векторов б содержит в себе и скачок нормальных составляющих вектора б, и закон сохранения электрического заряда. Существование поверхностной плотности а сторонних электрических зарядов на границе раздела сред может оказаться важным в определенных физических ситуациях.
СПИСОК ЛИТЕРАТУРЫ
1. Топтыгин И. Н. Современная электродинамика. Ч.2. Теория электромагнитных явлений в веществе: Учеб. пособие. - Москва-Ижевск: НИЦ "Регулярная и хаотическая динамика". - 2005. - 848 с.
2. С и в у х и н Д. В. Оптика: Учеб. пособие для вузов. - 2-е изд.- М.: Наука, Физматлит, 1985. - 752 с.
3. Ахманов С. А., Никитин С. Ю. Физическая оптика: Учебник. - М.: Изд-во Моск. ун-та. - 1998. - 656 с.
4. Никольский В. В., Никольская Т. И. Электродинамика и распространение радиоволн: Учеб. пособие для вузов. - М.: Физматлит, 1989. - 544 с.
5. М а т в е е в А. Н. Электродинамика. Учеб. пособие для вузов. - 2-е изд. - М.: Высш. шк., 1980.- 383 с.
6. Кугушев А. М., Голубева Н. С., Митрохин В. Н. Основы радиоэлектроники. Электродинамика и распространение радиоволн. Учеб. пособие для вузов. - М.: Изд-во МГТУ им. Н.Э. Баумана. - 2001. - 368 с.
7. К о ч и н Н. Е. Векторное исчисление и начала тензорного исчисления. - М.: Наука, 1965.-426 с.
Статья поступила в редакцию 15.01.2008