Главная страница

Курсовая рх. Федеральное агентство по образованию гоу впо Уральский государственный технический университет упи


Скачать 390.5 Kb.
НазваниеФедеральное агентство по образованию гоу впо Уральский государственный технический университет упи
АнкорКурсовая рх
Дата19.04.2023
Размер390.5 Kb.
Формат файлаdoc
Имя файла№13.doc
ТипДокументы
#1074620
страница2 из 3
1   2   3

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



m (Pu-239) - m (Pu-239,tK) = 1022-537,3 = 484,7 кг.
Количества плутония-239 на замену урана-235 (деление) и образование плутония-240 (захват) определим, разделив полученное количество израсходованного плутония-239 в отношении, равном отношению сечений деления и захвата:

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

израсходовано на захват Pu-239 (на образование Pu-240) – 128,09 кг
Количество Pu-240, дожившего до конца кампании:




N240,t = 2,776·1026 шт.

= 110,7 кг,
где М (Pu-240) = 240,0538 г/моль – молярная масса Pu-240[4].
Общее количество израсходованного Pu-240:
m (Pu-240) - m (Pu-240,tK) = 128,09 – 110,7 = 17,39 кг.


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

израсходовано на захват Pu-240 (на образование Pu-241) – 17,39 кг;

то есть, почти весь Pu-240 расходуется на образование Pu-241.
Количество Pu-241, дожившего до конца кампании:
N
241(tK) = 2,749·1025 шт.


= 11,01 кг,
где М (Pu-241) = 240,0538 г/моль – молярная масса Pu-241[4].
Общее количество израсходованного Pu-241:
m (Pu-241) - m (Pu-241,tK) = 17,39 – 11,01 = 6,38 кг.

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

израсходовано на захват Pu-241(на образование Pu-242) – 1,706 кг
Результаты расчета эволюции изотопов Pu приведены в таблице 14.
Таблица 14

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

Радионуклид

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

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

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

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

1

2

3

4

5

Pu-239

1022

356,61

128,09

537,3

Pu-240

128,09

1,799·10-3

17,39

110,7

Pu-241

17,39

4,674

1,706

11,01




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


Примем, что на 1 кг сгоревшего U-235 образуется 1 кг продуктов деления. Тогда масса продуктов деления, с учетом того, что на 1 кг загруженного U-235 сгорает 537 г U-235 (из fifa) и что масса загрузки U-235 составляет 1505,2 кг, равна 1219,21 кг.


      1. 1219,21 / 48,1 = 25,35 г ПД / кг UО2.

      2. 1219,21 / 42,4 = 28,75 кг ПД / т U-мет.




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


В общем случае коэффициенты очистки и времена выдержки облученного топлива зависят от многих факторов, таких как тип реактора, тип топливного цикла, дальнейшее использование топлива, степень выгорания и ряда других. Тем не менее есть общие "моменты": выдержка облученного топлива желательна для полного завершения эволюции цепочки U-238 (n,, 2)... или Th-232(n,, 2)...; распада U-237 в урановом топливе; возможно более заметного снижения удельной активности облученного материала с целью избежать в последующем заметных радиационно химических осложнений при переработке топлива. Очистка дает возможность либо получить новое ядерное топливо для повторного использования, либо направить его (топливо) на другие цели (оборона, использование в народном хозяйстве).

Реактор на тепловых нейтронах, гетерогенный, на природном (или слабообогащенном) топливе, с невысокой плотностью потока нейтронов и умеренным временем кампании (т.е. с малой степенью выгорания). Время выдержки лимитируется главным образом необходимостью исключить в будущем работу с ураном-237. Это ведет к временам выдержки порядка 80-100 дней. Коэффициент очистки определяется из основной посылки: активность топлива после радиохимической переработки не должна превышать активности "первичного" урана (на один - 2 распада ближайших дочерних, что дает удельную активность по -2,44·10+4 Бк/г).[1]

При активности облученного ядерного топлива (Тохл =96 сут), равной 6,68 Ku/г U-мет (п. 2.4.1.2.) коэффициент очистки равен:
Kоч = 6,68·3,7·1010 / 2,44·104 = 1,013·107.

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


Цепочку с А = 93 выписываем из [1]:
9336Kr → 9337Rb → 9338Sr → 9339Y → 9340Zr → 93m41Nb → 9341Nbстаб
Таблица 15

Характеристики радионуклидов цепочки с А=93

Р/нуклид

ηнез,доли

Т1/2

λ,с-1

1

2

3

 

4

Kr

5,38E-03

1,289

с

5,377E-01

Rb

3,13E-02

6,12

с

1,133E-01

Sr

2,51E-02

439

с

1,579E-03

Y

1,78E-03

9,58

ч

2,010E-05

Zr

2,90E-06

1,53E+06

лет

1,437E-14

mNb

2,58E-11

13,6

года

1,616E-09


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

  1. Y: Λ*з = з·ср = 2,185·1013·0,078·10-24 = 1,704·10-12 с-1, так как для Y λ больше Λ*з, то реакцию захвата для него не учитываем.

  2. Zr: Λ*з = з·ср = 2,185·1013·2,52·10-24 = 5,506·10-11 с-1, так как для Zr λ меньше Λ*з, то реакцию захвата учитываем и необходимо рассчитать обобщенную константу распада: Λ = λ + Λ*з = 5,507·10-11 с-1, т.к. она нам пригодится для дальнейших расчетов.

  3. Sr: для стронция необходимо поправить значение независимого выхода, учитывая его приход из более тяжелой цепочки за счет β-распада с испусканием запаздывающего нейтрона 9437Rb: ηSrΣ = ηSr + ηRb*0,106 = 2,51·10-2 + 1,55·10-2·0,106 = 2,67·10-2, это значение мы в дальнейших расчетах будем принимать за независимый выход стронция.

  4. 9240Zrстаб: рассчитаем Λ*з для стабильного изотопа циркония из более легкой цепочки: Λ*з = з·ср = 2,185·1013·0,26·10-24 = 5,681·10-12 с-1, она нам также пригодится в дальнейших расчетах.

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


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

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

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

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

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


Расчетный интервал времени, учитывая что Т1/2(Rb) = 6,12 с выбираем от 1с до 1,2 мин. К этому времени скорость поставки рубидия сравнится со скоростью его распада и в дальнейшем меняться не будет. Этот интервал разбиваем на 10 точек.

Результаты расчета приведены в таблице 17.
Накопление Sr

Так как период полураспада стронция больше 10-ти периодов полураспада криптона и рубидия, то примем, что скорость поставки стронция постоянна и состоит из собственного независимого выхода и суммы независимых выходов его предшественников.
Расчетная формула:


Расчетный интервал времени – от 1,5 мин до 75 мин.

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


Накопление Y

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


Расчетный интервал времени – от 80 мин до 96 ч.

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

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



Накопление циркония за счет прихода из 92-й цепочки рассчитаем по двучленной формуле, взяв в качестве независимого выхода стабильного циркония-92 суммарный выход 92-й цепочки, так как все предшественники стабильного циркония короткоживущие:

где ηΣ – суммарный выход 92-й цепочки

ΛZr – обобщенная константа распада Zr-93

Λ*з,Zr-92 – модифицированная константа радиационного захвата Zr-92

ηΣZr-93 – суммарный независимый выход циркония-93(по своей цепочке).

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



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

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

Накопление Nb

Эволюция ниобия целиком зависит от его долгоживущего предшественника циркония, расчет также необходимо вести по двум цепочкам: по двучленной формуле по своей 93-й цепочке и по трехчленной формуле по цепочке 92Zr→93Zr→93mNb. Остальные рассуждения аналогичны рассуждениям при расчете циркония-93-го.

По своей 93-й цепочке:



По второй цепочке:



Пояснения к формулам те же, что и к формулам для циркония.

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



Расчетный интервал тот же, что и для циркония.

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


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


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

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

Р/нуклид

λ,с-1

N0,шт

Kr

5,377E-01

4,487E+17

Rb

1,133E-01

1,453E+19

Sr

1,579E-03

1,800E+21

Y

2,010E-05

1,454E+23

Zr

1,437E-14

1,9286E+26

mNb

1,616E-09

1,43E+20


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

N1  N2  N3  Nj  Ni ,

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

Распад Кr
Кr – короткоживущий радионуклид. Его распад считаем по простейшей формуле:



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

Результаты расчета приведены в таблице 23.
Распад Rb
Расчетная формула:

В виду малого отличия в периодах полураспада распад рубидия считаем по двучленной формуле:



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

Результаты расчета приведены в таблице 24.
Распад Sr
Так как Sr – более долгоживущий, чем его предшественники, то формулу Бейтмана можно упростить, просто просуммировав начальные числа ядер предшественников и стронция.
Расчетная формула:


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

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

Распад Y

Для него рассуждения аналогичны рассуждениям для стронция.

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


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

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

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

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



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

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

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

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



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

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


Таблица 16.Накопление Kr

t, с

N, шт

A, Бк

Lgt

LgA

1

1,866E+17

1,004E+17

0

17,002

1,7

2,688E+17

1,446E+17

0,23

17,160

3

3,593E+17

1,932E+17

0,48

17,286

7

4,383E+17

2,357E+17

0,85

17,372

15

4,486E+17

2,412E+17

1,18

17,382

50

4,487E+17

2,413E+17

1,7

17,383

500

4,487E+17

2,413E+17

2,7

17,383

10000

4,487E+17

2,413E+17

4

17,383

200000

4,487E+17

2,413E+17

5,3

17,383

5000000

4,487E+17

2,413E+17

6,7

17,383


Таблица 17.Накопление Rb

t, с

N, шт

A, Бк

Lgt

LgA

1

1,380E+18

1,563E+17

0

17,194

4

4,996E+18

5,659E+17

0,6

17,753

10

9,665E+18

1,095E+18

1

18,039

25

1,364E+19

1,544E+18

1,4

18,189

100

1,452E+19

1,645E+18

2

18,216

500

1,453E+19

1,645E+18

2,7

18,216

2500

1,453E+19

1,645E+18

3,4

18,216

25000

1,453E+19

1,645E+18

4,4

18,216

300000

1,453E+19

1,645E+18

5,48

18,216

5500000

1,453E+19

1,645E+18

6,74

18,216


Таблица 18.Накопление Sr

t, с

N, шт

A, Бк

Lgt

LgA

90

2,385E+20

3,766E+17

1,95

17,576

150

3,797E+20

5,994E+17

2,18

17,778

250

5,872E+20

9,271E+17

2,4

17,967

450

9,157E+20

1,446E+18

2,65

18,160

1200

1,530E+21

2,415E+18

3,08

18,383

3500

1,793E+21

2,831E+18

3,54

18,452

13000

1,800E+21

2,843E+18

4,11

18,454

70000

1,800E+21

2,843E+18

4,85

18,454

600000

1,800E+21

2,843E+18

5,78

18,454

7000000

1,800E+21

2,843E+18

6,85

18,454



Таблица 19.Накопление Y

t, с

N, шт

A, Бк

Lgt

LgA

4800

1,337E+22

2,688E+17

3,68

17,429

10000

2,647E+22

5,321E+17

4

17,726

20000

4,813E+22

9,673E+17

4,3

17,986

35000

7,345E+22

1,476E+18

4,54

18,169

80000

1,163E+23

2,337E+18

4,9

18,369

200000

1,428E+23

2,870E+18

5,3

18,458

450000

1,454E+23

2,922E+18

5,65

18,466

1500000

1,454E+23

2,922E+18

6,18

18,466

5000000

1,454E+23

2,922E+18

6,7

18,466

20000000

1,454E+23

2,922E+18

7,3

18,466


Таблица 20.Накопление Zr

t, с

N, шт

A, Бк

Lgt

LgA

349200

1,020E+24

1,466E+10

5,54

10,166

3000000

8,767E+24

1,260E+11

6,48

11,100

15000000

4,382E+25

6,296E+11

7,18

11,799

66096000

1,928E+26

2,771E+12

7,82

12,443


Таблица 21.Накопление Nb

t, с

N, шт

A, Бк

Lgt

LgA

349200

3,81E+15

6,16E+06

5,54

6,789

3000000

2,31E+17

3,74E+08

6,48

8,573

15000000

5,99E+18

9,68E+09

7,18

9,986

66096000

1,43E+20

2,31E+11

7,82

11,364


Таблица 23.Распад Kr

t, с

N, шт

A, Бк

Lgt

LgA

1

2,621E+17

1,409E+17

0

17,149

1,5

2,003E+17

1,077E+17

0,176091

17,032

2

1,531E+17

8,231E+16

0,30103

16,915

3

8,940E+16

4,808E+16

0,477121

16,682

4

5,222E+16

2,808E+16

0,60206

16,448

5,5

2,331E+16

1,253E+16

0,740363

16,098

7

1,040E+16

5,595E+15

0,845098

15,748

9

3,549E+15

1,909E+15

0,954243

15,281

11

1,211E+15

6,511E+14

1,041393

14,814

13

4,130E+14

2,221E+14

1,113943

14,347


Таблица 24.Распад Rb

t, с

N, шт

A, Бк

Lgt

LgA

1

1,315E+19

1,489E+18

0

18,173

2

1,184E+19

1,341E+18

0,30103

18,127

4

9,529E+18

1,079E+18

0,60206

18,033

7

6,818E+18

7,722E+17

0,845098

17,888

12

3,876E+18

4,390E+17

1,079181

17,643

20

1,567E+18

1,775E+17

1,30103

17,249

30

5,048E+17

5,718E+16

1,477121

16,757

40

1,627E+17

1,842E+16

1,60206

16,265

50

5,241E+16

5,936E+15

1,69897

15,773

62

1,346E+16

1,525E+15

1,792392

15,183


Таблица 25.Распад Sr

t, с

N, шт

A, Бк

Lgt

LgA

1

1,812E+21

2,862E+18

0

18,457

3

1,807E+21

2,853E+18

0,477121

18,455

8

1,793E+21

2,830E+18

0,90309

18,452

24

1,748E+21

2,760E+18

1,380211

18,441

70

1,625E+21

2,566E+18

1,845098

18,409

190

1,345E+21

2,123E+18

2,278754

18,327

500

8,243E+20

1,302E+18

2,69897

18,114

1300

2,331E+20

3,680E+17

3,113943

17,566

2700

2,556E+19

4,035E+16

3,431364

16,606

4390

1,773E+18

2,799E+15

3,642465

15,447


Таблица 26.Распад Y

t, с

N, шт

A, Бк

Lgt

LgA

1

1,472E+23

2,959E+18

0

18,471

6

1,472E+23

2,959E+18

0,778151

18,471

30

1,471E+23

2,957E+18

1,477121

18,471

130

1,468E+23

2,951E+18

2,113943

18,470

700

1,452E+23

2,918E+18

2,845098

18,465

4000

1,359E+23

2,730E+18

3,60206

18,436

20000

9,849E+22

1,980E+18

4,30103

18,297

70000

3,606E+22

7,247E+17

4,845098

17,860

190000

3,233E+21

6,497E+16

5,278754

16,813

344880

1,438E+20

2,890E+15

5,537668

15,461



Таблица 27.Распад Zr

t, с

N, шт

A, Бк

Lgt

LgA

1

1,929E+26

2,771E+12

0,00

12,443

100000

1,929E+26

2,771E+12

5,00

12,443

1000000

1,929E+26

2,771E+12

6,00

12,443

8294400

1,929E+26

2,771E+12

6,92

12,443


Таблица 28.Распад Nb

t, с

N, шт

A, Бк

Lgt

LgA

1

1,43E+20

2,31E+11

0,00

11,364

10000

1,43E+20

2,31E+11

4,00

11,364

1500000

1,47E+20

2,37E+11

6,18

11,374

3500000

1,51E+20

2,45E+11

6,54

11,388

8294400

1,63E+20

2,63E+11

6,92

11,420
1   2   3


написать администратору сайта