Развитие моделей поведения ядерного топлива в условиях повышенного выгорания, переходных режимов и при пенальном контроле герметичности оболочек твэлов тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат физико-математических наук Зборовский, Вадим Гарольдович

• Зборовский, Вадим Гарольдович
• кандидат физико-математических наук
• 2012

Оглавление диссертации кандидат физико-математических наук Зборовский, Вадим Гарольдович

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

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

1.1. Формирование рим-структуры в облучаемом и02-топливе.

1.2. Модели поведения точечных дефектов в облучаемом 1Юг.

1.3. Модели поведения газовых продуктов деления в и02-топливе.

1.4. Модели роста зерна при облучении иОг-топлива в условиях переходных режимов.

1.5. Модели дрейфа газовых пузырьков в градиенте температуры.

1.6. Основные физические модели процессов в твэле при пенальном КГО на остановленном реакторе.

Глава 2. Моделирование пороговых условий для формирования рим-структуры

2.1. Аналитическое исследование неустойчивости Киношита.

2.2. Моделирование поведения точечных дефектов вблизи дислокации.

2.3. Самосогласованное моделирование точечных дефектов и атомов ГПД вблизи дислокации.

2.4. Численные методы для моделирования кинетики точечных дефектов и атомов ГПД вблизи дислокации.

2.5. Численная реализация модели поведения точечных дефектов и ГПД вблизи дислокации.

2.6. Выводы к главе 2.

Глава 3. Моделирование выхода газовых продуктов деления и микроструктурных изменений при облучении топлива в условиях переходных режимов.

3.1. Модель роста зерна во внутриреакторных условиях для топлива с крупным начальным зерном и легирующими добавками.

3.2. Верификация моделей роста зерна на экспериментах ШбиЗ и НВЕР.

3.3. Модели ускоренного выхода ГПД в условиях переходных режимов.

3.4. Верификация моделей выхода ГПД в переходных режимах.

3.5. Выводы к главе 3.

Глава 4. Моделирование выноса активности из дефектного твэла при контроле герметичности оболочек на остановленном реакторе кодом РТОП-КГО.

4.1. Физические модели кода РТОП-КГО.

4.2. Численные схемы, реализованные в коде РТОП-КГО.

4.3. Валидация и верификация кода РТОП-КГО.

4.4. Выводы к главе 4.

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

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

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

Для увеличения эффективности использования ядерного топлива в настоящее время повышается проектное выгорание топлива, и внедряются эксплуатационные режимы с маневрированием мощностью [2]. В рамках проводимых работ является важным развитие ряда моделей. К таковым относится моделирование 1Ю2 топлива высокого выгорания. На периферии топливных таблеток в твэлах легководных реакторов при средних по таблетке выгораниях более 40 МВт-сут/кгИ формируется так называемая рим-структура [3,4], характеризуемая дроблением исходных зерен и значительным ростом пористости. Изменения микроструктуры влияют на теплофизические и механические свойства топлива в области рим-слоя, а также на выход активности из твэла в случае разгерметизации. Самосогласованная физическая модель процессов, определяющих формирование рим-слоя в 1Ю2-топливе, до сих пор отсутствует, несмотря на развитую феноменологию данного явления [4]. Поэтому задача о получении критерия для порога формирования рим-слоя в зависимости от различных параметров топлива и условий облучения, является актуальной.

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

Для повышения радиационной безопасности на АЭС необходимым является развитие методик контроля герметичности оболочек (КГО) твэлов на остановленном реакторе. В настоящее время для российских атомных станций разработана и внедрена усовершенствованная методика пенального КГО с цитированием давления [5,6]. Данная методика обладает повышенной чувствительностью по сравнению со штатной методикой КГО, а измерение кинетики выноса активности дает возможность оценивать величину эффективного гидравлического диаметра дефекта в оболочке негерметичного твэла. Знание размера дефекта позволяет прогнозировать активность в первом контуре в случае загрузки негерметичной TBC. Для обоснования новой методики было необходимым разработать расчетный код, описывающий физические процессы, которые определяют кинетику выноса активности из негерметичного твэла при пеналь-ном КТО.

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

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

Научная новизна результатов, представленных в диссертации:

1. Разработана расчетно-теоретическая модель поведения точечных дефектов и атомов ГПД вблизи краевой дислокации в условиях низкотемпературного облучения иОг.

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

3. Определены условия формирования газовых пузырьков разного размера вблизи краевых дислокаций в облучаемом иС^-топливе.

4. Развиты расчетно-теоретические модели, описывающие ускоренный выход ГПД из иОг-топлива в переходных режимах: модели роста зерна во'внутриреакторных условиях и модель дрейфа газовых пузырьков в градиенте температуры.

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

6. Создан расчетный код РТОП-КГО, моделирующий вынос активности из дефектного твэла при проведении пенального КГО.

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

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

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

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

Достоверность результатов, полученных в диссертационной работе, подтверждается сравнением расчетов с экспериментальными данными. Модели микроструктурных изменений в топливе и выхода ГПД в переходных режимах верифицированы на реакторных экспериментах. Код РТОП-КГО верифицирован на экспериментах на стенде ПЕНАЛ, результаты его применения на АЭС подтверждаются послереакторными исследованиями облученных TBC. Проведено сравнение результатов расчетов с аналитическими тестами, и обоснована сходимость использованных численных схем. Программный код РТОП-КГО аттестован в Федеральной службе по экологическому, технологическому и атомному надзору РФ, получен аттестационный паспорт программного средства №221 от 21.02.2007.

На защиту выносятся следующие основные положения:

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

2. Формирование пространственно-неоднородных распределений точечных дефектов и атомов ГПД вблизи дислокации может определять пороговые условия формирования рим-структуры.

3. Развитые модели роста зерен и ускоренного выхода ГПД хорошо описывают поведение иОг-топлива с крупным зерном в условиях переходных режимов.

4. Созданный расчетный код РТОП-КГО решает задачу переноса жидкости и растворенных радионуклидов под оболочкой дефектного твэла и определения кинетики выноса активности из негерметичного твэла при пенальном КГО.

Апробация работы. Результаты работы докладывались и обсуждались на 6й, 7й, 8й, 9й международных конференциях "International Conference on WWER Fuel Performance, Modelling and Experimental Support" в Болгарии (2005, 2007, 2009, 2011 гг); на конференциях серии Top Fuel: "Water Reactor Fuel Performance Meeting" (Япония, 2005 г), "International Meeting on LWR Fuel Performance" (Испания, 2006 г, США, 2010 г); на IX Российской конференции по реакторному материаловедению (г. Димитровград, 14-18 сентября 2009 г); на международном семинаре "NXO International workshop: Radiation Effect and Self-Recovery in Materials" (Япония, 2004 г); на российской научно-технической конференции ОАО «ТВЭЛ» (ОАО ВНИИНМ, 2010 г); на российской межотраслевой научно-технической конференции «Исследовательскому комплексу ИВВ-2М - 45 лет» (Заречный, 2011 г); на российском семинаре «Теория и многоуровневое моделирование дефектов, явлений и свойств материалов ядерной техники - ТММ-2008» (ОАО ВНИИНМ, 2008 г); на российском семинаре «Физика радиационных повреждений материалов атомной техники» (Обнинск, 2004, 2005, 2006, 2009, 2010 гг); семинаре «Физическое моделирование изменения свойств реакторных материалов в номинальных и аварийных условиях» (Димитровград, 2006,2009 гг; Троицк, 2005 г).

Публикации. Материалы, вошедшие в диссертацию, опубликованы в 12 печатных работах, из них 5 статей в рецензируемых журналах, 7 докладов, опубликованных в сборниках трудов конференций.

Личный вклад автора состоит в разработке, реализации и сравнении с экспериментами расчетно-теоретических моделей поведения точечных дефектов и газовых продуктов деления в облучаемом и02-топливе (в условиях формирования рим-структуры и в условиях переходных режимов). Автором разработан модуль конвективно-диффузионного переноса и растворения радионуклидов, а также гидравлический модуль (последний совместно с Афанасьевой Е.Ю.) кода РТОП-КГО. Автор проводил верификацию разработанных расчетно-теоретических моделей и кода РТОП-КГО в целом.

Структура и объем диссертации. Настоящая диссертация состоит из введения, четырех глав, заключения, списка цитируемой литературы. Работа включает в себя 144 страницы, 52 иллюстрации, 5 таблиц и 163 цитирования литературы. Во введении излагаются предпосылки, послужившие выбору темы диссертации, цель проведенного исследования, формулируются научная новизна и практическая значимость представленной работы, положения, выносимые на защиту, список публикаций и личный вклад автора. В первой главе излагаются основные представления, опубликованные в литературе, и использованные в настоящей работе. Во второй главе моделируются пороговые условия формирования рим-структуры в облучаемом иОг-топливе, для этого численно решается задача о кинетике точечных дефектов и распределении ГПД вблизи дислокации. Полученные результаты позволяют проанализировать модели неустойчивости Киношита и формирования крупных пересжатых пузырьков в топливе. В третьей главе описаны расчетные модели ускоренного выхода ГПД из иОг-топлива в переходных режимах облучения и приведены результаты их применения для условий экспериментов. В четвертой главе излагаются физические и численные модели кода РТОП-КГО, которые описывают массоперенос

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Зборовский, Вадим Гарольдович

4.4. Выводы к главе 4

В данной главе рассмотрен код РТОП-КГО, моделирующий вынос активности из негерметичного твэла при КГО на остановленном реакторе. Представлены физические модели и их численные реализации.

В коде используются обоснованные численные схемы, которые являются консервативными и имеют порядок сходимости не ниже первого по пространству и по времени, а также 2й порядок сходимости по пространству во многих случаях. Это подтверждаются численной оценкой скорости сходимости и вали-дацией на аналитических тестах. Использование ТУЭ-схемы для уравнения переноса позволило получать монотонные решения при отсутствии значительной схемной диссипации. Код верифицирован на большом числе экспериментов, выполненных на стенде ПЕНАЛ, и аттестован в надзорных органах РФ [163].

Результаты, полученные в данной главе, имеют практическое значение. Код РТОП-КГО успешно применяется для решения прямой и обратной задач определения кинетики выноса активности из дефектного твэла при КГО на остановленном реакторе. Методика пенального КГО с циклированием давления апробирована на российских АЭС, в частности, результаты ее применения подтверждаются послереакторными исследованиями герметичных и негерметичных TBC. С использованием кода РТОП-КГО разработаны атласы карт кинетических зависимостей для топлива с центральным отверстием после 1-3 лет эксплуатации и топлива без центрального отверстия после 1 года эксплуатации. В настоящее время методика с циклированием давления включена в Руководящие Документы [127,128] и внедрена на все российские энергоблоки с реакторами ВВЭР.