Hindawi Publishing Corporation Science and Technology of Nuclear Installations Volume 2015, Article ID 104739, 23 pages http://dx.doi.org/10.1155/2015/104739
Research Article Transuranics Transmutation Using Neutrons Spectrum from Spallation Reactions Maurício Gilberti,1,2 Claubia Pereira,2,3 Maria Auxiliadora F. Veloso,2,3 and Antonella Lombardi Costa2,3 1
Eletrobr´as Termonuclear S.A.-ELETRONUCLEAR, Departamento GEI.T, BR 101 Sul-Rodovia Governador M´ario Covas Km 517, 23948–000 Angra dos Reis, RJ, Brazil 2 Programa de P´os-graduac¸a˜ o em Ciˆencias e T´ecnicas Nucleares (PCTN), Departamento de Engenharia Nuclear, Escola de Engenharia, Universidade Federal de Minas Gerais, Avenida Antˆonio Carlos 6627, 31270-901 Belo Horizonte, MG, Brazil 3 Instituto Nacional de Ciˆencia e Tecnologia de Reatores Nucleares Inovadores/CNPq, Brazil Correspondence should be addressed to Claubia Pereira;
[email protected] Received 14 October 2014; Revised 20 March 2015; Accepted 20 March 2015 Academic Editor: Alejandro Clausse Copyright © 2015 Maur´ıcio Gilberti et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The aim is to analyse the neutron spectrum influence in a hybrid system ADS-fission inducing transuranics (TRUs) transmutation. A simple model consisting of an Accelerator-Driven Subcritical (ADS) system containing spallation target, moderator or coolant, and spheres of actinides, “fuel,” at different locations in the system was modelled. The simulation was performed using the MCNPX 2.6.0 particles transport code evaluating capture (𝑛, 𝛾) and fission (𝑛, 𝑓) reactions, as well as the burnup of actinides. The goal is to examine the behaviour and influences of the hard neutron spectrum from spallation reactions in the transmutation, without the contribution or interference of multiplier subcritical medium, and compare the results with those obtained from the neutron fission spectrum. The results show that the transmutation efficiency is independent of the spallation target material used, and the neutrons spectrum from spallation does not contribute to increased rates of actinides transmutation even in the vicinity of the target.
1. Introduction In the last 50 years, the nuclear systems have produced amounts of highly radioactive waste. The spent fuel consists of a wide variety of elements and isotopes. Most of the long-lived isotopes come from few transuranic elements as plutonium, neptunium, americium, and curium and some fission products such as technetium and iodine. Every year, a LWR of 1 GWh(e) discharges about 21 tons of radioactive fuel with the following inventory: 20 tons of uranium containing 0.9% (180 kg) of 235 U, 200 kg of plutonium, 21 kg of minor actinides (10 kg of neptunium, 10 kg of americium, and 1 kg of curium), and 760 kg of fission products (18 kg of 99 Tc, 16 kg of 93 Zr, 9 kg of 135 Cs, 5 kg of 107 Pd, and 3 kg of 129 I) that are radioactive isotopes with a long half-life [1]. Transmutation of Pu, minor actinides, and long half-life fission products is
a promising concept to reduce the radioactive waste and its long-term radiotoxicity. Many different technologies have been investigated for transmutation, including a variety of reactors, fusion-fission systems, and ADSs according to the long-term programme on Partitioning and Transmutation P&T of the Nuclear Energy Agency. The programme focuses on the possible benefits for waste management and used specific fuel cycle strategies, covering plutonium recycling and additional burning of minor actinides in dedicated reactor systems [2]. To understand and quantitatively assess the advantages and drawbacks of different plutonium and minor actinides (MA) burning strategies, it is important to investigate the mechanisms involved in the transmutation. The systems or strategies proposed are primarily distinguished according to their fast or thermal neutron spectra and/or used combined
2
Science and Technology of Nuclear Installations 225.00 200.00 Neutron yield
175.00 150.00 125.00 100.00 75.00 50.00 25.00 U-Cm 3%
U-Am 3%
U-Pu 3%
U-Np 3%
U-3%
U
Th
Hg
PbBi
Bi
Pb
W
0.00
Target Leakage Capture
Total
Figure 1: Neutrons produced per incident proton in target, beam energy 2.0 GeV. Image adapted from [13], spallation simulation in MCNPX 2.6.0, Bertini Model. Z
330 cm Y
20 cm
∅ 330 cm
Transmutation cells Target 300 cm
100 cm X
15 cm
45∘
X
120 cm
Moderator/coolant
Proton beam (a)
(b)
Figure 2: Radial (right) and axial (left) view of the proposed model, according to the Cartesian planes 𝑋𝑍 and 𝑋𝑌, respectively. The origin of the Cartesian coordinate system is located in the centre of the array (image obtained from the MCNPX Vised, version X 22S).
Figure 3: Arrangement view in three dimensions showing the external vessel, the target (blue), transmutation cells (black), and tube of beam (pink) (image obtained from the MCNPX Vised, version X 22S).
Science and Technology of Nuclear Installations
3
(U, 0, Am)
0.07
(U, 0, Cm)
0.09 0.08
0.06
0.07 0.05 0.06 0.04
0.05
B
B 0.04
0.03
0.03 0.02 0.02 0.01 0.00
0.01 0
15
30
45
60
75
90
0.00
105 120 135 150
0
15
30
45
Distance (cm)
75
90
105 120 135 150
Distance (cm)
(U, 0, Np)
0.14
60
(U, 0, Pu)
0.18 0.16
0.12
0.14 0.10 0.12 0.08
0.10
B
B 0.08
0.06
0.06 0.04 0.04 0.02 0.00
0.02 0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
0.00
0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
Figure 4: Relative number of reactions 𝐵 in (U, 0, and TRU) as a function of the distance to the TRU. The proton beam energy of 2.0 GeV and 400,000 particles without moderator/coolant.
burning schemes [3–7]. Transuranics transmutation can be planned in both thermal and fast reactors. However, a critical reactor with solid fuel and liquid coolant, thermal and fast, in which neutron production and loss are in balance, may contain a quite limited amount of particular TRU mixture components like Pu, Np, Am, or Cm isotopes, due to safety constraints. This is one of the main reasons why the accelerator-driven systems (ADS) are currently studied worldwide for nuclear waste burning. ADS reactors have been proposed for many applications such as energy production, fertile-to-fissile transmutation, and conversion of long-lived
radioisotopes to stable or much shorter-lived isotopes [8– 11]. The transmutation technology to burnup the radioactive isotopes with long half-lives using an accelerator-driven subcritical system is one of the best solutions. ADSs can be designed to have a fast neutron energy spectrum, and they are subcritical being safer than a fast reactor. Another advantage of ADSs is that they can burn mixtures of material that would not maintain criticality in a reactor [12]. This paper is a contribution to best comprehension of mechanisms involved in the transuranic transmutation in ADS-fission hybrid systems. The goal is to compare
4
Science and Technology of Nuclear Installations (Hg, 0, Am)
0.06
(Hg, 0, Cm)
0.07 0.06
0.05
0.05 0.04 0.04 B 0.03
B 0.03
0.02 0.02 0.01
0.00
0.01
0
15
30
45
0.00
60 75 90 105 120 135 150 Distance (cm) (Hg, 0, Np)
0.10
0
15
30
45
(Hg, 0, Pu)
0.14
0.09
60 75 90 105 120 135 150 Distance (cm)
0.12
0.08 0.10
0.07 0.06
0.08
B 0.05
B 0.06
0.04 0.03
0.04
0.02 0.02
0.01 0.00
0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
0.00
0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
Figure 5: Relative number of reactions 𝐵 in (Hg, 0, and TRU) as a function of the distance to the TRU. The proton beam energy of 2.0 GeV and 400,000 particles without moderator/coolant.
the influence of the neutrons spectrum from spallation reactions without the contribution of neutron spectrum from fission reactions and to compare the results with the standard fission source. Building a model to assess only the effects of neutrons spectrum from spallation in the TRU transmutation, the simulations evaluate the behaviour of the neutron flux emitted by the target (source of neutrons), focusing on the transuranic elements—Pu and minor actinides (MA)— considering graphite as the moderator and lead-bismuth and sodium, as coolant. The simulation involved the use of the Los Alamos High Energy Transport Code (LAHET) into Monte Carlo code MCNPX 2.6. The LAHET allows simulating
the transport of a wide variety of particles and to analyse the production of neutrons through interaction of high energy protons with several targets material, in any form of geometry. The neutrons from spallation reactions are used to drive the subcritical nuclear assembly for calculations of longlived TRU isotopes transmutation. The targets chosen, energy and geometry, were defined based on optimization developed in [13, 14]. These targets are also commonly reported in the literature. The analysis will be focused on the capture (𝑛, 𝛾) and fission (𝑛, 𝑓) reactions in the layers of transmutation allowing verifying the influences of the spallation source submitted
Science and Technology of Nuclear Installations (PbBi, 0, Am)
0.07
5 (PbBi, 0, Cm)
0.08 0.07
0.06
0.06
0.05
0.05 0.04 B
B 0.04 0.03 0.03 0.02
0.02
0.01 0.00
0.01
0
15
30
45
60 75 90 105 120 135 150 Distance (cm) (PbBi, 0, Pu)
0.14 0.12
0.00
0
15
30
45
60 75 90 105 120 135 150 Distance (cm) (PbBi, 0, Np)
0.12
0.10
0.10 0.08 0.08 B
B 0.06 0.06 0.04 0.04 0.02
0.02 0.00
0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
0.00
0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
Figure 6: Relative number of reactions 𝐵 in (PbBi, 0, and TRU) as a function of the distance to the TRU. The proton beam energy of 2.0 GeV and 400,000 particles without moderator/coolant.
in the empty, moderator, or coolant medium. A second extended analysis was generated comparing neutrons spectrum from spallation reactions versus neutrons spectrum from fission reactions. Reactions rates and the burning and depletion of actinides were analysed.
2. Calculation Procedure and Modelling 2.1. Methodology. The methodology was based on the model used by Hashemi-Nezhad and collaborators [15, 16], in which some modifications and new aspects were considered, such as the target and energy beam optimization, according to results
obtained in [13]. Modifications in the geometry of the system allowed a significant reduction of the relative error, for values of the ≈1% in all the cells. Also, it opted for more realistic applications of P&T, based on the evaluation of a full isotopic composition of TRU. The isotopic composition of output by mass of fuel burnt in a power plant of 1000 MWe PWR after 5 years cooling is shown in Table 1, and it is used as a reference for defining the percentage of the isotopic composition used in the transmutation layer. The arrangement is simple and maintains a wide margin of subcriticality in order to avoid computation of the influence of fission neutrons generated in a subcritical medium.
6
Science and Technology of Nuclear Installations 1.6E − 03
2.5E − 03
Neutron flux relative (1/cm2 )
Neutron flux relative (1/cm2 )
3.0E − 03
2.0E − 03 1.5E − 03 1.0E − 03 5.0E − 04
1.4E − 03 1.2E − 03 1.0E − 03 8.0E − 04 6.0E − 04 4.0E − 04 2.0E − 04 0.0E + 00 1.00E − 04 0.02636 0.05909 0.09182 0.34545 0.67273 1 4.2727 7.5455 18.182 50.909 83.636 263.64 590.91 918.18 1272.7 1636.4 2000
1.00E − 03 0.03454 0.06727 0.1 0.42727 0.75455 1.8182 5.0909 8.3636 26.364 59.091 91.818 345.45 672.73 1000 1545.5
0.0E + 00
Energy (MeV)
Energy (MeV)
W
U 2.0E − 03
2.5E − 03
Neutron flux relative (1/cm2 )
2.0E − 03 1.5E − 03 1.0E − 03 5.0E − 04
1.6E − 03 1.4E − 03 1.2E − 03 1.0E − 03 8.0E − 04 6.0E − 04 4.0E − 04 2.0E − 04 0.0E + 00
1.00E − 03 0.03454 0.06727 0.1 0.42727 0.75455 1.8182 5.0909 8.3636 26.364 59.091 91.818 345.45 672.73 1000 1545.5
0.0E + 00
1.8E − 03
1.00E − 04 0.02636 0.05909 0.09182 0.34545 0.67273 1 4.2727 7.5455 18.182 50.909 83.636 263.64 590.91 918.18 1272.7 1636.4 2000
Neutron flux relative (1/cm2 )
3.0E − 03
Energy (MeV) PbBi
Energy (MeV) Hg
Figure 7: Neutrons emission spectrum for spallation reaction in different targets. 2.0 GeV proton beam, Bertini model.
It is necessary to ensure that only the influence of the hard neutron flux of yield target is computed, by the processes of intra/internuclear cascade, evaporation, and (𝑛, 𝑥𝑛) reactions, involved in the spallation reaction. The configuration for the target, beam energy, and transport particles were based on a study presented in [14]. The simulation consists of a parabolic proton beam of 2.0 GeV, focusing on the top of the surface of a cylindrical target. For the MCNPX 2.6.0 simulations, it used Bertini intranuclear cascade model, continuous ENDF/B-VI neutrons microscopic cross sections at 293.6 K for all material, and LANL/T data library for beam proton. As shown in [13], temperature has little significance in the results presented here. In spallation reaction, it used 2×106 incident particles— protons (totalling > 107 neutron histories). It was considered a wide range of transport particles—muons, pions, alpha, deuterons, tritium, and heavy ions that take part in the secondary cascade reactions [13, 17].
2.2. Modelling and Calculations. The system consists of a subcritical arrangement with target immersed in a moderator or coolant with transmutation material. For the analysis of spallation neutron flux, it considers three different target materials—eutectic mixture of PbBi, Hg, and natural U. These materials and others were selected because they are commonly used in nuclear power systems and some of them are proposed as target ADS systems. In [13] their properties as spallation target are widely discussed. Figure 1 and Table 2 show the neutron emissivity results (neutron leakage, total neutrons produced, and neutrons captured, all in the target) for some target materials. The hard neutron flux from spallation is analysed and compared in three conditions: (a) target without moderator or coolant and with transmutation material—named Reference System; (b) target immersed in a coolant with transmutation material— Fast System; and (c) target immersed in a moderator with transmutation material—Thermal System. For compatibility
Science and Technology of Nuclear Installations 2.0E − 02
Table 1: Mass of each isotope of actinides in the fuel burned by a 1000 MWe PWR plant after 5 years cooling (adapted from [22]).
1.8E − 02
Relative neutron flux (1/cm2 )
7
1.6E − 02
Radionuclide 234 U
1.4E − 02 1.2E − 02 1.0E − 02 8.0E − 03 6.0E − 03 4.0E − 03 2.0E − 03 0.00 0.01 0.03 0.04 0.06 0.08 0.09 0.18 0.35 0.51 0.67 0.84 1.00 2.64 4.27 5.91 7.55 9.18 18.18 34.55 50.91 67.27 83.64 100.00 263.64 427.27 590.91 754.55 918.18 1,090.90 1,272.70 1,454.50 1,636.40 1,818.20 2,000.00
0.0E + 00
Uranium target Fission source
Energy (MeV) Evaporation source
Figure 8: Neutrons emission spectrum for spallation reaction in U target (30 cm diameter, 2.0 GeV proton beam, Bertini Model) compared to MCNPX 2.6.0 standard evaporation and fission source.
Activity (Ci) 1.94𝐸01
235
U
2.15𝐸02
4.61𝐸 − 01
236
U
1.14𝐸02
7.22
237
U
9.15𝐸 − 07
7.47𝐸01
238
U
2.57𝐸04
Total
2.60E04
237
Np
2.04𝐸01
239
Np
2.05𝐸 − 06
8.56 𝛼 3.56E01 𝛽 7.47E01 1.44𝐸01 4.78𝐸02 𝛼 1.44E01
Total
2.04E01
236
2.51𝐸04
𝛽 4.78E02 1.34𝐸02
Pu
238
Pu
5.99
1.01𝐸05
239
Pu
1.44𝐸02
8.82𝐸03
240
Pu
5.91𝐸01
1.30𝐸04
Pu
2.77𝐸01
2.81𝐸06
9.65
3.76𝐸01
241
the same nomenclature used in Hashemi-Nezhad et al. [15] was maintained, naming the systems as (target, medium, and TRU). The analysed model consists of an outer vessel of 330 cm diameter and 300 cm in length, filled with moderator/coolant homogeneous, occupying a volume (∼26 m3 ) with the target and spheres embedded. The target cylinder 20 cm in diameter and 1 meter long is centred in the 𝑋𝑌 plane at positions 𝑍 = −30 cm and 𝑍 = 70 cm, with the spheres filled with the material transmutation radially arranged in the 𝑋𝑌 plane and 𝑍 = 0. The spheres (transmutation layer) are each filled with transuranic elements (Np, Pu, Am, and Cm) according to the percentage of the isotope given in Table 1. Proton beam energy of 2.0 GeV is introduced into the system along the 𝑧-axis in a hollow tube of length 120 cm, internal diameter of 6 cm, and thickness of 1.5 cm. The tube material is composed of alloy steel H9. The proton beam falls directly on the window of the tube toward the target crossing axially. This beam has a parabolic distribution with 1.4 cm of diameter [13]. Preliminary simulations have shown that, for obtaining a relative error below 10% calculated in all cells, it is necessary to transport more than 𝑛𝑝𝑠 > 106 protons per simulation, requiring considerable computational time. Due to this long execution time, some alternative geometry is considered such as increasing the number of spheres and diameter, still maintaining a large margin of subcriticality as shown in Figures 2 and 3. The modification was made increasing the diameter of the balls from 1.5 cm to 3.0 cm, increasing the total number of spheres from 20 to 80, constituting 8 rows with 10 spheres spaced in 15 cm from each other and away from the centre, and occupying the initial distance of 15 until 150 cm. It was considered the sum of records (events, traces) of neutrons that reach on all cells (rings) equidistant from the centre of the array. This analysis is more coherent with
Mass (kg) 3.14
242
Pu
Total 241
Am
242m 243
Am
Am
Total
2.46E02 1.32
𝛼 1.23E05 𝛽 2.81E06 4.53𝐸03
1.19𝐸02
1.16𝐸02
2.48
4.77𝐸02
3.81 1.33𝐸 − 01
𝛼 5.01E03 𝛽 1.16E02 4.40𝐸05
242
Cm
243
Cm
1.96𝐸 − 03
9.03𝐸01
244
Cm
9.11𝐸 − 01
7.38𝐸04
245
Cm
5.54𝐸 − 02
9.79
246
Cm
6.23𝐸 − 03
1.92
Total
1.11
Total
2.63𝐸04
𝛼 5.14E05 𝛼 6.42𝐸05 𝛽 2.81𝐸06
the practice, where the composition of actinides follows a reprocessing matrix fuel instead of the rates isolated analysis for each isotope. The results obtained are compared with the typical pattern of a conventional critical system. The results are shown in terms of the rates of reactions (𝑛, 𝛾) and (𝑛, 𝑓) in the cells (spheres) filled with transmutation material. These reactions are the more important involving the TRU transmutations [15, 18, 19]. The results are obtained through a special treatment of MCNPX card event record that allows number counting of captures in a given nuclide or a specific combination of nuclides at the end of each history, converting the pulses created in a power distribution of a detector in neutron capture. The number
8
Science and Technology of Nuclear Installations (U, Na, Am)
(U, Na , Cm)
0.35
0.60
0.30
0.50
0.25 0.40 0.20 B 0.30
B 0.15
0.20
0.10
0.10
0.00
0.05
0
15
30
45
60 75 90 105 120 135 150 Distance (cm) (U, Na, Np)
0.80
0.00
0
15
30
45
60 75 90 105 120 135 150 Distance (cm) (U, Na, Pu)
1.00 0.90
0.70
0.80 0.60
0.70
0.50
0.60
B 0.40
B 0.50 0.40
0.30
0.30 0.20 0.20 0.10 0.00
0.10 0
15
30
45
B-fission B-capture
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
0.00
0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
Figure 9: Relative number of reactions 𝐵 in (U, Na, and TRU) as a function of the distance to the TRU. For proton beam energy of 2.0 GeV and 400,000 particles.
of records of fission events occurring in a given cell is obtained indirectly, by converting the amount of average energy deposited in the cell by fission events. This conversion is made considering the material density, bulk, and 𝑄 value (average energy generated by each fission event ∼200 MeV).
3. Analysis of Results 3.1. Reference System. The system without moderator is used as a reference model for comparison with the fast sodium or fast lead-bismuth systems and thermal graphite system. Figure 4 shows the rate of TRU transmutation via reactions (𝑛, 𝑓) and (𝑛, 𝛾) as a function of the parameter 𝐵 and the distance 𝑥 from the centre of the system (U, 0, and TRU),
Figure 5 for (Hg, 0, and TRU), and Figure 6 for (PbBi, 0, and TRU). The parameter 𝐵 is defined as the relative number of interactions, while 𝐵 will represent the relative number of 𝐵 interactions normalized to 1 g of material. To systems without any moderation and scattering, the transmutation mechanism is more effective since all nuclides have much larger (𝑛, 𝑓) fission reaction rates than (𝑛, 𝛾) capture reaction rates. These behaviours can best be represented by the average value of the ratio between the capture and fission reactions 𝑅 = (𝐵(𝑛,𝑓) /𝐵(𝑛,𝛾) ) as shown in Table 3. The ratio 𝑅 is independent of the number of transport particles (i.e., the number of events) and material mass. This relation quantifies the effect in transmutation rates of the different neutrons spectrum emitted of each target.
Science and Technology of Nuclear Installations
9
(Hg, Na, Am)
0.30
(Hg, Na, Cm)
0.25
0.25
0.20
0.20 0.15 B
B 0.15
0.10 0.10 0.05
0.05 0.00
0
15
30
45
0.00
60 75 90 105 120 135 150 Distance (cm) (Hg, Na, Np)
0.45
0
15
30
45
60 75 90 105 120 135 150 Distance (cm) (Hg, Na, Pu)
0.60
0.40 0.50 0.35 0.30
0.40
0.25
B
B 0.30
0.20 0.15
0.20
0.10 0.10 0.05 0.00
0
15
30
45
0.00
60 75 90 105 120 135 150 Distance (cm)
B-fission B-capture
0
B-total (cap + fiss)
15
30
45
60 75 90 105 120 135 150 Distance (cm)
B-fission B-capture
B-total (cap + fiss)
Figure 10: Relative number of reactions 𝐵 in (Hg, Na, and TRU) as a function of the distance to the TRU. For proton beam energy of 2.0 GeV and 400,000 particles.
Table 2: Neutrons produced per incident proton in target, beam energy 2.0 GeV, and spallation simulation in MCNPX 2.6.0.
Leakage Capture Total
W 39.2 22.6 66.9
Pb 60.5 0.4 60.8
Bi 58.8 0.2 59
PbBi 60 0.3 60.3
Hg 53.4 10.3 63.7
Th 66.2 19.8 85.9
U 91.3 50.7 142
U-3% 110.1 59.9 170
U-Np 3% 93.4 64.1 157.5
U-Pu 3% 127.5 73.7 201.1
U-Am 3% 104.2 32.7 136.9
U-Cm 3% 112.9 32.4 145.3
Table 3: Average value of the ratio of the capture and fission reactions 𝑅 = (𝐵(𝑛,𝑓) /𝐵(𝑛,𝛾) ) for TRU fuels analysed in reference systems. System (Hg, 0, TRU) (PbBi, 0, TRU) (U, 0, TRU)
Material rate-𝑅 = (𝐵(𝑛,𝑓) /𝐵(𝑛,𝛾) ) 𝑅Am
𝑅Cm
𝑅Np
𝑅Pu
2.13 4.99 1.3
5.82 11.06 3.88
2.25 4.77 1.64
14.3 22.43 11.9
10
Science and Technology of Nuclear Installations (PbBi, PbBi, Cm)
(PbBi, PbBi, Am)
0.90
0.80
0.80
0.70
0.70
0.60
0.60 B
0.50
0.50
B 0.40
0.40
0.30
0.30
0.20
0.20
0.10
0.10 0.00
0
15
30
45
0.00
60 75 90 105 120 135 150 Distance (cm) (PbBi, PbBi, Np)
1.40
0
15
30
45
60 75 90 105 120 135 150 Distance (cm) (PbBi, PbBi, Pu)
1.80 1.60
1.20
1.40 1.00 1.20 0.80
1.00
B
B 0.80
0.60
0.60 0.40 0.40 0.20 0.00
0.20 0
15
30
B-fission B-capture
45
60 75 90 Distance (cm)
105 120 135 150
B-total (cap + fiss)
0.00
0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
Figure 11: Relative number of reaction 𝐵 in (PbBi, PbBi, TRU) as a function of the distance to the TRU. For proton beam energy of 2.0 GeV and 400,000 particles.
The results show that Pu followed by Cm is the materials that have the highest ratio fission/capture in a hardened neutron spectrum. The data only quantify the expected result, since all major isotopes among these actinides have fission cross sections above capture for a range of the spectrum above 1 MeV. This information will be important when analysing the same conditions, but with the presence of moderator or coolant. The data only quantify the expected result, since all major isotopes among these actinides analysed have fission cross sections above capture for a range of the spectrum above 1 MeV. This information will be important when analysing the same conditions, but with the presence of moderator or coolant.
The ratio 𝑅 also shows that the spectrum of neutrons emitted by the target of PbBi is harder than other targets. It is due to the low PbBi microscopic scattering cross section [13, 14]. Figure 7 shows the neutron emission spectrum to different targets. These results are consistent with the literature [2, 20, 21] corresponding to the typical energies expected for the process of evaporation, characterized by nuclear temperatures of 𝑇 = 2 ⋅ ⋅ ⋅ 8 MeV that follow Maxwellian distributions [20, 21]. Figure 7 shows also two high peaks: one of them is close to 2 MeV, approximately the value of the fission energy of neutron emission. The second wide peak is due to the process of capture/scattering of the target material. The PbBi target has only one peak centred close to 2 MeV because of
Science and Technology of Nuclear Installations
11
(U, C, Am)
5.00
(U, C, Cm)
9.00
4.50
8.00
4.00
7.00
3.50
6.00
3.00 5.00 B
B 2.50
4.00 2.00 3.00
1.50 1.00
2.00
0.50
1.00
0.00
0
15
30
45
60 75 90 105 120 135 150 Distance (cm) (U, C, Np)
6.00
0.00
0
15
30
45
60 75 90 105 120 135 150 Distance (cm) (U, C, Pu)
10.00 9.00
5.00
8.00 7.00
4.00
6.00 B 3.00
B 5.00 4.00
2.00
3.00 2.00
1.00
1.00 0.00
0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
0.00
0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
Figure 12: Relative number of reactions 𝐵 in (U, C, and TRU) as a function of the distance to the TRU. For proton beam energy of 2.0 GeV and 400,000 particles.
its low scattering and capture cross sections. Figure 8 shows the neutron emission spectrum of the MCNPX standard evaporation and fission source, compared with the spallation spectrum of U target simulated using Bertini model. Due to the material properties, the expected theoretical distribution of neutrons emission spectrum in the temperatures involved is slightly modified with peaks resonant emission, although within the energy range expected. Other reactions such as (𝑛, 𝑥𝑛) were not considered in this study, due to the low emissivity of the target at energies above 4 MeV as shown in Figures 7 and 8. 3.2. Sodium Cooled Fast System. The sodium cooled system and eutectic lead-bismuth cooled system are shown
in Figures 9–11, and Table 4 presents the average value of the ratio 𝑅 in transmutation materials for all systems. The analysis of these figures and ratio 𝑅 show that: (i) Comparing with the reference system, the Na and PbBi systems show opposite behaviour, increasing (𝑛, 𝛾) reactions to Am, Cm, and Np. Only Pu has the fission reactions as the main mechanism of transmutation. (ii) The behaviour of the three systems is very similar and weakly dependent on the selected target, being the transmutation a direct function of moderator and coolant. For target distances larger than 30 cm, it is not possible to observe any distinction between
12
Science and Technology of Nuclear Installations (Hg, C, Am)
(Hg, C, Cm)
3.00
2.00 1.80
2.50
1.60 1.40
2.00
1.20 B 1.50
B 1.00 0.80
1.00
0.60 0.40
0.50
0.20 0.00
0
15
30
45
0.00
60 75 90 105 120 135 150 Distance (cm) (Hg, C, Np)
2.00
0
15
30
45
(Hg, C, Pu)
3.50
1.80
60 75 90 105 120 135 150 Distance (cm)
3.00
1.60 2.50
1.40 1.20
2.00
B 1.00
B 1.50
0.80 0.60
1.00
0.40 0.50
0.20 0.00
0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm)
0.00
B-total (cap + fiss)
0
15
30
45
60 75 90 105 120 135 150 Distance (cm)
B-fission B-capture
B-total (cap + fiss)
Figure 13: Relative number of reactions 𝐵 in (Hg, C, and TRU) as a function of the distance to the TRU. For proton beam energy of 2.0 GeV and 400,000 particles.
Table 4: Average value of the ratio between fission and capture reactions for TRU fuels analysed in Na and PbBi systems. Material rate-𝑅 = (𝐵(𝑛,𝑓) /𝐵(𝑛,𝛾) )
System (Hg, Na, TRU) (PbBi, PbBi, TRU) (U, Na, TRU)
𝑅Am
𝑅Cm
𝑅Np
𝑅Pu
0.034 0.03 0.029
0.96 1.026 0.92
0.024 0.022 0.016
2.06 2.091 2.029
the 𝑅 values in the three systems, indicating that despite of the low power of moderation, Na and PbBi have more influence than the emission spectrum of the target used. Although these results indicate this direction, this conclusion is not as trivial as shown
forward, when performing the same analysis, but using a default neutron fission source. (iii) The eutectic lead-bismuth has behaviour very similar to sodium. Again, despite of the low power of moderation, sodium and lead-bismuth, when compared
Science and Technology of Nuclear Installations
13
(PbBi, C, Am)
4.50
(PbBi, C, Cm)
8.00
4.00
7.00
3.50
6.00
3.00
5.00
2.50 B 4.00
B 2.00
3.00
1.50
2.00
1.00
1.00
0.50 0.00
0
15
30
45
60
75
90
0.00
105 120 135 150
0
15
30
45
Distance (cm) (PbBi, C, Np)
4.50
B
8.00
3.50
7.00
3.00
6.00
2.50
B
5.00
2.00
4.00
1.50
3.00
1.00
2.00
0.50
1.00
0.00
0
15
30
B-fission B-capture
45
60 75 90 105 120 135 150 Distance (cm) B-total (cap + fiss)
(PbBi, C, Pu)
9.00
4.00
60 75 90 105 120 135 150 Distance (cm)
0.00
0
50 B-fission B-capture
100 Distance (cm)
150
200
B-total (cap + fiss)
Figure 14: Relative number of reaction 𝐵 in (PbBi, C, and TRU) as a function of the distance to the TRU. For proton beam energy of 2.0 GeV and 400,000 particles.
to the Reference System, significantly influence the neutrons spectrum and the reactions involved in the transmutation. This behaviour does not indicate significant neutronic differences between the metals in transmutation, although the values of Moderating Power and Moderating Ratio of Na and PbBi are different, as shown in Table 5. 3.3. Graphite Moderated System. The behaviour of the graphite-moderated system (U, C, and TRU) is shown in Figures 12–14. Table 6 shows the ratio 𝑅 for the graphite moderated systems.
The results show the following. (i) The ratios 𝑅 of all targets have values closer than that in the case of Na and PbBi coolants. These data are consistent with expectations due to strong moderation related to the graphite and show that the transmutation of TRU in an ADS thermal system is independent of the chosen target. As it is evidenced in charts, at 30 cm from the target, there is only influence of the moderator on the transmutation. (ii) For Am and Np, for the three systems, the transmutation occurs mainly via (𝑛, 𝛾) reaction, due to
Science and Technology of Nuclear Installations
4.00
0.16
3.50
0.14
3.00
0.12 R = [B(n,f) /B(n,𝛾) ]
R = [B(n,f) /B(n,𝛾) ]
14
2.50 2.00 1.50
0.10 0.08 0.06
1.00
0.04
0.50
0.02 0.00
0.00 15
30
45
60
75 90 105 Distance (cm)
120
135
150
15
R-(U, C, Cm) R-(U, Na, Cm)
30
45
60
75 90 105 120 135 150 Distance (cm)
R-(U, C, Am) R-(U, Na, Am) 0.23
4.00 3.50
0.18 R = [B(n,f) /B(n,𝛾) ]
R = [B(n,f) /B(n,𝛾) ]
3.00 2.50 2.00 1.50
0.13
0.08
1.00 0.03 0.50 −0.02
0.00 15
30
45
60
75 90 105 Distance (cm)
120
135
150
15
R-(U, C, Pu) R-(U, Na, Pu)
30
45
60
75 90 105 120 135 150 Distance (cm)
R-(U, C, Np) R-(U, Na, Np)
Figure 15: The ratio between the rates of fission and capture to (U, C, and TRU) and (U, Na, and TRU) systems at different locations.
Table 5: Some properties of the moderator/coolant used, with Σ𝑠 and Σ𝑎 being the scattering and absorption macroscopic cross sections, respectively. Moderator C Pb Bi PbBi Na H2 O He (a)
Average log. energy decrement 𝜉
Moderating power 𝜉Σ𝑠 (cm−1 )
Moderating ratio 𝜉Σ𝑠 /Σ𝑎
Mean free path 1/Σ𝑡 (cm)
𝑁1 (a)
𝑁2 (b)
0.158 0.0096 0.0095 0.0095 0.0845 0.948 0.425
0.0608 0.0035 0.0024 0.0029 0.0086 3.27 8.93𝐸 − 06
190 0.58 2.4 1.58 0.66 149 45
2.60 2.71 3.91 3.37 8.7 0.29 47619
58 959 969 965 109 9.7 22
115 1894 1914 1905 215 19 43
Number of neutrons collisions to cross the energy range 10 keV–1 eV. Number of neutrons collisions to cross the energy range 2 MeV–0.025 eV.
(b)
Science and Technology of Nuclear Installations
15 0.25
4.00 3.50
0.20 R = [B(n,f) /B(n,𝛾) ]
R = [B(n,f) /B(n,𝛾) ]
3.00 2.50 2.00 1.50 1.00
0.15
0.10
0.05
0.50 0.00
0.00 15
30
45
60
75
90
105
120
135
150
15
30
45
60
Distance (cm)
75
90
105
120
135
150
105
120
135
150
Distance (cm)
R-(Hg, Na, Cm) R-(Hg, C, Cm)
R-(Hg, Na, Am) R-(Hg, C, Am)
0.33
4.50
0.28
4.00
R = [B(n,f) /B(n,𝛾) ]
R = [B(n,f) /B(n,𝛾) ]
3.50 0.23 0.18 0.13 0.08
3.00 2.50 2.00 1.50 1.00
0.03
0.50
−0.02
0.00 15
30
45
60
75
90
105
120
135
150
15
30
45
60
Distance (cm)
75
90
Distance (cm)
R-(Hg, Na, Np) R-(Hg, C, Np)
R-(Hg, Na, Pu) R-(Hg, C, Pu)
Figure 16: The ratio between the fission and capture rates for systems (Hg, C, and TRU) and (Hg, Na, and TRU) at different locations.
Table 6: Average value of the ratio between fission and capture reactions for TRU fuels analysed in C systems. Material ratio-𝑅 = (𝐵(𝑛,𝑓) /𝐵(𝑛,𝛾) )
System (Hg, C, TRU) (PbBi, C, TRU) (U, C, TRU)
𝑅Am
𝑅Cm
𝑅Np
𝑅Pu
0.114 0.114 0.114
3.544 3.554 3.544
0.001 0.001 0.0005
2.06 2.085 2.087
the low cross sections of the main isotopes. On the other hand, the Cm exhibits a different behaviour having higher rates of reactions than the Pu, inversely as occurred for Na and PbBi. The Pu behaves in all systems reactions as the main mechanism of transmutation, direct function of the high cross sections of its isotopes.
(iii) The results show that the rates of transmutation are independent of target material used. The TRU transmutation depends directly on moderator or coolant used in the system. These results are valid for any target and any energy proton beam incident [13, 14]. Under the conditions analyzed, the PbBi and Na coolants have very similar neutronic behaviours.
16
Science and Technology of Nuclear Installations 0.25
4.00 3.50
0.20
R = [B(n,f) /B(n,𝛾) ]
R = [B(n,f) /B(n,𝛾) ]
3.00 2.50 2.00 1.50
0.15
0.10
1.00 0.05 0.50 0.00
0.00 15
30
45
60
75 90 105 Distance (cm)
120
135
150
15
R-(PbBi, PbBi, Cm) R-(PbBi, C, Cm)
45
60
75 90 105 Distance (cm)
120
135
150
120
135
150
R-(PbBi, PbBi, Am) R-(PbBi, C, Am)
0.38
4.50
0.33
4.00 3.50 R = [B(n,f) /B(n,𝛾) ]
0.28 ] R = [B(n,f) /B(n,𝛾)
30
0.23 0.18 0.13 0.08
3.00 2.50 2.00 1.50 1.00
0.03
0.50
−0.02
0.00 15
30
45
60
75 90 105 Distance (cm)
120
135
150
R-(PbBi, PbBi, Np) R-(PbBi, C, Np)
15
30
45
60
75 90 105 Distance (cm)
R-(PbBi, PbBi, Pu) R-(PbBi, C, Pu)
Figure 17: The ratio between the fission and capture rates for systems (PbBi, C, and TRU) and (PbBi, PbBi, and TRU) at different locations.
A comparison between the 𝑅 for the U, Hg, and PbBi target, in the fast and thermal systems, is shown in Figures 15–17 for all TRU, while Figure 18 shows the values of total reaction rates 𝐵𝑇 (𝐵𝑇 = 𝐵𝑓 + 𝐵𝛾 ). The results show that: (i) Pu—it presents a higher rate of (𝑛, 𝑓) than (𝑛, 𝛾) reactions for all systems. It is due to the significant amount of fissile isotopes present, ≈69%, in the composition of Pu reprocessed, mainly the isotopes 239 Pu and 241 Pu (Table 1).
(ii) Np—its main transmutation mechanism is (𝑛, 𝛾) reactions. The main isotope, 237 Np, has a fission cross section for the thermal spectrum very low ∼1 b, against a high capture cross section of 10,000 b. Such as quantified in charts, its reaction rates are small in all conditions with low transmutation. (iii) Am—it has the lowest values for the overall reaction rates 𝐵𝑇 . Its main transmutation mechanism is via (𝑛, 𝛾) reactions, which is due to high capture cross section in the thermal spectrum of the main
Science and Technology of Nuclear Installations
17 (U, C, TRU) 10
0.9
9
0.8
8
0.7
7 B-total (cap + fiss)
B-total (cap + fiss)
(U, Na, TRU) 1
0.6 0.5 0.4
6 5 4
0.3
3
0.2
2
0.1
1 0
0 15
30
Pu Np
45
60
75 90 105 120 135 150 Distance (cm) Cm Am
15
30
45
60
75 90 105 Distance (cm)
Pu Np
120
135
150
Cm Am
Figure 18: 𝐵𝑇 TRU for the Na and C systems, as a function of the distance.
the isotopes 241 Am (∼34%) and 243 Am (∼65%), whose values are 50,000 b and 5,000 b, respectively. (iv) Cm—it has features favourable to transmutation in a thermal system; the value of 𝐵𝑇 is lesser for all actinides in a system without moderator, and higher than for Pu in a graphite-moderated thermal system. The Cm has three main isotopes with representative mass in the spent fuel, 242 Cm (∼12%), 244 Cm (∼82%), and 245 Cm (∼5%). 242 Cm and 245 Cm are fissile with high fission and capture cross section for the thermal spectrum, while the 244 Cm (with the largest mass percentage) has low fission cross section (∼80 b) and high capture cross section (∼1,000 b). 3.4. Fission Source. Until now the results were obtained based on simulation of the source of protons focusing on a target, generating neutrons by spallation reaction and analyzing how the neutron spectrum emitted interferes with transmuting layer. In this section it considers only the neutrons generated by a typical distribution of fission energy (fission source default available in MCNPX 2.6.0 for criticality calculation), allowing us to compare quantitatively the importance and difference between these two coexisting types of sources in a hybrid ADS-fission system. This comparison is performed by means of fission and capture reaction rates, in the uranium target system. The relative error was found to be less than 3% for all cells in all systems examined. The comparison between sources is better represented by the average value of the ratio between the capture and fission reactions 𝑅, which is independent of the number of transport
particles (i.e., the number of events) and material mass. This relation 𝑅 allows us to infer the effect in transmutation rates of fission source in each TRU as shown in Figures 19–22. The results show the spacing between the curves of the fission source and spallation reaction in the reference system, showing that despite the similarity of the neutron emission spectrum of these sources, the behavior is significantly different when analyzing interactions with actinides. This behavior is observed in Cm and Pu in the Na system, while the Am and Np approach its curves moving away from the center of source origin. In graphite moderate systems, the curves are very close to each other due to the strong influence of moderator on the neutron spectrum. The results indicate that not considering the spallation source in the simulations involving the depletion of fuel in fast systems can lead to significant errors in the results. The result shows that the reaction of spallation, with a hardening neutron spectrum, does not favor the increase in transmutation rates of TRU. In all systems examined, and in all TRU matrices, the ratios 𝑅 are higher for the standard fission source, which indicates that the transmutation mechanism of direct process—fission reaction—is more statistically favorable to occur in a typical critical system, as summarized in Table 7. It is important to note that despite the 𝐵 reaction rates being higher than in the moderate system, it is incorrect to infer that thermal systems have higher reaction rates and thus greater efficiency as mentioned in [15]. The behavior displayed in the graphics for the reaction rates only reflects that the cross sections are higher in the thermal spectrum. Fast systems compensate the lower cross section in the spectrum
18
Science and Technology of Nuclear Installations 0.45
3.00
0.4 2.50
2.00
] R = [B(n,f) /B(n,𝛾)
R = [B(n,f) /B(n,𝛾) ]
0.35
1.50
1.00
0.3 0.25 0.2 0.15 0.1
0.50 0.05 0
0.00 15
30
45
60
75 90 105 Distance (cm)
120
135
150
15
30
45
60
75 90 105 Distance (cm)
120
135
150
kcode -(U, Na, Am) sdef -(U, Na, Am)
kcode -(U, 0, Am) sdef -(U, 0, Am)
0.14 0.12
] R = [B(n,f) /B(n,𝛾)
0.1 0.08 0.06 0.04 0.02 0 15
30
45
60
75 90 105 Distance (cm)
120
135
150
kcode -(U, C, Am) sdef -(U, C, Am)
Figure 19: The ratio between fission and capture rates for systems (U, 0, and Am), (U, Na, and Am), and (U, C, and Am), at the different locations in 𝑥-axis. Comparison between fission source and spallation reaction models.
Table 7: Comparison of the average value of the ratio of the fission and capture reactions 𝑅 = (𝐵(𝑛,𝑓) /𝐵(𝑛,𝛾) ) for TRU fuels analysed in uranium systems. System
Material rate-𝑅 = (𝐵(𝑛,𝑓) /𝐵(𝑛,𝛾) ) 𝑅Am
𝑅Cm
𝑅Np
𝑅Pu
2.14 1.3
7.51 3.88
3.06 1.64
23.12 11.9
(U, Na, TRU)Fission (U, Na, TRU)Spallation
0.5 0.029
2.85 0.92
0.05 0.016
6.43 2.029
(U, C, TRU)Fission (U, C, TRU)Spallation
0.11 0.114
3.8 3.544
0.0003 0.0005
2.46 2.087
(U, 0, TRU)Fission (U, 0, TRU)Spallation
Science and Technology of Nuclear Installations
19
10.00
4.5
9.00
4
8.00
3.5 R = [B(n,f) /B(n,𝛾) ]
R = [B(n,f) /B(n,𝛾) ]
7.00 6.00 5.00 4.00 3.00
3 2.5 2 1.5
2.00
1
1.00
0.5 0
0.00 15
30
45
60
75 90 105 120 135 150 Distance (cm)
15
30
45
60
75 90 105 Distance (cm)
120
135
150
kcode -(U, Na, Cm) sdef -(U, Na, Cm)
kcode -(U, 0, Cm) sdef -(U, 0, Cm)
5 4.5 4 ] R = [B(n,f) /B(n,𝛾)
3.5 3 2.5 2 1.5 1 0.5 0 15
30
45
60
75 90 105 Distance (cm)
120
135
150
kcode -(U, C, Cm) sdef -(U, C, Cm)
Figure 20: The ratio between fission and capture rates for systems (U, 0, and Cm), (U, Na, and Cm), and (U, C, and Cm), at the different locations in 𝑥-axis. Comparison between fission source and spallation reaction models.
and leakage with higher flows of neutrons by increasing the fuel enrichment. A fully closed loop based on thermal spectrum reactors presents difficulties P & T because there is a limit for the number of fuel reprocessing cycles due to the percentage increase of nonfissile isotopes. As an example, Figure 23 shows the burnup of the Pu matrix as fuel for sodium and reference systems. There is a drop in all Pu isotopes; however, the graphite moderated system exhibits a growth of 240 Pu and 242 Pu, fact related to very low fission cross sections for these nuclides to thermal neutrons, less than 1 b.
4. Conclusion This study evaluated how the spectrum of neutrons emitted by different spallation targets, Hg, U, and PbBi, interacts with a layer composed of transuranic elements. The goal was to analyse the rates of reactions of fission and capture involved in the purpose of transmutation. The results allowed us to quantify and consider the mechanisms involved when neutron fluxes from spallation target, with or without the presence of a moderation medium, interact with the transmutation layer.
Science and Technology of Nuclear Installations 4.50
0.9
4.00
0.8
3.50
0.7
3.00
0.6
R = [B(n,f) /B(n,𝛾) ]
R = [B(n,f) /B(n,𝛾) ]
20
2.50 2.00 1.50
0.5 0.4 0.3
1.00
0.2
0.50
0.1
0.00
0 15
30
45
60
75 90 105 120 135 150 Distance (cm)
15
30
45
60
75 90 105 Distance (cm)
120
135
150
kcode -(U, Na, Np) sdef -(U, Na, Np)
kcode -(U, 0, Np) sdef -(U, 0, Np)
0.06
R = [B(n,f) /B(n,𝛾) ]
0.05
0.04
0.03
0.02
0.01
0 15
30
45
60
75 90 105 Distance (cm)
120
135
150
kcode -(U, C, Np) sdef -(U, C, Np)
Figure 21: The ratio between fission and capture rates for systems (U, 0, and Np), (U, Na, and Np), and (U, C, and Np), at the different locations in 𝑥-axis. Comparison between fission source and spallation reaction models.
In the presence of a medium with low or high power moderation, there is no significant difference in the behaviour of the TRU transmutation using different targets. Transmutation in ADS system is independent of the target being a direct function of the medium, moderator/coolant used. However, in the simulations of burning fuel, the insertion of a spallation spectrum in addition to fission spectrum has relevance especially in fast systems. Spallation simulations involving PbBi and Na suggest that these coolants are neutron equivalent. The results do not indicate that the distribution
of PbBi neutron spectrum has advantages front to the Na coolant. Capture reactions in Am and Np, in environments with moderator/coolant, are the main mechanism of transmutation of these elements, and they presented, on average, lower overall reaction rate, making it difficult to reduce their inventory. In turn, the transmutation of the Cm is favoured in the moderated system having the fission as the main transmutation mechanism. About the Pu, it is easily converted in all systems having high reaction rates, with
Science and Technology of Nuclear Installations
21
25.00
12
10
R = [B(n,f) /B(n,𝛾) ]
R = [B(n,f) /B(n,𝛾) ]
20.00
15.00
10.00
5.00
8
6
4
2
0.00
0 15
30
45
60
75
90
105 120 135 150
15
30
45
60
Distance (cm)
75
90
105
120
135
150
Distance (cm)
kcode -(U, 0, Pu) sdef -(U, 0, Pu)
kcode -(U, Na, Pu) sdef -(U, Na, Pu) 3
R = [B(n,f) /B(n,𝛾) ]
2.5
2
1.5
1
0.5
0 15
30
45
60
75 90 105 Distance (cm)
120
135
150
kcode -(U, C, Pu) sdef -(U, C, Pu)
Figure 22: The ratio between fission and capture rates for systems (U, 0, and Pu), (U, Na, and Pu), and (U, C, and Pu), at the different locations in 𝑥-axis. Comparison between fission source and spallation reaction models.
fission reactions being higher than the capture reactions. However, mechanisms of transmutation via capture reaction may not be the best solution for reducing the isotopic inventory. The conversion of these isotopes with lower fission cross section may make the implementation of closed fuel cycles based on thermal systems difficult. The efficiency of the reactor is most closely connected to the fuel burnup, the amount of energy generated per mass of fuel. The higher the reactor burnup, higher rates of reaction involved, thus increasing the chances of interaction between the neutron flux and low cross sections TRU isotopes along the burnup cycle. This is due to the greater probability
of direct reactions (𝑛, 𝑓) or by indirect reactions by (𝑛, 𝛾) successive absorptions and decay chains. The most hardening spallation spectrum in a ADS-fission hybrid system does not seem to have advantages transmutation in relation to a standard critical system. The advantages of recycling in ADS systems appear to be limited security issues, the flexibility of fuel combinatons, and the application dedicated to inventory reduction. As shown in [2], systems that combine closed fuel cycle and fast reactors have greater reduction in the inventory of spent fuel as well as economic advantages. From the point of view analyzed in this work, the results are corroborated herewith.
22
Science and Technology of Nuclear Installations
545
Pu-238
Pu-239
13100 13050
540 13000 12950 Mass (g)
Mass (g)
535 530 525
12900 12850 12800 12750
520 12700 515
5420
0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360
0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360
12650 Time (days)
Time (days)
Pu-240
5400
2500
5380
2450
5360
Mass (g)
Mass (g)
Pu-241
2550
5340
2400 2350
5320
5280
2250 0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360
2300
0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360
5300
Time (days)
Time (days) Pu-242
894 892 890
Mass (g)
888 886 884 882 880 878 0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340 360
876
(U,C,Pu) (U,Na,Pu)
Time (days) (U,0,Pu)
Figure 23: Variation of the mass to Pu isotopes, for 360 days of burnup, in the uranium target systems analyzed. MCNPX 2.6.0 standard fission source.
Science and Technology of Nuclear Installations
23
Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments The authors are grateful to the Brazilian Research Funding Agencies, CNEN (Comiss˜ao Nacional de Energia Nuclear), CNPq (Conselho Nacional de Desenvolvimento e Pesquisa), CAPES (Coordenac¸a˜o de Aperfeic¸oamento de Pessoal de N´ıvel Superior), and FAPEMIG (Fundac¸a˜o de Amparo a` Pesquisa do Estado de Minas Gerais), for the support.
[12]
[13]
[14]
References [1] R. L. Garwin and G. Charpak, Megawatts and Megatons a Turning Point in the Nuclear Ages? Alfred A. Knopf, New York, NY, USA, 2001. [2] Nuclear Energy Agency Organisation for Economic Cooperation and Development, Accelerator-Driven Systems (ADS) and Fast Reactors (FR) in Advanced Nuclear Fuel Cycles: A Comparative Study, Nuclear Energy Agency Organisation for Economic Co-operation and Development, Paris, France, 2002. [3] G. P. Barros, C. Pereira, M. A. F. Veloso, and A. L. Costa, “Thorium and reprocessed fuel utilization in an acceleratordriven system,” Annals of Nuclear Energy, vol. 1, pp. 1–7, 2015. [4] R. V. Sousa, L. P. A. R. Tanure, D. F. Costa, C. Pereira, M. A. F. Veloso, and A. H. De Oliveira, “A preliminary neutronic evaluation of the high temperature nuclear reactor (HTTR) using reprocessed fuel,” Annals of Nuclear Energy, vol. 65, pp. 232–238, 2014. [5] A. Fortini, F. B. A. Monteiro, M. E. Scari et al., “Recent advances on the use of reprocessed fuels and combined thorium fuel cycles in HTR systems,” Progress in Nuclear Energy, 2014. [6] C. E. Velasquez, G. P. de Barros, C. Pereira, M. A. F. Veloso, and A. L. Costa, “Axial neutron flux evaluation in a Tokamak system: a possible transmutation blanket position for a fusion-fission transmutation system,” Brazilian Journal of Physics, vol. 42, no. 3-4, pp. 237–247, 2012. [7] F. C. Silva, C. Pereira, M. A. F. Veloso, and A. L. Costa, “Shifting study of a VHTR using reprocessed fuel with various TRISO packing fractions,” Nuclear Engineering and Design, vol. 248, pp. 42–47, 2012. [8] C. Rubbia, J. A. Rubio, S. Buono et al., “Conceptual design of a fast neutron operated high power energy amplifier,” Tech. Rep. CERN/AT/95-44 (ET), European Organization for Nuclear Research (CERN), Meyrin, Switzerland, 1995. [9] C. M. Van Atta, J. D. Lee, and W. Heckrotte, “The electronuclear conversion of fertile to fissile material,” Tech. Rep. UCRL-52144, Lawrance Livermore Laboratory, 1970. [10] C. D. Bowman, E. D. Arthur, P. W. Lisowski et al., “Nuclear energy generation and waste transmutation using an accelerator-driven intense thermal neutron source,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 320, no. 1-2, pp. 336–367, 1992. [11] M. Salvatores, I. Slessarev, A. Tchistiakov, and G. Ritter, “The potential of accelerator-driven systems for transmutation or power production using thorium or uranium fuel cycles,”
[15]
[16]
[17]
[18] [19] [20]
[21]
[22]
Nuclear Science and Engineering, vol. 126, no. 3, pp. 333–340, 1997. A. Rineiski and W. Maschek, Dynamics and Reactivity Control in Accelerator Driven Systems for Nuclear Waste Burning, Nuclear Energy, Forschungszentrum (FZK), Institute for Nuclear and Energy Technologies, Karlsruhe, Germany, 2006. M. Gilberti, Contribuic¸o˜es ao Desenvolvimento de Sistemas SubCr´ıticos Dirigido por Aceleradores Aplicados a` Transmutac¸a˜ o de Rejeitos Altamente Radioativos [Master thesis], Departamento de Engenharia Nuclear, UFMG, Belo Horizonte, Brazil, 2011, (Portuguese). M. Gilberti, C. Pereira, M. A. F. Veloso, A. L. Costa, and C. A. M. Silva, “Evaluation of spallation source for neutrons production by Monte Carlo methods,” in Proceedings of the 8th International Conference on Mathematical Problems in Engineering, Aerospace and Sciences, pp. 1–8, S˜ao Jos´e dos Campos, Brazil, 2010. S. R. Hashemi-Nezhad, R. Brandt, W. Westmeier et al., “Monte Carlo calculations on transmutation of trans-uranic nuclear waste isotopes using spallation neutrons: difference of lead and graphite moderators,” Nuclear Instruments and Methods in Physics Research A, vol. 481, no. 1–3, pp. 547–557, 2002. S. R. Hashemi-Nezhad, R. Brandt, W. Westmeier et al., “Monte Carlo studies of accelerator driven systems: energy and spatial distribution of neutrons in multiplying and non-multiplying media,” Nuclear Instruments and Methods in Physics Research A, vol. 482, no. 1–3, pp. 537–546, 2002. F. Goldenbaum, The Physics of Spallation Processes: Theory, Experiments and Applications, Institut F¨ur Kernphysik: Forschungszentrum J¨ulich Gmbh, 2004. J. J. Duderstadt and L. J. Hamilton, Nuclear Reactor Analysis, John Wiley & Sons, New York, NY, USA, 1976. H. Nifenecker, “Hybrid nuclear reactors,” Progress in Particle and Nuclear Physics, no. 43, pp. 683–827, 1999. F. Goldenbaum, The Physics of Spallation Processes: Theory, Experiments and Applications, Institut F¨ur Kernphysik, Forschungszentrum J¨ulich Gmbh, 2004. K. Tucek, Neutronic and Burnup Studies of Accelerator-driven Systems Dedicated to Nuclear Waste Transmutation, Department of Physics, Royal Institute of Technology, Stockholm, Sweden, 2004. Committee on Separations Technology and Transmutation Systems Board on Radioactive Waste Management Commission on Geosciences and Environment and Resources National Research Council, Nuclear Wastes: Technologies for Separations and Transmutation, National Academy Press, Washington, DC, USA, 1996.
Journal of
Wind Energy
Hindawi Publishing Corporation http://www.hindawi.com
Journal of
International Journal of
Rotating Machinery
The Scientific World Journal Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Energy
Advances in
Mechanical Engineering Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014
Journal of
Journal of
Industrial Engineering
Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Petroleum Engineering
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014
Journal of
Solar Energy
Submit your manuscripts at http://www.hindawi.com
Engineering Journal of
Hindawi Publishing Corporation http://www.hindawi.com
Advances in
High Energy Physics Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Volume 2014
Journal of
Structures Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
International Journal of
Volume 2014
Advances in
Advances in
Volume 2014
Science and Technology of
Tribology
Hindawi Publishing Corporation http://www.hindawi.com
Hindawi Publishing Corporation http://www.hindawi.com
Renewable Energy
Power Electronics Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Journal of
Nuclear Energy
Hindawi Publishing Corporation http://www.hindawi.com
Photoenergy
Fuels
Combustion Volume 2014
International Journal of
Journal of
Journal of
Volume 2014
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Nuclear Installations Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
Hindawi Publishing Corporation http://www.hindawi.com
Volume 2014
