Электродинамика сплошных сред
..pdfДругой пример черенковской волны – образование ударной волны (резкий скачок плотности, давления и температуры воздуха) при сверхзвуковом полете самолета. При этом возможно и проявление аномального эффекта в виде возбуждения источника – явление флаттер – вибрации самолета.
Излучение Вавилова-Черенкова нашло разнообразное применение в экспериментальной ядерной физике и физике элементарных частиц. На нем основано действие так называемых черенковских счетчиков, т.е. детекторов релятивистских заряженных частиц, излучение которых регистрируется с помощью фотоумножителей. Основное назначение черенковских счетчиков – разделение релятивистских частиц с одинаковыми импульсами, но различными скоростями. Пусть, например, пучок, состоящий из релятивистских протонов и мезонов, проходит через однородное поперечное магнитное поле. Направления траекторий прошедших частиц будут определяться только их импульсами, но не будут зависеть от их скоростей. С помощью диафрагм можно выделить протоны и мезоны с одинаковыми импульсами. Из-за различия масс скорости мезонов окажутся несколько больше скоростей протонов. Если полученный пучок направить в газ и подобрать показатель преломления газа, то мезоны будут давать излучение Вавилова-Черенкова, а протоны – нет. Таким образом, счетчик будет регистрировать только мезоны, но не будет регистрировать протоны.
1.9.4. Переменное электромагнитное поле в проводнике
Рассмотрим распределение электромагнитного поля по сечению проводника. Возьмем операцию ротора от уравнения Максвелла (1.148)
rot rot E = grad div E − 4E = µµ0 |
∂ rot H |
. |
||
∂t |
||||
Тогда уравнение будет выглядеть следующим образом: |
|
|||
∂E |
|
|
||
4E = σµµ0 |
|
. |
|
(1.205) |
∂t |
|
Аналогичное уравнение можно получить и для магнитного поля. Решим уравнение для случая переменного тока, текущего в омическом проводнике, заполняющем полупространство, z > 0. Направление тока вдоль
101
оси OX и его зависимость от времени заданы – jx(z, t) = j(z) exp ( i ωt); jy = jz = 0. В соответствии с законом Ома ищем электрическое поле, удовлетворяющее уравнению (1.205), в виде Ex(z, t) = E(z) exp ( i ωt);
Ey = Ez = 0. Получаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
∂2E(z) |
= i σµµ0ωE(z). |
|
(1.206) |
||||||||
|
∂z2 |
|
||||||||||
Общим решением будет E(z) = C 1 exp (−αz) + C 2 exp (αz), где |
||||||||||||
α2 = i ωσµµ0, |
|
|
|
|
i + 1 |
|
|
|||||
|
|
|
α = |
|
√ |
|
ωσµµ0. |
|
||||
|
|
|
2 |
|
||||||||
Тогда общее решение выглядит следующим образом: |
z!. |
|||||||||||
E(z) = C 1 exp − |
δs |
z! + C 2 exp |
δs |
|||||||||
|
|
i + 1 |
|
|
|
|
|
|
i + 1 |
|
||
Здесь введена величина |
|
|
|
|
|
|
|
|
|
|
|
|
|
δs = s |
|
2 |
|
, |
|
(1.207) |
|||||
|
|
|
||||||||||
|
ωσµµ0 |
|
которую называют глубиной проникновения поля в среду, или толщиной скин-слоя. При удалении на бесконечность, z → ∞, поле должно иметь конечное значение, поэтому константа C 2 = 0. Тогда выражение для поля будет выглядеть следующим образом:
Ex = E0 exp −δs ! exp "− i |
δs |
− ωt!#. |
(1.208) |
|
|
z |
z |
|
|
Формула (1.208) показывает, что напряженность электрического поля убывает экспоненциально вглубь проводника. Эффективное уменьшение напряженности (уменьшение в «e» раз) происходит на расстоянии δs от поверхности.
Зная распределение электрического поля в проводнике, можно найти магнитное поле H = (0, Hy, 0) следующим образом:
|
− i µµ0ωHy = rot yE = |
∂Ex |
|
|
i + 1 |
|
|
|
||||||||||||
|
|
|
|
= − |
|
|
|
|
Ex. |
|
|
|||||||||
|
|
∂z |
|
δs |
|
|
|
|||||||||||||
Отсюда |
y |
µµ0 |
ωδs |
−δs |
! |
|
"− |
|
δs |
− |
|
!# |
|
|||||||
|
|
|
|
|
||||||||||||||||
H |
|
= |
E0(1 |
− i ) |
exp |
|
z |
|
exp |
|
i |
|
z |
|
|
ωt |
. |
(1.209) |
||
|
|
|
|
|
|
|
|
|
|
|
|
102
Магнитное поле оказывается перпендикулярным электрическому. Оно убывает вглубь по тому же закону, что и электрическое. По абсолютной величине |Hy| ≈ |Ex|/ωµµ0δs.
Таким образом, электромагнитное поле и, соответственно, весь ток в проводнике оказываются локализованными в тонком поверхностном слое толщиной порядка δs. При ω → 0, т.е. при переходе к постоянному полю, δs → ∞ и скин-эффект исчезает. Чем выше электропроводность, тем меньше глубина проникновения магнитного поля в среду. В предельном случае (для сверхпроводника) внешнее поле в него не проникает. Все джоулево тепло выделяется в области скин-слоя.
Результат, полученный нами для простого случая, имеет совершенно общий характер. При любой геометрической конфигурации проводников поле в них оказывается локализованным в скин-слое. Скинэффект играет большую роль в технике переменных токов: он позволяет применять полые кабели или кабели, покрытые слоем металла с особенно высокой проводимостью. Также скин-эффект применяют для экранирования аппаратуры от переменного магнитного поля. На скин-эффекте основан метод фокусировки пучка электрических зарядов в металлической трубе: при приближении пучка к стенке трубы в ней наводятся токи, поля которых отталкивают пучок от стенки. Таким образом, пучок отталкивается трубой «со всех сторон» и устойчиво держится внутри трубы. Такая фокусировка позволяет транспортировать достаточно интенсивный пучок, например плазменный, по изогнутой трубе и, в частности, удерживать его в кольцевой трубе.
1.9.5. Прохождение электромагнитной волны сквозь диэлектрик
Рассмотрим распространение электромагнитной волны в среде. Диэлектрическая проницаемость среды будет зависеть от частоты поля, что определяется динамикой атомных электронов. Характерной особенностью этой динамики являются гармонические колебания электронов с разными частотами, что следует из вида атомных спектров. Поэтому
103
в качестве простейшей модели диэлектрика можно взять набор различных осцилляторов с собственными частотами ωj и плотностью соответствующих осцилляторов N j. Уравнение движения электрона для каждого осциллятора будет следующим:
∂2r j |
∂r j |
+ ω2j r = |
e¯ |
|
||
|
+ λj |
|
|
E exp [− i ωt] , |
(1.210) |
|
∂t2 |
∂t |
m e¯ |
где λj – коэффициенты, характеризующие потери энергии в диэлектрике. Решая уравнение (1.210) и суммируя по всем осцилляторам, найдем вектор поляризации
Xj |
|
Xj |
N j e¯2E/m e¯ |
|
|
|
|
ω2j |
− ω2 − i λjω |
|
|
||
P = |
e¯N jr j = |
|
|
|
. |
(1.211) |
Для случая малой плотности среды действующее электрическое поле равно среднему. Тогда выражение для диэлектрической проницаемости будет выглядеть таким образом:
Xj |
ω2j |
(ωp)2j |
|
(1.212) |
− ω2 − i λjω |
|
|||
ε(ω) = 1 + |
|
|
, |
где «плазменная» частота осцилляторов j определяется выражением
(ωp)2j = |
N j e¯2 |
(1.213) |
|
|
. |
||
|
|||
|
m e¯ |
|
Теперь найдем зависимость диэлектрической проницаемости от частоты поля для однородного и изотропного диэлектрика, помещенного в однородное магнитное поле с напряженностью H. Запишем уравнение движения атомного электрона, выбрав декартову систему координат с осью Z, параллельной вектору напряженности поля H.
∂2 x |
|
|
|
∂y |
|
|
|
e¯ |
|
|||||
|
|
+ ω02 x = |
−λ |
|
+ |
|
|
|
|
Ex exp [− i ωt] , |
|
|||
∂t2 |
∂t |
m e¯ |
|
|||||||||||
|
∂2y |
+ ω02y = λ |
∂x |
|
|
|
e¯ |
(1.214) |
||||||
|
|
|
+ |
|
|
Ey exp [− i ωt] , |
||||||||
|
∂t2 |
∂t |
m e¯ |
|||||||||||
|
|
|
|
∂2z |
+ ω02z = |
e¯ |
|
|||||||
|
|
|
|
|
|
Ez exp [− i ωt] . |
|
|||||||
|
|
|
|
∂t2 |
m e¯ |
|
104
Здесь ω0 – собственная частота атома-осциллятора. Потеря энергии вследствие действия силы Лоренца характеризуется величиной
λ = |
e¯H |
. |
(1.215) |
|
m e¯c
Решение системы (1.214) записывается в виде
x= aEx + i bEy,
y= aEy − i bEx,
z = |
e¯ Ez |
, |
|
|||||
|
|
|
|
|||||
m e¯ Ω |
|
|||||||
где |
|
|
||||||
|
|
|
|
|
||||
a = |
e¯ |
|
|
Ω |
|
, |
||
|
|
|
||||||
m e¯ |
Ω2 − ω2λ2 |
|||||||
b = |
e¯ |
|
ωλ |
, |
||||
m e¯ |
|
Ω2 − ω2λ2 |
Ω = ω20 − ω2.
Компоненты вектора электростатической индукции связаны с компонентами вектора поля тензорной зависимостью
Di = Ei + Pi = εi jE j,
где εi j – компоненты тензора диэлектрической проницаемости. Для их нахождения подставим в выражение для Pi значения x, y, z аналогично (1.211). В результате получим
|
εi j = |
|
|
i δ ε |
0 |
, |
||||
|
|
− |
ε |
i δ |
0 |
|
||||
|
0 |
0 |
ε0 |
|||||||
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ω2pΩ |
|
|||
|
ε = 1 + |
|
, |
|||||||
|
Ω2 − ω2λ2 |
|||||||||
|
δ = |
|
|
|
ω2pωλ |
|
, |
|
||
|
Ω2 − ω2λ2 |
|
||||||||
|
ε0 |
= 1 + |
ω2p |
. |
|
|
||||
|
Ω |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
105
Теперь рассмотрим, что произойдет с плоскостью поляризации плоской линейно поляризованной монохроматической волны, распространяющейся в диэлектрике вдоль слабого однородного магнитного поля. Волновое уравнение (1.170) использовать нельзя, поэтому для плоской монохроматической волны при µ = 1 используем выражение
|
|
|
|
|
|
∂2E 1 ∂2D |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
|
. |
|
||
|
|
|
|
|
|
∂z2 |
|
c2 |
∂t2 |
|
|||||||||
Для компонент Ex и Ey запишем |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
∂2Ex |
|
|
|
|
ω2 |
|
|||||||||||
|
|
|
|
|
= − |
|
|
|
(εEx + i δEy), |
|
|||||||||
|
|
∂z2 |
c2 |
|
|||||||||||||||
|
∂2Ey |
|
|
|
ω2 |
|
|||||||||||||
|
|
|
|
= − |
|
(− i δEx + εEx). |
|
||||||||||||
|
∂z2 |
c2 |
|
||||||||||||||||
Умножим второе уравнение на i |
|
и сложим с первым. Получим уравне- |
|||||||||||||||||
ние относительно величины E = Ex + i Ey, |
|
||||||||||||||||||
|
|
|
∂2E |
|
|
ω2 |
(1.216) |
||||||||||||
|
|
|
|
= |
|
− |
|
(ε + δ)E. |
|||||||||||
|
|
|
∂z2 |
|
c2 |
Теперь представим нашу линейно поляризованную волну в виде суперпозиции двух волн с круговой поляризацией:
Elin(z, t) = E0 ( exp [ i (kRz − ωt)] ± exp [− i (kLz − ωt)]). 2
Подставляя правую и левую составляющие в уравнение (1.216), получим выражение для волнового вектора
|
|
kR = |
ω √ |
|
|
|
|
kL = |
ω √ |
|
|
|||
|
|
ε + δ, |
|
ε − δ. |
||||||||||
|
|
c |
|
|
c |
|
||||||||
Для слабого магнитного поля, когда λ {ω, ω0}, величина ε ≈ ε0, |
||||||||||||||
√ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
а корень 1 ± δ/ε0 ≈ (1 ± δ/2ε0). Тогда |
|
|
|
|
||||||||||
|
|
kR,L |
≈ c |
√ε0 |
1 ± 2ε0 ! ≡ k0 ± 4k. |
|||||||||
|
|
|
|
|
ω |
|
|
δ |
|
|
|
|
Таким образом, на расстоянии l от входа в диэлектрик поле волны будет описываться выражением
Elin(z, t) ≈ E0 ( exp [ i (k0l − ωt)] ± exp [− i (k0l − ωt)])(1 + i 4kl). 2
106
Это выражение означает, что вектор Elin повернулся относительно вектора магнитного поля H на угол α,
α = arctg [4kl] = VH H · l, |
(1.217) |
где постоянная Верде
e¯ ω2pω2
VH = 2m e¯c2 pε0Ω2 .
Описанное явление поворота плоскости поляризации электромагнитного излучения, распространяющегося в среде вдоль магнитных силовых линий, называется эффектом Фарадея. Этот магнитооптический эффект был первым доказательством наличия прямой связи между магнетизмом и светом. Эффект применяется в различных устройствах радиосвязи, а также в астрофизических исследованиях. С помощью него изучаются космические магнитные поля.
107
2. Взаимодействие сред с электрическими и магнитными полями
2.1.МГД-течение в плоском слое электропроводящей жидкости
2.1.1. Уравнения магнитной гидродинамики
Магнитная гидродинамика (МГД) – это наука о движении электропроводных газов или жидкостей во взаимодействии с магнитным полем. При движении электропроводящей среды, находящейся в магнитном поле, в ней индуцируются электрические поля и токи, на которые действует магнитное поле и которые сами могут повлиять на магнитное поле. Таким образом, возникает сложная картина взаимодействия магнитных и гидродинамических явлений, которая должна рассматриваться на основе совместных уравнений гидродинамики и электродинамики.
Основными уравнениями магнитной гидродинамики являются: уравнение Навье-Стокса с учетом электромагнитных сил
ρw |
" ∂t |
+ (VO)V# = −OP + η4V + j × B + ρwg, |
(2.1) |
|
|
|
∂V |
|
|
уравнение неразрывности
∂ρw |
+ div (ρwV) = 0, |
(2.2) |
|
∂t |
|||
|
|
закон Ома
j = σ(E + V × B) |
(2.3) |
и уравнения Максвелла
108
rot E = − |
∂B |
, |
(2.4) |
|
|
∂t |
|||
rot B = µµ0j, |
(2.5) |
|||
div B = 0, |
|
|
(2.6) |
|
div E = 0, |
|
|
(2.7) |
|
div j = 0. |
|
|
(2.8) |
Токами смещения в данном случае пренебрегают вследствие их малости по сравнению с рассматриваемыми токами. Также считается, что плотность заряда равна нулю. Как было показано ранее, можно с помощью законов Ома и Ампера вывести уравнение переноса поля
∂B |
= |
1 |
4B + rot [V × B]. |
(2.9) |
∂t |
σµµ0 |
Также было указано, что уравнение записывают в безразмерном виде:
∂B |
− rot [V × B] = |
1 |
4B, |
(2.10) |
∂t |
Rem |
с использованием магнитного числа Рейнольдса Rem = V0µ0L0σ, которое характеризует отношение индуцированного движением магнитного поля к величине первоначального поля. Уравнение Навье-Стокса (2.1) тоже записывают в безразмерном виде с помощью гидродинамическо-
го числа Рейнольдса Reh |
= V0L0/ν следующим образом: |
|
|||||||
∂V |
|
|
1 |
|
|
|
|
||
|
|
+ (VO)V = −OP + |
|
4V + |
|
|
|||
|
∂t |
Reh |
|
|
|||||
|
|
+ |
Ha 2 |
[[V × B] × B] + |
N1 Ha |
[E × B]. |
(2.11) |
||
|
|
Reh |
Reh 2 |
Возможно использовать другой тип обезразмеривания, в котором в качестве масштаба скорости выбирается V0 = ν/L0. Тогда получается другая версия уравнения:
∂V |
+ (VO)V = −OP + 4V + N2[j × B]. |
(2.12) |
∂t |
Здесь параметры МГД-взаимодействия выражаются таким образом:
|
|
|
E |
L |
|
√ |
|
|
|
|
|
|
|
|
j0 B0L3 |
|
||
|
|
0 |
σρ |
w |
|
|
|
|
|
|||||||||
N |
|
= |
0 |
|
|
|
|
|
|
, |
N |
|
= |
0 |
. |
(2.13) |
||
1 |
|
p |
|
|
|
|
|
2 |
|
|||||||||
|
|
|
η3 |
|
|
|
|
|
ρwν2 |
109 |
Можно оценить порядок электромагнитных сил как fe σB20V0, вязкостных сил (учитываемых слагаемым ρwν4V) как fv = ρwνV0/L02 и сил инерции (учитываемых слагаемым ρw(VO)V) как fi = ρwV02/L0. Тогда порядок отношения электромагнитных сил к вязкостным силам описывается квадратом числа Гартмана,
|
|
f |
e |
|
σB2V0L2 |
|
σB2L2 |
|
||
Ha 2 |
= |
|
= |
0 0 |
= |
|
0 0 |
. |
(2.14) |
|
fv |
ρwνV0 |
|
||||||||
|
|
|
|
|
η |
|
Порядок отношения электромагнитных сил к силам инерции описывается числом Стюарта,
St = |
fe |
|
Ha 2 |
(2.15) |
|
|
= |
|
. |
||
fi |
|
||||
|
|
Reh |
|
Порядок отношения инерциальных сил к вязкостным силам описывается гидродинамическим числом Рейнольдса Reh . Следует упомянуть о числе Альфвена, характеризующем отношение магнитной энергии к кинетической,
Al = |
St |
|
|
B02 |
(2.16) |
|||
|
= |
|
|
|
|
. |
||
Rem |
µ |
ρ |
w |
V2 |
||||
|
|
0 |
|
0 |
|
|
Достаточно важным является магнитное число Прандтля, равное отношению гидродинамической и магнитной вязкостей,
Prm |
= |
Rem |
= |
ν |
. |
(2.17) |
|
Reh |
νm |
||||||
|
|
|
|
|
2.1.2. МГД-течения Куэтта и Гартмана
Рассмотрим течение электропроводной жидкости между двумя плоскопараллельными горизонтальными поверхностями, которые движутся в противоположные стороны со скоростями V0 под действием постоянного перепада давления P0 = −∂x P в присутствии поперечного внешнего магнитного поля B0 = (0, 0, B0) (рис. 2.1). Наша цель – полу-
чить профили скорости течения при Rem 1. |
|
Будем рассматривать стационарный установившийся |
процесс, |
в котором все производные по времени равны нулю, ∂t = 0. |
|
Рассматриваемое течение однородно в плоскости ZX, |
поэтому |
{∂xV; ∂yV} = 0. Если подставить это в уравнение неразрывности (2.2),
110