Расчет электромагнитного поля дипольных источников с использованием векторных потенциалов Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

Расчет электромагнитного поля дипольных источников с использованием векторных потенциалов

Кадников С.Н., д-р техн. наук, Веселова И.Е., ст. преп.

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

Ключевые слова: расчет поля, электромагнитное поле диполя, анизотропная среда, векторный потенциал, калибровка.

Calculation of electromagnetic field of dipole source using vector potential

Kadnikov S.N., Veselova I.E.

The design procedure of an electromagnetic field of an electric dipole in a homogeneous environment with anisotropy of electric and magnetic properties is considered. For vector potential expressions with elementary functions were found. It is shown that the result of calculation of a field with use of vector potential does not depend on a way of calibration.

Keywords: field calculation, electromagnetic field of dipole, anisotropic medium, vector potential, calibration.

Методика расчета электромагнитного поля диполей электрического типа в однородной одноосной анизотропной среде основана на использовании векторных потенциалов. Данная методика является естественным развитием схемы расчета электромагнитного поля в однородной изотропной среде [1]. Целесообразность использования векторных потенциалов обусловлена сложностью правых частей уравнений второго порядка относительно Еи Н , которые возникают при прямом решении уравнений Максвелла. Кроме того, их применение повышает точность расчета магнитных потоков и индуктивностей, поскольку не требует предварительного вычисления поля. С теоретической точки зрения расчет векторных потенциалов представляет интерес при исследовании эффекта Ааронова-Бома [2].

В изотропной среде фактически единственным эффективным (с расчетной точки зрения) способом введения связи между векторным и скалярным потенциалом является калибровка Лоренца. При этом возникают, как известно, три отдельных уравнения второго порядка относительно декартовых компонент векторного потенциала. В отличие от этого в анизотропной среде связь между векторным и скалярным потенциалами можно установить, вообще говоря, различными способами. В одноосной среде выбор формул связи (калибровки) является вполне очевидным и зависит от ориентации вектора магнитного диполя. Тем не менее, при наличии различных способов калибровки возникают как минимум два вопроса: 1) является ли очевидный с физической точки зрения выбор способа калибровки действительно оптимальным с вычислительной точки зрения; 2) зависит ли электромагнитное поле, т.е. функции Е и Н , от способа калибровки.

Рассмотрим анизотропную среду с анизотропией электрических и магнитных связей (квазистационарное синусоидальное электромагнитное поле). Исходная система уравнений имеет следующий вид:

г&Н = уЕ + р8; (1)

ОЕ = -] юД аН, (2)

где р - вектор плотности электрического тока (постоянный вектор); 8 = f (х2у2 z)- дельта-

функция.

В выбранной системе координат

y o o 1 o o

Y = o y o . 11 = o 1 o

o o Y z o o 1 z

где y*Yz.1*1z . Определим,

(2), индукцию

согласно

В = ]1аН формулой В = г&Л , где А - векторный магнитный потенциал. Тогда, согласно (1), Е = -]юА-Уф. Подставляя эти выражения, а

также H = 1 JrotA в (1), получим следующую систему уравнений:

д

дф

А1 Ax - j ®1o1z YAx = diviA + 1o1z Yt- -

-1o1zPx S.

dx

d

dx

дф

А1 Ay - j®1o1zYAy = — diviA + 1o1z Y— -

-1o1zPy S.

dy

d

dy

дф

А1 Az - j® 1o1YzAz =^T diviA + 1o1Yz -r- -

dz

dz

(З)

(4)

(5)

1o1zpzS.

д2 д2

д

где Аіі - + -Т^ + Т^Г;

дх2 ду2

дА дАу

х у + ^„

аЦ. А - + —

и дх ду

ид2 2 дА,

X., -

и,

Рх . Ру . Р, -

д, г и

компоненты постоянной вектора плотности тока в начале координат.

Рассмотрим сначала вариант вертикального диполя рх - ру - 0, р, ф 0 . Тогда, по-

1

лагая ф- —

.0., У

а'м.. а,

I-1

(6)

получим вместо (4), (5) два однородных уравнения (во всем пространстве), решения которых Ах = Ау = 0. Уравнение (5) с учетом (6) примет вид

д2 А, д2 А, .

1А +-дА + ^

ду 2

дА

дг2

А, .0.,Р,8 ,

(7)

У 2

где к2 = уюц0цу ; К=— .

У

После замены г = это уравнение

приводится к известному уравнению Гельмгольца, решение которого после обратной замены г., = дается формулой

А .0.Рэ Є

4п

~К~

(8)

где рэ - момент электрического диполя с размерностью Ам; = ^х2 + у2 + Х-1г2 .

Отметим, что приведенный способ калибровки не является единственным. Если положить

ф = -

1

.0. У;

-<а/у„ а,

I-1

то уравнение (5) примет вид

д2 А,

(9)

(10)

дх2 ду2 г д,2

После подстановки решения этого уравнения в (3), (4) и при использовании нормировки (9) получится система, состоящая из двух уравнений относительно Ах, Ау, в правых частях которой появятся производные решения (10) (известные функции). Однако такой способ решения, очевидно, лишен практического смысла, поскольку электромагнитное поле будет выражаться через три компоненты векторного потенциала, в то время как в первом способе достаточно одной компоненты А2, с помощью которой легко определяется магнитное поле, имеющее только азимутальную составляющую, а через него с использованием (1) - электрическое.

Рассмотрим теперь вариант диполя с

моментом рх ф 0 , Ру - р, - 0 (горизонтальный

диполь). После использования нормировки (6) уравнения (3) и (4) приведем к виду

А. Ах - к2х Ах --ц^Рх8. (11)

(12)

где кх - у®UoUzУ .

Последнее уравнение (5) переходит в следующее:

д2 А, д2 А, .

—2z + —^ + X, дх 2

ду 2

!У, ■

д,

- к\А, - (1 -X)

д 2 Ах

д,дх ’

(13)

где Х = ; к2 = у юцрцу 7 ■

Ц 7 У

Чтобы решить это уравнение, используем спектральный метод [1]. Прямое и обратное преобразование Фурье векторного потенциала представим в виде

і (хх+уу + ,,

)хауа,, (14)

+0)

А(х.у,z)--^-т ГГГА(х,у,Z)eіXx+yУ+zZ'}dXdydZ. (15) (2п) ^

Используя формулу (15), подставим составляющую Ах векторного потенциала в (7) и получим

А -__________u0uzpэ______

Ах - 2 - 2 - 2 і 2 ’

х + у + Ху •, + кх

где рэ - момент электрического диполя с размерностью Ам.

Тогда из уравнения (12) следует

________И^э (1 - X) х,____________

-

(16)

Подставляя это выражение в (14), получаем решение (12) в виде

Н-рН-7РЭ(1 — Х)

А, -

(2п)3

х,е

ахауа,.

(17)

Интеграл по 7 в этой формуле вычислим с помощью теории вычетов [3]:

ге

, (х + у2 + Ху,2 +

Х.(1 -Х)

в/дп,

к2 2 ^ у2 + Xuz + кх)

e-“xlz e~az\z

а, -

-2 -2 -2 -2 х + 'у х + у

( <.2

где ах -

2

х + у ,2

—— + к X,,

( ,-,2

2

х + ,2

—— + к

X.

к - У^.0! У ■

Далее нужно вычислить двойной интеграл вида

+ад

.. хеу(хх+уу )е-ах,| хеу(хх+уу )е-аАА

Г1^--------------------------С+у*) аха,у■ (18)

После замены декартовых координат цилиндрическими г, а и использования формулы [4]

ад

Г

^jrr cos(a-a) =

= 1 j

nejn(a-a

Jn (rr),

где J (гг) - функция Бесселя порядка п, интеграл (17) дает следующее окончательное решение для Л2:

A = 1o1 Рэ cos а z

z 4п r

kj

k, К, r2

К1r 2

I

(19)

где ц - составляющая тензора магнитной проницаемости по оси х (или у); рэ - момент электрического диполя; Ху=у z/y, Хц=ц2/ц ;

г = ^х2 + у2 ; cos a = x/r .

При = 1, т.е. для однородной в магнитном отношении среды, формула (19) совпадает с выражением, полученным в [5], где оно выведено методом разделения переменных Фурье, причем исходное представление для векторного потенциала введено, по сути, эвристически. Кроме того, при таком подходе возникает необходимость дополнительной процедуры по вычислению момента диполя, что отмечается как недостаток в [6]. Использованный метод, основанный на преобразованиях Фурье (14), (15), не требует дополнительных процедур и сразу дает готовое решение. Кроме того, такой подход позволяет, что более существенно, сравнительно просто проанализировать различные способы нормировки. В частности, в данном случае (горизонтального диполя) возможен вариант нормировки следующего типа:

ф = --

1

1o1Y;

-div,. A.

I-1

(20)

Уравнение (5) после этого приобретает вид

А. А, - к,2 А, - 0. (21)

Решение этого уравнения, очевидно, равно нулю. Из (4), (5) с использованием (20) возникает следующая система:

К

d2 Ax d2 Ax „

—x+^т+К1

dx2 ду

= -1o1zPxS’

d2 Ax dz 2

5% dx 2

-К-

d2 Ay

dy 2

6%

1 dz2

d2 Av 2

______j - k2 A =

dxdy x x (22)

d2Ax і 2 л л ----x - ky Ay = o.

dy dx

(23)

Решение этой системы дает две ненулевые компоненты Ах, Ау. Обозначим ненулевые компоненты потенциалов в задаче о горизонтальном диполе: для первого варианта решения - А., А., ф1, для второго - А2, А^, ф2 .

Выпишем выражения для составляющих векторов поля в обоих вариантах:

1-1 ■ л1 дф1.

Ex = - j и Ax —"

dx

E1 = ^ф1 ■

Ey dy ’

i-1 ■ л 1 дф1.

Ez = -j и Az д ■

dz

B1 =dAz; B1 _a_■ B1 = dAx

Bx dy ’ By dz dx ’ Bz dy '

1-2 ■ Л 2 дф2

e2 =-j и Ax -^r- ■ dx

2

31 5a1 5a1 _ _1

5A1

E2y =-j и A' --^- ■

2 ■

dy

Ez2 =_

ф1 =■

дф2 ■ dz ■ dAj dz 1

R* ^у ■ R2 dAx

Bx = _— ■ By =~dz

bZ = dAy _ dAx z dx dy

-div1 A1, ф2 =----------------------------------------1— div11 A2.

НоНг У НоНУг

По форме записи эти выражения существенно различаются. Тем не менее они представляют одно и то же решение, т.е.

Е1 = Е2 = Е, В1 = В2 = В. Это можно показать, используя данные выражения как представления для произвольных векторов. Если, например, взять ротор от вектора Е1, то, согласно формулам, для вектора В1 получим уравнение (2). Если же взять ротор от вектора Н1 = ц0Д-1В1, то при условии, что уравнения (11), (13) выполняются, получим формулы для вектора у Е1, т.е. придем к системе уравнений Максвелла (однородность этих уравнений в данном случае принципиального значения не имеет). То же самое можно проделать и с формулами для векторов Е2, В2, т.е. снова получить уравнения Максвелла. Поскольку, как известно, решения уравнений Максвелла единственны, равносильны будут и два рассмотренные способа нормировки.

В заключение отметим, что полученные выражения для векторных потенциалов являются в некоторой степени универсальными. В частности, полагая н = н7 (однородная магнитная среда) и заменяя действительный тензор у на комплексный у + у'ю80ё, можно получить выражения для векторных потенциалов в анизотропной диэлектрической среде с потерями [6].

Поскольку полученное выражение для векторного потенциала является довольно простым, с его помощью легко можно найти выражения для векторов электромагнитного поля диполя, которые являются, как известно, фундаментальными решениями уравнений Максвелла, имеющими практическое значение [7].

П=—і

З

Список литературы

1. Марков Г.Т., Чаплин А.Ф. Возбуждение электромагнитных волн. - М.: Радио и связь, 1983.

2. Чирков А.Г., Агеев А.Н. О природе эффекта Ааро-нова-Бома // Журнал технической физики. - 2001. - Т.1. -Вып. 2. - С. 16-22.

3. Лаврентьев М.А., Шабат Б.В. Методы теории функций комплексного переменного. - М.: Наука, 1965.

4. Лебедев Н.Н. Специальные функции и их приложения. - М.: ГИФМЛ, 1963.

5. Кауфман А.А., Качанский А.М. Индукционный метод изучения поперечного сопротивления в скважинах. -Новосибирск: Наука, 1972.

6. Савченко А.О., Савченко О.Я. Электромагнитное поле диполя в анизотропной среде // Журнал технической физики. - 2005. - Т.75. - Вып.10. - С. 118-121.

7. Тозони О.В., Маергойз И.Д. Расчет переменных электромагнитных полей. - Киев: Техника, 1975.

Кадников Сергей Николаевич,

Ивановский государственный энергетический университет,

доктор технических наук, профессор кафедры теоретической электротехники и электротехнологии, e-mail: kadnikovsn@mail.ru

Веселова Ирина Евгеньевна,

Ивановский государственный энергетический университет имени Ленина, ассистент кафедры высшей математики, e-mail: iveselova@math.ispu.ru