Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

DEVELOPMENT AND VALIDATION OFMODEL FOR ASSESSMENTOF CORIUM – SODIUM COOLANT THERMAL INTERACTION PARAMETERS

Melikhov V.I. 1 Melikhov O.I. 1 Parfenov Y.V. 1 Rtishchev N.A. 1 Tarasov A.E. 1
1 National Research University «Moscow Power Engineering Institute»
The paper is devoted to development of the model for assessment of melt-liquid sodium thermal interaction parameters and to its validation on experimental data FARO [8] and THINA [10]. During the liquid metal fast breeder reactor severe accident, direct contact of liquid sodium with melt is possible. As a result – their thermal interaction occurs. Under certain conditions, this interaction has explosive character, and may be dangerous for reactor integrity. Brief description of melt–sodium coolant thermal interaction model, based on mass, momentum and energy conservation laws, is presented in the paper. This model has been validated for two types of interaction: melt injection from below into sodium pool (THINA), and from above (FARO). It has been acceptable agreement of calculations results with experimental ones: interaction zone pressure, and energy conversion ratios.
thermal interaction
melt
sodium coolant
mathematical model
validation
1. Buksha Yu.K., Bagdasarov Yu.E., Kuznetsov I.A. Issledovanie vzaimodejstvija rasplavlennogo topliva s natriem v aktivnoj zone bystrogoreaktora–AtomayaEnegiya. 1976. T. 41, Vyp. 1, рр. 9–14.
2. Kirillov P.L., Jurev Yu.S., Bobkov V.P. Spravochnik po teplogidravlicheskim raschetam. M., Energoatomizdat, 1990. 360 p.
3. Tarasov A.E., Rtishchev N.A. Ocenka udarnyh vozdejstvij na konstrukcii pri parovom vzryve. Sborniktrudov XIII nauchnoj shkoly molodyhuch jonyh IBRAE RAN.– (Preprint / Institut problem bezopasnogorazvitijaatomnojjenergetiki RAN, aprel 2012, no. IBRAE-2012-02). M. IBRAE RAN, 2012. 146 p.
4. Caldalora L.A. Theoretical Model for the Molten Fuel-Sodium Interaction in a Nuclear Fast Reactor – Nuclear Engineering and Design, Vol. 22, 1972, рр. 175–211.
5. CronenbergA.W., Benz R.Vapourexplosion phenomena with respect to nuclear reactor safety assessment.– Advance in Nuclear Science and Technology, Editors: J. Lewins and M. Becker, Vol. 12, 1980, рр. 247–321.
6. Experimental Facilities for Sodium Fast Reactor Safety Studies. Task Group on Advanced Reactors Experimental Facilities (TAREF). NEA/CSNI/R(2010)12. OECD 2011.
7. Farahat M.M.K., Eggen Donald T. Pool Boiling in Subcooled Sodium at Atmospheric Pressure – Nuclear Science and Engineering, Vol. 53, 1974, рр. 240–253.
8. Kaiser A., Huber F., Wilhelm D. Experiments on the behaviour of a hot melt injected into sodium. Experimental Thermal and Fluid Science Vol. 18, 1998, рр. 48–69.
9. Kolev N.I. Film boiling on vertical plates and spheres // Experimental Thermal and Fluid Science, Vol. 18, 1998, рр. 97–115.
10. Magallon D., Hohmann H., Schins H. Pouring of 100-kg-Scale Molten UO2 into Sodium. Nuclear Technology Vol.98, 1992, рр. 79–90.
11. Witte L.G., Cox J.E. Thermal explosion hazard–Advances in Nuclear Science and Technology, Vol. 7, 1973,рр. 329–364.
12. Zyszkowski W. Thermal Explosion Hazards in (Fast) Nuclear Reactors – Atomic Energy Review. Vol. 16, 1978, no. 1, pp. 3–86.

Термическое взаимодействие высокотемпературного расплавленного топлива с охладителем может иметь место в реакторах с натриевым теплоносителем при реактивностных авариях с быстрым ростом мощности и при авариях с нарушением охлаждения активной зоны (аварии типа ULOF). Происходящее при этом взаимодействии высвобождение большого количества механической энергии представляет собой угрозу для безопасности АЭС и в связи с этим данное явление является предметом многочисленных исследований [5, 6, 12]. В настоящей работе предложена математическая модель, позволяющая оценивать параметры термического взаимодействия расплава с жидким натрием. Эта модель была валидирована на опытных данных, полученных на установке THINA (Германия), где расплав подавался в нижнюю часть экспериментального сосуда с натрием, и на установке FARO (Италия), где расплав подавался сверху (рис. 1). Было получено хорошее совпадение расчётных и экспериментальных данных как по давлению в зоне взаимодействия, так и по коэффициентам конверсии.

Цель работы: разработка и валидация модели термического взаимодействия расплава с натриевым теплоносителем.

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

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

1) интенсивным тепловым взаимодействием между расплавом и охладителем за счёт прямого контакта между ними;

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

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

pic_20.tif

Рис. 1. Экспериментальные установки THINA и FARO

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

Вторая фаза процесса термического взаимодействия начинается с момента закипания натрия. С этого момента теплопередача характеризуется:

1) ухудшением теплопередачи между расплавом и охладителем из-за образования плёнки пара натрия на внешней поверхности фрагментов расплава;

2) большой площадью поверхности взаимодействия.

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

Описанная модель основана на [1, 4, 11]. Дополнительно к этим исследованиям в настоящей работе для расчёта теплообмена частиц расплава при плёночном кипении были использованы зависимости [2, 4, 7, 9]. Стоит отметить, что в моделируемых экспериментах газовая полость над уровнем натрия была заполнена аргоном. Для учёта процессов в газовой полости математическая модель была дополнена уравнениями сохранения массы и энергии для аргона и уравнением теплообмена газа со стенками сосуда. Подробное описание математической модели изложено в [3].

Результаты исследования и их обсуждение

Расчёт эксперимента THINA

Установка THINA (см. рис. 1) состоит из двух основных частей: экспериментального сосуда и системы для приготовления и инжекции расплава. Около половины экспериментального сосуда заполнено жидким натрием, над поверхностью которого находится газ (аргон) для компенсации объёмных расширений. Малый диаметр сосуда (0,3 м) был выбран для возможности использования рентгеновского излучения с целью диагностики исследуемых процессов. Около 5 кг расплава подготавливалось в тигле с помощью экзотермической реакции алюминий/оксид железа согласно уравнению:

2Al + Fe2O3 = 2Fe + Al2O3 + 4117 Дж/г.

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

Предложенная модель была валидирована на всех четырёх экспериментах [10]. В данной работе приводятся результаты расчётов эксперимента TH561, основные параметры которого приведены в таблице:

Основные параметры экспериментов

 

TH 561 (THINA)

T2 (FARO)

Температура расплава, К

2780

3273

Время подачи расплава, с

0,190

-

Состав расплава

4 кг Al2O3 + 1,5 кг стали

112 кг UO2 + 28кг ZrO2

Масса натрия, кг

155

130

Температура натрия, К

790

673

Высота столба натрия, м

2,67

2,5

Начальное давление аргона, бар

1,1

1,3

С помощью разработанной математической модели были получены следующие результаты. Рассчитанный пик давления в зоне взаимодействия составил 7,4 бар, что хорошо согласуется с экспериментально полученным значением в 8,0 бар. Временные зависимости основных параметров, полученных в эксперименте и в расчёте, приведены на рис. 2. Коэффициент конверсии, полученный в эксперименте для данного теста, составил 0,13 %, что несколько выше полученного значения с помощью разработанной математической модели – 0,07 %. На рис. 2, 3 наблюдается качественное соответствие расчётных и экспериментальных данных. Различия вызваны принятой схематизацией исследуемого процесса (интегральный подход для формулировки законов сохранения) и неполнотой информации о параметрах экспериментальной установки и теплофизических свойств высокотемпературного расплава. Расчёты трёх других экспериментов на установке THINA [10] дали аналогичные результаты.

pic_21.tif

Рис. 2. Временные зависимости давления натрия в зоне взаимодействия, давления аргона

pic_22.tif

Рис. 3. Временная зависимость высоты зоны взаимодействия

Расчёт эксперимента FARO

Установка FARO [10] состоит из экспериментального сосуда TERMOS и системы подачи кориума.

TERMOS представляет собой сосуд высокого давления, внешним диаметром 800 мм и высотой 4,0 м. Контейнер с натрием находится внутри сосуда TERMOS и имеет диаметр 470 мм. Экспериментальный сосуд подсоединен к устройству для подготовки расплава через канал и три межсекционных клапана. Объём инертного газа, находящегося над уровнем натрия, составлял 0,15 м3. Температура натрия равна 673,15 К, уровень жидкого натрия – 2,5 м. Начальные условия эксперимента приведены в таблице.

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

Выводы

Предложена модель термического взаимодействия расплава с натрием, основанная на законах сохранения массы и энергии для зоны взаимодействия и импульса для столба холодного натрия. Данная модель была валидирована для двух типов взаимодействия: при подаче расплава в нижнюю часть экспериментального сосуда (THINA) и в верхнюю (FARO). Сравнение проводилось по нескольким параметрам: максимальное давление в натрии и аргоне, высоты подъёма столба натрия и коэффициентов конверсии. Было получено приемлемое соответствие расчётных и экспериментальных результатов.

Работа была выполнена при финансировании ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009–2013 гг. (госконтракт № 14.В37.21.0151), НИР НИУ«МЭИ» «Исследование теплофизических процессов в перспективном энергетическом оборудовании» рег. № 7.3163.2011, проекта РФФИ № 11-08-00410_а.

Рецензенты:

Семенов В.Н., д.ф.-м.н., заведующий отделом моделирования АЭС с реакторными установками нового поколения, ФГБУН «Институт проблем безопасного развития атомной энергии» РАН, г. Москва;

Якуш С.Е., д.ф.-м.н., ведущий научный сотрудник лаборатории термогазодинамики и горения, ФГБУН «Институт проблем механики им. А.Ю. Ишлинского» РАН, г. Москва.

Работа поступила в редакцию 12.03.2013.