Электромагнитное поле, излучаемое движущимся диполем

Язык труда и переводы:
УДК:
621.3.095.3
Дата публикации:
23 декабря 2022, 12:20
Категория:
Перспективные направления исследования необратимых физических процессов
Авторы
Аннотация:
Рассмотрена задача об излучении электромагнитной волны в вакууме движущимся диполем. Проведены вычисления электромагнитных полей при двух видах движения заряда — равноускоренное движение (закон движения отрицательного заряда Z(t)=β∙t^4), колебания по закону синуса (закон движения отрицательного заряда Z(t)=α∙sin(t)). Исследовано распределение электромагнитных полей в пространстве, полученное при вычислении. Проанализированы зависимости полей движущегося диполя, а также зависимости от расстояния до точки наблюдения, а также поляризация электромагнитных волн, излучаемых зарядом.
Ключевые слова:
электромагнитное поле диполя, движение диполя, дипольный момент, напряженность магнитного поля, скорость и ускорение движения заряда, поляризация электромагнитных волн
Основной текст труда

Электромагнитное поле, излучаемое движущимся диполем

Рассмотрим задачу об излучении электромагнитной волны в вакууме движущимся диполем в постановке, идея которого принадлежит Р. Фейнману [1]. Выберем систему координат Oxyz , в которой положительный заряд Q движется по закону {\overrightarrow {r}}^{+}\left(x,y,z,t\right) , а отрицательный заряд -Q  по закону {\overrightarrow {r}}^{-}\left(x,y,z,t\right) , постоянно находясь от положительного заряда на расстоянии l . Такая пара зарядов образует электрический диполь с “собственным” дипольным моментом, определяемым вектором {\overrightarrow {p}}_{0}=Q{\overrightarrow {l}} .

Через векторный потенциал A  и скалярный потенциал \varphi выводится формула вектора магнитной индукции одиночного заряда {\overrightarrow {B}}={\frac {\left\lbrack {\overrightarrow {\dot {p}}}\left(t-{\frac {r}{c}}\right),{\overrightarrow {n}}\right\rbrack }{cr^{2}}}+{\frac {\left\lbrack {\overrightarrow {\ddot {p}}}\left(t-{\frac {r}{c}}\right),{\overrightarrow {n}}\right\rbrack }{c^{2}r}},(1)

где {\overrightarrow {r}}  задает вектор от заряда до точки наблюдения P\left(x,y,z\right) , соответственно, {\overrightarrow {n}}={\frac {\overrightarrow {r}}{r}} — единичный вектор, задающий направление.

Производные дипольного момента введены в соответствии с методом Р. Фейнмана. Помещая пару фиктивных зарядов Q  и -Q в начало координат, равных по величине рассматриваемому заряду, рассматривается диполь, составленный из рассматриваемого заряда одного фиктивного с противоположным знаком. Компоненты поля в проекции на оси координат раскрываются в следующем виде: 

{B_{x}={\frac {z}{cr^{3}}}\left\lbrack {\dot {p}}_{y}+{\frac {r}{c}}{\ddot {p}}_{y}\right\rbrack -{\frac {y}{cr^{3}}}\left\lbrack {\dot {p}}_{z}+{\frac {r}{c}}{\ddot {p}}_{z}\right\rbrack },

{B_{y}={\frac {x}{cr^{3}}}\left\lbrack {\dot {p}}_{z}+{\frac {r}{c}}{\ddot {p}}_{z}\right\rbrack -{\frac {z}{cr^{3}}}\left\lbrack {\dot {p}}_{x}+{\frac {r}{c}}{\ddot {p}}_{x}\right\rbrack },(2) {B_{z}={\frac {y}{cr^{3}}}\left\lbrack {\dot {p}}_{x}+{\frac {r}{c}}{\ddot {p}}_{x}\right\rbrack -{\frac {x}{cr^{3}}}\left\lbrack {\dot {p}}_{y}+{\frac {r}{c}}{\ddot {p}}_{y}\right\rbrack .} Производные дипольных моментов пары зарядов, образующих диполь, запаздывают на разное время. Дипольный момент положительного заряда зависит от времени {\overrightarrow {p}}={\overrightarrow {p}}\left(t-{\frac {r^{+}}{c}}\right) , в то время как дипольный момент отрицательного заряда — {\overrightarrow {p}}={\overrightarrow {p}}\left(t-{\frac {r^{-}}{c}}\right) . Магнитное поле, создаваемое диполем, как единой системой, находится как сумма магнитных полей зарядов, образующих его {\overrightarrow {B}}_{\text{dipol}}={\overrightarrow {B^{+}}}+{\overrightarrow {B^{-}}} . Для этого требуется преобразовать поле, например, отрицательного заряда таким образом, чтобы оно запаздывало на то же время, что и положительный. Используем разложение дипольного момента в ряд Тейлора до первой производной: 

p\left(t-{\frac {r^{-}}{c}}\right)=p\left(t-{\frac {r^{+}}{c}}\right)-{\dot {p}}\left(t-{\frac {r^{+}}{c}}\right){\frac {r^{-}-r^{+}}{c}},(3)

тогда x — компонента вектора магнитной индукции отрицательного заряда и диполя: 

B_{x}^{-}={\frac {1}{cR^{3}}}\left\lbrack z{\dot {p}}_{y}-y{\dot {p}}_{z}\right\rbrack +{\frac {r}{c^{2}R^{3}}}\left\lbrack z{\ddot {p}}_{y}-y{\ddot {p}}_{z}\right\rbrack -{\frac {R-r}{c^{3}r^{2}}}\left\lbrack z{\stackrel {...}{p}}_{y}-y{\stackrel {...}{p}}_{z}\right\rbrack .(4)

С учетом того, что дипольные моменты зарядов имеют разные знаки, находим магнитное поле диполя: 

B_{x}={\frac {R^{3}-r^{3}}{cr^{2}R^{3}}}{\frac {\left\lbrack z{\dot {p}}_{y}-y{\dot {p}}_{z}\right\rbrack }{r}}+{\frac {R^{3}-r^{3}}{c^{2}rR^{3}}}{\frac {\left\lbrack z{\ddot {p}}_{y}-y{\ddot {p}}_{z}\right\rbrack }{r}}+{\frac {\left(R-r\right)r}{c^{3}R^{2}}}{\frac {\left\lbrack z{\stackrel {...}{p}}_{y}-y{\stackrel {...}{p}}_{z}\right\rbrack }{r}}.(5)

Вектор магнитной индукции диполя в момент времени t-{\frac {r^{+}}{c}}

{\overrightarrow {B}}_{dipol}={\frac {R^{3}-r^{3}}{cr^{2}R^{3}}}\left\lbrack {\overrightarrow {\dot {p}}},{\overrightarrow {n}}\right\rbrack +{\frac {R^{3}-r^{3}}{c^{2}rR^{3}}}\left\lbrack {\overrightarrow {\ddot {p}}},{\overrightarrow {n}}\right\rbrack +{\frac {\left(R-r\right)r}{c^{3}R^{2}}}\left\lbrack {\overrightarrow {\stackrel {...}{p}}},{\overrightarrow {n}}\right\rbrack .(6)

Расстояние от начала координат до положительного заряда можно определить следующими способами:

  • R=r+l\cos {\theta \ } — приближенное выражение, применимо для нахождения поля в сильно удаленной от диполя точке.
  • R={\sqrt {r^{2}+l^{2}-2rl\cos \theta }}   — точное аналитическое значение.

Движение диполя вдоль оси

Исследуем подробно электромагнитное поле, излучаемое диполем, движущимся вдоль оси Oz  со скоростью v , как показано на рис. 1.

Рис. 1. Диполь, движущийся вдоль оси Oz [2]

Получим для x  составляющей вектора индукции B_{x}^{+} магнитного поля, излучаемого положительным электрическим зарядом: B_{x}^{+}=-{\frac {\mu _{0}}{4\pi }}{\frac {y}{r^{2}}}\lbrack {\frac {1}{r}}{\dot {p}}_{z}^{+}\left(t-{\frac {r}{v}}\right)+{\frac {1}{v}}{\ddot {p}}_{z}^{+}\left(t-{\frac {r}{v}}\right)\rbrack ,(7a) где в соответствии с методом Р. Феймана введено обозначение p_{z}^{+}=Qz .

Аналогичным образом определим x  составляющую вектора индукции {\overrightarrow {B}}^{-}  магнитного поля, излучаемого отрицательным электрическим зарядом: 

B_{x}^{-}=-{\frac {\mu _{0}}{4\pi }}{\frac {y}{\left(r+l\cos {\theta \ }\right)^{2}}}\lbrack {\frac {1}{r+l\cos {\theta \ }}}{\dot {p}}_{z}^{-}\left(t-{\frac {r+l\cos {\theta \ }}{v}}\right)+{\frac {1}{v}}{\ddot {p}}_{z}^{+}\left(t-{\frac {r+l\cos {\theta \ }}{v}}\right)\rbrack ,(7b)

где p_{z}^{-}=-Qz

Разложим в ряд Тейлора входящие в выражение (7b) функции {\frac {1}{\left(r+l\cos {\theta \ }\right)^{2}}} , {\frac {1}{\left(r+l\cos {\theta \ }\right)^{3}}} , {\dot {p}}_{z}^{-}\left(t-{\frac {r+l\cos {\theta \ }}{v}}\right)  и {\ddot {p}}_{z}^{-}\left(t-{\frac {r+l\cos {\theta \ }}{v}}\right) . Учитывая, что излученные диполем электромагнитные волны определяются на большом расстоянии от зарядов, т. е. при r\gg \ l\cos {\theta \ } , получим: 

{\frac {1}{\left(r+l\cos {\theta \ }\right)^{2}}}\approx {\frac {1}{r^{3}}}(1-2lcos\theta ),(8a)

{\frac {1}{\left(r+l\cos {\theta \ }\right)^{3}}}\approx {\frac {1}{r^{3}}}(1-3lcos\theta ),(8b)

{\dot {p}}_{z}^{-}\left(t-{\frac {r+l\cos {\theta \ }}{v}}\right)\approx {\dot {p}}_{z}^{-}\left(t-{\frac {r}{v}}\right)-{\ddot {p}}_{z}^{-}\left(t-{\frac {r}{v}}\right){\frac {l\cos {\theta \ }}{v}},(8c)

{\ddot {p}}_{z}^{+}\left(t-{\frac {r+l\cos {\theta \ }}{v}}\right)\approx {\ddot {p}}_{z}^{+}\left(t-{\frac {r}{v}}\right)-{\stackrel {...}{p}}_{z}^{-}\left(t-{\frac {r}{v}}\right){\frac {l\cos {\theta \ }}{v}},(8d)

Учитывая выражения (8a), (8b), (8c) и (8d), получим выражения для вектора индукции {\overrightarrow {B}} магнитного поля, излучаемого движущимся диполем: 

{\overrightarrow {B}}={\overrightarrow {B}}^{+}+{\overrightarrow {B}}^{-}={\frac {\mu _{0}l\cos {\theta }}{4\pi r^{2}v^{2}}}\lbrack {\stackrel {...}{p}}_{z}^{+}\left(t-{\frac {r}{v}}\right),{\overrightarrow {r}}\rbrack ,(9)

где учтено, что p_{z}^{-}=-p_{z}^{+} .

С помощью формулы (9) определим среднюю интенсивность W_{\Sigma }  электромагнитного поля, излучаемого движущимся диполем с относительной диэлектрической и магнитной проницаемостями \varepsilon  и \mu

W_{\Sigma }=\oint _{\Sigma }^{}{\left({\overrightarrow {S}},d{\overrightarrow {\Sigma }}\right)={\sqrt {\frac {\mu }{\varepsilon }}}}{\frac {{\stackrel {...}{p}}^{2}}{30\varepsilon _{0}\pi ^{2}v^{5}}}\int _{0}^{\pi }{\int _{0}^{2\pi }{\sin ^{3}\theta \cos ^{3}\theta d\theta d\varphi =}}{\sqrt {\frac {\mu }{\varepsilon }}}{\frac {{\stackrel {...}{p}}^{2}}{30\varepsilon _{0}\pi v^{3}}},(10)

где {\overrightarrow {S}} — вектор Пойтинга, d{\overrightarrow {\Sigma }}=\ r^{2}\sin \theta {\text{dθdφ}}{\frac {\overrightarrow {r}}{r}} — ориентированный элемент сферической поверхности радиуса \ r , определенный полярным углом \theta и азимутальным углом \varphi ; \int \limits _{0}^{\pi }{\sin ^{3}\theta \cos ^{3}\theta \ d\theta }=\ -\int \limits _{0}^{\pi }{\left(1-\ \cos ^{2}\theta \right)\cos ^{2}\theta \ d}\cos \theta ={\frac {4}{15}} .

При расчетах использовано точное расстояние от зарядов до точки наблюдения P\left(x,y,z\right) .

  • От отрицательного заряда: r\left(t\right)={\sqrt {x^{2}+y^{2}+\left(Z\left(t\right)+l-z\right)^{2}}} ,
  • От положительного заряда: R\left(t\right)={\sqrt {x^{2}+y^{2}+\left(Z\left(t\right)-z\right)^{2}}} ,

где Z\left(t\right) — закон движения отрицательного заряда, известный изначально.

Математическое моделирование для равноускоренного движения 

Результаты математического моделирования для равноускоренного движения (закон движения отрицательного заряда  Z(t)=\beta \cdot t^{4} ) показаны на рис. 2—4.

 
Рис. 2. Распределение поля в плоскости Оху
Рис. 3. Графики поля в зависимости от расстояния между зарядами: синий цвет — величина магнитного поля отрицательного заряда; красный цвет — величина магнитного поля положительного заряда положительного; зеленый цвет — суммарное поле диполя
Рис. 4. Графики компонент магнитного поля:  1 — компонента со второй производной дипольного момента; 2 — компонента с третьей производной

Математическое моделирование для колебаний по закону синуса 

Результаты математического моделирования для колебаний по закону синуса [3] (закон движения отрицательного заряда Z(t)=\alpha \cdot \sin(t) ) показаны на рис. 5–7.

Рис. 5. Распределение поля в плоскости Оху
Рис. 6. Графики поля в зависимости от расстояния между зарядами вблизи начала координат: синий цвет — величина магнитного поля отрицательного заряда; красный цвет — величина магнитного поля положительного заряда положительного; зеленый цвет — суммарное поле диполя
Рис. 7. Графики компонент магнитного поля: 1 — компонента со второй производной дипольного момента; 2 — компонента с третьей производной

Оценка графических результатов

На рис. 4. компонента со второй производной и компонента с третьей производной при сопоставлении формы графической зависимости имеют практически одинаковый характер.

На рис. 3 и 6 представлена графическая зависимость поля (компонента вектора магнитной индукции) от расстояния между зарядами, на которых представлено не только суммарное влияние, но и влияние, оказываемое каждым зарядом. При анализе разных типов движения амплитуда поля имеет различие в один порядок. Таким образом, расположение заряда относительно точки наблюдения и расстояние между зарядом оказывает существенное влияние на распределение поля в пространстве.

На рис. 7 изображены компоненты магнитного поля, которые при расчете используют составляющие второй и третьей производной. Характер движения рассматриваемого диполя находит отражение и в форме получаемых графиков для компонент магнитного поля. Данный график третьей производной ассиметричен, а в случае второй производной график имеет симметрию, отражающую периодичный характер движения исходного диполя. Вклад составляющей, рассчитанной по третьей производной при данных условиях, оказывается на девять порядков ниже, чем вторых производных, что свидетельствует о более быстром затухании электромагнитного поля в зависимости от расстояния, однако величина такой составляющей остается существенной для использования в последующих расчетах.

Анализ численных результатов

В случае равноускоренного движения магнитная индукция: B_{p}=10^{-12}{\text{Тл}} , а в случае колебаний по закону синуса: B_{k}=10^{-13} Тл.

Запишем формулу для магнитной индукции: {\overrightarrow {B}}=\mu \mu _{0}{\overrightarrow {H}}. В дальнейшем значения относительной диэлектрической ( \varepsilon ) и магнитной ( \mu ) проницаемостей примем равные единице: \varepsilon =1;\mu =1.(11) Напряженность магнитного поля: H={\frac {B}{\mu _{0}}}.(12) Равенство, связывающее амплитуды колебаний электрической и магнитной компонент [2, 4]: E={\sqrt {\frac {\mu \mu _{0}}{\varepsilon \varepsilon _{0}}}}H.(13)

Учитывая соотношения (11), (12) и (13), найдем напряженность электрического поля для каждого из видов движений и \varepsilon _{0}=8,85\cdot 10^{-12}{\frac {\text{Ф}}{\text{м}}} , \mu _{0}=4\pi \cdot 10^{-7}{\frac {\text{Гн}}{\text{м}}} (соответственно электрическая и магнитная постоянные): E_{p}=3\cdot 10^{-4}{\frac {\text{В}}{\text{м}}};E_{k}=3\cdot 10^{-5}{\frac {\text{В}}{\text{м}}},(14)

вектор Пойнтинга: {\overrightarrow {S}}=\left\lbrack {\overrightarrow {E}},{\overrightarrow {H}}\right\rbrack ,(15) Учитывая соотношения (12), (14), (15) и значения магнитной индукции для двух типов движения: S_{p}=2,4\cdot 10^{-10}{\frac {\text{Вт}}{{\text{м}}^{2}}};S_{k}=2,4\cdot 10^{-12}{\frac {\text{Вт}}{{\text{м}}^{2}}},

соответственно, интенсивность поля: I_{p}=1,2\cdot 10^{-10}{\frac {\text{Вт}}{{\text{м}}^{2}}};I_{k}=1,2\cdot 10^{-12}{\frac {\text{Вт}}{{\text{м}}^{2}}}. При исследовании характера распределения электромагнитного поля с учетом расположения оказалось, что вблизи оказывает влияние компонента третьей производной, а по мере удаления имеет место классическое излучение, то есть на больших расстояниях данной частью можно пренебречь.

Заключение

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

Литература
  1. Фейнман Р., Лейтон Р., Сендс М. Феймановские лекции по физике. Электричество и магнетизм. Т. 3. Москва, АСТ, 2019, 304 с.
  2. Литвинов О.С., Горелик В.С. Электромагнитные волны и оптика: учеб. пособие для вузов. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2006, 446 с.
  3. Хармут Х.Ф. Несинусоидальные волны в радиолокации и связи. Москва, Радио и связь, 1985, 375 с.
  4. Горелик Г.С. Колебания и волны. Москва; Ленинград, Гостехтеоретиздат, 1950, 572 с.
Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.