Массоперенос продуктов коррозии с учетом химического взаимодействия в системе натрий - конструкционный материал - примеси тема диссертации и автореферата по ВАК РФ 05.14.03, кандидат технических наук Кондратьев, Александр Сергеевич

• Кондратьев, Александр Сергеевич
• кандидат технических наук
• 2011

Оглавление диссертации кандидат технических наук Кондратьев, Александр Сергеевич

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

Введение.

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

Список литературы диссертационного исследования кандидат технических наук Кондратьев, Александр Сергеевич, 2011 год

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.