Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ЭМ волны в материальных средах

.pdf
Скачиваний:
30
Добавлен:
23.02.2015
Размер:
579.82 Кб
Скачать

Свойства вектора 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 ,ω eikr dk ;

(2.14)

 

−∞

 

(

 

)

 

H (r ,ω)=

H

 

k ,ω eikr dk ;

(2.15)

 

 

−∞

 

(

 

 

)

 

j (r ,ω)=

j

 

k ,ω eikr dk ;

(2.16)

(

 

−∞

 

 

 

)

 

ρ(r ,ω)=

ρ k ,ω eikr dk .

(2.17)

Подставляя E (r,ω), H (r ,ω), j (r ,ω), ρ(r ,ω)

из (2.14)-(2.17)

в уравнения Мак-

свелла (2.10)-(2.13) и учитывая соотношения

 

 

 

 

 

 

 

rot (E (ω, k )eikr )= eikr , E (ω, k ) = −i k , E (ω, k ) eikr ; div(D(ω, k )eikr )= ( eikr , D(ω, k ))= −i(k , D(ω,k ))eikr ,

где D(r,ω)=εE (r,ω),

получаем уравнения Максвелла в алгебраической форме [1]:

k , H (ω,k ) = −ωcε E (ω, k )+ 4cπ j st ;

(2.1а)

21