Анализ влияния ступенчатых нерегулярностей внутреннего проводника на частотные характеристики кусочно-однородной брэгговской СВЧ структуры в коаксиальном волноводе
В.А. Ефимов, М.Э. Валъяров, Р.Ф. Шакиров, А.Р. Насыбуллин, Р.В. Фархутдинов, Т.М. Ишкаев
Казанский национальный исследовательский технический университет им. А.Н. Туполева - КАИ, Казань
Аннотация: Работа посвящена анализу частотных характеристик кусочно-однородной коаксиальной брэгговской СВЧ структуры с нерегулярностями в виде ступенчатого изменения радиуса внутреннего проводника. В статье рассматриваются результаты одноволнового и полноволнового моделирования периодических структур с целью анализа влияния параметров структуры на точность расчета. Показана необходимость применения эквивалентных схем нерегулярностей в виде сосредоточенных емкостей в одноволновой математической модели коаксиальной брэгговской СВЧ структуры. При приближении к условиям, благоприятным для распространения первой высшей моды, необходимо включать в расчет эффекты ее распространения.
Ключевые слова: брэгговская СВЧ структура, реактивность, диэлектрическая проницаемость, коаксиальный волновод
Брэгговские СВЧ структуры в коаксиальном волноводе (далее БСВЧСвКВ) представляют собой систему периодически повторяющихся каскадно-соединенных волноводных отрезков с различными волновыми сопротивлениями. Электрические размеры отрезков являются резонансными, то есть соизмеримыми с длиной волны возбуждающего электромагнитного колебания. Такая структура схожа по физическим принципам с брэгговскими зеркалами и фотонными кристаллами в оптическом диапазоне. Практическое применение брэгговских структур не ограничивается использованием в частотно-селективных схемах, антеннах и линиях задержки. В последнее время брэгговские СВЧ структуры активно внедряются в измерительные устройства в качестве их чувствительных элементов [1-3]. Одним из таких приложений служит измерение диэлектрических параметров материалов и веществ. БСВЧСвКВ нашли
применение для измерения диэлектрической проницаемости жидкостей и твердых веществ [4-6]. Сущность метода заключается во введении измеряемого образца во внутреннее пространство волновода, содержащего дополнительно диэлектрические слои с известными характеристиками. Для формирования брэгговского резонанса требуется значительная разница в диэлектрических свойствах сопрягаемых диэлектриков, поэтому в случае близости их диэлектрических проницаемостей необходимо вводить дополнительную нерегулярность в сечениях между диэлектриками. Эта задача может быть решена с помощью ступенчатого изменения внутреннего проводника коаксиального волновода.
Для математического описания частотных характеристик БСВЧСвКВ можно воспользоваться упрощенными методами, основанными на одноволновом приближении, либо полноволновыми методами, использующими численный аппарат решения уравнений электродинамики. Одноволновые методы просты с вычислительной точки зрения и удобны для анализа и синтеза устройств, особенно когда требуется произвести расчет большого количества вариантов. Полноволновые требуют больших вычислительных затрат, но обладают более высокой точностью расчета. В статье ставится задача сравнения двух подходов к решению задачи определения частотной характеристики коэффициента отражения БСВЧСвКВ со ступенчатыми нерегулярностями внутреннего проводника. В качестве одноволнового метода применяется метод трансформации импеданса, а многоволновый реализован в программе электромагнитной симуляции, использующей метод конечных разностей во временной области. В методе трансформации импеданса периодическая структура рассматривается как последовательное соединение отрезков однородных линий передачи (рис. 1), в которых распространяется ТЕМ-волна. Периодически повторяющаяся ячейка состоит из двух волноводных секций,
и
обладающих волновыми сопротивлениями 201 и 2)2, геометрическими длинами Ь1 и Ь2 и комплексными коэффициентами распространения у± и у2.
Рис. 1. - Периодическая структура в виде каскадного соединения кусочно-
однородных отрезков
Коэффициент распространения линии с ТЕМ-волной зависит от диэлектрических параметров заполняющего линию материала и рассчитывается по формуле:
у=а+/р
2л/У17
1+1+
(1)
где/- частота, Гц, с - скорость света, м/с, в' и в" - реальная и мнимая части относительной комплексной диэлектрической проницаемости диэлектрика, заполняющего внутреннее пространство линии.
Волновые сопротивления коаксиальных отрезков рассчитывались, как:
1
¿п =
1-1 п (") ,
1 " а
2ТТу!£'-)£" л/ £О
где - электрическая постоянная, - магнитная постоянная,
(2)
электрическая постоянная, Ь - радиус внешнего проводника, а - радиус внутреннего проводника.
Расчет производится рекуррентно, начиная с конца периодической структуры, то есть отрезка с волновым сопротивлением 201 (рис. 1 ). Для
и
нечетного участка линии передачи (/ - нечетное) входное сопротивление определяется по формуле:
2 вх —2 О 1 ' ( I - 1 ) , (3)
Для четного участка линии передачи (/ - четное) входное сопротивление определяется по формуле:
2 вх —2 О 2~~~Ж~~~), (4)
где /=1 ,2...Ы; N - количество секций; 2Н — 2 (0) — 2 01 - сопротивление нагрузки линии.
Последний отрезок линии передачи с номером N+1 имеет длину Ь1, волновое сопротивление 201 и коэффициент распространения у1. После
вычисления входного сопротивления последнего отрезка , можно
найти коэффициент отражения Г кусочно-однородной линии передачи по формуле
7(л+1) 7
т-<_ -вх -0 1
1 (N+1) . (5)
¿вх 01
Модель исследуемой линии передачи в программе электромагнитного моделирования представлена на рис. 2. Радиус внешнего проводника Ь составляет 5 мм. В качестве проводника выбран материал в виде идеального проводника. Длины секций составляли Ь1=Ь2=12 мм с одинаковым материалом проводника, но разным диэлектриком и разными диаметрами внутреннего проводника, разница между диэлектриками составляла 0.25, 0.5 и 1 мм. В нечетных секциях радиус внутреннего проводника выбирался из ряда 0,75, 1 и 1 ,5 мм и в качестве материала заполнения выступал фторопласт (е'=2,1, е'=0.00042), а в секциях с четным номером радиус внутреннего проводника составлял 0,5 мм, а диэлектрическая проницаемость выбиралась из ряда е'=5, 10, 15. Для каждого варианта получена частотная характеристика коэффициента отражения S 1 1 в диапазоне частот от 0 до 1 0
ГГц в программе электромагнитного моделирования и методом
трансформации импеданса.
а
Рис. 2. - Модель коаксиальной БСВЧСвКВ с неоднородным диэлектрическим заполнением и ступенчатыми нерегулярностями
На рис. 3-5 изображены результаты сравнения частотных характеристик БСВЧСвКВ, полученных методом трансформации импеданса
и рассчитанных в программе электромагнитного моделирования.
2-10® д»10* 4.10® 8«10® 1»101в 0 2хЮ9 4хЮ9 6х109 8х109 1хЮ10
Частота (Гц) Частота (Гц)
а б
Рис. 3. - Коэффициент отражения при в=5 и разницей между радиусами
внутреннего проводника: а) 0,25 мм; б) 1 мм
М Инженерный вестник Дона, №7 (2022) ivdon.ru/ru/magazine/arcliive/n7y2022/7811
а б в
Рис. 4. - Коэффициент отражения при £= 1 0 и разницей между радиусами
внутреннего проводника: а) 0,25 мм; б) 0,5 мм; в) 1 мм
О 2x10® 4хЮ® 6x10® 8x10® 1*Ю10 0 2x10* *хЮв 6хЮ® 8«10® 1хЮ10
Частота (Гц) Частота (Гц)
-Одноволновый расчет
1 - - - • Полноводное моделирование
Г\
О 2x10® 4x10® 6x10® 8x10® 1хЮ10
Частота (Гц)
а б в
Рис. 5. - Коэффициент отражения при £=15 и разницей между радиусами
внутреннего проводника: а) 0, 5 мм; б) 0,5 мм; в) мм
Как можно заметить из рис.3-5, нерегулярность внутреннего проводника значительно влияет на коэффициент отражения. При увеличении диэлектрической проницаемости и разнице между радиусами внутреннего проводника наблюдается частотный сдвиг полосы заграждения и изменение ее ширины при сравнении одноволнового и полноволнового методов. Это говорит о том, что при одноволновом приближении необходимо учитывать реактивности, возникающие в местах изменения внутреннего проводника волновода.
и
В связи с этим было принято решение использовать математическую модель с учетом влияния реактивностей посредством введения в одноволновую модель эквивалентной схемы неоднородности в виде сосредоточенной емкости [7,8].
Эквивалентная емкость может быть описана [9] эмпирическим выражением:
4
Си = 2пЬЕ
1007Г\ а
ч—15
+1.11*10 * (1 - а) * (т - 1)
10^ Ф/м,
(6)
где т = — и а =
Ъ-а? Ь-а
'-, аг - меньший радиус внутреннего проводника, а2
диэлектрическая
больший радиус внутреннего проводника, е -проницаемость в отрезке волновода с меньшим радиусом.
Частотную зависимость емкости [10] можно учесть с помощью эффекта возникновения первой высшей моды ТМ в коаксиальном волноводе:
с(Л =
Са
(7)
/2а2/у£у
На рис. 6-8 отображается результат применения эквивалентной емкости в одноволновой модели БСВЧСвКВ. Как можно заметить включение сосредоточенных емкостей в эквивалентную схему значительно улучает точность расчета характеристик при малых значениях разницы в радиусах ступенчатой нерегулярности.
а
б
М Инженерный вестник Дона, №7 (2022) ivdon.ru/ru/magazine/arcliive/n7y2022/7811
Рис. 6. - Коэффициент отражения при £=5 и разницей между радиусами
внутреннего проводника: а) 0, 5 мм; б) мм
.10® 6-10* Я«Ю® 1«10,с о Зх'о' 4x10* 6.10» 1x10® 1х|0,в 0 2-10® 4«10® 6«10® ««10^ 1x10й
Частота (Гц) Частота (Гц) Частота (Гц)
а б в
Рис. 7. - Коэффициент отражения при £= 1 0 и разницей между радиусами
внутреннего проводника: а) 0,25 мм; б) 0,5 мм; в) 1 мм
Частота (Гц)
а б в
Рис. 8. - Коэффициент отражения при £'=15 и разницей между радиусами
внутреннего проводника: а) 0, 5 мм; б) 0,5 мм; в) мм
Влияние сосредоточенных емкостей, включенных в эквивалентную схему кусочно-однородной коаксиальной брэгговской СВЧ структуры, значительно улучшает частотные характеристики одноволнового расчета, приближая их к результатам более точного полноволнового расчета. При большой разнице в радиусах ступенчатой нерегулярности центрального проводника сказывается близость к условиям распространения высшей ТМ моды. Следовательно, при приближении к условиям, благоприятным для
распространения первой высшей моды, необходимо включать в расчет эффекты ее распространения.
Литература
1. Насыбуллин А.Р., Ишкаев Т.М. СВЧ диэлектрический волновод круглого сечения с брэгговскими периодическими неоднородностями // Инженерный вестник Дона, 2 0 1 7, №2. URL: ivdon.ru/ru/magazine/archive/N2y2017/4232.
2. Насыбуллин А.Р., Ишкаев Т.М., Шаронов Д.Е. Диэлектрических волновод с периодическими неоднородностями как способ измерения характеристик диэлектриков // Инженерный вестник Дона, 20 1 8, №4. URL: ivdon.ru/ru/magazine/archive/n4y2018/5277.7.
3. Усанов Д.А. СВЧ коаксиальный фотонный кристалл с нарушением периодичности в виде диэлектрической структуры // Взаимодействие сверхвысокочастотного, терагерцового и оптического излучения с полупроводниковыми микро- и наноструктурами, метаматериалами и биообъектами. СНИГУ им. Н.Г. Чернышевского. С.: Издательство «Саратовский источник», 2 0 1 9. С. 1 7-20.
4. Фархутдинов Р.В., Насыбуллин А.Р., Морозов О.Г., Вазиев Т.О., Ишкаев Т.М., Садчиков В.В. Брэгговская СВЧ-структура в коаксиальном волноводе как датчик контроля диэлектрических параметров жидких сред // Физика волновых процессов и радиотехнические системы. - 2019. - №2 2 (4-2). - С. 1 1 4-120.
5. Усанов Д.А., Никитов С.А., Скрипаль А.В., Пономарев Д.В., Рузанов О.М. и Тимофеев И.О. Использование сверхвысокочастотной коаксиальной брэгговской структуры для измерения параметров диэлектриков // Радиотехника и электроника. - 2020. - № 5 (65). - С. 495-503.
6. Севастьянов А.А. Брэгговские СВЧ-структуры в коаксиальном кабеле для систем контроля уровня жидких сред // КНИТУ-КАИ, 2015. 136 с.
7. Whinnery J.R., Jamieson H.W., Robbins T.E. Coaxial-Line Discontinuities. Proceedings of the I.R.E. 1944. P. 697-700.
8. Whinnery J.R., Jamieson H.W. Equivalent Circuits for Discontinuities in Transmission Lines. Proceedings of the I.R.E. 1944. P. 98-114.
9. Solmo P.I. The Computation of Coaxial Line Step Capacitatnces. Transactions on microwave theory and texniques. 1967. P. 48-53.
10. Brian C. Wadell. Transmission line design handbook / Brian C. Wadell. Artech House, Inc. 1991. P. 517.
References
1. Nasybullin A.R., Ishkaev T.M. Inzhenernyj vestnik Dona, 2 0 1 7, №2 URL : ivdon.ru/ru/magazine/archive/N2y2017/4232.
2. Nasybullin A.R., Ishkaev T.M., Sharonov D.E. Inzhenernyj vestnik Dona, 2 0 1 8, №4 URL: ivdon.ru/ru/magazine/archive/n4y20 1 8/5277.7.
3. Usanov D.A. Vzaimodejstvie sverhvysokochastotnogo, teragercovogo i opticheskogo izluchenija s poluprovodnikovymi mikro- i nanostrukturami, metamaterialami i bioob'ektami. SNIGU im. N.G. Chernyshevskogo. S.: Izdatel'stvo «Saratovskij istochnik», 20 1 9. pp. 17-20.
4. Farhutdinov R.V., Nasybullin A.R., Morozov O.G., Vaziev T.O., Ishkaev T.M., Sadchikov V.V. Fizika volnovyh processov i radiotehnicheskie sistemy. 2019. №22 (4-2). pp. 114-120.
5. Usanov D.A., Nikitov S.A., Skripal' A.V., Ponomarev D.V., Ruzanov O.M. i Timofeev I.O. 2020. № 5 (65). pp. 495-503.
6. Sevast'janov A.A. Brjeggovskie SVCh-struktury v koaksial'nom kabele dlja sistem kontrolja urovnja zhidkih sred [Bragg microwave structures in a coaxial cable for liquid level control systems]. KNITU-KAI, 2015. 136 p.
7. Whinnery J.R., Jamieson H.W., Robbins T.E. Coaxial-Line Discontinuities. Proceedings of the I.R.E. 1944. pp. 697-700.
8. Whinnery J.R., Jamieson H.W. Equivalent Circuits for Discontinuities in Transmission Lines. Proceedings of the I.R.E. 1944. pp. 98-114.
9. Solmo P.I. The Computation of Coaxial Line Step Capacitatnces. Transactions on microwave theory and texniques. 1967. Pp. 48-53.
10. Brian C. Wadell. Transmission line design handbook. Brian C. Wadell. Artech House, Inc. 1991. pp. 517.