Конспект лекций Электродинамика и РРВ
.pdfБаланс энергии будет активным, когда РΣ > 0, т.е. отдача энергии со внешнего пространства преобладает; при этом dW/dt + P < 0. В случае чистого излучения может оказаться, что внутренний запас энергии остается постоянным: W = const, тогда, как видно РΣ = -P. Поскольку РΣ > 0, то P < 0: излучение создается сторонними силами в V. (P0 = Pп’+ Pст и в случае отсутствия потерь РΣ = -P). Но возможно также, что P = 0 (нет ни сторонних сил, ни внутренних потерь, либо, они взаимно уравновешены, тогда РΣ = -dW/dt, а поскольку РΣ > 0, то dW/dt < 0. Это означает, что излучение обусловлено убыванием запаса энергии в V.
Варианты (в), (г) и (д) соответствуют условию: РΣ = 0, это нейтральный баланс энергии. Поток энергии в данном случае может проходить насквозь (в), так, что число входящих линий вектора Пойнтинга равно числу входящих; он также может не входить в область V (г) или вообще отсутствовать (д). Наконец, возможен пассивный баланс, когда поглощение преобладает над излучением (е, ж). При чистом поглощении (е) и постоянстве внутреннего запаса энергии РΣ = -P . Если же Р = 0, то РΣ = -dW/dt . Поскольку РΣ < 0, то dW/dt > 0: поглощение внешнего излучения приводит к росту запаса энергии.
Энергия электромагнитного поля
Исходя из равенства
dW |
= |
V∫ |
|
∂B |
+ E |
∂D |
(1) |
|
H |
|
dV |
||||
dt |
|
|
∂t |
|
∂t |
|
можно путем интегрирования определить энергию поля.
Если параметры среды ε и µ не зависят от времени, и среда изотропна,
E ∂D ∂t = ε0εE ∂E ∂t = ∂(ε0εE22) ∂t , H ∂B ∂t =µ0µH ∂H ∂t = ∂(µ0µH22)∂t.
Это значит, что операцию дифференцирования по времени в (1) можно вынести за знак интеграла.
В результате знак энергии в объеме V выражается следующим образом:
W = |
1 ∫(ε0εE2 +µ0µΗ2 )dV = |
1 ∫(ED + HB)dV . |
(2) |
|
2 V |
2 V |
|
как видим, энергия слагается из двух частей, одна из которых связана с электрическим полем, а другая – с магнитным. Поэтому пишут W = WЭ + WМ, разла-
гая электрическую энергию
W Э = |
ε0 |
∫εE2dV = |
1 ∫E DdV |
(3) |
|
2 |
V |
2 V |
|
и магнитную энергию
W М = |
µ0 |
∫µH2dV = |
1 ∫H BdV . |
(4) |
|
2 |
V |
2 V |
|
Подынтегральное выражение в (2) – это плотность энергии электро-
магнитного поля:
W = lim |
∆W |
= |
1 |
(ε0εE2 |
+µ0µH2 )= |
1 |
(E D + H B). |
(5) |
∆V →0 |
∆V |
|
2 |
|
|
2 |
|
|
Соотношение имеет смысл плотной электрической и магнитной энергии:
W = WЭ + WМ.
Локальный баланс и движение энергии
Если допустить, что поток вектора Пойнтинга РΣ через любую, не только замкнутую поверхность, представляет собой поток энергии через эту поверхность, то П следует понимать как плотность потока энергии
Π = lim |
∆PΣ π0 . |
(6) |
∆S→0 |
∆S |
|
здесь π0 - единичный вектор, указывающий направление движения энергии, ∆S
– ортогонально проецированная площадка, ∆РΣ - количество энергии, проходящей через ∆S за единицу времени.
На основании интегрального соотношения
∫ ΠdV = −∫∂w dV − ∫j EdV |
|
||||
V |
V ∂t |
V |
|
||
можем записать |
|
∂w |
|
|
|
Π + |
|
+ p = 0. |
(7) |
||
|
∂t |
||||
|
|
|
|
Если P = 0, то (7) совпадает по структуре с дифференциальной формулировкой закона сохранения заряда. Полной аналогии нет, потому что в отличии от заряда q, энергия электромагнитного поля не сохраняется: она переходит в другие виды энергии, порождается ими.
Равенство (7) есть уравнение баланса энергии в дифференциальной фор-
ме. Оно характеризует баланс энергии. Если в исчезающе малой окрестности некоторой точки баланс активен, то dW/dt + P < 0 и П > 0. При пассивном балансе dW/dt + P > 0 и П < 0, а при нейтральном dW/dt + P = 0 и П = 0. Таким образом, при активном балансе рассматриваемая точка является истоком линий вектора Пойнтинга, при пассивном балансе – стоком, а при нейтральном – лежит на некоторой линии вектора Пойнтинга.
Сторонние силы в уравнениях Максвелла
Описание неэлектромагнитных факторов, как говорят, сторонних сил в большинстве случаев сводится к изменению вида материальных уравнений
j = σE.
Используем одну из следующих
j = σ(E + Eст ),
(8)
j = σE + jст .
здесь Ест и jст при решении электродинамических задач являются заранее заданными. Величина Ест называется напряженностью сторонних сил (или про-
сто сторонней напряженностью), а jст – сторонним током. Выражение для плотности мощности р = jЕ можем теперь представить в следующем виде:
p = σ-1j2 − jEст; p =σE2 + jстE. |
(9) |
Таким образом, можно записать: |
|
p = pп + pст , |
(10) |
где рп = σ-1j2 = σE2 - характеризует потери энергии электромагнитного процесса, pст = -jEст = jстE - характеризует действие сторонних сил.
Сторонние силы обычно локализованы и равны нулю вне объема Vи –
области источников.
Лекция № 5. Гармонические колебания в электродинамике
Гармонические колебания и комплексные амплитуды
Если векторная величина u(t) изменяется во времени по закону |
|
u(t ) = um cos(ωt +ϕ), |
(1) |
то говорят, что происходят гармонические колебания этой величины. При этом um называется амплитудой; ω – круговой частотой, а аргумент косинуса ωt + ϕ - фазой (полной фазой); ϕ - начальной фазой колебаний.
Гармоническое колебание – переходной процесс. Периодом Т называется наименьший отрезок времени, обладающий тем свойством, что u(t + T) = u(t). Очевидно
ω(t +T )+ ϕ = ωt + ϕ+ 2π; T = |
2π |
= |
1 |
, |
(2) |
|
ω |
f |
|||||
|
|
|
|
где f – частота колебаний, число периодов в секунду.
В теории электромагнетизма встречаются скалярные и векторные функции координат и времени, описывающие гармонические колебания. Векторная функция V(r, t) в общем случае разлагается на три скалярных в выбранной системе координат. Например:
V |
( |
r,t |
) |
0 mx ( |
r |
) |
|
|
x ( |
|
) |
0 my ( |
r |
) |
|
y ( |
) |
|
|||
|
|
= x V |
|
|
cos ωt + ϕ |
|
|
r + y V |
|
cos ωt + ϕ |
|
r + |
(3) |
||||||||
|
|
|
|
+ z V |
|
r |
) |
cos ωt +ϕ |
z ( |
r . |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
0 mz ( |
|
|
|
|
) |
|
|
|
|
|
|
|
Здесь амплитуды и начальные фазы скалярных компонент являются функциями координат.
Если компоненты вектора имеют одинаковые начальные фазы, то эта запись координат имеет вид:
( |
r,t |
) |
m ( |
) |
|
( |
) |
(4) |
V |
|
=V |
r |
cos ωt +ϕ |
r , |
где Vm = x0Vmx + y0Vmy + z0Vmz; ϕ = ϕx = ϕy = ϕz .
В теории гармонических колебаний обычно применяется метод комплексных амплитуд, суть которого состоит в том, что вместо тригонометрических функций в выражениях типа (1), (2), (4) употребляются экспоненциальные. При этом получаются комплексные представления величин. Например, вместо
(1) запишем:
u = umei(ωt+ϕ) = umeiωt . |
(5) |
Точка над величиной говорит о ее комплексном представлении. Величина u = umeiϕ , несущая информацию об амплитуде и начальной фазе, называется
комплексной амплитудой. В силу известной формулы Эйлера физическая величина и есть вещественная часть ее комплексного представления:
u = Reu = Reumeiωt . |
(6) |
Отметим, что если как в (4) амплитуда и фаза являются функциями координат, то комплексное представление (5) есть произведение функции координат um (r) и функции времени exp(iωt).
Из формулы Эйлера вытекает следующее соотношение
u = 12 (u +u ), |
(7) |
где звездочка означает комплексное сопряжение.
Метод комплексных амплитуд значительно упрощает технику преобразований при получении решения дифференциальных уравнений в частных производных.
Все члены линейного дифференциального уравнения оказываются умноженными на exp(iωt). Опуская этот множитель, получаем уравнения относительно комплексной амплитуды, не зависящей от времени. Если в результате решения комплексная амплитуда отрицательна, то для получения искомой физической величины надо лишь умножить комплексную амплитуду на exp(iωt) и отделить вещественную часть.
Для периодической функции от t средним значением называется деленный на Т интервал от 0 до Т. Очевидно, что среднее значение гармонической функции и равно нулю. Среднее от квадрата гармонически колеблющийся величины есть
___ |
1 |
|
T |
2 |
T |
|
|
|
|
|
|
|
2 |
1 |
T |
|
||||
u2 = |
|
u |
2dt = um |
cos2 (ωt +ϕ)dt = um |
cos2 |
(ωt +ϕ)+1 dt = |
||||||||||||||
|
|
|
|
|||||||||||||||||
|
|
T |
∫0 |
T |
∫0 |
|
|
|
|
|
|
T 2 |
∫0 |
|
||||||
|
u2 |
|
|
1 |
|
|
|
|
|
|
|
|
1 |
|
|
1 |
|
|
(8) |
|
|
|
sin (x) |
|
2(ωT +ϕ) |
|
T |
|
2 |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|||||||||||||
= |
m |
|
|
|
|
|
2ϕ |
+ x |
0 |
|
= |
|
um |
= |
|
umum . |
|
|||
|
|
|
|
|
|
|
||||||||||||||
|
2T 2ω |
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Введем функцию v = vmcos(ωt + ψ) и найдем среднее от произведения uv:
uv = 1 |
T |
1 |
T |
|
∫uvdt = |
∫(u +u )(v + v )dt = |
|||
___ |
|
|
|
|
|
T |
0 |
2T |
0 |
41T T∫(umvmei2ωt + umvme−i2ωt + umvm + umvm )dt.
0
Первые два члена, выражающие гармонические колебания с частотами 2ω и -2ω дают при интегрировании нуль. В результате
uv = 14 (umvm + umvm )= 14 (umeiϕvme−iψ + ume−iϕvmeiψ )=
= |
1 umvm (ei(ϕ−ψ) + e−i(ϕ−ψ) )= |
1 umvm cos(ϕ− ψ)= |
(9) |
|
4 |
2 |
|
12 Re(umvm )= 12 Re(umvm ).
Подобные формулы для векторных гармонических величин v и w имеют
вид:
___ |
1 v v ; |
|
|
|
|
v2 = |
|
|
|
||
|
2 |
|
|
|
|
______ |
= 1 Re(vm wm )= |
1 Re(vm wm ); |
|||
v w |
|||||
|
|
2 |
|
2 |
|
______ |
|
1 Re |
(vm ×wm ) |
|
1 Re(vm ×wm ). |
v ×w = |
= |
||||
|
|
2 |
|
|
2 |
Рассмотрим разложение некоторой периодической функции в ряд Фу-
рье:
|
a0 |
∞ |
∞ |
|
|
u (t )= |
+ ∑an cos(nωt )+ ∑bn sin (nωt ), |
(10) |
|||
|
|||||
2 |
n=1 |
n=1 |
|
которое представим в комплексной форме
∞ |
1 |
T / 2 |
|
|
u (t )= ∑ cneinωt , cn = |
∫ u (t )e−iωt dt , |
(11) |
||
T |
||||
n=−∞ |
−T / 2 |
|
где
c−n = cn = 12 (an −ibn ).
Видно, что коэффициенты ряда Фурье (11) – не что иное, как комплексные амплитуды, а члены – комплексные представления гармонических колеба-
ний с частотами nω (…, -2ω, -ω, 0, ω, 2ω, …).
В случае произвольной временной зависимости запишем разложение в интеграл Фурье
u (t )= |
1 |
∞∫u (ω)eiωt dω; u (ω)= |
1 |
∞∫u (t )e−iωt dt. |
(12) |
|
|
||||
|
2π −∞ |
2π −∞ |
|
Комплексные проницаемости
Уравнения Максвелла для гармонически меняющихся полей имеет вид:
×H |
m |
= iωε |
εкE |
m |
+ jст; |
|
|
|
0 |
|
m |
(13) |
|||
×Em = iωµ0µкHm , |
|||||||
|
где εк = ε −i ωεσ0 - комплексная диэлектрическая проницаемость.
В рамках метода комплексных амплитуд любой из параметров уравнения Максвелла должен рассматриваться в пределах комплексной плоскости, что приводит к расширению физического содержания некоторых понятий. В первую очередь это относится к проницаемостям, которые будем обозначать в виде
′ ′′ |
′ ′′ |
(14) |
ε = ε −iε ; |
µ =µ −iµ . |
Таким образом теперь может быть описана инерционность поляризации диэлектриков. При гармонических колебаниях достаточно ввести фазовое за-
паздывание D по отношению к E, т.е. Dm = ε0εEme−iα . Но это означает, что в
данном случае диэлектрическая проницаемость – комплексная величина ε(cosα- isinα).
В дополнении к (14) вводят следующие обозначения:
tg(∆) = ε′′ ε′, tg(∆м )=µ′′ µ′, |
(15) |
где ∆ называют углом электрических потерь (или просто углом потерь, а ∆м - углом магнитных потерь).
Ввиду этого выражениям (15) можно придать другую форму
′ |
( |
−i tg |
( |
)) |
|
|
e |
−i∆ |
, |
|
||||
ε = ε 1 |
∆ = |
ε |
|
|
(16) |
|||||||||
′ |
( |
|
( |
)) |
|
|
|
|
|
−i∆ |
||||
1−i tg ∆ = |
µ |
e |
. |
|||||||||||
µ =µ |
|
|
|
Получим уравнения электродинамики второго порядка в комплексной форме
×(ε−1 ×Hm )− ω |
2 |
µHm = ×(ε−1jстm ), |
(17) |
||
c |
|
|
|
|
|
×(µ−1 ×Em )− ω |
2 |
εEm = −iωµ0 jстm . |
(18) |
||
c |
|
|
|
|
Пусть среда однородна (ε = const, µ = const). Вынося обратные проницаемости за знак операций дифференцирования и применяя слева тождество
× × v = ( v)− 2 v,
получаем
2 |
|
ω 2 |
ст |
, |
(19) |
|
|
Hm + |
|
εµHm = − × jm |
|||
|
|
c |
|
|
|
|
2 |
|
ω 2 |
i |
ст |
ст |
|
||
|
Em + |
|
εµEm = |
|
( jm )+ iωµ0µjm . |
(20) |
||
ωε0ε |
||||||||
|
|
c |
|
|
|
|
Это так называемые неоднородные уравнения Гельмгольца.
Векторный потенциал
Поскольку B = 0 (четвертое уравнение Максвелла), то вектор магнитной индукции можно представить в виде ротора некоторого вспомогательного вектора А
B = × A или H = (µ0µ)−1 × A. |
(21) |
Подстановка этого выражения Н во второе уравнение Максвелла, приводит к равенству
|
∂A |
|
= 0. |
(22) |
× E + |
|
|||
|
∂t |
|
|
|
Мы видим, что векторная функция в скобках является потенциальной. Приравнивая эту функцию градиенту ϕ, получаем
E = − ϕ − |
∂A |
. |
(23) |
|
|||
|
∂t |
|
Таким образом, напряженности Е и Н выражены при помощи соотно-
шений (21) и (23) через электродинамические потенциалы А и ϕ.
В комплексной форме соотношения (22) и (23) записываются в виде:
Hm = (µ0µ)−1 × Am , |
(24) |
Em = − ϕm −iωAm . |
(25) |
Внося их в первое уравнение Максвелла, записываем в случае однородной среды:
|
|
|
ω 2 |
|
|
|
|
εµ |
|
ст |
|
× × Am − εµAm = −iω |
c |
2 |
ϕm +µ0µjm . |
||||||||
|
|
|
c |
|
|
|
|
|
|
|
|
2 |
|
ω 2 |
|
|
|
εµ |
ϕm |
|
ст |
||
|
Am + εµAm |
= iω |
c |
2 |
+ Am |
−µ0µjm . |
|||||
|
|
c |
|
|
|
|
|
|
|
|
|
Налагая дополнительное условие |
|
|
|
|
|
|
|
||||
|
|
|
iω |
εµ |
ϕm + Am = 0, |
(26) |
|||||
|
|
|
|
c2 |
|
|
|
|
|
|
|
которое иногда называют лоренцевой калибровкой, получаем следующие уравнения Гельмгольца для электродинамических потенциалов
2 |
|
|
ω 2 |
ст |
|
|
|
|
|
||||
|
Am + |
εµAm = −µ0µjm . |
|
|
|
||||||||
|
|
|
c |
|
|
|
|
|
|
|
(27) |
||
|
|
ω |
2 |
|
|
|
|
|
|||||
2 |
ϕm + |
|
|
−1 |
ст |
|
|||||||
|
|
|
εµϕm = −i(ωε0ε) |
|
|
jm . |
|
||||||
|
|
c |
|
|
|
|
|
|
|
|
|||
Посредством (26) можно исключить из (25) скалярный потенциал, в ре- |
|||||||||||||
зультате чего |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
2 |
|
|||
|
|
|
|
|
|
( Am )+ ω |
|
|
|
||||
Em = −i |
|
|
|
|
|
εµAm . |
(28) |
||||||
|
ωεµ |
|
|||||||||||
|
|
|
|
c |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, комплексные амплитуды напряженностей поля выражены при помощи формул (26) и (28) только через векторный потенциал.
Лекция № 6. Гармонические колебания в электродинамике (продолжение)
Баланс энергии при гармонических колебаниях
Поскольку гармонические колебания электромагнитных полей, представляющие интерес для радиоэлектроники, являются весьма быстрыми, обычно имеют дело с их усредненными во времени энергетическими характеристиками.
Средняя плотность энергии в предположении безынерционной среды
w |
= |
1 |
(ε0εEEm +µ0µHHm ). |
(1) |
|
|
4 |
|
|
Среднее значение плотности мощности
p = Re p; p = |
1 j Em . |
(2) |
|
2 |
|
Величина p называется плотностью комплексной мощности, а сама комплексная мощность P , есть интеграл от p по объему.
___
Среднее значение Π вектора Пойнтинга
___ |
|
1 |
(Em ×Hm ). |
|
Π |
= ReΠ; Π = |
(3) |
||
|
|
2 |
|
|
Промежуточная величина Π называется комплексным вектором Пойн-
тинга. Поток Π через некоторую поверхность S называют комплексным пото-
ком энергии.
Выведем уравнение для среднего баланса энергии электромагнитного поля в некоторой области V.
Будем исходить из комплексной формы уравнений Максвелла, записывая первое из них комплексно-сопряженным:
×Hm = −iωε0ε*Em + (jстm ) , |
(4) |
|
×Em = −iωµ0µHm . |
||
|
||
Все члены первой строки умножим на Em , а второй – на Hm . Проведем |
||
вычитание соответствующих частей |
|
|
Em ×Hm − Hm ×Em = |
(5) |
|
= −iω(ε0ε*Em Em −µ0µHm Hm )+ (jстm )* Em |
||
|
||
Используя известное соотношение векторного анализа, получим |
|
|
− (Em ×Hm )= −iω(ε0ε*Em Em −µ0µHm Hm )+ (jстm )* Em . |
|
Изменив знаки и используя соотношения (1) – (3), перепишем это выражение в виде
|
|
Π = i |
ω(ε0ε*Em Em −µ0µHm Hm )− pcт . |
|
|
|
(6) |
|||||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
Это – комплексный аналог полученного ранее уравнения баланса энер- |
||||||||||||||||
гии в дифференциальной форме для линейных значений поля |
|
|
|
|
||||||||||||
|
|
|
Π = − |
|
∂D |
+ H |
∂B |
− j E. |
|
|
|
|
|
|||
|
|
|
E |
∂t |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
∂t |
|
|
|
|
|
|
|
|
Внесем |
в |
(6) |
|
представления |
комплексных |
|
|
проницаемостей |
||||||||
′ ′′ |
′ |
′′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ε =ε −iε ; µ =µ |
- iµ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Π = i |
ω |
|
′ |
′′ |
|
|
′ |
′′ |
|
|
|
cт |
. |
||
|
2 ε0 |
(ε + iε )Em Em −µ0 (µ −iµ |
)Hm Hm |
− p |
|
|||||||||||
и выделим в уравнении действительную и мнимую части |
|
|
|
|
|
|||||||||||
|
|
ReΠ = − |
ω(ε0ε′′Em Em +µ0µ′′Hm Hm )− Re pcт , |
|
|
|||||||||||
|
|
|
|
ω |
2 |
|
|
|
|
|
|
|
|
|
|
(7) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Im Π = |
′ |
|
|
|
′ |
|
− Im p |
cт |
. |
|
|
|
||
|
|
2 |
(ε0ε Em Em |
−µ0µ Hm Hm ) |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Проведем интегрирование по некоторому объему V с границей S и перейдем к следующим соотношениям:
Re ∫Π dS= − |
ω∫(ε0ε′′Em Em +µ0µ′′Hm Hm )dV − Re Pcт , |
(8) |
||||||||
S |
|
2 V |
|
|
|
|
|
|
|
|
Im ∫Π dS = |
ω |
|
′ |
|
′ |
|
cт |
, |
(9) |
|
|
|
∫(ε0ε Em Em |
−µ0µ Hm Hm )dV − Im P |
|
||||||
S |
2 V |
|
|
|
|
|
|
|
||
где Pст - интеграл от рст |
|
по V, выражающий комплексную мощность источни- |
||||||||
ков. |
|
|
|
|
|
|
|
|
|
|
Обсудим смысл первого из полученных равенств. В левой части – вещественная часть комплексного потока энергии PΣ . Это – средний поток энергии
___
через S: PΣ = Re PΣ . Последний член справа дает среднюю мощность источни-
____
ка: Pст = Re Pст . Таким образом, рассматриваемое равенство, которому удобно придать вид
___ |
____ |
|
PΣ = − |
ω∫(ε0ε′′Em Em +µ0µ′′Hm Hm )dV − Pст , |
(10) |
|
2 V |
|
является уравнением среднего баланса энергии при гармонических колебаниях.
____
Пусть источники в среднем отдают энергию поля Pст < 0 . Если проницаемости ε и µ вещественные (ε” = µ” = 0), то объемный интеграл в (10) исчезает. При этом в среднем вся мощность источников идет на излучение