Science and Technology of Nuclear Installations
Volume 2016 (2016), Article ID 2612459, 11 pages
http://dx.doi.org/10.1155/2016/2612459

Research Article

Time-Dependent Neutronic Analysis of a Power-Flattened Gas Cooled Accelerator Driven System Fuelled with Thorium, Uranium, Plutonium, and Curium Dioxides TRISO Particles

Gizem Bakır,¹ Gamze Genç,² and Hüseyin Yapıcı²

¹Cumhuriyet Üniversitesi Teknoloji Fakültesi, 58140 Sivas, Turkey
²Erciyes Üniversitesi Mühendislik Fakültesi, 38039 Kayseri, Turkey

Received 14 July 2016; Accepted 3 August 2016

Academic Editor: Eugenijus Ušpuras

Copyright © 2016 Gizem Bakır 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.

Abstract.

This study presents the power flattening and time-dependent neutronic analysis of a conceptual helium gas cooled Accelerator Driven System (ADS) loaded with TRISO (tristructural-isotropic) fuel particles. Target material is lead-bismuth eutectic (LBE). ThO₂, UO₂, PuO₂, and CmO₂ TRISO particles are used as fuel. PuO₂ and CmO₂ fuels are extracted from PWR-MOX spent fuel. Subcritical core is radially divided into 10 equidistant subzones in order to flatten the power produced in the core. Tens of thousands of these TRISO fuel particles are embedded in the carbon matrix fuel pebbles as five different cases. The high-energy Monte Carlo code MCNPX 2.7 with the LA150 library is used for the neutronic calculations. Time-dependent burnup calculations are carried out for thermal fission power () of 1000 MW using the BURN card. The energy gain of the ADS is in the range of 99.98–148.64 at the beginning of a cycle. Furthermore, the peak-to-average fission power density ratio is obtained between 1.021 and 1.029 at the beginning of the cycle. These ratios show a good quasi-uniform power density for each case. Furthermore, up to 155.1 g ²³³U and 103.6 g ²³⁹Pu per day can be produced. The considered system has a high neutronic capability in terms of energy multiplication, fissile breeding, and spent fuel transmutation with thorium utilization.

1. Introduction

Currently, Light Water Reactors (LWRs) and Canada Deuterium Uranium (CANDU) reactors, which use uranium sources, generate the most of nuclear electricity. Commercial nuclear reactors produce highly radioactive materials as high-level wastes. They mainly contain transuranic isotopes (Np, Pu, Am, and Cm) and long-lived fission products. Most countries prefer to bury these wastes in concrete containers in the sea bed. On the other hand, transmutation of these wastes by driving high-energetic neutron and/or proton source is an improved approach. Lawrence [1] is the first to transmute thorium to ²³³U by releasing fast neutrons from a spallation target bombarded with high-energetic protons. The potentials of nuclear fuel transmutation in various fusion-fission hybrid reactors fuelled with several spent fuels extracted from conventional nuclear reactors are investigated in our previous studies (Yapıcı et al. [2, 3]).

Many studies on transmutation of nuclear waste and fissile fuel breeding in ADSs have been performed in recent years. Barros et al. investigated potentials of fissile breeding and transmutation for a lead cooled ADS loaded with thorium and reprocessed fuel [4, 5]. Their results show that the use of ThO₂ and reprocessed fuel combination enabled uranium (²³³U) production without the initial ²³³U enrichment. Vu and Kitada [6] propose a conceptual design of ADS utilizing thorium fuel and reprocessed fuel for transmutation of nuclear waste and production of energy. Abánades and Pérez-Navarro [7] examine the transmutation of nuclear waste in a gas cooled ADS moderated with graphite. TRISO is formed from nuclear wastes and they show that 95% of plutonium can be transmuted (except ²⁴²Pu). García et al. [8] evaluate the real number of pebbles fitting in a cylindrical ADS core in detail. The transmutation of MA in a lead-bismuth cooled ADS is analyzed by Takizuka et al. [9] and Tsujimoto et al. [10]. Takizuka et al. indicate that the waste transmutation of 250 kg/year is obtained by 80% plant factor. According to Tsujimoto et al., their ADS design productively transmutes and burns the MA when the effective neutron multiplication factor is 0.97. A uranium spallation target in an ADS loaded with MA is considered by Ismailov et al. [11]. Lead-bismuth (PbBi) target and uranium target are compared in their study and results show that the limited geometrical size of uranium target has better neutron multiplication than that of lead-bismuth (PbBi) target. The criticality of the targets made of pure ²⁴¹Am and ²⁴³Am is investigated by applying the Monte Carlo code for ADS [12]. The neutronic data of several infinite target mediums irradiated with a proton source of 1000 MeV is studied by Yapıcı et al. [13] to attain a real ADS design. Yapıcı et al. [14] and Bakir et al. [15] bring out that a high neutronic performance in nuclear waste transmutation, fissile breeding, and energy generation can be obtained for various considered configurations and fuel compositions of an ADS. Martínez et al. [16] claim that high-energetic neutrons are necessitated for transmutation of transuranium (TRU) and long-lived fission products in ADS. There are several other studies of various types of ADS for energy production and transmutation of radioactive wastes (Adam et al. [17]; Artisyuk et al. [18]; Brolly and Vértes, [19]; Haeck et al. [20]; Mukaiyama et al. [21]; Park et al. [22]; Seltborg and Wallenius [23]; Tsujimoto et al. [24]; Wade et al. [25]; Wallenius and Eriksson [26]; Westlén and Wallenius [27]).

The main cause of temperature and radiation damage gradients in an ADS is the nonuniform fission power density. Therefore, power flattening can help cooling of fuel core and reduce material stresses. To obtain a uniform fission power density profile, peak-to-average fission power density ratio () is expected to be reduced to 1.00. Power flattening of the helium-cooled DT driven blanket in the Prometheus-H (heavy ion) breeder reactor fuelled with various mixed fuels (UC–ThC, UO₂–ThO₂, UC–C, UO₂–C, and ²⁴⁴CmO₂–UO₂) and nuclear waste actinide are studied by Yapıcı and Übeyli [30]. According to their results, high amounts of energy, fissile fuel, and self-sufficient tritium production for (D, T) fusion reaction are obtained. Fissile fuel breeding and a flat fission power density are examined by Yapıcı [31] in a blanket-driven ICF (inertial confinement fusion) neutron and based on /SiC composite material. In order to accomplish a flat fission power density, the blanket is fuelled with ThO₂ and UO₂ mixed by several different mixing methods and cooled with natural lithium, (LiF)₂BeF₂, Li₁₇Pb₈₃, and ⁴He for the nuclear heat transfer. Peak-to-average fission power density ratio of the blanket is decreased to almost 1.1. Bakır et al. [32] have investigated the medical radioisotope production performance of a power-flattened ADS fuelled with UO₂ and PuO₂. Their results show that a good quasi-uniform power density is achieved and good neutronic performance in terms of energy production, radioisotope production, and the transmutation of spent fuel is also achieved.

In this study, power flattening and time-dependent neutronic analysis of a conceptual helium gas cooled ADS loaded with TRISO fuel particles are presented. The paper is organized as follows. In Section 2, the computational model of a conceptual gas cooled ADS is explained. Calculation procedure is outlined in Section 3. Numerical results and conclusions are presented in Sections 4 and 5, respectively.

2. Conceptual Gas Cooled Accelerator Driven System

Figure 1 demonstrates the layout of 1000 MeV linear accelerator (LINAC) for an ADS. As is apparent from this figure, a proton of 50 keV is accelerated to 1000 MeV in four steps: ion source (IS) of 50 keV is accelerated firstly to 3 MeV by radio frequency quadrupole (RFQ), secondly to 40 MeV by drift tube LINAC (DTL), thirdly to 100 MeV by cavity coupled DTL (CCDTL), and finally to 1000 MeV by superconducting linear accelerator (SC LINAC).

Figure 1: Layout of the 1 GeV LINAC for ADS (IS: ion source; RFQ: radio frequency quadrupol; DTL: drift tube LINAC; CCDTL: cavity coupled DTL; SC LINAC: superconducting linear accelerator).

Figures 2 and 3 show the cross-sectional view of the investigated conceptual cylindrical ADS and TRISO coated fuel particles in the fuel pebble sphere, respectively. There are four parts in the considered ADS as follows: spallation neutron target (SNT), subcritical core (SC), reflector zone (RZ), and shielding zone (SZ). In order to obtain flattened power profiles in SC, this core is radially divided into 10 equidistant subzones (see Figure 2). The isotopic percentages and densities of the materials used in this ADS are given in Table 1. Geometric characteristics of a TRISO particle (Kim et al. [29]) are given in Table 2.

Table 1: Isotopic percentages and densities of the materials used in the investigated ADS.

Material
Density [g/]
Nuclide
Percentage

LBE
11.344
Pb
44.5

9.8
Bi
55.5

He
0.1786
He
100

ThO₂
9.88
²³²Th
100

UO₂
10.54
²³⁵U
7.5, 15, 22.5

²³⁸U
92.5, 85, 77.5

PuO₂
11.50
²³⁸Pu
3.53535

²³⁹Pu
45.0154

²⁴⁰Pu
26.3505

²⁴¹Pu
15.9640

²⁴²Pu
9.13483

CmO₂
10.55
²⁴²Cm
3.91520

²⁴³Cm
0.04721

²⁴⁴Cm
85.5422

²⁴⁵Cm
9.54125

²⁴⁶Cm
0.95413

Graphite
2.10
¹²C
100

B₄C
2.52
¹⁰B
18.431

¹¹B
81.569

Discharged PWR-MOX fuel with plutonium recycle, 1000 MWe reactor, 80% capacity factor, 33 MWd/kg, 2.5% thermal efficiency, and 150 days after discharge [28, p. 370, Table 8.5].

Table 2: Characteristics of a TRISO particle [29].

Layer
Size [m]
Material
Density [g/]

Fuel kernel
250 (radius)
UO₂
10.55

PyC buffer
95 (thickness)
C/C
1.05

Inner PyC
40 (thickness)
C/C
1.90

SiC layer
35 (thickness)
SiC
3.18

Outer PyC
40 (thickness)
C/C
1.90

Figure 2: Cross-sectional view of the investigated ADS (① SNT: spallation neutron target; ② SCZ: Subcritical Core Zone; ③ RZ: reflector zone; ④ SZ: shielding zone; dimensions are in cm).

Figure 3: Fuel pebble sphere and TRISO coated fuel particle (PyC: pyrocarbon; SiC: silicon carbide).

Spallation Neutron Target. The target is liquid lead- (Pb-) bismuth (Bi) eutectic, (LBE: 44.5% Pb-55.5% Bi eutectic). Due to its good neutron release characteristic and thermal and chemical properties, the LBE is the most preferable target material for ADS designs among many target materials in the literature [33]. The radius of the target is optimized as 55 cm in terms of spallation neutrons. A continuous uniform proton beam of 1000 MeV and a source radius of 4 cm bombards the target to release a few tens of high-energy spallation neutrons. In this work, the target radius is gradually increased to optimize neutron leakage. It is determined as 55 cm and neutron leakage is calculated as about 30 neutrons per proton. The release neutrons penetrate through SC to activate fission and breeding reactions.

Subcritical Core. Tristructural-isotropic (TRISO) fuel is originated from a kernel made of fuel such as uranium carbide/uranium oxide and surrounded by carbon and ceramic layer. The TRISO fuel is an attractive fuel for high temperature nuclear reactors due to the fact that it has a quite high neutronic performance and good burn-up ability in high temperatures. As is apparent from Figure 3 and Table 2, a microspherical tristructural-isotropic (TRISO) fuel particle is designed as five layers from the inside out as follows: () fuel kernel, () Porous Carbon Buffer, () Inner Pyrolytic Carbon (IPyC), () silicon carbide (SiC), and () Outer Pyrolytic (OPyC) (see Figure 3) (Kim et al. [29]). Dimensions of these layers are given in Table 2. These TRISO composite particles are embedded inside a spherical carbon matrix fuel pebble with a determined packing fraction which can be up to 32% (Conway and Sloane [34]). Spherical fuel pebble is made of carbon matrix having a 60 mm diameter (see Figure 3). These carbon matrix fuel pebbles are located in a cylindrical SC with a determined packing fraction which can be up to 74% (Conway and Sloane [34]). In this study, to produce fissile fuel and energy and to transmute PWR-MOX spent fuel, tens of thousands of six different TRISO fuel particles (ThO₂, average 9% enrichment UO₂, 15% enriched UO₂, 22.5% enriched UO₂, PuO₂, and CmO₂) are embedded in the carbon matrix fuel pebbles with a packing factor of 29% as five different cases as follows and the embedding percentages of different TRISO fuel particles in the carbon matrix fuel pebbles in each subzone are given in Table 3.

Table 3: Enrichment percentages of UO₂ TRISO particles (Case 1(a)) and embedding percentages of 15% enriched UO₂, 22.5% enriched UO₂, PuO₂, and CmO₂ TRISO particles in Cases 1(b), 1(c), 2, and 3 in the fuel pebbles in the fuel subzones at the beginning of cycle [ MW].

Case
Number of the subzones

1
2
3
4
5
6
7
8
9
10
1–10

1(a)
8.52
9.32
9.82
10.25
10.75
11.25
11.25
9.92
6.47
3.20
7.5

1(b)
70.00
78.50
83.00
86.50
89.00
93.00
94.50
93.00
63.50
28.00
71.00

1(c)
47.00
56.00
60.50
63.50
66.50
68.50
69.50
66.00
43.00
18.25
46.00

2
20.30
26.00
27.90
28.85
29.76
31.40
33.65
35.55
26.80
6.76
26.50

3
7.90
7.80
7.80
7.90
8.10
8.31
8.31
7.42
6.01
4.04
6.45

Constant enrichment percentage or constant embedding percentages in all subzones.

Case 1. (a) Only UO₂ TRISO particles with different enrichment percentages (3.20%–11.25%) are embedded in the carbon matrix fuel pebbles and tens of thousands of these fuel pebbles are located in the ten equidistant layers according to their enrichment percentages.
(b) ThO₂ and 15% enriched UO₂ TRISO particles are separately embedded in the carbon matrix fuel pebbles with different embedding percentages and tens of thousands of these fuel pebbles are located in the ten equidistant layers according to their embedding percentages.
(c) ThO₂ and 22.5% enriched UO₂ TRISO particles are separately embedded in the carbon matrix fuel pebbles with different embedding percentages and tens of thousands of these fuel pebbles are located in the ten equidistant layers according to their embedding percentages.

Case 2. ThO₂ and PuO₂ TRISO particles are separately embedded in the carbon matrix fuel pebbles with different embedding percentages and tens of thousands of these fuel pebbles are located in the ten equidistant layers according to their embedding percentages.

Case 3. ThO₂ and CmO₂ TRISO particles are separately embedded in the carbon matrix fuel pebbles with different embedding percentages and tens of thousands of these fuel pebbles are located in the ten equidistant layers according to their embedding percentages.

Packing factor of the carbon matrix fuel pebbles in SC is 60%. This core is cooled with helium gas whose percentage is 40%. PuO₂ and CmO₂ fuels mentioned above are extracted from PWR-MOX spent fuels (Manson et al. [28]), fuel with plutonium recycle, 1000 MWe reactor, 80% capacity factor, 33 MWd/kg, 32.5% thermal efficiency, and 150 days after discharge.

Neutron Safety. The last two zones (RZ and SZ) serve for neutron safety. Reflector part is the zone made of graphite reflecting neutrons leaking from the SC zone to increase fission and breeding reactions. The graphite is selected because its scatter cross section is much greater than its absorption cross section. Additionally, it is a good neutron moderator and is high temperature resistant. Because of these properties, the graphite is widely preferable in nuclear applications. Shielding part is the last zone made of boron carbide (B₄C) absorbing the neutrons leaking from RZ. Due to the fact that boron has a quite high absorption cross section and B₄C has good thermomechanical properties, B₄C is usually used in nuclear reactors.

3. Calculation Procedure

The numerical calculations have been computed by using the high-energy Monte Carlo code MCNPX 2.7 (Pelowitz et al. [35]) with the LA150 library (Chadwick et al. [36]). “The library consists of evaluated reaction cross-sections and emission spectra up to 150 MeV for incident neutrons and protons, for over 40 target isotopes important in SNTs, structural materials, and shielding” (Yapıcı et al. [14]). Model for the intranuclear cascade of spallation reactions is selected as Bertini INC model (Bertini, [37]). The literature and our previous studies (Yapıcı et al. [13, 14]) bring out that the energy gain () is at the maximum level at proton energy () of 1000 MeV. Therefore, source proton energy is assumed as 1000 MeV (see Figure 2). The BURN card of the MCNPX 2.7 code (Pelowitz et al. [35]) is used for time-dependent calculations and these calculations are performed for the power-flattened cases.

4. Numerical Results

4.1. Effective Neutron Multiplication Coefficient

The effective neutron multiplication coefficient is the ratio of one generation of neutrons to the next generation. This coefficient is less than 1 in the subcritical systems. The compositions of the fuels in this study are determined in all cases so that is 0.98 at the beginning of the burn cycle. Thermal fission power is considered as 100, 300, 500, 700, and 1000 MW for the burn-up calculations. Figures 4–6 demonstrate the variations of with the operation time (a) and the burn cycle times calculated according to the decrease of from 0.95 to 90 (b) for all fuel cases.

(a)

(b)

Figure 4: (a) Variations of the effective neutron multiplication in the various cases of thermal power depending on the operation time. (b) The operation time at which effective neutron multiplication reached various values (0.90–0.95) versus the thermal power.

(a)

(b)

Figure 5: (a) Variations of the effective neutron multiplication in the various cases of thermal power depending on the operation time. (b) The operation time at which effective neutron multiplication reached various values (0.90–0.95) versus the thermal power.

(a)

(b)

Figure 6: (a) Variations of the effective neutron multiplication in the various cases of thermal power depending on the operation time. (b) The operation time at which effective neutron multiplication reached various values (0.90–0.95) versus the thermal power.

As is apparent from (a) of Figures 4–6, the profiles of gradually decline with the operation time in all ’s and fuel cases. These profiles indicate that fissile isotopes are consumed by burning with the operation time. For all values of (0.90–0.95), the operation times shorten exponentially with the increase of (see (b) of Figures 4–6). For  MW, the times that decreases from 0.98 to 0.95 are 0.75, 1.25, and 1.25 days in Cases 1(a), 1(b), and 1(c) (including natural, 15% and 22.5% enriched UO₂ TRISO particles), respectively. These times are 30 days and half a day in Cases 2 and 3 (PuO₂ and CmO₂ spent fuel TRISO particles). To prevent energy gain () from decreasing significantly, the fuel pebbles are refreshed when the value of decreases from 0.98 to 0.95.

4.2. Fission Power Density

Large numbers of energetic spallation and fission neutrons are released in SNT and SC of an ADS, and in turn, more fission reactions occur in SC. From inner side to outer side of SC in radial direction, profile of the fission power density declines exponentially because of the decrease of neutron fluxes in this direction. This nonuniformity of fission power density causes radiation damage and large temperature gradients that result in cooling and material issues. This situation is usually observed in ADSs and similar nuclear power systems. On the contrary, a uniform fission power density would help preventing of these cooling and material damage problems. A measure of uniformity of fission power density profile is determined with the peak-to-average fission power density ratio. This ratio must be 1.00 for an exact uniform fission power density profile.

In this work, in order to achieve a uniform power profile in SC, this core is radially divided into 10 equidistant subzones (see Figure 2), and the embedding percentages of different TRISO fuel particles in the carbon matrix fuel pebbles in each subzone of SC are adjusted in radial direction. The embedding percentages are optimized after many computational trials of percentages (see Table 3). The fission power density profiles obtained in the cases without percentage adjustments (constant percentages in all subzones) and in the cases with percentage adjustments (varying percentages in all subzones) are plotted in Figures 7–9 for all fuel cases. Dashed lines in these figures indicate the cases without percentage adjustments. In general, uniformity of power profiles of the cases without percentage adjustments breaks down (the profiles sharply increase) in the last subzones because of neutrons reflecting from RZ. On the other hand, in the optimized adjustment cases, quasi-uniform power profiles are achieved for all fuel cases. At the beginning of cycle (BOC), the values of are reduced to between 1.021 and 1.029 in Cases 1(a)–1(c), and as to Cases 2 and 3, these values are 1.029 and 1.024, respectively. It is important to note that these values are lower with respect to the values obtained from Yapıcı and Übeyli [30] ( 1.051–1.069) and Yapıcı [31] ( 1.131–1.403). The uniformity of power profile does not deteriorate significantly until decreases to 0.95. These values and power profiles bring out that the power of the investigated ADS is flattened for all fuel cases.

Figure 7: Variations of fission densities in the fuel core at the beginning and at the end of the cycle versus the core radius (solid lines indicate the flattened power case).

Figure 8: Variations of fission densities in the fuel core at the beginning and at the end of the cycle versus the core radius (solid lines indicate the flattened power case).

Figure 9: Variations of fission densities in the fuel core at the beginning and at the end of the cycle versus the core radius (solid lines indicate the flattened power case).

4.3. Fuel Burn-up

Burn-up (BU) is defined as the total energy generated per initial unit mass of initially loaded fuel, and its unit is GWd/MTU (GWd: gigawatt days; MTU: per metric ton of uranium). So it is directly related to the fission reactions and calculated as follows:

where is operation time and is time interval.

The calculations emphasize that the values of BU increase quasi-linearly in all ’s with operation time. In order to find an empirical relation for BU, curve-fitting method is used. Curve-fitting method gives the relation between independent variables and dependent variable by applying the best fit to data. The coefficients of the BU profiles obtained by this method are given in Table 4 for all fuel cases. The operation time limits in this BU() expression are between the BOC and the end of cycle (EOC) which varies depending on the fuel cases (see footnote of Table 5). The longest EOC time is in Case 2 and it is 30 days for . In this operation time, the BU value is 8.58 GWd/MTU for  MW.

Table 4: Increase of burn-up with the thermal power and time.

⁢ [GWd/MTU]

⁢ = 100–1000 MW

Case
1(a)
1(b)
1(c)
2
3

α
2.79
2.87
2.83
2.86
2.97

Table 5: Neutronic data at the beginning and at the end of the cycle.

= 1000 MW
Case

1(a)
1(b)
1(c)
2
3

Fission
BOC
692.17
743.24
675.64
615.68
499.92

EOC
263.56¹
234.39³
248.73⁵
179.48⁷
78.38⁹

EOC
109.81²
106.88⁴
107.48⁶
94.01⁸
47.34¹⁰

BOC
138.43
148.64
135.13
123.14
99.98

EOC
52.71
46.88
49.75
35.90
15.68

EOC
21.96
21.38
21.50
18.80
9.47

PE [MW]
BOC
7.22
6.73
7.40
8.12
10.00

EOC
18.97
21.33
20.10
27.85
63.77

EOC
45.54
46.77
46.51
53.19
105.60

PF × 10¹⁷ [p/s]
BOC
0.451
0.420
0.462
0.507
0.624

EOC
1.184
1.331
1.255
1.738
3.981

EOC
2.843
2.919
2.903
3.320
6.592

⁢, , and = 0.90

¹After 0.75 days. ²After 25 days.

³After 1.25 days. ⁴After 30 days.

⁵After 1.25 days. ⁶After 30 days.

⁷After 30 days. ⁸After 165 days.

⁹After 0.5 days. ¹⁰After 2.5 days.

4.4. Gain

The energy gain proportional to fission reactions () occurring in SC is one of the most major parameters of an ADS. It can be defined as the ratio of the total fission energy production in the core to E p and calculated as follows:

or

where is the energy per fission (200 MeV) and PE is the proton beam power.

The gain values of the investigated ADS, calculated for  MW at the BOC and the EOC, are given in Table 5. While these values vary in the range of 99.98–148.64 at the BOC depending on the fuel case, this range is 15.68–52.71 at the EOC (). Case 1(b) (including 15% enriched UO₂ TRISO particles) among Cases 1(a)–1(c) is the best fuel case in terms of the energy gain. As for the spent fuel cases, Case 2 (including PuO₂ TRISO particles) is the best. In order to maintain  MW, the value of proton flux (PF) is varied in the range of protons/s at the BOC depending on the fuel case; and during the operation, PF is increased to the range of 1.184 to 3.981 × 10¹⁷ protons/s. 10¹⁷ protons of 1000 MeV of energy per second correspond to a PE of 16.02 MW.

4.5. Thorium Utilization

As known, thorium utilization means production of fissile ²³³U. This production reaction is as follows:

Masses of isotopes which are loaded to SC at the BOC and produced in this core at the EOC are given in Tables 6–8 for all fuel cases. Although the highest thorium load (3114 kgHM) is in Case 3 (including CmO₂ TRISO particles), the best thorium utilization is realized in Case 2. As for the cases with uranium, Case 1(c) (including 22.5% enriched UO₂ TRISO particles) is the best case among Cases 1(a)–1(c) due to the fact that it has more enrichment uranium. Fissile fuel ²³³U of 155.1 g per day (4.653 kg/30 days) and 6.1 g per day ( kg/1.25 days) is produced by utilizing thorium for in Cases 2 and 1(c), respectively.

Table 6: Masses of actinide isotopes [kgHM] at the beginning and at the end of the cycle (Cases 1(a)–1(c)).

Case
1(a)
1(b)
1(c)

Isotopes
BOC
EOC
EOC
BOC
EOC
EOC
BOC
EOC
EOC

²³¹Th
—
—
—
—
7.448E − 5
1.597E − 4
—
1.680E − 4
2.989E − 4

²³²Th
—
—
—
8.015E + 2
8.012E + 2
7.941E + 2
1.519E + 3
1.519E + 3
1.507E + 3

²³³Th
—
—
—
—
5.405E − 3
5.611E − 3
—
9.076E − 3
9.442E − 3

²³³Pa
—
—
—
—
2.861E − 1
4.904E + 0
—
4.844E − 1
8.354E + 0

²³³U
—
—
—
—
4.488E − 3
2.002E + 0
—
7.619E − 3
3.456E + 0

²³⁴U
—
7.634E − 6
8.355E − 4
—
4.174E − 4
2.331E − 1
—
6.433E − 4
3.748E − 1

²³⁵U
3.120E + 2
3.110E + 2
2.807E + 2
4.051E + 2
4.034E + 2
3.663E + 2
4.372E + 2
4.355E + 2
3.978E + 2

²³⁶U
—
2.194E − 1
6.952E + 0
—
3.892E − 1
8.978E + 0
—
3.954E − 1
9.181E + 0

²³⁷U
—
3.845E − 4
5.459E − 2
—
5.932E − 4
6.547E − 2
—
4.674E − 4
6.445E − 2

²³⁸U
3.276E + 3
3.275E + 3
3.252E + 3
2.325E + 3
2.324E + 3
2.304E + 3
1.525E + 3
1.524E + 3
1.509E + 3

²³⁹U
—
2.151E − 2
2.274E − 2
—
1.574E − 2
1.675E − 2
—
1.175E − 2
1.240E − 2

²³⁸Pu
—
—
1.053E − 3
—
—
1.629E − 3
—
—
1.494E − 3

²³⁹Pu
—
6.525E − 2
1.768E + 1
—
1.295E − 1
1.615E + 1
—
9.697E − 2
1.224E + 1

²⁴⁰Pu
—
7.266E − 4
9.220E − 1
—
1.359E − 3
7.479E − 1
—
9.769E − 4
5.428E − 1

²⁴¹Pu
—
—
1.733E − 1
—
—
1.435E − 1
—
—
1.008E − 1

Total
3.588E + 3
3.586E + 3
3.559E + 3
3.532E + 3
3.529E + 3
3.498E + 3
3.481E + 3
3.480E + 3
3.448E + 3

After 0.75 days (). ²After 25 days ().

³After 1.25 days (). ⁴After 30 days ().

⁵After 1.25 days (). ⁶After 30 days ().

Table 7: Masses of actinide isotopes [kgHM] at the beginning and at the end of the cycle (Case 2).

Isotopes
BOC
EOC
EOC

²³¹Th
—
6.711E − 4
6.574E − 4

²³²Th
2.468E + 3
2.451E + 3
2.374E + 3

²³³Th
—
1.231E − 2
1.297E − 2

²³³Pa
—
1.100E + 1
1.995E + 1

²³³U
—
4.653E + 0
5.359E + 1

²³⁴U
—
4.039E − 1
5.973E + 0

²³⁵U
—
8.541E − 3
6.315E − 1

²³⁸Pu
3.632E + 1
3.520E + 1
3.045E + 1

²³⁹Pu
4.644E + 2
4.360E + 2
3.256E + 2

²⁴⁰Pu
2.730E + 2
2.696E + 2
2.516E + 2

²⁴¹Pu
1.661E + 2
1.648E + 2
1.535E + 2

²⁴²Pu
9.543E + 1
9.525E + 1
9.465E + 1

²⁴³Pu
—
3.420E − 2
3.485E − 2

²⁴⁴Pu
—
7.213E − 4
4.449E − 3

²⁴¹Am
—
6.296E − 1
2.790E + 0

²⁴²Am
—
1.944E − 3
3.276E − 2

²⁴³Am
—
3.129E + 0
1.345E + 1

²⁴⁴Am
—
4.735E − 3
1.914E − 2

Total
3.503E + 3
3.472E + 3
3.326E + 3

After 30 days ().

²After 165 days ().

Table 8: Masses of actinide isotopes [kgHM] at the beginning and at the end of the cycle (Case 3).

Isotopes
BOC
EOC
EOC

²³¹Th
—
4.644E − 4
1.414E − 3

²³²Th
3.114E + 3
3.113E + 3
3.110E + 3

²³³Th
—
3.370E − 2
3.608E − 2

²³³Pa
—
7.145E − 1
3.730E + 0

²³³U
—
4.390E − 3
1.170E − 1

²³⁴U
—
7.117E − 4
2.635E − 2

²³⁵U
—
—
1.304E − 4

²⁴²Cm
9.990E + 0
9.961E + 0
9.848E + 0

²⁴³Cm
1.210E − 1
1.258E − 1
1.453E − 1

²⁴⁴Cm
2.201E + 2
2.198E + 2
2.187E + 2

²⁴⁵Cm
2.465E + 1
2.428E + 1
2.284E + 1

²⁴⁶Cm
2.475E + 0
2.553E + 0
2.863E + 0

Total
3.371E + 3
3.370E + 3
3.368E + 3

After 0.5 days (). ²After 2.5 days ().

In addition to these, fissile fuel ²³⁹Pu of 103.6 g per day ( kg/1.25 days) and 77.6 g per day ( kg/1.25 days) is produced by utilizing uranium for in Cases 1(b) and 1(c), respectively. Reaction of ²³⁹Pu produced from ²³⁸U is as follows:

5. Conclusion

One way of transmutation of nuclear wastes is using ADSs by driving high-energetic proton source. In this study, a conceptual LBE target-ADS has been investigated to transmute spent fuel with thorium and to flatten the fission power density of the fuel core, as well as energy production. Many time-dependent computational trials have been performed to obtain optimum neutronic data by using MCNPX 2.7 computer code. The major results are, in brief, as follows:

(i)The optimized radius (55 cm) of LBE target is the best size in terms of spallation neutrons (about 30 neutrons).(ii)The longest and shortest times when decreases to 0.95 are in Case 2 (including ThO₂ and PuO₂ TRISO particles) and Case 3 (including the ThO₂ and CmO₂ TRISO particles), which are 0.5 days and 30 days, respectively. After these times, the fuel pebbles should be refreshed.(iii)A quasi-uniform fission power density profile is obtained for each fuel case ( 1.021–1.029 at the BOC).(iv)Burn-up profiles increase quasi-linearly in all ’s with operation time. The highest value is 8.58 GWd/MTU.(v)Case 1(b) (including ThO₂ and 15% enriched UO₂ TRISO particles) is the best fuel case in terms of the energy gain ().(vi)²³³U and ²³⁹Pu fissile fuels can be produced up to 155.1 g and 103.6 g per day, respectively.

In conclusion, the ADS flattened fission power density has a good neutronic performance in terms of effective utilizations of thorium and spent fuel and energy production. Furthermore, this work would light the way for obtaining a uniform fission power density profile for similar future works.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

This study is supported by the Research Fund of the Erciyes University, Project no. FDK-2015-5811.

References

E. O. Lawrence, “AEC research and development report, facilities for electronuclear (MTA) program,” Report LWS-24736, 1953.

H. Yapıcı, N. Demir, and G. Genç, “Neutronic analysis for transmutation of minor actinides and long-lived fission products in a Fusion-Driven Transmuter (FDT),” Journal of Fusion Energy, vol. 25, no. 3-4, pp. 225–239, 2006.

H. Yapıcı, G. Genç, and N. Demir, “Transmutation-incineration potential of transuraniums discharged from PWR-UO₂ spent fuel in modified PROMETHEUS fusion reactor,” Fusion Engineering and Design, vol. 81, no. 18, pp. 2093–2108, 2006.

G. P. Barros, C. Pereira, M. A. F. Veloso, and A. L. Costa, “Thorium and reprocessed fuel utilization in an accelerator-driven system,” Annals of Nuclear Energy, vol. 80, pp. 14–20, 2015.

G. D. P. Barros, C. Pereira, M. A. F. Veloso, and A. L. Costa, “Study of an ADS loaded with thorium and reprocessed fuel,” Science and Technology of Nuclear Installations, vol. 2012, Article ID 934105, 12 pages, 2012.

T. M. Vu and T. Kitada, “Transmutation strategy using thorium-reprocessed fuel ADS for future reactors in Vietnam,” Science and Technology of Nuclear Installations, vol. 2013, Article ID 674638, 5 pages, 2013.

A. Abánades and A. Pérez-Navarro, “Engineering design studies for the transmutation of nuclear wastes with a gas-cooled pebble-bed ADS,” Nuclear Engineering and Design, vol. 237, no. 3, pp. 325–333, 2007.

L. García, J. Pérez, C. García, A. Escrivá, J. Rosales, and A. Abánades, “Calculation of the packing fraction in a pebble-bed ADS and redesigning of the Transmutation Advanced Device for Sustainable Energy Applications (TADSEA),” Nuclear Engineering and Design, vol. 253, pp. 142–152, 2012.

T. Takizuka, K. Tsujimoto, T. Sasa, K. Nishihara, and H. Takano, “Design study of lead-bismuth cooled ADS dedicated to nuclear waste transmutation,” Progress in Nuclear Energy, vol. 40, no. 3-4, pp. 505–512, 2002.

K. Tsujimoto, T. Sasa, K. Nishihara, H. Oigawa, and H. Takano, “Neutronics design for lead-bismuth cooled accelerator-driven system for transmutation of minor actinide,” Journal of Nuclear Science and Technology, vol. 41, no. 1, pp. 21–36, 2004.

K. Ismailov, M. Saito, H. Sagara, and K. Nishihara, “Feasibility of uranium spallation target in accelerator-driven system,” Progress in Nuclear Energy, vol. 53, no. 7, pp. 925–929, 2011.

Y. Malyshkin, I. Pshenichnov, I. Mishustin, and W. Greiner, “Monte Carlo modeling of spallation targets containing uranium and americium,” Nuclear Instruments and Methods in Physics Research, Section B: Beam Interactions with Materials and Atoms, vol. 334, pp. 8–17, 2014.

H. Yapıcı, G. Genç, and N. Demïr, “Neutronic limits in infinite target mediums driven by high energetic protons,” Annals of Nuclear Energy, vol. 34, no. 5, pp. 374–384, 2007.

H. Yapıcı, G. Genç, and N. Demir, “A comprehensive study on neutronics of a lead-bismuth eutectic cooled accelerator-driven sub-critical system for long-lived fission product transmutation,” Annals of Nuclear Energy, vol. 35, no. 7, pp. 1264–1273, 2008.

G. Bakır, B. S. Selçuklu, G. Genç, and H. Yapıcı, “Neutronic analysis of LBE-uranium spallation target accelerator driven system loaded with uranium dioxide in TRISO particles,” Acta Physica Polonica A, vol. 129, no. 1, pp. 30–32, 2016.

A. H. Martínez, Y. Kadi, and G. Parks, “Transmutation of nuclear waste in accelerator-driven systems: thermal spectrum,” Annals of Nuclear Energy, vol. 34, no. 7, pp. 550–563, 2007.

J. Adam, K. Katovský, A. Balabekyan et al., “Transmutation of ¹²⁹I, ²³⁷Np, ²³⁸Pu, ²³⁹Pu, and ²⁴¹Am using neutrons produced in target-blanket system ‘Energy plus Transmutation’ by relativistic protons,” Pramana, vol. 68, no. 2, pp. 201–212, 2007.

V. Artisyuk, A. Chmelev, M. Saito, M. Suzuki, and Y. Fujii-E, “²⁴⁴Cm transmutation in accelerator-driven system,” Journal of Nuclear Science and Technology, vol. 36, no. 12, pp. 1135–1140, 1999.

Á. Brolly and P. Vértes, “Concept of a small-scale electron accelerator driven system for nuclear waste transmutation: part 2. Investigation of burnup,” Annals of Nuclear Energy, vol. 32, no. 4, pp. 417–433, 2005.

W. Haeck, E. Malambu, V. P. Sobolev, and H. Aït Abderrahim, “Assessment of americium and curium transmutation in magnesia based targets in different spectral zones of an experimental accelerator driven system,” Journal of Nuclear Materials, vol. 352, no. 1–3, pp. 285–290, 2006.

T. Mukaiyama, H. Takano, T. Ogawa, T. Takizuka, and M. Mizumoto, “Partitioning and transmutation studies at JAERI both under OMEGA program and high-intensity proton accelerator project,” Progress in Nuclear Energy, vol. 40, no. 3-4, pp. 403–413, 2002.

W. S. Park, T. Y. Song, B. O. Lee, and C. K. Park, “A preliminary design study for the HYPER system,” Nuclear Engineering and Design, vol. 219, no. 3, pp. 207–223, 2003.

P. Seltborg and J. Wallenius, “Proton source efficiency for heterogeneous distribution of actinides in the core of an accelerator-driven system,” Nuclear Science and Engineering, vol. 154, no. 2, pp. 202–214, 2006.

K. Tsujimoto, T. Sasa, K. Nishihara, T. Takizuka, and H. Takano, “Accelerator-driven system for transmutation of high-level waste,” Progress in Nuclear Energy, vol. 37, no. 1–4, pp. 339–344, 2000.

D. C. Wade, W. S. Yang, and H. Khalil, “ATW neutronics design studies,” Progress in Nuclear Energy, vol. 40, no. 3-4, pp. 497–504, 2002.

J. Wallenius and M. Eriksson, “Neutronics of minor-actinide burning accelerator-driven systems with ceramic fuel,” Nuclear Technology, vol. 152, no. 3, pp. 367–381, 2005.

D. Westlén and J. Wallenius, “Neutronic and safety aspects of a gas-cooled subcritical core for minor actinide transmutation,” Nuclear Technology, vol. 154, no. 1, pp. 41–51, 2006.

B. Manson, H. P. Thomas, and W. L. Hans, Nuclear Chemical Engineering, McGraw-Hill, New York, NY, USA, 1981.

H.-C. Kim, S. Y. Kim, J. K. Kim, and J. M. Noh, “Monte Carlo benchmark calculations for 400MWth PBMR core,” in Proceedings of the 13th International Conference on Emerging Nuclear Energy Systems (ICENES '07), pp. 498–502, Istanbul, Turkey, June 2007.

H. Yapıcı and M. Übeyli, “Power flattening in Prometheus breeder reactor using nuclear fuel and waste actinide,” Annals of Nuclear Energy, vol. 30, no. 2, pp. 159–173, 2003.

H. Yapıcı, “Power flattening of an inertial fusion energy breeder with mixed ThO₂–UO₂ fuel,” Fusion Engineering and Design, vol. 65, no. 1, pp. 89–108, 2003.

G. Bakır, S. B. Selçuklu, and H. Yapıcı, “Medical radioisotope production in a power-flattened ADS fuelled with uranium and plutonium dioxides,” Science and Technology of Nuclear Installations, vol. 2016, Article ID 5302176, 11 pages, 2016.

S. T. Mongelli, J. R. Maiorino, S. Anéfalos, A. Deppman, and T. Carluccio, “Spallation physics and the ADS target design,” Brazilian Journal of Physics, vol. 35, no. 3, pp. 894–897, 2005.

J. H. Conway and N. J. A. Sloane, Sphere-Packings, Lattices and Groups, Springer, New York, NY, USA, 3rd edition, 1999.

D. B. Pelowitz, “MCNPX, user’s manual, version 2.6.0,” Tech. Rep. LA-CP-07-1473, Los Alamos National Laboratory, 2008.

M. B. Chadwick, P. G. Young, S. Chiba et al., “Cross-section evaluations to 150 MeV for accelerator-driven systems and implementation in MCNPX,” Nuclear Science and Engineering, vol. 131, no. 2-3, pp. 293–328, 1999.

H. W. Bertini, “Low-energy intranuclear cascade calculation,” Physical Review, vol. 131, no. 4, p. 1801, 1963.