Влияние транспорта водорода в первой стенке термоядерных реакторов на рециклинг

Маренков, Евгений Дмитриевич

Влияние транспорта водорода в первой стенке термоядерных реакторов на рециклинг тема автореферата и диссертации по физике, 01.04.08 ВАК РФ

Маренков, Евгений Дмитриевич АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Москва МЕСТО ЗАЩИТЫ

2013 ГОД ЗАЩИТЫ

01.04.08 КОД ВАК РФ

Диссертация по физике на тему «Влияние транспорта водорода в первой стенке термоядерных реакторов на рециклинг»

Автореферат диссертации на тему "Влияние транспорта водорода в первой стенке термоядерных реакторов на рециклинг"

На правах рукописи

Маренков Евгений Дмитриевич

Влияние транспорта водорода в первой стенке термоядерных реакторов на рециклинг

Специальность 01.04.08 — Физика плазмы

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Автор: /і/і^Р 2 8 МАР 2013

и «

Москва 2013

005051007

Работа выполнена в Национальной исследовательском ядерном университете «МИФИ»

Научный руководитель:

Официальные оппоненты:

Ведущая организация:

кандидат физико-математических наук, доцент кафедры «Физика плазмы» НИЯУ МИФИ Цветков Игорь Владимирович Хвесюк Владимир Иванович, доктор технических наук, профессор, заведующий кафедрой теплофизики МГТУ им. Н.Э. Баумана Цвентух Михаил Михайлович, кандидат физико-математических наук, научный сотрудник Отдела физической электроники Физического института им. П.Н.Лебедева РАН (ФИАН)

Национальный исследовательский центр «Курчатовский институт», Москва

Зашита состоится 17 апреля 2013 г. в 16:30 на заседании диссертационного совета Д 212.130.05 в НИЯУ МИФИ по адресу: Москва, 115409, Каширское шоссе, 31, корпус «К», ауд. К-608, тел: 324-84-96, 323-95-26. С диссертацией можно ознакомится в библиотеке НИЯУ МИФИ.

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

Автореферат разослан « 16 » марта 2013 года.

Ученый секретарь диссертационного совета

Р.С. Стариков

Общая характеристика работы

Актуальность работы

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

Можно выделить следующие три проблемы, играющие ключевую роль во взаимодействии плазмы с материалом первой стенки:

1. время жизни обращенных к плазме материалов;

2. образование пыли;

3. накопление трития в материалах вакуумной камеры.

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

в камере трития должно быть настолько минимальным, насколько это возможно, с тем, чтобы избежать загрязнения окружающей среды и эвакуации местного населения из-за его случайного высвобождения. В соответствии с этим, для международного экспериментального термоядерного реактора ИТЭР (ITER —International Thermonuclear Experimental Reactor) установлены следующие «административные» пределы, которые гарантируют, что предельные с точки зрения безопасности значения не будут достигнуты даже с учетом погрешности измерений: 670 кг пыли и 700 г накопленного в камере трития.

Понимание важности процессов взаимодействия плазмы с первой стенкой токамаков стимулировало большое количество как экспериментальных, так и теоретических работ в этой области. Разработаны компьютерные коды, позволяющие моделировать процессы переноса в пристеночной плазме (UEDGE, ERO, SOLPS и т.п.); достигнут существенный прогресс в понимании нестационарных (off-normal) процессов в токамаке, таких как срывы, MARFE, ЭЛМы (ELM), возникновение «горячих точек» (hot spots) и т.д. Накопление и транспорт изотопов водорода в материалах, которые планируется использовать в первой стенке (бериллий, вольфрам и СС-композит), также являются объектом интенсивного исследования. Анализ экспериментальных данных показывает, что транспорт водорода в чистых материалах определяется процессами диффузии, захвата/выхода из дефектов различных типов, накоплением в порах и поверхностными процессами. В условиях первой стенки, однако, большую роль играет также эрозия и перенапыление эродированных частиц, что может существенно изменить свойства материала.

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

яния могут служить резкая десорбция большого количества водорода во время MARFE и срывов, ЭЛМов, появление на стенке «горячих точек» (hot spots) и т.п. В тоже время, большинство экспериментальных работ по изучению взаимодействия плазмы со стенкой делается в лабораторных условиях не на токамаках, а на модельных установках, в которых материал подвергается воздействию плазмы или облучению ионным пучком с параметрами, близкими к ожидаемым в термоядерных реакторах. Результаты, полученные в таких условиях, не могут учесть взаимного влияния процессов, происходящих в пристеночной плазме токамаков и процессов, отвечающих за транспорт водорода на стенке. В большинстве теоретических работ такое взаимное влияние также не учитывается, т.к., из-за своей сложности, даже расчет этих процессов независимо друг от друга наталкивается на существенные трудности. Например, известно, что хорошее удержание плазмы в Н-режиме связано с установлением полоидального профиля плотности плазмы с резким градиентом на периферии («пьедестал плотности»). Во время вспышки ЭЛМа этот пьедестал разрушается, а после — восстанавливается. На настоящий момент механизм этого процесса не ясен, хотя понятно, что процессы рециклинга могут оказывать влияние на изменение плотности периферийной плазмы и, следовательно, на восстановление пьедестала.

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

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

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

Несмотря на большое количество различных процессов, участвующих в транспорте водорода в твердом теле, на практике транспорт чаще всего описывается с помощью решения уравнений переноса (уравнение диффузии для водорода в поле ловушек). Существует несколько кодов для решения этих уравнений, из которых код ТМАР является одним из наиболее часто используемых. Вместе с тем, в коде ТМАР 7 может рассматриваться не более трех различных типов ловушек, отличающихся энергиями связи с ними водорода. Однако, некоторые эксперименты могут быть удовлетворительно объяснены только при учете значительно большего количества типов ловушек (например, термодесорбционные спектры водорода из вольфрама с примесью карбида титана; вольфрама, подвергшегося нейтронному облучению и т.п.). В связи с этим разумно предположить, что в материалах со сложной структурой (например, перенапыленные слои на первой стенке) более адекватной моделью транспорта является введе-

ние непрерывного распределения ловушек по энергиям связи с ними водорода.

В соответствии с указанными проблемами была сформулирована цель данной работы.

Цель

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

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

2. Исследовать роль первой стенки в восстановлении пьедестала после вспышки ЭЛМа.

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

4. Разработать модель транспорта изотопов водорода в материалах с широким распределением ловушек по энергии выхода.

На защиту выносятся следующие результаты, содержащие научную новизну:

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

2. Нульмерная модель реакции первой стенки токамака на вспышки ЭЛМа и условия, в которых газовыделение из первой стенки токамака может играть роль в восстановлении пьедестала после вспышки.

3. Компьютерный код, основанный на разработанной модели транспорта изотопов водорода в материалах с учетом теплопроводности, транспорта и возникновения дефектов, взаимодействия водорода с дефектами и междоузельными атомами, диффузии растворенного водорода. Результаты компьютерного моделирования реакции стенки на вспышку ЭЛМа, оценка влияния бе-риллиевой стенки на восстановление пьедестала Н-моды.

4. Модель проницаемости водорода в многослойных материалах в стационарном режиме с учетом поверхностных процессов.

5. Анализ качественных изменений в газовыделении из материала, вызванных наличием дефектов с непрерывным распределением по энергиям выхода водорода.

Научная и практическая значимость работы

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

Апробация работы

Основные результаты работы докладывались на следующих конференциях:

1. 16 конференция «Взаимодействие плазмы с поверхностью», Москва, НИЯУ МИФИ, 2013

2. 13th International Workshop on Plasma-Facing Materials and Components for Fusion Applications and 1st International Conference on Fusion Energy Materials Science (PFMC), Rosenheim, Germany, 2011

3. 20th international conference PSI-2012. Aachen, Germany, 2012

4. 15 конференция «Взаимодействие плазмы с поверхностью». Москва, НИЯУ МИФИ, 2012

5. 20 international conference Ion-Surface Interactions (ISI)-2011, Zvenigorod, Russia, 2011

6. 19 international conference IonSurface Interactions ISI-2009, Zvenigorod, Russia

7. 18 international conference IonSurface Interactions ISI-2007, Zvenigorod, Russia

8. International conference of atomic collisions on surfaces ICACS-24, Krakow, Poland, 2012

9. Четвертая международная конференция «Взаимодействие изотопов водорода с конструкционными материалами. IHISM-2008»

10. Третья международной конференция «Взаимодействие изотопов водорода с конструкционными материалами. IHISM-2007»

11. Научная сессия МИФИ-2007, Москва, НИЯУ МИФИ

12. Учебно-научная конференция-конкурс по физике УНК-2007, Москва, ФИАН

По результатам исследований опубликовано 19 печатных работ, из них 5 — в рецензируемых журналах, рекомендованных ВАК (см. список публикаций).

Объем и структура диссертации

Диссертация состоит из введения, пяти глав и заключения. Полный объем диссертации составляет 114 страниц с 45 рисунками и 3 таблицами. Список литературы содержит 116 наименований.

Содержание работы

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

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

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

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

мущений были получены дисперсионные уравнения, а из них — условия развития неустойчивости. В случае десорбции первого порядка, когда предполагается, что выход водорода происходит из ловушек одного типа, характеризующихся энергией выхода Еи этот критерий имеет вид:

+ (1) 1 - /г \ТЕ Т\у /

Здесь те — характерное время удержания энергии в плазме, тр — характерное время удержания частиц в плазме, а тцг = N/Т есть отношению количества водорода в стенке N к потоку водорода из стенки Г в равновесии и представляет собой, поэтому, время удержания водорода в стенке. Особую роль играет параметр

Е1 То — Тг

£ = 2—--=-£г

Т0 То

характеризующий, в конечном итоге, реакцию стенки. Здесь То, Тг — температуры обращенной к плазме и внешней поверхностей стенки в стационарном режиме. Величина £г = п2И/Н определяет отношение стационарных потерь на излучение и перезарядку п2Н (п — плотность плазмы, Я — коэффициент пропорциональности), к эффективной мощности нагрева Н. Параметр 1 — Н равен отношению потока из стенки только за счет десорбции, N ехр (—Е^/То) /т4, к потоку из плазмы на стенку, Г в равновесии, где т{ — характерное время выхода из ловушек.

Таким образом, возможность возникновения неустойчивости определяется безразмерными параметрами £, К и временами удержания энергии и частиц. При этом характерные значения е ~ Юч-ЗО , времена те,р обычно порядка 1 с, а К и т\у могут меняться в очень широких пределах из-за их экспоненциальной зависимости от температуры. Из критерия (1) видно, что неустойчивость вероятна для больших е , что соответствует, прежде всего, большой доле потерь £г (температура и

энергия входят также в правую часть неравенства и поэтому их влияние сложнее).

На практике выполнение условия (1) возможно, только если температура стенки остается в определенном диапазоне, зависящем от параметров взаимодействия плазмы со стенкой. При типичных для токамаков параметров (поток на стенку порядка 1015 см~2с-1, времена удержания энергии и частиц тЕ ~ тр ~ 1 с, энергия десорбции Et ~ 1 эВ) этот диапазон лежит в районе 600 К. Вместе с тем, на практике выход водорода часто определяется не одним, а несколькими типами ловушек, отличающихся энергиям связи Ег, Так как каждая из этих энергий дает свой температурный диапазон развития неустойчивости, реальные температуры стенки, при которых возникает неустойчивость, могут лежать в гораздо более широких пределах.

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

Аналогичный анализ был приведен для случая, когда десорбция водорода из стенки определяется диффузией. Как и в предыдущем приближении, возникновение неустойчивости определяется параметром е, в котором вместо энергии связи с ловушками Et стоит энергия активации диффузии Еа. Кроме того, большую роль играет величина 5 = ттГ/птахХг,пр, где тт — характерное время теплопроводности, Г — поток водорода из плазмы на стенку, хгтр — глубина внедрения водорода в стенку, а птах — максимальная плотность водорода, достижимая в материале стенки. Развитие неустойчивости маловероятно для значений 5 » 1. Для случая ё «С 1 критерий развития неустойчивости может быть записан в довольно простой форме:

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

£) = 3.5 х Ю-3 ехр(—0.39 эВ/Т) см2/с

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

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

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

. = ^ехр(-Е^Т0(Ь)) +Т/Мт. (2)

где Г — поток частиц из плазмы на стенку (зависящий от времени), 7\Гт — максимальное накопление водорода в стенке (на единицу площади), V ~ 1013 с-1 — характерная частота колебаний атомов. Если изменение потока частиц из плазмы на стенку, вызванное вспышкой ЭЛМа, мало,

|<*1пГ(4)/сЙ| ^ '

поток водорода из стенки в плазму просто равен потоку из плазмы на стенку, и стенка не может играть роли в восстановлении пьедестала. В противоположенном случае

|dlnr(í)/df| <<; (3)

влияния стенки на восстановление пьедестала можно ожидать только

когда термодесорбция водорода достаточно эффективна, то есть

K <4>

(для Et = 1.5 эВ, Nm = 1018 см"2, v = 1Ü13 с"1 и Г = 1017 см"2^1 ). Эти условия приводят к тому, что существенное влияние десорбции водорода возможно лишь для определенного температурного диапазона, который для указанных параметров лежит в районе 500 - 800 К.

Эти качественные выводы были проверены численным моделированием реакции бериллиевой стенки в коде FACE (автор Р.Д. Смирнов, UCSD), специально модифицированном для решения этой задачи. Взаимодействие водорода со стенкой в этом коде описывается одномерной моделью, учитывающей транспорт растворенного водорода, вакансий и междоузельных атомов бериллия, созданных в результате воздействия на стенку потока плазмы, а также возможные реакции между этими компонентами. Температура стенки вычисляется из решения одномерного уравнения теплопроводности.

На рис. 1 показаны зависимости потока газовыделения из стенки Гout от входного потока T¿n от времени для трех различных температур стенки, полученные моделированием (показан один из импульсов). Для каждой температуры рассмотрены два режима: когда стационарная температура определяется излучением (rad) или потоком частиц (flux). Из рисунка видно, что при выбранных параметрах взаимодействия водорода со стенкой можно ожидать существенного

t, ms

Рис. 1. Зависимости потока газовыделения из стенки rout и входного потока Tin от времени для трех различных температур стенки в случаях, когда температура контролируется излучением (rad) или потоком частиц (flux). Оба потока нормировании на максимальную величину потока на стенку

Гятгах-

(а) (Ъ)

Рис. 2. Распределение по глубине концентраций растворенного водорода Н, вакансий V, водорода, захваченного в вакансии УН, междоузельных атомов бериллия I, и междоузельных атомов бериллия, связанных с водородом 1Н для То = 400 К (а) и То = 800 К (б), температура контролируется излучением. Сплошные маркеры соответствуют началу вспышки, пустые —

окончанию.

влияния десорбции газа на восстановление пьедестала, только если значение этой температуры мало (около 400 К). При температуре порядка 800 К поток десорбции газа из стенки практически совпадает с потоком водорода из плазмы на стенку, т.е. система находится в квазистационарном режиме, и реакция стенки на вспышку ЭЛМа не может оказывать существенного влияния на восстановление пьедестала.

Анализ изменения профилей распределения частиц по глубине показывает, что газовыделение определяется, в основном, диффузией растворенного водорода (рис. 2, где показано изменение профилей для двух температур 400 К и 800 К; поток на стенку равен 2,64 х Ю20 см~2с-1, мощность излучения 2 МВт для 400 К и 6 МВт для 800 К).

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

В главе «Проницаемость водорода через многослойные мембраны» получены некоторые общие закономерности проницаемости водорода через многослойные материалы из газа.

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

1 - — —

V" К + К (5)

где и \г2 — параметры проницаемости входной и выходной сторон,

равные:

г=1

(7)

Здесь Р — давление газа. Di, Si, h — коэффициенты диффузии, растворимости и толщины слоя г, К,,л и К„п — коэффициенты абсорбции на входной и выходной сторонах. Для случая п = 1 параметры V\ и Vn сводятся к известным параметрам проницаемости через несимметричную однослойную мембрану.

Если параметр 1, проницаемость водорода ограничена диф-

фузией в слоях мембраны, а поверхностные процессы не играют роли (DLR режим). Если же V <С 1, поверхностные процессы доминируют (SLR режим).

На практике большое значение имеет зависимость проникающего потока через мембрану от температуры. Анализ уравнений проницаемости показывает, что эта зависимость во многом определяется соотношением между «энергиями активации транспорта» Ер = Ed + Qs в различных слоях, где Ed — энергия активации диффузии, a Qs — теплота растворимости. Выделяются следующие случаи:

1. если все Ер > 0, то проникающий поток Jp является возрастающей функцией температуры;

2. если все Ер < 0, то JP(T) имеет максимум при температуре, когда V = 1;

3. если некоторые из Ер положительны, а некоторые отрицательны и выполняется условие

Рис. 3. Зависимости проникающего потока через двухслойную мембрану от температуры для различных энергий активации транспорта Ер. Квадратные маркеры: Ер\ = 0.1 эВ, Ер2 = 0.7 эВ. Треугольные маркеры: Ер\ = 0.3 эВ, ЕР2 = -0.4 эВ. Круглые маркеры: Ер1 = —0.1 эВ, Ер2 — —0.7 эВ.

то JP(T) возрастает. Если же условие (8) не выполнено, JP(T) имеет максимум.

Эти результаты проиллюстрированы на рис. 3, где показаны возможные зависимость JP{T) для двухслойной мембраны для различных энергий активации транспорта Ер. Использованные значения Ер указаны в подписи к рисунку. Величины li/DoiSoi = 104 см'"2, ¿2/A)2<So2 = 1 см-2. (Такие параметры выбраны, чтобы сделать иллюстрацию ясной; вместе с тем, они разумны по порядку величины для реальных материалов).

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

зависит от свойств поверхности).

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

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

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

менах плотность водорода становится много меньше плотности ловушек, можно использовать для анализа транспорта водорода модель случайных блужданий по решетке с распределением времен ожидания прыжка из узла в узел т, Рт (т). Характерное время выхода из ловушки, связанное с температурными флуктуациями, равно Те = то ещ>(Е/Т), где То — характерное время. Поэтому распределение по временам прыжков связано с распределением ловушек по энергиям связи с ними водорода Е, Ре(Е/Т). В наиболее общем виде эту связь можно записать как

здесь / — функция, удовлетворяющая условиям нормировки и определяющая детали связи Рт(т) с Ре(Е/Т). Нетрудно видеть, что конкретный вид функции / не влияет на асимптотику Более того, асимптотика будет, в основном, определяться хвостом распределения Ре(Е), которое должно достаточно быстро спадать с ростом Е. Приведем здесь результаты для случая, когда

Ре(Е) = 2(1 -р)5(Е/Т) + арехр(-аЕ/Т), (9)

где р — вероятность встретить ловушку, а — параметр. При выполнении некоторых условий на р, выход водорода из такой системы будет состоять из двух участков, разделенных характерным временем

¿р-то/р1^1-0).

Для £ <С получается закон классической диффузии:

Зоиг ~ 1 ¡\Я.

Для 4 1р мы приходим к «аномальной» зависимости,

Эои1 ~ 1/г(1~а/2). (10)

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

Результаты этого пункта проиллюстрированы на рис. 4, где показаны результаты компьютерного моделирования блуждания частицы методом Монте-Карло для случая р = Ю-3 и а = 1/3. На рисунке показана зависимость полного количества частиц С}(Ь), вышедших к

Рис. 4. Зависимость полного количества вышедшего водорода Q(t) от времени.

моменту времени t, от времени. Так как jout = dQ/dt, в данном случае получается jout

Заметим, что такая зависимость потока газовыделения действительно иногда наблюдается в экспериментах. Так, измерения зависимости потока газовыделения из стенок от времени в токамаках Тоге Supra и JET показывают, что выход водорода происходит по асимптотическому закону 1/i0'7 и 1 /£2/3 вместо классической диффузионной зависимости 1/i0-5. Это может быть связано с тем, что значительная доля водорода накапливается в перенаправленных на стенках пленках, имеющих сложную внутреннюю структуру, которая предположительно может быть описана именно наличием дефектов с непрерывным распределением по энергиям выхода водорода.

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

Поэтому мы модифицируем уравнения диффузии в поле ловушек, чтобы учесть непрерывное распределение ловушек по энергиям. Изменения сводятся к тому, что теперь вместо различных типов ловушек необходимо ввести функцию распределения ловушек по энергиям выхода пе{Е, ж), так что пе(Е, х) <1Е есть плотность числа частиц, находящихся в ловушках с энергией выхода от Е до Е + йЕ (считаем, что энергия захвата Е^-ар = Ед одинакова для всех ловушек). Тогда уравнения, определяющие транспорт частиц в такой системе, запишутся в виде:

Здесь Ne — плотность распределения всех ловушек (с водородом и без), введенная аналогично пе, va 11 vt —характерные частоты выхода и захвата в ловушки, соответственно, D — коэффициент диффузии, ns — плотность растворенного водорода, N — плотность материала.

В случае малой концентрации водорода большинство ловушек остаются свободными. В этом приближении, считая температуру постоянной, а образец — полубесконечным, можно решить уравнения (11) методом Лапласа. Тогда, подставляя конкретные виды распределения ловушек по энергиям, мы приходим к тем же результатам, что были получены из модели случайных блужданий.

Основные результаты работы

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

оо

(И)

дпЕ

~дГ

-vde Е/Тпв + щп.

Ne - пЕ N

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

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

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

Основные результаты диссертации опубликованы в 19 печатных работах. Из них 5 в ведущих рецензируемых научных журналах, определенных ВАК.

Публикации в журналах из перечня ВАК:

1. Е. Marenkov, S. Krasheninnikov, A. Pisarev, I. Tsvetkov. On the tokamak first wall response to the ELM bursts. Pliys. Plasmas 19,

092501 (2012); doi: 10.1063/1.4742982

2. E. D. Marenkov, S. I. Krasheninnikov, A. Yu. Pigarov, A. A. Pisarev, and I. V. Tsvetkov. On thermal instability caused by plasma-wall coupling. Phys. Plasmas 18, 092502 (2011); doi: 10.1063/1.3626832 (6 pages)

3. Е.Д. Маренков, С.И. Крашенинников, A.A. Писарев, И.В. Цветков, О термической неустойчивости, вызванной взаимодействием плазмы со стенкой. Физика плазмы, 2012, том. 38, №3, с. 1-7

4. Е.Д. Маренков, И. В. Цветков, А. А. Писарев, Проницаемость изотопов водорода через многослойные мембраны. Известия РАН. Серия физическая, 2010, том 74, 2, с. 316-322

5. Е.Д. Маренков, И. В. Цветков, А. А. Писарев, Проницаемость изотопов водорода через мембраны с неоднородными по толщине свойствами. Поверхность. Рентгеновские, синхротронные и нейтронные исследования. 2008, №5, с. 1-5

Прочие публикации:

1. Маренков Е.Д., Крашенинников С.П., Гаспарян Ю.М., Писарев A.A. Транспорт водорода в материалах с широким распределением ловушек по энергиям захвата. Материалы 16 конференции «Взаимодействие плазмы с поверхностью». М.: НИЯУ МИФИ, 2013. Сборник тезисов

2. Маренков Е.Д., Смирнов Р.Д., Крашенинников С.П., Писарев A.A. О реакции первой стенки токамаков на ЭЛМ. Материалы 16 конференции «Взаимодействие плазмы с поверхностью». М.: НИЯУ МИФИ, 2013. Сборник тезисов

3. E. Marenkov, S. Krasheninnikov, A. Pisarev, I. Tsvetkov. On the tokamak first wall response to the ELM bursts. Bokk of reports of 20th international conference PSI-2012. Aachen, Germany, 2012

4. E. D. Marenkov, S. I. Krasheninnikov, A. A. Pisarev, and I. V. Tsvetkov. On thermal instability caused by plasma-wall coupling. Book of reports of 20 international conference Ion-Surface Interactions (ISI)-2011, Zvenigorod, Russia, 2011, Vol. 2, p. 275

5. E.D. Marenkov, A.A. Pisarev, I.V. Tsvetkov, Hydrogen permeation through sandwich membranes. Books of reports of 19 international conference IonSurface Interactions ISI-2009, Zvenigorod, Russia. Vol.2, p. 369.

6. E.D. Marenkov, A.A. Pisarev, I.V. Tsvetkov, Hydrogen permeation through the membrane with non uniform properties. Books of reports

of 18 international conference IonSurface Interactions ISI-2007, Zvenigorod, Russia. Vol.1, p. 303.

7. E. D. Marenkov, S. I. Krasheninnikov, A. Yu. Pigarov, A. A. Pisarev, and I. V. Tsvetkov. On the thermal instability caused by plasma-wall coupling, тезисы Book of Abstracts of 13th International Workshop on Plasma-Facing Materials and Components for Fusion Applications and 1st International Conference on Fusion Energy Materials Science, 2011, p. 119 1

8. Е.Д. Маренков, И. В. Цветков, А. А. Писарев. Транспорт водорода через мембраны с неоднородными по толщине свойствами, Труды восемнадцатой международной конференции «Взаимодействие ионов с поверхностью — 2007»

9. Е.Д. Маренков, И. В. Цветков, А. А. Писарев. Проницаемость водорода через трехслойные структуры, тезисы. Тезисы докладов

Третьей международной конференции «Взаимодействие изотопов водорода с конструкционными материалами. IHISM-2007»

10. Е.Д. Маренков, И. В. Цветков, А. А. Писарев. Проницаемость водорода через трехслойные структуры, тезисы Труды научной сессии МИФИ-2007 Т.4 Лазерная физика. Физика плазмы. Сверхпроводимость и физика наноструктур №48 13986

11. Е.Д. Маренков, И. В. Цветков, А. А. Писарев. Проницаемость изотопов водорода через мембраны с неоднородными по толщине свойствами, тезисы учебно-научной конференции-конкурса по физике УНК-2007.

12. Е.Д. Маренков, И. В. Цветков, А. А. Писарев. Проницаемость водорода через многослойные мембраны, тезисы докладов Четвертой международной конференции «Взаимодействие изотопов водорода с конструкционными материалами. IHISM-2008

13. Е.Д. Маренков, И. В. Цветков, А. А. Писарев. Проницаемость водорода через трехслойные структуры, тезисы Труды Третьей международной конференции «Взаимодействие изотопов водорода с конструкционными материалами. IHISM-2007

14. A.A. Писарев, И.В.Цветков, Е.Д. Маренков, С.С. Ярко. Проницаемость водорода через металлы. Учебное пособие. М.: МИФИ, 2008. - 144 с.

Подписано в печать 14.03.2013. Заказ № 140. Тираж 120 экз.

Типография НИЯУ МИФИ. Каширское шоссе, 31.

Текст научной работы диссертации и автореферата по физике, кандидата физико-математических наук, Маренков, Евгений Дмитриевич, Москва

НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ЯДЕРНЫЙ УНИВЕРСИТЕТ

МИФИ

На правах рукописи

04201355768

МАРЕНКОВ Евгений Дмитриевич

«У*

Влияние транспорта водорода в первой стенке термоядерных

реакторов на рециклинг

Специальность 01.04.08 — «Физика плазмы»

Диссертация па соискание учёной степени кандидата физико-математических наук

Научный руководитель: канидидат физико-математических наук

Цветков И.В.

Москва - 2013

Оглавление

Введение 4

1 Обзор литературы 10

1.1 Накопление и рециклинг водорода в материалах, контактирующих с плазмой..........................................................................10

1.1.1 Внедрение................................................................10

1.1.2 Диффузия водорода в металлах......................................11

1.1.3 Проницаемость через простую мембрану............................19

1.1.4 Проницаемость через многослойную мембрану ....................22

1.1.5 Захват водорода в металлах..........................................24

1.2 Н-режим и ЭЛМы..............................................................27

1.2.1 Основные понятия......................................................27

1.2.2 МГД теория ............................................................30

1.3 MARFE, срывы ................................................................35

1.4 Длительный разряд............................................................40

1.5 Выводы..........................................................................44

2 Термическая неустойчивость 47

2.1 Введение ........................................................................47

2.2 Десорбция первого порядка....................................................47

2.3 Случай диффузии..............................................................57

2.4 Выводы..........................................................................63

3 Реакция на ЭЛМ 64

3.1 Введение ........................................................................64

3.2 Нульмерная модель............................................................65

3.2.1 Код FACE................................................................72

3.3 Результаты моделирования....................................................75

3.4 Выводы..........................................................................80

4 Проницаемость водорода через многослойные мембраны 82

4.1 Введение........................................................................82

4.2 Основные уравнения............................................................82

4.3 Проникающий поток............................................................85

4.4 Накопление......................................................................91

4.5 Выводы..........................................................................94

5 Транспорт водорода в материале с дефектами, имеющими непрерывное распределение по энергиям выхода 96

5.1 Введение .................................... 96

5.2 Случайные блуждания в решетке с ловушками.............. 96

5.3 Уравнение с распределенными ловушками.................102

5.4 Выводы.....................................103

Заключение 105

Список литературы 114

Введение

Взаимодействие плазмы со стенкой представляет собой сложную проблему физики термоядерных реакторов. Процессы взаимодействия важны для выбора материалов стенки термоядерного реактора (хотя для проекта ITER такой выбор уже сделан [1]), а так же для его правильной эксплуатации, учитывающей ограничения, накладываемые взаимодействием со стенкой. Активные исследования в этой области, хотя и прояснили большое число вопросов, связанных со взаимодействием плазмы с твердыми телами, все же оставляют и ие меньше еще не понятых проблем. Состояние исследований в этой области отражены в большом количестве книг и обзоров (например, [2,3]). Особенно трудно дается понимание явлений, связанных со взаимным влиянием процессов, происходящих в простеночной плазме, и процессов рециклинга изотопов водорода в первой степке.

В обзоре [4], посвященном основным проблемам взаимодействия со стенкой для проекта ITER, выделены следующие три критических проблемы

1. время жизни обращенных к плазме материалов;

2. образование пыли из-за эрозии обращенных к плазме материалов;

3. накопление трития (Т) в материалах вакуумной камеры.

Последние два пункта определяют пределы безопасной эксплуатации установки по следующим причинам:

• опасность случайного выброса в окружающую среду токсичных (Be) и радиоактивных (из-за накопления в них трития) материалов;

• возможное активное выделение водорода из-за реакции с паром при возникновении случайной течи воды;

• возможный взрыв пыли, следующий за выделением водорода и случайным прорывом воздуха.

1Любо11ытно, что жесткие условия, в которых находятся материалы первой стенки токамака, приводят также к появлению технических проблем, напрямую не связанных с взаимодействием плазмы с веществом; например, в работе [5] обсуждается механизм самопроизвольного раскручивания болтов, крепящих тайлы в токамаке EAST из-за тепловых нагрузок, которым подвергается первая стенка.

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

Приняты следующие максимальные уровни по содержанию пыли и накопления трития, полученные из анализа безопасности эксплуатации реактора ITER (см. также [6]):

• один килограмм подвижного (mobilisable) трития в камере (предел «эвакуации»);

• одна тона подвижной пыли (mobilisable dust) в камере во время D-D и D-T фазы;

• 6 кг углеродной, 6 кг вольфрамовой и 6 кг бериллиевой пыли на горячих поверхностях. Этот предел обусловлен риском производства водорода. Значения соответствуют максимально допустимой концентрации водорода (2,5 кг) в камере, которое обусловлено требованием сохранности камеры в случае взрыва выделившегося водорода.

Вместе с тем, точное измереиие этих величин наталкивается на определенные трудности, в связи с чем для ITER установлены также следующие «административные» пределы, которые гарантируют, что предельные с точки зрения безопасности значения не будут достигнуты даже с учетом погрешности измерений:

• 670 кг пыли в камере;

• 700 г накопленного в камере трития.

Известны также оценки для основных параметров воздействия плазмы на стенку [7,8]. В диверторе в области ударных точек (strike-point) потоки ионов и нейтралов достигают значений больших, чем 1024м~2с-1 (что приводит к дозе > 1026м-2 на каждый импульс ITER). Плотность диверторной плазмы составляет при этом величину порядка 1021м-3 при температуре около 3 эВ. что соответствует для ионов D+ энергии падения на поверхность около 15 эВ из-за ускорения в пристеночном потенциале плазмы. Температура поверхности элементов дивертора в районе ударной точки ожидается порядка 1300 К для внешнего дивертора и около 800 К на внутреннем. Моделирование дает также ожидаемую мощность излучения в диверторе около 60 - 70 МВт. что значительно превосходит параметры современных токамаков [9]

Оценки для потоков на первую стенку менее точны. Результаты моделирования показывают, что поток нейтрального D на первую стенку будет достигать значений 1019 - 1021м-2с-1 при энергиях 8 - 300 эВ

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

работ в этой области. Разработаны компьютерные коды, позволяющие моделировать процессы переноса в пристеночной плазме (UEDGE, ERO, SOLPS и т.п.); достигнут существенный прогресс в понимании нестационарных (off-normal) процессов в тока-маке, таких как срывы, MARFE, ЭЛМы. возникновение «горячих точек» (hot spots) и т.д. Накопление и транспорт изотопов водорода в материалах, которые планируется использовать в первой стенке (бериллий, вольфрам и CFC), также являются объектом интенсивного исследования. Анализ экспериментальных данных показывает, что транспорт водорода в чистых материалах определяется процессами диффузии, захвата/выхода из дефектов различных типов, накоплением в порах и поверхностными процессами. В условиях первой стенки, однако, большую роль играет также эрозия и перенапыление эродированных частиц, что может существенно изменить свойства материала.

Вместе с тем, эксперименты на существующих в настоящее время токамаках показывают, что процессы, обычно ассоциируемые с чисто «плазменными» эффектами, оказывают влияиие на транспорт изотопов водорода в стенке, и наоборот. Примерами такого взаимного влияния могут служить резкая десорбция большого количества водорода во время MARFE и срывов [10]. Основное же количество экспериментальных работ по изучению взаимодействия плазмы со стенкой делается в лабораторных условиях ие на токамаках, а на модельных установках, в которых материал подвергается воздействию плазмы или облучению ионным пучком с параметрами, близкими к ожидаемым в термоядерных реакторах. Результаты, полученные в таких условиях, не могут учесть взаимного влияния процессов, происходящих в пристеночной плазме токамаков и процессов, отвечающих за транспорт водорода на степке. В большинстве теоретических работ и в существующих компьютерных кодах, посвященных моделированию пристеночной плазмы (UEDGE [11]. ERO [12] и другие) такое взаимное влияние также не учитывается, т.к., из-за своей сложности, даже расчет этих процессов независимо друг от друга наталкивается па существенные трудности.

Большинство конструкционных элементов токамаков, контактирующих с плазмой, представляет собой многослойные структуры, и расчет накопления и транспорта в них водорода должен учитывать это. Кроме того, использование защитных покрытий является одним из возможных способов снижения накопления водорода в контактирующих с плазмой материалах. Так. в работе [13] рассматривается применение вольфрамовых покрытий, в работе [14] — напыление а-алюмииия на сталь Eurofer для снижения проницаемости через нее водорода, в работе [15] — покрытие стали оксидом эрбия и т.д. Поэтому задача проницаемости водорода через многослойные материалы привлекает к себе особое внимание [16-18].

Задача проницаемости из газа через перегородки («мембраны») имеет также другое применение Любое моделирование поведения водорода в материале требует знания большого числа параметров (коэффициенты диффузии, растворимости, реком-

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

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

Несмотря на большое количество различных процессов, участвующих в транспорте водорода в твердом теле, па практике транспорт чаще всего описывается с помощью решения уравнений переноса (уравнение диффузии для водорода в поле ловушек). Существует несколько кодов для решения этих уравнений, из которых код ТМАР является одним из наиболее часто используемых. Вместе с тем, в коде ТМАР7 может рассматриваться не более трех различных типов ловушек, отличающихся энергиями связи с ними водорода. Однако, некоторые эксперименты могут быть удовлетворительно объяснены только при учете значительно большего количества типов ловушек (например, термодесорбционные спектры водорода из вольфрама с примесыо карбида титана; вольфрама, подвергшегося нейтронному облучению и т.п.). В связи с этим разумно предположить, что в материалах со сложной структурой (например, перепапыленпые слои па первой стенке) более адекватной моделью транспорта является введение непрерывного распределения ловушек по энергиям связи с ними водорода.

Исходя из вышеперечисленного, была сформулирована цель данной работы.

Цель

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

2. Исследовать роль первой стенки в восстановлении пьедестала после вспышки ЭЛМа.

На защиту выносятся следующие результаты, содержащие научную новизну:

3. Компьютерный код, основанный на разработанной модели транспорта изотопов водорода в материалах с учетом теплопроводности, транспорта и возникновения дефектов, взаимодействия водорода с дефектами и междоузельпыми атомами, диффузии растворенного водорода. Результаты компьютерного моделирования реакции стенки на вспышку ЭЛМа, оценка влияния бериллиевой стенки на восстановление пьедестала Н-моды.

Научная и практическая значимость работы

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

Апробация работы

Основные результаты работы докладывались па следующих конференциях:

1. 16 конференция «ВЗАИМОДЕЙСТВИЕ ПЛАЗМЫ С ПОВЕРХНОСТЬЮ», Москва, НИЯУ МИФИ, 2013

2. 13th International Workshop on Plasma-Facing Materials and Components for Fusion Applications and 1st International Conference on Fusion Energy Materials Science (PFMC), Rosenheim, Germany, 2011

3. 20th international conference PSI-2012. Aachen, Germany, 2012

4. 15 конференция «Взаимодействие плазмы с поверхностью». Москва. НИЯУ МИФИ, 2012

5. 20 international conference Ion-Surface Interactions (ISI)-2011, Zvenigorod, Russia, 2011

6. 19 international conference IonSurface Interactions ISI-2009, Zvenigorod, Russia

7. 18 international conference IonSurface Interactions ISI-2007, Zvenigorod, Russia

8. International conference of atomic collisions on surfaces ICACS-24, Krakow, Poland, 2012

11. Научная сессия МИФИ-2007, Москва, НИЯУ МИФИ

12. Учебно-научная конференция-конкурс по физике УНК-2007. Москва, ФИАН

Объем и структура диссертации

Глава 1

Обзор литературы

1.1 Накопление и рециклинг водорода в материалах, контактирующих с плазмой

1.1.1 Внедрение

В токамаках материалы, обращенные к плазме, бомбардируются ионами и нейтральными частицами, обладающими энергией до нескольких килоэлектронвольт (кэВ). Некоторые из падающих частиц отражаются от материала с сохранением существенной доли их изначальной энергии. Большая часть частиц при отражении нейтрализуется. Доля частиц, которые отражаются от материала, выше для падающих частиц с более низкой энергией и для материалов-мишеней с более высоким атомным числом Z и может превышать 50% [19]. Частицы, которые не отражаются, проникают в глубину материала, теряя свою кинетическую энергию в результате упругих и неупругих соударений вплоть до их полной остановки. Распределение этих внедренных частиц по глубине зависит от энергии и сорта падающих частиц, а также от материала мишени.

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