Процесс развития проектной аварии с блокировкой проходного сечения отдельной ТВС при работе реактора на быстрых нейтронах с натриевым теплоносителем (БН) на мощности представляет интерес с точки зрения распространения повреждений на соседние ТВС и перемещения вниз расплавленных компонентов. Анализ распространения аварии в радиальном направлении важен для разработки в проекте реакторной установки мер по предотвращению повреждения тепловыделяющих элементов (твэлов) в объеме за пределами аварийной ТВС. Важны также последствия такой аварии для нижележащих конструкций – будет ли расплав удержан внутри аварийных ТВС или проплавит их до хвостовиков и выйдет на плиту напорной камеры, причинив повреждения и плите. Для предотвращения развития аварии в отдельной ТВС за ее пределы в активной зоне и соответственно перехода аварии в разряд запроектных необходима своевременная остановка реактора средствами воздействия на реактивность (аварийной защитой). Чтобы избежать распространения аварии за пределы аварийной ТВС, реактор должен быть остановлен до достижения кипящего бассейна в этой сборке [6]. Для доказательства эффективности средств обнаружения аварий и аварийной защиты – доказательства своевременности срабатывания аварийной защиты необходим детальный анализ аварии, в результате которого должны быть определены временные параметры ее развития.
Цель исследования
Целью анализа является обоснование эффективности третьего уровня глубокоэшелонированной защиты, связанного с предотвращением развития проектной аварии в запроектную [4] – распространения повреждения и расплавления твэлов от одной аварийной ТВС на окружающие ТВС активной зоны сверх проектного предела повреждения твэлов. Проектный предел повреждения твэлов для проектной аварии с блокировкой отдельной ТВС для реактора БН состоит в разрушении твэлов аварийной ТВС с непревышением пределов безопасной эксплуатации по повреждению твэлов во всем остальном объеме активной зоны.
Методика исследования
Для численного моделирования процессов, необходимость рассмотрения которых возникает при анализе аварий с плавлением материалов активной зоны в быстрых реакторах, был разработан комплекс расчетных моделей [1]. Расчетные модели для анализа проектной аварии:
– модель разогрева твэла, плавления оболочки и топливного сердечника с приближенным учетом выпаривания натрия и с учетом перемещения расплавленной оболочки при прекращении расхода теплоносителя в ТВС быстрого реактора, реализованная в программах TRAMS-TWEL и MOTE [2];
– модель проплавления чехлов аварийной и окружающих ТВС расплавленным топливом при проектной аварии с блокировкой отдельной ТВС быстрого реактора, реализованная в программе TRAMS-DUCT [3];
– модель опускного продвижения расплавленной тепловыделяющей массы из активной зоны в нижнюю часть корпуса быстрого реактора при авариях с плавлением ТВС в активной зоне, реализованная в программе TRAMS‑FDOWN [8].
В результате расчетов по программе TRAMS-TWEL для ряда высотных участков определяются времена выпаривания натрия по высоте отдельного твэла, расплавления оболочки и топливного сердечника, полной или частичной блокировки межтвэльного пространства затвердевающим материалом, а также повторного расплавления слоя затвердевшего материала (твердой корки) за счет теплопереноса в радиальном направлении. Времена начала и окончания перемещения расплавленного материала поверх твердой оболочки твэла в рассчитываемых высотных участках выше области расплавления оболочки, начала такого перемещения ниже области расплавления и заполнения расплавом всего сечения ячейки в нижней части активной зоны определяются в результате расчетов аксиального перемещения расплавленного материала оболочки по программе MOTE. По данной программе определяются также время смены восходящего движения на нисходящее, длина проникновения расплавленного материала, возможно, в верхнюю и в нижнюю зоны воспроизводства до момента полной блокировки последней, а затем время образования и высота слоя из расплавленного материала оболочки над блокировкой каналов под активной зоной. Расчетная модель теплообмена в твэле с твердой коркой связана по итерационной схеме с моделью перемещения расплавленной оболочки. В результате расчетов по программе TRAMS-DUCT определяются времена полного проплавления либо глубина частичного проплавления стенок чехлов аварийной и соседней ТВС. В результате расчетов по программе TRAMS‑FDOWN определяются глубина и время опускного продвижения расплавленной тепловыделяющей массы из активной зоны через нижнюю зону воспроизводства и последующие конструкции с учетом отвода тепла от тепловыделяющей массы в радиальном направлении к соседним неповрежденным ТВС.
Результаты исследования
В результате расчетов по программе TRAMS-TWEL получено, что температура натрия достигает точки насыщения сначала в центре по высоте активной зоны примерно через 0,7 с после прекращения расхода в ТВС, далее зона кипения расширяется вверх и вниз. Расплавление оболочки твэла наступает первоначально в центре через 3,1 с после прекращения расхода.
Восходящее движение увлекаемого парами натрия, испаряющегося в нижней части активной зоны, жидкого материала оболочки начинается с момента ее расплавления на наиболее горячем участке. В момент 4,4 с происходит изменение направления движения расплава. При этом верхняя граница стальной пленки достигает высоты z = 78,6 см, далеко не доходя до верха активной зоны. Следовательно, блокировок каналов для прохода теплоносителя затвердевающим расплавом оболочек в верхней части активной зоны и в области необогреваемого участка твэла не произойдет. При этом источники запаздывающих нейтронов будут выноситься к выходу из аварийной ТВС непосредственно после разрыва либо расплавления оболочек твэлов.
Нижняя граница движущейся стальной пленки начинает опережать границу зоны расплавления оболочки с высоты z » 13 см. В момент 5,9 с жидкий материал достигает зоны воспроизводства и затвердевает. В результате затвердевания происходит полная закупорка проходного сечения в нижней зоне воспроизводства на глубину 1,8 см. Продолжающееся опускное движение расплава в активной зоне приводит ко времени 6,4 с после прекращения расхода к заливу сталью канала над блокировкой на высоту 29,4 см. Перемещение оболочки в данном случае завершается при частично расплавленном материале внутри топливного сердечника. В нижней части активной зоны образование твердой корки поверх твердой оболочки начинается с высоты z » 4 см.
Плавление топливного сердечника в центральной плоскости активной зоны начинается в момент 5,4 с и заканчивается через 7,8 с после прекращения расхода. В момент 16,6 с топливо расплавляется до верха активной зоны.
Ко времени расплавления топлива в центральной области активной зоны на участке, составляющем примерно одну треть ее высоты, в области стального бассейна будет расплавлена лишь малая доля топлива. Следовательно, жидкое топливо из средней и верхней частей активной зоны будет собираться над стальным бассейном при преимущественно твердых пеньках топливных сердечников в нижней части активной зоны. При этом образуется конфигурация, при которой около двух третей топлива собирается в компактную массу в центральной части по высоте активной зоны, вносящую наибольший вклад в реактивность. Такое перемещение расплава может быть использовано для обнаружения аварии системой вычисления баланса реактивности.
Расчеты проплавления стенок чехлов аварийной и окружающих ТВС расплавом топлива проведены по программе TRAMS-DUCT. В результате расчетного анализа выявлены основные параметры, определяющие конечный результат и продолжительность плавления стенок чехлов. Получено, что с момента контакта начинается плавление стенки чехла аварийной ТВС и затвердевание расплавленного топлива у стенки. Через 0,9 с вскипает натрий в межпакетном зазоре. Это приводит к расплавлению стенки чехла аварийной ТВС. Расплавление соседней стенки произойдет при кризисе теплоотдачи к теплоносителю соседней ТВС. Плавление стенки чехла предотвращается срабатыванием системы аварийной защиты (АЗ) до начала ее плавления.
Время проплавления стенки чехла аварийной ТВС тонким неподвижным слоем топлива при сохранении прослойки из расплава стенки в месте плавления составит около 8 с, в случае вытекания расплава плавление стенки происходит быстрее – примерно за 5 с.
Расчеты, проведенные по двум схемам, показали, что время плавления стенок чехлов зависит от массы топлива у стенки и интенсивности теплоотдачи от расплавленного топлива к стенкам. Время проплавления стенки чехла аварийной ТВС первоначально большой неподвижной массой расплавленного топлива составит ~7 с при сохранении прослойки из жидкой стали, а при ее вытекании – ~4 с. Минимальное время проплавления стенки в соответствии со второй схемой составит ~5 с при сохранении жидкой стали, а в случае полного ее оттока – всего ~3 с.
Расчеты опускного продвижения тепловыделяющей массы в аварийной и окружающих ТВС выполнены по программе TRAMS‑FDOWN. Рассматривается медленное продвижение при остаточном тепловыделении, при котором радиальное распространение прекратилось, т.е. предполагаемая конфигурация не может нарушиться. Расчеты показали значительное замедление продвижения тепловыделяющей массы в ограниченном числе ТВС по сравнению с перемещением в значительном числе ТВС, характеризующемся малой долей боковых потерь тепла. Получено, что в семи ТВС тепловыделяющая масса проплавляет лишь около 1/3 части зоны воспроизводства, после чего происходят прекращение движения и охлаждение расплава с сопутствующей конденсацией паров топлива из вышележащего парового объема. Для одной ТВС получено, что при принятых предпосылках тепловыделяющая масса практически не будет плавить нижележащие конструкции. Такой результат для случая одной ТВС качественно согласуется с результатами экспериментов серии BE+ программы SCARABEE [7], в которых не было движения расплава вниз после остановки реактора.
Выполненный анализ проектной аварии с блокировкой отдельной ТВС реактора БН показал, что проектный предел при данной аварии выполняется. Проектный предел повреждения твэлов для проектной аварии с блокировкой отдельной ТВС для реактора БН состоит в разрушении твэлов аварийной ТВС с непревышением пределов безопасной эксплуатации по повреждению твэлов во всем остальном объеме активной зоны. Время срабатывания АЗ по сигналу системы контроля герметичности оболочек твэлов (КГО) после выхода источников из твэлов в соответствии с оценками специалистов составит примерно 5–10 с [5]. Тогда время срабатывания АЗ после расплавления оболочек в центральной плоскости активной зоны составит 8–13 с от момента исходного события аварии (рис. 1). В результате, топливо будет локализовано пределами одной аварийной ТВС – проектным пределом повреждения твэлов. В маловероятном случае слабого сигнала системы КГО после расплавления оболочек твэлов со свежим топливом произойдет срабатывание АЗ по несколько более позднему сигналу системы КГО после выноса источников при плавлении топлива. В случае проплавления стенки чехла аварийной ТВС топливо с остаточным тепловыделением затвердеет на стенке чехла соседней ТВС. В результате, топливо также фактически будет локализовано пределами одной аварийной ТВС. Срабатывание АЗ остановит формирование бассейна расплавленного топлива в аварийной ТВС и, следовательно, проплавление стенок чехлов ТВС. Получение достаточно сильного сигнала системы КГО для срабатывания АЗ должно быть подтверждено экспериментально для штатных ТВС конкретного проекта реактора БН. Тепловыделяющая масса при аварии будет удержана внутри аварийной ТВС и в высотном направлении, не причинив повреждений нижележащим конструкциям, в частности плитам напорной камеры.
Рис. 1. Последовательность событий при проектной аварии с блокировкой отдельной ТВС.
Заключение
Выполнен анализ проектной аварии с блокировкой отдельной ТВС реактора БН в обоснование эффективности его защиты. Анализ показал, что проектный предел при данной аварии выполняется. Срабатывание АЗ по сигналу системы КГО ограничит аварию пределами одной аварийной ТВС – проектным пределом повреждения твэлов.
Рецензенты:
Сухарев Ю.П., д.т.н., главный специалист, Открытое Акционерное Общество «Опытное Конструкторское Бюро Машиностроения им.И.И.Африкантова» (ОАО «ОКБМ Африкантов»), г. Н. Новгород;
Радионов А.А., д.т.н., профессор, зав. кафедрой общей и ядерной физики, Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Нижегородский государственный университет им. Р.Е. Алексеева», г. Н.Новгород.