автореферат диссертации по энергетике, 05.14.03, диссертация на тему:Массоперенос продуктов коррозии с учетом химического взаимодействия в системе натрий - конструкционный материал - примеси

На правах рукописи УДК 621.039.534.63

Кондратьев Александр Сергеевич

МАССОПЕРЕНОС ПРОДУКТОВ КОРРОЗИИ С УЧЕТОМ ХИМИЧЕСКОГО ВЗАИМОДЕЙСТВИЯ В СИСТЕМЕ НАТРИЙ - КОНСТРУКЦИОННЫЙ МАТЕРИАЛ - ПРИМЕСИ

Специальность 05.14.03 - Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

484433'

Обнинск-2011

4844337

Работа выполнена в Государственном научном центре Российской Федерации -Физико-энергетическом институте имени А. И. Лейпунского (ГНЦ РФ-ФЭИ)

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

доктор технических наук Алексеев Виктор Васильевич

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

доктор технических наук, профессор Арнольдов Михаил Николаевич

кандидат технических наук Ромадова Елена Леонардовна

Ведущая организация: Открытое акционерное общество «Опытное конструкторское бюро машиностроения им. И.И. Африкантова» (ОАО «ОКБМ АФРИКАНТОВ»), г. Нижний Новгород.

Защита состоится 25 |и/У1<Х 2011 г. в 10 часов на заседании диссертационного совета Д 201.003.01 при ГНЦ РФ-ФЭИ в конференц-зале по адресу: 249033, г. Обнинск Калужской области, пл. Бондаренко, 1.

С диссертацией можно ознакомиться в научно-технической библиотеке ГНЦ РФ-ФЭИ. Автореферат разослан 11 г.

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

доктор технических наук

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

Цель диссертационной работы состоит в создании метода расчета массопе-

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

Задачи исследования:

» представить математическое описание физико-химических процессов массопе-реноса продуктов коррозии с учетом химического взаимодействия в системе натрий - конструкционный материал - примеси;

• создать расчетную программу для моделирования в одномерном приближении массопереноса конструкционных материалов в натриевых контурах для различных условий эксплуатации ЯЭУ;

• получить физико-химические константы, характеризующие массоперенос в системе натрий - конструкционный материал - примеси;

• разработать трехмерную модель сопряженного тепломассопереноса (переноса примесей);

• на основании выполненных расчетов проанализировать и обобщить полученные данные, на основе которых сформулировать научно обоснованные рекомендации для ЯЭУ.

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

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

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

• получены уточненные физико-химические константы, характеризующие массоперенос в системе натрий - конструкционный материал - примеси. Получен ряд простых формул для расчета термодинамических и транспортных теплофи-зических свойств натрия;

• разработана новая трехмерная модель сопряженного тепломассопереноса (переноса примесей).

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

• на детальном анализе процессов переноса продуктов коррозии в циркуляционном натриевом контуре и на сопоставлении результатов анализа с известными экспериментальными данными;

• на реализации методов трехмерного моделирования для расчетов процессов массопереноса, результаты которых находятся в хорошем соответствии с литературными данными;

• на системном подходе к проведенным исследованиям, в ходе которых один и тот же результат получен различными методами;

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

Практическая значимость диссертационной работы состоит в том, что:

• разработанная методика и вычислительная программа расчета позволяет определять массоперенос продуктов коррозии в натриевых контурах ЯЭУ в реальных условиях с учетом химического взаимодействия в системе натрий -конструкционный материал - примеси;

• определены экстремальные потоки примесей на стенки каналов для различных параметров натриевого теплоносителя первого контура реактора БН-600 в условно поставарнйный период, а также для высокотемпературных ЯЭУ с натриевым теплоносителем;

• на основе анализа полученных данных по массопереносу в натриевых контурах, сформулированы научно обоснованные рекомендации для ЯЭУ.

Основные положения, выносимые на защиту:

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

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

• методика расчета массопереноса примесей в холодной ловушке в трехмерной динамической постановке.

Апробация работы: основные положения диссертационной работы докладывались и обсуждались на Всероссийских конференциях и школах: "Реакторы на быстрых нейтронах" (Обнинск, 2009); "Теплофизика-09" (Обнинск, 2009); V Студенческая конференция «Студенчество - будущее атомной энергетики» (Обнинск, 2009); "Реакторы на быстрых нейтронах" (Обнинск, 2010).

На международных конференциях и семинарах: «Безопасность АЭС и подготовка кадров», Обнинск, Россия, 29 сентября - 2 октября, 2009; International Conference on Fast Reactors and Related Fuel Cycles: Challenges and Opportunities FR09, Kyoto, Japan, 07 -11 December 2009.

Публикации: основное содержание диссертации изложено в трех статьях в реферируемом журнале, в двух препринтах ГНЦ РФ-ФЭИ и двенадцати публикациях в сборниках тезисов докладов и трудах конференций.

Структура и объем диссертации: диссертация состоит из введения, четырех глав, заключения, списка литературы. Материал изложен на 130 страницах, содержит 41 рисунок, 19 таблиц, список литературы из 72 наименований.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Массопереиос хрома с учетом его химического взаимодействия с кислородом

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

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

В контурах из хромоникелевой нержавеющей стали при определенном содержании кислорода в натрии и температуре образуются наряду с оксидом натрия Na20 соединения кислорода и натрия с элементами конструкционных материалов Na4Fe03 и NaCr02.

Известна зависимость для определения границы стабильности NaCr02 на поверхности стали в виде Ige = -3326/Т+4,51 -0,5lg(a(C7j), где с - концентрация кислорода в натрии, млн"1; a(Cr) - активность хрома в стали; Т - температура, К. В соответствии с этой зависимостью для стали AISI 316 (a(Cr) = 0,4) соединение NaCrCb стабильно не выше температуры 600°С при с - 8 млн"'.

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

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

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

Cr + 2Na20 = 3Na + NaCr02. (1)

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

струкционный материал - кислород в значительном количестве (для хромоникеле-вых сталей порядка 10-20 мае. %).

Константа равновесия реакции (1) выражается через активности соответствующих компонентов

, = а(КаСг02)-а(]Ма)3 а(Сг)а(Ыа20)2

Кинетическое уравнение этой реакции, если рассматривать изменение концентрации ИаСгОг, выражается зависимостью

(А^аСЮ^ =^^с(Сг)с(1Ча20)2_^^с(ЫаС1.02)с(Ка)з. (3)

V л /р

Индекс «р» здесь указывает на изменение концентрации сложного оксида за счет реакции (1).

Константа равновесия реакции (1) К = К~*/К*~. С учетом конвективного переноса растворенного сложного оксида в натрии получим изменение его концентрации в теплоносителе в мольных долях (активность натрия равна единице)

ММаСг02) = ^^с(Сг)с(Ма20)2 - + 4Хпл., (4)

где Ок - диаметр канала, м;

коэффициент массоотдачи сложного оксида от стенки канала, м/с; индекс «и'» характеризует пристеночные значения концентраций.

Изменение концентрации хрома в натрии в дифференциальном виде с учетом его конвективного переноса в канале получим также в мольных долях

с1с{Сг) = Г</с(НаСЮ2П +4х с(Сг),„ - с(Сг) (5)

с1х \ с/т )р Ок

Аналогично для изменения концентрации кислорода (в форме оксида) в натрии получено выражение

Л I ¿х )р 0 Вк

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

с(№20)„,=С(№20). (7)

Пристеночная концентрация хрома с(Сг)„, рассчитывается на основании баланса потоков хрома:

1) поступающего из стали через слой сложного оксида на ее поверхности в натрий

J _ ^ух(аст~а\у) (8)

* 6

и 2) выносимого в поток натрия за счет конвективного массообмена

^=Хх(аСг\,-с(Сг)), (9)

где Кт - коэффициент проницаемости хрома в двойном оксиде, кг/(м-с); 8 - толщина слоя двойного оксида хрома, м; аС1 - активность хрома в стали;

а„. - активность хрома в натрии в пристеночной области.

Совместно решая (8) и (9), с учетом соотношения

с(Сг),=С(Сг)5аи„ (10)

получим

с(Сг) ^(Сг^+МСг)^ (11)

" д+рс(С г)у

где (12)

о

Р-ХхРыа —. (13)

^Ыа

цСг и |Лма - молярные массы хрома и натрия соответственно, г/моль; с(Сг), - концентрация насыщения хрома в натрии, м. д.

Концентрация двойного оксида хрома в пристеночном слое определяется на основании уравнения (2), записанного для соответствующих условий

^аС^^^-^аСЮ,),^^-)2 (И)

с(Сг), [с(Иа2ои

где с(На20)1 - концентрация насыщения оксида натрия в натрии, м. д.

Коэффициент х определяется исходя из известной критериальной зависимости тепломассообмена

№/¿=5 + 0,025-Ре%*, (15)

о и °к

где

X/= ' ■ (16)

ик

где Dj - коэффициент диффузии примеси в теплоносителе (здесь индекс г характеризует вид примеси).

Константа равновесия (2)

к =ехр

АО0

ят

(17)

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

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

Задание исходных данных и определение констант

Расчеты массопереноса хрома проводятся для первого контура реактора БН-600, схема которого представлена на рисунке 1. Для построения расчетной модели:

- Выделен гидравлический тракт с циркуляцией основной массы теплоносителя (не менее 94 % от максимального расхода по контуру).

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

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

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

Рисунок 1 - Схема гидравлического тракта первого контура реактора (БН-600): 1 - активная зона (0,52 с)*; 2 - зона воспроизводства (1,89 с); 3 - хранилище (9,12 с);

4 - верхняя камера (5,28 с) + смесительная камера (6,05 с) + коллектор ПТО (0,63 с);

5 - ПТО (6,63 с); 6 - коллектор ПТО (0,63 с) + сливная камера ПТО (18,03 с) + всасывающая камера насоса (7,34 с); 7 - напорный трубопровод (1,10 с); 8 - напорный коллектор (1,59 с).

*'В скобках указано время пребывания теплоносителя в соответствующем участке

- Определено время прохождения теплоносителя в каждом участке контура, геометрические параметры, гидродинамические и температурные режимы на каждом участке.

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

Значения констант сведены в таблицу 1.

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

Считается, что вся поверхность контура, омываемая натрием, покрыта пленкой двойного оксида хрома №СЮ2 толщиной 2,5 мкм. Данное значение принято, исходя из опыта работы с оборудованием натриевых контуров в ГНЦ РФ-ФЭИ.

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

Таблица 1 - Значения констант для расчета массопереноса хрома в натриевом контуре

Параметры Обозначение Размерность Значение параметра или расчетная зависимость

Константа скорости реакции образования двойного оксида хрома КГ> 1/(м.дЛс) 108-ехр(0,477/10-(Т-823))

Концентрация насыщения двойного оксида хрома в натрии с(№СЮ2)л м.д. . 5-Ю"'2

Концентрация насыщения хрома в натрии с(Сг)5 млн"' ехр(16,2-20746/Т)

Активность хрома в стали с'ст(Сг) - 0,6

Коэффициент диффузии хрома в натрии А м2/с Ю-8

Изобарно-изотермический потенциал реакции образования двойного оксида хрома ДО0 Дж/моль -35539-76,9 Т

Коэффициент диффузии хрома в двойном оксиде Д» м2/с ю-'8

Константа скорости кристаллизации или селективного растворения Р* м/с 310"\

Коэффициент, учитывающий шероховатость поверхности Кг - 2

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

Аналогично представленной выше модели рассматривается массоперенос железа в натрии с учетом образования системы натрий-кислород-железо и протекания следующей реакции в теплоносителе

Ре + ЗИа20 = 2ТЧа + №4Ре03. (18)

Полагаем, что кислород в натрии находится в форме На20. Другими реакциями кислорода пренебрегаем. Железо присутствует в системе натрий-кислород-конструкционный материал в значительном количестве (до десятков процентов в конструкционном материале). Соединение №4Ре03 наиболее устойчивое в системе натрий-кислород-железо.

В главе 3 представлены результаты расчетов переноса хрома и железа для первого контура реактора БН-600, в том числе при гипотетическом увеличении температуры и концентрации кислорода в натрии.

Расчеты массопереноса хрома проведены для первого контура реактора БН-600 с концентрациями кислорода в натрии: 2, 10, 50, 100, 500 млн"1.

В таблице 2 приведены результаты расчета потока двойного оксида хрома (ЫаСгСЬ) на стенку активной зоны в зависимости от концентрации кислорода в натрии при различной температуре.

Таблица 2 - Поток двойного оксида хрома (ИаСгСЬ) на стенку, мкм/год,

активная зона

Концентрация, млн-1 550°С 650°С 750°С 850°С 950°С

2 0,641 7,95 43,1 73,4 127

10 3,54 40,7 83 112 532

50 29,0 83,5 115 355 3820

100 45,7 99,7 128 821 5890

500 83,9 124 238 - 5800 15400

Из данной таблицы видно, что при увеличении концентрации кислорода в натрии с 2 до 10 млн'1, скорость роста отложений в виде двойного оксида хрома (NaCr03) в экстремальной точке активной зоны (на выходе) возрастает в 6 раз, что при номинальных режимах (550 °С) не опасно для эксплуатации. При гипотетиче-

ском увеличении температуры теплоносителя на выходе из активной зоны на 100 °С, при концентрации кислорода в натрии 50 млн-1, скорость роста отложений составляет 83,5 мкм в год.

Относительный шаг расположения твэлов в сборке в БН-600 я/йМ,17 при диаметре твэла ¿N>,9 мм, при этом зазор между твэлами в узком сечении составляет 1173 мкм. Из таблицы 2 видно, что расчетная величина годовых отложений заметно превышает величину межтвэльных зазоров при некоторых условиях (при значительных превышениях концентрации кислорода и температуры выше номинальных). Кроме того, с учетом двухстороннего роста отложений в зазорах табличные данные необходимо удвоить.

На рисунках 2-5 представлено распределение потоков хрома и двойного оксида хрома ЫаСЮ2 в виде взвеси на стенку по длине контура (характеризуется временем прохождения теплоносителя по контуру).

Рисунок 2 - Распределение потоков хрома по длине первого контура реактора БН-600 в номинальном режиме работы

1.00Е-11

-1.00Е-И

■ -З.ООЕ-11

-5.00Е-11

— -7.00Е-11

-9,00Е-П -1

Ч 2! 2 2.с .4 2Й.6 28.5 2 :< 2 25 .......

/

---------

Время, с

Рисунок 3 - Распределение потоков хрома по длине первого контура реактора БН-600 в номинальном режиме работы (фрагмент рис. 2)

1,80Е-)0

и

<ч Е 1.60Е-10

X

Л 1.40Е-10

И Й

е 1.20Е-10

я

X

и 1.00Е-10

аэ

33 и 8.00Е-И

я а 6.00Е-11

гч

9

и 4.00Е-11

Я

X

* 2 2.00Е-11

Ь

С

О.ООЕ+ОО

15 20 25 30 Время от точки отсчета, с

Рисунок 4 - Распределение потоков двойного оксида хрома (№СгСЬ) по длине первого контура реактора БН-600 в номинальном режиме работы

1.Я0Е-10

28 28,2 28,4 28,6 28,8 29 29,2 29,4 29,6 29,8

Время, с

Рисунок 5 - Распределение потоков двойного оксида хрома (№СЮ2) по длине первого контура реактора БН-600 в номинальном режиме работы (фрагмент рисунка 4)

-Расчет -Эксперимент

б.ООЕ-11

4.00Е-Ц

о гм 2.00Е-П

.2

м 0.00Е+00

£ -2.00Е-11

а

X -4.00Е-11

к й

й -6.00Е-11

и

й -8.00Е-11

С

-1.00Е-10

-1,20Е-1П

-1.40Е-10

1

0 1 -кГ*— 5 2 ) 2

" Г 1 5 3 4 !) 4 5

у

1

Время, с

Рисунок 6 - Распределение суммарного потока компонентов стали по длине первого контура реактора БН-600 в номинальном режиме работы

Расчетные данные получены для случая отсутствия оксидной пленки на поверхности стали. Как видно из приведенного сравнения (рис. 6), в основном по контуру данные, полученные по расчетам, отличаются от данных работы Н.Д. Краева не более чем на 5 %. В активной зоне, где скорость теплоносителя достигает 8 м/с, расчетные данные отличаются от экспериментальной формулы на 55 %. Эти различия можно объяснить тем, что экспериментальная зависимость не учитывает изменение скорости коррозии стали при скоростях теплоносителя больше 4 м/с.

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

Таблица 3 - Средний размер частиц взвеси двойного оксида хрома №Сг02

Температура, °С Концентрация кислорода, млн"' Средний размер частиц, мкм

550 2 0,04

500

650 2 0,09

500 0,2

750 2 0,2

500 0,4

На рисунке 7 представлена зависимость выхода продуктов коррозии от толщины пленки в виде двойного оксида хрома при номинальном режиме работы реактора БН-600.

Полученные данные по выходу продуктов коррозии в натрии первого контура реактора БН-600 (при толщине пленки двойного оксида хрома равной 2,5 мкм) соизмеримы с данными расчетов, полученными без учета оксидной пленки на поверхности конструкционного материала.

120.00 100,00

§

и.

ю

^ 80.00 §

г*

о о.

§ 60,00 ю

о 40.00 &

а

¿8 ю,оо

0,00 -

0123456789 10

Толщина пленки двойного оксида хрома, мкм

Рисунок 7 - Зависимость выхода продуктов коррозии от толщины пленки двойного оксида хрома при номинальном режиме работы реактора БН-600

Сравнение данных по плотности потока массы продуктов коррозии конструкционных материалов на стенки каналов показывает: при отсутствии химической реакции образования двойного оксида хрома на выходе из активной зоны плотность потока равна 5-10"11 кг/(м2-с), что в пересчете на скорость роста отложений соответствует 0,2 мкм/год, а плотность потока при учете химической реакции составляет 1,45-10"'° кг/(м2-с) или 0,6 мкм/год.

На рисунках 8-9 представлено распределение потоков железа в виде взвеси на стенку по длине контура (характеризуется временем прохождения теплоносителя по контуру).

При параметрах, соответствующих номинальному режиму работы реактора, на выходе из активной зоны поток железа в виде взвеси на стенку составляет 7,43-10~7кг/(м2'С).

8.00е-07

й 7.00Е-07 £

1 6.00Е-07

В

и

в 5.00Е-07

и

и га

3 4.00Е-07

в!

ш

2 3,0<)е-07

М

| 2.00Е-07 N

| 1,ООЕ-07 о

к

О.ООЕ+ОО

10

1

! ; I

:

: ;

1 :

( :

| |

15

20 25 30 35

Время от точки стчеяа, с

45

Рисунок 8 - Распределение потоков железа в виде взвеси на стенку по длине первого контура реактора БН-600 в номинальном режиме работы

Врем я,-с

Рисунок 9 - Распределение потоков железа в виде взвеси на стенку по длине первого контура реактора БН-600 в номинальном режиме работы (фрагмент рисунка 8)

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

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

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

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

Уравнение сохранения массы

|(р) + У(ри) = 0. (19)

Уравнение сохранения импульса

^• + У.(р1/.(7) = -У.р + У.[(р + ц7')УС7] + р^. (20)

Турбулентная вязкость получается на основе двухпараметрической модели турбулентности. При использовании стандартной к-£ -модели турбулентности ве-

личина р - коэффициента турбулентной динамической вязкости определяется по

формуле Колмогорова

(21)

Запись уравнения сохранения энергии отличается для теплоносителя и стенки для принятой сопряженной постановки задачи теплообмена. Для стенки

Для теплоносителя

с}Т

(22)

(23)

Турбулентная теплопроводность X1 определяется на основе по (21). Используется соотношение

[\ + Хт)=/

1 +

и

г \

Рг Рг7

(24)

В (24) число Прандтля Рг— — ', удельная температ7ропроводность а -

а р Ср

турбулентная температуропроводность а1 = ^

рСр

Уравнения переноса концентрации примеси записаны для размерной величины С, [кг/м3] - концентрации /-ой компоненты примеси в теплоносителе или стальных конструкциях. Для принятой сопряженной постановки задачи переноса концентрации примеси получим. Для стенки

дС,

Эт

Для теплоносителя

^+V ■ (ис,■)=V ■ [( О/ + Ьг) V • С, ]■+ .

(25)

(26)

На основе ц' определяется Эт [м2/с] - коэффициент турбулентной диффу-

зии теплоносителя. Используется соотношение

Г Л-Л

+ о7 = о/

1 +

Бс1'

(27)

Здесь число Шмидта 5с = —; V [м /с] - коэффициент кинематической вязкости натрия; 5с7 =0,9 - турбулентное число Шмидта; О/ [м2/с] - коэффициент молекулярной диффузии /-ой компоненты примеси в теплоносителе.

Для определения всех неизвестных системы, кроме описанных выше уравнений сохранения (уравнение сохранения импульса решается в виде трех скалярных уравнений для компонент скорости), необходимо иметь дополнительно еще три уравнения. Это уравнение состояния и два уравнения модели турбулентности. При учете температурной зависимости использовалась полученная с участием автора система простых формул для описания термодинамических свойств натрия.

Учет температурной зависимости предела растворимости примеси от температуры приводит к необходимости записи в уравнении переноса (26) источниково-го компонента Jj [кг/(м3 с)].

Выполнен расчет переноса примесей в холодной ловушке (ХЛ) в трехмерной динамической постановке. Сетка в расчетной области выполнена упрощенно, дно и верх ловушки выполнены плоскими, три фильтра на выходе для улавливания уходящих частиц не моделировались (рис. 10).

Результаты расчетов переноса примесей в холодной ловушке

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

Рисунок 10 - Расчетная сетка рабочей полости холодной ловушки

Таблица 4 - Исходные данные

Примесь Оксид натрия Хром

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

на входе в ХЛ 483,15 483,15

на стенках ХЛ 423,15 423,15

Скорость натрия на входе в ХЛ, м/с 2,65 2,65

Концентрация насыщения, млн"1 j Q6,2571 -2444,5/Т ехр(16,2-20746/Т)

Концентрация на входе, кг/м3 3,6-10"2 2,4-

Плотность, кг/м'' 2270 7180

На рисунке 11 представлено установившееся векторное поле скорости. Осевая скорость на входе в XJ1 задана 2,65 м/с, в то время как выходная (отрицатель' ная) скорость порядка -4 мм/с на этой шкале практически не видна. Несмотря на

I

малую величину средней скорости для всей области течения, наблюдается смешанное турбулентное вынужденное (Re~104>2300) и турбулентное свободно-конвективное (GrPr- 10|3> 109) течение теплоносителя. Температура натрия снижалась от 210 до 150 °С.

I I

23

|___ _ _ _ _ __ _ _ _ _ _ _

Рисунок 11 - Векторное поле скорости в рабочей области холодной ловушки и в зоне выхода из подающего патрубка

Сложное течение определяется турбулентной естественной конвекцией и турбулентным вынужденным течением. В нижней части при проникновении струи в зону отстоя наблюдается слабое растекание струи и обратное течение. В охлаждаемой части по высоте наблюдаются периодические структуры течения, характерные для естественной конвекции при неизотермическом течении. На охлаждаемой поверхности стенок в верхней части рабочей полости движение натрия направлено вниз против основного потока, а вдоль поверхности подводящей трубки, имеющей более высокую температуру, скорость натрия направлена вверх. В объеме отстойника образуется зона рециркуляции натрия с малой скоростью около (1-2) мм/с.

На рисунках 12-13 представлено распределение потоков оксида натрия и хрома.

Рисунок 12 - Поток оксида натрия на стенку рабочей полости холодной ловушки

Э.Ов-Л

°'°0.0 0.3 1.0 (Л 2.0 2.5 3.0 3.3 4.0 4.5 5.0 5.5 в!о 0.3

Рисунок 13 - Поток хрома на стенку рабочей полости холодной ловушки

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

Проведенные расчеты демонстрируют возможность применения модифицированного пакета программ ОрепРОАМ для решения реальных практических задач анализа сложного тепломассопереноса в оборудовании реакторов типа БН.

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

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

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

ОСНОВНЫЕ ВЫВОДЫ И РУЗУЛЬТАТЫ РАБОТЫ

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

2. В соответствии с предложенной моделью разработана расчетная программа на алгоритмическом языке «Паскаль». Выполнены расчеты и получены новые данные в обоснование безопасности натриевых контуров быстрых реакторов:

- По переносу хрома: при наличии оксидной пленки на поверхности стали для номинального температурного режима теплоносителя первого контура реактора БН-600 с концентрацией кислорода в натрии (2 млн-') выход хрома в пересчете на изменение толщины стенки канала в экстремальной точке активной зоны (на выходе) составит 0,34 мкм в год, на входе в промежуточный теплообменник -0,06 мкм в год. Увеличение толщины пленки двойного оксида хрома (МаСЮ2) в каналах на выходе из активной зоны составляет 0,64 мкм в год, на входе в промежуточный теплообменник - 0,07 мкм в год, что не приведет к изменению эксплуатационных характеристик первого контура.

- По переносу э/селеза: в условиях загрязнения поверхности первого контура двойным оксидом железа (Ыа4РеОз) для реактора БН-600 при номинальном температурном режиме теплоносителя скорость роста отложений двойного оксида железа в виде взвеси на выходе из активной зоны может достигать 9,9 мкм в сутки.

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

3. Химическое взаимодействие хрома с кислородом в натрии приводит к увеличению взвешенной фазы в натриевом контуре. При отсутствии химической реакции образования двойного оксида хрома в натрии концентрация частиц равна 4,5-10~7 кг/м3, при учете химической реакции концентрация взвеси составляет 3,77-10~3 кг/м3 (номинальный режим).

4. Получены данные по среднему размеру частиц взвеси двойного оксида хрома NaCrC>2, при температуре натрия 550 °С и концентрации кислорода в натрии 2 млн"1 средний размер частиц взвеси NaCr02 составляет 0,04 мкм, при увеличении концентрации кислорода средний размер частиц не изменяется. При температуре натрия 650 °С и концентрации кислорода в натрии 2 млн"1 средний размер частиц взвеси NaCrCb составляет 0,09 мкм, при концентрации 500 млн"1 средний размер частиц взвеси равен 0,2 мкм.

5. Сравнение данных по распределению суммарного потока компонентов стали в натрий по длине первого контура реактора БН-600 в номинальном режиме работы показывают, что в основном по контуру данные, полученные по расчетам, отличаются от известных экспериментальных данных не более, чем на 5 %. В отдельных зонах с высокой скоростью потока натрия (8 м/с) различие составляет 55 %.

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

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

данные гидродинамики и массопереноса для системы натрий - конструкционный

материал - примеси применительно к процессам осаждения примесей в холодной

ловушке при заданных температурах поверхности омываемой натрием.

СПИСОК РАБОТ, ОПУБЛИКОВАННЫХ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Казанцев A.A., Кондратьев A.C. Простые формулы для определения термодинамических свойств натрия // Изв. вузов. Ядерная энергетика. - 2008. - № 4. -С. 86-93.

2. Алексеев В.В., Козлов Ф.А., Орлова Е.А., Кондратьев A.C., Торбенкова И.Ю. Моделирование массопереноса продуктов коррозии в контурах ЯЭУ с натриевым теплоносителем / Материалы XI-й Межд. конф. «Безопасность АЭС и подготовка кадров-2009». - Обнинск, 29 сентября - 2 октября 2009. - С. 139 - 141.

3. Казанцев A.A., Кондратьев A.C., Левченко Ю.Д. Расчет CFD кодом гидравлического сопротивления прямоугольного канала при изменении расстояния между перегородками / Материалы XI-й Межд. конф. «Безопасность АЭС и подготовка кадров-2009». - Обнинск, 29 сентября - 2 октября 2009. - С. 154 - 155.

4. Алексеев В.В., Козлов Ф. А., Орлова Е.А., Кондратьев A.C., Торбенкова И.Ю. Моделирование массопереноса продуктов коррозии в контурах ЯЭУ с натриевым и свинцовым теплоносителем / Материалы Научной школы «Реакторы на быстрых нейтронах». - Обнинск, 26 - 30 октября 2009 г. - С. 1-12.

5. Алексеев В.В., Козлов Ф.А., Орлова Е.А., Кондратьев A.C., Торбенкова И.Ю. Моделирование процессов массопереноса и коррозии сталей в ядерных энергетических установках со свинцовым теплоносителем (часть 2: Разработка одномерной модели массопереноса): Препринт ФЭИ-3154. - Обнинск: ГНЦ РФ-ФЭИ, 2009 г.-27 с.

6. Алексеев В.В., Козлов Ф.А., Орлова Е. А., Кондратьев A.C., Торбенкова И.Ю. Моделирование массопереноса продуктов коррозии в контурах ЯЭУ с натриевым и свинцовым теплоносителем / Материалы Межведомственного семинара «Технология щелочных жидкометаллических теплоносителей. Теплофизика-2009». -Обнинск, 28 - 30 октября 2009. - С. 109-111.

7. Alexeev V.V., Kozlov F.A., Orlova Е. A., Kondratiev A.S., Torbenkova I.Yu. The modeling of corrosion products mass transfer in circuits of LMFDR with sodium and lead coolant // Int. Conf. on Fast Reactors and Related Fuel Cycles: Challenges and Opportunities FR09. - P. 343 - 345.

8. Варсеев Е.В., Кондратьев A.C. Демонстрация возможностей кода OpenFOAM в теплогидравлике быстрых реакторов // Материалы V Студенческой конф. физико-энергетического ф-та «Студенчество - будущее атомной энергетики». - Обнинск, 18 декабря 2009 г. - С. 7 - 8.

9. Алексеев В.В., Козлов Ф.А., Орлова Е.А., Кумаев В.Я., Кондратьев A.C., Торбенкова И.Ю. Моделирование процессов массопереноса и коррозии сталей в ядерных энергетических установках со свинцовым теплоносителем (часть 3: Уточнение предварительных расчетов. Разработка трехмерной модели): Препринт ФЭИ-3179. - Обнинск: ГНЦРФ-ФЭИ, 2010 г. -24 с.

10. Алексеев В.В., Козлов Ф.А., Орлова Е.А., Кондратьев A.C., Торбенкова И.Ю. Моделирование процессов массопереноса и коррозии сталей в ядерных энергетических установках со свинцовым теплоносителем // МНТК-2010 «Безопасность эффективность и экономика атомной энергетики». - Москва, 26 - 27 мая 2010. -С. 356-357.

11. Алексеев В.В., Кондратьев A.C. Моделирование массопереноса продуктов коррозии в контурах ЯЭУ с натриевым теплоносителем / Материалы Научной школы «Реакторы на быстрых нейтронах». Обнинск, 13-17 сентября 2010 г. - С. 42 - 43.

12. Казанцев A.A., Кондратьев A.C., Варсеев Е.В. Численное моделирование тепло-гидравлики жидких металлов с помощью открытого пакета OpenFOAM // Материалы Научной школы «Реакторы на быстрых нейтронах». - Обнинск, 13 -17 сентября 2010 г. - С. 25-26.

13. Алексеев В.В., Козлов Ф.А., Орлова Е.А., Кондратьев A.C., Торбенкова И.Ю. Моделирование массопереноса продуктов коррозии в контурах со свинцовым теплоносителем / Сб. д-дов Межведомственного семинара «Тяжелые жидкометал-лические теплоносители в быстрых реакторах. Теплофизика-2010». - Обнинск, 20-22 октября 2010.-С. 141-152.

14. Казанцев A.A., Кондратьев A.C., Левченко Ю.Д. Расчет CFD кодом гидравлического сопротивления прямоугольного канала при изменении расстояния между перегородками // Изв. вузов. Ядерная энергетика. - 2010. - № 4. - С. 193 - 198.

15. Алексеев В.В., Кондратьев A.C. Моделирование массопереноса продуктов коррозии в контурах ЯЭУ с натриевым теплоносителем // Изв. вузов. Ядерная энергетика. - 2010.-№ 4. - С. 162 -171.

Подписано к печати 16.02.2011 г. Заказ №.56. Формат 60x84 '/16. Усл. п. л. 0,9. Уч.-изд. л. 1,1. Тираж 50 экз.

Отпечатано в ОНТИ методом прямого репродуцирования с оригинала автора. 249033, Обнинск Калужской обл., пл. Бондаренко, 1 ГНЦ РФ - Физико-энергетический институт имени А.И. Лейпунского.

(

--301

Список сокращений.

Введение.

Глава 1. Состояние вопроса.

1.1. Физико-химические процессы в системе натрий - конструкционный материал -примеси.И

1.1.1. Физико-химические и коррозионные свойства натрия.

1.1.2. Уравнения коррозии для системы натрий — конструкционный материал - примеси

1.2. Массоперенос в натриевом теплоносителе.

1.2.1. Массообмен одиночной частицы.

1.2.2. Массообменные процессы в дисперсных потоках.

1.2.3. Образование и состав частиц взвесей в натрии.

1.3. Гидродинамические эффекты.

1.3.1. Гидродинамика течения жидких металлов в каналах.

1.3.2. Гидродинамические особенности поведения частиц в дисперсных потоках.

Выводы к главе 1.

Глава 2. Модели массопереноса.

2.1. Обоснование модели процесса переноса хрома.

2.1.1. Модель циркуляционного контура.

2.1.2. Задание исходных данных-и определение констант для хрома.

2.2. Обоснование модели процесса переноса железа.

2.2.1. Задание исходных данных и определение констант для железа.

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

Глава 3. Результаты расчета в одномерном приближении.

3.1. Результаты расчета переноса хрома.

3.1.1. Результаты расчета переноса хрома для первого контура БН-600 при номинальном режиме работы, концентрация кислорода в натрии 2 млн"1.

3.1.2. Сравнение с экспериментальными данными.

3.1.3. Результаты расчета переноса хрома в зависимости от концентрации кислорода в натрии: 2, 10, 50, 100, 500 млн"1 и температуры: 550, 650, 750, 850, 950 °С.

3.1.4. Результаты расчета переноса хрома в зависимости от температуры натрия: 550, 650, 750, 850, 950 °С, для концентрации кислорода в натрии: 2, 10, 50, 100, 500 млн"1.

3.2. Результаты расчета переноса железа.

3.2.1. Результаты расчета переноса железа для первого контура БН-600, концентрация кислорода в натрии 2 млн"1.

3.2.2. Результаты расчета переноса железа для температур: 550, 650, 750, 850, 950 °С при различных концентрациях кислорода.

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

Глава 4. Трехмерная модель.

4.1. Физическая постановка задачи.

4.2. Обоснование выбора кода.

4.2.1. Структура пакета программ OpenFOAM.

4.3. Расчет CFD кодом гидравлического сопротивления прямоугольного канала при изменении расстояния между перегородками.

4.4. Численное моделирование тепломассопереноса жидких металлов с помощью пакета программ OpenFOAM.

4.4.1. Особенности моделирования жидкометаллических теплоносителей.

4.4.2. Адаптация OpenFOAM для задач расчета ЯЭУ с ЖМТ.

4.4.3. Моделирование процессов переноса в ПТО реактора БН-600.

4.4.4. Решение задач гидродинамики и сопряженного тепломассообмена.

4.4.5. Результаты расчета переноса примесей Сг и №гО в холодной ловушке.

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

Выводы.

Рост энергетических потребностей в развивающемся мире требует разработки новых технологических процессов, направленных как на использование новых источников энергии, так и на более рациональное и экономичное использование тех энергетических ресурсов, которыми в данное время располагает человек. Исследования и разработки, ведущиеся в этом направлении, весьма обширны. Многие из них связаны с использованием жидких металлов. Об этом убедительно свидетельствуют многочисленные конференции, на которых обсуждаются технологии жидких металлов, используемых при производстве энергии. Анализ трудов этих конференций показывает, что большое внимание в настоящее время уделяется щелочным металлам, а среди них натрию. Области их применения различны, но наибольший практический выход получен при использовании жидких металлов как теплоносителей ядерных энергетических установок (ЯЭУ), и в первую очередь атомных электростанций с реакторами на быстрых нейтронах. О достигнутых здесь успехах убедительно свидетельствует эксплуатация промышленных АЭС БН-600 в России, «Феникс» во Франции и РБЯ в Англии, а также работы, ведущиеся в этих и других странах по созданию более мощных АЭС этого типа.

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

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

Цель диссертационной работы состоит в разработке метода расчета массопереноса продуктов коррозии в натриевых контурах с учетом химического взаимодействия в системе натрий—конструкционный материал-примеси. Обоснование безопасной работы натриевых контуров ЯЭУ на основании полученных данных.

Задачи исследования:

• представить математическое описание физико-химических процессов массопереноса продуктов коррозии с учетом химического взаимодействия в системе натрий - конструкционный материал -примеси;

• разработать расчетную программу для моделирования в одномерном приближении массопереноса конструкционных материалов в натриевых контурах для различных условий эксплуатации ЯЭУ;

• получить физико-химические константы, характеризующие массоперенос в системе натрий - конструкционный материал - примеси;

• разработать трехмерную модель сопряженного тепломассопереноса (переноса примесей);

• выполнить расчеты. Проанализировать и обобщить полученные данные, на основе которых сформулировать научно обоснованные рекомендации для ЯЭУ.

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

• создана методика расчета массопереноса компонентов конструкционных материалов в натриевых контурах с учетом их химического взаимодействия для различных условий эксплуатации ЯЭУ;

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

• получены физико-химические константы, характеризующие массоперенос в системе натрий-конструкционный материал-примеси. Получен ряд простых формул для расчета термодинамических и транспортных теплофизических свойств натрия;

• разработана трехмерная модель сопряженного тепломассопереноса (переноса примесей).

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

• на детальном анализе процессов переноса продуктов коррозии в циркуляционном натриевом контуре и на сопоставлении результатов анализа с известными экспериментальными данными;

• на реализации методов трехмерного моделирования для расчетов процессов массопереноса, результаты которых находятся в хорошем соответствии с литературными данными;

• на системном подходе к проведенным исследованиям, в ходе которых один и тот же результат получался различными методами;

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

Практическая значимость диссертационной работы состоит в том, что:

• разработана методика и вычислительная программа расчета массопереноса продуктов коррозии в натриевых контурах с учетом химического взаимодействия в системе натрий-конструкционный материал-примеси;

• определены экстремальные потоки примесей на стенки каналов для различных условий работы первого контура реактора БН-600, а также для перспективных ЯЭУ;

• получен ряд простых формул для расчета теплофизических свойств натрия;

• проанализированы и обобщены полученные данные по массопереносу в натриевых контурах, сформулированы научно обоснованные рекомендации для ЯЭУ.

Основные положения, выносимые на защиту:

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

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

• результаты расчета массопереноса примесей в холодной ловушке в трехмерной динамической постановке.

Апробация работы: основные положения диссертационной работы докладывались и обсуждались на Всероссийских конференциях и школах: "Реакторы на быстрых нейтронах" (Обнинск, 2009); "Теплофизика-09" (Обнинск, 2009); V Студенческая конференция «Студенчество - будущее атомной энергетики» (Обнинск, 2009); "Реакторы на быстрых нейтронах" (Обнинск, 2010).

На международных конференциях и семинарах: «Безопасность АЭС и подготовка кадров», Обнинск, Россия, 29 сентября-2 октября, 2009; International Conference on Fast Reactors and Related Fuel Cycles: Challenges and Opportunities FR09, Kyoto, Japan, 07-11 December 2009.

Публикации: основное содержание диссертации изложено в трех статьях в реферируемых журналах, в двух препринтах ГНЦ РФ-ФЭИ и двенадцати публикациях в сборниках тезисов докладов и трудах конференций.

Структура и объем диссертации: диссертация состоит из введения, четырех глав, заключения, списка литературы. Материал изложен на 130 страницах, содержит 41 рисунок, 19 таблиц, список литературы из 72 наименований.

Выводы

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

2. В соответствии с предложенной моделью разработана расчетная программа на алгоритмическом языке "Паскаль". Выполнены расчеты и получены новые данные в обоснование безопасности натриевых контуров быстрых реакторов:

По переносу хрома: для номинального режима работы реактора БН-600 с концентрацией кислорода в натрии (2 млн"1) выход хрома в пересчете на изменение толщины стенки канала в экстремальной точке активной зоны (на выходе) составит 0,34 мкм в год, на входе в промежуточный теплообменник 0,06 мкм в год. Увеличение толщины пленки двойного оксида хрома (МаСг02) в каналах на выходе из активной зоны составляет 0,64 мкм в год, на входе в промежуточный теплообменник 0,07 мкм в год, что не приведет к изменению эксплуатационных характеристик первого контура.

В условиях загрязнения поверхности контура двойным оксидом хрома (№Сг02), при достижении уровня температуры теплоносителя 650 °С, при концентрациях кислорода 50 млн"1 и выше скорость роста отложений двойного оксида хрома (ЫаСгСЬ) на выходе из активной зоны 83,5 мкм в год.

По переносу железа: в условиях загрязнения поверхности контура двойным оксидом железа (№4Ре03) для реактора БН-600 при выходе на номинальный режим работы скорость роста отложений двойного оксида железа в виде взвеси на выходе из активной зоны может достигать

9,9 мкм в сутки.

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

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

7 3 концентрация частиц равна 4,5*10" кг/м , при учете химической реакции концентрация взвеси составляет 3,77-10" кг/м (номинальный режим).

4. Получены данные по среднему размеру частиц взвеси двойного оксида хрома NaCr02, при температуре натрия 550 °С и концентрации кислорода в натрии 2 млн"1 средний размер частиц взвеси NaCr02 составляет 0,04 мкм, при увеличении концентрации кислорода средний размер частиц не изменяется. При температуре натрия 650 °С и концентрации кислорода в натрии 2 млн"1 средний размер частиц взвеси NaCr02 составляет 0,09 мкм, при концентрации 500 млн"1 средний размер частиц взвеси равен 0,2 мкм.

5. Сравнение данных по распределению суммарного потока компонентов стали по длине первого контура реактора БН-600 в номинальном режиме работы показывают, что в основном по контуру данные, полученные по расчетам, отличаются от известных экспериментальных данных не более, чем на 5 %. В отдельных зонах с высокой скоростью потока натрия (8 м/с) различие составляет 55 %.

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

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

1. Никитин В.И. Физико-химические явления при воздействии жидких металлов на твердые. М.: Атомиздат, 1967.

2. Жидкометаллические теплоносители. Перев. с англ./ Под ред. А.Е. Шейндлина. М.: Изд. иностр. лит., 1958.

3. Tyzack С. In.: Steel reactor pressure circuits. London, 1961.

4. Horsley G.W. J. Iron and steel inst., 1956. V. 182. - № 1. - P. 43.

5. Kovacina T.A., Miller R.R. Nucl. sei. and engng., 1961. V. 10. - № 2. - P. 163.

6. Williams D.D. et. al. J. Phys. chem. 1959. - V. 63. - № 1. - P. 68.

7. Невзоров Б.А. и др. Коррозия конструкционных материалов в жидких щелочных металлах. М.: Атомиздат, 1977.

8. Kolster В.H. Mechanism of Fe and Cr transport by liquid sodium in nonisothermal loop systems. J. of Nucl. Mat., 1975. V. 55. - P. 155-168.

9. Невзоров Б.А. Коррозия конструкционных материалов в натрии. М.: Атомиздат, 1968.

10. Azad A.M., Sreedharan О.М., Gnanamoothy J.B. J. Nucl. Mater., 1988.-V. 151.-P. 293.

11. Ивановский M.H. и др. О термодинамической устойчивости сложных оксидов в щелочных металлах: Теплофизика высоких температур. 1983. -Т. 21.-№5.-С. 909.

12. Бескоровайный Н.М., Иолтуховский А.Г. Конструкционные материалы и жидкометаллические теплоносители. М.: Энергоатомиздат, 1983.

13. Краев Н.Д., Зотов В.В., Старков О.В. Влияние скорости потока и содержания кислорода в натрии на коррозию сталей //Kernenrgie. 1978. 21. Heft 8.

14. Загорулько Ю.И., Козлов Ф.А., Воробьева Т.А. Дополнительные данные по растворимости индивидуальных веществ в натрии: Препринт ФЭИ-2412. Обнинск, 1997.

15. Краев Н.Д., Умняшкин Е.В., Старков О.В. и др. Влияние физико-химических параметров на скорость массопереноса в жидких щелочных металлах: Препринт ФЭИ-122. Обнинск, 1981.

16. Corrosion by liquid metals. New York: Plenum Press, 1970.

17. Proceeding of the 1st Int. Conf. On Liquid Metal Technology in Energy Production. Champion (PA), USA, May 1976. Conf-760503.

18. Proceeding of the 2nd Int. Conf. On Liquid Metal Technology in Energy Production. Richland (WA), USA, 1980. ERD of Conf-800401.

19. Proceeding of the 4th Int. Conf. On Liquid Metal Engineering and Technology. Avignon, France, Oct. 1988. Ed. SEEN F 75724, Paris.

20. Yong R.S., Lockhart R.W. Sodium component development programmer -General Electric Sodium mass-transfer programmer//Chicago Conf. 650620, 1965.

21. Weeks J.R. //Fall meeting of AIME, Detroit, Oct. 1971.

22. Iizava K. et al. Calculation model and code for corrosion products behaviour in primary circuits of LMFBRs //Karlsrue, Germany, 1987. /Vienna: IAEA, IWGFR/64. KFK Report 4279. P. 191-225.

23. Матусевич JT.H. Кристаллизация из растворов в химической промышленности. М.: Химия, 1968. - С. 60, 84, 96.

24. Batchelor G.K. Mass transfer from small particles suspended in turbulent fluid // J. Fluid Mech. 1980. - V.98, №3. - P. 609-623.

25. Уббелоде А. Плавление и кристаллическая структура. Пер. с англ. М.: Мир, 1969.-420 с.

26. Мелихов И.В., Меркулова М.С. Сокристаллизация. М.: Химия, 1975. -С. 16.

27. Каргельский И.В., Виноградова И.Э. Коэффициенты трения. М.: Машгиз, 1962. - 220 с.

28. Турчин Н.М., Дробышев А.В. Экспериментальные жидкометаллические стенды. -М.: Атомиздат, 1978. С. 27.

29. Glaeser W.A. Wear and friction characteristics of structural materials in liquid sodium//Reactor Technol. 1972. - V.15, №1. - P. 1-10.

30. Алексеев B.B., Козлов Ф.А., Загорулько Ю.И. Расчет критической скорости транспортировки взвесей в потоке теплоносителя // Атомная энергия. 1984. - Т.57, №7. - С. 14-18.

31. Алексеев В.В., Козлов Ф.А., Загорулько Ю.И., Лихарев В.А., Краев Н.Д. Исследование динамики поведения взвесей в натриевых контурах быстрых реакторов: препринт ФЭИ-2576. Обнинск, 1996. — 16 с.

32. Кириллов П.Л., Юрьев Ю.С. и др. Справочник по теплогидравлическим расчетам, 2-е изд. М.: Энергоатомиздат, 1990.

33. Романков П.Г., Курочкина М.И. Гидромеханические процессы химической технологии. М.: Химия, 1974. - 288 с.

34. Rouhiainen P.O., Stachiewicz J.W. On the deposition of small particles from turbulent streams // J-. of Heat Transfer, Trans. ASME, Series C. 1970. - V.92, №1.- P. 169-177.

35. Морозова И.К. и др. Вынос и отложения продуктов коррозии реакторных материалов. М.: Атомиздат, 1975. - 280 с.

36. Юфин А.П. Движение неоднородной жидкости по горизонтальным стальным незаиленным трубам // Изв. АН СССР, ОТН. 1949. - №8. - С. 1146-1159.

37. Пирумов А.И. Обеспыливание воздуха. М.: Стройиздат, 1974. - С. 33.

38. Шифрин К.С., Гордон И.З., Файнштейн М.Г. Экспериментальное определение константы турбулентной диффузии внутри замкнутого сосудапри малых градиентах температуры. Изв. АН СССР. Сер. геогр. и геофиз., 1949.-№3. -С. 300-310.

39. Руцков А.П. Краткий курс коллоидной химии. М.: Госхимиздат, 1958.

40. Левин В.Г. Физико-химическая гидродинамика. М.: Физматгиз, 1959. -699 с.

41. Nicholas M.G., Cavell I.W. Materials Development Division, A.E.R.E. // Proc. Second Int. Conf. on Liquid Metal Technology in Energy Production, 20-24 April 1980, Richland/CONF-800401, p. 3-35 3-41.

42. Алексеев В.В. Массоперенос трития и продуктов коррозии конструкционных материалов в контурах с натриевым теплоносителем: диссертация на соискание ученой степени доктора технических наук. — ГНЦ-РФ ФЭИ, Обнинск, 2002.

43. Алексеев В.В., Кондратьев А.С. Моделирование массопереноса продуктов коррозии в контурах ЯЭУ с натриевым теплоносителем // Известия вузов. Ядерная энергетика. 2010. - № 4. - С. 162-171.

44. Hooper A.J. The Sodium Compatibility of 9Cr-lMo Ferritic Steel // Proc. Second Int. Conf. on Liquid Metal Technology in Energy Production, 20-24 April 1980, Richland/CONF-800401, p. 7-35 7-44.

45. Shaiv В .J., Wu P.C.S., Chiotti P. Thermodynamic Properties of the Double Oxides of Cr, Ni and Fe // J. Nucl. Mater. 1977. - V. 67. - P. 13-23.

46. Краев Н.Д., Коррозия и массоперенос конструкционных материалов в натриевом и натрий-калиевом теплоносителях // Известия вузов. Ядерная энергетика. 1999. - № 3. - С. 40-48.

47. Henric Rusche. Computational Fluid Dynamics of Dispersed Two-Phase Flows at High Phase Fractions. PhD Thesis, University of London, Imperial College of Science, Technology and Medicine, Department of Medical Engineering, 2002.

48. Кумаев В.Я., Леончук М.П. Численное моделирование трехмерного течения теплоносителя в пучках твэлов: Атомная энергия, 1986. — Т. 60. — Вып. 6.-С. 382-386.

49. OpenFOAM User Guide Version 1.6. 2009.

50. OpenFOAM Programmer's Guide Version 1.6. 2009.

51. The Open Source Integration Platform for Numerical Simulation URL: http://www.salome-platform.org/ (дата обращения 20 сентября 2010).

52. Gmsh: 3-D finite element mesh generator URL: http://geuz.org/gmsh/ (дата обращения 16 июля 2010).

53. Ефанов А.Д., Левченко Ю.Д., Федотовский B.C., Щукин Н.М. Гидравлические потери в контуре с взаимным влиянием местных сопротивлений // Теплоэнергетика. 1997. - Т. 3. - С. 8-13.

54. Казанцев А.А., Кондратьев А.С., Левченко Ю.Д. Расчет CFD кодом гидравлического сопротивления прямоугольного канала при изменении расстояния между перегородками // Известия вузов. Ддерная энергетика. — 2010.-№ 4.-С. 193-198.

55. Eckert, S., Gerbeth, G., "Local velocity measurements in lead-bismuth and sodium flows using the Ultrasound Doppler Velocimetry". NURETH-10, Seoul (Korea), Oct. 5-9, 2003.

56. Казанцев A.A., Кондратьев A.C. Простые формулы для термодинамических свойств натрия. // Известия вузов. Ядерная энергетика. -2008. № 4. - С. 86-94.

57. Open Source CFD URL http://www.os-cfd.ru/ (дата обращения 16 июля 2010).

58. Персональный сайт Варсеева Евгения URL: http://varseev.heliohost.org (дата обращения 16 июля 2010).

59. Кириллов П.Л., Терентьева М.И., Денискина Н.Б. Теплофизические свойства материалов ядерной техники: Учебное справочное пособие длястудентов / под общ. ред. проф. П.Л. Кириллова; 2-е изд., перераб. и доп. -М.: ИздАтб, 2007. 200 с.

60. Кирилов П.Л., Галин Н.М., Грошев А.И., Слободчук В.И. О влиянии свойств стенки на теплоотдачу при турбулентном течении жидких металлов в трубах // Теплоэнергетика. 1984. - № 3. - С. 63-64.

61. Грошев А.И., Слободчук В.И. Численный расчет нестационарного турбулентного теплообмена в круглой трубе с учетом растечек теплоты в стенке // Теплоэнергетика. 1986. - № 9. - С. 56-59.

62. Александров A.A. Григорьев Б.А. Стандарт воды IAPWS-IF97. http://www.iapws.org/ (дата обращения: 25.06.2010).

63. Таблицы теплофизических свойств воды и водяного пара: Справочник. Рек. Гос. Службой стандартных справочных данных. ГСССД Р-776-98 М.: Издательство МЭИ, 1999. - 168 с.

64. Software for the Industrial-Formulation IAPWS-IF97 for Water und Steam http://www.ruhr-uni-bochum.de/thermo/Software/Seiten/IAPWS-IF97-eng (дата обращения: 25.06.2010).

65. Large eddy simulation (LES) http://www.cfd-online.com/Wiki/LargeEddySimulation (дата обращения: 10.10.2010).

66. Боришанский В.М., Кутателадзе С.С., Новиков И.И., Федынский О.С. Жидкометаллические теплоносители./Изд. 3-е. — М.: Атомиздат, 1976. 328 с.

67. Субботин В.И., Козлов Ф.А. В сб. "Alkali metal coolants", IAEA, Vienna, 1967 (SM-85/48).

68. Жидкометаллические теплоносители ЯЭУ. Очистка от примесей и их контроль / Под ред. д-ра техн. наук Козлова Ф.А. / Козлов Ф.А., Волчков Л.Г., Кузнецов Э.К., Матюхин В.В. -М.: Энергоатомиздат, 1983.

69. Баклушин Р.П., Козлов Ф.А. О схемах очистки натрия первого и второго контуров быстрых реакторов. — Атомная энергия, 1978, т. 44, вып. 3, с. 224228.