Abstract

This study presents the medical radioisotope production performance of a conceptual accelerator driven system (ADS). Lead-bismuth eutectic (LBE) is selected as target material. The subcritical fuel core is conceptually divided into ten equidistant subzones. The ceramic (natural U, Pu)O2 fuel mixture and the materials used for radioisotope production (copper, gold, cobalt, holmium, rhenium, thulium, mercury, palladium, thallium, molybdenum, and yttrium) are separately prepared as cylindrical rods cladded with carbon/carbon composite (C/C) and these rods are located in the subzones. In order to obtain the flattened power density, percentages of PuO2 in the mixture of UO2 and PuO2 in the subzones are adjusted in radial direction of the fuel zone. Time-dependent calculations are performed at 1000 MW thermal fission power () for one hour using the BURN card. The neutronic results show that the investigated ADS has a high neutronic capability, in terms of medical radioisotope productions, spent fuel transmutation and energy multiplication. Moreover, a good quasiuniform power density is achieved in each material case. The peak-to-average fission power density ratio is in the range of 1.02–1.28.

1. Introduction

Medical radioisotopes are used for the diagnosis and treatment of several illnesses. Radioisotopes can be produced in cyclotrons, nuclear reactors, or radioisotope generators. Accelerator driven systems (ADSs) can also be used for production of radioisotopes as well as fissile fuel breeding and energy generation. ADSs operate in subcritical mode; therefore, they are safer than conventional reactors which operate in critical mode. Large numbers of high energetic spallation neutrons can be produced in ADSs by means of target bombarding. Furthermore, these neutrons multiply in subcritical fuel core via fission reactions. For example, in Europe, MYRRHA (Multipurpose Hybrid Research Reactor for High-Tech Applications) is currently being studied. MYRRHA can produce neutron-rich radioisotopes due to its neutron flux characteristics [1]. ADSs can also be used for transmutation of nuclear wastes. Transmutation of the 99Tc, 129I, and 135Cs isotopes is studied in our previous work [2]. This study investigates production potentials of some medical radioisotopes in a conceptual cylindrical lead-bismuth eutectic (LBE) accelerator driven system (ADS) via neutron-gamma reactions.

Various methods of medical radioisotope production are studied by many researchers. Starovoitova et al. [3] investigate medical radioisotope production with photoneutron and photoproton reactions in linear accelerators. They consider 100Mo(γ,n)99Mo and 68Zn(γ,p)67Cu photonuclear capture reactions and compare results of Monte Carlo simulations with experimental data obtained with an electron accelerator. Richards et al. [4] examine the production of from 100Mo2C targets in a medical cyclotron accelerator. Webster et al. [5] study the production of 89Y(p,n)89Zr, 64Ni(p,n)64Cu, and 103Rh(p,n)103Pd reactions in a compact and low energy accelerator system. They suggest that the considered accelerator system has sufficient production quantities for the medical applications and this system can be used for production of many other medical isotopes. Abbas et al. [6] develop an accelerator driven neutron activator based on a modified version of the Adiabatic Resonance Crossing (ARC) concept for medical radioisotope production. Kin et al. [7] study new production routes for copper isotopes (64Zn(n,p)64Cu, 67Zn(n,p)67Cu, and 68Zn(n,x)67Cu) by using accelerator neutrons and claim that 64Zn(n,p)64Cu reaction is an encouraging route to produce 64Cu. Tárkányi et al. [811] examine the production of the therapeutic radioisotopes 165Er, 169Yb, and . Tárkányi et al. [12] later report new cross sections for 197Au(d,xn) and 197Au(d,x) nuclear reactions and discuss the production of the medically relevant isotopes 198Au and . They also provide the comparison with other charged particle induced production routes and the possible use of the 197Au(d,x) and reactions for monitoring deuteron beam parameters. Liem et al. [13] propose a conceptual design of a homogeneous solution reactor for production of medical radioisotope. Lebedev et al. [14] consider possible ways for production of 173Lu and usages of other radioactive nuclides as radiation sources for crack detectors in gamma defectoscopes. Artun and Aytekin [15] calculate the excitation functions for production of medical radioisotopes with proton, alpha, and deuteron induced reactions. Nortier et al. [16] measure integral excitation functions for the production of 16 radioisotopes of Cs, Xe, and I in the bombardment of with protons up to 100 MeV. Okuducu et al. [17] compute the nuclear level density parameters of some deformed radioisotopes of target nuclei (W, Hg) utilized on an ADS.

Power flattening can help in cooling of fuel core and reduce material stresses in an ADS because the nonuniform fission power density is the main cause of temperature and radiation damage gradients. Yapıcı and Übeyli [18] examine power flattening of the DT driven blanket in the Prometheus-H (heavy ion) breeder reactor cooled with helium and fuelled with different mixed fuels (UC–ThC, UO2–ThO2, UC–C, UO2–C, and 244CmO2–UO2) and nuclear waste actinide. Their calculations show that the breeder reactor has high neutronic performance and can produce significant amount of energy, fissile fuel, and tritium required for (D, T) fusion reaction. Yapıcı [19] also investigates the neutronic performance of a blanket-driven ICF (inertial confinement fusion) neutron and based on SiCf/SiC composite material for fissile fuel breeding and a flat fission power density. The blanket is loaded with ThO2 and UO2 mixed by various mixing methods to achieve a flat fission power density and cooled with natural lithium, (LiF)2BeF2, Li17Pb83, and 4He for the nuclear heat transfer. Peek-to-average fission power density ratio (г) of the blanket is reduced to ~1.1 which is expected to be reduced to ~1.00 to obtain a uniform fission power density profile. Yapıcı [20] considers the transmutation of transuranium (TRU) discharged from PWR spent fuel and the possibility of a flat fission power (FFP) generation along the transmutation process in the force-free helical reactor (FFHR). In this study, potential of a conceptual cylindrical proton accelerator driven system (ADS) fuelled with mixture of ceramic (U, Pu)O2 is investigated for medical radioisotope production as well as production of fissile fuel and energy. The paper is organized as follows. In Section 2, the computational model of a conceptual accelerator driven system is explained. Calculation procedure is outlined in Section 3. Numerical results and conclusions are presented in Sections 4 and 5, respectively.

2. Computational Model of Conceptual Accelerator Driven System

The geometric model of the considered ADS is plotted in Figure 1(a). Linear accelerator (LINAC) is used as proton accelerator. The proton acceleration process consists of a 50 keV Ion Source (IS), a 3 MeV Radio Frequency Quadrupole (RFQ), a 40 MeV Drift Tube LINAC (DTL), a 100 MeV Cavity Coupled DTL (CCDTL), and superconducting linear accelerator (SC LINAC) to accelerate the beam to 1 GeV. As is apparent from Figure 1(a), the considered ADS contains four different zones: (i) spallation neutron target (SNT), (ii) subcritical core (SC), (iii) reflector zone (RZ), and (iv) shielding zone (SZ). The densities of used materials in the investigated ADS are given in Table 1. Furthermore, the power of the considered ADS is flattened by varying PuO2 percentage in the mixture throughout the subcritical core.

(i) Spallation Neutron Target. This zone includes liquid lead- (Pb-) bismuth (Bi) eutectic (LBE: 44.5% Pb-55.5% Bi eutectic). Although there are many materials in literature as target material, the LBE is still the most attractive for ADS designs because of its good neutronic, chemical, and thermal properties. The SNT is bombarded with high energetic proton particles, which in turn releases a few tens of high-energy spallation neutrons depending on the energy of proton. These neutrons diffuse through the SC to make transmutation and breeding reactions. In LBE targets, 210Po, which is an alpha emitter, is produced. The accelerator components must be isolated to prevent contamination of 210Po and other possible radioactive wastes [21].

(ii) Subcritical Core. The ceramic (natural U, Pu)O2 fuel mixture and the materials used for radioisotope production (copper, gold, cobalt, holmium, rhenium, thulium, mercury, palladium, thallium, molybdenum, and yttrium) are separately prepared as cylindrical rods cladded with carbon/carbon composite (see Table 1). The considered PuO2 is extracted from PWR-MOX spent fuel (Manson et al. [22], fuel with plutonium recycle, 1000 MWe reactor, 80% capacity factor, 33 MWd/kg, and 32.5% thermal efficiency, 150 days after discharge). In order to achieve a flattened core power, the subcritical core is radially divided into ten equidistant subzones (these subzones are shown in Figure 1(a) as dashed lines) and the percentages of PuO2 in the fuel mixture are varied in each subzone. Nonetheless, the volumetric fractions (VFs) are kept the same. The cylindrical rods containing either the fuel or the materials are separately placed in the subzones of the subcritical fuel core in hexagonal order as shown in Figure 1(b) (for every 12 rods, 5 of them are fuel rods and 7 of them are material rods used for radioisotope production). Light water is used as the core coolant and it also serves as neutron moderator. According to the dimensions and arrangements of the rods, the VFs of the fuel mixture, the materials used for radioisotope production, the clad, and the coolant are calculated as 25%, 35%, 10%, and 30%, respectively.

(iii) Reflector Zone. The function of this zone is to reflect the neutrons escaping from the SC zone to enhance transmutation reactions. Therefore, graphite (carbon) is selected as a reflective material due to the fact that its scatter cross section is much greater than its absorption cross section. Moreover, the graphite is a good neutron moderator and it also has a high-temperature-resistance property. The graphite is widely used in nuclear implementation as an effective neutron reflector and moderator.

(iv) Shielding Zone. The role of this last zone made of boron carbide (B4C), which has a very high absorption cross section and excellent thermomechanical properties, is to absorb the neutrons leaking from the RZ. B4C is usually preferred in nuclear reactors as a neutron absorber.

3. Calculation Procedure

The neutronic computations have been carried out with the high-energy Monte Carlo code MCNPX 2.7 [23] by using the LA150 library [24]. “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 the SNTs, structural materials, and shielding” [2]. Bertini INC model [25] is used for the intranuclear cascade of spallation reactions. Time-dependent calculations have been performed at 1000 MW thermal fission power () for one hour using the BURN card option of the MCNPX2.7 code [23]. The burn cycle times are considered as one hour in all cases due to the fact that the produced medical isotopes have short half-lives.

In the literature and our previous studies [2, 26], it is found that the gain () reaches the maximum value when proton energy () is 1000 MeV. Hence, energy of one source proton is assumed as 1000 MeV in this study (see Figure 1(a)). A continuous uniform proton source bombards on the target material and the source radius is 4 cm.

4. Numerical Results

4.1. Neutron Flux

The spatial variations of neutron fluxes in the case of copper and in the case of copper and thulium are plotted in Figure 2. Generally, it is observed that the fluxes of neutrons, which are produced with spallation and fission reactions, decrease by deeper penetration in the target and subcritical zones. In particular, the number of neutrons having energy of less than 2–4 MeV decreases due to collisions with hydrogen and other atoms in the subcritical zones.

4.2. Flattened Fission Power Density

The fuel zone of an ADS is a region containing highly energetic spallation and fission neutrons. The fission power density in this zone decreases exponentially from inner side to outer side in radial direction due to the decrease of neutron fluxes. Therefore, the fission power density in an ADS is nonuniform in radial direction. The nonuniform fission power density is the main cause of temperature and radiation damage gradients. This characteristic is generally observed in fast ADS and similar power systems. On the other hand, a flattened fission power density would help in cooling of fuel core and reduce material stresses. The peak-to-average fission power density ratio (г) is a measure of fission power density uniformity. To obtain a uniform fission power density profile, this ratio must be reduced to ~1.00. In this study, in order to obtain the flattened power density, percentages of PuO2 in the mixture of UO2 and PuO2 in the subzones are adjusted in radial direction of the fuel zone (see Table 2). All values of г in all investigated cases are also given in Table 2. Values of г vary in the range of 1.02–1.28. The best uniformity is acquired in the case of thallium (г = 1.02). These values when compared to the values of Yapıcı and Übeyli [18] (г = 1.051–1.069) and Yapıcı [19] (г = 1.131–1.403) show that a good quasiuniform power density is achieved in each material case of this study.

Figure 3 shows variations of nonuniform and quasiuniform fission densities in the fuel core in the cases of mercury, palladium, thallium, molybdenum, and yttrium. As it is apparent from this figure, in the cases of constant PuO2 fraction in the fuel mixture, the fission power profiles decrease exponentially. However, these profiles increase toward the outer zone due to neutron reflection from the reflector zone. The profile of flattened fission power is lower than that of nonuniform fission power for all considered material cases and the highest flattened fission power profile is in the yttrium case (about 6.2·10−5 fissions/cm3).

In the cases of other considered materials (copper, gold, cobalt, holmium, rhenium, and thulium), the variations of nonuniform and quasiuniform fission densities in the fuel core are plotted in Figure 4. The first subplot of this figure shows the case where the copper rods are placed in all ten subzones. In other cases, the rods containing copper that has low capture cross section are placed into the first five subzones, and the rods including one of the other considered materials (gold, cobalt, holmium, rhenium, and thulium) having higher capture cross section than copper are located into the next five subzones. The reason that gold, cobalt, holmium, rhenium, and thulium are used with copper is to keep in the range of 0.97–0.98. Furthermore, for the same reason in the cases copper-rhenium and copper-thulium, the VFs of the fuel mixture and the materials have to be changed from 25% and 35% to 30% and 30%, respectively. Similar to the cases shown in Figure 3, in the cases of constant PuO2 fraction in the fuel mixture, the fission power profiles decrease rapidly and the profile of flattened fission power is lower than that of nonuniform fission power for all considered material cases. The highest flattened fission power profile is in the copper-holmium case (about 7.4·10−5 fissions/cm3).

4.3. Gain

The energy gain, , is the ratio of the total fission energy production in the fuel core to the energy of the proton beam and it is calculated as follows:where is the number of fission reactions and is the energy per fission (200 MeV). The gain can also be calculated aswhere is the thermal power and it is assumed as 1000 MW, and PE is the proton beam power.

The energy gain is one of the most important outputs of an ADS and the gain values acquired in this study are given in Table 2. These results show that gain values of power-flattened case (30.13–72.91) are lower than the values of the case of fuel mixed with a constant PuO2 fraction (65.43–115.04). The highest gain values of both cases are reached when copper-holmium is used:The value of proton flux (PF) varies in the ranges of 0.5·1017–0.96·1017 protons/s in the constant PuO2 fraction case and 0.86·1017–1.5·1017 protons/s in the flattened power case, depending on the PE value. 1017 protons (having a 1000 MeV of energy) per second correspond to a 16.02 MW of PE.

4.4. Medical Radioisotope Production

Medical radioisotopes are used for diagnosis, treatment, and therapy. Usage areas and production reactions of the radioisotopes considered in this study are shown in Table 3. Medical radioisotopes are generally produced artificially because most of the naturally found radioisotopes have long half-lives and they are mostly harmful for human body. Radioisotopes are usually produced in cyclotrons, nuclear reactors, or radioisotope generators depending on the target nucleus or energy of radiator particles. The conceptual system of this study is an accelerator driven system (ADS). A general radioisotope production reaction that takes place in an ADS is a neutron capture reaction as follows:Activities of the radioisotopes of this study at the end of one hour for the cases of fuel mixed with a constant PuO2 fraction and the cases of flattened power are presented in Tables 4 and 5, respectively. Activities of produced radioisotopes are generally higher in the flattened power cases except for 203Hg and 103Pd. The production of 198Au, 60Co, 166Ho, 186Re, and 170Tm increases almost threefold in the flattened power cases compared to the cases of fuel mixed with a constant PuO2 fraction. The reason of this increase is the adjustment of the PuO2 fraction in the fuel mixture; that is, the PuO2 fraction is decreased in the inner zones and increased in the outer zones in the flattened power cases (see Table 2). The production of other radioisotopes is almost the same in both cases. 197Au production (28.780 g) is the highest and the lowest production is 103Pd (0.192 g) among all the considered radioisotopes. These values are quite high with respect to medical isotope production in low energy proton accelerators (75 MeV–150 MeV) because the ADSs are bombarded with high energetic protons (PE > 500 MeV) and have a subcritical fuel core that also multiply neutrons.

5. Conclusions

A conceptual cylindrical ADS is examined for medical radioisotope production, spent fuel transmutation, and energy production. The main results of this study are given briefly as follows:(i)A good quasiuniform power density is achieved in each material case (г is in the range 1.02–1.28). The best uniform power density profile is obtained in the thallium case.(ii)The production of radioisotopes is in the range of 0.192–31.810 g at the end of one-hour cycle. The case of copper-thulium is the best case [31.810 g 170Tm ( Ci)].(iii)The highest flattened fission power profile is in the copper-holmium case (about 7.4·10−5 fissions/cm3).(iv)The energy gain is in the range of 30.13–72.91. The best energy multiplication is obtained in the case of copper-holmium.

In conclusion, the investigated ADS has a good neutronic performance in terms of energy production, radioisotope production, spent fuel transmutation, and management of nuclear waste.

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 Erciyes University, Project no. FDK-2015-5811.