УДК 621.039.743
Ю.Е. Карякин, А.А. Плетнев, Е.Д. Федорович
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ТЕПЛОМАССООБМЕНА
ПРИ АВАРИЙНОЙ СИТУАЦИИ В ХРАНИЛИЩЕ ОТРАБОТАВШЕГО ЯДЕРНОГО ТОПЛИВА
Yu.E. Karyakin, A.A. Pletnev, E.D. Fedorovich
NUMERICAL SIMULATION OF THE HEATMASSTRANSFER PROCESSES IN EMERGENCY SITUATION AT A STORAGE FOR SPENT NUCLEAR FUEL
В статье представлена математическая модель и результаты расчетов по определению локальных значений температур в хранилище отработавшего ядерного топлива в различные моменты времени как в штатной, так и в аварийной ситуации, связанной с прекращением циркуляции воды через бассейн. Предложен способ управления аварийной ситуацией в бассейне с использованием гидроаккумулятора.
АТОМНАЯ ЭЛЕКТРОСТАНЦИЯ. ОТРАБОТАВШЕЕ ЯДЕРНОЕ ТОПЛИВО. ХРАНИЛИЩЕ «МОКРОГО» ТИПА. АВАРИЙНАЯ СИТУАЦИЯ. ГИДРОАККУМУЛЯТОР. ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ.
In the paper are presented the mathematical model and the results of calculations of local temperatures in «wet» type storage for spent nuclear fuel. There are considered normal and accident non-stationary situations, connected with the stopping of water circulation through pool. It is suggested a mode of control with use of hydroaccumulators.
NUCLEAR POWER PLANT. SPENT NUCLEAR FUEL. «WET» TYPE STORAGE. EMERGENCY SITUATION. HYDRAULIC ACCUMULATOR. NUMERICAL SIMULATION.
В настоящее время большая часть отработавшего ядерного топлива (ОЯТ) атомных электростанций (АЭС) находится в условиях так называемого «мокрого» хранения в хранилищах отработавшего ядерного топлива (ХОЯТ). Это бассейны, через которые прокачивается вода, нагреваемая теплотой остаточного энерговыделения ОЯТ и охлаждаемая затем в теплообменниках внешней системы охлаждения. На Ленинградской АЭС охлаждающей средой теплообменников служит вода Финского залива.
ХОЯТ «мокрого» типа располагаются как на территории атомной электростанции, так и на территориях предприятий ядерного топливного цикла, предназначенных для последующей химической переработки ОЯТ, либо для его «окончательного» захоронения. Одновременно в одном ХОЯТ может храниться до нескольких десятков тысяч отработавших тепловыделяющих сборок (ОТВС). Мощность
тепловыделения одной ОТВС относительно невелика, например для ХОЯТ Ленинградской АЭС — от 30 до 80 Вт в зависимости от промежутка времени, прошедшего после выгрузки из реактора. Однако при наличии большого числа ОТВС общая мощность тепловыделения в хранилище оказывается достаточной, чтобы при прекращении принудительной циркуляции воды через отсеки бассейна температура воды спустя несколько суток достигла высоких значений, при которых интенсивность испарения становится недопустимо большой.
Именно такая ситуация возникла на АЭС «Фукусима-1» (Япония) в 2011 году [1]. Как известно, ее последствия были крайне тяжелыми: полная потеря (выпаривание) воды из бассейна, перегрев циркониевых оболочек тепловыделяющих элементов, пароциркониевая реакция с выделением водорода, образованием гремучей смеси и взрывом. Известны также случаи утечек
и прекращения циркуляции воды через бассейны хранения ОЯТ атомного предприятия в Сел-лафилде (Англия) [2].
Один из возможных способов охлаждения отработавшего топлива в условиях прекращения циркуляции воды через бассейн — использование предварительно запасенной в специальных баках-гидроаккумуляторах относительно холодной воды, подаваемой автоматически в бассейн самотеком по сигналу противоаварийной системы, действующей даже в условиях полного длительного обесточивания объектов АЭС.
Для условий Российской Федерации возможные причины потери электроснабжения объектов АЭС и потерь воды из бассейнов ХОЯТ отличаются от зарубежных. Например, для российских АЭС возникновение цунами маловероятно, но причиной аварии могут быть иные природные или техногенные катастрофы: землетрясения, смерчи, наводнения, падение самолетов с горящим топливом или атаки террористов, вероятность которых, к сожалению, нельзя исключить. Во всех этих случаях может возникнуть длительное прекращение циркуляции воды через бассейны ХОЯТ и ее усиленное испарение.
Данное обстоятельство делает актуальной разработку сценариев возникновения и развития гипотетических аварий, которые могут быть в результате длительного обесточивания объектов АЭС и прекращения циркуляции воды через бассейн хранения отработавшего ядерного топлива. Наиболее полная информация в этой области может быть получена с помощью методов математического моделирования. Этому вопросу и посвящена данная статья.
Расчетная модель ХОЯТ
В центре «Техническая диагностика и надежность АЭС и ТЭС» Санкт-Петербургского государственного политехнического университета разработана математическая модель, предназначенная для определения локальных значений температур элементов ХОЯТ в различные моменты времени как в штатной, так и в аварийной ситуации, связанной с прекращением циркуляции воды через бассейн.
Расчетная модель ХОЯТ представляет собой приведенную к одной ОТВС элементарную осе-симметричную ячейку хранилища ОЯТ, содержащую все участвующие в теплообмене элементы
натурного хранилища отработавшего ядерного топлива: тепловыделяющие элементы, пеналы, воду пеналов и бассейна, торцевое и боковое ограждения бассейна и т. д.
В модели все вертикальные размеры, размеры пенала и его содержимого сохранены натурными (рис. 1), а межпенальное пространство смоделировано более простым для рассмотрения, но подобным в отношении конвективной теплопередачи кольцевым каналом. При этом радиальные размеры модели выбраны так, чтобы в расчете на одну ОТВС в модели и в натурном ХОЯТ были равны следующие площади: а) поперечного сечения межпенального пространства; б) поперечного сечения боковых
I
Рис. 1. Вертикальное сечение двумерной осесимметричной модели элементарной ячейки ХОЯТ:
1 — центральный (несущий) стержень; 2, 3 — твэлы внутреннего и внешнего ряда; 4 — пространство между твэлами в пенале; 5 — стенка пенала; 6 — межпенальное пространство; 7 — внешнее бетонное ограждение; 8 — днище бассейна
бетонных ограждений; в) поверхности испарения внутри пенала. Таким образом, в расчете на одну ОТВС в модели и в натурном хранилище соблюдено равенство объемов, масс и тепло-емкостей рассматриваемых элементов ХОЯТ.
Условия теплопередачи через внешние ограждения бассейна определены с использованием критерия краевого подобия Био, представляющего собой отношение внутреннего и внешнего термических сопротивлений, В1 = а8/К , где а — коэффициент конвективной теплоотдачи, 5 — толщина стенки, X — коэффициент ее теплопроводности.
Поскольку в реальном хранилище лишь небольшая часть пеналов непосредственно примыкает к боковым стенам бассейна, внешний периметр кольцевого бетонного ограждения модели в расчете на одну ОТВС оказывается значительно больше, чем в натуре, а его толщина соответственно меньше. По этим причинам «эффективный» (расчетный) коэффициент теплопроводности эквивалентного бетонного ограждения оказывается на несколько порядков меньше, чем у железобетона. Необходимое условие Вмодел = В/натур выполнено путем соответствующего выбора коэффициента конвективной теплоотдачи а.
Аналогично [3] предполагалось, что источник тепловыделения обусловленный поглощением гамма-квантов и электронов, испускаемых продуктами деления, сосредоточен в топливных элементах, тепловая мощность которых известна и одинакова для всех ОТВС. Система дифференциальных уравнений тепломассообмена для элементарной ячейки ХОЯТ имеет вид [4]
эт 1 э (. эт ттт
ср-=--1 гк--сригт
Э t г Э г I Э г
(
Э у
к—-срУТ Э у
+Ог;
Э с 1 Э
.Э с
- = -—I гБ---игС
Э t г Э г I Э г
(
Э у
Б—-УС Э у
(
Эи ттЭи ?и
-+и-+У-
Э t Э г Э у
Эр . 1 (Эи и
= -_д +2Ц I д---
Э г г I Э г г
2 ~ , Эи
+2 н Ц-^-
Э г I Э г
Э у
(
Ц
- V
Эи + ЭУ Э у Э г
+р/г ;
(1)
ЭУ ТТЭУ „ЭУ
-+и-+У-
Э t Э г Э у
ЭР
= -—+ 2Ц-Э у г
(ЭУ Эи
+ -
Э г Э у
+ 2— Э у
(
ЭУ
Э у
Э г
(
Ц
-
Эи + ЭУ Э у Э г
+р/у ;
+^ = 0 ; Г=?|!(Т„-Т).
г Э г Э у
Здесь: t — время; г — радиальная координата; у — вертикальная координата; Т — температура; С — концентрация водяного пара; и, У — компоненты вектора скорости соответственно вдоль осей г, у; р — давление; с — удельная теплоемкость; г — плотность; Б — коэффициент диффузии; т — коэффициент динамической вязкости; О — удельная мощность источника тепловыделения; Ь — коэффициент теплового объемного расширения; g — вектор ускорения силы тяжести; Та—«опорная» температура для свободной конвекции; Г — вектор силы Архимеда.
Система дифференциальных уравнений тепломассообмена (1) для элементарной ячейки ХОЯТ представлена в ламинарной постановке [4]. В случае турбулентного движения авторами использовалась низкорейнольдсовая ^-е-модель турбулентности [5].
Численное решение системы дифференциальных уравнений (1) с соответствующими граничными и начальными условиями выполнено неявным конечно-разностным методом.
Граничные условия для решения задачи заданы на основе экспериментальных данных эксплуатации ХОЯТ Ленинградской АЭС; начальные условия определены из стационарного решения исходной системы уравнений в штатном режиме эксплуатации ХОЯТ.
В расчетах физические свойства паро-воз-душной смеси зависят от температуры и относительной влажности, а физические свойства воды и твердых сред были приняты постоянными.
Высота уровня воды в межпенальном пространстве в момент возникновения аварии принята 10,7 м, в пенале — 10,6 м. Начальная сред-необъемная температура воды в пенале — 40 °С, что соответствует штатным условиям эксплуатации бассейнов ХОЯТ. Все расчеты выполнены
при заполнении отсеков ОЯТ на 80 % от максимальной проектной вместимости (7100 ОТВС на один отсек) для тепловой мощности одной ОТВС, равной 100 Вт, что примерно соответствует ее остаточному тепловыделению после 10-летней выдержки.
Результаты расчетов
Согласно результатам численного моделирования на начальном этапе развития аварии через несколько часов после прекращения циркуляции теплоносителя тепловое состояние бассейна ХОЯТ имело следующие особенности. Вследствие перемешивания за счет действия свободной конвекции при весьма малой плотности теплового потока на поверхности твэлов вода в пенале и межпенальном пространстве имеет практически одинаковую температуру во всем объеме.
Темп разогрева воды составляет 0,43 °С/час, темп испарения — 1,6 г/час в пенале и 1,9 г/час в межпенальном пространстве. С учетом разницы площадей зеркала испарения (7,4110-3 м2 в пенале и 12,6510-3 м2 в межпенальном пространстве) скорость понижения уровня воды в пенале
а)
несколько выше скорости снижения уровня воды, заполняющей бассейн. Так как вентиляция отсутствует, пар удаляется из-под настила только за счет диффузии (скорость вблизи настила равна нулю, концентрация постоянна).
Для моделирования более поздних стадий развития аварии высота уровня воды в пенале была задана равной 2,3 м, средняя температура воды — 92 °С. Межпенальное пространство не рассматривалось, так как на основании результатов предыдущих расчетов была зафиксирована плотность теплового потока на внешней границе пенала, равная 20 Вт/м2 в области твэла и 5 Вт/м2 между настилом и твэлом. Температура настила бассейна — 50 °С.
Максимальная температура твэла (рис. 2), зафиксированная в области между уровнем воды и верхним торцом твэла, составляет 120 °С и увеличивается с темпом приблизительно 3,5 °С/час. Верхний торец твэла имеет более низкую температуру (96 °С) благодаря конвективному теплообмену с влажным воздухом, циркулирующим в пространстве между твэлом и верхним настилом. Температура воды, имеющейся в пенале, не увеличивается. Это свидетельствует
б)
Рис. 2. Расчетные поля температур (а), осредненной скорости и концентрации (б) в пенале при высоте уровня воды в пенале 2,4 м
о том, что вся теплота, поступающая к воде от твэла, расходуется на ее испарение и частично на теплоотдачу от внешней поверхности пенала. Темп испарения составляет 16,2 г/час, что в 10 раз больше, чем в начальном тепловом состоянии при температуре воды в бассейне 40 °С. Можно утверждать, что при дальнейшем снижении уровня воды в пенале темп испарения будет уменьшаться вследствие уменьшения площади контакта воды с теплоотдающей поверхностью твэла, а также из-за увеличения гидравлического сопротивления выходу пара из пенала.
Если расчетная концентрация пара вблизи поверхности испарения или около ограничивающих поверхностей превышала параметры насыщения, то в этом случае моделировался процесс конденсации: концентрация пара в соответствующих узлах уменьшалась до значений насыщения (при температуре в узле) и вычислялась масса конденсата, которая достигала 20 % от массы испарившейся жидкости.
Для оценки влияния температуры воды в бассейне на тепловое состояние ХОЯТ была выполнена серия из трех расчетов, в которых средняя температура воды в ХОЯТ была при-
близительно равна 40, 60 и 80 °С, а температура верхнего настила принята равной соответственно 30, 40 и 50 °С.
Исходные данные, одинаковые для всех расчетов серии: высота уровня воды в пенале — 9,1 м; высота уровня воды в бассейне — 9,2 м; относительная влажность вблизи настила — 50 %.
Вентиляция пространства под настилом отсутствовала, люки верхнего настила предполагались закрытыми.
Результаты серии расчетов обобщены в таблице.
Поскольку значения тепловых потоков и темпа испарения, указанные в таблице, относятся к элементарной ячейке ХОЯТ, содержащей один пенал, параметры для всего отсека ХОЯТ могут быть получены путем умножения на число пеналов в отсеке.
Расчеты показали, что темп роста температуры воды в бассейне на начальном этапе развития аварии составил 0,43 °С/час; время полного испарения воды из пеналов и бассейна — не менее 48 суток. До момента полного испарения воды из ХОЯТ температура твэлов не превысит 200 °С.
Таблица 1
Название параметра Значения параметров воды в бассейне и пенале при разных температурах
Т ~ * ср 40 °С Т ~ * ср 60 °С Т ~ * ср 80 °С
Пенал Бассейн Пенал Бассейн Пенал Бассейн
Максимальная температура, °С 42,7 41,2 60,2 58,7 79,5 77,8
Средняя по объему температура воды, °С 41,1 40,6 58,5 58,0 77,6 77,1
Тепловой поток через ограждения, Вт 0,5 3,0 1,9 8,0 3,6 11,6
Мощность на нагрев воды, Вт 23,8 64,1 19,7 53,0 9,2 22,9
Мощность на нагрев ОТВС, Вт 5,8 - 5,4 - 2,9 -
Мощность на испарение, Вт 1,2 1,4 5,2 6,8 14,6 35,2
Тепловой поток за счет испарения, кВт/м2 0,16 0,11 0,7 0,54 2,0 2,8
Темп испарения, г/час 1,75 2,05 7,7 10,1 21,6 52,1
Темп снижения уровня воды, мм/час 0,24 0,16 1,04 0,80 2,92 4,13
Темп роста ср. температуры воды, °С/час 0,44 0,44 0,34 0,33 0,17 0,15
Темп роста тах температуры ТВЭЛ, °С/час 0,47 - 0,39 - 0,19 -
Сравнительные результаты расчетов аварийных режимов ХОЯТ, выполненных для трех значений температуры воды в бассейне
Использование гидроаккумуляторов
Для предотвращения тяжелых последствий аварии, связанной с прекращением циркуляции воды через бассейн ХОЯТ, авторами настоящей работы предложено использовать автономную энергонезависимую систему подачи воды в бассейн. Основной элемент системы — резервуар (бак-гидроаккумулятор), заполненный водой и расположенный выше уровня воды в бассейне. Автоматический клапан открывает подачу воды, если ее уровень в бассейне опустится ниже установленного. Клапан остается открытым до тех пор, пока уровень воды не восстановится, т. е. не поднимется до заданной верхней отметки. Долив воды в бассейн происходит порциями. Система функционирует, пока имеется вода в гидроаккумуляторе (считается, что целостность корпуса бассейна не нарушена, т.е. аварийной течи из него нет).
Выполняя тепловой расчет бассейна ХОЯТ в аварийном режиме, мы имели цель определить: установившиеся (квазистационарные) тепловые параметры при работающей системе автономной подачи воды; время достижения такого квазистационарного теплового состояния;
среднечасовой расход воды из гидроаккумулятора в установившемся тепловом режиме.
Как показали результаты численного моделирования, из-за относительно малой удельной мощности остаточного тепловыделения ОЯТ (менее 4 кВт/м3 после 10 лет хранения) градиенты температур в бассейнах выдержки малы и не превышают 2 °С как в штатном режиме, так и в режиме аварийного нагрева (см. таблицу). Поэтому для расчета тепловых параметров ХОЯТ при работающей системе автономной подачи воды был выбран метод интегрального теплового баланса, дифференциальное уравнение которого в рассматриваемом случае имеет вид
dt _ £>оят бконв °исп1 °исп2 Обет
dt
C
, (2)
басс
мощность, расходуемая на нагрев бетонных ограждений бассейна и теплопередачу через указанные ограждения к окружающему воздуху.
Полная теплоемкость бассейна за вычетом внешних ограждений складывается из теплоемкости ОТВС, пеналов и теплоемкости воды, т. е.
Сбасс = СОТВС + твсв .
Изменение массы воды в пеналах (нижний индекс «1») и бассейне (нижний индекс «2») выразим так:
m*1 - _jiSMcni+Gi;
dт
dm,
в2
dт
- -j2SMcn2 +G2-
(3а)
(3б)
где Т — средняя по объему ХОЯТ температура воды; I — время; Сбасс — полная теплоемкость бассейна с ОЯТ; 0ОЯТ — полная мощность остаточного тепловыделения ОЯТ; (2конв — конвективный поток тепла с поверхности испарения; 0исп1 — мощность, расходуемая на испарение воды из пеналов; бисп2 — мощность, расходуемая на испарение воды из бассейна; 0бет — тепловая
Здесь: тв — масса воды в пеналах (mBi) или бассейне (тв2), G — темп поступления воды (массовый расход воды из гидроаккумулятора); j — массовый поток пара с единицы поверхности испарения; £исп — площадь поверхности испарения.
Удельный поток испарения равен j - Рм (Рп " Рппв), где Рм — коэффициент массо-отдачи; рп — плотность водяного пара вблизи поверхности испарения (равна плотности насыщения при температуре воды); рп пв — плотность водяного пара в воздухе помещения.
Коэффициент массоотдачи Рм обычно находят опытным путем, но если соответствующие экспериментальные данные отсутствуют, то его приближенное значение можно найти из аналогии тепло- и массообмена [6]:
в -а.
Рм cp
Коэффициент теплоотдачи а был рассчитан по критериальным соотношениям для свободной конвекции на плоской горизонтальной поверхности, обращенной вверх [7]:
Nu = 0,96Ra1/6 при 1 < Ra < 200;
Nu = 0,54Ra1/4 при 200 < Ra < 8-106.
В качестве характерного размера при вычислении критериев Нуссельта и Рэлея взято отношение площади отсека бассейна к его периметру.
Обе величины — jk и Gk — изменяются во времени. При этом зависимость G(t) имеет вид 5-функции конечной ширины, что соответствует двум возможным положениям клапана подачи воды — «открыт» или «закрыт».
а) б)
Рис. 3. Зависимость массы воды от времени в отсеке бассейна ХОЯТ, заполненного ОЯТ 10-летней (а) и 30-летней (б) выдержки: 1 — воды, испаренной из отсека (бассейн и пеналы); 2 — воды, испаренной из пеналов; 3 — воды, испаренной из бассейна; 4 — воды, добавленной в отсек из бака-аккумулятора
Суммарная масса воды mисп, испарившейся из бассейна за время t, есть результат интегрирования
t
тисп
0
Уравнения (2) и (3) были решены численно методом конечных разностей.
Расчеты выполнены для одного отсека бассейна ХОЯТ Ленинградской АЭС. Принятое
Рис. 4. Зависимость температуры воды от времени в отсеке бассейна ХОЯТ, заполненного ОЯТ: 1-10-летней выдержки, 2-30-летней выдержки
(о — момент открытия клапана подачи воды в бассейн; □ — момент закрытия клапана подачи воды в бассейн)
в расчете суммарное тепловыделение для ОЯТ 10-летней выдержки равно 610 кВт/отсек, для ОЯТ 30-летней выдержки — 370 кВт/отсек.
Расход (темп долива) воды в бассейн и пеналы из гидроаккумулятора, одинаковый для бассейна и пеналов, был принят 2 л/с. Температура доливаемой воды — 15 °С. Аварийный клапан открывает подачу воды при снижении уровня воды в бассейне и пеналах на 0,4 м.
Влажность воздуха в здании ХОЯТ не вычислялась, она была принята постоянной и равной начальному значению — 60 %.
Результаты расчетов для ОЯТ 10-летней и 30-летней выдержки представлены на рис. 3, 4.
Для принятых исходных данных во время каждого срабатывания клапана из гидроаккумулятора в отсек бассейна ХОЯТ добавляется около 60 тонн воды (рис. 3, а, кривая 4). Масса единовременно добавляемой воды составляет менее 5 % от массы воды в отсеке на момент открытия клапана, в связи с чем колебания температуры воды в бассейне в установившемся тепловом режиме не превысили 4 °С.
Время выхода аварийного отсека ХОЯТ на квазистационарный тепловой режим и температура воды в таком режиме зависят от количества пеналов и суммарной мощности остаточных тепловыделений ОЯТ. В исследованном диапазоне ООЯТ = 400-600 кВт/отсек время достижения установившегося теплового режима составило от 15 до 25 суток, а максимальная
температура воды в ХОЯТ — соответственно 93 °С и 79 °С, (см. рис. 4).
Автономность системы аварийной подачи воды в ХОЯТ должна быть не менее 10 суток. Исходя из расчетного темпа испарения воды из отсека с ОЯТ 10-ти летней выдержки следует рекомендовать иметь запас воды в гидроаккумуляторе около 100 тонн на отсек.
Методом численного моделирования на основе решений уравнений тепломассопереноса и уравнения теплового баланса выполнен расчет теплового состояния бассейна хранения отработавшего ядерного топлива в аварийной
ситуации, вызванной прекращением принудительной циркуляции воды через бассейн (на примере Ленинградской АЭС).
Темп роста среднеобъемной температуры воды в бассейне на начальном этапе развития аварии составил 0,43 °С/час. Расчетное время полного испарения воды из пеналов и бассейна — не менее 48 суток.
Показано, что система автономной подачи воды в бассейн из гидроаккумулятора в состоянии стабилизировать температурный режим ХОЯТ в аварийной ситуации и поддерживать его на безопасном уровне до тех пор, пока имеется запас воды в гидроаккумуляторе.
СПИСОК ЛИТЕРАТУРЫ
1. Катастрофа на АЭС «Фукусима» [Электронный ресурс].— Режим доступа: http://nuclphys.sinp.msu.ru/ ecology/fukushima/index.html.
2. Bolter, H. Inside Sellafield [Текст] / H. Bolter.— London: Quartet Books, 1996.— 298 p.
3. Фромзель, В.Н. Методика определения эффективной теплопроводности сборки твэлов и расчет температурного поля в сборках, размещенных в вертикальных контейнерах [Текст] / В.Н. Фромзель, Л.В. Фромзель, Н.В. Вдовец // Процессы тепломассообмена и гидродинамики в системах безопасности АЭС с ВВЭР-640: Сб. трудов.— СПб: Изд-во АООТ «НПО ЦКТИ», 1997.— С. 139-150.
4. Astafieva, V. Modeling of heatmasstransfer in «wet» and «dry» storages of spent nuclear fuel— normal and accident situations [Текст] / V. Astafieva, E. Fedorovich,
Yu. Karyakin, A. Pletnev, A. Trifonov // 2011 Baltic Heat Transfer Conference— 6th BHTC (August 24-26, 2011, Tampere, Finland).— Extended Abstracts.— Tampere: Univ. of Technology, 2011.— P. 81-82.
5. Launder, B.E. On the computation of convective heat transfer in complex turbulent flows [Текст] / B.E. Launder // Transaction of the ASME: Journal of Heat Transfer.— 1988. Vol. 110, №4 (B).— P. 11121128.
6. Гухман, А.А. Применение теории подобия к исследованию процессов тепло-массообмена. (Процессы переноса в движущейся среде) [Текст] / А.А. Гухман.— М.: Высшая школа, 1974.— 328 С.
7. Мартыненко, О.Г. Свободно-конвективный теплообмен [Текст]: Справочник / О.Г Мартыненко, Ю.А. Со-ковишин.— Минск: Наука и техника, 1982.— 400 с.
REFERENCES
1. Katastrofa na AES «Fukusima» [Elektronnyi resurs].— Rezhim dostupa: http://nuclphys.sinp.msu.ru/ ecology/fukushima/index.html.
2. Bolter H. Inside Sellafield [Text] / H. Bolter.— London: Quartet Books, 1996.— 298 p.
3. Fromzel' V.N., Fromzel' L.V., Vdovets N.V. Meto-dika opredeleniia effektivnoi teploprovodnosti sborki tvelov i raschet temperaturnogo polia v sborkakh, razmesh-chennykh v vertikal'nykh konteinerakh [Tekst] // Protsessy teplomassoobmena i gidrodinamiki v sistemakh bezopas-nosti AES s VVER-640: Sb. trudov.— SPb: AOOT «NPO TsKTI», 1997.— S. 139-150. (rus.)
4. Astafieva V., Fedorovich E., Karyakin Yu., Pletnev A., Trifonov A. Modeling of heatmasstransfer in «wet» and «dry» storages of spent nuclear fuel— normal and accident
situations [Text] // 2011 Baltic Heat Transfer Conference- 6th BHTC (August 24-26, 2011, Tampere, Finland).— Extended Abstracts.— Tampere Univ. ofTechnol-ogy, 2011.— P. 81-82.
5. Launder, B.E. On the computation of convective heat transfer in complex turbulent flows [Текст] / B.E. Launder // Transaction of the ASME: Journal of Heat Transfer.— 1988.— Vol. 110, №4 (B).— P. 1112-1128.
6. Gukhman, A.A. Primenenie teorii podobiia k issle-dovaniiu protsessov teplo-massoobmena. (Protsessy perenosa v dvizhushcheisia srede) [Tekst] / A.A. Gukhman.— M.: Vysshaia shkola, 1974.— 328 s. (rus.)
7. Martynenko O.G., Sokovishin Yu.A. Svobodno-konvektivnyi teploobmen [Tekst]: Spravochnik.— Minsk: Nauka i tekhnika, 1982.— 400 s. (rus.)
СВЕДЕНИЯ ОБ АВТОРАХ /AUTHORS
КАРЯКИН Юрий Евгеньевич — доктор физико-математических наук профессор директор Центра технической диагностики и надежности АЭС и ТЭС Санкт-Петербургского государственного политехнического университета; 195251, ул. Политехническая, 29, Санкт-Петербург, Россия; e-mail: yu_kar@mail.ru
KARYAKIN Yuriy E. — St.-Petersburg State Polytechnic University; 195251, Politekhnicheskaya Str. 29, St.-Petersburg, Russia; e-mail: yu_kar@mail.ru
ПЛЕТНЕВ Александр Александрович — кандидат технических наук старший преподаватель Санкт-Петербургского государственного политехнического университета; 195251, ул. Политехническая, 29, Санкт-Петербург, Россия; e-mail: aapletnev@yandex.ru
PLETNEV Alexandr A. — St.-Petersburg State Polytechnic University; 195251, Politekhnicheskaya Str. 29, St.-Petersburg, Russia; e-mail: aapletnev@yandex.ru
ФЕДОРОВИЧ Евгений Данилович — доктор технических наук профессор Санкт-Петербургского государственного политехнического университета; 195251, ул. Политехническая, 29, Санкт-Петербург, Россия; e-mail: efed1@yandex.ru
FEDOROVICH Evgeny D. — St.-Petersburg State Polytechnic University; 195251, Politekhnicheskaya Str. 29, St.-Petersburg, Russia; e-mail: efed1@yandex.ru
© Санкт-Петербургский государственный политехнический университет, 2013