# Medical Radioisotope Production in a Power-Flattened ADS Fuelled with Uranium and Plutonium Dioxides

Authors: Gizem Bakır; Saltuk Buğra Selçuklu; Hüseyin Yapıcı
Journal: Science and Technology of Nuclear Installations
 (2016)

---

## 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 (P
t
h) 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.


---

## Body

## 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 T
c
99
m 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. [8–11] examine the production of the therapeutic radioisotopes 165Er, 169Yb, and H
o
166
g. Tárkányi et al. [12] later report new cross sections for 197Au(d,xn)H
g
197
m
,
197
g
,
195
m
,
195
g
,
193
m
,
193
g and 197Au(d,x)A
u
198
m
,
198
g
,
196
m
,
196
g
,
195,194 nuclear reactions and discuss the production of the medically relevant isotopes 198Au and H
g
195
m
,
195
g
,
197
m
,
197
g. They also provide the comparison with other charged particle induced production routes and the possible use of the 197Au(d,x)H
g
197
m
,
197
g
,
195
m
,
193
m and A
u
196
m
,
196
g reactions for monitoring deuteron beam parameters. Liem et al. [13] propose a conceptual design of a homogeneous solution reactor for production of M
o
99
/
T
c
99
m 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 I
122–125 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 X
e
n
a
t 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 Figure1(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.Table 1
Isotopic fractions and densities of the materials used in the investigated ADS.
Material
Density [g/cm3]
Nuclide
Fraction [%]
LBE
11.344
Pb
44.5
9.80
Bi
55.5
H2O
1.00
—
100
C/C
1.85
12C
100
UO2
10.54
235U
0.7
238U
99.3
PuO2
11.50
238Pu
3.53535
239Pu
45.0154
240Pu
26.3505
241Pu
15.9640
242Pu
9.13483
For radioisotope production
Cu
8.92
63Cu
68.499
65Cu
31.501
Au
18.880
197Au
100
Co
8.900
59Co
100
Ho
8.795
165Ho
100
Re
20.530
185Re
37.148
187Re
62.852
Tm
9.321
169Tm
100
Hg
13.546
196Hg
14.6
198Hg
9.869
199Hg
16.763
200Hg
23.028
201Hg
13.225
202Hg
30.004
204Hg
6.965
Pd
12.02
102Pd
0.977
104Pd
10.877
105Pd
22.013
106Pd
27.199
109Pd
26.830
110Pd
12.104
Tl
11.850
203Tl
29.316
205Tl
70.684
Mo
10.20
92Mo
14.230
94Mo
9.053
95Mo
15.748
96Mo
16.673
97Mo
9.646
98Mo
24.623
100Mo
10.028
Y
4.469
89Y
100
Graphite
2.10
12C
100
B4C
2.52
10B
18.431
11B
81.569Figure 1
(a) Cross-sectional view of the investigated ADS (① SNT: spallation neutron target, ② SCZ: subcritical core zone (divided into ten equidistant subzones), ③ RZ: reflector zone, and ④ SZ: shielding zone; dimensions are in cm). (b) Top cross-sectional view of sample two hexagonal arrangements of the fuel mixture rods (shaded) and the material rods used for the radioisotope production in the subcritical core zone (inner and outer radiuses (r
i and r
o) of the rods are 0.51 cm and 0.55 cm, resp., pitch length P = 1.2 cm, and the volume fractions of fuel mixture (total 5 rods), materials used for isotope production (total 7 rods), clad, and coolant are 25%, 35%, 10%, and 30%, resp; the dimensions are not in scale).
(a)
(b)(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 (P
t
h) 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 (G) reaches the maximum value when proton energy (E
p) 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 Figure2. 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.Figure 2
Neutron fluxes [SCZ = subcritical zone, energy = (E
i
+
1
-
E
i)/Lethargy, and Lethargy = ln
⁡
(
E
i
+
1
/
E
i
)].
(a)
(b)


### 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.Table 2
Neutronic data in the cases of0.97
<
k
e
f
f
<
0.98.
Material
Percentages of PuO2 in the mixture of UO2 and PuO2 in the subzones
Neutronic data
1
2
3
4
5
6
7
8
9
10
1–10a
г
G
Hg
25.50
33.00
38.00
40.00
42.00
42.00
43.00
47.00
51.00
44.00
41.50
1.66a
1.06
78.48a
40.48
Pd
26.50
40.00
48.50
51.50
54.00
55.00
55.50
59.00
63.00
45.00
53.00
1.63
1.14
69.38
30.13
Tl
5.75
10.00
14.50
17.00
17.75
18.75
20.25
22.25
22.50
6.75
16.00
1.48
1.02
67.30
53.77
Mo
11.5
16.5
23.0
26.0
28.0
28.5
30.5
32.5
34.0
17.0
26.4
1.56
1.10
72.47
40.64
Y
2.2
3.7
6.0
8.1
10.1
11.9
12.9
13.5
12.3
3.0
2.2
1.36
1.03
78.26
62.00
Cu
6.00
10.00
15.00
17.00
18.00
19.50
21.50
23.50
25.50
10.50
17.50
1.59
1.03
65.43
41.48
Cu and Au
8.00
14.00
19.00
24.00
33.00
55.00
65.00
68.00
74.00
81.00
26.00
2.64
1.11
106.95
58.78
Cu and Co
6.50
10.70
16.40
22.10
27.80
37.80
42.40
46.50
50.00
47.50
24.50
2.41
1.11
103.09
45.98
Cu and Ho
9.00
15.00
20.50
25.50
38.00
76.50
90.00
95.00
99.00
97.50
27.00
2.77
1.16
115.04
72.91
Cu and Re
7.00
14.00
18.00
20.00
30.00
95.00
100.00
100.00
100.00
100.00
23.30
2.83
1.28
106.63
53.65
Cu and Tm
5.00
10.00
14.00
16.00
26.00
76.00
84.00
85.50
86.00
82.00
23.00
2.81
1.07
97.27
42.11
aIn the cases of fuel mixed with a constant PuO2 fraction.Figure3 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).Figure 3
Variations of fission densities in the fuel core versus the core radius (solid line indicates the flattened power case).
(a)
(b)
(c)
(d)
(e)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 Figure4. 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 k
e
f
f 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).Figure 4
Variations of fission densities in the fuel core versus the core radius (solid line indicates the flattened power case).
(a)
(b)
(c)
(d)
(e)
(f)


### 4.3. Gain

The energy gain,G, 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:1
a
G
=
R
f
E
f
E
p
,where R
f is the number of fission reactions and E
f is the energy per fission (200 MeV). The gain can also be calculated as1
b
G
=
P
t
h
P
E
,where P
t
h 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 Table2. 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:1
c
P
E
=
P
t
h
G
.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 Table3. 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:2
i
s
o
t
o
p
e
n
,
γ
A
i
s
o
t
o
p
e
A
+
1
.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.Table 3
Usage areas and reactions of the produced radioisotopes.
Radioisotope
Half-life
Usage areas
64Cu
63
C
u
+
n →
γ 64
C
u
12.7 hours
Radiopharmaceuticals for PET (positron emission tomography) imaging of cancer
198Au
197
A
u
+
n →
γ 198
A
u
2.69 days
Treatment of cancer of prostate and bladderfor reduction of fluid accumulation secondary to cancer
60Co
59
C
o
+
n →
γ 60
C
o
5.27 years
Modern radiation therapy and industrial radiography
166Ho
165
H
o
+
n →
γ 166
H
o
26.77 hours
Medical radiotherapeutic applications, radioembolization for the treatment of patients with liver metastases
170Tm
169
T
m
+
n →
γ 170
T
m
128.6 days
X-ray source for cancer therapy with brachytherapy
186Re
185
R
e
+
n →
γ 186
R
e
3.7183 days
Radiotherapy, bone pain palliation
197Hg
196
H
g
+
n →
γ 197
H
g
64.128 hours
Scanning kidneys
103Pd
102
P
d
+
n →
γ 103
P
d
16.991 days
Treatment of prostate cancer
204Tl
203
T
l
+
n →
γ 204
T
l
3.78 years
Myocardial perfusion scans
99Mo
98
M
o
+
n →
γ 99
M
o
65.94 hours
99
mTc generators
99
mTc
99
M
o → 99
m
T
c
+
β
-
+
ν
e
6 hours
Imaging of the skeleton, heart muscle, brain, thyroid, lungs, liver, spleen, kidney, gall bladder, bone marrow, salivary, lacrimal glands, heart blood pool, and infection
90Y
89
Y
+
n →
γ 90
Y
64.053 hours
Cancer brachytherapy and as silicate colloid for the relieving of the pain of arthritis in larger synovial joints.Table 4
Activities of the radioisotopes at the end of one hour in the case of fuel mixed with a constant PuO2 fraction (P
t
h
=
1000 MW).
Materials
Isotopes
Activity (A) [Ci]
Mass [g]
Cu
29064
8.229
E
+
06
2.132
Cu and Au
29064
5.868
E
+
06
1.520
79198
2.857
E
+
06
11.680
Cu and Co
29064
5.581
E
+
06
1.446
27060
3.127
E
+
03
2.764
Cu and Ho
29064
5.978
E
+
06
1.549
67166
6.782
E
+
06
9.624
Cu and Re
29064
5.188
E
+
06
1.344
75186
9.486
E
+
05
5.102
75188
6.167
E
+
06
6.281
Cu and Tm
29064
5.220
E
+
06
1.353
69170
6.530
E
+
04
10.930
Hg
80197
1.422
E
+
06
5.728
80203
1.190
E
+
04
0.862
Pd
46103
1.433
E
+
04
0.192
46107
1.183
E
-
03
2.300
46109
1.051
E
+
07
5.002
46111
2.030
E
+
07
0.280
46112
6.087
E
+
04
0.046
Tl
81202
9.066
E
+
02
0.017
81204
4.343
E
+
03
9.366
Mo
42099
1.075
E
+
06
2.237
Y
39090
3.402
E
+
05
0.626Table 5
Activities of the radioisotopes at the end of one hour in the flattened power case (P
t
h
=
1000 MW).
Materials
Isotopes
Activity [Ci]
Mass [g]
Cu
29064
8.574
E
+
06
2.222
Cu and Au
29064
2.909
E
+
06
0.754
79198
7.040
E
+
06
28.780
Cu and Co
29064
2.809
E
+
06
0.728
27060
6.662
E
+
03
5.887
Cu and Ho
29064
2.930
E
+
06
0.759
67166
1.809
E
+
07
25.680
Cu and Re
29064
2.656
E
+
06
0.688
75186
2.707
E
+
06
14.560
75188
1.821
E
+
07
18.550
Cu and Tm
29064
2.285
E
+
06
0.592
69170
1.901
E
+
05
31.810
Hg
80197
4.453
E
+
05
5.912
80203
1.551
E
+
03
0.850
Pd
46103
1.869
E
+
03
0.192
46107
1.562
E
-
04
2.298
46109
1.302
E
+
06
4.978
46111
2.610
E
+
06
0.279
46112
7.791
E
+
03
0.046
Tl
81202
8.069
E
+
02
0.015
81204
4.524
E
+
03
9.756
Mo
42099
1.092
E
+
06
2.272
Y
39090
4.296
E
+
05
0.790


## 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 Figure2. 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.Figure 2
Neutron fluxes [SCZ = subcritical zone, energy = (E
i
+
1
-
E
i)/Lethargy, and Lethargy = ln
⁡
(
E
i
+
1
/
E
i
)].
(a)
(b)


## 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.Table 2
Neutronic data in the cases of0.97
<
k
e
f
f
<
0.98.
Material
Percentages of PuO2 in the mixture of UO2 and PuO2 in the subzones
Neutronic data
1
2
3
4
5
6
7
8
9
10
1–10a
г
G
Hg
25.50
33.00
38.00
40.00
42.00
42.00
43.00
47.00
51.00
44.00
41.50
1.66a
1.06
78.48a
40.48
Pd
26.50
40.00
48.50
51.50
54.00
55.00
55.50
59.00
63.00
45.00
53.00
1.63
1.14
69.38
30.13
Tl
5.75
10.00
14.50
17.00
17.75
18.75
20.25
22.25
22.50
6.75
16.00
1.48
1.02
67.30
53.77
Mo
11.5
16.5
23.0
26.0
28.0
28.5
30.5
32.5
34.0
17.0
26.4
1.56
1.10
72.47
40.64
Y
2.2
3.7
6.0
8.1
10.1
11.9
12.9
13.5
12.3
3.0
2.2
1.36
1.03
78.26
62.00
Cu
6.00
10.00
15.00
17.00
18.00
19.50
21.50
23.50
25.50
10.50
17.50
1.59
1.03
65.43
41.48
Cu and Au
8.00
14.00
19.00
24.00
33.00
55.00
65.00
68.00
74.00
81.00
26.00
2.64
1.11
106.95
58.78
Cu and Co
6.50
10.70
16.40
22.10
27.80
37.80
42.40
46.50
50.00
47.50
24.50
2.41
1.11
103.09
45.98
Cu and Ho
9.00
15.00
20.50
25.50
38.00
76.50
90.00
95.00
99.00
97.50
27.00
2.77
1.16
115.04
72.91
Cu and Re
7.00
14.00
18.00
20.00
30.00
95.00
100.00
100.00
100.00
100.00
23.30
2.83
1.28
106.63
53.65
Cu and Tm
5.00
10.00
14.00
16.00
26.00
76.00
84.00
85.50
86.00
82.00
23.00
2.81
1.07
97.27
42.11
aIn the cases of fuel mixed with a constant PuO2 fraction.Figure3 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).Figure 3
Variations of fission densities in the fuel core versus the core radius (solid line indicates the flattened power case).
(a)
(b)
(c)
(d)
(e)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 Figure4. 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 k
e
f
f 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).Figure 4
Variations of fission densities in the fuel core versus the core radius (solid line indicates the flattened power case).
(a)
(b)
(c)
(d)
(e)
(f)


## 4.3. Gain

The energy gain,G, 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:1
a
G
=
R
f
E
f
E
p
,where R
f is the number of fission reactions and E
f is the energy per fission (200 MeV). The gain can also be calculated as1
b
G
=
P
t
h
P
E
,where P
t
h 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 Table2. 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:1
c
P
E
=
P
t
h
G
.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 Table3. 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:2
i
s
o
t
o
p
e
n
,
γ
A
i
s
o
t
o
p
e
A
+
1
.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.Table 3
Usage areas and reactions of the produced radioisotopes.
Radioisotope
Half-life
Usage areas
64Cu
63
C
u
+
n →
γ 64
C
u
12.7 hours
Radiopharmaceuticals for PET (positron emission tomography) imaging of cancer
198Au
197
A
u
+
n →
γ 198
A
u
2.69 days
Treatment of cancer of prostate and bladderfor reduction of fluid accumulation secondary to cancer
60Co
59
C
o
+
n →
γ 60
C
o
5.27 years
Modern radiation therapy and industrial radiography
166Ho
165
H
o
+
n →
γ 166
H
o
26.77 hours
Medical radiotherapeutic applications, radioembolization for the treatment of patients with liver metastases
170Tm
169
T
m
+
n →
γ 170
T
m
128.6 days
X-ray source for cancer therapy with brachytherapy
186Re
185
R
e
+
n →
γ 186
R
e
3.7183 days
Radiotherapy, bone pain palliation
197Hg
196
H
g
+
n →
γ 197
H
g
64.128 hours
Scanning kidneys
103Pd
102
P
d
+
n →
γ 103
P
d
16.991 days
Treatment of prostate cancer
204Tl
203
T
l
+
n →
γ 204
T
l
3.78 years
Myocardial perfusion scans
99Mo
98
M
o
+
n →
γ 99
M
o
65.94 hours
99
mTc generators
99
mTc
99
M
o → 99
m
T
c
+
β
-
+
ν
e
6 hours
Imaging of the skeleton, heart muscle, brain, thyroid, lungs, liver, spleen, kidney, gall bladder, bone marrow, salivary, lacrimal glands, heart blood pool, and infection
90Y
89
Y
+
n →
γ 90
Y
64.053 hours
Cancer brachytherapy and as silicate colloid for the relieving of the pain of arthritis in larger synovial joints.Table 4
Activities of the radioisotopes at the end of one hour in the case of fuel mixed with a constant PuO2 fraction (P
t
h
=
1000 MW).
Materials
Isotopes
Activity (A) [Ci]
Mass [g]
Cu
29064
8.229
E
+
06
2.132
Cu and Au
29064
5.868
E
+
06
1.520
79198
2.857
E
+
06
11.680
Cu and Co
29064
5.581
E
+
06
1.446
27060
3.127
E
+
03
2.764
Cu and Ho
29064
5.978
E
+
06
1.549
67166
6.782
E
+
06
9.624
Cu and Re
29064
5.188
E
+
06
1.344
75186
9.486
E
+
05
5.102
75188
6.167
E
+
06
6.281
Cu and Tm
29064
5.220
E
+
06
1.353
69170
6.530
E
+
04
10.930
Hg
80197
1.422
E
+
06
5.728
80203
1.190
E
+
04
0.862
Pd
46103
1.433
E
+
04
0.192
46107
1.183
E
-
03
2.300
46109
1.051
E
+
07
5.002
46111
2.030
E
+
07
0.280
46112
6.087
E
+
04
0.046
Tl
81202
9.066
E
+
02
0.017
81204
4.343
E
+
03
9.366
Mo
42099
1.075
E
+
06
2.237
Y
39090
3.402
E
+
05
0.626Table 5
Activities of the radioisotopes at the end of one hour in the flattened power case (P
t
h
=
1000 MW).
Materials
Isotopes
Activity [Ci]
Mass [g]
Cu
29064
8.574
E
+
06
2.222
Cu and Au
29064
2.909
E
+
06
0.754
79198
7.040
E
+
06
28.780
Cu and Co
29064
2.809
E
+
06
0.728
27060
6.662
E
+
03
5.887
Cu and Ho
29064
2.930
E
+
06
0.759
67166
1.809
E
+
07
25.680
Cu and Re
29064
2.656
E
+
06
0.688
75186
2.707
E
+
06
14.560
75188
1.821
E
+
07
18.550
Cu and Tm
29064
2.285
E
+
06
0.592
69170
1.901
E
+
05
31.810
Hg
80197
4.453
E
+
05
5.912
80203
1.551
E
+
03
0.850
Pd
46103
1.869
E
+
03
0.192
46107
1.562
E
-
04
2.298
46109
1.302
E
+
06
4.978
46111
2.610
E
+
06
0.279
46112
7.791
E
+
03
0.046
Tl
81202
8.069
E
+
02
0.015
81204
4.524
E
+
03
9.756
Mo
42099
1.092
E
+
06
2.272
Y
39090
4.296
E
+
05
0.790


## 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 g170Tm (1.901
E
+
05 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.



---

*Source: 5302176-2016-05-04.xml*