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

Кураев АА, Попкова ТЛ, Синицын АК 2004

.pdf
Скачиваний:
383
Добавлен:
15.06.2014
Размер:
5 Мб
Скачать
☆

V

= 1

. Коэффициент при первых производных по времени определяет

ф

εa µa

 

 

 

пространственное затухание волнового процесса в среде с джоулевыми потерями.

Таким образом, из УМ непосредственно следует, что электромагнитные процессы имеют волновой характер и, следовательно, любое ЭМП представимо как некоторое наложение электромагнитных волн.

5.2. Электродинамические потенциалы

Установим далее формулы, выражающие напряженность электрическо-

го поля E иr индукцию магнитного поля B через электродинамические по-

тенциалы A (векторный) и Ф (скалярный). Для этого обратимся к системе уравнений Максвелла.

Четвертое уравнение Максвелла утверждает, что поле B чисто вихревое:

r

(5.8)

div B = 0.

Из (5.8) следует, что поле B следует определить как

 

B = rot A ,

(5.9)

так как div rot A = 0 .

 

Воспользуемся теперь вторым уравнением Максвелла

 

r

 

rot E = −

∂B

.

(5.10)

 

 

∂t

 

Подставляя в (5.10) выражение B (5.9) и имея в виду, что координаты и время в используемой лабораторной системе отсчета независимы, получим

r

 

 

 

 

 

 

∂A

= 0 .

(5.11)

rot E +

 

 

 

∂t

 

 

Из (5.11) с учетом того, что rot grad Ф = 0,следует

 

r

r

 

 

 

∂A

− gradФ.

(5.12)

E = −

 

∂t

 

 

 

31

Введенные таким образом функции координат и времени Ar(rr,t), Ф(rr,t)

(сrтрадиционным следованием знаков), как известно, и являются векторным r r

(A)и скалярным (Ф) потенциалами электромагнитного поля E, B . Заметим, что такое определение потенциалов является очевидным

образом неоднозначным. Действительно, поле E, B не изменяется, если r

вместо исходных A , Ф использовать A′, Ф′, определенные следующим образом:

r′

r

′

=Ф −

∂f (x, y, z,t)

.

(5.13)

A

= A + grad f (x, y, z,t), Ф

 

∂t

 

 

 

 

 

 

Таким образом, векторный потенциал определен с точностью до градиента произвольной функции f (x, y, z,t), а скалярный – с точностью до

производной той же функции по времени ∂f (x, y, z,t). Это означает, что

∂t

электродинамические потенциалы могут быть доопределены путем наложения на них любого условия связи, ненарушающего (5.13). Это весьма важный момент, позволяющий существенно улучшить вычислительную процедуру при решении электродинамических задач путем соответствующего выбора указанного условия.

Налагаемые на электродинамические потенциалы условия или связи называются условиями калибровки потенциалов. Из (5.13) очевидно, что

′

 

всегда можно f (x, y, z,t) выбрать так, что Ф (x, y, z,t)= 0 , (т.е. чисто

переменное электромагнитное поле E, B , когда ∂f

≠ 0 , может быть в

∂t

r′

 

принципе описано полностью одним векторным потенциалом A ).

Сформулируем теперь уравнения, определяющие электродинамические

 

 

 

 

r

 

 

 

 

потенциалы A , Ф. Для этого используем первое и третье уравнения

 

Максвелла

 

 

 

 

 

 

1

 

 

 

 

r

r

 

 

r

 

∂E

 

 

 

rot B = εa

 

+δ ,

(5.14)

 

µa

∂t

 

 

 

 

 

 

 

 

 

r

 

ρ

 

 

 

 

div E =

 

 

.

 

 

(5.15)

 

εa

 

 

 

 

 

 

 

 

 

 

Здесь εa , µa - соответственно диэлектрическая и магнитная проницаемости пространства, которые при записи (5.14), (5.15) считались постоянными, δ - плотность электрического тока, ρ - плотность электрического заряда.

32

Подставляя в уравнение (5.14) выражения (5.9) и (5.12), а в уравнение

(5.15) – (5.12), получим уравнения, определяющие A , Ф:

 

 

 

 

r

 

∂

 

 

 

r

 

 

∂2 A

 

 

r

 

rot rot A + εa µa

 

 

 

= −εa µa

 

 

grad Ф + µaδ ,

(5.16)

∂

t2

∂

t

 

r

 

 

 

 

 

 

 

 

2

∂A

 

ρ

 

 

 

 

 

Ф = −div

 

 

−

 

.

 

 

 

∂ t

 

εa

 

 

 

Наиболее простой вид принимает уравнения (5.16), (5.17), если использовать условие калибровки Лоренца:

divAr +εaµa ∂∂Фt = 0 .

При условии (5.18) получим

 

 

 

 

 

→

 

 

 

 

 

2

→

µ

∂2 А

= −µ

 

→

 

А−ε

a ∂t2

 

δ

 

 

a

 

 

 

a

 

 

2

 

 

 

∂2Ф

 

ρ

 

 

Ф−εaµa

 

= −

 

.

∂t2

εa

При записи (5.19) использовано векторное тождество:

→

→

→

rot rot A = grad div A− 2

A .

(5.17)

(5.18)

(5.19)

(5.20)

Однако в случае, когда необходимо выделение содержащей разрыв на границахr источников квазистатической составляющей электрического поля

E , можно воспользоватся кулоновской калибровкой потенциалов:

 

r

 

 

 

 

 

 

(5.21)

div A = 0.

 

 

 

 

 

В этом случае уравнения (5.16), (5.17) принимают вид

 

 

2 r

 

r

 

 

∂Ф

 

 

 

µa

∂2 A

 

r

r

(5.22)

A −εa

∂t2

= −µa

δ −εa

 

= −µаδ′,

 

 

 

 

 

∂t

 

 

2Ф = −

ρ

.

 

 

 

 

(5.23)

 

 

 

 

 

εа

 

 

 

 

33

Соответственно квазистатическая часть электрического поля Е1

определяется только через Ф, а динамическая Е2

r

- только через А :

r

r

∂Аr r r

r

 

E1 = − Ф,

Е2 = −

 

, E = E1

+ E 2 .

(5.24)

∂t

 

 

 

 

 

В дальнейшем мы будем использовать калибровку Лоренца (5.18), поскольку уравнения (5.19), (5.20) в этом случае наиболее просто связывают

электродинамические потенциалы с их источниками - δ и ρ .

5.3.Электричеcкий вектор Герца или поляризационный потенциал

Определим векторный потенциал через электрический вектор Герца Пe следующим образом:

A = εa µa ∂∂Пtе .

Скалярный потенциал Ф определим через Пe так, чтобы выполнялось условие калибровки Лоренца (5.18) :

εaµa ∂∂Фt = −div Ar = −εaµa ∂∂t div Пrе,

тогда Ф = −div Пrе.

Подчеркнем, что указанные определения возможны только при ∂∂t ≠ 0,

т.е. для чисто переменных полей. В общем случае, когда имеются статиче-

ские и стационарныеr

компоненты ЭМП, нужны все четыре составляющие

потенциалов, т.е. и Ar, и Ф.

 

r

определяются только через один вектор

При введении Пe поля E, H

Пre :

 

 

 

 

 

 

r

1

r

 

∂Пrе

 

 

H =

 

rot A = εarot

 

,

(5.25)

µa

 

 

 

∂t

 

 

34

r

r

 

 

 

 

 

2

 

e

r

 

∂A

− grad Ф = −ε

 

µ

 

∂

П

(5.26)

E = −

a

a

 

 

+ grad div Пе .

 

∂t

 

 

∂t2

 

 

 

Выясним физический смысл Пe . Для этого установим его источник. Обратимся к понятию поляризованности среды и вектору поляризации p .

r

r

= −qсвяз = −∫ ρсвязdV .

(5.27)

∫ pdS

S

 

V

 

Здесь qсвяз - связанный заряд, пересекающий S при поляризации.

Уравнение (5.27) можно переписать, используя теорему Остроградского – Гаусса:

r

∫div pdV = −∫ρсвязdV .

V V

Ввиду произвольности V получаем дифференциальный эквивалент

(5.27):

div prcm = −ρсвяз.

Если действуют сторонние источники, создающие рrcm и ними, очевидно, должна существовать идентичная (5.28) связь:

div prcm = −ρcm .

(5.28)

ρст , то между

(5.29)

Дифференцируя по времени (5.29) и, используя уравнение непрерывности, имеем

r

 

 

∂ρcm = divδrcm .

 

 

div ∂p

= −

(5.30)

 

∂t

 

 

 

 

∂t

 

 

 

 

Из уравнения (7.6) с учетом тождества div rot a ≡ 0 получаем

 

r

r

 

 

 

 

r

r

r

 

 

∂p

cm

 

 

 

 

δcm =

 

+ rot a

= δcmP

+δcm3 .

 

 

 

 

 

 

 

∂t

 

 

 

 

 

 

 

 

 

 

rP

 

 

 

r

 

 

 

 

 

 

 

 

∂pcm

 

 

3

 

Здесь δcm

=

 

 

 

- разомкнутая плотность тока, образующая ρcm , δст

-

 

∂t

 

 

 

 

 

 

 

 

 

 

замкнутая, вихревая плотность тока, для которой divδcm3

= 0.

 

Рассмотрим далее процессы, связанные только с δстP . Для этого случая

 

УМ принимают вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

35

r rot H

r rot E

 

 

r

 

 

r

 

= εa

∂E

+

 

∂Pcm

,

(5.31)

∂t

 

 

 

 

∂t

 

 

 

 

r

 

 

 

 

= −µa

∂H

.

 

(5.32)

∂t

 

 

 

 

 

 

 

Рассматриваемый случай может соответствовать, например, возбуждению ЭМП с помощью системы диполей, образованных разомкнутыми проводниками с возбужденным в них от источника сторонним током.

Подставим в (5.31) выражения H , E (5.25) и (5.26):

 

 

 

 

 

 

 

∂Пr

е

 

 

∂

 

 

 

 

 

∂2 Пr е

r

 

 

∂Рr

ε

a

rot rot

 

 

 

= ε

a

 

 

−ε

a

µ

a

 

 

+ grad div П

е

+

ст

.

 

 

 

 

2

 

 

 

∂t

 

 

 

 

 

 

 

∂t

 

 

 

∂t

 

 

 

 

 

 

 

∂t

 

 

 

 

 

 

 

 

Используя в левой части записанного уравнения тождество rot rot Пе = grad div Пе − 2 Пе , получаем

∂

 

 

 

 

 

∂2

r

е

 

r

 

 

2

r

е

 

П

 

pcm

 

 

−

 

П

 

+εа µа

 

 

 

−

 

 

= 0 .

 

 

 

∂t2

 

εa

∂t

 

 

 

 

 

 

 

 

В случае, который мы рассматриваем ∂∂t ≠ 0 , и поэтому окончательно имеем

 

 

r е

 

r

 

 

r

 

 

 

 

2

 

∂2 П

е

 

pcm

 

 

 

 

 

П

−εа µа

 

 

= −

 

.

 

(5.33)

 

∂t

 

εa

 

 

 

 

 

 

 

 

 

 

Полученное волновое уравнение для

Пe

(5.33) по структуре не отлича-

 

 

 

 

r

 

r

 

 

 

 

ется от уравнения для E,

H , A, Ф. Источником же Пe , как следует из (5.33),

является сторонняя поляризация pcm , откуда и название этого вектора – поляризационный потенциал. Простая связь Пe и pcm весьма полезна при решении конкретных задач. Например, если prcm = Z0 pcm (т.е. prcm имеет только Z – компоненту) из (7.9) очевидно, что для описания возбуждаемого ЭМП достаточно использовать только одну компоненту Пе − Z0 Пze . Т.е. в этом случае векторная задача сводится к скалярной! Заметим, что в этом случае

ЭМП имеет только пять компонент: в H отсутствует в соответствии с (5.25) продольная компонента Hz. Такое поле можно назвать поперечно – магнитным.

36

5.4. Фиктивные магнитные точки и заряды. Перестановочная двойственность УМ.

Магнитный вектор Герца

Обратимся теперь к замкнутым токам. Как известно из решения задач магнитостатики, магнитное поле рамки с током i соответствует полю маг-

нитного диполя, т.е. магнитный момент M замкнутого тока определяется как

Mr = iµa Snr,

где Snr - вектор-площадка, охватываемая током i. С другой стороны, можно представить, что

Mr = iµa Snr = rlqm ,

где qm – магнитные заряды диполя длиной l.

В свою очередь, im = dqdtm , где im – разомкнутый магнитный ток.

Таким образом, замкнутый ток i создает поле, подобное тому, которое возбуждает разомкнутый магнитный ток im. Такие эквивалентные qm и im называются фиктивными магнитными зарядами и токами. В ряде случаев такая замена эффективна, поскольку позволяет не решать непосредственно поставленную задачу, а воспользоваться уже решенной с соответствующими переобозначениями.

При введении фиктивных магнитных токов и зарядов УМ переписываются следующим образом:

r

 

 

r

 

 

 

 

∂E

 

 

 

 

rot H =εa

,

 

(5.34)

∂t

 

 

 

 

 

 

r

 

 

∂Hr

r

 

rotE

= −µa

 

 

 

 

−δm ,

(5.35)

 

∂t

 

 

 

 

 

 

 

r

 

 

 

 

 

 

(5.36)

div D = 0 ,

 

 

 

 

 

r

= ρm .

 

 

 

div B

 

 

(5.37)

Очевидно от системы (5.34)...(5.37) можно перейти к исходным УМ путем следующих переобозначений:

r r

r r

r

ρm → −ρ .

E → H ,

H → E, εa → −µa ,

µa → −εa , δm → −δ ,

37

Это свойство называется перестановочной двойственностью УМ. Это свойство позволяет, как указывалось выше, не решать ряд задач непосредственно, а воспользоваться переобозначениями в решении, относящемся к “двойственной” задаче. В частности, очевидно, что задачу возбуждения ЭМП замкнутыми токами решать нет смысла как самостоятельную, т.к. она формулируется как задача возбуждения ЭМП разомкнутыми магнитными токами, т.е. в виде двойственной уже решенной задаче возбуждения ЭМП разомкнутыми электрическими токами.

Для разомкнутых магнитных токов, как и в электрическом случае, имеет место связь

rP

∂Jrcm

 

 

 

 

 

 

δm =

 

,

 

 

 

 

 

∂t

 

 

 

 

 

 

 

 

 

 

 

 

где Jrcm - намагниченность среды, создаваемая δmP .

 

Второе УМ при этом записывается в виде

 

r

 

 

r

 

r

 

 

 

∂H

 

∂Jcm

 

 

rot E

= −µa

−

.

(5.38)

∂t

 

 

 

 

 

∂t

 

ца Пr тПо: аналогии с “электрической” задачей введем магнитный вектор Гер-

r

 

 

 

 

∂Пr т

 

 

E

= −µ

a

rot

 

 

,

(5.39)

 

 

 

 

 

∂t

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

∂

2 r т

 

r

 

H = −εa µa

П

+ grad div Пт .

(5.40)

 

∂t 2

 

 

 

 

 

 

 

 

Подставляя (5.39), (5.40) в (5.38), получаем при ∂∂t ≠ 0 следующее вол-

новое уравнение для Пr т :

 

2

r

т

 

∂2 Пr т

 

Jcm

 

 

 

 

П

 

−εа µа

 

= −

 

.

(5.41)

 

 

∂t 2

µa

 

 

 

 

 

 

 

 

Структура уравнения (5.41) полностью идентична (5.33), что непосредственно может быть использовано при решении двойственных задач путем переобозначений.

38

Опять можно отметить, что если специфика задачи такова, что при Jrcm = Zr0 Jcm , достаточно использовать только одну компоненту Пr т =Z0 Пzm . Компонентой Пzm определяются в соответствии с (5.39) ЭМП, в которых отсутствует компонента Ez. Такие поля называются поперечно-электрическими.

5.5. Граничные условия для Пzm и Пze на идеально

проводящих продольных и поперечных поверхностях

Для большой группы задач, относящихся к регулярным волноводам и резонаторам, ЭМП может быть представлено в виде суперпозиции попереч- но-электрических (ТЕ) и поперечно-магнитных (ТМ) полей (Самарский А.А., Тихонов А.Н. «О представлении полей в волноводе в виде суммы полей ТЕ и ТМ, ЖТФ 1948 – т.18, №7 – с.959-970»). Поскольку эти поля выражаются со-

ответственно через Пzm и Пze, важнейшее значение приобретает формули-

ровка граничных значений (ГУ) для этих функций на идеально проводящих продольных и поперечных стенках указанных устройств.

Введем обобщенно-цилиндрическую систему координат q1,q2 ,Z , где q1,q2 - обобщенные криволинейные в общем случае координаты в попереч-

ной к Z плоскости Q3.

Координатные поверхности в этой системе определяются следующим образом :

Q1 : q1 = Const; Q2 : q2 = Const; Q3 : Z = Const .

Соответственно направляющие вектора системы определяются как

er1 dqr1 Q1, er2 dqr2 Q 2 , Z0 .

ГУ для Пze .

Будем исходить из того, что тангенциальная составляющая E на идеально проводящей поверхности S равна нулю. Таким образом,

Ez

Q , Q

2

= 0, Eq , q

2

Q

= 0.

 

1

1

 

3

Запишем E , выраженное через Z0 Пze :

r

 

µ

 

∂2 Пe

r

 

 

∂Пe

E = −ε

a

a

 

z

Z

0

+ grad

z .

 

 

∂t

2

 

 

 

 

 

 

 

 

 

 

 

∂ z

39

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂

 

 

 

 

 

 

 

 

 

Полагая процессы гармоническими

 

 

 

 

 

= jω , имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r

 

 

 

 

 

∂ Пe

k 2

=ω2ε

 

µ

 

 

,

 

 

 

 

 

 

E = k 2Z

0

Пe + grad

 

 

z ,

a

a

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂Пe

r 1

 

∂

 

 

∂Пe

 

r

1 ∂

 

 

 

 

∂Пe

r

 

∂2 Пe

 

grad

 

 

z

= e

 

 

 

 

 

z

+ e

 

 

 

 

 

 

 

 

 

z

+ Z

 

 

z

,

 

 

 

 

∂q

 

2 h

 

∂q

 

 

 

 

 

 

∂ z

1 h

 

 

 

∂ z

 

 

 

 

2

 

 

∂ z

 

0

∂ z

2

 

 

 

 

 

1

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

h1, h2 – соответственно метрические коэффициенты Ламэ обобщенных

координат q1, q2.

Таким образом, имеем

 

 

 

 

 

 

∂2

 

Пe

 

 

 

 

 

 

 

 

 

E

z

Q1

, Q2

= k 2

+

 

 

 

Q1 , Q2

= 0 .

 

2

 

 

 

∂ z

 

z

 

 

 

 

 

 

 

 

 

 

 

 

Это условие должно выполняться при любых Z, q1 или Z, q2, следовательно, необходимо потребовать

Пze

 

Q , Q

= 0 .

(5.42)

 

 

 

1

2

 

Иначе говоря, на идеально проводящих боковых стенках Q1, Q2 Пze об-

ращается в ноль.

Условие на поперечной стенке Q3 имеет вид

 

 

 

 

 

 

 

 

1

 

 

 

 

∂

 

∂Пe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Eq1 , q2

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

Q3 = 0.

 

 

Q3

= h

 

 

 

∂q

 

∂ z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1, 2

 

 

 

 

1, 2

 

 

 

 

 

Таким образом, если

∂

 

 

≠ 0 , получаем

∂q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1, 2

 

 

 

 

 

∂Пe

 

 

 

 

=

∂Пe

 

 

 

= 0 .

 

 

 

 

 

 

(5.43)

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

Q

 

rz

Q

 

 

 

 

 

 

∂ z

 

3

 

∂ n

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ГУ для Пzm .

 

 

 

 

 

 

 

 

 

 

 

 

 

v

Будем исходить из ГУ для нормальной составляющей H на идеально

проводящей поверхности:

 

 

 

 

 

 

 

Hq

 

Q

= 0, Hq

2

 

Q

= 0, H z

 

Q

= 0 .

 

 

 

1

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

3

40