Расчет характеристик ядерного реактора и эволюции радионуклидов за время кампании и время выдержки облученного ядерного топлива

  • Вид работы:
    Курсовая работа (т)
  • Предмет:
    Физика
  • Язык:
    Русский
    ,
    Формат файла:
    MS Word
    384,24 Кб
  • Опубликовано:
    2013-03-11
Вы можете узнать стоимость помощи в написании студенческой работы.
Помощь в написании работы, которую точно примут!

Расчет характеристик ядерного реактора и эволюции радионуклидов за время кампании и время выдержки облученного ядерного топлива















Расчет характеристик ядерного реактора

1. Исходные данные

Исходные данные из [1] приведены в таблице 1.

Таблица 1. Исходные данные.

Наименование

Принятое обозначение

Численное значение

Мощность, МВт Тепловая: Электрическая:

 W Q

 3070 930

Обогащение, %

q

3,95

Загрузка, т U-мет

m(U-мет)

67,5

Выгорание, МВт-сут / т U-мет

УВ

4,2∙104

Время выдержки ТВЭЛов, сут.

Тохл

1020

Массовое число цепочки

А

107

Коэффициент воспроизводства

Кв

0,752

Заряд ядра нуклидов семейства изобар

Z

35


Тип реактора - ВВЭР-1000 (на тепловых нейтронах).

Химическая формула топлива - UО2.

2. Расчетная часть

.1 Определение удельного выгорания

Поскольку удельное выгорание в МВт·сут/ т U-мет. уже задано, пересчет в другие размерности труда не представляет. Единственное, что надо учесть, это измененный, по сравнению с природным, изотопный состав исходного материала. Поэтому требуется уточнение значений атомной и молярной масс урана и топлива с применением изотопных масс (в углеродной шкале) и заданной степени обогащения топлива U - 235.

Таблица 2.Содержание изотопов урана в природном и обогащенном топливе.

Изотоп

Молярная масса,г/моль*

Природный уран,масс %*

Уран, обогащенный до 3,95 масс % по U-235, масс %

U-238

238,0508

99,2739

96,01875

U-235

235,0439

0,7204

3,95000

U-234

234,0410

0,0057

0,03125

* - [1]

Массовый процент обогащенного урана был найден следующим образом:

В процессе обогащения содержание U-235 увеличилось от 0,7204 масс % до 3,95 масс%, то есть в 5,5 раз. Содержание U-234, который сопровождает U-235 в технологических схемах, также увеличилось в 5,5 раз по сравнению с природным.

%U-234(в топливе) = [q /% U-235(природный)]·U-234(природный) = [3,95/0,7204]·0,0057 = =0,03125 масс%.

Содержание U-238 в топливе было найдено как разность:

%U-238(в топливе) = 100% - %U-234(в топливе) - %U-235(в топливе) =100-0,03125-3,95= = 96,01875 масс%.

Находим молярную массу U-мет по формуле:

М(U-мет) = М(U-238)·%U-238(в топливе)+ М(U-235)·%U-235(в топливе)+ М(U-234)·%U-234(в топливе);

где М(U-238), М(U-235), М(U-234) - молярные массы изотопов U [4]

М(U-мет)=[238,0508·0,9601875]+[235,0439·0,0395]+[234,0410·0,0003125]=

=237,9308 г/моль.

Молярная масса топлива равна:

М(UО2) = 269,9206 г/моль,

где М(О) = 15,9949 г/моль.

Масса топлива в пересчете на UО2:

m(UО2) = [m(U-мет)(U-мет)]·М(UО2),

m(UО2) =[67,5/237,9308] ·269,9206= 76,575 т.

Масса U-235 в загрузке топлива с учетом обогащения:

m(U-235) = m(U-мет)·q = 67,5·0,0395 = 2,6663 т.

2.1.1 Определение удельного выгорания в МВт·сут/т топлива

Удельное выгорание (УВ) топлива меньше удельного выгорания U-мет во столько раз, во сколько масса топлива UО2 больше массы загрузки U-мет:

УВ(U-235) = [m(U-мет)/ m(UО2)]·УВ(U-мет) = [67,5/ 76,575]·4,2·104 =

= 3,7·104 МВт·сут/т топлива.

2.1.2 Определение удельного выгорания в МВт·сут/т U-235

УВ(U-235) = УВ(U-мет)/q = 3,7·104 /0,0395 = 9,37·105 МВт·сут/т U-235.

2.1.3 Определение удельного выгорания в кг U-235 сгоревшего/ кг U-235 загруженного

Можно рассчитать как отношение удельного выгорания U-235 к предельному выгоранию:

УВ = УВ(U-235)/ПВ(U-235) = 9,37·105 /9,4·105 =

=0,9968 кг U-235 сгоревшего/ кг U-235 загруженного;

где ПВ=9,4·105 МВт·сут/т U-235 - предельная величина выгорания.

.1.4 Определение степени выгорания

Степень выгорания - это отношение реального удельного выгорания U-235 к предельному выгоранию:

СВ = 9,37·105 ·100/9,4·105 =99,68%.

.1.5 Определение fifa

Fifa (fissions per initial fussionablе atom) - число актов деления ядер на один первоначально загруженный делящийся атом.

Первоначально загружено ядер U-235:

N(U-235) = [m(U-235)·NA] /M(U-235) = [2,6663·106·6,02·1023] /235,0439 = 6,829·1027 ядер,

где NA=6,02·1023 моль -1 -постоянная Авогадро;

В реакторе при заданном удельном выгорании U-мет делится ядер:

,2·104·67,5·86400·3·1016 = 7,35·1027 ядер,

где 86400 с = 1 сут, 3·1016 дел/с = 1МВт

fifa = 7,35·1027 /6,829·1027 = 1,076

.1.6   Определение fima

Fima (fissions per initial metal atom) - число актов деления ядер на один первоначально загруженный атом тяжелого металла.

Fima меньше fifa в n раз, где:

n = (N35 + N38 )/ N35 >> 1= (0,9601875+0,0395) / 0,0395 = 25,31 раз= fifa / n = 1,076 / 25,31 = 0,0425.

2.1.7 Определение глубины выгорания

Глубина выгорания - доля выгоревших атомов топлива или тяжелого металла.

ГВ = fifa (1 + sЗ/ sД)35.

где sЗ = 97,4 - сечение захвата U-235;

sД = 583,5 - сечение деления U-235 [1]

ГВ = 1,076·(1+97,4 / 583,5) = 0,92208 или 92,21%.

.2 Определение скорости выгорания

Скорость выгорания - это удельная мощность, т.е. мощность, отнесенная к единице массы топлива, металла или U-235.

2.2.1 Определение скорости выгорания в МВт / т топлива

V = W / m(UО2) = 3070/ 76,575= 40,09 МВт / т топлива.

.2.2 Определение скорости выгорания в МВт / т U-мет

V = W / m(U-мет) = 3070 / 67,5 = 45,48 МВт / т U-мет.

.2.3 Определение скорости выгорания в МВт / т U-235

V = W / m(U-235) = 3070 / 2.6663 = 1151,408 МВт / т U-235.

.3 Определение времени кампании

УВU-мет = (W·tK) / mU-мет ;

Следовательно, tK = УВ U-мет ·m U-мет / W,

tK = 4.2·104 ·67,5 / 3070 = 923 сут.

.4 Определение суммарной b- активности продуктов деления с учетом и без учета короткоживущих радионуклидов

Для реакторов на тепловых нейтронах и уранового топлива (любой степени обогащения) набор продуктов деления и их независимые выходы являются постоянными. Это дает возможность табулировать значения суммарной b-активности продуктов деления как функцию времени кампании при неизменной (как правило - единичной) удельной мощности реактора для неизменной (тоже, как правило - единичной) массы топлива любой степени обогащения.

Таблицы для определения b-активности часто имеют достаточно большой шаг по времени кампании и полученное время кампании может оказаться лежащим между двумя табулированными величинами. В этом случае необходимая величина активности может быть найдена путем графической линейной интерполяции. Если время кампании выходит за пределы табулированных величин, то возможно два варианта. Первый - поискать справочник, в котором будет интересующее значение времени кампании. Второй - использовать графическую нелинейную экстраполяцию, имея при этом в виду, что чем больше интервал времени, для которого проводится экстраполяция, тем больше шансов получить заметную погрешность в конечном результате. Далее полученная величина активности умножается на удельную мощность, на массу топлива и окончательный результат выдается в кюри и в беккерелях[1].

2.4.1 Определение суммарной b- активности продуктов деления без учета короткоживущих радионуклидов

2.4.1.1        К концу облучения (Тохл=0)

Справочные данные для заданного реактора из [2] приведены в таблице 3.

Таблица 3. Полная активность смеси продуктов деления 235U, Ки / кг в зависимости от времени кампании tк и времени выдержки Tохл для единичной мощности реактора.

Тохл, сут

tK, сут


130

180

200

360

720

0

673

707

717

763

859

920

720

7,47

9,68

10,6

17,0

27,4

165

1080

4,35

5,68

6,23

10,3

17,6

148


Рисунок 1. Зависимость активности продуктов деления от времени кампании.

По графику находим, что для tK = 923 сут. активность для единичной мощности реактора равняется 880,0 (Ku/кг) / (Вт/г).

Тогда удельная активность для полной мощности реактора:

AWуд = A0уд·V,

где V - удельная мощность (скорость выгорания), МВт / т U-мет.

AWуд = 880 (Ku/кг) / (Вт/г)*45,48 кВт / кг U-мет = 4,0*104 Ku / кг

Полная активность на всю загрузку, Ku:

A = AWуд·m(U-мет)кг = 4,0*104 ·67,5·1000= 2,70·109 Ku.

Полная активность на всю загрузку, Бк:

AS = ASKu) ·3,7·1010 = 2,70·109 ·3,7·1010 ≈ 10,0·1019 Бк,

где 1Ku = 3,7·1010 Бк.

Итоговые результаты расчетов b- активности продуктов деления для Тохл = 0 без учета короткоживущих радионуклидов представлены в таблице 4.

Таблица 4. Удельная и полная b- активность ПД на конец кампании без учета короткоживущих радионуклидов.

Процесс

A0уд , (Ku/кг) / (Вт/г)

АWуд, Ku / кг

А, Ku

А, Бк

Работа реактора - облучение ТВЭЛов

880,0

4,0*104

2,70·109

10,0·1019


2.4.1.2    К концу выдержки (Тохл = 1020 суток)

По данным таблицы 3 строим зависимость активности ПД от времени охлаждения и находим искомую величину активности для единичной мощности реактора при tK = 923 сут, Тохл = 1020 дней.

Рисунок 2. Зависимость активности продуктов деления от времени охлаждения.

Из уравнения у=-0,0272х+47 следует, что при Тохл = 1020 суток, tK =720 дней, активность равна 19,233 (Ku/кг) / (Вт/г).

А для tK = ∞ активность равняется -0,0472*1020+199=150,836 (Ku/кг) / (Вт/г).

По полученным данным можно построить кривую, которая будет характеризовать изменение активности ПД во время выдержки.

Рисунок 3. Зависимость b- активности долгоживущих продуктов деления от времени компании.

По графику находим, что для tK = 923 сут., величина активности ПД для единичной мощности реактора на конец времени выдержки будет равняться 120(Ku/кг) / (Вт/г).

Результаты расчетов в - активности продуктов деления без учета короткоживущих радионуклидов при Тохл = 1020 сут приведены в таблице 5.

Таблица 5.Удельная и полная в- активность на конец выдержки.

Процесс

А°уд, (Ku/кг) / (Вт/г)

АWуд,Ku / кг

А∑,Ku

А∑,Бк

Выдержка ТВЭЛов

120

5,46·103

3,684·108

13,63·1018


2.4.2 Определение суммарной в- активности продуктов деления с учетом короткоживущих радионуклидов

.4.2.1 К концу облучения (Тохл = 0).

Справочные данные для заданного реактора из [2] приведены в таблице 6.

Таблица 6.Удельная в- активность (в кюри) смеси продуктов деления U-235 на 1 кг металла в зависимости от времени кампании tK и времени выдержки Тохл при мощности реактора W = 1 Вт/ U-мет.

Тохл, сут.

tK, сут.


720

0

2307

2368

720

27,4

165

1080

17,6

148


По данным таблицы 6, аналогично пункту 2.4.1.1. находим искомую величину активности для единичной мощности реактора, tK= 923 сут, Тохл= 0 сут (рис.4):

Рисунок 4. Зависимость в- активности смеси продуктов деления от времени кампании.

По данным таблицы 6, аналогично пункту 2.4.1.1. находим искомую величину активности для единичной мощности реактора, tK= 923 дня, Тохл= 0:

А°уд = 2348 Ku·г/кг·Вт.

Результаты расчетов b- активности продуктов деления при Тохл = 0 представлены в таблице 7.

Таблица 7. Удельная и полная в- активности продуктов деления на конец кампании.

Процесс

А°уд, (Ku/кг) / (Вт/г)

АWуд, Ku / кг

А∑, Ku

А∑, Бк

Работа реактора - облучение ТВЭЛов

2,35·103

1,07·105

7,22·109

2,67·1020


2.4.2.2        К концу выдержки (Тохл = 1020 сут).

Значения b- активности ТВЭЛов на конец выдержки с учетом короткоживущих радионуклидов совпадают с соответствующими значениями без их учета, т.к. за время охлаждения Тохл = 1020 сут. короткоживущие радионуклиды распадаются и не вносят вклад в значение активности.

.5 Определение суммарного г - эквивалента продуктов деления

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

2.5.1  Определение суммарного г - эквивалента продуктов деления без учета короткоживущих радионуклидов

2.5.1.1        К концу облучения (Тохл = 0)

Справочные данные для заданного реактора из [2] приведены в таблице

Таблица 8.Удельный г-эквивалент долгоживущих продуктов деления 235U (в г-экв Ra) на 1 кг металла в зависимости от времени кампании tK и времени выдержки Тохл при мощности реактора W = 1 Вт/г 235U.

Тохл, сут.

tK, сут.


720

0

240

251

720

1,38

23

1080

1,08

22,2


По данным таблицы 8, аналогично пункту 2.4.1.1. находим искомую величину активности для единичной мощности реактора, tK= 923 дня, Тохл= 0 (рис.5): А°уд = 242 г-экв Ra / кг.

Рис. 5 Зависимость г - эквивалента долгоживущих продуктов деления от времени кампании.

Результаты расчета г - эквивалента продуктов деления без учета короткоживущих радионуклидов для Тохл = 0 сут приведены в таблице 9.

Таблица 9. Удельный и полный г- эквивалент на конец кампании.

Процесс

А°уд,г-экв Ra / кг

АWуд, г-экв Ra / кг

А∑, Ku

Работа реактора - облучение ТВЭЛов

242

1,10·104

7,43·108


К концу выдержки (Тохл = 1020 сут)

По данным таблицы 8, аналогично пункту 2.4.1.2., находим искомую величину активности для единичной мощности реактора, tK = 923 сут,Тохл= 1020 сут (рис.6 и 7)

Рис.6. Зависимость г-эквивалента долгоживущих продуктов деления от времени выдержки.

Из уравнения у=-0,02722х+47 следует, что при Тохл = 1020 суток, tK =720 дней, активность равна 19,233 (г-экв Ra/кг).

А для tK = ∞ активность равняется -0,04722*1020+199=150,836 (г-экв Ra/кг).

Рис.7. Зависимость г-эквивалента долгоживущих продуктов деления от времени компании.

Результаты расчета г- эквивалента продуктов деления без учета короткоживущих радионуклидов для Тохл = 1020 сут приведены в таблице 10.

Таблица 10. Удельный и полный г- эквивалент на конец выдержки.

Процесс

А°уд, г -экв Ra /кг

АWуд, г -экв Ra/ кг·

А∑,г -экв Ra

Выдержка ТВЭЛов

38,00

1,73·103

1,17·105


2.5.1  Определение суммарного г - эквивалента продуктов деления с учетом короткоживущих радионуклидов

2.5.1.1        К концу облучения (Тохл = 0).

Справочные данные для расчета заданного реактора из [2] приведены в таблице 11.

Таблица 11. Удельный г -эквивалент смеси продуктов деления 235U (в г-экв Ra) на 1 кг металла в зависимости от времени кампании tK и времени выдержки Тохл при мощности реактора W = 1 Вт/г 235U с учетом короткоживущих изотопов.

Тохл, сут.

tK, сут.


720

0

632

643

720

1,38

23,0

1080

1,08

22,2


По данным таблицы 11, аналогично пункту 2.4.1.1. находим искомую величину активности для единичной мощности реактора, tK= 923 сут, Тохл= 0 сут (рис.8):

Рис. 8 Зависимость г-эквивалент смеси продуктов деления от времени кампании.

А°уд = 635 (г-экв Ra/кг)/(Вт/г).

Результаты расчета г - эквивалента продуктов деления с учетом короткоживущих радионуклидов для Тохл = 0 приведены в таблице 12.

Таблица 12. Удельный и полный г- эквивалент на конец кампании.

Процесс

А°уд,г -экв Ra / кг

АWуд,г -экв Ra/ кг·

А∑,г -экв Ra

Работа реактора - облучение ТВЭЛов

635

2,89·104

19,5·108


2.4.2.2    К концу выдержки (Тохл = 1200 сут).

Значения г-эквивалента ТВЭЛов на конец выдержки с учетом короткоживущих радионуклидов совпадают с соответствующими значениями без их учета, т.к. за время охлаждения Тохл = 1020 сут короткоживущие радионуклиды распадаются и не вносят вклад в значение активности.

2.5   Определение плотности потока нейтронов и количества образовавшихся Pu-239, 240, 241 с учетом всех путей их расходования

Расход плутония идёт в трех направлениях: радиоактивный распад, деление, захват. Каждое из направлений характеризуется своей константой, а следовательно характеризуется количественно. Общее количество наработанного плутония за вычетом количества, дожившего до конца кампании, дает общее число израсходованного плутония, которое в свою очередь распределяется пропорционально сечениям деления и захвата между реакциями деления и захвата. Плутоний, претерпевший деление, вносит свой вклад в создании мощности реактора, который не учитывался в расчете пункта 2.1, поэтому реальное количество выгоревшего U-235 будет ниже заявленного - его необходимо пересчитать.

Плутоний претерпевший захват нейтрона, вносит свой вклад в поставку более тяжелого изотопа плутония.

2.5.1  Определение плотности потока нейтронов

рср = 2,2·1010 ·W / m(U-235) [1]

рср - средняя плотность потока нейтронов.

рср = 2,2·1010 ·3070/ 2,6663 = 2,533·1013 см-2·с-1.

.5.2   Определение количества образовавшихся Pu-239, 240, 241 с учетом всех путей их расходования

Цель этого пункта - дать оценку вклада образующегося плутония-239 в работу реактора и определить количество плутония-239 и его более тяжелых изотопов, накопившихся к концу кампании. Это позволит исправить степень выгорания U-235, полученную в п.2.1 по чисто формальному признаку (через энерговыделение), и получить более объективную оценку глубины выгорания.[1]

Рассмотрим цепочку превращений Pu в реакторе:


Эволюция Pu-239 в реакторе на тепловых нейтронах зависит от мощности W, потока нейтронов р, времени кампании tK.

Пусть, W = const, то есть скорость поставки Pu-239 Р239 = const. При этом, р = const, если загрузка U-235 постоянна.

Данные для расчета эволюции Pu представлены в таблице 13.

Таблица 13.Данные для расчета эволюции плутония-239.

Радионуклид

Т1/2,лет

l, с-1

sд, барн

sз, барн

L*з-1

L*д-1

L,с-1

Т°,с

1

2

3

4

5

6

7

8

Pu-239

2,4·104

9,16·10-13

 744

267,2

6,49·10-9

1,81·10-8

2,46·10-8

2,82·107

Pu-240

6,54·103

3,36·10-12

0,03

290

7,04·10-9

7,29·10-13

7,05·10-9

9,83·107

Pu-241

1,45·101

1,52·10-9

1009

368

8,94·10-9

2,45·10-8

3,50·10-8

1,98·107


, 5 - 1 барн = 10-24 см; [1]

6, 7 - L* = s·pср;

8 - L= l + L*з +L*д ;

- Т°= ln2 / L.

Скорость деления ядер:

F = 3,1·1016·W = 3,1·1016·3070 = 9,21·1019 дел/с

Скорость поставки Pu:

P239 = F·KB·(1+sз/sд) = 9,21·1019·0,752·(1+ 97,4/583,5) = 8,082·1019 дел/с.

Общее количество наработанного за время кампании Pu-239:

N239 = P239 ·tK = 8,082·1019·923·86400 = 6,445·1027 ядер

m (Pu-239) = N239 ·M(Pu-239)/NA = 6,445·1027·239,0522/6,02·1023 = 255,94*104г =2559,4 кг ≈ 2,6 т.

где М (Pu-239) = 239,0522 г/моль - молярная масса Pu-239[1].

2.5.2.1        Определение количество образовавшегося Pu-239

Количество Pu-239, дожившего до конца кампании, определим по формуле, тождественной случаю накопления дочернего радионуклида при вековом равновесии.

Вместо постоянной активности материнского нуклида используем постоянную скорость поставки Pu-239 при работе реактора с постоянной мощностью:

2,823*1027 шт.

1121 кг.

 

Общее количество израсходованного Pu-239:


m (Pu-239) - m (Pu-239,tK) = 2559,4 - 1121 = 1438,4 кг.

Рu-39 в реакторе расходуется по трем путям - реакция радиационного захвата, реакция деления и собственный радиоактивный распад. Последний путь расходования при работающем реакторе не вносит значительной вклад, т.к. период полураспада Pu-39 (T1/2 = 24 тыс. лет) значительно превосходит время кампании реактора (tк = 923 сут.), следовательно, определим количества Pu-39, израсходованные только на деление и захват, разделив полученное количество израсходованного Рu-39 через отношение сечений деления и захвата:


т.о. израсходовано на деление Pu-239 - 1057,87 кг

израсходовано на захват Pu-239(на образование Pu-240) - 380,53 кг.

2.5.2.2 Определение количество образовавшегося Pu-240

Количество Pu-240, дожившего до конца кампании:


N240,t =7,5899*1026 шт.

= 302,65 кг.

где М (Pu-240) = 240,0538 г/моль - молярная масса Pu-240.

Общее количество израсходованного Pu-240:

m (Pu-240) - m (Pu-240,tK) = 380,53 - 302,65 = 77,88 кг.


т.о. израсходовано на деление Pu-240 - 8·10-3 кг

израсходовано на захват Pu-240 - 77,872 кг;

то есть, почти весь Pu-240 расходуется на образование Pu-241.

Из полученных значений видно, что практически весь Рu - 40 расходуется на образование Рu - 41 по реакции захвата, т.к. его сечение захвата значительно превосходит сечение деления.

.5.2.3 Определение количество образовавшегося Pu-241

Количество Pu-241, дожившего до конца кампании:

N241(tK) = 8,603·1025 шт.

= 34,45 кг.

где М (Pu-241) = 241,0569 г/моль - молярная масса Pu-240.

Общее количество израсходованного Pu-241:

m (Pu-241) - m (Pu-241,tK) = 77,872 - 34,45 = 43,43 кг.


т.о. израсходовано на деление Pu-241 - 31,81 кг.

израсходовано на захват Pu-241 - 11,61 кг.

Результаты расчета эволюции изотопов Pu приведены в таблице 14.

Таблица 14. Эволюция изотопов плутония во время кампании

Радионуклид

Наработано за время кампании, кг

Израсходовано на деление, кг

Израсходовано на захват, кг

Дожило до конца кампании, кг

Pu-239

2559,40

1057,87

380,53

1121,00

Pu-240

380,53

8·10-3

77,872

302,65

Pu-241

77,87

31,81

11,61

34,45


2.5.3 Определение количества реально сгоревшего U-235

С учетом полученных данных уточним количество реально сгоревшего U-35.

В реакторе при заданном удельном выгорании U-мет делится ядер:

,2·104·67,5·86400·3,1·1016=7,593·1027 ядер,

где 86400 с = 1 сут, 3,1·1016 = 1МВт,

Тогда массу сгоревшего U-235 равна:

m(U-235) = [N(U-235) · M(U-235)] / NA = [7,593·1027·235,0439] /6,02·1023 = 2964,6 кг,

где NA=6,02·1023 моль -1 -постоянная Авогадро;

M(U-235)=235,0439 г/моль - молярная масса U-235.

Учитывая, что часть Pu-239, расходуемая на деление, идет на замену U-235, получим, что реально сгорает U-235:

=

= 2964,6 - 1057,87 - 0,008 - 31,81 = 1874,912 кг.

 кг U-235 / кг топл.

Тогда степень выгорания U-235:

СВреал. = mU-35,сгор/ mU-35,загр. = 1874,912/2666,3 = 70,32%,

где ПВ=9,4·105МВт·сут/т U-235 - предельная величина выгорания.

Уточним значения fifa, fima, глубины выгорания:

в реакторе делится ядер U-235:

N = 1874,912 ·1000·6,02·1023 / 235,0439 = 4,802·1027 шт.

fifa = 4,802·1027 / 6,829·1027 = 0,703

где 6,829·1027 шт.- число, первоначально загруженных ядер U-235.

fima = 0,703/25,31 = 0,0278

ГВ = 0,703·100%·(1+97,4 / 583,5) = 82,04%.

2.6    Определение весового содержания продуктов деления (ПД)

Примем, что на 1 кг сгоревшего U-235 образуется 1 кг ПД. Тогда масса ПД с учетом того, что на 1 кг загруженного U-235 сгорает 703г U-235 (из ftfa) и что масса загрузки U-235 составляет 2666,3 кг, равна 1874,41 г.

24,48 гПД/кгUО2.

43,92 кг ПД / т U-мет.

.7 Определение коэффициента очистки топлива

В общем случае коэффициенты очистки и времена выдержки облученного топлива зависят от многих факторов, таких как тип реактора, тип топливного цикла, дальнейшее использование топлива, степень выгорания и ряда других. Тем не менее есть общие "моменты": выдержка облученного топлива желательна для полного завершения эволюции цепочки U-238 (n,g, 2b)... или Th-232(n,g, 2b)...;распада U-237 в урановом топливе; возможно более заметного снижения удельной активности облученного материала с целью избежать в последующем заметных радиационно химических осложнений при переработке топлива. Очистка дает возможность либо получить новое ядерное топливо для повторного использования, либо направить его (топливо) на другие цели (оборона, использование в народном хозяйстве). Реактор на тепловых нейтронах, гетерогенный, на природном (или слабообогащенном) топливе, с невысокой плотностью потока нейтронов и умеренным временем кампании (т.е. с малой степенью выгорания). Время выдержки лимитируется главным образом необходимостью исключить в будущем работу с ураном-237. Это ведет к временам выдержки порядка 80-100 дней. Коэффициент очистки определяется из основной посылки: активность топлива после радиохимической переработки не должна превышать активности "первичного" урана (на один a- 2b распада ближайших дочерних, что дает удельную активность по b-2,44·10+4 Бк/г).[1] При активности облученного ядерного топлива (Тохл =1020 сут), равной АWуд=5,46 Ku/г U-мет (п. 2.4.1.2.) коэффициент очистки равен:оч = 5,46·3,7·1010 / 2,44·104 = 8,28·106.

Определение суммарной концентрации ПД после растворения ТВЭЛов при соотношении Т:Ж = 1:10. Ранее найдено, что на 1 кг UО2 образуется 24,48 г ПД (п. 2.7.1.). После растворения ТВЭЛов концентрация ПД составит:

 кг ПД/ т жидк.фазы;

где  т.

Определение изменения активности для временных в абсолютных единицах (Бк) для цепочки с заданным А.

Цепочку с А = 107 выписываем из [1]:


Таблица 15. Характеристики радионуклидов цепочки с А = 107.

Радионуклид

hнез, доли

Т1/2

l, с-1

Мо107

1,07E-03

3,5

c

1,98E-01

Tc107

1,82E-04

29,0

с

2,39E-02

Ru107

1,69E-04

4,2

мин

2,75E-03

Rh107

1,10E-06

21,7

мин

5,32E-04

mPd107

4,17E-10

22,0

c

3,15E-02

Pd107

1,19E-10

6,50E+06

лет

3,38E-15

mAg107

0

44,3

с

1,56E-02


пояснение к таблице:

з - независимый выход, доли;

T1/2 - период полураспада, с;

л - постоянная распада, с-1;

Примечания: ветви, для которых приход или уход в соседние цепочки с большим(меньшим) А составляет меньше 5% и для которых л больше Л*з не учитываются.

1.     
106Ru: Л*з = sз·pср = 2,533·1013 ·0,15·10-24 = 3,8·10-12 с-1, так как для Ru л= =2,75·10-3 больше Л*з, то реакцию захвата для него не учитываем.

2.             107Pd: Л*з = sз·pср = 2,533·1013 ·10,1·10-24 = 2,58·10-10 с-1, так как для Pd л меньше Л*з, то реакцию захвата учитываем и необходимо рассчитать обобщенную константу распада: Л = л + Л*з = 2,58·10-10 с-1, т.к. она нам пригодится для дальнейших расчетов.

3.             10646Pdстаб: рассчитаем Л*з для стабильного изотопа палладия из более легкой цепочки: Л*з = sз·pср = 2,533·1013 ·0,24·10-24 = 6,08·10-12 с-1, она нам также пригодится в дальнейших расчетах.

За время кампании (tK = 923 сут)

Накопление Mo

Родоначальник цепочки - Mo является короткоживущим. Его активность уже через 35 секунд достигнет своей предельной величины, равной скорости его поставки, и после этого меняться не будет.

АMo max = PMo = F·hMo = 9,21·1019 ·1,07·10-3 = 9,86·1016 Бк

lg АMo max = 16,994;

Расчетная формула:


расчетный интервал времени 1-35с

Результаты расчета приведены в таблице 16.

Таблица 16.Эволюция Mo за время кампании.

t,сек

N,шт

A,Бк

lgt

lgA

1

8,94E+16

1,77E+16

0,000

16,248

10

4,29E+17

8,49E+16

1,000

16,929

15

4,72E+17

9,35E+16

1,176

16,971

20

4,88E+17

9,67E+16

1,301

16,985

25

4,94E+17

9,78E+16

1,398

16,991

30

4,94E+17

9,83E+16

1,477

16,992

35

4,94E+17

9,85E+16

1,544

16,993

79747200

4,94E+17

9,85E+16

7,902

16,994


пояснение к таблице:- время от начала кампании, с;- число ядер данного изотопа;

A - активность данного изотопа, А=N*л, Бк;

(примечание: последнее значение времени - время кампании);

Накопление Tc

Так как периоды полураспада молибдена и технеция отличаются не сильно, то накопление технеция рассчитываем по двучленной формуле.

Расчетная формула:


Расчетный интервал времени, учитывая что Т1/2(Tc) = 29 с выбираем от 36с до 290 c. К этому времени скорость поставки технеция сравнится со скоростью его распада и в дальнейшем меняться не будет. Этот интервал разбиваем на 10 точек.

Таблица 17. Изменение активности Tc за время кампании.

t,сек

N,шт

A,Бк

lgt

lgA

36

2,55E+18

6,08E+16

1,556

16,78

56

3,41E+18

8,15E+16

1,748

16,91

76

3,95E+18

9,44E+16

1,881

16,97

96

4,28E+18

1,02E+17

1,982

17,01

126

4,56E+18

1,09E+17

2,100

17,04

146

4,66E+18

1,11E+17

2,164

17,05

186

4,76E+18

1,14E+17

2,270

17,06

200

4,78E+18

1,14E+17

2,301

17,06

246

4,81E+18

1,15E+17

2,391

17,06

290

4,82E+18

1,15E+17

2,462

17,06

79747200

4,82E+18

1,15E+17

7,902

17,06

 

Накопление Ru

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

Расчетная формула:


Расчетный интервал времени - от 300 с до 2500 с.

Результаты расчета приведены в таблице 18.

Таблица 18.Эволюция Ru за время кампании.

t,секN,штA,БкlgtlgA





300

2,44E+19

6,70E+16

2,477

16,83

600

3,74E+19

1,03E+17

2,778

17,01

800

4,17E+19

1,15E+17

2,903

17,06

1000

4,42E+19

1,22E+17

3,000

17,08

1300

4,61E+19

1,27E+17

3,114

17,10

1500

4,67E+19

1,29E+17

3,176

17,11

1800

4,72E+19

1,30E+17

3,255

17,11

2000

4,74E+19

1,30E+17

3,301

17,11

2200

1,31E+17

3,342

17,12

2400

4,75E+19

1,31E+17

3,380

17,12

2500

4,75E+19

1,31E+17

3,398

17,12

79747200

4,76E+19

1,31E+17

7,902

17,12


Накопление Rh

Расчетная формула:


Расчетный интервал времени - от 2550 с до 13020 с (217 мин).

Результаты расчета приведены в таблице 19.

Таблица 19.Эволюция Rh за время кампании.

t,секN,штA,БкlgtlgA





2550

1,68E+20

8,92E+16

3,407

16,95

3000

1,84E+20

9,81E+16

3,477

16,99

3500

1,99E+20

1,06E+17

3,544

17,02

4500

2,18E+20

1,16E+17

3,653

17,07

6000

2,34E+20

1,24E+17

3,778

17,09

8000

2,42E+20

1,29E+17

3,903

17,11

9000

2,44E+20

1,30E+17

3,954

17,11

10000

2,45E+20

1,30E+17

4,000

17,11

11500

2,46E+20

1,31E+17

4,061

17,12

13020

2,46E+20

1,31E+17

4,115

17,12

79747200

2,46E+20

1,31E+17

7,902

17,12

 

Накопление Pdm

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

Расчётная формула:


Расчетный интервал времени - от 1 с до 220 с.

Результаты расчета приведены в таблице 20.

Таблица 20.Эволюция Pdm за время кампании.

t,секN,штA,БкlgtlgA





1

3,78E+10

1,19E+09

0,000

9,08

20

5,70E+11

1,80E+10

1,301

10,25

55

1,00E+12

3,16E+10

1,740

10,50

75

1,10E+12

3,48E+10

1,875

10,54

100

1,17E+12

3,68E+10

2,000

10,57

125

1,20E+12

3,77E+10

2,097

10,58

150

1,21E+12

3,81E+10

2,176

10,58

175

1,21E+12

3,83E+10

2,243

10,58

200

1,22E+12

3,83E+10

2,301

10,58

220

1,22E+12

3,84E+10

2,342

10,58

79747200

1,22E+12

3,84E+10

7,902

10,58


Накопление Pd

Палладий накапливается за счет своей цепочки и за счет прихода из более легкой цепочки в результате реакции радиационного захвата стабильным 106-м палладием. Поэтому накопление палладия рассчитываем по двум цепочкам: по своей и по 106-й. По своей цепочке рассчитаем его по одночленной формуле, просуммировав все независимые выходы предшественников, а так как палладий расходуется по двум направлениям: распад и захват, то используем обобщенную константу распада:


Накопление палладия за счёт прихода из 106-й цепочки рассчитаем следующим образом: мы должны учесть реакцию радиационного захвата 106-м стабильным палладием, а т.к. его накопление определяется в основном эволюцией долгоживущего 106Ru, следовательно и рутений мы должны учесть тоже. Ниобий, молибден и технеций являются короткоживущими, поэтому мы приплюсуем их независимые выходы к рутению 106-му. Родием 106-м мы можем пренебречь в силу его незначительного независимого выхода и малого периода полураспада (29,9с), он получается связующим звеном между 106-м рутением и стабильным 106-м палладием, и практически не оказывает никакого влияния на эволюцию цепочки. Накопление палладия 107-го рассчитаем по трёхчленной формуле:


где

ЛPd-107 - обобщенная константа распада Pd-107.

Л*з,Pd-106 - модифицированная константа радиационного захвата Pd-106.

Расчетная формула:


Расчетный интервал времени от 4ч до 11ч, т.к. Pd долгоживущий, то ограничимся 4-мя точками.

Результаты расчета приведены в таблице 21.

Таблица 21. Накопление Pd за время компании.

t,сек

N,шт

A,Бк

lgt

lgA

13020

1,71E+21

5,78E+06

4,12

6,76

500000

6,55E+22

2,21E+08

5,70

8,35

1000000

1,31E+23

4,43E+08

6,00

8,65

10000000

1,31E+24

4,43E+09

7,00

9,65

79747200

1,06E+25

3,58E+10

7,90

10,55


Накопление Agm

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

За время выдержки (Тохл = 1020 сут).

К моменту выгрузки радионуклиды цепочки А = 107 будут иметь характеристики, представленные в таблице 22.

Таблица 22. Исходные данные для расчета эволюции цепочки с массовым числом А = 107 при выдержке ТВЭЛов.

Радионуклид

, с-1

Начальное число ядер,N0

Мо107

1,98E-01

4,94E+17

Tc107

2,39E-02

4,82E+18

Ru107

2,75E-03

4,76E+19

Rh107

5,32E-04

2,46E+20

mPd107

3,15E-02

1,22E+12

Pd107

3,38E-15

1,06E+25


Эволюцию цепочки во время выдержки рассчитаем по уравнению Бейтмана, в соответствии с которым для цепочки вида:

¾® N2 ¾® N3 ¾® Nj ¾® Ni ,


Если принять, что при t =0 N1 = N01, a N2 = N3 = ...Nj...= Ni = 0, то справедливо

Распад Mo

Mo - короткоживущий радионуклид. Его распад считаем по простой формуле:


Расчетный интервал времени от 1с до 125с.

Результаты расчета приведены в таблице 23.

 

Таблица 23. Распад Mo за время выдержки.

t,сек

N,шт

A,Бк

lgt

lgA

1

4,05E+17

8,02E+16

0

16,90

10

6,82E+16

1,35E+16

1

16,13

20

9,42E+15

1,86E+15

1,30

15,27

40

1,80E+14

3,55E+13

1,60

13,55

60

3,42E+12

6,78E+11

1,78

11,83

80

6,52E+10

1,29E+10

1,90

10,11

100

1,24E+09

2,46E+08

2,00

8,39

125

8,81E+06

1,74E+06

2,10

6,24


Распад Tc

Т.к. периоды полураспада Мо и Тс соизмеримы,то распад Тс считаем по двучленной формуле:


Расчетный интервал времени от 1 с до 1000 с.

Результаты расчета приведены в таблице 24.

Таблица 24. Распад Тс за время выдержки.

t,сек

N,шт

A,Бк

lgt

lgA

1

4,79E+18

1,14E+17

0

17,057

50

1,63E+18

3,90E+16

1,70

16,591

75

8,96E+17

2,14E+16

1,88

16,330

100

4,93E+17

1,18E+16

2,00

16,072

150

1,49E+17

3,56E+15

2,18

15,551

200

4,52E+16

1,08E+15

2,30

15,033

500

3,48E+13

8,31E+11

2,70

11,920

1000

2,25E+08

5,37E+06

3,00

6,730


Распад Ru

Т.к. периоды полураспада Тс и Ru соизмеримы, то распад Ru считаем по двучленной формуле:


Расчетный интервал времени от 1 с до 2520 с.

Результаты расчета приведены в таблице 25.

Таблица 25. Распад Ru за время выдержки.

t,сек

N,шт

A,Бк

lgt

lgA

1

4,76E+19

1,309E+17

0

17,117

300

2,35E+19

6,463E+16

2,48

16,810

600

1,03E+19

2,833E+16

2,78

16,452

900

4,51E+18

1,240E+16

2,95

1200

1,98E+18

5,445E+15

3,08

15,736

1500

8,66E+17

2,382E+15

3,18

15,377

1800

3,80E+17

1,045E+15

3,26

15,019

2100

1,66E+17

4,565E+14

3,32

14,659

2500

5,54E+16

1,524E+14

3,40

14,183

3000

1,40E+16

3,850E+13

3,48

13,585

9000

9,56E+08

2,629E+06

3,95

6,420


Распад Rh

Расчетная формула:


Расчетный интервал времени от 1 с до 45000 с.

Результаты расчета приведены в таблице 26.

Таблица 26. Распад Rh за время выдержки.

t,сек

N,шт

A,Бк

lgt

lgA

1

2,46E+20

1,31E+17

0

17,12

2000

1,07E+20

5,71E+16

3,30

16,76

3000

6,31E+19

3,36E+16

3,48

16,53

5000

2,18E+19

1,16E+16

3,70

16,06

7000

7,52E+18

4,00E+15

3,85

15,60

10000

1,52E+18

8,11E+14

4,00

14,91

11000

8,96E+17

4,76E+14

4,04

14,68

12000

5,26E+17

2,80E+14

4,08

14,45

13020

3,06E+17

1,63E+14

4,12

14,21

45000

1,25E+10

6,65E+06

4,65

6,82

 

Распад Pdm

Расчетная формула:


Расчетный интервал времени от 1 с до 220 с.

Результаты расчета приведены в таблице 27.

Таблица 27.Эволюция Pdm за время выдержки.

t,сек

N,шт

A,Бк

lgt

lgA

1

1,18E+12

3,72E+10

0

10,571

40

3,46E+11

1,09E+10

3,30

10,037

80

9,82E+10

3,09E+09

3,48

9,490

120

2,78E+10

8,77E+08

3,70

8,943

160

7,90E+09

2,49E+08

3,85

8,396

200

2,24E+09

7,06E+07

4,00

7,849

240

6,35E+08

2,00E+07

4,04

7,301

280

1,80E+08

5,68E+06

4,08

6,754

360

1,45E+07

4,57E+05

4,12

5,660


Распад Pd

Расчетная формула:

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


За небольшое по сравнению с периодом полураспада палладия время выдержки его активность не изменится и для построения графика можно ограничиться 4-мя точками.

Результаты расчета приведены в таблице 28.

Таблица 28.Эволюция Pd в целом за время выдержки.

t,сек

N,шт

A,Бк

lgt

lgA

500000

1,03E+25

3,47E+10

5,70

10,541

1000000

1,03E+25

3,47E+10

6,00

10,541

10000000

1,03E+25

3,47E+10

7,00

10,541

88128000

1,03E+25

3,47E+10

7,95

10,541


Графические зависимости lg A-lg tK, lg A-lg Tохл для радионуклидов цепочки А= 107 по данным таблиц 16 - 21, 23-28 представлены на рисунке 9.

Рисунок 9. Графические зависимости lg A-lg tK, lg A-lg Tохл для радионуклидов цепочки А= 107.

3. Определение содержания стабильного (конечного) члена цепочки

к концу компании:

N=F*знез*tk=9,21·1019*0,00142*923*86400=10,43*1024 шт,

где знез-общий независимый выход цепочки.

F-скорость деления чдер.

к концу выдержки:

Nост= ∑Nчл.цеп.

В этой формуле не будем учитывать те радионуклиды цепочки, у которых время уменьшения их активности в 10 раз относительно первоначальной не сравнимо с временем выдержки - 1020 сут, значит нет смысла для них рассчитывать числа оставшихся ядер на момент окончания выдержки, т.е. для Mo, Tc, Ru, Rh, Pdm, Agm. Для оставшегося Pd107 с Т1/2=6,5·106 лет рассчитаем число оставшихся ядер цепочки по формуле Бейтмана за время охлаждения.

Число оставшихся ядер для Pd107 рассчитываем по формуле:

NостPd= =1,03E+25шт.

Nстаб= N- Nост = 10,43*1024 - 1,03E+25=0,013*1025 шт- число оставшихся стабильных ядер цепочки на момент окончания выдержки.

mстаб=[Nстаб · MAg] / NA=0,013*1025*107,8682/6,02*1023=23,294 г- масса оставшихся стабильных ядер цепочки на момент окончания выдержки.

3.1   
Эволюция радионуклидов с предложенным Z = 35 в координатах “lg A-lg Tохл

Из [3] выписываем изменение активности радионуклидов элемента с Z = 35 (Br) во время охлаждения для реактора ВВЭР-1000 с обогащением 3,95%.

В [3] приведен только 1 радионуклид Br, у которого Т1/2=35,30 ч.

Для него из [3] выписываем значения удельной активности в зависимости от времени. Справочные данные приведены в таблице 29.

Таблица29.Изменение активности, расп/(с·т-мет) ПД во время охлаждения.

радионуклид

Т1/2

60 сут

90 сут

120 сут

180 сут

1 г

2,3 г

3 г

10 лет

Br-82

35,30 ч

8,08Е+1

-

-

-

-

-

-

-



Для уточнения данных воспользуемся источником [1], из которого выпишем все изотопы брома.

По периодам полураспада данных радионуклидов можно сказать, что время снижения их активности в 10 раз относительно первоначальной не сравнимо с временем охлаждения - 1020 сут., поэтому нет смысла строить для них графическую зависимость, т.к. они полностью распадутся.

Самый долгоживущий из них радионуклид с периодом полураспада равным 35,30 часа уже ко времени равном 90 суток распадется полностью. (Мы можем говорить об этом, опираясь на данные, выписанные из [3])

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

3.2 Определение скорости образования гремучего газа при растворении ТВЭЛов в л/мин·кг топл. и общего количества выделившегося газа за 6 часов

Время растворения ТВЭЛов:

сек;

Суммарная в -активность продуктов деления на конец времени охлаждения (п. 2.4.2.):

 Бк;

Начальный выход водорода (с.41[1]):

шт/100 эв;

Объем водорода за 6 часов растворения (с.41[1]):

л;

где: 2*106 - средняя энергия в -частиц, эВ;

Поскольку гремучий газ - это смесь Н2 и О2 в соотношении 2:1, то:

л;

Скорость образования гремучего газа (с.41[1]):

 л/мин*кг;

топливо ядерный реактор изотоп

Список используемой литературы

1.  Особенности эволюции радионуклидов в природных и техногенных системах/ Учебное пособие / Н.Д. Бетенеков, Ю.В. Егоров, Т.А. Недобух, В.Д. Пузако. Екатеринбург: ГОУ УГТУ-УПИ, 2001, 72с.

2.      Гусев Н.Г. и др. Гамма-излучение радиоактивных изотопов и продуктов деления.-М:НФМЛ,1958.-208 с.

.        Колобашкин В.М. и др. Радиационные характеристики облученного ядерного топлива: Справочник.-М.:Энергоатомиздат.1983.-384 с.

.        Курс лекций по курсу Радиохимии III часть, 2006г. Н.Д. Бетенеков.

Похожие работы на - Расчет характеристик ядерного реактора и эволюции радионуклидов за время кампании и время выдержки облученного ядерного топлива

 

Не нашли материал для своей работы?
Поможем написать уникальную работу
Без плагиата!