Способ определения максимальной температуры окисного топлива в сердечнике вентилируемого твэла термоэмиссионного элемента
Использование: атомная энергетика. Сущность изобретения: при испытаниях электрогенерирующего канала в реакторе измеряют тепловыделение (Q) в твэле и тепловыделение (q) конденсата топлива, вышедшего из твэла, в зависимости от времени, определяют по этим зависимостям для интересующего момента времени скорости изменения dQ/d и dq/d, измеряют давление P парогазовой смеси в межэлектродном зазоре ЭГЭ, а оценку максимальной температуры топлива в сердечнике твэла проводят по выражению Mо - начальная масса топлива в твэле; R - суммарное сопротивление газоотводного устройства; A и B - коэффициенты, зависящие от вида топливного материала; P - давление парогазовой смеси. 2 ил.
Изобретение относится к атомной энергетике, к созданию и наземной обработке твэлов, в частности электрогенерирующих элементов (ЭГЭ) термоэмиссионных электрогенерирующих каналов (ЭГК).
Главным этапом разработки вентилируемых твэлов являются петлевые испытания в исследовательских реакторах, где изучаются все специфические проблемы, связанные с созданием длительно работающих ЭГК. Одной из важнейших характеристик твэлов при испытаниях в исследовательском реакторе является максимальная температура (Т) топливного материала (ТМ), во многом определяющая работоспособность, стабильность и воспроизводимость характеристик ЭГК, свеллинговые процессы в твэлах и т.д. Известны косвенные методы определения Т, например, пользуясь значениями температуры на оболочке ЭГЭ и тепловыделения в твэле, можно рассчитать температурные поля в топливном сердечнике ЭГЭ [1] Наиболее близким к изобретению по технической сущности является способ определения максимальной температуры ТМ, описанный в [2] Он включает измерение тепловыделения qv и TЕ и оценки Т при известных и Rт по соотношению (1) для сплошного топливного цилиндра с источниками тепла: qv плотность объемного тепловыделения в твэле; Rт радиус топливного цилиндра; теплопроводность ТМ; TЕ температура эмиттера. Температуру эмиттера (ТЕ можно определить косвенными методами, например реперной точки, сравнения вольт-амперных характеристик для ЭГК, теплового баланса [3] Точность определения Т по выражению (1) зависит от точности определения l и ТЕ, которые во время испытаний определяются или известны с большой погрешностью. Так, теплопроводность l ТМ, используемая в соотношении (1), относится к структурно-чувствительным свойствам и, значит, l в зависимости от технологии изготовления, режимов работы и т.п. может отличаться в несколько раз от справочного значения [4] Аналогичным образом высока погрешность определения и ТЕ. Как показывают расчетные исследования [3] при определении ТЕ методом теплового баланса температура эмиттера зависит от приведенной степени черноты материалов эмиттера-коллектора, которая известна с большой погрешностью. В результате погрешность определения Т по (1) достаточно велика. Положительным результатом, получаем при использовании изобретения, является повышение точности определения максимальной температуры ТМ, поскольку не требуется предварительного определения температуры эмиттера ТЕ и результаты не зависят от исходного и рабочего состояния ТМ и теплофизических свойств материалов эмиттера-коллектора. Указанный технически результат достигается способом определения максимальной температуры (Т) окисного топлива в сердечнике вентилируемого твэла термоэмиссионного элемента, включающим измерение рабочих значений тепловой мощности и оценку максимальной температуры топлива в твэле, отличающимся тем, что в процессе испытаний измеряют зависимости тепловыделения (Q) в твэле и тепловыделения (q) конденсата топлива, вышедшего из твэла, от времени, определяют по этим зависимостям для интересующего момента времени t, скорости изменения dQ/d и dq/d,, измеряют давление (P) парогазовой смеси с межэлектродном зазоре ЭГЭ, а оценку максимальной температуры топлива в сердечнике твэла проводят по выражениюMo начальная масса топлива в твэле [г]
R суммарное сопротивление газоотводного устройства [1/мм]
A и B коэффициенты, зависящие от вида топливного материала, для UO2 A 1,681032, B 7,4108;
P давление парогазовой смеси [Па] [c]
На фиг. 1 представлен общий вид ЭГЭ, в котором может быть реализован данный способ определения максимальной температуры ТМ в сердечнике вентилируемого твэла, а на фиг. 2 график, поясняющий суть способа. На фиг. 1 обозначено 1 конденсат топлива, вышедшего из твэла, 2 - система вентиляции, 3 топливо, 4 эмиттерная оболочка, 5 коллектор ЭГЭ, 6 коллекторная изоляция, 7 чехловая труба ЭГК, 8 датчик тепловой мощности в твэле, 9 датчик тепловой мощности, выделяемой в конденсате топлива, 10 - изоляция, 11 жиклер. На фиг. 2 обозначено: NA характерная кривая вывода реактора на мощность; Q кривая изменения тепловой мощности в твэле; q кривая изменения тепловой мощности, выделяемой в конденсате топлива. Твэл (см. фиг. 1) содержит оболочку 4, заключающую ТМ 3, и систему вентиляции, выполненную в виде трубки 2 с капиллярным наконечником (жиклером) 11. Способ реализуется следующим образом. ЭКГ, представляющий собой последовательно соединенную сборку ЭГЭ, с системами регистрации 8, 9 тепловой мощности, выделяемой в твэле и конденсате ТМ 1, вышедшего из твэла, помещают в ячейку реактора. В процессе вывода реактора на мощность NA (фиг. 2) и затем в процессе работы реактора измеряют давление P парогазовой смеси в межэлектродном зазоре ЭГЭ, например, с помощью датчика давления. Измеряют изменение во времени тепловой мощности, выделяемой в твэле Q(), и тепловой мощности q(), выделяемой в конденсате ТМ. Для интересующего момента времени , определив по зависимостям Q() и q() (типа приведенных на фиг. 2) скорость изменения dQ/d и dq/d,, регистрируя суммарное давление парогазовой смеси в межэлектродном промежутке ЭГЭ (P), при известной начальной массе топлива в твэле Мо по выражению (2) оцениваем максимальную температуру в твэле. Выражение (2) получено из выражения для потока молекул ТМ (m) в случае диффузии их в среде паров Сs и газообразных продуктов деления из топливного сердечника ЭГЭ через систему вентиляции, используя первый закон Фика:
где w средняя скорость теплового движения молекул ТМ;
L средняя длина свободного пробега молекул ТМ в парогазовой смеси;
no концентрация пара ТМ в сердечнике ЭГЭ;
mтм масса молекулы ТМ;
R сопротивление газоотводного устройства.
где l1 и r1 длина и радиус отдельных участков газоотводного устройства соответственно. Выражение для no можно получить из зависимости для стехиометричной UO2 для равновесия между паровой и адсорбированной фазой [5]
lg P [мм рт.ст. 32258/Т + 12,183. (4)
Поток молекул ТМ (m) через систему вентиляции твэла можно определить также из соотношения
m = dv/d, (5)
где v объем конденсата топлива, вышедшего из ЭГЭ. Тепловыделение конденсата топлива определим из соотношения
q vqv, (6)
qv Q/V q/v (q + Q)/Vo, (7)
Vo и V начальный объем ТМ в твэле и текущее значение объема ТМ в твэле соответственно. Продифференцировав (6) и используя (7), получим выражение для dv/d, откуда после небольших преобразований (5) получим выражение для m:
Используя известные выражения для L и w из [6] и no из формулы (4), по соотношениям (3) и (8) получим выражение (2) для Т. В качестве примера рассмотрим использование способа определения Т для типичного ЭГЭ, где в качестве ТМ используется двуокись урана (UO2). Основное сопротивление газоотводного устройства 2 (фиг. 1) оказывает жиклер 11, сопротивление которого R l/r2. Здесь l длина жиклера, r радиус жиклера. Для типичного ЭГЭ примем R 1600 1/мм. Давление парогазовой смеси, регистрируемое в межэлектродном зазоре ЭГЭ (обычно 130-1500 Па), примем P 1330 Па. Из анализа имеющихся экспериментальных данных по испытаниям энергонапряженных ЭГК примем для типичных режимов испытаний для момента времени тепловыделение в твэле Q 900 Вт и тепловыделение в конденсате ТМ q 0,2 Вт. Определяем по измеренной зависимости dq/d = 710-8 Вт/с.. Зная первоначальную массу топлива в твэле Мo 33 г, определяем по выражению (2) Т 2511,3 К. В данном примере уравнение (2) решалось итерационным методом [7]
Можно решать уравнение (2) и графически с помощью обратной функции. Выражение (2) преобразуем к виду
где m определяется выражением (8). Построим зависимость m f(Т) для заданного R и P по выражению (9). Определив m по выражению (8), из графической зависимости m f(T) находим Т. Литература
1. Корнилов В.А. и др. Метод расчета температурных полей гетерогенного топливного сердечника термоэмиссионного электрогенерирующего элемента. Атомная энергия, 1980, т.49, вып.6, с.393-394. 2. Займовский А.С. и др. Тепловыделяющие элементы атомных реакторов, М. Атомиздат, 1966г, с.504. 3. Синявский В.В. Особенности определения температуры эмиттера при испытаниях термоэмиссионных преобразователей. ТВТ, т. 12, N6, с.1267-1271, 1974. 4. Колядин В.И. и др. Теплопроводность двуокиси урана. Атомная энергия, т.36, вып.1, 1974г, с.59-60. 5. Горбань Ю.А. и др. Исследование испарения двуокиси и карбидов урана. Атомная энергия, 1967, т.22, вып.6, с.465-467. 6. С. Дэшман. Научные основы вакуумной техники. Мир, М. 1964, с.68 и с. 84. 7. Г. Корн и Т.Корн. Справочник по математике для научных работников и инженеров. Издание второе, Москва, 1970, с.572.
Формула изобретения
где М0 начальная масса топлива в ТВЭЛе, г;
R суммарное сопротивление газоотводного устройства, мм-1;
A и B коэффициенты, зависящие от вида топливного материала, для UO2 A 1,68 1032, B 7,4 108;
Р давление парогазовой смеси, Па.
РИСУНКИ
Рисунок 1, Рисунок 2