Определение динамических нагрузок при термическом взаимодействии кориума с теплоносителем в ходе тяжёлых аварий на аэс
На правах рукописи
Ртищев Никита Александрович ОПРЕДЕЛЕНИЕ ДИНАМИЧЕСКИХ НАГРУЗОК ПРИ ТЕРМИЧЕСКОМ ВЗАИМОДЕЙСТВИИ КОРИУМА С ТЕПЛОНОСИТЕЛЕМ В ХОДЕ ТЯЖЁЛЫХ АВАРИЙ НА АЭС Специальность 05.14.03 – Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата технических наук
Москва - 2013
Работа выполнена в ФГБОУ ВПО «НИУ«МЭИ» на кафедре Атомных электрических станций
Научный консультант: доктор технических наук, начальник вычислительного отдела ОАО «ЭНИЦ», профессор кафедры АЭС НИУ «МЭИ» Мелихов Владимир Игорьевич
Официальные оппоненты: доктор технических наук, профессор кафедры Инженерной теплофизики НИУ «МЭИ» Кузма-Кичта Юрий Альфредович доктор физико-математических наук, начальник научно исследовательского сектора Института экспериментальной газодинамики и физики взрыва ФГУП «РФЯЦ-ВНИИЭФ» Дреннов Олег Борисович
Ведущая организация: Институт теплофизики Сибирского отделения РАН
Защита состоится «26» июня 2013 года в 14 часов, на заседании диссертационного совета Д 212.157.07 ФГБОУ ВПО «НИУ«МЭИ» по адресу:
111250, г. Москва, ул. Красноказарменная, д.14, Малый актовый зал МЭИ.
С диссертацией можно ознакомится в библиотеке ФГБОУ ВПО «НИУ«МЭИ»
Автореферат разослан «_» мая 2013 г.
Ученый секретарь диссертационного совета к.т.н. И.П.Ильина
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы. Требование обеспечения высокого уровня безопасности и надежности АЭС диктует необходимость выполнения углубленных оценок безопасности, включая рассмотрение запроектных аварий.
Тяжелая авария на АЭС, вызывающая разрушение и плавление активной зоны реактора, является одним из самых опасных по своим последствиям типов запроектной аварии. В ходе ее развития возможны ситуации, приводящие к прямому контакту высокотемпературного расплава материалов активной зоны (кориума) с теплоносителем и их термическому взаимодействию, которое при определенных условиях может осуществляться взрывным образом вследствие мощного когерентного вскипания теплоносителя и сопровождаться динамическим воздействиями на окружающие конструкции. Вследствие возникающих угроз целостности АЭС при таком развитии аварийных событий термическое взаимодействие кориума с теплоносителем и соответствующие динамические воздействия уже много лет является важным предметом исследований специалистов в области безопасности атомной энергетики. Ввиду исключительной сложности реализующихся при термическом взаимодействии многофазных неравновесных многоскоростных процессов, сопровождающихся неустойчивостями различной природы, исчерпывающая теория этого явления еще далека до своего завершения. Необходимы дальнейшие экспериментальные и теоретические исследования, разработка, усовершенствование и валидация расчетных кодов для анализов термического взаимодействия, квалифицированное проведение оценок динамических воздействий для АЭС различного типа с помощью этих кодов.
Диссертационная работа посвящена исследованиям по этим направлениям, что и определяет ее актуальность.
Цель работы состоит в определении динамических нагрузок при термическом взаимодействии кориума с теплоносителем в ходе тяжелой аварии на АЭС.
Основными задачами
исследования являются:
- аналитический обзор исследований термического взаимодействия кориума с теплоносителем при тяжёлых авариях на АЭС;
- разработка и усовершенствование математических моделей термического взаимодействия кориума с теплоносителем, доработка/создание соответствующих расчетных кодов;
- валидация этих кодов на экспериментальных данных по термическому взаимодействию кориума с теплоносителем;
- определение динамических нагрузок при термическом взаимодействии кориума с водой в корпусе реактора ВВЭР-1000 (внутрикорпусной паровой взрыв);
- определение динамических нагрузок при термическом взаимодействии кориума с водой в шахте реактора ВВЭР-1000 (внекорпусной паровой взрыв);
- определение динамических нагрузок при термическом взаимодействии кориума с натрием в корпусе реактора БН-600.
Научная новизна. В диссертационной работе:
1. Усовершенствована карта режимов течения, входящая в состав математической модели кода VAPEX-P, предназначенного для моделирования стадии перемешивания кориума с водой.
2. Выполнена валидация кода VAPEX-P на опытных данных экспериментов FARO L-14, FARO L-24 и FARO L-28.
3. Разработана интегральная модель термического взаимодействия кориума с натрием и на ее основе создан расчетный код CSI.
4. Выполнена валидация кода CSI на опытных данных экспериментов THINA, ПЛУТОН.
Достоверность результатов выполненных исследований подтверждается:
- применением общепризнанных методов и подходов механики многофазных сред при создании/доработки математических моделей;
- использование апробированных корреляций и соотношений для описания трения и тепло и массообмена между фазами;
- результатами валидации кодов VAPEX-P, VAPEX-D, CSI.
Практическая ценность проведенного исследования состоит в применении доработанных и валидированных расчетных кодов VAPEX-P и VAPEX-D для оценок динамических воздействий при внутрикорпусном и внекорпусном паровых взрывах на АЭС с ВВЭР и в применении созданного и валидированного расчетного кода CSI для оценок динамических воздействий при термическом взаимодействии кориума с натрием в ходе тяжелой аварии на АЭС с БН.
Все основные этапы исследования выполнялись по договору между Проектно-конструкторского филиалом ОАО «Концерн Росэнергоатом» и НИУ «МЭИ» № 2077110 от 04.04.2011 г. и Государственному контракту с Министерством образования и науки РФ ГК № 14.740.11.0093 от 08.09 2010 г.
Положения, выносимые на защиту:
1. Результаты валидации кода VAPEX-P на опытных данных экспериментов FARO L-14, FARO L-24 и FARO L-28.
2. Математическая модель и расчетный код CSI для анализа термического взаимодействия кориума с натрием.
3. Результаты валидации созданного расчетного кода CSI на опытных данных экспериментов THINA, ПЛУТОН.
4. Значения динамических нагрузок при внутрикорпусном паровом взрыве в корпусе реактора ВВЭР-1000.
5. Значения динамических нагрузок при внекорпусном паровом взрыве в шахте реактора ВВЭР-1000.
6. Значения динамических нагрузок при термическом взаимодействии кориума с натрием в корпусе реактора БН-600.
Личный вклад автора. Все этапы работы по доработке математической модели кода VAPEX-P, валидации кода VAPEX-P на опытных данных экспериментов FARO L-14, FARO L-24 и FARO L-28, созданию интегральной математической модели термического взаимодействия кориума с натрием и на ее основе кода CSI, валидации кода CSI на опытных данных экспериментов THINA, ПЛУТОН, определению динамических нагрузок при внутрикорпусном и внекорпусном паровых взрывах на АЭС с ВВЭР, а также при термическом взаимодействии кориума с натрием в ходе тяжелой аварии на АЭС с БН были выполнены непосредственно автором, либо проходили при его непосредственном участии.
Апробация работы. Результаты диссертационной работы докладывались и обсуждались на конференции–конкурсе молодых физиков (Москва, 13 февраля 2012 г.);
XIII школе молодых учёных ИБРАЭ РАН (26 – 27 апреля 2012 г.);
XVIII международной научно-технической конференции студентов и аспирантов "РАДИОЭЛЕКТРОНИКА, ЭЛЕКТРОТЕХНИКА И ЭНЕРГЕТИКА" (Москва, 1-2 марта 2012 г.);
Научно-технической конференции «Теплофизика реакторов на быстрых нейтронах (Теплофизика 2012) (Обнинск, 24-26 октября 2012 г.);
международном совещании участников международного проекта SERENA, посвященного исследованиям термического взаимодействия кориума с теплоносителем (Франция, Кадараш, 5 - 7 ноября 2012 г.);
XIХ международной научно-технической конференции студентов и аспирантов "РАДИОЭЛЕКТРОНИКА, ЭЛЕКТРОТЕХНИКА И ЭНЕРГЕТИКА" (Москва, февраля-1 марта 2013 г.).
Публикации. Основные результаты работы были изложены в трех статьях, опубликованных в журналах из списка ВАК "Фундаментальные исследования" "Вестник МЭИ", в одной статье в рецензируемом международном журнале "Science and Technology of Nuclear Installations" (импакт-фактор 0,562), включенном в реферативную базу Scopus, и в пяти докладах, опубликованных в трудах отечественных и международных конференций.
Структура и объем диссертации. Диссертация состоит из введения, четырех глав, заключения, списка литературы из 158 наименований.
Диссертация содержит 228 страниц текста, в том числе 13 таблиц и рисунков.
КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении показана актуальность проведенных в диссертации исследований, определены цели и задачи работы и дана общая постановка решаемых в диссертации задач.
В первой главе приводится аналитический обзор исследований термического взаимодействия кориума с теплоносителем при тяжёлых авариях на АЭС. В обзоре отмечается, что термическое взаимодействие расплавленных материалов активной зоны (кориум) с теплоносителем, возникающее в ходе тяжёлой аварии на АЭС с ВЭЭР и АЭС с БН, может вызывать существенные динамические воздействия на окружающие конструкции и создавать угрозу их целостности. Различие между типичными сценариями протеканий тяжёлых аварий на АЭС с ВВЭР и АЭС с БН приводит к разным начальным конфигурациям системы «кориум-теплоноситель», определяющим последующий характер протекания процесса термического взаимодействия кориума с водяным и с натриевым теплоносителями.
Наличие плёночного кипения в системе «кориум-теплоноситель» оказывает существенное влияние на процесс их термического взаимодействия.
В системе «кориум-вода» на поверхности крупных капель устанавливается устойчивый режим плёночного кипения, что способствует формированию значительной области перемешивания, содержащей большое количество тепловой энергии, высвобождение которой может привести к мощному паровому взрыву.
При первоначальном контакте кориума с натрием, для типичных условий тяжёлой аварии в реакторах с натриевым теплоносителем, устойчивое плёночное кипение на каплях кориума не реализуется. Поэтому при термическом взаимодействии кориума с натрием не происходит формирование протяжённой области перемешивания, а зона их взаимодействия определяется масштабом первоначального контакта. Передача тепла от кориума в натрий вызывает быстрое расширение кипящего натрия в зоне взаимодействия.
Расширяющаяся зона взаимодействия приводит в движение прилегающий к ней «холодный» натрий, что приводит к динамическим нагрузкам на окружающие конструкции.
В обзоре показано, что основной неразрешенной проблемой при изучении термического взаимодействия в системе «вода-кориум» на данный момент является отсутствие паровых взрывов в экспериментах KROTOS для системы натурный кориум - вода, хотя в экспериментах, выполненных в НПО «Луч», а также на экспериментальных установках FARO и TROI, паровые взрывы были зафиксированы. При этом для других систем (Al2O3-вода, олово вода) на установке KROTOS были осуществлены паровые взрывы.
Для оценки параметров термического взаимодействия в системе «кориум вода» и возникающих при этом динамических нагрузок созданы расчётные коды, валидированные на имеющихся опытных данных. В расчётах тяжёлых аварий на АЭС с ВВЭР/PWR по этим кодам событие, инициирующее паровой взрыв консервативно постулируется.
Основная проблема при изучении термического взаимодействия в системе «натрий-кориум» состоит в отсутствии механистических моделей начальной стадии перемешивания кориума с натрием и его фрагментации при разрушении оболочек твэл. Как правило, предполагается или мгновенное перемешивание, или упрощённое параметрическое описание.
Ввиду малой изученности процесса перемешивания кориума с натрием, для оценки параметров термического взаимодействия в этой системе и возникающих при этом динамических нагрузок используются (или дорабатываются) упрощённые математические модели и расчётные коды, в которых процессы перемешивания и фрагментации моделируются упрощённо.
Во второй главе излагаются разработанные математические модели для анализа взаимодействия кориума с теплоносителем.
На основе подходов механики многофазных сред сформулированы математические модели, описывающие пространственное перемешивание кориума с водой и распространения волны термической детонации с её последующим расширением.
Модель перемешивания кориума с водой рассматривает следующие фазы: вода, пар, струя кориума, капли кориума и пористая структура, образующаяся в результате падения и затвердевания кориума. Описание движения воды, пара и струи кориума основывается на эйлеровом (континуальном) подходе, а динамика капель кориума описывается с помощью лагранжева подхода. Используется цилиндрическая система координат в осе симметричном случае. Для описания межфазного трения, тепло и массобмена используется подробный набор корреляций и соотношений, основанный на применении расширенной автором диссертации карты режимов течения двухфазного пароводяного потока (пузырьковый, подразделяющийся на режим мелких пузырьков, переходный режим и режим крупных паровых снарядов;
эмульсионный режим и капельный). Для описания фрагментации капель кориума используется модель Пилча-Эрдмана. Фрагментация струи кориума описывается корреляцией Саито:
f 2 Lbrk FrJ 2. a DJ где: Lbrk - длина полного распада струи, DJ - диаметр струи при входе в теплоноситель, f - плотность расплава, a - плотность теплоносителя, UJ - число Фруда, U J - скорость струи при входе в теплоноситель.
FrJ gDJ Данная модель была реализована в виде кода VAPEX-P, предназначенного для анализа процессов перемешивания кориума с водой при тяжёлых авариях на АЭС с ВВЭР/PWR.
При разработке модели распространения волны термической детонации была применена концепция микровзаимодействий кориума с водой. Эта концепция предполагает, что образующиеся при дроблении крупных капель кориума мелкие фрагменты мгновенно смешиваются с находящейся рядом водой, которая сразу же испаряется, и образуется гомогенная равновесная среда (фрагмент +.окружающая его паровая область). Остальная часть воды, находящаяся «далеко» от фрагментов – не участвует в таком «микровзаимодействии». Модель распространения волны термической детонации рассматривает в рамках эйлерового подхода следующие фазы:
крупные капли кориума, мелкие капли кориума (фрагменты кориума), пар вблизи фрагментов кориума, вода вдали от фрагментов кориума. Скорость фрагментации крупных капель кориума в волне термической детонации описывается корреляцией Патела и Теофануса, основанной на концепции неустойчивости Релея-Тейлора. Особенностью модели микровзаимодействий является унос (захват) части «удаленной» воды в «близлежащую», находящуюся в фазе микровзаимодействий, в результате фрагментации материнской капли кориума, при этом скорость «захвата» определяется пропорционально скорости фрагментации крупных капель кориума.
Данная модель была реализована в виде кода VAPEX-D, предназначенного для анализа процессов распространения волны термической детонации в системе «кориум-вода» и последующего расширения.
Излагается разработанная автором интегральная модель термического взаимодействия кориума с натрием. Предполагается, что происходит мгновенное перемешивание фрагментов кориума с натрием. При этом процесс термического взаимодействия происходит через две последовательные фазы.
На первой фазе происходит прямой контакт натрия с частицами кориума и нагрев натрия. Как только натрий достигает температуры насыщения наступает вторая фаза процесса, характеризующаяся возникновением плёночного кипения натрия на частицах кориума и резким расширением зоны взаимодействия, приводящим к возникновению динамических нагрузок на окружающие конструкции.
Данная модель была реализована в виде кода CSI, предназначенного для анализа процессов термического взаимодействия кориума с натрием и расчёта возникающих динамических нагрузок.
В третьей главе описывается результаты валидации кодов VAPEX-P и CSI, описывающих термическое взаимодействие кориума с теплоносителем.
Валидация кодов VAPEX-P и VAPEX-D, описывающих термическое взаимодействие расплава с теплоносителем на стадии предварительного перемешивания и на последующих стадиях термической детонации и расширения, была выполнена в 1998 – 2008 без участия автора на большом объёме экспериментальных данных, а с участием автора были выполнены валидационные расчёты кодом VAPEX-P экспериментов FARO-L14, FARO L24, FARO-L28, в ходе разработки модуля VAPEX-M для интегрального системного одномерного кода СОКРАТ.
Крупномасштабная установка FARO предназначена для исследования термического взаимодействия расплавленного кориума с водой. Она состоит из пяти основных компонент: плавильной камеры, системы отсечки плавильной камеры от рабочего участка, резервуара для накопления и последующей подачи расплава в рабочий участок, собственно рабочего участка, где происходит взаимодействие расплава с водой, а также системы отвода пара и неконденсируемых газов, позволяющей сбрасывать давление в рабочем участке. В экспериментах использовалось до 180 кг расплава UO2-ZrO2(-Zr) при температуре около 3100 K.
Основным параметром, характеризующим интенсивность протекания процесса термического взаимодействия, является рост давления в системе. На рисунке 1 показаны результаты расчёта давления кодом VAPEX-P эксперимента FARO-L28. Видно, что код VAPEX-P хорошо качественно и количественно описывает временную зависимость давления в системе.
Рисунок 1 – Временная зависимость давления для эксперимента FARO-L Были выполнены валидационные расчёты кодом CSI экспериментов по термическому взаимодействия кориума с натрием, проведённых на установках THINA и ПЛУТОН.
На установке THINA исследовалось термическое взаимодействие расплава оксида алюминия и железа с натрием. Расплав, массой до 5 кг, подготавливался в тигле и подавался в нижнюю часть рабочего участка, заполненного натрием. С помощью кода CSI были проведены валидационные расчёты экспериментов TH561, TH562, TH564, TH567. В расчётах варьировался параметр отношения массы расплава к массе натрия в зоне взаимодействия, а также использовались различные модели плёночного кипения. Анализ результатов проведённых расчётов показал, что наилучшее совпадение с экспериментальными данными достигается при использовании модели плёночного кипения, предложенной L. Caldarola, и отношении массы расплава = 0,83. Расчётная временная зависимость давления натрия в к массе натрия зоне взаимодействия на начальном этапе процесса достаточно хорошо совпадает с экспериментальной, рисунок 2.
Рисунок 2 –Расчётные и временные зависимости давления натрия для эксперимента TH Стенд ПЛУТОН предназначен для экспериментального определения коэффициентов конверсии термического взаимодействия и перемещения материалов при термическом взаимодействии имитаторов кориума с натрием. В экспериментах в качестве имитатора кориума использовалась термитная смесь массой 70 г, которая выливалась в сосуд, заполненный натрием, недогретым до температуры насыщения при атмосферном давлении. В таблице 1 приведено сопоставление экспериментальных значений коэффициентов конверсии с расчётными значениями, полученными кодом CSI с использованием модели плёночного кипения L. Caldarola и отношении массы расплава к массе натрия.
Таблица 1 – Сравнение коэффициентов конверсии Экспериментальное значение Расчётное значение Эксперимент Н2 0,125 0, Н4 0,41 0, Проведённые валидационные расчёты показали, что, несмотря на интегральный характер математической модели кода CSI, этот код позволяет адекватно воспроизводить основные параметры термического взаимодействия кориума с натрием.
Четвертая глава посвящена определению динамических нагрузок на конструкции при термическом взаимодействии кориума с теплоносителем в ходе тяжёлой аварии на АЭС.
В первом разделе этой главы выполнена оценка последствий внутрикорпусного парового взрыва для ВВЭР-1000. Рассматривается следующая ситуация: в центральной части активной зоны формируется область, заполненная расплавленными материалами активной зоны. Кориум вытекает через отверстие в виде струи в нижнюю часть корпуса реактора, заполненную водой. Геометрические и теплофизические параметры приведены в таблице 2. В расчётах варьировались диаметр струи кориума и уровень воды в нижней части сосуда.
Таблица 2 – Основные параметры расчётов Масса расплава, т Температура кориума, К Состав кориума 80%UO2 + 20%ZrO Высота падения струи кориума, м Диаметр струи кориума, см 10/20/ Давление, МПа 0, Температура воды, К 393 (температура насыщения) Уровень воды, м 1/ Радиус сосуда, м Основной целью расчетов перемешивания кориума с водой являлось определение массы жидких (еще не успевших затвердеть) капель кориума, образующихся в результате фрагментации струи, окруженных пароводяной смесью с низким паросодержанием, поскольку именно в таких областях («капли кориума - вода»), возможна очень быстрая передача тепловой энергии кориума к пароводяной смеси (паровой взрыв). Исследования взрывного взаимодействия капель кориума с водой свидетельствуют, что способность к возникновению и распространению значительных волн давления резко уменьшается, если капли расплава находятся в пароводяной смеси с паросодержанием () больше 0,5-0,6. В расчётах консервативно полагается, что паровой взрыв возможен, если капли кориума находятся в пароводяной смеси с объёмным паросодержанием меньшим 0,6. Расчёты выполнялись с помощью кода VAPEX-P.
Было проведено девять расчетов кодом VAPEX-P для различных диаметров вытекающего расплава (10 см, 20 см, 40 см) и различных начальных уровней воды в нижней части корпуса реактора (1 м, 2 м, 3 м). Для каждого варианта была рассчитана масса жидких капель расплава, окруженная достаточным количеством воды (объемное паросодержание меньше 0,6).
Предполагается, что в этом случае данная часть расплава способна участвовать во взрывном взаимодействии. Масса таких капель оказалась существенно меньше всей массы расплава и составила 6 – 140 кг.
Предполагается, что после удара капель кориума о дно корпуса реактора в области жидких капель расплава, окруженных пароводяной смесью с объемным паросодержанием меньшим 0,6, начинается взрывное взаимодействие кориума с водой по механизму тонкой фрагментации с образованием мелких частиц (0,1 – 0,2 мм).
Для оценки кинетической энергии, способной выделиться в ходе взрывного взаимодействия кориума с водой, используется коэффициент конверсии:
K E где K – кинетическая энергия движущегося потока, образующегося в результате парового взрыва, E – тепловая энергия, запасенная в массе кориума, взаимодействующего с водой.
Максимально возможное значение коэффициента конверсии, полученное с помощью термодинамических оценок, составляет ~ 30%. Проведенные многочисленные экспериментальные исследования показали, что реальное значение коэффициента конверсии минимум на порядок меньше (~3%), а наиболее вероятное значение составляет доли процента (~0,3-0,5%).
В таблице 3 приведены значения кинетической энергии, способной выделится при взрывном взаимодействии кориума с водой, при трех значениях коэффициента конверсии: 30%, 3% и 0,3%. Оценки этой величины сделаны для трех случаев: 1) во взаимодействии участвует вся масса вытекающего кориума, 2) во взаимодействии участвуют все образующиеся капли кориума, 3) во взаимодействии участвуют только жидкие капли кориума, окруженные достаточным количеством воды (объемное паросодержание меньше 0,6).
На основе анализа литературы, в качестве минимального значения кинетической энергии потока, способной привести к деформации корпуса реактора, выбрано значение 1 ГДж.
Из таблицы 3 следует, что соизмеримое с опасной величиной кинетической энергии значение ~ 0,8 ГДж получается при взаимодействии всей массы вылившегося расплава (1639 кг) с водой при максимальном коэффициенте конверсии 30% для варианта «уровень воды 3м, диаметр струи 10 см». Однако, детальный расчет процесса перемешивания расплава с водой показывает, что в паровом взрыве может участвовать только примерно в десять раз меньшая масса расплава (140 кг) и соответственно выделяется в 10 раз меньше кинетической энергии. Это, даже при консервативном значении коэффициента конверсии 30%, приводит к безопасным для целостности корпуса реактора значениям кинетической энергии ~ 0,6 ГДж. При определении выделяющейся кинетической энергии для варианта «уровень воды 3м, диаметр струи 10 см» с более реалистическими значениями коэффициента конверсии и рассчитанной кодом VAPEX-P массой жидких капель, находящихся в области малых паросодержаний, получаются еще меньшие значения (0,6-6 МДж).
Таблица 3 – Сводные результаты расчёта внутрикорпусного парового взрыва Уровень воды 1 м Уровень воды 2 м Уровень воды 3 м Диаметр струи, см 10 20 40 10 20 40 10 20 Полная масса 107,5 416,5 1596 110,5 426 1600 117 437,5 расплава, кг Кин.энергия (=30%), 49 189 720 50 193 730 54 200 МДж Кин.энергия (=3%), 4,9 18,9 72 5 19,3 73 5,4 20 75, МДж Кин.энергия(=0,3%), 0,49 1,89 7,2 0,5 1,93 7,3 0,54 2 7, МДж Масса капель 9 20,5 59 35,5 100 264 83 253,5 расплава, кг Кин.энергия (=30%), 4 9 26,5 16 45 119 37 114 МДж Кин.энергия (=3%), 0,4 0,9 2,65 1,6 4,5 11,9 3,7 11,4 33, МДж Кин.энергия(=0,3%), 0,04 0,09 0,265 0,16 0,45 1,19 0,37 1,14 3, МДж Масса жидких капель расплава в 6 15 45 6 22,5 81 8 42 пароводяной смеси с 0,6, кг Кин.энергия (=30%), 3 7 20 3 10 36 4 19 МДж Кин.энергия (=3%), 0,3 0,7 2 0,3 1 3,6 0,4 1,9 6, МДж Кин.энергия(=0,3%), 0,03 0,07 0,2 0,03 0,1 0,36 0,04 0,19 0, МДж Таким образом, детальный расчет процесса перемешивания расплава с охладителем выявил, что лишь небольшая часть капель расплава находится в области с достаточно большим содержанием воды. Это приводит к существенному снижению взрывного потенциала системы кориум–охладитель.
Во втором разделе этой главы были определены динамические нагрузки при термическом взаимодействии кориума с водой в шахте реактора ВВЭР 1000. Рассматривается следующая задача: кориум вытекает из разрушенного корпуса реактора в шахту, частично заполненную водой. Были выполнены два расчета, соответствующие различным начальным уровням воды в шахте реактора. Геометрические и теплофизические параметры расчётов приведены в таблице 4.
Таблица 4 – Основные параметры расчётов Температура кориума, К Состав кориума 80%UO2 + 20%ZrO Высота, с которой вытекает струя кориума, м Диаметр струи кориума, см Давление, МПа 0, Недогрев воды, К Уровень воды, м 1;
Диаметр шахты, м 5, Проведенные расчеты термического взаимодействия кориума с водой в шахте реактора показали, что оно может представлять опасность для целостности шахты.
При наличии достаточного количества воды в шахте реактора в результате распада струи кориума, истекающей из корпуса, и смешения с водой образуется взрывоопасная смесь высокотемпературных капель кориума с водой. Под действием инициирующего события возникает волна термической детонации, проводящая к динамическому нагружению стенок шахты до значений, при которых возможно разрушение шахты.
Сопоставительные расчеты при различных начальных уровнях воды в шахте выявили, что наиболее опасным с точки зрения воздействия на шахту будет термическое взаимодействие при высоком уровне воды. Импульсы давления на боковые стенки отличаются более чем на порядок (180 кПа·с уровень 3 м, 13 кПа·с - начальный уровень 1 м). Следует отметить, что воздействие на основание шахты в обоих случаях примерно одинаковое - около 400 кПа·с. Это связано с тем, что волны давления, генерируемые из зоны взрыва, проходят вниз, практически не ослабляясь, в то время как их распространение в поперечном направлении существенно ослабляется для малых уровней воды из-за разгрузки на поверхности раздела вода-водяной пар.
Таким образом, выполненный анализ показал, что при проливе расплавленного кориума в шахту реактора возникающий паровой взрыв оказывает значительные нагрузки на стенки, создавая угрозу ее целостности.
В третьем разделе этой главы были определены динамические нагрузки при термическом взаимодействии кориума с натрием в ходе тяжёлой аварии на АЭС с БН-600. Рассматривалась авария типа ULOF (Unprotected Loss Of Flow) авария с потерей расхода теплоносителя. Данная авария реализуется при полной потере системного и автономного электроснабжения с отказом аварийной защиты реактора.
Для расчёта термического взаимодействия кориума с натрием в ходе тяжёлой аварии на АЭС с БН-600 использовались следующие исходные данные:
объем газа (аргона) в газовой полости реактора 45 м3;
средняя температура натрия в верхней смесительной камере реактора 800C;
температура натрия в зоне взаимодействия с топливом 900C;
начальное давление газовой подушки над свободным уровнем реактора 0.2 МПа;
масса расплавленного топлива 3500 кг;
температура топлива 3600C.
С помощью кода CSI моделировался выход кориума в верхнюю камеру реактора БН-600 и его последующее взаимодействие с натрием. Результаты валидации кода CSI показали, что наилучшее совпадение с экспериментальными данными получается при использовании отношения массы кориума к массе натрия в зоне взаимодействия в диапазоне 0,83 - 1, и использовании модели плёночного кипения, предложенная L. Caldarola.
Поэтому при расчётах термического взаимодействия в верхней камере реактора БН-600 также использовалась эта модель плёночного кипения, а значение 1 использовались с целью варьировалось в пределах 0,8 – 1,2 (значения получения консервативных результатов).
= 1,2 максимальное В наиболее консервативном расчёте, при значение давления в натрии составило 20 бар, а давление в газе – 11 бар.
Рассчитанные по коду CSI значения коэффициентов конверсии находились в диапазоне 0,16 – 0,18 %.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ 1. Термическое взаимодействие расплавленных материалов активной зоны (кориум) с теплоносителем, возникающее в ходе тяжёлой аварии на АЭС с ВЭЭР и АЭС с БН, может приводить к существенным динамическим воздействиям на окружающие конструкции и создавать угрозу их целостности.
2. Существенное различие между типичными сценариями протеканий тяжёлых аварий на АЭС с ВВЭР и АЭС с БН приводит к различным начальным конфигурациям системы «кориум-теплоноситель», определяющим последующий характер протекания процесса термического взаимодействия кориума с водяным и с натриевым теплоносителями и связанные с ним динамические нагрузки на конструкции.
3. На основе подходов механики многофазных сред сформулированы математические модели, описывающие пространственное перемешивание кориума с водой и распространения волны термической детонации с её последующим расширением.
4. Разработана интегральная модель термического взаимодействия кориума с натрием. Предполагается, что происходит мгновенное перемешивание фрагментов кориума с натрием. При этом процесс термического взаимодействия происходит через две последовательные фазы.
На первой фазе происходит прямой контакт натрия с частицами кориума и нагрев натрия. Как только натрий достигает температуры насыщения наступает вторая фаза процесса, характеризующаяся возникновением плёночного кипения натрия на частицах кориума и резким расширением зоны взаимодействия, приводящим к возникновению динамических нагрузок на окружающие конструкции.
5. Выполнены валидационные расчёты кодом VAPEX-P экспериментов FARO-L14, FARO-L24, FARO-L28, в ходе разработки модуля VAPEX-M для интегрального системного одномерного кода СОКРАТ. Код VAPEX-P адекватно воспроизводит основные черты термического взаимодействия кориума с водой. Основным параметром, характеризующим интенсивность протекания процесса термического взаимодействия, является рост давления в системе. Код VAPEX-P хорошо качественно и количественно описывает временные зависимости давления в системе.
6. Выполнены валидационные расчёты кодом CSI экспериментов на установках THINA и ПЛУТОН. Эти валидационные расчёты показали, что, несмотря на интегральный характер математической модели кода CSI, этот код позволяет адекватно воспроизводить основные параметры термического взаимодействия кориума с натрием.
7. С помощью кода VAPEX-P выполнены расчёты динамических нагрузок при термическом взаимодействии кориума с теплоносителем внутри корпуса реактора в ходе тяжёлой аварии на АЭС с ВВЭР-1000. На основе анализа результатов расчёта показано, что внутрикорпусное термическое взаимодействие не представляет угрозу для целостности корпуса реактора.
8. С помощью кодов VAPEX-P и VAPEX-D были выполнены расчёты динамических нагрузок при термическом взаимодействии кориума с теплоносителем в шахте реактора в ходе тяжёлой аварии на АЭС с ВВЭР-1000.
Установлено, что динамические нагрузки на шахту реактора существенно зависят от уровня воды в шахте. Показано, что для обеспечения целостности шахты реактора необходимо поддерживать низкий уровень воды в шахте (около 1 м).
9. С помощью кода CSI были выполнены расчёты термического взаимодействия кориума с натриевым теплоносителем в ходе тяжёлой аварии типа ULOF в реакторе БН-600. Получены оценки максимальных значений давления натрия в зоне взаимодействия и давления аргона в газовой полости при различных соотношениях массы кориума к массе натрия в зоне взаимодействия.
СПИСОК РАБОТ, ОПУБЛИКОВАННЫХ ПО ТЕМЕ ДИССЕРТАЦИИ 1. Мелихов В., Мелихов О., Якуш С., Ртищев Н. Валидация модели термического взаимодействия топлива с теплоносителем при тяжёлой аварии на АЭС (публикация на английском языке)// Science and Technology of Nuclear Installations, vol. 2011, Article ID 560157, 11 pages, 2011. doi:10.1155/2011/560157.
2. Мелихов В.И., Мелихов О.И., Парфенов Ю.В., Ртищев Н.А., Боровкова Е.М. Оценка ударного воздействия на корпус реактора вследствие внутрикорпусного парового взрыва // Вестник МЭИ – 2012. – №2. С.50–57.
3. Мелихов В.И., Ртищев Н.А. Разработка модели термического взаимодействия расплавленного кориума с натриевым теплоносителем. // Приложение к журналу Физическое образование в вузах. Т.18, №1, 2012. Труды конференции-конкурса молодых физиков. Москва, 13 февраля 2012 г. П21.
4. Тарасов А.Е., Ртищев Н.А. Мелихов О.И., Мелихов В.И. Оценка ударных воздействий на конструкции при паровом взрыве. Сборник трудов XIII научной школы молодых учёных ИБРАЭ РАН, проходившей 26-27 апреля 2012 г. – (Препринт / Институт проблем безопасного развития атомной энергетики РАН, апрель 2012, № IBRAE-2012-02). – М. ИБРАЭ РАН, 2012. – С.146.
5. Блинков В.Н., Мелихов В.И., Мелихов О.И., Давыдов М.В., Парфёнов Ю.В., Ртищев Н.А., Тарасов А.Е., Гудеменко Д.В., Климов П.С. Определение динамических нагрузок на контейнмент при внекорпусном паровом взрыве на АЭС с ВВЭР // Фундаментальные исследования, № 9 (часть 4), 2012, С.
889-893.
6. Ртищев Н.А., Мелихов В.И. Разработка модели термического взаимодействия кориума с натрием // Восемнадцатая международная научно-техническая конференция студентов и аспирантов. 1 - 2 марта 2012 г., Москва. Тезисы докладов, 2012, том 4, МЭИ, С.25.
7. Мелихов О.И., Мелихов В.И., Ртищев Н.А. Разработка модели термического взаимодействия расплавленного топлива с натриевым теплоносителем. Теплофизические экспериментальные и расчётно теоретические исследования в обоснование характеристик и безопасности ядерных реакторов на быстрых нейтронах (Теплофизика-2012) / Сборник тезисов докладов на научно технической конференции «Теплофизика-2012». Обнинск. ГНЦ РФ ФЭИ. 2012. C.55.
8. Мелихов В.И., Мелихов О.И., Парфенов Ю.В., Ртищев Н.А., Тарасов А.Е. Разработка и валидация математической модели для оценки параметров термического взаимодействия расплава с натриевым теплоносителем // Фундаментальные исследования, №4 (часть 5), 2013, С. 1079-1083.
9. Ртищев Н.А., Тарасов А.Е., Мелихов В.И., Мелихов О.И.
Разработка модели термического взаимодействия расплавленного топлива с натриевым теплоносителем // Девятнадцатая международная научно-техническая конференция студентов и аспирантов. 28 февраля - 1 марта 2013 г., Москва. Тезисы докладов, 2013, том 4, МЭИ, С.18.