автореферат диссертации по радиотехнике и связи, 05.12.07, диссертация на тему:Исследование резонансных ускоряющих систем с использованием численных методов электродинамического моделирования

П П ЕДЕ 11И E

1. ЛИНЕЙНЫЕ РЕЗОНАНСНЫЕ УСКОРЯЮЩИЕ СТРУКТУРЫ

1.1 пр1шцш1 действия линейных ускорителей попов

1.2 Ускоряющая структура с ПФФ

1.3 Ускоряющая структура с НОКФ

1.4 Радиотехнические характеристики ускоряющих структур

1.5 Этапы проектирования ускоряющих структур ti ныводы но главе

2. МЕТОДЫ ЭЛЕКТРОДИНАМИЧЕСКОГО МОДЕЛИРОВАНИЯ УСКОРЯЮЩИХ СТРУКТУР

2.1 Метод конечных элементов

2.2 метод ПУЛ)

2.2.1 Формулировка метода FDTD

2.2.2 Моделирование резонаторов методом FDTD 39 t* 2.2.3 Разработка программного обеспечения для электродинамического моделирования методом FDTD выводы по главе

3. ОБРАБОТКА ИМПУЛЬСНОЙ РЕАКЦИИ УСКОРЯЮЩЕЙ СТРУКТУРЫ В ЧАСТОТНОЙ ОБЛАСТИ

3.1 Эффект растекания спектра

3.2 ог.зор методов частотного анализа

3.3 пол и гармонический модовый анализ 63 выводы по главе

4. СОБСТВЕННАЯ И НАГРУЖЕННАЯ ДОБРОТНОСТИ УСКОРЯЮЩЕЙ СТРУКТУРЫТ

4.1 Учет потерь в металле методом FDTI)

4.2 Моделирование устройства высокочастотного ввода 82 Выводы но г лаве

5. СТРУКТУРА С ПЕРЕМЕННО-ФАЗОВОЙ ФОКУСИРОВКОЙ

5.1 Исследование короткого макета структуры с ПФФ

5.2 Электродинамическое моделирование иолноразмернои структуры

5.3 Результаты натурного эксперимента 103 Выводы но главе

ЗА КЛ104 ЕНИ Е

Линейные резонаторные ускорители заряженных частиц находят все более широкое применение - как в научных целях, так и для решения прикладных задач. Спектр проблем, решаемых при помощи ускорителей прикладного назначения, чрезвычайно широк. Это лечение раковых опухолей, просвечивание металлических контейнеров на таможне, дефектоскопия, стерилизация медицинского инструмента и многое другое. Качественные характеристики ускорителей при использовании в подобных задачах, как правило, заметно превосходят характеристики альтернативных методов. Линейные ускорители, используемые в промышленных целях, имеют относительно небольшую выходную энергию пучка (от 2 до IOiYIdB), что позволяет делать их небольшими по габаритам (до 4-х метров). Преимущества линейных резонаторных ускорителей делает их использование привлекательным и часто незаменимым.

Требования к технико-экономическим характеристикам ускорителям постоянно растут, поэтому для сохранения конкурентоспособности в рыночных условиях мировой экономики необходимо постоянно совершенствовать высокочастотные системы ускорителей. Необходимо добавить, что производство ускорителей заряженных частиц является весьма наукоемким и требует значительных капиталовложений, поэтому вложение средств в исследования с целыо создания новых, более совершенных ускоряющих структур является вполне оправданным.

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

В рабочем режиме в ускоряющей структуре должен быть возбужден только один — рабочий тип колебаний, остальные типы - паразитные - даже при небольшой интенсивности значительно снижают эффективность ускорения. Ускорение ионов требует мощностей порядка нескольких мегаватт. Необходимо также обеспечить отвод тепла от структуры и компактность источника питания. Поэтому ВЧ-питание ускорителей прикладного назначения осуществляется в импульсном режиме. Паразитные типы колебаний, как правило, имеют частоту, отличающуюся на величину всего лишь порядка нескольких процентов от частоты рабочего типа, поэтому задача оптимизации узла ввода ВЧ-мощности с целыо максимального подавления соседних паразитных типов колебаний является достаточно важной.

При расчете структур применяют различные численные методы. Выбор метода моделирования, а также его качественные характеристики непосредственно влияют на объем и достоверность информации о моделируемой структуре. Наибольшее распространение получил метод конечных элементов, который обладает рядом преимуществ, среди них -возможность задания оптимальной пространственной сетки, сравнительно небольшие затраты вычислительных ресурсов. Поскольку метод конечных элементов имеет дело с собственными функциями, потери в металле и способ возбуждения учитываются в нем лишь косвенно. Среди недостатков также можно отметить - невозможность получить информацию о характере переходного процесса, нестабильность модели с большим количеством узлов.

Метод конечных разностей во временной области (FDTD) дает полную информацию о переходном процессе (результатом FDTD-моделирования являются значения напряженности полей на каждом временном отсчете), равно как и возможность учета особенностей узла связи. Преимущество метода FDTD — его универсальность и возможность эволюционного развития. Метод FDTD позволяет изучить проблемы, связанные с ВЧ питанием ускоряющей структуры и моделированием переходного процесса. Однако, применение метода FDTD, как и любого другого метода моделирования электромагнитных полей во времени, имеет ряд особенностей.

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

Учет потерь в резонаторах из металлов с конечной проводимостью -задача, не имеющая хорошего решения в классической формулировке метода FDTD. На первый взгляд, можно представить реальный металл материалом с потерями, задав соответствующую проводимость. Однако в проводниках с высокой проводимостью наблюдается скин-эффект, иоле проникает в такой проводник очень неглубоко, и для корректного представления поля в проводнике с потерями потребуется нереально малый шаг сетки.

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

Цслыо работы является разработка, практическое воплощение и верификация современных эффективных методов анализа и оптимизации ВЧ-систем резонансных ускорителей ионов.

Для достижения поставленной цели в работе решаются следующие задачи:

1. Построение универсального программного комплекса, реализующего метод FDTD.

2. Моделирование полноразмерной ускоряющей структуры с переменно-фазовой фокусировкой с использованием метода FDTD. Нахождение распределений ускоряющего поля, значений добротности и резонансных частот различных типов колебаний

3. Реализация методики учета потерь в резонаторах из металлов с конечной проводимостью при использовании метода FDTD.

4. Создание электродинамической модели устройства ВЧ-ввода. Исследование влияния способа возбуждения на характер переходного процесса

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

Методы исследовании

Решение поставленных задач моделирования потребовало привлечения математического аппарата, включающего в себя методы численного решения уравнений Максвелла в частных производных и методы частотного-временного анализа. Экспериментальная проверка результатов моделирования осуществлялась в ходе сравнений с результатами, полученными с помощью метода конечных элементов и натурного эксперимента. Научная новизна работы состоит в следующем:

1. Впервые построена имитационная электродинамическая модель полноразмерной ускоряющей структуры с переменно-фазовой фокусировкой с использованием метода FDTD.

2. Разработана методика непосредственного учета омических потерь в реальных металлах в рамках метода FDTD, позволяющая получить полную информацию о переходном процессе с учетом потерь, определять добротность, сдвиг резонансной частоты и другие характеристики.

3. Создана электродинамическая модель устройства ВЧ-ввода. Произведено сравнение качественных характеристик работы ускорителя при различных способах возбуждения. Выработаны рекомендации по оптимизации способа возбуждения.

4. Разработана методика обработки отклика на импульсное возбуждение в резонаторе ускорителя, позволяющая определять амплитуды и фазы полей в ускоряющих зазорах на модах колебаний. Алгоритм позволяет повысить точность определения собственных частот типов колебаний, использование данного алгоритма незаменимо при определении декремента затухания и распределения поле й мод колебаний сложной ускоряющей структуры, у которой резонансные частоты расположены достаточно близко.

На защиту выносится:

1. Методика и алгоритм учета омических потерь в реальных металлах при использовании метода FDTD.

2. Имитационная модель узла ВЧ-ввода мощности, позволяющая моделировать ускоряющие структуры во временной области с учетом особенносей ВЧ-питания.

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

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

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

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

Выводы по главе 5

Проведено моделирование реальных ускоряющих структур с помощью метода FDTD. Полученные значения резонансных частот, добротностей и распределения ускоряющего поля на различных модах соответствуют результатам, полученным в ходе эксперимента, что подтверждает корректность имитационной модели. Таким образом, показана возможность использования метода FDTD для моделирования ускоряющих структур. Моделирование ускоряющих структур методом FDTD, результаты которого представлены в главе 5, осуществлялись автором. Результаты моделирования методом конечных элементов предоставлены Силевым С.А. Натурный эксперимент осуществлялся в НИИЭФА им. Д.В. Ефремова. Материалы главы 5 отражены в публикациях автора [52,58].

ЗАКЛЮЧЕНИЕ

Резонаторные ускоряющие структуры, применяемые в ускорителях прикладного назначения, имеют сложную форму границ. Зачастую это совокупность связанных резонаторов или сравнительно протяженные, нагруженные электродами структуры, в которых, кроме рабочего типа колебаний, представлены многочисленные паразитные типы, в том числе имеющие близкие к рабочему собственные частоты. Ужесточение требований к технико-экономическим характеристикам промышленных ускорителей делает поиск эффективных методов расчета и теоретически и практически важными. Макетирование таких устройств представляется достаточно дорогостоящим и трудоемким, кроме того изменение параметров ускоряющей структуры может производиться лишь в небольших диапазонах, а зачастую такая возможность вообще может быть недоступна. Вследствие этого моделирование ускоряющих структур имеет большое значение.

Для электродинамического моделирования ускоряющих структур построен программный комплекс, реализующий метод конечных разностей во временной области (FDTD). С помощью программного комплекса впервые проведено электродинамическое моделирование полноразмерной структуры с переменно-фазовой фокусировкой (ПФФ) во временной области.

Отличительной особенностью программного комплекса является предложенная методика учета потерь в металлах, основанная на приближенных граничных условиях Щукина-Леонтовича, реализованных в методе FDTD. Методика позволяет непосредственно учесть влияние потерь и получить полную информацию о переходном процессе, определять добротность, сдвиг резонансной частоты и другие характеристики ускоряющей структуры, изготовленной из реального металла.

Предложен и программно реализован метод определения типов колебаний - полигармонический модовый анализ (ПМА), позволяющий определять амплитуды и фазы полей в ускоряющих зазорах на модах колебаний. Метод позволяет повысить точность определения собственных частот типов колебаний, использование данного алгоритма незаменимо при определении декремента затухания и распределения нолей мод колебаний сложной ускоряющей структуры, у которой резонансные частоты расположены достаточно близко. Особенностью предложенного метода является возможность определять распределение поля сразу нескольких мод по продольной координате за один акт моделирования.

Выполнено сопоставление результатов расчетов частот, добротностей и распределения ускоряющего поля на различных модах с результатами натурного эксперимента. Амплитуды полей получены с точностью не хуже чем 4%, собственные частоты - 2%, добротности - 20%. Таким образом, подтверждена корректность имитационной электродинамической модели и показана возможность использования метода FDTD для моделирования ускоряющих структур.

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

1. Каретников Д.В., Шембель, И. П., Сливков Б.К. и др. Линейные ускорители ионов-М.: Госатомиздат, 1962 г. С. 5-10.

2. Линейные ускорители. Теория и эксперимент. Под ред. О. Л. Вальднера. М. Энергоатомиздат. 1991. С. 2-4.

3. Лебедев А.Н., Шальнов А.В. Основы физики и техники ускорителей. т.1 "Ускорители заряженных частиц". М.: Энергоиздат 1981.

4. Лебедев А.Н., Шальнов А.В. Основы физики и техники ускорителей. т.З "Линейные ускорители". М.: Энергоиздат 1983.

5. The Stanford two-mile accelerator. Editor R.B. Neal. New York 1968.

6. Linear accelerators. Pierre M. Lapostolle. Amsterdam 1970.

7. Линейные ускорители ионов, под. ред. Мурина Б.П. т.2. Основные системы. М.: Атом издат 1978.

8. Новиков Д.Л. Линейные ускорители ионов. М.: ЦНИИАтоминформ 1975.

9. Овсянников Д.А., Свистунов Ю.А. Моделирование пучков заряженных частиц в ускорителях. СПб.: НИИ химии СпбГУ, 1998. 60 с.

10. Капчинский И.М. Теория линейных резонансных ускорителей. М.: Энергоиздат. 1982.

11. Линейные ускорители, под ред. Вальднера О.А. М.: МИФИ Энергоатомиздат 1987.

12. Minaev S.A., Svistunov Yu. A., Silaev S. A. Modelling and Testing of APF Cavity RF Field // Workshop BDO-95. p-131-136.

13. Laverne S. Williams. Designing of an RFQ linear accelerator for the use as a mass spectrometer. 1997 // Summer internship in science and technology (SIST). May 19, 1997-August 1, 1997.

14. Ворогушин М.Ф., Малышев B.H. Динамическое согласование в системах ВЧ-питания резонаторных ускорителей // М.: ЦНИИатоминформ, 1984. 12с.

15. Ворогушин М.Ф., Малышев В.Н. Высокочастотное питание резонаторных ускорителей прикладного назначения // М.: Энергоатомиздат, 1989.

16. Милованов О.С., Пятнов Е.Г., Собенин Н.П. Высокочастотные системы линейных ускорителей электронов. // М.: МИФИ 1988.

17. Линейные ускорители. Теория и эксперимент. // Сборник научных трудов. М.: МИФИ 1988.

18. Мальцев И. Г., Тепляков В.А. Система ввода ВЧ-мощности в резонаторы линейного ускорителя протонов УРАЛ-30. // Серпухов, 1976.

19. Техника заряженных ускорителей заряженных частиц. // Сб. науч. трудов, под ред. Спока В.К. М.: 1982.

20. Николаев В.М., Румянцев В.В. Мощные источники СВЧ-питания в линейных ускорителях электронов. Ленинград: НИИЭФА 1979.

21. Линейные ускорители для народного хозяйства. / под ред. Вальднера О.А. М.: Энергоатомиздат 1989.

22. Безногих Ю.Д., Говоров А.И., Зиновьев Л.П. , Мончинский В.А., Семешошкин И.Н. Линейные ускорители ионов на энергию 10 МэВ/нуклон.-Дубна 1983.

23. Mullett L.B., Clay R.E., Madden R.J.B. Mutipactor effect in linear accelerators and other evacuated RF systems, and a new cathode valve // Atomic Energy Research Establishment. Harwell, Berkcsire 1957.

24. Комар Е.Г. Основы ускорительной техники. М. Атомиздат 1975. стр. 368.

25. Будтов А.А., Бахрушин Ю.П., Ворогушин М.Ф. и др. Линейные ВЧ ускорителей прикладного назначения. // Труды двенадцатого совещания но ускорителям заряженных частиц т.2. М. 1990. С. 110-119.

26. Демченко П.А., Ковальчук И.К., Митроченко В.В. и др. Экспериментальные исследования динамики протонов в ускорителях с переменно-фазовой фокусировкой. 1985.

27. Вальднер О.А. Линейные ускорители электронов, (основы расчета ЛУЭ на малые энергии). М. Атомиздат 1965. стр. 186.

28. Ворогушин М.Ф., Давыдов М.С., Смирнов Ю.В. и др. Малогабаритный линейный ускоритель на стоячей волне для дефектоскопии. // Труды третьего совещания по применению ускорителей заряженных частиц. Л. 1979. стр. 239.

29. Малышев В. II., Пикулин Е.В. Характеристики рабочего типа колебаний ускоряющей структуры с переменно-фазовой фокусировкой // Материалы конференции 58 IITK НТО ГЭС. С-Пб. апрель 2003. С. 3-4.

30. Малышев В.II., Пикулин Е.В. Расчет резонансных ускоряющих структур методом FDTD // Тезисы докладов 57 НТК НТО РЭС. С-Пб. 2002. С. 6.

31. Силин Р.А.,. Сазонов В.П. Замедляющие системы «Советское радио» 1966. стр. 253

32. Вычислительные методы в электродинамике / Под редакцией Р. Миттры. Москва, "Мир", 1997.

33. Силаев С. А. Решение трехмерных внутренних задач электродинамики методом конечных элементов // Электронная техника. Сер. Электровакуумные и газоразрядные приборы. 1992. Вып.З. С. 25-30.

34. Норри Д., де Фриз Ж. Введение в метод конечных элементов: Пер. с англ. -М.: Мир, 1981.-304с.

35. Сабоннадьер Ж.-К., Кулон Ж.-Л. Метод конечных элементов и САПР: Пер. с франц. М.: Мир, 1989. - 190 с.

36. Численные методы анализа электрических машин / Под ред. Я.Б. Данилевича. Л.: Наука, 1988.-222 с.

37. Силаев С. А. Разработка методов моделирования электромагнитных полей и устройств и их применение для анализа электровакуумных приборов. С-Петербург. 1997.

38. Silaev S. Л. Isoparametric finite element analysis of time-harmonic electromagnetic fields in three dimensions // Nucl. Instr. and Meth. 1993 Vol. Л328. P. 535-541.

39. Lynch D. R., Paulsen K. D., Origin of Vector Parasites in Numerical Maxwell Solutions // IEEE Trans. Microwave Theory and Tech., vol. MTT-39, March 1991,P. 383-394.

40. Paulsen K. D., Lynch D. R. Elimination of Vector Parasites in Finite Element Maxwell Solutions // IEEE Trans. Microwave Theory and Tech., vol. MTT-39, March 1991,P. 395-404.

41. Силаев С.Л. Исключение ложных видов колебаний при решении трехмерных спектральных задач электродинамики методом конечных элементов // Радиотехника и электроника. 1996. т. 41, Л«4. С. 404-411

42. Yee К. S. Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media // IEEE Transactions of Antennas and Propagation 1966. Vol. AP-14. P. 302-307.

43. Mur G. Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic field equations // IEEE Trans. Electromagn. Compat. 1981. Vol. EMC-22. P. 377-382.

44. Данилин A.A., Малышев В. II. и др. Электродинамическое моделирование методом конечных разностей во временной области (FDTD) СПб.: СПбГЭТУ "ЛЭТИ". 2000.

45. Victor N. Malyshev, Igor A. Lepikhin, Eugene V. Pickulin. Program module for electromagnetic computations using FDTD method // Proceedings 2001 of St. Petersburg IEEE Chapters. SPb ETU "LETI" Publishing Mouse. P. 75-78.

46. Костылева В.В., Лепихин FLA., Малышев В.П., Николаев В.Н. Пикулин Е.В., Смирнов А.Б. FDTD моделирование импульсных излучателей // Материалы конференции 58 НТК НТО РЭС. С-Пб. апрель 2003. С. 5-6.

47. Лепихин И. Л., Малышев В.Н., Пикулин Е. В. Применение метода FDTD для расчета антенно-фидериых устройств // Материалы международной конференции по телекоммуникациям IEEE/ICC2001/St.Petersburg 11-15 июня 2001.С. 242-245.

48. Pursel J.D., Coggans P.M. Л finite-difference time-domain method for solving electromagnetic problems with bandpass-limited sources // IEEE Transactions on antennas and propagation, vol. 47 Л»1, January 1999.

49. Сергиенко Л.Б. Цифровая обработка сигналов. СПб.: Питер, 2002. С 608.

50. Марпл С.А.-мл. Цифровой спектральный анализ и его приложения. М. Мир. 1990г. 584с., с ил.

51. Welch, P.D, "The Use of Fast Fourier Transform for the Estimation of Power Spectra: Л Method Based on Time Averaging Over Short, Modified Periodograms, "IEEE Trans. Audio Electroacoustics, Vol. AU-15 (June 1967), pp. 70-73.

52. Пикулин E.B. Моделирование ВЧ-систем ускорителей ионов // Материалы конференции 59 НТК НТО РЭС. С-Пб. апрель 2004. С. 12-14.

53. Пименов Ю.В., Вольман В.И., Муравцов А.Д. Техническая электродинамика. М.: "Радио и связь" 2000. 526с.

54. Петров Б.М. Электродинамика и распространение радиоволн. М.: Горячая линия Телеком, 2003. 558с.

55. Федоров Н.Н. Основы электродинамики: / Учеб. пособие для вузов. -М.: "Высш. школа", 1980. 399 е., ил.

56. Малышев В.II., Пикулин Е.В. Моделирование резонансных систем с потерями методом FDTD // Тезисы докладов 56 НТК НТО РЭС. С-Пб. 2001. С. 4-5.

57. Eric Jacobsen, Richard Lions. The sliding DFT. // IEEE Signal Processing magazine. March 2003, vol. 20, №2.

58. Маттей Д.Л., Янг JI., Джонс Е.М.Т. Фильтры СВЧ, согласующие цепи и цепи связи, т.2. М.: "Связь", 1972.