ЭМ волны в материальных средах
.pdf
Свойства вектора M определяем, исходя из правильного описания физических явлений. Так, например, токи, протекающие в прямом и обратном направлениях через сечение диэлектрического тела произвольно расположенной плоскостью, равновероятны. Поэтому из физических соображений должно выполняться следующее соотношение:
|
|
|
|
|
|
∫ |
j navndS = ∫ |
ρυ |
ndS = 0 . |
(1.16) |
|
S |
|
S |
|
||
Здесь n – нормаль к рассматриваемому сечению диэлектрического тела.
Сформулируем условия, которым должен удовлетворять вектор M . Подставим в
(1.16) вместо вектора ρυ его представление из (1.15). В соответствии с теоремой Стокса, должно выполняться равенство
c ∫ nrotMdS = c ∫ Mdl , |
(1.17) |
|
S |
L |
|
где L – контур, ограничивающий поперечное сечение S , а контурный интеграл в (1.17) берется по контуру, прилегающему к диэлектрическому телу с внешней стороны. Вектор
M вне диэлектрического тела, где отсутствует ток ρυ , равен нулю. Поэтому контурный интеграл в правой части равенства (1.17) обращается в ноль, что обеспечивает выполнение равенства (1.16) и правомерность представления (1.15).
В выражении (1.15) вектор M и наведенный ток связаны между собой прямой зависимостью. Из (1.14) следует: чем больше приложенное поле B , тем больше наведенный ток, значит, тем больше длина вектора M . Иными словами, векторы M и B тоже связаны прямой зависимостью. В линейном приближении можно считать, что вектор M прямо пропорционален вектору B . С учетом (1.15) уравнение (1.1а) в статическом случае, когда
∂ |
= 0 , перепишется в виде |
|
∂t |
|
|
rot (B − 4πM )= 0 . |
|
|
|
(1.18) |
|
|
Обозначим содержимое круглых скобок в (1.18) символом H : |
|
|
B − 4πM = H . |
(1.19) |
Поскольку, как следует из (1.19), векторы B и H пропорциональны друг другу, для симметрии уравнений Максвелла положили M ~ H . В качестве коэффициента пропор-
циональности выбрано число χm – магнитная восприимчивость вещества.
M = χmH . |
(1.20) |
12
Подставляя в (1.19) вектор M , представленный в виде (1.20), получаем |
|
B = (1 + 4πχm )H . |
(1.21) |
Введем обозначение |
|
μ =1 + 4πχm , |
(1.22) |
где μ – магнитная проницаемость среды, измеряемая величина. В вакууме и любой не-
магнитной среде μ =1. Из формулы (1.22) следует, что в вакууме χm = 0 , следователь-
но, M = χm H = 0 . Это лишний раз подтверждает, что вектор M , представленный в виде (1.20), обеспечивает выполнение условий (1.17), необходимых для адекватного описания равной вероятности наведенных токов двух взаимно противоположных направлений.
Свяжем между собой векторы B и H . В соответствии с (1.21) и (1.22)
B = μH . |
(1.23) |
Из (1.23) легко понять физический смысл введенного вспомогательного вектора H – это напряженность магнитного поля в вакууме или в любой немагнитной среде, где μ =1.
Расписывая подробно введенное обозначение (1.14) для ρυ с учетом (1.20) и (1.23), получаем следующее выражение для среднего значения плотности наведенного тока в статическом случае:
|
|
cχm |
rot B . |
|
|
ρυ = |
(1.24) |
||||
μ |
|||||
|
|
|
|
||
Эта часть наведенного тока отражает явление искривления в магнитной среде траекторий движущихся заряженных частиц под действием магнитного поля. В немагнитной
среде, где χm = 0 , такой ток отсутствует. Данное явление рассмотрено в статическом случае.
1.2.5. Усредненное значение плотности наведенного тока ρυ в динамическом случае
В динамическом случае, когда ∂∂t ≠ 0 , к наведенному току (1.24) должна доба-
виться часть, обусловленная высокочастотным изменением направления вектора электрического поля E . Найдем эту добавку. Временно введем в уравнения Максвелла сторонние токи и заряды. Запишем уравнение Максвелла (1.3) при наличии сторонних токов j st и зарядов ρst :
13
divD = 4πρst . |
(1.25) |
||
Учтем, что ток j st и заряд ρst связаны между собой уравнением непрерывности |
|
||
∂ρst |
= −divj st . |
(1.26) |
|
∂t |
|||
|
|
||
Продифференцируем уравнение (1.25) по времени t и перепишем результат с учетом
(1.26):
div ∂D = 4π ∂ρst = −4π div j st ,
∂t ∂t
откуда получаем
|
∂D |
+ 4π j |
st |
= 0 . |
div |
∂t |
|
||
|
|
|
|
Из равенства нулю дивергенции следует, что в круглых скобках стоит произвольный вихревой вектор. Мы можем выбрать его, как нам удобно. Полагаем
|
∂D |
+ 4π j st |
= c rot H . |
|
||||||
|
∂t |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
«Выведем из игры» сторонние токи и заряды, полагая, что |
j st = 0 . Тогда получим новое |
|||||||||
уравнение |
|
|
|
|
|
|
|
|
|
|
|
|
rotH = |
1 |
∂D |
. |
|
|
(1.27) |
||
|
|
|
|
|
||||||
|
|
|
c |
|
∂t |
|
|
|||
Перепишем еще раз уравнение (1.2): |
|
|
|
|
|
|
|
|
|
|
|
|
rotE = − |
1 |
∂B |
. |
|
(1.2) |
|||
|
|
|
|
|||||||
|
|
|
|
c ∂t |
|
|
||||
Исключив из (1.27), (1.2) вектор E , приходим к уравнению Даламбера |
||||||||||
|
|
H − |
εμ ∂2 H |
= 0 . |
(1.28) |
|||||
|
|
c2 ∂t2 |
||||||||
|
|
|
|
|
||||||
14
В вакууме, где ε =1, μ =1 это уравнение принимает вид
H − 1 ∂2 H = 0 , c2 ∂t2
где c имеет смысл скорости света. Сравнивая это уравнение с предыдущим, приходим к
выводу о том, что |
|
c |
есть фазовая скорость волны в среде. Перепишем еще раз урав- |
|||||||
|
εμ |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
нения (1.1) и (1.27): |
|
|
|
|
|
|
|
|
|
|
|
|
|
rotB = 1 ∂E |
+ |
4π |
|
, |
|
||
|
|
|
ρv |
(1.1) |
||||||
|
|
|
c ∂t |
|
|
c |
|
|
|
|
|
|
|
rotH = |
1 |
∂D |
. |
|
|
(1.27) |
|
|
|
|
|
|
|
|||||
|
|
|
|
c |
∂t |
|
|
|
||
Из уравнения (1.1) вычтем уравнение (1.27):
rot (B − H )= −1c ∂∂t (D − E)+ 4cπ ρv ,
откуда следует с учетом полученных ранее соотношений (1.12) и (1.23), что
rotH (μ −1)= −1 |
∂ |
E (ε −1)+ 4π |
|
. |
|
||||
ρv |
(1.29) |
||||||||
|
|||||||||
c ∂t |
|
|
c |
|
|
|
|||
Подставляя в (1.29) выражения (1.11) и (1.22) для ε и μ , находим |
|
||||||||
4πχmrotH + 1 ∂E 4πχe = |
4π |
|
. |
|
|
|
|||
ρv |
|
|
|
||||||
c ∂t |
c |
|
|
|
|||||
Отсюда следует:
ρυ = ∂∂t (χe E)+ cχμm rotB .
В соответствии с (1.8), χeE = P ; где вектор P , параллельный вектору E , назван вектором поляризации. Окончательно получаем
j nav = |
|
= ∂P |
+ cχm rot B . |
|
ρv |
(1.30) |
|||
|
|
∂t |
μ |
|
Первое слагаемое в правой части равенства (1.30) обусловлено изменением длины |
||||
и направления (от положительного к отрицательному) вектора E в пределах периода вы- |
||||
сокой |
|
частоты. |
||
15
Оно называется током поляризации. Заметим, что вихревая часть тока (второе слагаемое в правой части) тождественно совпадает с результатом, полученным для статического случая.
Таким образом, проведенное исследование позволило получить усредненное зна-
чение тока ρυ в виде двух слагаемых: тока поляризации
jpol = |
∂P |
= χe ∂E |
(1.31) |
|
∂t |
∂t |
|
и вихревого тока
j = |
cχm |
rotB . |
(1.32) |
|
|||
rot |
μ |
|
|
|
|
|
К этим двум наведенным токам следует добавить так называемый ток потерь, временно исключенный при нахождении токов (1.31) и (1.32).
|
|
jσ =σE . |
(1.33) |
|
Таким образом, окончательно получаем |
|
|
||
|
|
nav = χe ∂E |
+ cχm rotB +σE . |
|
|
j |
(1.34) |
||
|
|
∂t |
μ |
|
Соотношение (1.34) получило название обобщенного дифференциального закона Ома.
1.3.Разные формы записи уравнений Максвелла
Втеории электромагнитного поля можно привести примеры таких типичных задач, важных с практической точки зрения.
1.В вакууме имеются сторонние источники электромагнитного поля в виде элек-
трических или магнитных токов, распределенных с заданной плотностью j e st или j m st . Необходимо найти поле этих источников.
2. В неограниченной материальной среде с макроскопическими параметрами ε , μ ,σ , зависящими от координат пространства (ε =ε (r ), μ = μ(r ), σ =σ (r )), рас-
пространяется электромагнитная волна, порожденная сторонними источниками, удаленными на бесконечность. Необходимо найти закономерности распространения этой волны в обозримой области пространства.
16
3. В безграничной области в вакууме имеется неоднородность, представляющая собой диэлектрическое или магнитодиэлектрическое тело, занимающее ограниченный объем V . На него падает электромагнитная волна. Необходимо найти поле, рассеянное этой неоднородностью.
Для решения первой задачи мы записываем неоднородные уравнения Максвелла в вакууме при наличии сторонних токов и зарядов:
rotH = 1 ∂E |
+ |
4π |
j st ; |
(1.1б) |
||
c ∂t |
|
|
c |
|
|
|
rotE = −1 |
∂H |
; |
(1.2б) |
|||
∂t |
||||||
|
c |
|
|
|
||
divE = 4πρst ; |
|
(1.3б) |
||||
divB = 0. |
|
|
|
|
(1.4б) |
|
Вуравнениях (1.1б)-(1.4б) фигурируют векторы поля E и H . Параметры ε , μ ,σ отсутствуют, вернее, полагаются такими: ε =1, μ =1, σ = 0 .
Впроцессе решения второй задачи удобно иметь дело с уравнениями Максвелла,
вкоторых материальные параметры ε (r ), μ(r ), σ (r ) фигурируют в явном виде. Для
этого выполним следующие преобразования. Запишем усредненные уравнения Максвелла (1.1)-(1.4) с учетом полученных выше выражений для усредненного наведенного тока (1.34). Уравнение (1.1) с учетом всех наведенных токов принимает вид
rotB = |
1 |
∂E |
+ |
4π |
(χe ∂E |
+ cχm rotB +σE) . |
(1.35) |
|
c |
∂t |
|
c |
∂t |
μ |
|
Объединим первое слагаемое в правой части с первым слагаемым в скобках, а второе слагаемое в скобках перенесем в левую часть уравнения. При этом получаем
|
|
4πχm |
|
1 + 4πχe ∂E |
|
4π |
|
|||
rotB 1 |
− |
|
|
= |
|
|
|
+ |
|
σE . |
μ |
|
c |
∂t |
c |
||||||
|
|
|
|
|
|
|
||||
Учитывая введенные ранее обозначения (1.11), (1.12), (1.22) и (1.23), перепишем первые два уравнения Максвелла в виде
rotH = |
ε (r )∂E |
|
+ 4π |
σ (r )E ; |
(1.36) |
|||
c |
|
∂t |
||||||
|
|
|
c |
|
|
|||
rotE = − |
μ(r ) |
|
∂H |
, |
|
(1.37) |
||
|
|
|
||||||
|
|
|
c |
|
∂t |
|
|
|
17
Запишем эти уравнения для комплексных амплитуд, когда зависимость от времени описывается функцией exp{iωt} . Уравнение (1.36) принимает вид:
rotH = |
iωε (r ) |
E + |
4π |
σ (r )E . |
(1.38) |
|
c |
c |
|||||
|
|
|
|
Его можно свести к более компактной форме путем введения комплексной диэлектрической проницаемости, учитывающей потери в среде:
ε =ε′+iε′′,
где ε′ =ε (r ), ε′′ = − |
4πσ (r ) |
, i – мнимая единица. |
|
|
|
|||
ω |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
При этом уравнение (1.38) сводится к виду |
|
|
|
|||||
|
|
rotH = |
iωε (r ) |
E . |
(1.39) |
|||
|
|
|
|
|||||
|
|
|
|
|
|
c |
|
|
Аналогичную форму принимает и уравнение (1.36): |
|
|
|
|||||
|
|
rotE = − |
iωμ(r ) |
|
H . |
|
(1.40) |
|
|
|
c |
|
|||||
|
|
|
|
|
|
|||
При отсутствии потерь в среде ε =ε (r ). Решение задач распространения элек-
тромагнитных волн в неоднородных средах, как правило, начинают именно с этого случая. Исходными являются уравнения (1.39), (1.40) при условии, что ε =ε (r ).
Решение третьей задачи начинается с того, что мы ищем в вакууме поле, рассеянное неоднородностью, занимающую объем V . Это поле порождается наведенными токами, возникающими в магнитодиэлектрическом теле под действием электромагнитного поля падающей на него волны. Для простоты полагаем, что потери в объеме V отсутст-
вуют (σ = 0 ). Первое из усредненных уравнений Максвелла для векторов полей E и B в
точках пространства r V при условии σ = 0 имеет вид |
|
|
|||||||
rotB = |
1 ∂E |
+ |
4π ∂P |
+ |
4π cχm |
rotB . |
(1.41) |
||
c |
∂t |
c ∂t |
c μ |
||||||
|
|
|
|
|
|||||
Здесь в правой части уравнения фигурируют два электрических тока: ток электрической поляризации jpole = ∂∂Pt = χe ∂∂Et и вихревой ток магнитного происхождения
18
jrot = cχμm rotB . Во втором уравнении Максвелла (1.2) относительно векторов E и B
токи отсутствуют. Уравнения (1.41) и (1.2) образуют несимметричную систему уравнений Максвелла. От этой системы можно перейти к симметричной системе усредненных урав-
нений относительно векторов E и H , оставив в первом уравнении только электрический ток поляризации j e = jpole и введя во втором уравнении (1.2) магнитный ток. Процедура исключения вихревого тока из уравнения (1.41) описана в процессе преобразования уравнения (1.35). Слагаемое, содержащее jrot , перенесено из правой части в левую часть
уравнения, которая превратилась в rotH . Преобразуем теперь уравнение (1.2). Для этого представим его в виде
rotE = − |
1 ∂B |
= − |
μ ∂H |
= − |
μ −1 +1 ∂H |
= − |
1 ∂H |
− |
4π |
|
μ −1 ∂H |
|
|
|
. |
||||||||||||
|
c ∂t |
|
c ∂t |
|
c |
∂t |
|
c ∂t |
|
c |
|
4π ∂t |
|
Множитель в скобках в правой части уравнения назовем |
j m . Как следует из (1.22), |
||||||||||||
μ −1 |
= χm . Поэтому χm ∂H |
= |
∂M = jpolm = j m естественно назвать наведенным током |
|||||
4π |
∂t |
|
∂t |
|
|
|
|
|
магнитной поляризации аналогично тому, как |
jpole |
= |
∂P |
= χe ∂E = j e |
назван наведен- |
|||
|
|
|
|
|
|
∂t |
∂t |
|
ным током электрической поляризации. |
|
|
|
|
|
|||
|
В результате всех проведенных преобразований получены симметричные уравне- |
|||||||
ния Максвелла: |
|
|
|
|
|
|
|
|
|
|
|
rotH = 1 ∂E |
+ 4π |
j e , |
|
(1.42) |
|
|
|
|
c ∂t |
c |
|
|
|
|
|
|
|
rotE = −1 ∂H − |
4π |
j m . |
(1.43) |
||
|
|
|
c ∂t |
|
c |
|
|
|
В этих уравнениях токи j e и j m при определении рассеянного поля следует рас-
сматривать как сторонние токи. Форма записи этих уравнений такая же, как и при решении первой задачи, с той лишь разницей, что эти квазисторонние токи являются наведенными. Задача рассеяния решается в два этапа. На первом этапе наведенные токи определяются путем решения интегральных уравнений макроскопической электродинамики [1]. На втором этапе найденные токи используют для определения поля, рассеянного диэлектрической или магнитодиэлектрической неоднородностью.
19
2. УРАВНЕНИЯ МАКСВЕЛЛА В БЕЗГРАНИЧНОЙ СРЕДЕ ПРИ
НАЛИЧИИ СТОРОННИХ ТОКОВ И ЗАРЯДОВ
2.1.Уравнения Максвелла для комплексных амплитуд
Вуравнении Максвелла (1.39), (1.40) сторонние токи и заряды ( j st , ρst ) отсутствуют.
Если они отличны от нуля, уравнения Максвелла в однородной магнитодиэлектрической среде принимают вид:
rotH = 1 |
ε ∂E |
+ 4π j st ; |
(2.1) |
||||
c |
|
∂t |
|
|
c |
|
|
rotE = −1 |
μ |
∂H |
; |
(2.2) |
|||
∂t |
|
||||||
|
c |
|
|
|
|
||
div(εE )= 4πρst ; |
(2.3) |
||||||
div(μH )= 0 . |
|
|
(2.4) |
||||
Сторонние токи и заряды связаны уравнением непрерывности |
|
||||||
divj st |
+ |
∂ρst |
= 0 . |
(2.5) |
|||
∂t |
|||||||
|
|
|
|
|
|||
Векторы электрического и магнитного поля E (r ,t ), |
H (r ,t ), плотность тока |
j st (r ,t ) |
|||||
и плотность заряда ρst (r,t ) являются произвольными функциями времени |
t и про- |
||||||
странственных координат r . Уравнения (2.1)-(2.5) можно упростить, сводя их к уравнениям относительно компонент Фурье. Временную зависимость для векторных и скалярных функций, присутствующих в уравнениях (2.1)-(2.5), представим в виде интегралов Фурье:
∞ |
E (r ,ω)eiωt dω; |
(2.6) |
E (r ,t )= R e ∫ |
||
−∞ |
|
|
∞ |
H (r,ω)eiωt dω; |
(2.7) |
H (r ,t )= Re ∫ |
||
−∞ |
|
|
∞ |
j (r ,ω)eiωt dω; |
(2.8) |
j (r,t )= R e ∫ |
||
−∞ |
|
|
∞ |
ρ(r,ω)eiωt dω. |
(2.9) |
ρ(r ,t )= R e ∫ |
||
−∞
20
Используем представления (2.6)-(2.9) в уравнениях (2.1)-(2.4). Легко показать, что уравне-
ния относительно E (r ,ω), H (r,ω) принимают известную форму, тождественно сов-
падающую с уравнениями Максвелла, записанными с использованием метода комплексных амплитуд:
rotH (r ,ω)= iωε E (r,ω)+ |
4π |
j st (r ,ω); |
(2.10) |
||
c |
c |
|
|
||
rotE (r ,ω)= − |
iωμ |
H (r,ω); |
|
(2.11) |
|
|
|
||||
|
c |
|
|
|
|
div(εE (r ,ω))= 4πρst (r ,ω); |
|
(2.12) |
|||
div(μH (r ,ω))= 0 . |
|
|
(2.13) |
||
2.2. Уравнения Максвелла в алгебраической форме
Комплексные амплитуды E (r ,ω), H (r ,ω) и др., в свою очередь, можно разложить в интеграл Фурье по плоским волнам:
−∞ |
|
( |
|
|
) |
|
|
E (r,ω)= |
∞ |
E k ,ω e−ikr dk ; |
(2.14) |
||||
∫ |
|||||||
|
−∞ |
|
( |
|
) |
|
|
H (r ,ω)= |
∞ |
H |
|
k ,ω e−ikr dk ; |
(2.15) |
||
∫ |
|
||||||
|
−∞ |
|
( |
|
|
) |
|
j (r ,ω)= |
∞ |
j |
|
k ,ω e−ikr dk ; |
(2.16) |
||
∫ |
( |
||||||
|
−∞ |
|
|
|
) |
|
|
ρ(r ,ω)= |
∞ |
ρ k ,ω e−ikr dk . |
(2.17) |
||||
∫ |
|||||||
Подставляя E (r,ω), H (r ,ω), j (r ,ω), ρ(r ,ω) |
из (2.14)-(2.17) |
в уравнения Мак- |
|||||
свелла (2.10)-(2.13) и учитывая соотношения |
|
|
|
|
|
|
|
rot (E (ω, k )e−ikr )= e−ikr , E (ω, k ) = −i k , E (ω, k ) e−ikr ; div(D(ω, k )e−ikr )= ( e−ikr , D(ω, k ))= −i(k , D(ω,k ))e−ikr ,
где D(r,ω)=εE (r,ω),
получаем уравнения Максвелла в алгебраической форме [1]:
k , H (ω,k ) = −ωcε E (ω, k )+ 4cπ j st ; |
(2.1а) |
21
