# Study of an ADS Loaded with Thorium and Reprocessed Fuel

Authors: Graiciany de Paula Barros; Claubia Pereira; Maria A. F. Veloso; Antonella L. Costa
Journal: Science and Technology of Nuclear Installations
 (2012)

---

## Abstract

Accelerator-driven systems (ADSs) are investigated for long-lived fission product transmutation and fuel regeneration. The aim of this paper is to investigate the nuclear fuel evolution and the neutronic parameters of a lead-cooled accelerator-driven system used for fuel breeding.  The fuel used in some fuel rods wasT232hO2 for U233 production. In the other fuel rods was used a mixture based upon Pu-MA, removed from PWR-spent fuel, reprocessed by GANEX, and finally spiked with thorium or depleted uranium. The use of reprocessed fuel ensured the use of T232hO2 without the initial requirement of U233 enrichment. In this paper was used the Monte Carlo code MCNPX 2.6.0 that presents the  depletion/burnup capability, combining an ADS source and kcode-mode (for criticality calculations). The multiplication factor (keff) evolution, the neutron energy spectra in the core at BOL, and the nuclear fuel evolution during the burnup were evaluated. The results indicated that the combined use of T232hO2 and reprocessed fuel allowed U233 production without the initial requirement of U233 enrichment.


---

## Body

## 1. Introduction

In recent years great interest has been given to accelerator-driven systems (ADSs). This is mainly because of their inherent safety features, their waste transmutation potential, and their capability to breed the required233U when the thorium fuel is used. ADSs are useful for recycling of americium, curium, neptunium, and plutonium. A great number of works on the ADS and the relative neutronics have been reported in the scientific literature [1–7]. Pioneers in this revival have been Furukawa et al. [8], Bowman et al. [9, 10], and Rubbia et al. [11, 12]. Very similar idea was first presented almost 60 years ago [13]. In an accelerator-driven system, an accelerator is coupled to a subcritical core loaded with nuclear fuel. The particles accelerated are injected into a spallation target that produces neutrons, which are used in the subcritical core for the fission chain maintenance.Both critical reactors and subcritical accelerator-driven systems are potential candidates as dedicated transmutation systems. Critical reactors, however, loaded with fuel containing large amounts of MA have safety problems caused by unfavorable reactivity coefficients and small delayed neutron fraction. Nevertheless, subcritical systems present criticality and safety advantages: they operate with a neutron source rate that can be increased to compensate for the negative reactivity or any variation during the operation. Therefore, research interest in spent nuclear fuel transmutation has focused on both accelerator-driven systems and fusion-driven systems. At DEN/UFMG we began about 2006 and some results were obtained [14–16]. In this work, the focuses are 233U production and transuranic transmutation possibility in accelerator-driven systems.The main characteristic of ADS (i.e., subcriticality) is particularly favorable and allows a maximum transmutation rate while operating in a safe manner. An advantage of accelerator-driven systems is that, since there is no criticality condition to satisfy, almost any fuel composition can be used in the system [17]. Moreover, ADS could transmute 99.9% of transuranics light-water-reactor-  (LWR-) spent nuclear fuel [18].Another fuel proposed for using in ADS is the thorium-based fuel. There are many reasons for the resurgence of interest in the thorium fuel cycle nowadays. The thorium’s abundance is about three times more than uranium abundance. Thorium-fueled reactor is also an attractive tool to produce long-term nuclear energy with low radiotoxicity waste [19]. The main issue verified in the usage of this fuel in subcritical systems is the need of initial enrichment (233U), since the use of natural thorium (232Th) is not feasible due to the very low values of achieved criticality.Based on those questions, this work presents a preliminary study of ADS for simultaneous use of thorium fuel and reprocessed fuel. The simultaneous usage of thorium fuel and spent fuel allows one to get appropriate values of criticality without the need of enrichment of thorium and still can reduce radiotoxicity of spent fuel by transmutation. The proposal is the simulation of a core with some fuel elements of232Th and other elements based on fuel spent from LWRs. There was simulated a cylindrical core of 12.0 m3 filled with a hexagonal lattice formed by 156 fuel rods, the fuel used in 36 rods was 232ThO2, and the fuel used in the other 120 rods was reprocessed fuel spiked with thorium or depleted uranium. The reprocessed fuel was obtained from spent fuel discharged from the Brazilian PWR ANGRA-I, reprocessed by GANEX processes, and spiked with thorium or depleted uranium. Different percentages of ThO2 or UO2 were added to the reprocessed fuel and analyzed, as the use of only reprocessed fuel generates high criticality values. The neutron spectra were evaluated at BOL (Begin Of Life) and the criticality and the fuel evolution were investigated during ten years with the system operating at steady state, at 515 MW of thermal power. During this time, the 233U production and the reprocessed fuel evolution were quantified.MCNPX 2.6.0 code [20] was used to simulate the geometrical and operational characteristics of the system. For the simulation there was used a combination of spallation source (external source-SDEF) and kcode-mode for the calculation of initial keff and flux, describing in this way the real behavior of an ADS; that is, the initial criticality is the sum of neutrons produced by fission in the fuel and the neutrons produced by spallation in the target. During the burnup the code does not take into account the flux from external source, so the keff and fuel evolution results obtained are just approximations. Moreover, the keffresults that were obtained during the burnup (kcode-mode) provide an insight of the necessary contribution of the external source during the operation period.The isotopic composition of the reprocessed fuel used in this simulation was obtained using GANEX process (Group ActiNide EXtration) which is currently in testing phase [21]. The GANEX process developed by CEA for the reprocessing of Generation IV spent nuclear fuels is composed of two extraction cycles following the dissolution of the spent fuel. Once the uranium is selectively extracted from the dissolution solution by a monoamide solvent, the transuranic elements (Np, Pu, Am, and Cm) are separated from the fission products in a second cycle prior to the co-conversion step [22].The GANEX first cycle recovers more than 99.99% of the total amount of uranium. In a second cycle neptunium, plutonium, americium, and curium are recovered altogether in one liquid flow (actinide product) and the losses of transuranics in the different outputs and in the solvent were estimated at a value lower than 0.5% (neptunium essentially) at the end of the test, corresponding to a recovery yield of actinides higher than 99.5%. Nevertheless, the decontamination factors versus some lanthanides (especially Nd, Sm, and Eu) were much lower than expected and the mass of lanthanides in the actinide product was around 5% at the end [23].


## 2. Methodology

### 2.1. Computational Tool

The MCNPX 2.6.0 code was used to simulate the geometrical and operational characteristics of the system. Such version is quite interesting for the ADS evaluation because it describes the nuclear fuel evolution during the operation. For the calculation of the initialkeffand flux there was used a combination of spallation source (external source-SDEF) and kcode-mode. The code does not take into account the flux from external source during the burnup, so the results obtained are just approximations.The depletion/burnup capability is based on CINDER90. MCNPX depletion is a linked process involving steady-state flux calculations by MCNPX and nuclide depletion calculations by CINDER90. The code runs a steady-state calculation to determine the system eigenvalue, 63-group fluxes, energy-integrated reaction rates, fission multiplicity, and recoverable energy per fission (Q values). CINDER90 then takes those MCNPX-generated values and performs the depletion calculation to generate new number densities for the next time step. MCNPX takes these new number densities and generates another set of fluxes and reaction rates. The process repeats itself until after the final time step specified by the user [24].


### 2.2. System Parameters

Figures1 and 2 show schematic views of the simulated ADS. The basic geometry includes the spallation target, a subcritical core, and the reflector. The accelerator tube has a radius of 1.5 cm, and the axial position is in the center of the target. The spallation source is represented by a neutrons source with a spectrum characteristic of spallation reactions. Such spectrum was generated, in a previous simulation, using a beam of 1-GeV protons with a parabolic spatial profile. There was used Bertini intranuclear cascade model for the transport of protons, neutrons, and charged pions. The parameters of this simulation were described in [15].Figure 1
Horizontal cross section of the ADS simulated.Figure 2
Vertical cross section of the ADS simulated.Due to its high neutron yield by spallation reactions there was used lead as spallation target. The Pb cylindrical target has 9.5 cm radius and 39 cm height. Lead was also used as coolant and as reflector. The use of lead as coolant offers many advantages like convective cooling, passive safety and small neutron absorption cross section [25].Another advantage of the lead coolant is that it is not a neutron moderator. This is important because the protactinium effect, which limits the achievable values ofkeff, is less severe for harder spectra. For solid fuels, systems without moderator and based on thorium,  smaller values of capture cross sections of fission products will reduce the keff variation and produce a hardening (shift to higher energies) in the neutron energy distribution. So, the inventory of 233U is much larger in fast reactors (about 7 times), with the associated larger breeding times and inventory radiotoxicity [26].The simulated core is a cylinder of 12.0 m3 filled with a hexagonal lattice formed by 156 fuel rods. The fuel rod diameter is 6 cm, the pitch is 12 cm, and rod length is 400 cm. The fuel used in 36 rods is 232ThO2. The fuel used in the other 120 rods is spent fuel discharged from reactor ANGRA-I, with initial enrichment of 3.1%. This fuel was burned in ORIGEN 2.1 code [27] for three cycles, with the burnup of approximately 11.000 MW d/t in each cycle, following the ANGRA-I power historic of real cycles 1, 2, and 3. The composition of the spent fuel used is shown in Table 1 [28], after remaining 5 years in a cooling pool.Table 1
UO2 fuel composition after 33.000 MW d/T.
ActinidesFission ProductsTotalweightfraction=0.9785Totalweightfraction=0.0215NuclideWeight fractionNuclideWeight fractionNuclideWeight fraction234U1.546E−04H1.525E−06Ag2.386E−03235U8.046E−03Li5.395E−09Cd3.267E−03236U4.113E−03Be4.384E−09In6.975E−05238U9.781E−01C7.707E−10Sn2.686E−03237Np4.759E−04Co2.672E−17Sb8.428E−04238Pu1.851E−04Ni9.566E−15Te1.413E−02239Pu4.847E−03Cu8.460E−14I7.037E−03240Pu1.657E−03Zn1.186E−09Xe1.527E−01241Pu1.558E−03Ga6.593E−10Cs8.153E−02242Pu5.888E−04Ge1.929E−05Ba4.294E−02243Am1.126E−04As5.897E−06La3.566E−02Others1.620E−04Se1.635E−03Ce7.546E−02Total 1.000E+00Br6.294E−04Pr3.211E−02OBS.: “Others” include:Kr1.058E−02Nd1.120E−014He1.202E−06Rb9.934E−03Pm3.342E−03230Th3.407E−09Sr2.512E−02Sm2.246E−02233U2.069E−09Y1.339E−02Eu4.776E−03237U5.860E−06Zr1.031E−01Gd2.583E−03238Np7.797E−07Nb6.179E−04Tb8.146E−05239Np4.923E−05Mo9.596E−02Dy4.039E−05241Am8.292E−05Tc2.262E−02Ho4.368E−06242Am1.526E−07Ru6.843E−02Er1.724E−06242Cm2.582E−05Rh1.275E−02Tm1.812E−09244Cm2.960E−05Pd3.906E−02Yb3.676E−10245Cm1.030E−06Total      1.000E+00After cooling by five years, the spent nuclear was submitted to GANEX reprocessing and the isotopic composition after the reprocessing is presented in Table2. The amount of uranium after the reprocessing is 0.01% of the total amount of uranium in the burned fuel. In terms of actinides, neptunium, plutonium, americium, and curium are recovered with a loss of only 0.5% (neptunium essentially). The lanthanides contamination in the reprocessed fuel is around 5%. The isotopic composition presented in Table 2 was adjusted (with a computational program) for several percentages of thorium or depleted uranium added to the reprocessed fuel and the final composition (ThO2orUO2+reprocessedfuel) was normalized. The ratio between the number of thorium fuel rods and reprocessed fuel rods and the positions of the thorium fuel rods were determined to maximize the amount of reprocessed fuel and achieve keff values close to 1.Table 2
Fuel composition (normalized) after reprocessing.
ActinidesWeight fraction234U9.29E−07235U4.83E−05236U2.47E−05238U5.88E−03237Np2.72E−02238Pu1.11E−02239Pu2.91E−01240Pu9.95E−02241Pu9.36E−02242Pu3.54E−02243Am6.76E−03233U1.20E−11237U3.52E−08238Np4.45E−05239Np2.81E−03241Am4.98E−03242Am9.17E−06242Cm1.55E−03244Cm1.78E−03245Cm6.19E−05Nd3.36E−01Sm6.75E−02Eu1.43E−02The isotopic compositions used in the simulations with several percentages of thorium or depleted uranium added to the reprocessed fuel are shown in Tables3 and 4.Table 3
Fuel composition (normalized) for the reprocessed fuels spiked with thorium.
ActinidesWeight fraction fuel spiked with 60% of thoriumWeight fraction fuel spiked with 61% of thoriumWeight fraction fuel spiked with 62% of thorium234U3.7145E−0073.6216E−0073.5287E−007235U1.9331E−0051.8848E−0051.8365E−005236U9.8820E−0069.6349E−0069.3879E−006238U0.00235000.00229130.0022325237Np0.0108620.0105910.010319238Pu0.00444730.00433610.0042249239Pu0.116460.113540.11063240Pu0.0398110.0388160.037821241Pu0.0374330.0364970.035561242Pu0.0141470.0137930.013439243Am0.00270540.00263770.0025701233U4.8052E−0124.6851E−0124.5650E−012237U1.4079E−0081.3727E−0081.3375E−008238Np1.7797E−0051.7352E−0051.6907E−005239Np0.00112370.00109560.0010675241Am0.00199230.00194240.0018926242Am3.6664E−0063.5747E−0063.4831E−006242Cm0.000620360.000604850.00058934244Cm0.000711180.000693400.00067562245Cm2.4747E−0052.4128E−0052.3510E−005Nd0.134550.131180.12782Sm0.0269810.0263070.025632Eu0.00573750.00559400.0054506232Th0.200000.203330.20667O0.400000.406670.41333Table 4
Fuel composition (normalized) for the reprocessed fuels spiked with depleted uranium.
ActinidesWeight fraction fuel spiked with 62% of depleted uraniumWeight fraction fuel spiked with 63% of depleted uraniumWeight fraction fuel spiked with 64% of depleted uranium234U3.5287E−0073.4359E−0073.3430E−007235U1.8365E−0051.7882E−0051.7398E−005236U9.3879E−0069.1408E−0068.8938E−006238U0.00223250.00217380.0021150237Np0.0103190.0100480.0097761238Pu0.00422490.00411370.0040025239Pu0.110630.107720.10481240Pu0.0378210.0368260.035830241Pu0.0355610.0346250.033690242Pu0.0134390.0130860.012732243Am0.00257010.00250250.0024348233U4.5650E−0124.4449E−0124.3247E−012237U1.3375E−0081.3023E−0081.2671E−008238Np1.6907E−0051.6462E−0051.6017E−005239Np0.00106750.00103940.0010113241Am0.00189260.00184280.0017930242Am3.4831E−0063.3914E−0063.2998E−006242Cm0.000589340.000573830.00055832244Cm0.000675620.000657840.00064006245Cm2.3510E−0052.2891E−0052.2272E−005Nd0.127820.124460.12109Sm0.0256320.0249580.024283Eu0.00545060.00530720.0051637238U0.206250.209580.21291235U0.000413330.000420000.00042667234U2.0667E−0062.1000E−0062.1333E−006O0.413330.420000.42667


### 2.3. Calculational Procedure

For this simulation there was used a combination of ADS source (SDEF) and kcode-mode for calculation of initialkeff and flux. The SDEF source was positioned on the center of the target. The burnup calculations were performed in kcode-mode using 53 time steps; the total simulated time was 10 years and the thermal power of operation throughout this period was 515 MW t. were performed simulations using reprocessed fuel spiked with thorium in different percentages and simulations using reprocessed fuel spiked with different percentages of depleted uranium.


## 2.1. Computational Tool

The MCNPX 2.6.0 code was used to simulate the geometrical and operational characteristics of the system. Such version is quite interesting for the ADS evaluation because it describes the nuclear fuel evolution during the operation. For the calculation of the initialkeffand flux there was used a combination of spallation source (external source-SDEF) and kcode-mode. The code does not take into account the flux from external source during the burnup, so the results obtained are just approximations.The depletion/burnup capability is based on CINDER90. MCNPX depletion is a linked process involving steady-state flux calculations by MCNPX and nuclide depletion calculations by CINDER90. The code runs a steady-state calculation to determine the system eigenvalue, 63-group fluxes, energy-integrated reaction rates, fission multiplicity, and recoverable energy per fission (Q values). CINDER90 then takes those MCNPX-generated values and performs the depletion calculation to generate new number densities for the next time step. MCNPX takes these new number densities and generates another set of fluxes and reaction rates. The process repeats itself until after the final time step specified by the user [24].


## 2.2. System Parameters

Figures1 and 2 show schematic views of the simulated ADS. The basic geometry includes the spallation target, a subcritical core, and the reflector. The accelerator tube has a radius of 1.5 cm, and the axial position is in the center of the target. The spallation source is represented by a neutrons source with a spectrum characteristic of spallation reactions. Such spectrum was generated, in a previous simulation, using a beam of 1-GeV protons with a parabolic spatial profile. There was used Bertini intranuclear cascade model for the transport of protons, neutrons, and charged pions. The parameters of this simulation were described in [15].Figure 1
Horizontal cross section of the ADS simulated.Figure 2
Vertical cross section of the ADS simulated.Due to its high neutron yield by spallation reactions there was used lead as spallation target. The Pb cylindrical target has 9.5 cm radius and 39 cm height. Lead was also used as coolant and as reflector. The use of lead as coolant offers many advantages like convective cooling, passive safety and small neutron absorption cross section [25].Another advantage of the lead coolant is that it is not a neutron moderator. This is important because the protactinium effect, which limits the achievable values ofkeff, is less severe for harder spectra. For solid fuels, systems without moderator and based on thorium,  smaller values of capture cross sections of fission products will reduce the keff variation and produce a hardening (shift to higher energies) in the neutron energy distribution. So, the inventory of 233U is much larger in fast reactors (about 7 times), with the associated larger breeding times and inventory radiotoxicity [26].The simulated core is a cylinder of 12.0 m3 filled with a hexagonal lattice formed by 156 fuel rods. The fuel rod diameter is 6 cm, the pitch is 12 cm, and rod length is 400 cm. The fuel used in 36 rods is 232ThO2. The fuel used in the other 120 rods is spent fuel discharged from reactor ANGRA-I, with initial enrichment of 3.1%. This fuel was burned in ORIGEN 2.1 code [27] for three cycles, with the burnup of approximately 11.000 MW d/t in each cycle, following the ANGRA-I power historic of real cycles 1, 2, and 3. The composition of the spent fuel used is shown in Table 1 [28], after remaining 5 years in a cooling pool.Table 1
UO2 fuel composition after 33.000 MW d/T.
ActinidesFission ProductsTotalweightfraction=0.9785Totalweightfraction=0.0215NuclideWeight fractionNuclideWeight fractionNuclideWeight fraction234U1.546E−04H1.525E−06Ag2.386E−03235U8.046E−03Li5.395E−09Cd3.267E−03236U4.113E−03Be4.384E−09In6.975E−05238U9.781E−01C7.707E−10Sn2.686E−03237Np4.759E−04Co2.672E−17Sb8.428E−04238Pu1.851E−04Ni9.566E−15Te1.413E−02239Pu4.847E−03Cu8.460E−14I7.037E−03240Pu1.657E−03Zn1.186E−09Xe1.527E−01241Pu1.558E−03Ga6.593E−10Cs8.153E−02242Pu5.888E−04Ge1.929E−05Ba4.294E−02243Am1.126E−04As5.897E−06La3.566E−02Others1.620E−04Se1.635E−03Ce7.546E−02Total 1.000E+00Br6.294E−04Pr3.211E−02OBS.: “Others” include:Kr1.058E−02Nd1.120E−014He1.202E−06Rb9.934E−03Pm3.342E−03230Th3.407E−09Sr2.512E−02Sm2.246E−02233U2.069E−09Y1.339E−02Eu4.776E−03237U5.860E−06Zr1.031E−01Gd2.583E−03238Np7.797E−07Nb6.179E−04Tb8.146E−05239Np4.923E−05Mo9.596E−02Dy4.039E−05241Am8.292E−05Tc2.262E−02Ho4.368E−06242Am1.526E−07Ru6.843E−02Er1.724E−06242Cm2.582E−05Rh1.275E−02Tm1.812E−09244Cm2.960E−05Pd3.906E−02Yb3.676E−10245Cm1.030E−06Total      1.000E+00After cooling by five years, the spent nuclear was submitted to GANEX reprocessing and the isotopic composition after the reprocessing is presented in Table2. The amount of uranium after the reprocessing is 0.01% of the total amount of uranium in the burned fuel. In terms of actinides, neptunium, plutonium, americium, and curium are recovered with a loss of only 0.5% (neptunium essentially). The lanthanides contamination in the reprocessed fuel is around 5%. The isotopic composition presented in Table 2 was adjusted (with a computational program) for several percentages of thorium or depleted uranium added to the reprocessed fuel and the final composition (ThO2orUO2+reprocessedfuel) was normalized. The ratio between the number of thorium fuel rods and reprocessed fuel rods and the positions of the thorium fuel rods were determined to maximize the amount of reprocessed fuel and achieve keff values close to 1.Table 2
Fuel composition (normalized) after reprocessing.
ActinidesWeight fraction234U9.29E−07235U4.83E−05236U2.47E−05238U5.88E−03237Np2.72E−02238Pu1.11E−02239Pu2.91E−01240Pu9.95E−02241Pu9.36E−02242Pu3.54E−02243Am6.76E−03233U1.20E−11237U3.52E−08238Np4.45E−05239Np2.81E−03241Am4.98E−03242Am9.17E−06242Cm1.55E−03244Cm1.78E−03245Cm6.19E−05Nd3.36E−01Sm6.75E−02Eu1.43E−02The isotopic compositions used in the simulations with several percentages of thorium or depleted uranium added to the reprocessed fuel are shown in Tables3 and 4.Table 3
Fuel composition (normalized) for the reprocessed fuels spiked with thorium.
ActinidesWeight fraction fuel spiked with 60% of thoriumWeight fraction fuel spiked with 61% of thoriumWeight fraction fuel spiked with 62% of thorium234U3.7145E−0073.6216E−0073.5287E−007235U1.9331E−0051.8848E−0051.8365E−005236U9.8820E−0069.6349E−0069.3879E−006238U0.00235000.00229130.0022325237Np0.0108620.0105910.010319238Pu0.00444730.00433610.0042249239Pu0.116460.113540.11063240Pu0.0398110.0388160.037821241Pu0.0374330.0364970.035561242Pu0.0141470.0137930.013439243Am0.00270540.00263770.0025701233U4.8052E−0124.6851E−0124.5650E−012237U1.4079E−0081.3727E−0081.3375E−008238Np1.7797E−0051.7352E−0051.6907E−005239Np0.00112370.00109560.0010675241Am0.00199230.00194240.0018926242Am3.6664E−0063.5747E−0063.4831E−006242Cm0.000620360.000604850.00058934244Cm0.000711180.000693400.00067562245Cm2.4747E−0052.4128E−0052.3510E−005Nd0.134550.131180.12782Sm0.0269810.0263070.025632Eu0.00573750.00559400.0054506232Th0.200000.203330.20667O0.400000.406670.41333Table 4
Fuel composition (normalized) for the reprocessed fuels spiked with depleted uranium.
ActinidesWeight fraction fuel spiked with 62% of depleted uraniumWeight fraction fuel spiked with 63% of depleted uraniumWeight fraction fuel spiked with 64% of depleted uranium234U3.5287E−0073.4359E−0073.3430E−007235U1.8365E−0051.7882E−0051.7398E−005236U9.3879E−0069.1408E−0068.8938E−006238U0.00223250.00217380.0021150237Np0.0103190.0100480.0097761238Pu0.00422490.00411370.0040025239Pu0.110630.107720.10481240Pu0.0378210.0368260.035830241Pu0.0355610.0346250.033690242Pu0.0134390.0130860.012732243Am0.00257010.00250250.0024348233U4.5650E−0124.4449E−0124.3247E−012237U1.3375E−0081.3023E−0081.2671E−008238Np1.6907E−0051.6462E−0051.6017E−005239Np0.00106750.00103940.0010113241Am0.00189260.00184280.0017930242Am3.4831E−0063.3914E−0063.2998E−006242Cm0.000589340.000573830.00055832244Cm0.000675620.000657840.00064006245Cm2.3510E−0052.2891E−0052.2272E−005Nd0.127820.124460.12109Sm0.0256320.0249580.024283Eu0.00545060.00530720.0051637238U0.206250.209580.21291235U0.000413330.000420000.00042667234U2.0667E−0062.1000E−0062.1333E−006O0.413330.420000.42667


## 2.3. Calculational Procedure

For this simulation there was used a combination of ADS source (SDEF) and kcode-mode for calculation of initialkeff and flux. The SDEF source was positioned on the center of the target. The burnup calculations were performed in kcode-mode using 53 time steps; the total simulated time was 10 years and the thermal power of operation throughout this period was 515 MW t. were performed simulations using reprocessed fuel spiked with thorium in different percentages and simulations using reprocessed fuel spiked with different percentages of depleted uranium.


## 3. Results

### 3.1. Neutronic Evaluation

In  Figure3 is shown the normalized neutron energy distribution (for the systems spiked with 62% of thorium and 62% of depleted uranium) in the BOL (beginning of life). Both spectra in the thorium fuel present their maximum around 150 keV and around 250 keV in the reprocessed fuel. Therefore, the reprocessed fuel neutron spectrum is harder than thorium fuel spectrum. That is result of the plutonium presence that emits fission neutrons with a harder spectrum. The spectra for the system loaded with fuel spiked with thorium are very similar to spectra for the fuel spiked with depleted uranium. The spectrum in the lead of the core (coolant) presents a similar neutrons energy distribution to fuel spectra, but the number of neutrons available in the fuel is larger due to neutron production by fission reactions.Neutron energy distributions in the (BOL).
(a)(b)Figure4 shows the multiplication factor (keff) as a time function for the systems loaded with 232ThO2 (24 rods) and reprocessed fuel (132 rods) spiked with thorium (Figure 4(a)) and spiked with depleted uranium (Figure 4(b)). The percentages of thorium and depleted uranium in the reprocessed fuel were varied and the keff values obtained are represented in Figure 4. The percentages of thorium or depleted uranium were chosen in such a way that the values for the initial keff were close to 1, since a large variability on those percentages could lead to considerably low initial keff values.Multiplication factor (keff) evolution.
(a)(b)From Figure4 can be observed that after some years of operation the decrease in the keff  value was smaller in the system which used reprocessed fuel spiked with thorium. That behavior indicates that the use of reprocessed fuel spiked with thorium allows an extension of the burnup without reloading fuel.For the same percentage of thorium and depleted uranium (62%), thekeff values for the system spiked with depleted uranium (keff variation of 0.991 to 0.761) were significantly higher than these for the system spiked with thorium (keff variation of 0.964 to 0.782), but the decrease in the kefffor the system with depleted uranium was also higher.The drop of thekeffvalue observed in two cases is mainly due to the increased poisoning caused by the accumulation of fission products having large neutron capture cross sections and, of course, the constant energy generation by the system.


### 3.2. Thorium Fuel Evolution

Figure5 describes the 232Th, 233Pa, and 233U concentrations in the thorium fuel rods during the 10 years of operation for the systems spiked with 62% of thorium and 62% of depleted uranium. It can be verified that for the system spiked with thorium the concentration of the isotope 232Th was reduced by about 17% (623 kg), for the system spiked with depleted uranium that was about 18% (644 kg). It is the result of the considerable capture cross section of that isotope, which allows the 233U production. 233U is formed when 232Th captures a neutron, and it soon undergoes two beta decays:(1)232Th+n⟶Th233→β-Pa233→β-U233.
A relevant question in the neutron spectrum choice is the protactinium effect. Protactinium captures neutrons (due to its large capture cross section) and, thus, decreases the reactivity and fuel regeneration. From Figure 5 can be observed that this isotope is formed in small scale in both systems. Furthermore, the protactinium effect is less severe for harder spectra like that.232Th, 233Pa, and 233U mass variations.
(a)(b)(c)Figure5 also presents the 233U concentration during the burnup, which is formed in large scale (~282 kg) for the system with reprocessed fuel spiked with thorium and (~284 kg) for the system spiked with depleted uranium. That 233U production was due to 232Th regeneration. Then these results show that the combined use of 232ThO2 and reprocessed fuel allowed 233U production without the initial requirement of 233U enrichment.


### 3.3. Reprocessed Fuel Evolution

The long-term potential radiotoxicity of the spent nuclear fuel is mainly associated with transuranics (TRU) particularly Pu and MA: Am, Cm, and Np. These constitute a significant radiological source term over a very long period within a spent fuel [29]. Figures 6–9 show the concentration of those isotopes during the burnup.Plutonium mass variations.
(a)(b)(c)(d)(e)Americium mass variations.
(a)(b)Curium mass variations.
(a)(b)(c)Neptunium mass variations.
(a)(b)In a standard uranium-uranium(235Ux,238U1-x) fuel used in nuclear reactors, the isotope of neptunium, 237Np, originates from 235U, and the isotopes of americium, 241Am and 243Am, and curium, 242Cm, 244Cm, and 245Cm, from 238Uas a result of nuclear reactions [30].Plutonium is the major element for radiotoxicity and mass in the storage [31]. 238Pu, which is generated from237Np, was produced in large scale in the system simulated (Figure 6), while 239Pu and241Pu mass were reduced mainly due to fission reactions. 240Pu is formed when 239Pucaptures a neutron that justifies the high production of that isotope due to the reduced in 239Pu. 242Pu was formed in small scale. The greater plutonium production in the system spiked with depleted uranium was due to the greater uranium availability in this system that allows the 237Np formation and consequently plutonium production.Long-lived isotopes of americium produced are241Am and 243Am. Americium is produced by β-decay of 241Pu (241Pu→241Am) and (n,γ) reaction of 242Pu(242Pu→243Am) during the operation. The mass variation of americium is showed in Figure 7; the difference between the americium productions in the two cases was very small.Long-lived isotopes of curium produced are242Cm, 244Cm, and 245Cm. At that, isotope 242Cm has a relatively short lifetime justifying its low production (Figure 8). The radioactive decay of 242Cm proceeds as follows:(2)242Cm(T1/2=163days;α)⟶Pu238(T1/2=87.7years;α)⟶U234(T1/2=2.46×105years;α).
Curium is the most intensely radioactive of the actinides for both neutron emission and α-activity [32]. The most abundant of its isotopes is 244Cm, which decays with a half-life of 18 years to form 240Pu. Radioactive decay of 244Cm occurs as follows:(3)Cm244(T1/2=18.1years;α)⟶Pu240(T1/2=6.56×103years;α).
The radioactive decay of 245Cm to the nearest long-lived daughter nucleus proceeds as follows:(4)Cm245(T1/2=8.5×103years;α)⟶Pu241(T1/2=14.4years;β)⟶Am241(T1/2=4.33×102years;α).From Figure8 can be observed increase in the masses of 245Cmand 244Cm. At long-term burnup of plutonium, there is the increase of radiotoxicity determined by accumulation of 244Cm. The burnup of plutonium should be done by much shorter lifetimes with intermediate processing and addition of new plutonium. Since 244Cm half-life makes 18.1 years, one can organize a controllable storage of extracted curium [32]. The Cm mass evolution in the two cases was very similar.Neptunium is produced by(n,γ) reaction of 235U, (n,2n) reaction of 238U, and α-decay of 241Am. The most important Np isotope is 237Np (half-life 2.14 × 106 years). The decrease observed in the 237Np mass (Figure 9) is mainly due to 238Np production by neutron capture. The very low concentration of 238Np (Figure 9) is justified by the quick decay (2.117 days) to 238Pu through β decay. Once more was not observed significant difference between the mass variations for each isotope in the two cases.The large amount of232Th in the fuel spiked with 232ThO2 results in a high 233U production (~258 kg), as can be observed in Figure 10. That fuel regeneration potential appears with an advantage of the reprocessed fuel spiked with thorium. Figure 10 also shows the 235U and238Umass variations.Thorium and uranium mass variation.
(a)(b)(c)(d)


## 3.1. Neutronic Evaluation

In  Figure3 is shown the normalized neutron energy distribution (for the systems spiked with 62% of thorium and 62% of depleted uranium) in the BOL (beginning of life). Both spectra in the thorium fuel present their maximum around 150 keV and around 250 keV in the reprocessed fuel. Therefore, the reprocessed fuel neutron spectrum is harder than thorium fuel spectrum. That is result of the plutonium presence that emits fission neutrons with a harder spectrum. The spectra for the system loaded with fuel spiked with thorium are very similar to spectra for the fuel spiked with depleted uranium. The spectrum in the lead of the core (coolant) presents a similar neutrons energy distribution to fuel spectra, but the number of neutrons available in the fuel is larger due to neutron production by fission reactions.Neutron energy distributions in the (BOL).
(a)(b)Figure4 shows the multiplication factor (keff) as a time function for the systems loaded with 232ThO2 (24 rods) and reprocessed fuel (132 rods) spiked with thorium (Figure 4(a)) and spiked with depleted uranium (Figure 4(b)). The percentages of thorium and depleted uranium in the reprocessed fuel were varied and the keff values obtained are represented in Figure 4. The percentages of thorium or depleted uranium were chosen in such a way that the values for the initial keff were close to 1, since a large variability on those percentages could lead to considerably low initial keff values.Multiplication factor (keff) evolution.
(a)(b)From Figure4 can be observed that after some years of operation the decrease in the keff  value was smaller in the system which used reprocessed fuel spiked with thorium. That behavior indicates that the use of reprocessed fuel spiked with thorium allows an extension of the burnup without reloading fuel.For the same percentage of thorium and depleted uranium (62%), thekeff values for the system spiked with depleted uranium (keff variation of 0.991 to 0.761) were significantly higher than these for the system spiked with thorium (keff variation of 0.964 to 0.782), but the decrease in the kefffor the system with depleted uranium was also higher.The drop of thekeffvalue observed in two cases is mainly due to the increased poisoning caused by the accumulation of fission products having large neutron capture cross sections and, of course, the constant energy generation by the system.


## 3.2. Thorium Fuel Evolution

Figure5 describes the 232Th, 233Pa, and 233U concentrations in the thorium fuel rods during the 10 years of operation for the systems spiked with 62% of thorium and 62% of depleted uranium. It can be verified that for the system spiked with thorium the concentration of the isotope 232Th was reduced by about 17% (623 kg), for the system spiked with depleted uranium that was about 18% (644 kg). It is the result of the considerable capture cross section of that isotope, which allows the 233U production. 233U is formed when 232Th captures a neutron, and it soon undergoes two beta decays:(1)232Th+n⟶Th233→β-Pa233→β-U233.
A relevant question in the neutron spectrum choice is the protactinium effect. Protactinium captures neutrons (due to its large capture cross section) and, thus, decreases the reactivity and fuel regeneration. From Figure 5 can be observed that this isotope is formed in small scale in both systems. Furthermore, the protactinium effect is less severe for harder spectra like that.232Th, 233Pa, and 233U mass variations.
(a)(b)(c)Figure5 also presents the 233U concentration during the burnup, which is formed in large scale (~282 kg) for the system with reprocessed fuel spiked with thorium and (~284 kg) for the system spiked with depleted uranium. That 233U production was due to 232Th regeneration. Then these results show that the combined use of 232ThO2 and reprocessed fuel allowed 233U production without the initial requirement of 233U enrichment.


## 3.3. Reprocessed Fuel Evolution

The long-term potential radiotoxicity of the spent nuclear fuel is mainly associated with transuranics (TRU) particularly Pu and MA: Am, Cm, and Np. These constitute a significant radiological source term over a very long period within a spent fuel [29]. Figures 6–9 show the concentration of those isotopes during the burnup.Plutonium mass variations.
(a)(b)(c)(d)(e)Americium mass variations.
(a)(b)Curium mass variations.
(a)(b)(c)Neptunium mass variations.
(a)(b)In a standard uranium-uranium(235Ux,238U1-x) fuel used in nuclear reactors, the isotope of neptunium, 237Np, originates from 235U, and the isotopes of americium, 241Am and 243Am, and curium, 242Cm, 244Cm, and 245Cm, from 238Uas a result of nuclear reactions [30].Plutonium is the major element for radiotoxicity and mass in the storage [31]. 238Pu, which is generated from237Np, was produced in large scale in the system simulated (Figure 6), while 239Pu and241Pu mass were reduced mainly due to fission reactions. 240Pu is formed when 239Pucaptures a neutron that justifies the high production of that isotope due to the reduced in 239Pu. 242Pu was formed in small scale. The greater plutonium production in the system spiked with depleted uranium was due to the greater uranium availability in this system that allows the 237Np formation and consequently plutonium production.Long-lived isotopes of americium produced are241Am and 243Am. Americium is produced by β-decay of 241Pu (241Pu→241Am) and (n,γ) reaction of 242Pu(242Pu→243Am) during the operation. The mass variation of americium is showed in Figure 7; the difference between the americium productions in the two cases was very small.Long-lived isotopes of curium produced are242Cm, 244Cm, and 245Cm. At that, isotope 242Cm has a relatively short lifetime justifying its low production (Figure 8). The radioactive decay of 242Cm proceeds as follows:(2)242Cm(T1/2=163days;α)⟶Pu238(T1/2=87.7years;α)⟶U234(T1/2=2.46×105years;α).
Curium is the most intensely radioactive of the actinides for both neutron emission and α-activity [32]. The most abundant of its isotopes is 244Cm, which decays with a half-life of 18 years to form 240Pu. Radioactive decay of 244Cm occurs as follows:(3)Cm244(T1/2=18.1years;α)⟶Pu240(T1/2=6.56×103years;α).
The radioactive decay of 245Cm to the nearest long-lived daughter nucleus proceeds as follows:(4)Cm245(T1/2=8.5×103years;α)⟶Pu241(T1/2=14.4years;β)⟶Am241(T1/2=4.33×102years;α).From Figure8 can be observed increase in the masses of 245Cmand 244Cm. At long-term burnup of plutonium, there is the increase of radiotoxicity determined by accumulation of 244Cm. The burnup of plutonium should be done by much shorter lifetimes with intermediate processing and addition of new plutonium. Since 244Cm half-life makes 18.1 years, one can organize a controllable storage of extracted curium [32]. The Cm mass evolution in the two cases was very similar.Neptunium is produced by(n,γ) reaction of 235U, (n,2n) reaction of 238U, and α-decay of 241Am. The most important Np isotope is 237Np (half-life 2.14 × 106 years). The decrease observed in the 237Np mass (Figure 9) is mainly due to 238Np production by neutron capture. The very low concentration of 238Np (Figure 9) is justified by the quick decay (2.117 days) to 238Pu through β decay. Once more was not observed significant difference between the mass variations for each isotope in the two cases.The large amount of232Th in the fuel spiked with 232ThO2 results in a high 233U production (~258 kg), as can be observed in Figure 10. That fuel regeneration potential appears with an advantage of the reprocessed fuel spiked with thorium. Figure 10 also shows the 235U and238Umass variations.Thorium and uranium mass variation.
(a)(b)(c)(d)


## 4. Conclusions

The main conclusions indicated by the simulation presented in this work were as follow.(i)
It is noteworthy that the code does not take into account the flux from external source during the burnup, so thekeff and fuel results obtained are just approximations. Nevertheless, such approximated results are reasonable ones as long as the external source would not be sufficient to drive the system towards criticality, because the keff values obtained were substantial low during the operation period. The initial keff obtained (keff(initial)=neutronsfission+neutronsspallation) were appropriate for an ADS and the keff obtained during the burnup (keff=neutronsfission) provide an insight of the necessary contribution of the external source on the operation period.(ii)
The insertion of reprocessed fuel in a232ThO2 system leads to harder neutron spectra due to the plutonium presence, facilitating the 233U production. This is result of the protactinium effect reduction by hardening of the neutron spectrum.(iii)
The reprocessed fuel evolution was very similar when that fuel was spiked with thorium or depleted uranium. However, the233U production appears with an advantage of the reprocessed fuel spiked with thorium. Furthermore, the decrease in the criticality was smaller in the system spiked with thorium that allows an extension in the burnup period.



---

*Source: 934105-2012-05-13.xml*