Ядерные энергетические реакторы
.pdf
Количество ядер |
урана |
|
в |
единице объема определяется в |
||
соответствии с заданной величиной “a” по формуле: |
||||||
|
N |
|
|
N |
3 |
, |
|
|
|||||
|
|
|
|
|
||
|
|
|
а |
|
||
|
|
|
|
|
||
где а – количество ядер замедлителя, приходящихся на одно ядро урана. Уран в топливной композиции состоит из двух изотопов: урана-235
(N5) и урана-238 (N8). Следовательно, количество каждого изотопа будет зависеть от обогащения урана-238 ураном-235 (Х).
Х |
N |
5 |
100 % ; |
|
|||
N |
|
||
|
|
|
Х
N5 100 N ,
N8 N N5.
Коэффициент использования тепловых нейтронов рассчитывается по
формуле: |
|
|
|
|
|
|
|
|
|
|
а |
|
|
|
|
1 |
|
|
, |
|
|
|
|
а3 |
|
||||
|
|
1 |
|
N3 |
|||||
|
а |
а3 |
|
|
|
|
|
||
|
|
а5 |
N5 |
а8 N8 |
|||||
|
|
|
|
|
|||||
где σа3 – микроскопическое |
эффективное сечение поглощения ядер |
||||||||
замедлителя; σа5 – микроскопическое эффективное сечение поглощения ядер урана-235;
σа8 – микроскопическое эффективное сечение поглощения ядер урана-238; σа3, σа5, σа8 – принимаются из таблицы (Приложение 1).
Количество быстрых нейтронов деления η, приходящихся на одно поглощение теплового нейтрона ураном, находится по зависимости:
2,47 |
|
f 5 |
|
, |
|
а5 |
а8 |
N8 |
|
||
|
N5 |
|
|
||
где: σf5 – микроскопическое эффективное сечение деления урана-235 принимается из табл. 4.
Определение величины θ и η ведется для трех значений обогащения урана Х. Поэтому необходимо предварительно задаться этими значениями в пределах от Х = 2 % до Х= 5 %. Три значения величины в указанном интервале студент принимает самостоятельно (например 2, 3, 5 или 2,5; 3,5; 4,5).
Таким образом, имеющиеся значения величин ε, φ, θ, η позволяют рассчитать коэффициент К∞ по формуле (4.3) для трех значений обогащений урана.
191
–Значение возраста нейтронов τ, входящее в уравнение (4.1), принимается по замедлителю нейтронов τ = τзам.
–Квадрат длины диффузий нейтронов L2 определяется по формуле:
2 |
2 |
1 , |
|
L |
L |
зам |
|
|
|
|
|
где: Lзам – длина диффузии нейтронов для принятого в расчете замедлителя.
θ– коэффициент использования тепловых нейтронов.
–Имея значения К∞,B2 , τ и L2 по формуле (4.1) рассчитывается
величина эффективного коэффициента размножения нейтронов в реактора. Результат этого расчета представляется в виде графика К
активной зоне |
||
эф |
f Х |
. |
|
|
|
Пересечение горизонтальной линии, соответствующей величине К 0эф (см. табл. 1) с полученной зависимостью Кэф = f(Х) покажет необходимое значение Хо.
– В заключение следует установить влияние отражателя нейтронов на величину обогащения Хо. Для этого необходимо рассчитать значение Кэф при другой величине геометрического параметра В2. В этом случае в формуле (1) будет изменяться только величина В2, которую необходимо рассчитать с учетом отражателя нейтронов:
|
|
2,4 2 |
|
2 |
||
|
В2 |
|
|
|
|
, |
|
|
|
||||
|
|
R |
|
Н |
||
где R' = R + Δ, H' = H + 2Δ – эффективные радиусы и высота активной зоны. |
||||||
Значение эффективной добавки принимается самостоятельно для |
||||||
принятого замедлителя нейтронов: |
|
|
|
|
|
|
- для воды – = 8 – 12,см ; |
|
|
|
|
|
|
- для тяжелой воды – = 20 – 40, см ; |
|
|
|
|
|
|
- для графита – |
= 80 – 120, см; |
|
|
|
|
|
- для бериллия – |
= 150 – 250, см. |
|
|
|
|
|
Зависимость Кэф = f(Х) для новой величины В2 пойдет несколько выше прежней. Из кривых будет видно, что наличие отражателя нейтронов снижает величину обогащения для требуемой величины К0эф.
КОНТРОЛЬНАЯ РАБОТА № 2
Определение размеров реактора активной зоны реактора
– Определяется тепловая мощность атомной энергетической установки
Qт Qэл. 1 ,
у
192
где: ηу – коэффициент полезного действия атомной энергетической установки, %.
– Оцениваются параметры теплоносителя.
По давлению Р1 из таблицы табл. 5 находится температура кипения теплоносителя ТI', К.
Принимается температура теплоносителя на выходе из реактора:
Т |
вых |
Т |
Т |
нк |
, |
|
I |
|
|
Температура теплоносителя на входе а активную зону находится по зависимости:
Т |
вых |
Т |
вых |
Т |
р |
. |
|
|
|
|
По определенной тепловой мощности реактора Qт и параметрам теплоносителя определяется расход теплоносителя
G т |
|
Qт |
|
|
|
|
, |
|
Ср т Т |
вых |
Т |
вх |
|
||||
|
|
|||||||
|
|
|
|
|
|
|||
где Qт – тепловая мощность реактора, кВт;
Срт – теплоемкость теплоносителя, принимается по его средней температуре.
– Определяется сечение активной зоны для прохода теплоносителя
Fт |
Ст т |
|
, |
W т |
|
||
|
ср |
||
|
ср |
|
|
где: γтср, кг/м3 - плотность теплоносителя (таблица) определяется по его средней температуре.
– Выбор тепловыделяющей сборки (ТВС).
Приводится обоснование варианта ТВС, принятой для расчета. Перед выбором ТВС необходимо предварительно ознакомиться с тепловыделяющими сборками действующих реакторов. Обратить внимание на величину диаметра
193
тепловыделяющего элемента (ТВЭЛ), количество ТВЭЛ в одной ТВС. Руководствоваться при этом следует имеющимися рекомендациями. В частности диаметр ТВЭЛ для реакторов ВВЭР можно принимать в интервале dТВЭЛ = (6 ÷ 10)·10-3 м, для реакторов РБМК dТВЭЛ = (12 ÷ 15)·10-3 м .
Необходимо выполнить в масштабе эскиз и разрез тепловыделяющей сборки и ТВЭЛ с соблюдением подлинных размеров. Вычислить площадь сборки fсб. Основным параметром обоснования варианта ТВС должно быть значение для прохода теплоносителя в одной ТВС → fт. Например, если сечение сборки по внутреннему размеру корпуса ТВС обозначить fсб, то сечение для прохода теплоносителя одной сборки найдется по формуле:
|
|
f |
|
n |
|
Пd |
2 |
|
f |
|
|
|
ТВЭЛ |
, |
|||
т |
сб |
ТВЭЛ |
|
|||||
|
|
|
|
|
||||
|
|
|
|
4 |
|
|||
|
|
|
|
|
|
|
|
где: nТВЭЛ – число тепловыделяющих элементов в ТВС; dТВЭЛ – наружный диаметр ТВЭЛ.
Число тепловыделяющих сборок в активной зоне реактора определяется по формуле:
n ТВС Fт .
f т
Сечение одной ячейки Sяч., соответствующей ТВС реактора ВВЭР может быть принято по зависимости:
Sяч fсб fзам ,
где fзам – сечение замедлителя (воды), находящегося снаружи тепловыделяющейся сборки.
Обозначим длину периметра ТВС тепловыделяющейся сборки Псб, тогда сечение замедлителя найдется
fзам Псб зам ,
где: σзам – толщина замедлителя вокруг сборки, м.
194
Толщина замедлителя для реакторов ВВЭР принимается самостоятельно в интервале:
|
|
10 15 10 |
3 |
зам |
. |
||
|
|
|
Площадь сечения активной зоны найдется по зависимости:
S |
1,15 n |
ТВС |
S |
яч |
, |
аз |
|
|
|
где коэффициент 1,15 введен с целью увеличения Sаз для размещения приборов контроля и регулирующих стержней.
Диаметр активной зоны определяется из равенства:
|
|
|
|
|
ПД |
2 |
|
|
|
|
S |
|
|
аз |
, |
||||
|
аз |
|
|
||||||
|
|
|
|
|
|
|
|||
|
|
|
4 |
|
|
|
|||
|
|
|
|
|
|
|
|
||
тогда |
|
|
|
|
|
|
|
|
|
Д |
|
|
|
|
4 Sаз |
. |
|||
аз |
|
|
|||||||
|
|
|
|
П |
|
|
|||
|
|
|
|
|
|
|
|
||
Высота активной зоны принимается в зависимости от диаметра
Наз 0,8 1,2 Даз .
Первая часть контрольной работы считается выполненной, так как определены размеры активной зоны Даз, Наз.
Расчет температурного режима наиболее теплонапряженного тепловыделяющего элемента
Энерговыделение в активной зоне реактора соответствует изменению нейтронного потока. Принимаем синусоидальный закон изменения тепловыделения по длине ТВЭЛ
q v qmaxv sin ПHh ,
195
где qvmax – максимальное, удельное тепловыделение, приходящее на единицу объема горючего, кВт/м3;
h – текущее значение длины ТВЭЛ, м; Н' – эффективная длина ТВЭЛ, м.
Н
Н
2
.
Расчетная схема температурного режима ТВЭЛ показана на рис. 38, а. Расчет температурного режима ТВЭЛ ведется при использовании
зависимости
qh qmaxh sin ПHh ,
где qhmax – максимальное удельное тепловыделение, приходящееся на единицу
длины ТВЭЛ, кВт/м.
qh – удельное тепловыделение приходящееся на единицу длины для участка при координате h, кВт/м.
196
а) |
б) |
Рис. 38. Расчетная схема температурного режима ТВЭЛ |
|
Максимальный удельный тепловой поток определяется по зависимости.
qmaxh qсрh Кh ,
где Кh – коэффициент неравномерности тепловыделения по высоте активной
зоны. При расчете Кh – принять самостоятельно в интервале Кh = = 2,2 ÷ 2,5
Среднее значение удельного теплового потока, приходящегося на единицу длины ТВЭЛ зависит от мощности реактора и длины ТВЭЛ, определяется по формуле:
197
q |
ср |
|
|
Q |
т |
, |
|
|
|||||
h |
Н n |
|
n |
|||
|
|
|
|
|||
|
|
|
ТВС |
ТВЭЛ |
||
|
|
|
|
|
где Н – высота (длина) ТВЭЛ, м;
nТВС – число тепловыделяющих сборок в активной зоне; nТВЭЛ – число ТВЭЛ в сборке.
Использование уравнения позволяет рассчитать изменение температур по длине ТВЭЛ: теплоносителя, оболочки ТВЭЛ со стороны теплоносителя, оболочки ТВЭЛ со стороны горючего и самого горючего по формуле:
qh
где gт – расход теплоносителя,
dh gтСрт dt ,
приходящийся на один ТВЭЛ.
Его значение определяется по зависимости
gт |
Gr |
Кr, |
|
nТВС nТВЭЛ |
|||
|
|
где Gr – расход теплоносителя на реактор; Ср – теплоемкость теплоносителя;
Кr – коэффициент неравномерности распределения теплоносителя по радиусу реактора Кr 1,45.
Определение температуры теплоносителя. Решение уравнения (30) будет
dТ |
|
qhmax |
|
|
sin |
Пh |
dh |
, |
|
|
||||||||
|
g т Ср |
т |
|
H |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Н1 |
|
qhmax |
|
|
|
|
|
|
|
Пh |
|
|
|
|||||
t |
|
|
|
|
|
sin |
|
|
|
|
|
dh , |
|
|||||
|
|
т |
|
|
|
H |
|
|
||||||||||
gт С |
|
р |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
П |
|
|
|
Пh |
|
||||||
t т tвх |
А cos |
|
|
|
|
|
cos |
|
|
, |
||||||||
Н |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
H |
|
|||||||
где Тт – температура теплоносителя по высоте h, К;
198
Твх – температура теплоносителя на входе в реактор, К;– эффективная добавка, обусловленная влиянием отражателя нейтронов.
А |
q |
max |
H |
, |
|||
|
|
|
|||||
|
h |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
Пg |
|
Ср |
т |
|
||
|
т |
|
|
||||
|
|
|
|
|
|
|
|
где А – постоянный комплекс величин, характеризующих тепловую нагрузку,
эффективную высоту ТВЭЛ и расход теплоносителя.
Уравнение позволяет определить температуру теплоносителя в любом сечении по длине ТВЭЛ.
Температуру теплоносителя на выходе из тепловыделяющей сборки найдем по уравнению, если вместо h подставить значение Н′-
|
вых |
t |
|
|
2 |
|
q |
max |
H |
cos |
П |
t |
|
|
Q |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
t |
|
|
|
h |
|
|
H |
|
|
|
твэл |
. |
|||||||||
т |
вх |
П |
Ст |
т |
Ст |
вх |
Ст |
т |
Ст |
||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
– Определение температуры оболочки ТВЭЛ со стороны теплоносителя.
Температура стенки Т'cт, со стороны теплоносителя определяется на
основании решения уравнения теплоотдачи. |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
max |
|
Пh |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
sin |
|
|
|
|
|
|
|||||
|
2 П rст tст t т qh |
Н |
2 П rст tст t т |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Пh |
|
|
П |
|
Пh |
|
Пh |
|
|||
|
|
|
tвх А cos |
|
|
|
cos |
|
|
|
, |
||||
|
tст t т В sin |
H |
H |
|
В sin |
H |
|||||||||
|
|
|
|
|
|
|
H |
|
|
||||||
|
qhmax |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где В |
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2П rст |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где r'ст – наружный периметр ТВЭЛ, м.
α – коэффициент теплоотдачи оболочки ТВЭЛ к теплоносителю, кВт/м2 К. Как следует из (4.36) определение величины В' может быть
выполнено, если известен коэффициент теплоотдачи α. Расчет коэффициента α делается по известной зависимости конвективного теплообмена. В данном случае имеет место теплоотдача стенки оболочки ТВЭЛ к однофазной среде (воде). При движении теплоносителя внутри труб рекомендуется формула:
199
N |
0.021Re |
0.8 |
Pr |
0.43 |
, |
|
|
||||
u |
|
|
|
|
|
где критерий Нуссельта (Nu), Рейнольдса (Re) и Прандтля (Pr) определяются по следующим формулам:
Nu
d
;
Re
Wd
;
Pr
.
Значение коэффициента теплопроводности α, кинематической вязкости ν и температуропроводности α принимаются из таблиц по средней температуре теплоносителя.
Значения эквивалентного диаметра α, входящего в критерии Nu и Re для тесных решеток и тепловыделяющих сборок можно принимать как величину наружного диаметра тепловыделяющего элемента.
Скорость теплоносителя можно определить по зависимости:
Wgт nТВЭЛ ,
ср fт
где gт – массовый расход теплоносителя, приходящийся на один тепловыделяющий элемент, кг/с;
γср – плотность теплоносителя;
fт – сечение для прохода теплоносителя одной сборки; nтвэл – число ТВЭЛ в сборке.
После расчета критериев определяется коэффициент теплоотдачи
0,021Re0.8 Pr0.4 |
. |
|
|
Температура оболочки ТВЭЛ |
Т'ст меняется по длине. |
Максимальная температура Т'стmax устанавливается на определенной длине h'о. Значение этой длины находится методом решения на максимум уравнения (35),
|
|
|
|
|
|
|
|
|
полагая |
dTст |
0. В результате этого решения получим |
||||||
dh |
||||||||
|
|
|
|
|
|
|
||
|
|
hо |
H |
|
|
А |
||
|
|
|
arcctg |
|
|
. |
||
|
|
П |
|
|||||
|
|
|
|
|
В |
|
||
Подставляя значение h'о в уравнение, будем иметь
200
