CC BY
68
2. МатейД. Л. Фильтры СВЧ, согласующие цепи и цепи связи / Д. Л. Матей, Л. Янг, Е. М. Т. Джонс ; [пер. с англ. под ред. Л. В. Алексеева]. - М.: Связь, 1972. -Т. 2. - 495 с.
3. Марков Г. Т. Возбуждение электромагнитных волн / Г. Т. Марков, А. Ф. Чаплин. - М. ; Л. : Энергия, 1967. -367 с.
4. Миллер М. А. Исследование понятия поверхностного импеданса в теории поверхностных волн (обзор) / М. А. Миллер, В. И. Таланов // Изв. вузов. Радиофизика. - 1961. - Т. 4, № 5. - С. 795-830.
5. Арсеньев Е. М. Поверхностные волны над трапециевидно гофрированной поверхностью / Е. М. Арсеньев, В. А. Калошин // Журнал радиоэлектроники. -2008. - № 4. - С. 127-135.
6. Горобец Н. Н. Приближенное решение дисперсионного уравнения для цилиндрического волновода с импе-дансными границами / Н. Н. Горобец, А. П. Удовенко // Радиофизика и радиоастрономия. - 2005. - Т 10, № 2. -С. 172-177.
7. Петрова Ю. Ю. Решение уравнения Гельмгольца в многосвязных волноводных областях: дис. ... канд. физ.-мат. наук : 01.01.03 / Петрова Юлия Юрьевна. -М., 2006. - 114 с.
8. Моденов В. П. Проекционные методы в теории волноводов / В. П. Моденов // Радиотехника и электроника (Россия). - 2005. - Т. 50, № 2. - С. 203-207.
9. Богданов Ф. Г. Дифракция волны Н10 на симметричных диэлектрических стержнях конечной длины / Ф. Г. Богданов // Изв. вузов. Радиофизика. - 1983. -Т. 26, № 2. - С. 246-250.
Надійшла 02.06.2009 Після доробки 03.09.2009
Логачова Л. М, Куцак С. В., Бондарєв В. П.
ДИФРАКЦІЯ ХВИЛІ #10 В ПРЯМОКУТНОМУ ХВИЛЕВОДІ НА СТРИБКУ ІМПЕДАНСА ВУЗЬКИХ СТІНОК
На основі проекційного методу дається строге рішення задачі дифракції основної хвилі Н10 на стику регулярного і нерегулярного (вузькі стінки утворені вторинними хвилеводами) прямокутних хвилеводів. На основі чисельних результатів приводяться та аналізуються залежності дифракційних характеристик системи від її параметрів.
Ключові слова: дифракція, імпедансні стінки, проекційний метод, прямокутний хвилевід.
Logacheva L. M., Kutsak S. V., Bondaryev V. P.
Н10 WAVE DIFFRACTION IN A RECTANGULAR WAVEGUIDE BY NARROW WALLS IMPEDANCE BUMP
Rigorous solution of the principal H10 wave diffraction problem at the joint of regular and non-regular (narrow walls are formed by secondary waveguides) rectangular waveguides has been obtained using the method of projection matching. On the basis of numerical results the dependence of the system diffraction characteristics on its parameters is determined and analyzed.
Key words: diffraction, impedance walls, method of projection matching, rectangular waveguide.
УДК 621.372.81
Самойлик С. С.1, Бондарев В. П.2
1 Аспирант Запорожского национального технического университета 2Канд. физ.-мат. наук, доцент Запорожского национального технического университета
ЭНЕРГЕТИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПРОХОДНОГО РЕЗОНАТОРА НА ПРЯМОУГОЛЬНЫХ ВОЛНОВОДАХ С ЧАСТИЧНЫМ ДИЭЛЕКТРИЧЕСКИМ ЗАПОЛНЕНИЕМ
Расчет параметров проходного прямоугольного резонатора с частичным диэлектрическим заполнением при сильной связи выполнен методом частичных пересекающихся областей с использованием интегральных уравнений макроскопической электродинамики. Этот метод позволяет учитывать как параметры каждого из диэлектриков, такие как диэлектрическая проницаемость и тангенс угла потерь вещества, из которого изготовлена неоднородность, так и его геометрическую форму и место расположения в резонаторе. Проведенный численный анализ позволил вычислить мощности, поглощенные неоднородностями, амплитуды прошедшей через структуру и отраженной волн в зависимости от параметров сложного диэлектрического заполнения.
Ключевые слова: проходной резонатор, коэффициент отражения, коэффициент прохождения.
ВВЕДЕНИЕ ний, а необходимо математическое моделирование
Моделированию процессов, происходящих в СВЧ этих щювдссот. Одн°й из теркмых проблем при раз-
резонаторах, заполненных диэлектрическим матери- работке таких устройств является трудность формиро-
алом, посвящены исследования [1, 2]. Обзор работ по- вания заданных распределений поля в диэлектри-
казывает, что для детального анализа наблюдаемых ческом материале, помещенном в резонатор. Наиболее
явлений недостаточно экспериментальных исследова- предпочтительным является аналитическое решение
© Самойлик С. С., Бондарев В. П., 2010
задач СВЧ нагрева, которое удобно для нахождения зависимости распределения температуры от параметров среды и характеристик системы возбуждения электромагнитного поля и их оптимизации. Например, в работе [3] получено аналитическое решение задачи определения температурного поля и глубины проникновения электрического поля при нормальном падении электромагнитной волны на плоский диэлектрик с малыми потерями, заполняющий полупространство.
При моделировании электромагнитных и тепловых полей в СВЧ резонаторах открытого типа и, особенно, в резонаторах закрытого типа применяются, как правило, такие численные методы решения системы уравнений Максвелла, как метод конечных элементов (МКЭ), метод конечных разностей (МКР), либо метод частичных областей (МЧО). Недостаток методов МКР и МКЭ заключается в необходимости использования больших объемов памяти и времени счета на ЭВМ. Метод МКР также обладает одним очень серьезным недостатком, связанным с трудностью учета криволинейных границ диэлектрического материала. Метод частичных областей (МЧО) широко применяется при решении многих электродинамических задач, в том числе задач расчета резо-наторных электродинамических структур и волноводных трансформаторов. В МЧО исследуемая область разбивается на подобласти, затем в каждой подобласти электромагнитное поле раскладывается в ряд, как правило, по собственным функциям подобластей. Наложение граничных условий и использование свойств ортогональности собственных функций приводят к матричному уравнению относительно неизвестных коэффициентов разложения. В настоящее время существует много модификаций МЧО, отличающихся способами разбиения электродинамических структур на частичные области и сшивания полей на границах раздела.
В данной работе расчет параметров проходного прямоугольного резонатора с частичным диэлектрическим заполнением при сильной связи решается методом частичных пересекающихся областей с использованием интегральных уравнений макроскопической электродинамики [4, 5].
ПОСТАНОВКА ЗАДАЧИ И МЕТОД РЕШЕНИЯ
В качестве математической модели выбран проходной резонатор с сильной связью, содержащий диэлектрическое заполнение с заданными значениями комплексной диэлектрической проницаемости (рис. 1).
Рис. 1. Резонатор с диэлектрическим заполнением
Разобьем всю область определения электромагнитного поля на две пересекающиеся области (области пересечения обозначены пунктирной линией):
1-я область АБЄВ
—о < 2 < +о; 0 < х < а,
2-я область Е¥ОИ \0 < 2 < 1;
(резонансная область) | о < х < а + С.
Будем считать, что на входе данной конструкции вдоль положительного направления оси 2 распространяется волна ТЕ10. Область резонатора (вторая область) заполнена веществом с комплексной диэлектрической проницаемостью ¿2 и отделена от основного волновода диэлектрической пластиной с диэлектрической проницаемостью £ь Предположим, что основной волновод заполнен средой со скалярными проницаемостями £0 и ц0 и имеет идеально проводящую поверхность.
Электрическое поле запишем через Еу составляющую.
Воспользовавшись второй теоремой Грина и учитывая вторичные волны, переизлученные диэлектриками, а также граничные условия Дирихле, для первой и второй областей можно составить следующую систему интегральных уравнений:
Е (х, у) = Ест(х, 2) -
- |ЕП(X', 2')(Х,^,Х'’ 2/)С2'
(1)
ЕП( X, у) = | Е1 (х', 2')
,і, , ,^дО (х, 2; х', 2')
Э 2'
сх
-}е'(х, 2')Э° (хі2 х' >2'-Сх'
д 2'
+ ^ (к + graddiv) х
2' = І і =1
| Оп( х, 2; х', 2' )(е ,■ - 1) Еп( х', 2') ёУ. (2)
Электромагнитное поле в резонансной области определим из интегрального уравнения для £п( х, у):
х = а
а
г = 0
0
о
É\x, y ) =
J Ecтіx’, z') - JEIIіx", z")’ " x ’z )dzr'
x dG (x,z; x',z0dx dz
- J Eet(x', z' ) - JE11 іx", z") -f" X" ^ dz"
x ЭО (y; x/>z0dx
dz
dx"
к 2
+ у і к + graddiv ) x
z' = І i =1
х | О11 (х, 2; х', г')(е,■ - 1)Е11 (х', г')ёУ, (3)
V
где О!(х, г; х', г'), Оп(х, г; х', 2) - функция Грина для первой и второй области, соответственно.
Для первой области функция Грина имеет вид:
G1 іx, z ; x', z' ) =
1 - i s п Л • i s п Л -ії> - z'l
у---------sin — x sin — x ■ є .
*-ч ysa \ а / \ а J
(4)
Для второй области:
G і x, z ; x', z' ) =
-уу
, mп пп
sinl -------- x I sinl — z
\ а + d I і І
і а + d) І
,2 m п пп
к-Іа-пІ -17
, mп , . пп ,
x sinl----------- x I sinl — z
і а + d І і І
(5)
i,2 f s П , 2 п
где ys = Jk - [ aJ ; k = — - волновое число.
Уравнение (3) является интегральным уравнением Фредгольма второго рода с вырожденным ядром. Решение этого уравнения будем искать в виде разложения по собственным функциям второй области
£“( X, z) = £ £ Eae sin ^ aOf-d x) sin f ^ z |, (6)
а сторонний источник представим следующим образом:
's П А -‘V
Ест(x, z) = E0sin^—xJ • e
Подставим в уравнение (2) выражение (6), получим:
f f Л En(x,y) = J ECT(x', z')-J£EaPsinI ^x"|sin( ^z"J x
0 0 a, p ^ J
ЭGIіx', z'; x", z") , /
x ¡x"----------------) dz
dG іх- z; x\ z'-dx/
dz'
- J EcTі x', z' ) - |у EаßSІn і а+d x "I sin і
0 і 0 а, ß і 1
• i ß:2 111 ---7
z x
ЭGIіx', z'; x", z") , ,
x ¡x"----------------) dz'
Э G і x, z ; x', z'-dx
d z'
+ у (к2 + graddiv)| ОП(х, г; х', г')(е,■ - 1)Еп(х', г')ёУ'.
,■ =1 у1
(7)
Вычислив интегралы, входящие в (7), и введя следующие обозначения, получим систему неоднородных линейных уравнений [6]:
EIIі x, z ) = уу Emnsin і ^ xj sin ^ zl = = -у E0B 1„nI„*smі -.Цx! sin{n)-z I +
0 mn ms
а + d
І
і x-' T-' 7-» 1 • i ап Л / * -.s sп
+ у EaßB 1-nsin(.0+7!і-1)^ x
m, n, а, ß, s
т л . i m п Л . i n п Л , x Im^Ißssin^ а"—- x| sin^—z| +
+ у і -1 )nE0B 1.пЄ ‘^Оіпі -.Ц x I sin і '^z І -
а + d
І
x-1 7-' пі • i ап Л , і\s+n sп
£ E„, B 1.nsiniO+dO) і-1) —¡ x
m, n, а, ß, s
T T2 . i тп Л . i пп .
x I.sIßs sinx| sin^—z\ -
Eаß ■ B2mn ' Ai, mn, аß ■ к і^^i 1 ) x
i, m, n, а, ß
1 mп Л . i пп x sin------------ sin — z
і а + d І і І
(8)
где введены следующие обозначения:
B 1 =
mn
і а + d) І V І
пп
2пЛ2 i тп Л2 iппЛ2’
B2mn = B 1 mn^ і ^
/ пп
j r. i sп л ■ i m п л j ,
= Jsin 1,7x І 5,п{:^ x J •* =
s ^ , .s + ^ i m па Л —і-1 ) sin------------
а \а + d )
sпЛ2 i тп Л 2 а + d-
z = 0
z =0
z
mn
а
л г • (Вп "] -іу*2" -1 "
/р, = | 8ІП ^ 1у-2 | • Є С2 =
0
[ 1 + (-1 )р+^1 ] рп
Рп
/В, = | БІпГI--?2"| • Є
І
+ (у,)
‘и‘-2"’С2" = (-1 )Р +1/Р,,
А
і, тп, ар
- коэффициенты, отвечающие за геометрию /-го диэлектрика.
Решая полученную систему численным методом с использованием ЭВМ, найдем амплитуды напряженности поля во второй области.
Подставив выражение (6) с определенными амплитудами напряженности электрического поля в выражение (1), получим:
Е’( х, 2) = Е0 БІп( а х] • Є ^ - |££ Еар :
0 ар
, ап Л . (Рп Л дО (х, 2; х , 2 ) , ,
х 8ІП -------, х Sin 2 -------------‘-■г--—:----- С2
Vа + С І V І I дх'
. (9)
Выражение для прошедшей волны 2 > 2 принимает вид:
7^+1 / \ 7~" ’ (,П ] У,2 1 ,П
Е (х, 2) = Е0§ІП(^х| Є - Ш (-1 )ЕаР~а х
ар ,
[ 1 + (-1 )р+1ЄЬІ] ІІЕ . ап ] І . (,п ] -у,2
х 8ІП --------а ----------- -------------------8ІШ — х є .
а + С
рп
+ (у,)2
а
Выражение для отраженной волны 2 < 2' запишется в виде:
Е Ч( х,2) = -Ш(-1 У Еар--а
ар
1 ,п . ( ап
8Ш ----------а х
у, а а V а + С
[ 1 + (-1 )р+1е^1 ] х --------------- ---------------1- • 8ІП( ^х] Є*2.
Рп
+ (у,)2
ЧИСЛЕННЫЕ РЕЗУЛЬТАТЫ
По значениям амплитуд напряженности полей рассчитаны коэффициенты прохождения и отражения, а также мощность излучения из резонансной области отраженной (Ря) и прошедшей (РТ) волны, мощность потерь в диэлектриках (Р^). Мощность потерь в диэлектрике определяется выражением [7]:
Рсм = уЄ"Є0 | Е2СУ,
(10)
где є - мнимая часть диэлектрической проницаемости; юр - частота излучения; N - номер диэлектрической неоднородности N = 1, 2).
Мощность потерь в ^й неоднородности, с учетом (6) и (10), будет определяться выражениями:
п™ ь 1^
PсN = Ю-є" IЕ2СУ = 1 ює^ | | | \Еу\гСхСуС2 =
лт 0 т п1 г1
1
Л N) , (N)
П2 Ь 2
= 2 ю^ I I I
п(N) 0 т п1 г1
Ё ЕарЄІП(—х] • ЄІПV^2
ар = 1
І
СхСуС2.
После преобразований и операций интегрирования получим окончательные выражения для мощности потерь в каждой из диэлектрических неоднородностей:
Р СИ = 2 Ю£ЛТ£0 Т Т ЕаеЕт»7- • "^ир). (11)
а, р =1 т, п =1
Коэффициенты /^а, ■У'пв соответственно равны:
^ sin ({Ш—З-Ж ^п((т-+-а)п^)
/(N) = а
/та = -2---п--
(т - а)
(т + а)
^^н-^п^ 5іп V ^-.±а}2- п1N)||]
(т - а)
+-
(т + а)
да _
^пр
2п
(п - р)
(п + р)
sinV(n■.:::lШ.■Il1N))
(п - В)
(п + В)
Графики зависимости прошедшей, отраженной и поглощенной в диэлектрике мощности от длины волны представлены на рис. 2.
При этом были выбраны следующие параметры диэлектрического заполнения: жидкость внутри резонатора £2 = 10, tg5 = 5 • 10-2, отделенная от волновода кварцевой диэлектрической пластиной (£1 = 3, 8, tg5 = 1 • 10-4) с толщиной 0, 1 • а; геометрические размеры волновода 23х10 мм; толщина резонаторной области С = 0,5а.
На графике рис. 2 можно выделить три области, в которых коэффициент прохождения минимальный; соответственно в интервале частот 0,68.0,7,
0,74.0,77, 0,85.0,87 следует помещать диэлектрики, чтобы они поглощали максимальное количество
2
х = а
+
V
Рис. 2. Зависимость прошедшей, отраженной и поглощенной в диэлектрике мощности от длины волны
энергии. Соответственно, именно в этих областях будет происходить наиболее эффективный нагрев неоднородностей в резонансной области.
ВЫВОДЫ
Расчет проходного резонатора с диэлектрическими неоднородностями проведен с использованием строгой электродинамической модели. В работе в широком диапазоне частот рассчитаны зависимости прошедшей, отраженной и поглощенной в диэлектрике мощности. Показано, что в данной структуре возможно как явление резонанса (коэффициент отражения по мощности максимален), так и явление антирезонанса (коэффициент прохождения по мощности близок к единице). На частотах резонанса наблюдается максимальное поглощение энергии диэлектриком. Эффективность рассмотренного метода интегральных уравнений обусловливает перспективу его применения для расчета более сложных волноводных трансформаторов, содержащих произвольное число диэлектрических неоднородностей различной геометрической формы.
СПИСОК ЛИТЕРАТУРЫ
1. Афанасьев А. М. Математическое моделирование взаимодействия СВЧ излучения с влагосодержащими плоскими слоистыми средами. Ч. 1. / А. М. Афанасьев,
В. В. Подгорный, Б. Н. Сипливый, В. В. Яцышен // Изв. вузов. Электротехника. - 2001. - № 2. - С. 14-21.
2. Степанов В. В. Оптимизация равномерного распределения СВЧ энергии в частично заполненном резонаторе / Степанов В. В. // Материалы 10-й Международной Крымской конференции [«СВЧ-техника и телекомму-
никационные технологии»], (Украина, Крым, Севастополь, 11-15 сентября, 2000). - Севастополь : Севастопольский национальный технический университет, 2000. - С. 596.
3. Рикенглаз Л. Э. К теории распространения СВЧ электромагнитных полей в диэлектриках с малыми потерями // Журнал технической физики. - 1974. - Т. 44. -С. 1125.
4. ХижнякН. А. Интегральные уравнения макроскопической электродинамики / Хижняк Н. А. - К. : Наукова думка, 1986. - 279 с.
5. Прохода И. Г. Метод частично пересекающихся областей для исследования волноводно-резонаторных систем сложной формы. / И. Г. Прохода, В. П. Чумаченко // Изв. вузов. Радиофизика. - 1973. - Т. XVI, № 10. -
С. 1578-1582.
6. SamoylykS. S. Complex frequencies of the rectangular resonator with the arbitrary number of the dielectric heterogeneities / S. S. Samoylyk, V. P. Bondaryev // Proceedings of 6th International Conference [«Antenna Theory and Techniques »], (Ukraine, the Crimea, Sevastopol, September 17-21, 2007). - Sevastopol : Sevastopol National Technical University, 2007. - P. 152-154.
7. Капилевич Б. Ю. Волноводно-диэлектрические фильтрующие структуры : [шравочник] / Б. Ю. Капилевич, Е. Р Требехин. - М. : Радио и связь, 1990. - 272 с.
8. Kalesinskas V Analysis of Electromagnetic Waves Scattering by a Finite Size Dielectric and Metal Cylinder in a Rectangular Waveguide / A. Konstantinov, V. Shugurov // International Journal of Infrared and Millimeter Waves. -1998. - Vol. 19, № 12. - P. 1773-1781.
9. MelezhikP. N. Open Resonators with Conducting Cylindrical Inserts. 1. Two-Dimensional Model / V. S. Miroshnichenko, Ye. B. Senkevich // Radiophysics and Quantum Electronics. - 2005. - Vol. 48, № 7. - P 529-536.
10. Shakhbazov S. Yu. Measurement of the distribution of the temperature field of materials in microwave resonators of complex shape / M. V. Nefedov, V. E. Nikishin, D. A. Loik, A. O. Nikishev // Measurement Techniques. - 2008. -Vol. 51, № 6. - P. 671-674.
Надійшла 28.05.2009 Після доробки 31.08.2009
Самойлик С. С., Бондарєв В. П.
ЕНЕРГЕТИЧНІ ХАРАКТЕРИСТИКИ ПРОХІДНОГО ПРЯМОКУТНОГО РЕЗОНАТОРА НА ПРЯМОКУТНИХ ХВИЛЕВОДАХ ІЗ ЧАСТКОВИМ ДІЕЛЕКТРИЧНИМ ЗАПОВНЕННЯМ
Розрахунок параметрів прохідного прямокутного резонатора з частковим діелектричним заповненням при сильному зв’язку виконується методом часткових пересічних областей і використанням інтегральних рівнянь макроскопічної електродинаміки. Цей метод дозволяє враховувати як параметри кожного з діелектриків, такі як діелектрична проникність і тангенс кута втрат речовини, з якого виготовлена неоднорідність, так і його геометричну форму й місце розташування в резонаторі. Проведений чисельний аналіз дозволив обчислити потужність, поглинену неод-норідностями, амплітуди відбитої хвилі та хвилі, що пройшла крізь структуру, залежно від параметрів складного діелектричного заповнення.
Ключові слова: прохідний резонатор, коефіцієнт
відбиття, коефіцієнт проходження.
Samoylyk S. S., Bondaryev V. P.
ENERGY PARAMETERS OF TRANSMISSION RESONATOR BASED ON RECTANGULAR WAVEGUIDES PARTIALLY FILLED WITH DIELECTRIC
Parameters of a rectangular transmission resonator partially filled with dielectric and having strong coupling are calculated by the method of partially intersecting fields using integral equations of macroscopic electrodynamics. This method permits to consider the parameters of each dielectric such as permittivity and a tangent of loss angle of the obstacle material, as well as its geometrical shape and location in the resonator. The numerical analysis permitted to calculate power absorbed by the obstacles and the reflected and transmitted wave amplitude depending on the parameters of the complex dielectric filling material.
Key words: transmission resonator, reflection coefficient, transmission coefficient.
