Тепловой расчёт реактора — основные формулы и порядок расчёта

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

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

Что рассчитывают

В базовом тепловом расчёте определяют:

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

Для реального проекта расчёт выполняется специализированными кодами, но основные формулы полезно понимать вручную.

Тепловая мощность реактора

Тепловая мощность реактора — это количество тепла, выделяемое в активной зоне за единицу времени:

``text Q = N_f · E_f ``

где:

Q — тепловая мощность, Вт;
N_f — число делений в секунду;
E_f — энергия одного деления.

Для урана-235 энергия деления примерно:

``text E_f ≈ 200 МэВ ≈ 3,2 · 10^-11 Дж ``

Если известна тепловая мощность, можно оценить скорость делений:

``text N_f = Q / E_f ``

Например, при мощности 3000 МВт(тепл.) число делений составляет порядка 10^20 делений в секунду.

Связь тепловой и электрической мощности

АЭС превращает тепло в электричество не полностью. Часть энергии теряется в цикле турбины, конденсаторе и вспомогательных системах.

КПД энергоблока:

``text η = N_эл / Q_т ``

где:

η — КПД;
N_эл — электрическая мощность;
Q_т — тепловая мощность реактора.

Для водо-водяных реакторов КПД обычно около 0,32–0,36. Поэтому блок с электрической мощностью 1000 МВт имеет тепловую мощность примерно 3000 МВт.

Тепловой баланс активной зоны

Главное уравнение теплового расчёта:

``text Q = G · c_p · (t_вых - t_вх) ``

где:

Q — тепловая мощность, Вт;
G — массовый расход теплоносителя, кг/с;
c_p — теплоёмкость теплоносителя, Дж/(кг·К);
t_вых — температура теплоносителя на выходе;
t_вх — температура теплоносителя на входе.

Отсюда расход:

``text G = Q / (c_p · Δt) ``

где:

``text Δt = t_вых - t_вх ``

Это первое уравнение, которое используют при оценке расхода теплоносителя через активную зону.

Средний тепловой поток

Если известна поверхность теплообмена твэлов, можно оценить средний тепловой поток:

``text q''_ср = Q / F ``

где:

q''_ср — средний тепловой поток, Вт/м²;
F — суммарная поверхность теплоотдачи твэлов, м².

Для цилиндрического твэла поверхность боковой теплоотдачи:

``text F_твэл = π · d · L ``

где:

d — наружный диаметр оболочки;
L — активная длина твэла.

Если твэлов n, то:

``text F = n · π · d · L ``

Линейная тепловая нагрузка твэла

Линейная тепловая нагрузка показывает, сколько тепла выделяется на одном метре длины твэла:

``text q_l = Q_твэл / L ``

или для всей активной зоны:

``text q_l,ср = Q / (n · L) ``

где:

q_l — линейная нагрузка, Вт/м;
Q_твэл — мощность одного твэла;
n — число твэлов;
L — активная длина.

Это важный параметр, потому что температура топлива сильно зависит от линейной нагрузки.

Коэффициент неравномерности мощности

Мощность в активной зоне распределена неравномерно. В одних каналах тепловыделение выше, в других ниже. Поэтому вводят коэффициенты неравномерности:

``text K_q = q_max / q_ср ``

где:

q_max — максимальное локальное тепловыделение;
q_ср — среднее тепловыделение.

Максимальная линейная нагрузка:

``text q_l,max = K_q · q_l,ср ``

Именно максимальные, а не средние значения определяют запас до перегрева и кризиса теплообмена.

Перепад температуры через оболочку

Тепло от топлива проходит через оболочку твэла к теплоносителю. В простом приближении перепад температуры на теплоотдаче:

``text ΔT = q'' / α ``

где:

ΔT — разность температур между стенкой и теплоносителем;
q'' — тепловой поток;
α — коэффициент теплоотдачи.

Тогда температура наружной поверхности оболочки:

``text T_ст = T_ж + q'' / α ``

где T_ж — температура теплоносителя в данном сечении.

Температура в топливе

В топливной таблетке тепло выделяется по объёму и идёт к поверхности. Для цилиндрического топлива в упрощённой оценке перепад между центром и поверхностью зависит от теплопроводности:

``text ΔT_топл ∼ q_l / (4πλ) ``

где:

q_l — линейная тепловая нагрузка;
λ — теплопроводность топлива.

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

Кризис теплообмена

Кризис теплообмена — это режим, при котором резко ухудшается отвод тепла от поверхности твэла. Для водоохлаждаемых реакторов это опасное состояние: температура оболочки может быстро вырасти.

Запас до кризиса часто описывают отношением:

``text DNBR = q''_кр / q''_факт ``

где:

q''_кр — критический тепловой поток;
q''_факт — фактический локальный тепловой поток.

Если DNBR слишком мал, запас по теплообмену недостаточен.

Гидравлическая часть

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

``text Δp = ξ · ρ · w² / 2 ``

где:

Δp — потери давления;
ξ — коэффициент сопротивления;
ρ — плотность теплоносителя;
w — скорость потока.

Для канала активной зоны также используют число Рейнольдса:

``text Re = ρ · w · d_г / μ ``

и число Нуссельта:

``text Nu = α · d_г / λ ``

Эти критерии позволяют оценить коэффициент теплоотдачи.

Упрощённый порядок расчёта

Типовой учебный алгоритм:

Задать тепловую мощность реактора Q.
Задать температуры теплоносителя на входе и выходе.
Найти массовый расход G.
Определить суммарную поверхность теплообмена твэлов.
Найти средний тепловой поток.
Найти среднюю и максимальную линейную нагрузку.
Оценить температуру стенки оболочки.
Оценить температуру топлива.
Проверить запас до кризиса теплообмена.
Оценить потери давления и требуемую работу насосов.

Пример оценки расхода

Пусть тепловая мощность реактора:

``text Q = 3000 МВт = 3 · 10^9 Вт ``

Подогрев воды в активной зоне:

``text Δt = 30 К ``

Теплоёмкость воды:

``text c_p ≈ 5000 Дж/(кг·К) ``

Тогда расход:

``text G = 3 · 10^9 / (5000 · 30) ≈ 20000 кг/с ``

Это грубая оценка, но она показывает масштаб: через активную зону проходят десятки тысяч килограммов теплоносителя в секунду.

Что важно помнить

Тепловой расчёт реактора — это не только формула Q = G c_p Δt. Главная задача — убедиться, что тепло отводится с достаточным запасом во всех точках активной зоны, включая наиболее нагруженные каналы.

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

Коротко

Основные формулы теплового расчёта:

``text Q = G · c_p · Δt q'' = Q / F q_l = Q / (n · L) K_q = q_max / q_ср DNBR = q''_кр / q''_факт ``

Эти выражения дают основу для понимания того, как реакторная мощность связана с расходом теплоносителя, температурой топлива и запасами безопасности.