Численное моделирование физических процессов в инжекторах отрицательных ионов тема автореферата и диссертации по физике, 01.04.08 ВАК РФ

ВВЕДЕНИЕ.

ГЛАВА. I. ОБЗОР ПРОБЛЕМ, СВЯЗАННЫХ С СОЗДАНИЕМ И ТРАНСПОРТИРОВКОЙ ПУЧКОВ ОТРИЦАТЕЛЬНЫХ ИОНОВ И ИХ МТШАТИЧЕСКИМ МОДЕЛИРОВАНИЙ! .ю

§ 1.1 Требования,предъявляемые к параметрам пучка отрицательных ионов,инжектируемого в термоядерные установки,и способы получения таких пучков

§ 1.2 Обзор численных методов,применяемых при расчетах ионно-оптических систем.

§ 1.3 Обзор методов математического моделирования процессов компенсации собственного.-.заряда ионных пучков.

§ 1.4 Использование перезарядных Сверхзвуковых струй в качестве вакуумных затворов и насосов

ГЛАВА II. МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ДОУСКОРЕНИЯ ПУЧКОВ

ОТРИЦАТЕЛЬНЫХ ИОНОВ.

§ 2.1 Одномерный анализ.

§ 2.2 Моделирование переходного слоя между ускоряющим пространством и компенсированной областью пучка

§ 2.3 Численная модель системы доускорения в двумерной геометрии.

§ 2.4 Результаты численных экспериментов по моделированию системы доускорения

§ 2.5 Доускорение многокомпонентных пучков

ВЫВОДЫ.

ГЛАВА III. МОДЕЛИРОВАНИЕ ПРОЦЕССОВ КОМПЕНСАЦИИ ОБЪЕМНОГО ЗАРЯДА

ПУЧКОВ ОТРИЦАТЕЛЬНЫХ ИОНОВ.

§ 3.1 Особенности компенсации объемного заряда пучков отрицательных ионов и общие вопросы,касающиеся построения моделей пучково-плазменных систем

§ 3.2 Изучение компенсации положительного объемного заряда электронами.

§ 3,3 Моделирование компенсации объемного заряда пучка отрицательных ионов

§ 3.4 Анализ устойчивости пучково-плазменной системы,образующейся при компенсации объемного заряда пучка отрицательных ионов

ГЛАВА 1У. ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ПРОНИКНОВЕНИЯ ГАЗА ЧЕРЕЗ ;

СВЕРХЗВУКОВУЮ ПАРОМЕТАЛЛИЧЕСКУЮ СТРУЮ ПРИ ИЗУЧЕНИИ ЗАТВОРНЫХ СВОЙСТВ ПЕРЕЗАРЯДНОЙ СТРУИ В ИНЖЕКТОРАХ.

§ 4.1 Постановка задачи и методы решения.

§ 4.2 Описание алгоритма расчета и выбор вычислительных параметров

§ 4.3 Обсуждение результатов расчетов

ВЫВОДЫ.

I. Необходимость использования отрицательных ионов для инфекционного нагрева плазмы в перспективных схемах управляемого термоядерного синтеза (УТС).

В настоящее время в исследованиях по УТС с магнитным удержанием плазмы достигнуты большие успехи. Дальнейший прогресс в этой области связывается с разработкой эффективных методов нагрева плазмы до термоядерных температур. Последние экспериментальные работы по ин-жекционному нагреву плазма в магнитных ловушках [I] подтвердили перспективность нагрева с помощью инжекции в ловушку высокоэнергетичных нейтральных атомов. Сейчас созданы и успешно развиваются инжекторы нейтральных атомов, основанные на перезарядке интенсивных пучков положительных ионов. Эффективность перезарядки в энергетическом диапазоне 10ч-50 кэв составляет при этом 90-50/&. Энергия нейтралов в десятки кэв достаточна для инжекции в Т# установки, предназначенные для демонстрации осуществимости управляемой термоядерной реакции с положительным энергетическим выходом. В проектах термоядерных реакторов промышленного типа требуемые энергии нейтралов составляют сотни кэв [2] . Реализация ряда перспективных схем УТС также связана с использованием высокоэнергетичных нейтральных атомов. В частности,в амбипо-лярных открытых ловушках Сз] требуемая энергия нейтралов, инжектируемых в концевые ловушки, должна составлять около I Мэв. Термоядерные з реакторы будущего, основанные на реакциях I) +]) и Р+Не» также потребуют разработки эффективных инжекторов с энергиями нейтралов в сотни кэв.

С ростом энергии выход нейтралов при перезарядке положительных ионов в газовой обдирочной мишени быстро падает (кривая I на рис.1), тогда как выход нейтралов при обдирке отрицательных ионов остаётся дейтерия в зависимости от энергии,приходящейся на один дейтрон.

Рис,2 Зависимость энергетического к.п.д. получения нейтралов от энергии ионов. практически постоянным до Еесьма больших энергий (кривая 2 рис.1). При перезарядке на газовой мишени эффективность преобразования отрицательных ионов в нейтралы составляет 60%, а на плазменных мишенях эффективность может достигать 80% [4] .

Оценки энергетических затрат на получение одного отрицательного иона путем перезарядки положительных ионов показывают, что энергетическая цена отрицательного иона, полученного при перезарядке на цезии, составляет около 30 кэв, а в натрии - около 75 кэв [5] . Эти величины определяют область энергий, при которых получение нейтралов через отрицательные ионы становится энергетически выгодным. На рис.2 представлена зависимость энергетического к.п.д. инжектора от энергии ионов для разных способов получения нейтралов,полученная в работе[б]. Энергетический к.п.д. получения нейтралов непосредственно из положительных ионов падает с энергией (кривая 3), даже с учетом рекуперации энергии неперезарядившихся ионов (кривая 4). При энергиях свыше 200 кэв использование отрицательных ионов имеет явное преимущество, а при энергиях свыше 300 кэв становится единственно приемлимым. Поэтому разработка и создание интенсивных потоков высокоэнергетичных отрицательных ионов является одной из ключевых проблем, стоящих на пути реализации некоторых программ УТС.

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

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

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

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

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

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

Окончательные результаты используются при планировании реального физического эксперимента и сравнении с его результатами.

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

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

Это позволяет проводить экспериментальные исследования по отработке инжектора более эффективно с точки зрения использования имеющихся ресурсов.

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

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

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

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

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

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

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

В главе 1У приводятся результаты моделирования методом статистических испытаний / метод Монте-Карло/ проникновения газа через сверхзвуковую парометаллическую струю. Проведено исследование затворных свойств таких струй в широком диапазоне параметров исследуемых пар "газ-па]:

Результаты исследования,изложенные в главе 1У, получены в соавторстве с А.И.Крыловым а Б.В.Кузнецовым.

Постояноое внимание к работе со стороны научного руководителя Н.Н.Семашко в большой мере способствовало её продвижению.

Заключение.

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

В результате проведенной работы был создан комплекс программ, обеспечивающих численное моделирование следующих процессов;

- доускорение интенсивных пучков отрицательных ионов до требуемой энергии ( гл.П);

- компенсация пространственного заряда интенсивных пучков отрицательных ионов (гл.Ш);

- газокинетическая очистка пучка от сопутствующего газа с помощью сверхзвуковых парометаллических струй - мишеней (гл.1У).

В соответствии с поставленными задачами были разработаны следующие математические модели:

1, Аналитическая модель системы доускорения, основанная на анализе одномерного доускоряемого потока с последующим нахождением распределения потенциала в окружающем пучок пространстве путем решения задачи Коши для уравнения Лапласа. Такая модель является исходной при синтезе системы доускорения (раздел 2.1).

2. Одномерная модель переходного слоя между компенсированной областью пучка и пространством доускорения, полученная на основе численного решения уравнения Власова для положительных ионов и уравнения Пуассона (раздел 2.2.)

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

4. Численные модели процесса компенсации положительного пространственного заряда электронами на основе одномерного уравнения Власова для электронов и уравнения Пуассона и на основе метода "крупных частиц" (раздел 3.2.)

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

6. Статистическая модель диффузии сопутствующего пучку, газа через сверхзвуковую перезарядную мишень на основе метода "Монте-Карло" (раздел 4.1).

Разработанные средства математического моделирования позволили провести детальное излучение упомянутых процессов путем проведения большого объема численных экспериментов на ЭВМ БЭСМ-6 при различных значениях физических параметров.

1. £и$сигк H. ei "PLT mochtet? ¿ea*rv Лесы ъьби ¿/г- " а, PA*, 4öC6 Ccsit/ Gpfii //uc £east /^¿c Ъ*? * f P<u> с. 7/A, Тгюбёъиск, /978) v. J, JA EÄ, Vc'tinnc^, /919, f. /67.2. ^ 2). L. ^üyn- (¿ttiTesb -¿окси^аж.1.ion. ') /Vuc~£ FuiiorL. f и /v/^, /9??:

2. Димов Г.И.,Закайдаков В.В.,Кишеневский M.E. Термоядерная ловушка с двойными пробками. "Физика плазмы",1976,т.2,в.4.

3. Димов Г.И. и др. "Ядерный синтез",1975,т.15,с.551.5# Семашко H.H.,Владимиров А.Н.,Кузнецов В.В.,Кулыгин В.М., . Панасенков A.A. Инжекторы быстрых атомов водорода, М.,"Энергоиздат",1981.

4. Крылов А.И.»Кузнецов В.В.,Семашко H.H. Физические аспекты инжек-ционной системы для больших токамаков и открытых ловушек. -"Атомная энергия",1980,т.48,в.З.7."Атомные и молещглярные процессы", под ред.Бейтса, М.,ИЛ,1964.

5. Доусон Дж. Физика пучков заряженных частиц.,М.,"Мир", 1980.

6. Димов Г.И. Перезарядный метод инжекции протонов в ускорители и накопители. Препринт ИЯФ СО АН СССР ЖЮ4,Новосибирск, 1969.

7. Семашко H.H. Исследование и создание инжекторов быстрых ионови атомов водорода для стационарных магнитных ловушек. Докторская диссертация, М.,1974.

8. I^iyrxo-tf Q Т^ ; " Ръос. ¿^Kv^b . Ptoo(u.c/Lio ГЬ о^ AN£QCL-iiTr4 Hjcl*. Tont Л/7с/ /'

9. Косслет В. Введение в электронную оптику. М.,1950.- из

10. Кельман В.М.,Явор С.Я. Электронная оптика, М.-Л.,изд.АН СССР,1963

11. Стэррок П. Статистическая и динамическая электронная оптика., М., изд."ИЛ",1965.

12. Тихонов А.Н.Дрсенин В.Я. Методы решения некорректных задач.,М., "Наука", 1974.

13. Лаврентьев М.М. О задаче Коши для уравнения Лапласа. -"Известия АН СССР",Серия "Математика",т.20,1956.

14. Pte>*cce У-Q. „ Theot^ Оле/ De-s^n, o/S/&18. tt тАе Гбеоъу Рсвтсе

15. Gun. '. ^О^Ъгъа £ &*?<>/ /btr/ta^, к 4л/%/$$<19. сА>г&/, г£ зг, ^ ¿я^ /p/s.20. ¿а^тссс^ , тГ. 2Z, f>. 34?; /9^3.

16. Волков Б.И.,Свешников А.Г.,Семашко H.H. К задаче определения формы плазменного эииттера в ускоряющем электрическом поле.-"ДАН СССР",1971,т.201,в.4.

17. Болдасов B.C.»Волков Б.И.,Свешников А.Г.,Семашко H.H. К определен: формы плазменного эмиттера и ускоренного пучка в ионно-оптической системе. "ДАН СССР",1974,т.218,в.5.

18. Волков Б.И. Математическое моделирование движения заряженных частиц в плазменном инжекторе. Сб."Численные методы в физике плазмы", М.,"Наука",1977.

19. Волков Б.И.,Свешников А.Г.,Семашко H.H. Стационарное движение пучка заряженных частиц с учетом собственного пространственного заряда. "ДАН СССР",1969,т.189,в.5.

20. Волков Б.И. КЕМ и МФ,1969,т.9,М,с.961.26. Büßerncl*V 0.}

21. PAgt. гг //S, лз, f. ¿03, /2SIl

22. Масленников M.B.,Сигов Ю.С. ДАН СССР,1964,т.159,в.5,с.1013.- 114

23. Вычислительные методы в физике плазмы. М.,"Мир", 1974.

24. Анастасевич B.C. Теория компенсации ионных пучков,ЖТФ,1956,т.26,в.7 3D. ¿¿/го4ъ E.G., Hebru^VLb4 KG, ft Syoace сЛа,уе ¿ffiic*оъ Secant агьс/ 'си/ис/сЪ/п. ^ ¿ootc•¿¿ve ¿en> ¿zcyeptstp 1 yousuv. o^ 4/bfitf^^

25. X/}. ft of?-/Ал ca&c&'o^ejs jb&AsriCL -¿/¿¿ьт^С TCO^I- " Of? ¿Zctcb/f^1. Я 6, /963.32. Qusbsv&A., о/ о/е**'^

26. One/ o/c'di(tt'&c/t'oib cl eecLSbi-^esitstA'/e-ajD&Li^c^. ^ourub. о/ , v. /264.

27. Габович М.Д. и др. Газовая фокусировка пучка отрицательных ионов и её особенности. ЖТФ,1974,Hi.44,в.4.34. МсСоък.£^>1. Ci'dc.'C. e/ZecJ^bort" "f

28. P^atmA. PAgte'ct f v. /3; /#?</.

29. Габович М.Д. и др. Самодекомпенсация стабильного квазинейтрального ионного пучка за счет кулоновских соударений. "Физика плазмы",1975,т.I,в.2.

30. Me-inpffosi „¿fee, се сЛ^и^е lh он с/ -¿Аас'ъ ^ы^/и c/t'oe*.jDOic'-ltire сое* •¿bOytybtv^ J ¿fowbsu. OJ?4f>/>£

31. Габович М.Д.,Горецкий В.П.,Джаббаров Д.Г.,Козырев Ю.Н.,Найда А.П. Коллективные процессы в плотных пучках отрицательных ионов и декомпенсация этих пучков. Препринт $9,ИФ АН УССР,Киев, 1979.

32. Бредихин М.Ю.»Ильченко А.М.,Скибенко Е.Н.,Юферов В.Б. Исследование сверхзвуковых газовых струй в качестве вакуумного затвора. ЖТФ,1969,т.39,в.II.

33. Иванов В.П. Огазокинетической очистке ионов от нейтральных нолекул, ЖТФ,1974,т.44,в.2.

34. Дьячков Б.А.,Зиненко В.И. К вопросу использования вакуумных свойств газодинамической струи паров щелочных металлов. ЖТФ, 1980,т.50,в.II.

35. Крылов А.И.,Кузнецов В.В. Сверхзвуковая паровая струя как вакуумный затвор. Препринт ИАЭ-3330/7, М.,1980.

36. Дьячков Б.А.,Крылов А.И.»Кузнецов В.В., Семашко H.H. Натриевая мишень инжектора отрицательных ионов. "Атомная энергия",1980, т.49,в.4.

37. Самарский A.A.,Андреев В.Б. Разностные методы для эллиптических уравнений. М.,"Наука",1976.

38. Захаров А.В.,Самарский A.A.»Свешников А.Г. Расчет движения заряженного пучка методом больших частиц с учетом собственногопространственного заряда. ДАН СССР,т.197,№3,1971.

39. Дьячков Б.А.,Крылов А.И.,Кузнецов В.В.,Семашко H.H. Получение отрицательных ионов путем перезарядки водорода на натриевой мишени.

40. Препринт ИАЭ г 2523, М.,1975. , Л , ЛЛ ,

41. Габович М.Д.»Джаббаров Д.Г.,Найда А.П. Эффект декомпенсации плотного пучка отрицательных ионов. "Письма в ЖЭТФ",т.29,в.9, 1979.

42. Телегин В.И. Об одной разностной схеме для уравнения Власова ЖВМиМФ,т.16,165,1376.

43. Кулыгин В.М. »Телегин В.И. Компенсация электронами пространственно-ограниченного положительного заряда. "Атомная энергия",т.41, в.4,1976.

44. Свешников А.Г.,Телегин В.И. О численном решении уравнения Власова для ионов в пространственно-ограниченной области", Сб."Вычислительные методы и программирование",Изд-во МГУ,в.28,1978.- 116

45. Нулыгин В*М., Соколов В.Г. "О функции:распределения электронов) компенсирующих положительный объемный заряд".Доклад на III Всесоюзном семинаре по физике и технике интенсивных источников ионов и ионных пучков,Киев,1980.

46. Михайловский A.B. "Теория плазменных неустойчивостей. т.1, М., "Атомиздат",1970.

47. Корн Г.,Корн Н. Справочник по математике.,М.,"Наука",1973.

48. Бусленко Н.П.,Голенко Д.И.,Соболь И.М.,Срагович В.Г.,Шрейдер Ю.А. Метод статистических испытаний Сметод монте-карло).,Физматгиз, М-., 1962.

49. Хастед Дж. Физика атомных столкновений. М.,"Мир",1965.

50. Габович М.Д. и др. Возбуждение колебаний плазмы быстрым пучком отрицательных ионов. ЖЭТФ,1974,т.67,в.5(11)