Определение температур теплоносителя и твэла, запаса до кризиса теплообмена

Исходные данные, полученные в результате вычислений предыдущего раздела, сведены в таблицу 5. Основная задача дальнейших расчетов — определение изменений по высоте канала следующих величин: температур теплоносителя, составных частей твэла (оболочки, топлива) и запаса до кризиса теплобмена. Все вычисления данного раздела производились по методике для ВВЭР, описанной в [2, стр. 322−327]. Значения параметров теплоносителя определены согласно данным [4]. Перед основным расчетом произведены дополнительные геометрические вычисления для ТВС и твэла, отраженные в таблице 6.

Таблица 5 — Параметры теплоносителя на входе и в среднем по активной зоне.

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

Тем не менее, в рамках данной работы влияние закипания теплоносителя в каналах (как объемного в ТВС, так и поверхностного на стенке твэла) не учитывалось. Таким образом:

• ? линейный тепловой поток имеет косинусоидальное распределение по высоте (высотная координата z=0 в центральной плоскости активной зоны)
• ? при достижении точки насыщения расчет коэффициента конвективной теплоотдачи производился по формулам для однофазного теплоносителя, т. е. для воды в состоянии насыщения с нулевым содержанием пара.

Стоит отметить, что для снижения неравномерности энерговыделения в современных реакторных установках типа ВВЭР используются следующие технологические решения:

• ? профилирование топлива по радиусу активной зоны;
• ? профилирование расхода теплоносителя по радиусу;
• ? использование интенсификаторов теплообмена;
• ? использование бесчехловых ТВС.

Распределение линейного теплового потока по высоте среднеи максимально нагруженного твэлов изображено на рисунке 4. На рисунке 5 отражена зависимость энтальпии от высоты, на рисунке 6 — температуры теплоносителя.

Вычисления энтальпии по высоте твэла производилось следующим образом:

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

На рисунке 7, в соответствии с данными [2, стр. 169], отражена использованная в вычислениях зависимость коэффициента теплоотдачи контактного слоя между оболочкой твэла и топливным сердечником в зависимости от величины его первоначального газового зазора.

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

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

Оценка запаса до кризиса теплообмена производилась согласно.

[2, стр. 238] по следующей формуле:

Как видно из рисунка 9, коэффициенты запаса до кризиса теплообмена не опускаются ниже 7, но в виду того, что у формулы (4) рекомендуемые пределы применения не в полной мере соответствуют значениям величин, полученным в данной работе, было дополнительно определено граничное паросодержание перехода от пузырькового кипения к пленочному xгр=0,38 [2, стр. 236]. Данное паросодержание не достигается ни в одной из высотных координат. Таким образом, с достаточной степенью точности был произведен вывод об отсутствии кризиса теплообмена в активной зоне.

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

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

Рисунок 7 — Коэффициент теплоотдачи контактного слоя Таблица 7 — Определение теплогидравлических параметров по высоте.

Наименование параметра.	Формула / обозначение.	Ед. изм.	Значение.
Коэффициент неравномерности энерговыделения по объему.	kV (принимается).	Ї.	2,8.
Коэффициент неравномерности энерговыделения по высоте активной зоны.	Ї.	1,51.
Коэффициент неравномерности энерговыделения по радиусу активной зоны.	Ї.	1,86.
Максимальный линейный тепловой поток (в центральной плоскости активной зоны):	Вт/м.
для средненагруженного твэла.	1,27*104.
для максимально нагруженного твэла.	2,37*104.
Распределение линейного теплового потока.	Вт/м.	см. таблицу 8.
Тепловая нагрузка на единицу поверхности твэла.		см. таблицу 8.
Массовая скорость теплоносителя в ТВС.
Расход теплоносителя в расчете на один твэл.	кг/с.	0,614.
Энтальпия теплоносителя по высоте канала.	см. формулу 3.		см. таблицу 8.
Относительная энтальпия по высоте канала.	Ї.	см. таблицу 8.
Критерий Рейнольдса по высоте канала.	Ї.	см. таблицу 8.
Критерий Прандтля по высоте канала.	Ї.	см. таблицу 8.
Коэффициент A для определения критерия Нуссельта.	Ї.	3,06*10−2.
Критерий Нуссельта по высоте канала (для стержней в треугольной упаковке).	Ї.	см. таблицу 8.
Коэффициент теплоотдачи между теплоносителем и твэлом.		см. таблицу 8.
Температура наружной поверхности оболочки твэла.	C.	см. таблицу 8.
Коэффициент теплопроводности оболочки (цирконий) [2, стр. 320].		20,9.
Температура внутренней поверхности оболочки твэла.	єС.	см. таблицу 8.
Отношение толщины газового зазора к внутреннему диаметру оболочки твэла.	дз/d1.	%.	0,76.
Коэффициент проводимости контактного слоя.	(см. рисунок 7).		4,7*103.
Температура наружной поверхности топливного сердечника.	єС.	см. таблицу 8.
Коэффициент теплопроводности нитрида урана.	где — средняя по толщине температура топлива в К.		см. таблицу 8.
Температура внутренней поверхности топливного сердечника.	єС.	см. таблицу 8.
Критический тепловой поток.		см. таблицу 8.
Коэффициент запаса до кризиса теплообмена.	см. таблицу 8.
Граничное паросодержание (переход от пузырькового кипения к пленочному).	*.	Ї.	0,37.

*.