Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Ядерные энергетические реакторы

.pdf
Скачиваний:
2
Добавлен:
30.11.2025
Размер:
11.34 Mб
Скачать

Количество ядер

урана

 

в

единице объема определяется в

соответствии с заданной величиной “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) будет

 

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

 

 

 

твэл

.

т

вх

П

Ст

т

Ст

вх

Ст

т

Ст

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

– Определение температуры оболочки ТВЭЛ со стороны теплоносителя.

Температура стенки Т', со стороны теплоносителя определяется на

основании решения уравнения теплоотдачи.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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