УДК 517.95
DOI: 10.18101/2304-5728-2019-3-17-31
ЗАДАЧА ДИРИХЛЕ ДЛЯ УРАВНЕНИЙ ГЕЛЬМГОЛЬЦА В ГИРОТРОПНОЙ ЭЛЛИПТИЧЕСКОЙ ОБЛАСТИ ПРИ ПРОДОЛЬНОМ НАМАГНИЧИВАНИИ
© Ширапов Дашадондок Шагдарович
доктор физико-математических наук, профессор, Бурятский государственный университет имени Доржи Банзарова Россия, 670000, г. Улан-Удэ, ул. Смолина, 24а E-mail: [email protected]
© Итигилов Гарма Борисович
кандидат технических наук, доцент,
Восточно-Сибирский государственный университет технологий и управления Россия, 670013, г. Улан-Удэ, ул. Ключевская, 40В E-mail: [email protected]
© Юмов Игорь Бимбаевич
кандидат физико-математических наук, доцент, Бурятский государственный университет имени Доржи Банзарова Россия, 670000, г. Улан-Удэ, ул. Смолина, 24а E-mail: [email protected]
© Анахин Владимир Дмитриевич
доктор технических наук, профессор,
Бурятский государственный университет имени Доржи Банзарова Россия, 670000, г. Улан-Удэ, ул. Смолина, 24а E-mail: [email protected]
© Дамбаев Жаргал Гомбоевич
доктор технических наук, профессор,
Бурятский государственный университет имени Доржи Банзарова Россия, 670000, г. Улан-Удэ, ул. Смолина, 24а E-mail: [email protected]
Поставлена и решена задача Дирихле для уравнений Гельмгольца электромагнитных волн, распространяющихся в эллиптическом цилиндре, заполненном продольно намагниченным ферритом, который описывается тензором второго ранга. Предполагается, что цилиндр имеет бесконечно проводящую стенку. Для решения краевой задачи уравнений Гельмгольца для продольных компонент электромагнитных волн применяется метод укорочения исходного дифференциального уравнения и метод разделения переменных. Решение указанной краевой задачи в эллиптических координатах связано с использованием четных и нечетных обыкновенных и модифицированных функций Матье 1-го рода. Используя полученные результаты, определены все компоненты электромагнитных волн для четных и нечетных решений. Применив условие Дирихле к компонентам электромагнитных волн и решив систему линейных
однородных алгебраических уравнений, получены дисперсионные уравнения электромагнитных волн, которые имеют важное практическое значение и позволяют проводить исследования распространения гибридных волн в данной области.
Ключевые слова: эллиптический цилиндр; феррит; задача Дирихле; уравнение Гельмгольца; электромагнитная волна; продольное намагничивание; гиро-тропная область; поперечные компоненты электромагнитного поля; функции Матье; дисперсионное уравнение.
Для цитирования:
Ширапов Д. Ш., Итигилов Г. Б., Юмов И. Б., Анахин В. Д., Дамбаев Ж. Г. Задача Дирихле для уравнений Гельмгольца в гиротропной эллиптической области при продольном намагничивании // Вестник Бурятского государственного университета. Математика, информатика. 2019. № 3. С. 17-31.
Введение
Вопросам распространения электромагнитных волн в магнитоактивных средах начали уделять пристальное внимание во второй половине XX в. в связи с дальнейшим развитием радиотехники и совершенствованием технологии получения ферритов, когда стало возможным изготовление невзаимных устройств в области сверхвысоких частот. Большой вклад в развитие теории распространения электромагнитных волн в ферритовых средах внесли такие видные ученые, как Г. Сул, Л. Уокер [1], Б. Лакс, К. Баттон [2], в чьих работах определение векторов электромагнитного поля осуществляется прямым решением системы уравнений Максвелла в строгом виде. В этих и других работах, как правило, не рассматриваются эллиптические волноводы с ферритовым заполнением или рассматривается изотропный случай [3]. Поэтому является актуальным анализ распространения электромагнитных волн в эллиптических волноводах с ферритовым заполнением.
Рассматривается продольно намагниченный эллиптический цилиндр с бесконечно проводящей стенкой. Область цилиндра заполнена ферритом, диэлектрическая проницаемость е которого изотропна, а тензор магнитной проницаемости определяется выражением
( т & о ^ т о о 0 Ц\\
где ] — мнимое число, ц = №0 _ т0—2—^Г, т0 = 4^-10 Гн/м — магнитная постоянная, С0 = т0УН0 — частота ферромагнитного резонанса, У = 1,76 -1011 Кл/кг — гиромагнитное отношение, Н0 — напряженность
постоянного магнитного поля, wm = /u0YM0, M0 — намагниченность на-
сыщения феррита, k = m0
—2-- , w — циклическая частота, m » m0-
w -wn
Для исследования распространения электромагнитных волн (ЭМВ) в этой области необходимо знать в эллиптической системе координат (£ ф, г), как поперечные компоненты электрического (Е%, Еф) и магнитного (И, Нф) полей, так и продольные компоненты Ег, Нг.
Если в [4] были получены поперечные компоненты электромагнитного поля для данной области
2 jga 1
22 g+2 g2 ed
jg2 1
Ex=-
F = -j = g+g2 ed
8Fz 1 wm c2 8Hz jw2ek (8Fz g 8Hz}
-+-
g a2 8j
a
we
8Fz wm c2 8Hz jw2sk
+
ga
a
8Fz , g 8Hz 1
+
we
Hx =
Hj=2
jg2 1
22 g+ g- ed
we 8Fz 8HZ + jw2ek
g 8j 8X
a
we 8Fz + 8Hz
jg2 1
22 g+ g- ed
we 8Fz + 8Hz jw ek
g 8x 8 j 0 we 8E 8H„
(1)
g 8X 8j a2
g
то продольные компоненты Ez электрического поля и Hz магнитного поля неизвестны. Здесь у — постоянная распространения,
2 2.2 / 22 2 2 / 2 g+=W S/ + W sk -g , g-=W SjU-W sk -g ,
d = ch2^ - cos2 p = -J0,5(ch2£ - cos 2 p) — геометрический параметр,
2 2 2 2 2 2 a = w s/u-g , e — фокусное расстояние эллипса, c = w ju±s -g .
В статье ставится задача определения продольных компонент электромагнитного поля.
1 Задача Дирихле для уравнения Гельмгольца
Для решения поставленной задачи необходимо решить следующую задачу Дирихле для уравнений Гельмгольца ЕН-обыкновенной и НЕнеобыкновенной волн [5]
82 FZ 82 F.
2-Z + e2 d2 (w2em± - g2) FZ - je2 d2 gam«— Hz = 0,
8f ^ ' " m
8X2 8f
82HZ ■ 82HZ , e2d2
где m± =
8X2 8f m2 - k2 m
m
2 Hi 2
w emii —g v m 0
(2)
HZ 2 je2 d 2gws— Fz = 0,
m
z
z
z
При этом условие Дирихле для электрического поля на границе бесконечно проводящей эллиптической области будет иметь вид:
ЕЛ=Хо = ЕЛ=Хо = (3)
Отметим, что система дифференциальных уравнений (2) описывает распространение гибридных ЭМВ, возникающих из-за гиротропности области распространения [1; 2; 6].
Для решения системы (2) преобразуем ее, применяя метод укорочения исходного дифференциального уравнения [1; 7]. Для этого умножим второе уравнение (2) на ]л и сложим с первым
а2 (е2 + ]аиг )д2 (е + ]аи2)
ах2
+
- + в2 ё2
' 2 т 2 ^
ю ет--у
т 0
№2 +
V
(4)
+ в2ё2(а2е/и1 - у2)е2 -Лв2ё2уюе —Е2 -]в2ё2ую/л,,—И2 = 0.
т т
Введем обозначение
Ег + ]ЛИ2 =у. (5)
Подставив (5) в (4) и учитывая, что Е2 = ]ЛИ2, после компоновки относительно ¥ и Е2 , получим:
а2¥ а2¥ -+-
ах2 ар2
+в2ё2
—
—
22 ю ет -у -Луюе—
¥ =
т0
(
-Л2в2ё2 усе—Лв2ё2
т
2 П\ 2 2 ,2
ю ет — У -с ет +У
т
л
—
22
+в ё уют ~
0
т
(6)
Н.
В правой части (6) имеем квадратное уравнение относительно Л ь ( ,, \ ь
уюе—Л +
т
2 т\\ 2 2 ,2 ю ет —-у -ю ец,^ + У
т
Л - уют. — = 0. (7) 0 т
В предположении, что Л12 являются корнями уравнения (7), формула (6) принимает вид:
а2^ , а2¥
+-
2 + в2 а2
г
ах2 а(р2
—
Л
ю ет± - у - Л12уюе — Т12 = 0. (8)
т0 '
Подставляя корни Лх и Л2 в (8), получим два его решения ¥1 и ¥2 . Тогда из (5) будем иметь:
Е + ЛИ, =¥; [Е2 + ]Л 2 И2 =¥2.
Из (9) получим:
2
E = ^Л2 -YA z Л2 -Л Y -Y
Hz = j —1—Y2.
Z Л 2 -Л1
(10) (11)
При известных Y1 и Y2 уравнение (6) может быть записано
а 2y1 а 2—1 ^2 œ
-2j +-2- + e 2d2
ах2 dj2
<a2sju1 - g2 - Lggwe — Y1 = 0,
u)
а2— а2—, +
+ e2d2
ах2 аj2
Уравнения Гельмгольца (12) после подстановки
d2 = 1 (ch2X - cos2j)
w 2eUi - g2 - Л2gws — Y2 = 0.
u)
(12)
примут вид:
а2— а2— e
— (ch2X - cos 2 j) w 2eu1 - g2 - Л-gœs —
+
ах2 аj2
+
л
2
а2y2 + а2—2 e
U0 k ö
Y1 = 0,
(13)
+--(ch2X - cos 2 j) w 2eu1 - g2 - Л2gws —
—2 = 0.
д£2 дф2 2 ^ ¡и0
Далее уравнения (13) решаются методом разделения переменных [8]. Решения будем искать в виде
^ =4^1, Y2 = ^2, (14)
где Yx1 и Yx2 — функции, зависящие только от X , а Yj1 и Yj2 — только от ф .
Подставив (14) в уравнения (13) и разделив первое уравнение (13) на Yx1Yj1, а второе — на Y£2Yj2, получим:
1 дХ в2 / 2 2 . к}
и)
———^+"V ch2x < 2£Ui- g2 - ^^ ~
YX1 ах 2 è
+
+
1 а 2—л e
f
k
Л
—л аj2
1 а2—^, e2 œ
--cos2j w eu1-g -Л^усоа —
2 è U0
= 0,
y—аХ2+"Vch2x ®2sui-g2 -Л№£—
Y Х2 аХ 2 è
k ö U0
(15)
+
1 а— e
k
Yj2 аj2
--cos2 j w su1 - g - Л^оа —
2 è U0
2
2
2
2
В (15) выражения в первых квадратных скобках зависят только от £ , а во вторых квадратных скобках — только от р .
Уравнения системы (15) могут выполняться только в том случае, если выражения в квадратных скобках в каждом из уравнений по отдельности будут равны одной постоянной величине, но с разными знаками, т. е.
¥ _ £1 5£2
1 5 %
¥ _ р1 5р2
1 5 %
¥ 5£2
1 5 2¥ р2
1 5% е2 Г 2 2 Л _ ^ ■ +— еп2Х со ац^ - у - Луута —
ц0 _ ^ Ц0
= Ь,
2 2
- — С0Б2р Ю2£Ц1-у2 -Л1ус£ —
(16)
= -Ь.
_
+—сИ2£ с 2ец1 - у2 - Лууюе —
Ц0
= 5,
Р р2 2
е 2 2 к _ --С0Б2р с ац1-у -Лута —
\ Ц0
(17)
= -5.
Умножив первую формулу в (16) на , вторую — на¥р1, получим:
52Ч£
52¥
5р
£1 -(Ь - 2д1сИ2£)х¥£1 = 0, Р + (Ь - 2^1С082р)Рр1 = 0,
(18)
где
е2
(
_2 - у2 -Луоа
Л
ц
4
(19)
и _2 = с 28Ц± , = + 1ЛхН2£ , ¥р1 = Е2р + 3ЛхН2р .
Умножив первую формулу в (17) на , а вторую — на¥р2, полу-
чим:
5
£2
5Х
5
р2
5р2
-(^ - 2ч2СИ2£)Ч>Х2 = 0, + (5 - 2ч2С082р)чр,2 = 0,
(20)
где
и
в 2
12 2 л к
к^-у - Л2уа>8
Чг =-4-и (21)
и ^2 = Ех + ]Л2Нч , ^,2 = Еф + ]Л2Нф.
В системах (18), (20) вторые уравнения являются обыкновенными уравнениями Матье, а первые — модифицированными уравнениями Ма-тье [3; 9]
Л =-и- (к1 -у2 - % Л2 =-и- (к!-у2 - % (22)
уюЕк в уюЕк в
Далее, представляя решения (18) и (20) в виде Y = YxYj согласно [3; 9], получим частные решения Yj = свт (ф, 2) или sвm (ф, Ч12) с постоянным множителем, а для Yx = Свт (X, 2) или Бвт (X, Ч12 ). Здесь свт (ф, 2) и sвm (ф, Ч12) — четная и нечетная, соответственно, периодические обыкновенные функции Матье целого порядка т с действительными Ь (или 5) и Ч .
Тогда, учитывая принцип суперпозиции, получим общие решения (18), (20) [3; 9]
¥
^,2 = ^Л,2 = / Ст1,2Свт (X Чх,2 )свт 41,2 -У) +
(23)
,2
т=0
/ 8т1,2$вт Х 41,2 )™т (ф 41,2 )С0^ - у),
+ / 1 ит1,2 8вт
т=1
где Ст12, 8т12 — произвольные постоянные.
Для любого т имеется два типа решений (четные и нечетные) [3; 9; 10]
С^1,2 =С^Х 1,2 = Ст1,2Свт (Х, Ч1,2 )Свт (ф, Ч\,2 ),
8%,2 =3Х¥Х 1,2 8^1,2 = 8т1,2^т (Х, 41,2 )5вт (ф 41,2 ),
где СY12 — четные решения, а 2 — нечетные.
В связи с тем, что в эллиптической области волны делятся на четные и нечетные, то (10) и (11) после объединения примут вид:
YЛ-YЛ YЛ-YЛ
Е = С х С х 2^1 Е = 8 1 2 8 х21
С ^ Л2-Л, '8 ^ Л2-Л, '
21 21 (25)
Y— Y Y— Y
Я- СТ1 Ст 2 и 8Т 1 8Т 2
7 = ]-, 8И7 = ]-.
С 7 Л2-Л! 8 7 Л2-Л!
Выражения (25) с учетом (24) примут вид:
СЕ2 Л Л (С¥1Л2 С¥2Л1) Л Л (с¥£1 С¥р1Л2 С¥£2 С¥р2Л1)
= Л - Л [Л2Ст 1Сет Ъ )сет Р Ъ ) - Л1Сш2Сет Ъ )Сет Р Ъ )1 Л2 Л1
Е =_1_( - ¥ Л ) =_1_( ¥ ¥ Л - ¥ ¥ Л ) =
5^2 Л Л \5 х1^2 5 х гМ/ Л Л \5 Х£15 хр1^2 5 х£2 5 х р2^ М /
= Л - Л [Л25т15ет (£, Ч1 )5ет Р Ъ ) - Л15т25ет Ъ Кт Р Ъ )1 Л2 Л1
Н =_^_( ¥-¥) =_^_( ¥ ¥ - ¥ ¥ )=
СА12 Л Л 1 Л Л £1С р1 С"1 £2 С ±р2/
]
Л2 -Л1
[Ст1Сет (£, Ч1 )сет Р Ъ ) - Ст2Сет (Х, Ъ )сет Р Ъ )1
Н =_^_( ¥-¥) =_^_(¥ ¥ ¥ ):
^г Л Л \5 А1 5"1 2) Л Л \5 £1С р1 5"1 £2 С ур2)
]
Л - Л " [5т15ет (X, Ч1 )*ет (Р Ъ ) - 5т25ет (£, Ъ Р Ъ )] (26)
Л2 Л1
Выражения (26), определяющие продольные компоненты электромагнитного поля, являются общим решением дифференциальных уравнений (2).
Подставив (26) в (1) и добавив продольные компоненты, получим выражения для всех шести составляющих (четных) гибридных волн
СЕ£ -т-г х
£ яV (Л2 -Л1)
х [ ]Сш1Сеш (£, Ъ ) Сет Ъ ) {у 2с_ - уа2Л2 } --}Сш2Сет (£, Ч2 ) Сет Ч2 ) ^ - уа"Л1 } + +Ст1Сет Ъ ) Сет (Ф, Ъ ) {цСС ' - уС '£_Л2 } -
- Ст2Сет (£, Ч2 ) Сет (Ф, Ч2 ) {цСС" - УС 2£_Л1 }], (271 )
:[]Ст1Сет (£, Ъ ) Се>т (Ф, Ъ ){у"®_ - уа2Л2 }-
С Ф_ яV (Л2 -Л1)[
- ]Сш2Сеш (£, Ъ2 ) Сет (Ф, Ъ ){у2®_ - уа2Л1 } -"Ст1Сет Ъ ) Сет (Ф, Ъ ){цСС ' - у^Л2 } +
+ Ст2Сет (£, Ъ2 ) Сет (Ф, ^2 ){цСС' - У^Л1 } ], (272)
1
Е7 = ■ 1 ■ [Л2Ст1Свт (Х , Ч ) Свт (Ф, Ч ) - Л1Ст2Свт (Х , 4 ) Свт (Ф, 4 )] ,(273)
Л2 -Л1
СН X = 2 2 ,/А-Лч []Ст1Свт Ч )свт ( ф Ч ^^'Л 2 - У® ^
вё (Л 2 -Л1 )
Н =-
- ]Ст2Свт (X, Ч2 ) Свт (Ф, Ч2 ){®8а 'Л1 - } +
+Ст1Свт (^ Ч1 ) Свт (Ф, Ч ) {Уа" - ^кЛ2 } -
-Ст 2Свт (Х, Ч2 ) Свт (Ф Ч2 ){уа" - ® "к Л1 } ], (274)
"[ 3Ст1Св'т (X, Ч1 ) Свт (Ф, 4 ) {у® ^ - ®£а"Л2 } -
С Ф" g-вё (Л2 -Л1 )1 - ]Ст2Свт (х, Ч2 ) свт (Ф, Ч2 ) {у® ^ - Л1 } + +Ст1Свт (Х, Ч1 ) Свт (Ф, Ч ) {уа" - ®8 "кЛ2 } -
Ст 2 Свт (X, Ч2 ) свт (Ф, Ч2 ^ -®3£2кЛ1}], (275)
СН 7 = - [ Ст1Свт (4, Ч1 ) Свт (Ф, 4 )-Ст2Свт (I , Ч2 ) Свт (Ф, Ч2 )]. (276) Л2 Л1
Здесь Ст1, Ст2 — амплитудные коэффициенты, свт (ф, ^ 2) — обыкновенная функция Матье 1-го рода целого порядка т, Свт (X, Ч12 ) — присоединенные (модифицированные) функции Матье 1-го рода (с целым индексом) и свт (ф, 2), Св (X, 2) — производные функций Матье 1-го
рода. Параметры функций Матье , Ч2 и корни (7) Л12 определяются,
соответственно, по формулам (19), (21) и (22).
Выражения для нечетных волн получаются аналогично, только необходимы следующие замены
Ст1 ® 8т1, Ст2 ® 8т2, Св ® 8в, св ® яв.
т1 т1 т 2 т 2
Для определения произвольных постоянных Ст1, Ст2 в (27) используем граничные условия (3). Из (3) следует, что тангенциальные составляющие электрического поля равны нулю на стенках (X =X0), ограничивающих эллиптическую область, т. е.
I СЕ7 =СЕф = 0,
Е7 = Еф = 0 ф (28)
7 ф 13Ег =8Еф= 0. ( )
Применяя условие (28) к (27) для Ег и Ер, получим:
СЕФ = Ст1 []Сет Ъ )Се>т (Ф,Ъ ) ^ С_ - уа2 } -
-Сет (£0 , Ъ ) Сет (Ф, Ъ ) {с2 - уС 2£_Л2 }] -Ст2 []Сет (£,д2)сеи (Ф,д2){у2с_ - уа^ } - (29)
-Сет (£0 , Ъ2 ) Сет (Ф, Ъ2 ) {цСС' - уС ^Л1 }] = 0,
СЕг = Ст1Сет Ъ ) Сет (Ф, Ъ )Л2 - Ст2Сет Ъ )Сет (Ф, Ъ )Л1 = 0.
Равенства (29) являются системой линейных однородных алгебраических уравнений относительно амплитудных коэффициентов Ст1 и Ст2.
Для существования нетривиальных решений системы (29) ее определитель должен равняться нулю
Сет (£0 , Ъ )Сет (ф, Ъ ) Л1 []СеИ (£0, Ъ ) СеИ (Ф, Ъ ) {у С_ - уа'Л2 } -
- Сет (£0 , Ъ ) Сет (ф, Ъ ) {Ц®С2 - уС£_Л2 }] Сет (£0 , Ъ )СеИ (Ф, Ъ )Л2 Х
[ ■ г /J (30)
Х []Сет (£0 , Ъ ) Сет (ф, Ъ ) у С_ - уа2Л1 } --Сет (X0, Ъ ) Сет (ф, Ъ ) {цс С" - ус ^Л1 }] =
В (30) вынеся за скобки Сет (£0, ъ ) в первом слагаемом и Сет (£0, ) — во втором, полученный результат разделим на
(Р, Ъ ) Сет (Р, Ъ )
се
т
- Л, ^^ Ъ|) {цсС - уСекЛ2}+Л2 Cе-^(Xo, ) {цсС - ую2г_Л,}+ 1 Сея (£0, ъГ ' 2 (£0, Ъ У 1 (31)
+]СР {у2с_ -уа2Л2}\1 -МЫ {у2с_ -уа2Л1(Л2 = 0. "" (Р Ъ1) Сет (Р, Ъ )
се
Вначале преобразуем первые две слагаемые величины в левой части (31)
. Се (^qi), 2 2 -Л,——,-т{¡июе -gas
1 Cem (Х qi )
Л Ce (x0, qi) f = -Л -m,.—г 1 ma
. Cem (X0, q2 )i 2 2
кЛ 2 } + Л 2——-- {mac -ga2s
' t^ r> I > /1 \ >
Cem (Xo, q2 )
kLi}:
+Л.
1 Cem (X0, qi )
Cem (x q2)
f k ö r 2 2 л A
k±-g -gas—Л2
m
(32)
mw
( k > < 2 2 л a
k±-g -gas—Л_
m /
Cem (X0, Я2 )
В (32) выражения в круглых скобках выразим через q1 (19) и q2 (21), а затем вместо Л1 и Л2 подставим их значения (22). Тогда получим:
. Ce (x qi ), 2 2 ,А ) A Cem (Х ^2 )i 2 2 ,A )
-Л^ТТ^-\ {/wc -gw2skЛ2} + Л2——-r{u®c -gw2skЛ1} =
Cem (x qi)
Cem (X q2 )
= -fk2_g2 -ö_ml4qiCe (xqi).
I 1 e" Jgsk e" Cem (^ qi )
4q2 ö m2 4qi Cem q2 )
(33)
V 1 g e" 0gsk e" Cem (X0, q2)' Далее преобразуем третье и четвертое слагаемые выражения в (3i)
Jcf) {g2ak -ga 2 Л 2 }Li - jemf {g2ak -ga2 Л1}л 2 =
ce.
(ф qi) cem (ф q2)
Jcem (Ф, q2 ) Jcem (ф qi )f +
= ga2 Л1Л
cem (ф, q2 ) cem (ф, qi )
(34)
2 a I Jcem (ф, qi ) Л - jcem (ф, q2 )
+g ak
cem (ф qi ) 1 cem (ф qi )
Л,
gsk
Подставив (33), (34) в (31) и умножив полученное выражение на-.
m
имеем:
-fk2±_g2 - 4q_ö Ce (xo,qi)+fk2_g2 - 4^ ö4Ce (xq2)
+j
e" 0 e" Cem (X0, qi ) g2ska2Л1Л2 f cem (ф, q2 ) - cem (ф, qi ) "
m2 1 cem (ф q2) cem (ф, qi),
g3ask f cem (ф qi ) Л - cem (ф, q2 )Л
m2 1 cem (ф, qi) 1 cem (ф,q2) 2
e" 0 e" Cem (X q2)
(35)
Подставляя в (35) Л1 и Л2 из (22), получим
2 -
-Гк2 - у2 - Ч1^ (Х 41 ) Г , 2 2 - 492.1 Ч_ ^ (Х 42 )
I 1 У е2 0 е2 Сет (Хо,41) Г 1 у е2 1 е2 Сеи (X,42)
се.
ф 41)
^ - у2 - 4411 - ^ (Ф 42 ) 1^1 - у2 - 442 ^ [ Сет (ф 41 )Г 1 е" 1 Сет (Ф, 42 )Г 1 е"
сет (Ф, 41)/ (36) = 0.
Формула (36) представляет собой дисперсионное уравнение четных волн, полученное в результате решения задачи Дирихле (2), (3). Отметим, что (36) имеет вид, аналогичный полученному для кругового гиротропного цилиндра [6].
Чтобы получить дисперсионное уравнение для нечетных волн, в (36) надо сделать замены
,41,2) ® ^(Х, 412) Се(х0,4],2) ® 5е (х0,4],2), ^ 41,2 ) ® se(j, 41,2 ), се 41,2 ) ® sejP, 41,2 ) где 0 , Я-\2) и 5е(Х0,4и) — нечетные присоединенные (модифицированные) функции Матье1-го рода (с целым индексом) и их производные, se{j.,412) иse (р,412) — нечетные обыкновенные функции Матье 1 -го рода целого порядка т и их производные.
Заключение
Решение задачи Дирихле для уравнений Гельмгольца в гиротропном волноводе в строгой постановке приводит к большим затруднениям из-за громоздких выражений для компонент электромагнитного поля. В то же время в работе [11] для круглого продольно намагниченного ферритового волновода показано, что решения краевых задач методом укорочения исходного дифференциального уравнения полностью совпадают с решениями задачи в строгой постановке и образуют полные решения краевых задач для ферритовых волноводов.
В данной работе впервые осуществлена адаптация метода укорочения исходного дифференциального уравнения для случая эллиптического волновода с ферритовым заполнением при продольном намагничивании. Полученные результаты позволяют исследовать различные характеристики распространения ЭМВ в эллиптическом гиротропном волноводе при продольном намагничивании и поставить задачу идентификации параметров намагничивающего магнитного поля и управляемых переменных, характеризующих «свойства заполняющей пространство ферритовой среды», в заданном классе функций. В частности, провести исследование зависимости постоянной распространения ЭМВ от свойств заполняющей волновод ферритовой среды и параметров намагничивающего магнитного поля.
Основные выводы таковы:
1. Поставлена и решена задача Дирихле для уравнений Гельмгольца с целью нахождения компонент электромагнитного поля в гиротропном эллиптическом цилиндре с продольным намагничиванием.
2. На основе решения задачи Дирихле получены дисперсионные уравнения ЭМВ, позволяющие проводить исследования распространения гибридных в эллиптическом гиротропном волноводе, такие как установление зависимости постоянной распространения ЭМВ от свойств заполняющей волновод ферритовой среды и параметров намагничивающего магнитного поля.
Литература
1. Сул Г., Уокер Л. Вопросы волноводного распространения электромагнитных волн в гиротропных средах: пер. с англ. М.: Изд-во иностр. лит-ры, 1955. 192 с.
2. Лакс Б., Баттон К. Сверхвысокочастотные ферриты и ферримагнетики: пер. с англ. М.: Мир, 1965. 676 с.
3. Мак-Лахлан Н. В. Теория и приложение функций Матье / пер. с англ. В. А. Братановского. М.: Изд-во иностр. лит-ры, 1953. 475 с.
4. Итигилов Г. Б., Ширапов Д. Ш. Метод инвариантных преобразований для определения поперечных компонент электромагнитного поля в гиротропных ограниченных областях // Вестник Бурятского государственного университета. 2012. Вып. 9: Математика, информатика. С. 162-166.
5. Ширапов Д. Ш., Итигилов Г. Б. Обобщенные уравнения Гельмгольца ги-ротропных волноводов произвольной формы поперечного сечения // Современные проблемы дистанционного зондирования, радиолокации, распространения и дифракции волн: материалы II Всерос. науч. конф. (г. Муром, 26-28 июня 2018 г.). Муром, 2018. С. 209-219.
6. Микаэлян А. Л. Теория и применение ферритов на сверхвысоких частотах. Л.: Госэнергоиздат, 1963. 664 с.
7. Назаров А. В., Раевский С. Б. Электромагнитные волны в структурах, содержащих продольно намагниченные ферритовые слои // Физика волновых процессов и радиотехнические системы. 2007. Т. 10, № 1. С. 76-82.
8. Бицадзе А. В. Уравнения математической физики. 2-е изд., перераб. и доп. М.: Наука, 1982. 336 с.
9. Анго А. Математика для электро- и радиоинженеров. М.: Наука, 1967. 780 с.
10. Бейтмен Г., Эрдейи А. Высшие трансцендентные функции. Эллиптические аморфные функции. Функции Ламэ и Матье / пер. с англ. Н. Я. Виленкина. М.: Наука, 1967. 300 с.
11. Назаров А. В., Новоселова Н. А., Раевский С. Б. О полноте системы решений краевых задач для ферритовых волноводов, полученных методом укорочения дифференциального уравнения // Антенны. М.: Радиотехника, 2016. Т. 7(227). С. 63-66.
DIRICHLET PROBLEM FOR HELMHOLTZ EQUATIONS
IN GYROTROPIC ELLIPTICAL REGION WITH LONGITUDINAL
MAGNETIZATION
Dashadondok Sh. Shirapov Dr. Sci. (Phys. and Math.), Prof., Dorzhi Banzarov Buryat State University 24a Smolina St., Ulan-Ude 670000, Russia E-mail: [email protected]
Garma B. Itigilov
Cand. Sci. (Engineering), A/Prof.,
East-Siberian State University of Technology and Management 40v, bldg 1 Klyuchevskaya St., Ulan-Ude 670013, Russia E-mail: [email protected]
Igor B. Yumov
Cand. Sci. (Phys. and Math.), A/Prof., Dorzhi Banzarov Buryat State University 24a Smolina St., Ulan-Ude 670000, Russia E-mail: [email protected]
Vladimir D. Anakhin
Dr. Sci. (Engineering), Prof.,
Dorzhi Banzarov Buryat State University
24a Smolina St., Ulan-Ude 670000, Russia
E-mail: [email protected]
Zhargal G. Dambaev
Dr. Sci. (Engineering), Prof.,
Dorzhi Banzarov Buryat State University
24a Smolina St., Ulan-Ude 670000, Russia
E-mail: [email protected]
We have formulated and solved the Dirichlet problem for Helmholtz equations of electromagnetic waves, propagating in an elliptical cylinder filled with longitudinally magnetized ferrite which is described by a second-rank tensor. It is assumed that the cylinder has an infinitely conductive wall. To solve the boundary value problem of Helmholtz equations for longitudinal components of electromagnetic waves, we have used the method of shortening the initial differential equation and the method of variables separation. The solution of the above boundary value problem in elliptic coordinates is associated with the use of even and odd ordinary and modified Mathieu functions of the first kind. Using the results obtained, we have determined all the components of electromagnetic waves for even and odd solutions. Applying the Dirichlet condition to the components of electromagnetic waves and solving the system of linear homogeneous algebraic equations, we have obtained dispersion equations. The found dispersion equations of electromagnetic waves are of great practical importance and allow studying the propagation of hybrid waves in this region.
Keywords: elliptical cylinder; ferrite; Dirichlet problem; Helmholtz equation; electromagnetic wave; longitudinal magnetization; gyrotropic region; transverse components of the electromagnetic field; Mathieu functions; dispersion equation.
References
1. Suhl H. and Walker L. Waveguide Propagation of Electromagnetic Waves in Gyrotropic Media (Russ. transl.). Moscow: Inostrannya literatura, 1955.
2. Lax B., Button K. J. Microwave Ferrites and Ferrimagnetics. New York: Mac-Graw-Hill, 1962. 752 p.
3. McLachlan N. W. Theory and Application of Mathieu Functions. London: Oxford University Press, 1947.
4. Itigilov G. B., Shirapov D. Sh. Metod invariantnykh preobrazovanii dlya opredeleniya poperechnykh komponent elektromagnitnogo polya v girotropnykh ogranichennykh oblastyakh [Method of Invariant Transformations for Determining the Transverse Components of Electromagnetic Field in Gyrotropic Bounded Areas]. Vest-nik Buryatskogo gosudarstvennogo universiteta. Matematika, informatika. 2012. V. 9. Pp.162-166.
5. Shirapov D. Sh., Itigilov G. B. Obobshchennye uravneniya Gelmgoltsa giro-tropnykh volnovodov proizvolnoi formy poperechnogo secheniya [Generalized Helm-holtz Equations of Gyrotropic Waveguides with Arbitrary Cross-Section]. Sovremen-nye problemy distantsionnogo zondirovaniya, radiolokatsii, rasprostraneniya i di-fraktsii voln. Proc. 2nd All-Russ. conf. (June 26-28, 2018). Pp. 209-219.
6. Mikaelyan A. L. Teoriya i primenenie ferritov na sverkhvysokikh chastotakh [Theory and Application of Ferrites at Microwave Frequencies]. Leningrad: Gosener-goizdat Publ., 1963. 664 p.
7. Nazarov A. V., Raevskii S. B. Elektromagnitnye volny v strukturakh, soderz-hashchikh prodolno namagnichennye ferritovye sloi [Electromagnetic Waves in Structures Containing Longitudinally Magnetized Ferrite Layers]. Fizika volnovykh protses-sov i radiotekhnicheskie sistemy. 2007. V. 10. No. 1. Pp. 76-82.
8. Bitsadze A. V. Uravneniya matematicheskoi fiziki [Equations of Mathematical Physics]. 2nd rev. ed. Moscow: Nauka Publ., 1982. 336 p.
9. Ango A. Matematika dlya elektro- i radioinzhenerov [Mathematics for Electrical and Radio Engineers]. Moscow: Nauka Publ., 1967. 780 p.
10. Bateman H., Erdelyi A. Higher Transcendental Functions. New York: McGraw-Hill, 1953-1955.
11. Nazarov A. V. Novoselova N. A., Raevskii S. B. O polnote sistemy reshenii kraevykh zadach dlya ferritovykh volnovodov, poluchennykh metodom ukorocheniya differentsialnogo uravneniya [On Completeness of the System of Boundary Value Problems Solutions for Ferrite Waveguides Obtained by the Method of Differential Equation Shortening]. Telecommunications and Radio Engineering. 2016. V. 7 (227). Pp. 63-66.