# The Use of Th in HTR: State of the Art and Implementation in Th/Pu Fuel Cycles

Authors: Guido Mazzini; Eleonora Bomboni; Nicola Cerullo; Emil Fridman; Guglielmo Lomonaco; Eugene Shwageraus
Journal: Science and Technology of Nuclear Installations
 (2009)

---

## Abstract

Nowadays nuclear is the only greenhouse-free source that can appreciably respond to the increasing worldwide energy demand. The use of Thorium in the nuclear energy production may offer some advantages to accomplish this task. Extensive R&D on the thorium fuel cycle has been conducted in many countries around the world. Starting from the current nuclear waste policy, the EU-PUMA project focuses on the potential benefits of using the HTR core as a Pu/MA transmuter. In this paper the following aspects have been analysed: (1) the state-of-the-art of the studies on the use of Th in different reactors, (2) the use of Th in HTRs, with a particular emphasis on Th-Pu fuel cycles, (3) an original assessment of Th-Pu fuel cycles in HTR. Some aspects related to Thorium exploitation were outlined, particularly its suitability for working in pebble-bed HTR in a Th-Pu fuel cycle. The influence of the Th/Pu weight fraction at BOC in a typical HTR pebble was analysed as far as the reactivity trend versus burn-up, the energy produced per Pu mass, and the Pu isotopic composition at EOC are concerned. Although deeper investigations need to be performed in order to draw final conclusions, it is possible to state that some optimized Th percentage in the initial Pu/Th fuel could be suggested on the basis of the aim we are trying to reach.


---

## Body

## 1. Introduction

Since the beginning of  the nuclear age, the  possibility of using Thorium as a nuclear fuel in thermal reactors appeared a very promising alternative to uranium exploitation. First of all, Thorium natural resources are three times more abundant than Uranium ones. What’s more, the lower mass number of the fertile element (Th232)  entails a by far smaller build-up of high mass number TRU, which are as known responsible for the most of the long-term radiotoxicity of the nuclear waste. On the other hand, natural Thorium is composed of a single isotope (Th232) that is not fissile, but only fertile. By a neutron capture, Th232 transmutes into U233, which is a fissile nuclide characterized by an excellent neutronic behaviour in the thermal range  [1, 2]. Hence, Th-based fuels need some amounts of U235 or Pu acting as a driver. Conversely, this kind of fertilization process produces, by (n,2n) reaction, some U232, which is a very pernicious nuclide because of its strong γ-emitting decay daughter. Clearly, it results particularly problematic if SNF reprocessing is envisaged.The High-Temperature gas-cooled Reactor (HTR or HTGR) is a graphite moderated, He cooled nuclear reactor that was studied and realised in the past, and that it has been recovered today in the framework of the Generation IV Initiative. HTR shows unique features due to its peculiar fuel form (TRISO-coated particles embedded in cylindrical or spherical fuel elements made of graphite) as well as to its inert coolant (He). The gaseous coolant allows us the possibility to make the neutronics and the thermalfluiddynamics of the core substantially independent of each other.  What’s more, He is an inert coolant from the chemical point of view as well, and it does not absorb neutrons parasitically. The TRISO-coated particles permit  us the possibility to reach burn-ups that are by far longer than those of LWRs (i.e., in principle up to 800 GWd/tHM without leakages of Fission Products).Starting from these consideration and taking into account the unique characteristics of the HTR concept, it is possible to appreciate the potentialities of loading the HTR core with a Th-based fuel.  First of all, reaching very long burn-ups without needing to reprocess the fuel in-between in principle enables us to exploit Th without the drawbacks previously mentioned (strongγ-rays of its daughter). Moreover, the excellent neutron economy of Th-fuels is coupled with the excellent neutron economy of HTR itself (indeed, He and graphite are the main components of this core). Nevertheless, a driver fuel is necessary to sustain the chain reaction at BOC.As far as the driver fuel is concerned, two main options are available, that is, enriched U or Pu. Although enriched U undoubtedly would be the best choice from the neutronic point of view, Gen. IV aims suggest Pu exploitation.In the present paper, on the basis of some work performed by our group in the past, Pu (RG) is adopted as a driver fuel instead of HEU. Particularly, a pebble-bed HTR using Thorium along with 1st generation Plutonium has been considered, analogously to the fuel considered in the framework of the EU PUMA project [3]. This work aims at analysing the influence of the initial Pu/Th mass ratio in the fresh fuel as far as the main neutronic and burn-up parameters are concerned: reactivity at BOC, possibility to reach very long irradiation cycles, Pu isotopic composition versus burn-up, energy released. Hence, the following aspects will be analysed:(i)
the state of art about the studies focusing on the use of Th in different core concepts,(ii)
the use of Th in HTRs, with a particular emphasis on Th-Pu fuel cycles,(iii)
an original assessment of Th-Pu fuel cycles in HTR, assuming, as evaluating parameters,
(a)
the total energy produced,(b)
the energy produced per initial loaded Pu,(c)
the ratio between final (discharged) and initial (loaded) Pu mass.


## 2. Thorium Resources, Drawbacks of the Current Fuel Cycle, and Proliferation Concerns: A Brief Overview

With  more than 300 NPPs running worldwide, LWR is currently the most widespread nuclear reactor technology. Developed and applied successfully since the 1950s, this technology is safe, reliable, and well proven [4]. Its main drawback consists of the extremely small U resource exploitation (less than 1%), which implies what follows:(i)
concerns about the limitation of U resource availability,(ii)
proliferation risk due to Pu content of SNF (please remember that the end of Cold War raised concerns about the proliferation risk, due to the large stockpiles of Pu produced by civil and military applications),(iii)
long-term toxicity of the final waste, mainly due to TRUs built up in SNF.Indeed, an integral U exploitation would entail not only a huge availability of nuclear fuel resources  but also a nuclear waste mainly composed of FPs instead of HMs.The commonly mentioned design objectives of Th-based fuel cycles include the following:(i)
increasing nuclear fuel resources by breeding U233  from thorium,(ii)
improving fuel utilization in thermal reactors,(iii)
reducing significantly the U235  enrichment requirements,(iv)
reducing the build-up of Pu  and of other TRU in comparison with uranium-based fuel cycles,(v)
burning Pu (RG or WP) without reprocessing, then reducing existing Pu stockpiles,(vi)
increasing the possibility to breed U233 and to incinerate long lived radiotoxic isotopes by combining the thorium-based fuel cycle with accelerator driven systems (ADS) and/or hybrid fusion driven systems.A Th-based fuel cycle for HTR is not a new idea: the THTR experience represents and summarizes the effort made in that direction. What is more, as shown by some past work, adopting the Th fuel cycle in HTR using Pu as a driver would allow us what follows:(i)
increase of the efficiency in TRU reduction, as illustrated in Figure1,(ii)
achieving higher fuel burn-ups.Figure 1
Mass of actinides at EOC starting from 1 g of Pu [4].As previously mentioned, Th is by far more abundant than U in the earth’s crust. The primary natural Th resource[5] is monazite, which is a thorium phosphate mineral. Thorium deposits have been found in several countries [5] as shown in Figure 2. According to the US Geological Survey [5], the world thorium reserves and reserve base (resources) are shown in Figure 3.Figure 2
Thorium Deposits in the World [5].Figure 3
The World Thorium Reserves and Reserve Base [5].Moreover, Th extraction technology as well as processes to recover valuable material from Th-based SNF (e.g., THOREX process) had been studied since a long time [2].Basic researches and developments  on the Th fuel cycle[1] have been conducted in many countries, and  a considerable experience about Th-based   fuels has been gained in power reactors worldwide (i.e., determination of material data, fabrication tests on the laboratory scale, irradiation experiments in material test reactors, use of Th-based fuel in HTRs, LWRs (including WWERs), LMFBRs, and (potentially) in MSRs, postirradiation examinations of spent Th-based fuels, fabrication of Th-based fuel, both on pilot and semi-industrial scale (AVR)). Further activities are on-going worldwide, aiming at using Th to burn MA in both Gen-IV reactors and/or ADS.Due to the lack of data, it seems impractical to develop meaningful cost projections[5] for any nuclear energy systems using thorium. Historical examples give some ideas of the funds that may be required. As an example, in the 1970s, Germany spent the equivalent of 500 million  Euros in current money to develop a thorium fuel cycle  and 2.5 billion  Euros for the development of HTR. More recently, the GIF in its technology roadmap estimated that the only assessment of viability and performance of a nuclear system requires around one billion dollars (before any decision to develop and build a demonstrator, which of course would require large additional funding). The market of Th raw material [5], that is, ores or concentrates of oxide of thorium, is extremely small and limited to not-energy employments (such as high-temperature ceramics, crucibles, catalysts, welding electrodes, and some specific alloys). Additionally, the use of thorium in most of these products has been continuously decreasing, due to the burden associated with its  natural radioactivity. Then, a nearly constant commercial price in the range 30—35 US $/kg[5] was published in the 1960s and 1970s, but this publication was suspended due to the progressively decreasing demand. Thus, nowadays a commercial Th price is not available. What is more, fluctuations of Th price have been minimized by a supply that by far exceeds the demand due to Th by-product build-up. That means that many countries have fairly large stockpiles of thorium that they consider as a waste.If Th were exploited to produce electricity in the future, a Th market would develop. On the basis of what previously explained, it is reasonable to suppose that the raw material contribution to the cost of the electricity generated will remain low, comparable with or lower than that of the uranium cycle. So, from the economical point of view, the idea of using Th as a fertile material will virtually have no significant cost.


## 3. A Brief Review of Thorium-Based Fuel Cycles for Different Core Concepts

### 3.1. Thorium-Plutonium Fuel in LWR

The possibility of using Th-based fuels for LWR has been extensively analysed since a long time. Performances of Th-U and Th-Pu mixtures have been studied and compared with other  more common concepts. Just to set an example, among some interesting comparisons between MOX and Pu-Th fuels shown in  [6], the case of Korean 900 MWe  PWR can be taken into account. Its infinite multiplication factors when fuelled by Th-based and MOX fuel assembly,  respectively, were calculated by means of the HELIOS code [7]. The boron concentration was kept constant (500 ppm). Figure 4 shows  the results obtained. This example is useful to appreciate the very good neutronic behaviour of Th-based fuels, but we have to underline that there are not any advantages concerning the long-term radiotoxicity of the final waste.Figure 4
k-infinite as a function of burn-up [6].Indeed, Cm244   contribution at the beginning  and  Pu239   and Th229  contribution successively make the radiotoxicity of the considered Th-based fuels higher than that of conventional UO2 PWR SNF.


### 3.2. Thorium-Uranium Fuel in CANDU Reactors

The possibility of using Th as fertile material in CANDU has been largely considered[8]. Particularly, three main types of fuel mixtures have been analysed: mixed SEU and ThO2 bundles, gadolinium-doped mixed bundles, and high burn-up ThO2  bundles. The variation of lattice  kinf  as well as the change in fissile content as a function of the bundle average burn-up is shown in Figure 5, for each of these considered compositions. Physical properties of natural UO2 and natural ThO2 fuels, respectively, have been shown for comparison purposes.Figure 5
kinf   as a function of burn-up in CANDU [8].Please note that the initial fissile content of the high burn-up thorium bundles has been carefully chosen so that the depletion rate of the fissile material is almost the same as the conversion rate of the fertile Th232 into fissile U233. Consequently, the reactivity and the fissile content of the high burn-up Th bundles are almost constant throughout the entire lifetime of the bundles.


### 3.3. Thorium-Uranium Fuel in Molten-Salt Reactors

Although the MSR core  appears very different from all the other nuclear reactor concepts, the graphite-moderate MSR actually owns some features that are substantially similar to those belonging to High-Temperature Gas-Cooled Reactor (HTGR). Indeed, in both concepts the graphite moderator is at an average temperature between600—700°C, a Th-U mixture can be adopted as a fuel, and the moderation ratios are very close to each other. Therefore, a large part of the thermal reactor technology (and, particularly, of HTGRs) is directly applicable to this kind of MSR design. Therefore, nuclear data and methods for calculating neutron fluxes, reactivity temperature coefficients, and so on  developed for HTGR can be applied for the MSR design, even though not-negligible differences remain. Some examples of applications of Th-based fuel in MSR could be found in [9, 10].


## 3.1. Thorium-Plutonium Fuel in LWR

The possibility of using Th-based fuels for LWR has been extensively analysed since a long time. Performances of Th-U and Th-Pu mixtures have been studied and compared with other  more common concepts. Just to set an example, among some interesting comparisons between MOX and Pu-Th fuels shown in  [6], the case of Korean 900 MWe  PWR can be taken into account. Its infinite multiplication factors when fuelled by Th-based and MOX fuel assembly,  respectively, were calculated by means of the HELIOS code [7]. The boron concentration was kept constant (500 ppm). Figure 4 shows  the results obtained. This example is useful to appreciate the very good neutronic behaviour of Th-based fuels, but we have to underline that there are not any advantages concerning the long-term radiotoxicity of the final waste.Figure 4
k-infinite as a function of burn-up [6].Indeed, Cm244   contribution at the beginning  and  Pu239   and Th229  contribution successively make the radiotoxicity of the considered Th-based fuels higher than that of conventional UO2 PWR SNF.


## 3.2. Thorium-Uranium Fuel in CANDU Reactors

The possibility of using Th as fertile material in CANDU has been largely considered[8]. Particularly, three main types of fuel mixtures have been analysed: mixed SEU and ThO2 bundles, gadolinium-doped mixed bundles, and high burn-up ThO2  bundles. The variation of lattice  kinf  as well as the change in fissile content as a function of the bundle average burn-up is shown in Figure 5, for each of these considered compositions. Physical properties of natural UO2 and natural ThO2 fuels, respectively, have been shown for comparison purposes.Figure 5
kinf   as a function of burn-up in CANDU [8].Please note that the initial fissile content of the high burn-up thorium bundles has been carefully chosen so that the depletion rate of the fissile material is almost the same as the conversion rate of the fertile Th232 into fissile U233. Consequently, the reactivity and the fissile content of the high burn-up Th bundles are almost constant throughout the entire lifetime of the bundles.


## 3.3. Thorium-Uranium Fuel in Molten-Salt Reactors

Although the MSR core  appears very different from all the other nuclear reactor concepts, the graphite-moderate MSR actually owns some features that are substantially similar to those belonging to High-Temperature Gas-Cooled Reactor (HTGR). Indeed, in both concepts the graphite moderator is at an average temperature between600—700°C, a Th-U mixture can be adopted as a fuel, and the moderation ratios are very close to each other. Therefore, a large part of the thermal reactor technology (and, particularly, of HTGRs) is directly applicable to this kind of MSR design. Therefore, nuclear data and methods for calculating neutron fluxes, reactivity temperature coefficients, and so on  developed for HTGR can be applied for the MSR design, even though not-negligible differences remain. Some examples of applications of Th-based fuel in MSR could be found in [9, 10].


## 4. The Use of Thorium in HTR

### 4.1. Preliminary Considerations on TH-PU Cycle in HTR

The reduction of Pu stockpiles as well as the discharge of a nuclear waste that is very poor from the Pu content point of view   is among the main goals of the current nuclear research. Then, Th-Pu-based fuels have to be considered as promising candidates to reach these objectives. What is more, the possibility of very long irradiation cycles allowed by HGTR cores has on the other hand to be taken into account as well. These considerations justify the efforts to develop Th-U and Th-Pu fuel cycles in HTGR core concepts [11–15]. Among the others, a few researches on Th-Pu fuel cycle in HTR were performed by our group in the past [4, 16].The results shown in  [4] lead us to draw some considerations concerning the behaviour of highly-burned (i.e., up to 800 GWd/tHM) Thorium-Plutonium fuel. Figure 6 reports the Level Of Mine Balancing Time (LOMBT) trend as a function of burn-up: it is clear that the best choice in terms of  radiotoxicity reduction is to reach a burn-up that is as high as possible. This results is mainly due to a progressively larger Pu consumption.Figure 6
LOMBT versus Burn-up (2/3 Th – 1/3 Pu) [4].Considering the rather long LOMBT of a typical LWR SNF, the reduction of the final radiotoxicity and of the total amount of waste is of course a key issue. Looking at this latter, the actinide mass reduction in the final waste to be stored in a permanent repository is another very attractive characteristic of HTRs (Figure7).Figure 7
Masses of actinides versus fuel discharge burn-up (2/3 Th – 1/3 Pu) [4].Despite a larger production of curium isotopes (especially Cm244), it is clearly shown by Figure 8 that reaching high burn-ups allows us to strongly reduce the masses of dangerous actinides.Figure 8
MCNP model without cut CPs [16].The main results of  [4] are summarized in Table 1.Table 1
Summary of obtained results(k∞~1.05) [4].
Fuel/PebbleEFPDLOMBTEnergy/PebblePuout/PuinTRN*/Puin1 g Pu 1st gen.6004300051.8 GJ32.9%37.6%1.5 g Pu 2nd  gen.6453500055.7 GJ46.4%54.9%1 g Pu 1st + 2 g Th7053400060.9 GJ26.6%32.8%NOTE:*TRN =  Puout + Am + Cm.


### 4.2. The PUMA Project

Recently, Th as fuel for HTR has been considered, among the others, in the framework of the EU PUMA project [5]. Starting from the current legacy of nuclear waste  as well as from the (V)HTR good capabilities of Pu burning, the PUMA activities focus  on core physics investigations, in order  to further  demonstrate the potential benefits of using the HTR core as a Pu/MA transmuter (as known, annually a 1000 MWe PWR produces about 30 tons of SNF (burn-up around 30 GWD/tHM) with the following average composition [7]:  94% U238, 1% U235, 1% Pu, 0.1% MA (Np, Am and Cm), 3—4% FP). Moreover, as (V)HTR Pu/MA transmuters are foreseen to operate in a global system of reactor designs and fuel cycle facilities, fuel cycle studies and socioeconomic/environmental assessments have been carried out. The PUMA consortium gathers 16 organisations from 9 countries (BE, DE, FR, IT, NL, PL, SE, UK, USA), including research organisations, leading nuclear engineering and fuel cycle firms, a fuel manufacturer, a utility, universities and institutes, and consultancy SMEs. The work carried out in the framework of PUMA represents part of the EURATOM contribution to the GIF [9].Among the other  activities of PUMA project, a deep burn-up fuel cycle strategy using the Th-Pu fuel in PBMR reactor [8] has been investigated.Many of the data used in the following sections are based on PUMA project findings.


## 4.1. Preliminary Considerations on TH-PU Cycle in HTR

The reduction of Pu stockpiles as well as the discharge of a nuclear waste that is very poor from the Pu content point of view   is among the main goals of the current nuclear research. Then, Th-Pu-based fuels have to be considered as promising candidates to reach these objectives. What is more, the possibility of very long irradiation cycles allowed by HGTR cores has on the other hand to be taken into account as well. These considerations justify the efforts to develop Th-U and Th-Pu fuel cycles in HTGR core concepts [11–15]. Among the others, a few researches on Th-Pu fuel cycle in HTR were performed by our group in the past [4, 16].The results shown in  [4] lead us to draw some considerations concerning the behaviour of highly-burned (i.e., up to 800 GWd/tHM) Thorium-Plutonium fuel. Figure 6 reports the Level Of Mine Balancing Time (LOMBT) trend as a function of burn-up: it is clear that the best choice in terms of  radiotoxicity reduction is to reach a burn-up that is as high as possible. This results is mainly due to a progressively larger Pu consumption.Figure 6
LOMBT versus Burn-up (2/3 Th – 1/3 Pu) [4].Considering the rather long LOMBT of a typical LWR SNF, the reduction of the final radiotoxicity and of the total amount of waste is of course a key issue. Looking at this latter, the actinide mass reduction in the final waste to be stored in a permanent repository is another very attractive characteristic of HTRs (Figure7).Figure 7
Masses of actinides versus fuel discharge burn-up (2/3 Th – 1/3 Pu) [4].Despite a larger production of curium isotopes (especially Cm244), it is clearly shown by Figure 8 that reaching high burn-ups allows us to strongly reduce the masses of dangerous actinides.Figure 8
MCNP model without cut CPs [16].The main results of  [4] are summarized in Table 1.Table 1
Summary of obtained results(k∞~1.05) [4].
Fuel/PebbleEFPDLOMBTEnergy/PebblePuout/PuinTRN*/Puin1 g Pu 1st gen.6004300051.8 GJ32.9%37.6%1.5 g Pu 2nd  gen.6453500055.7 GJ46.4%54.9%1 g Pu 1st + 2 g Th7053400060.9 GJ26.6%32.8%NOTE:*TRN =  Puout + Am + Cm.


## 4.2. The PUMA Project

Recently, Th as fuel for HTR has been considered, among the others, in the framework of the EU PUMA project [5]. Starting from the current legacy of nuclear waste  as well as from the (V)HTR good capabilities of Pu burning, the PUMA activities focus  on core physics investigations, in order  to further  demonstrate the potential benefits of using the HTR core as a Pu/MA transmuter (as known, annually a 1000 MWe PWR produces about 30 tons of SNF (burn-up around 30 GWD/tHM) with the following average composition [7]:  94% U238, 1% U235, 1% Pu, 0.1% MA (Np, Am and Cm), 3—4% FP). Moreover, as (V)HTR Pu/MA transmuters are foreseen to operate in a global system of reactor designs and fuel cycle facilities, fuel cycle studies and socioeconomic/environmental assessments have been carried out. The PUMA consortium gathers 16 organisations from 9 countries (BE, DE, FR, IT, NL, PL, SE, UK, USA), including research organisations, leading nuclear engineering and fuel cycle firms, a fuel manufacturer, a utility, universities and institutes, and consultancy SMEs. The work carried out in the framework of PUMA represents part of the EURATOM contribution to the GIF [9].Among the other  activities of PUMA project, a deep burn-up fuel cycle strategy using the Th-Pu fuel in PBMR reactor [8] has been investigated.Many of the data used in the following sections are based on PUMA project findings.


## 5. Calculations and Results

### 5.1. Calculation Parameters and Models

In this research, on the basis of the advantages outlined in previous paragraphs, we investigated the main neutronic characteristics related to the utilization of Th as fertile material in PBMR loaded with a Pu/Th fuel, by means of MCNP5  [17] code.The model used to perform these analyses is described in the following lines.In our simulations Pu-RG (i.e., Plutonium produced during irradiation of UO2 fuel in a typical LWR) has been chosen as driver fuel. In principle,  Pu-RG is not suitable for proliferation purposes, due to its high content of even isotopes (Pu238, Pu240, Pu242), which are characterized by a high spontaneous fission rate and, at least in the case of Pu238, by a high decay heat as well. Th content in the  considered pebble has been varied from 0% to 100% of the HM constituting the fuel. Table 2 shows the main parameters of the analyzed PBMR fuel elements, while Table 3 shows the Pu-RG (1st   gen. Pu) isotopic vector.Table 2
(Pu-Th)O2 fuel general characteristics in accordance with PUMA Project [10, 11].
DescriptionValueFuel pebbleFuel pebble outer radius3.0 cmThickness of fuel free zone0.5 cmTotal heavy metal loading per fuel pebble2.0 gMatrix density1.74 g/cm3Packing fraction in pebble bed61%Coated particleKernel Coating MaterialC/C/SiC/CLayer thickness90/40/35/40μmLayer densities1.05/1.90/3.18/1.90 g/cm3Fuel kernelFuel kernel diameter200 micronKernel material type(Pu-Th)O2.0Th-oxide density9.55 g/cm3 (95% TD) [13]Pu-oxide density10.89 g/cm3 (95% TD) [12]Table 3
Composition of 1st  generation Pu [10, 11].
Pu isotopesFirst Generation (weight %)2382.5923953.8524023.6624113.132426.78Additionally, we fixed the following calculation parameters.(i)
An infinite lattice of stochastically arranged pebble has been considered for each fuel composition, by modelling a single sphere (with the He gap associated) with white boundary conditions.(ii)
Each of the considered pebbles contains 2 g HM (Pu/Th) [5].(iii)
Due to the lack of data and for the sake of simplicity, the density of (Pu,Th)O2 fuel (with Th content varying from 0% to 100% of HM) has been assumed as varying linearly from the density of PuO2 kernels (10.89 g/cm3  [12]) to the density of  ThO2 kernels (9.55 g/cm3  [13]), as a function of the Th fraction.(iv)
Geometric error reduction: the pebbles have been modelled as actual as possible (Figure15), also on the basis of the results obtained in [16].(v)
No temperature profile inside the pebble: the temperature is uniform everywhere and equal to 1200 K.(vi)
The cross section libraries used in this work are JEFF3.1 [18].(vii)
Helium Gap: according to the average packing fraction of the PBMR core (61% [8]) the helium gap has been explicitly modelled.Figure8 shows an x-y view of the pebble model adopted to perform these calculations. Please note that clipped CPs along the boundary of the fuelled zone have been eliminated.


### 5.2. Preliminary Results

The variation of thek-infinite at BOC as a function of Th percentage is shown in Figure 9. Analysing the trend of this curve, we can note that some trends can be recognised on the basis of Th content:(i)
0% ≤ Th% < 10%: the k-infinite value remains near 1.26 (flat trend);(ii)
10% ≤ Th% < 30%: a slightly positive gradient is shown (the k-infinite increases from 1.26 to 1.27);(iii)
30% ≤ Th% < 84.5%: the value of k-infinite increases linearly from ~1.27 and to ~1.39;(iv)
84.5% ≤ Th% < 100%: the k-infinite decreases from the peak at 84.5% (k∞=1.39274, σ=0.00154) to 1.Figure 9
Variation ofk-infinite at BOC as a function of Th content.In order to complete this analysis, the fuel behaviour as a function of burn-up has to be considered. All the following burn-up calculations has been  performed by means of the BGCore code [19].


### 5.3. BGCore Code

The burn-up calculations were performed with BGCore. BGCore is a software package for comprehensive computer simulation of nuclear reactor systems and their fuel cycles. It consists of a number of modules:(i)
3-dimensional coupled neutronic and thermal hydraulic module, which calculates the reactor core criticality, as well as the core power and temperature distributions; only steady state core analysis is considered at this stage of the BG-Core package development; no thermal feedback features were used in this analysis;(ii)
fuel depletion and decay module, which calculates the fuel isotopic composition;(iii)
auxiliary modules for basic data processing and management.


### 5.4.k-Effective versus k-Infinite Burn-Up Calculations Comparison

The “real’’k-effective in PBMR depends  not only from the material and the core geometry, as in LWR reactors, but also from a random contribution due to pebbles recirculation [20]. Then, a not-negligible approximation is introduced by calculating the  k-effective fixing the pebble arrangement in the core. In addition, even adopting this approximation, k-effective calculations are quite expensive form the computational point of view.Hence, in order to roughly estimate the differences betweenk-effective  model (i.e., fixed geometry) and k-infinite one, we made a preliminary evaluation, performing two burn-up calculations with a representative fuel composition. Particularly, a fuel composed of 50% Th and 50% Pu was adopted. This choice was also due to the fact that the results  obtained can be directly compared with an independent calculation performed by FZJ [15].Figure10 and Table 4 show our main results.Table 4
Puout/Puin for  k-infinite and  k-effective models.
Pu238out/Pu238inPu239out/Pu239inPu240out/Pu240inPu241out/Pu241inPu242out/Pu242inPutotout/Putotink-infinite43,76%5,46%23,25%48,63%189,23%28,64%k-effective42,57%9,05%53,65%55,51%185,46%38,38%Figure 10
k-infinite and k-effective trends as a function of burn-up (expressed in MWd/tHM).Thus, we can outline some considerations.(i)
If we adoptk-inf = 1 as exit conditions for the burn-up calculations, we obtain a reasonable discharge burnup (comparable to that one we can obtain in a real pebble-bed HTR for a discharged pebble with the same initial composition).(ii)
The shape ofk trends versus burn-up are quite similar.(iii)
The ratios (OUT/IN) among Pu isotopes have (in the two cases) the same order of magnitude.(iv)
The total ratioPuout/Puin in the two cases considered is comparable.(v)
Finally the discharge burn-up obtained in the FZJ calculation is quite similar (Δ~5%) to that we obtained assuming k=1 as exit condition for the k-infinite calculation (see Table 5).Table 5
Comparison betweenk-infinite model and FZJ [15] results.
Heavy metal load and composition per pebbleFZJ calculation—discharge burn-up  [MWd/kgHM]k-infinite calculation—discharge burn-up  [MWd/kgHM)]1 g Pu, 1.0 g Th-232332.38314.42For the reasons described above and because of the kind of study (preliminary parametric evaluation), we performed the calculations shown in the following paragraphs using a simplifiedk-infinite model instead of a k-effective one.


### 5.5. Burn-Up Results

As anticipated above, we analysed some Th/Pu compositions by means of BGCore, in order to assess the fuel behaviour as a function of burn-up. The basic assumptions for these calculations are the following:(i)
times steps of 25 days at full power,(ii)
power per pebble fixed at 1 kW,(iii)
uniform fuel temperature (1200 K) inside the pebble,(iv)
discharge burn-up is set at the time whenk-infinite becomes equal to 1, to give a preliminary estimation of the theoretical nuclear burn-up.(As already underlined, please note that the real PBMR will not be fuelled with pebbles at the same burn-up level, but with a mixture of pebbles at different burn-ups [21]: the real discharge burn-up in principle could be extremely high [21], i.e., up to twice the assumed discharge burnup. On the other hand, we calculated the k-infinite instead of k-effective (so we have not took into account leakages and effects related to the structures surrounding the core)).The results concerningk-inf are summarized in Figure 11. Please note that all the standard deviations are around 100 pcm, which is the value commonly accepted in this kind of calculation by means of MC codes.Figure 11
k-infinite variation as a function of burn-up (expressed in MWd/tHM).As anticipated in the introduction, we evaluated the fuel behaviour on the basis of the following parameters:(i)
the total energy produced,(ii)
the energy produced per initial Pu mass,(iii)
the ratio between final (discharged) and initial (loaded) masses for the Pu isotopes.The obtained results are shown in  Figures12–19.Figure 12
Produced energy (from BOC up to the time at whichk-infinite becomes equal to unity) versus Th concentration.Figure 13
Energy produced pergPuin  versus Th concentration.Figure 14
Pu238out/Pu238in versus Th concentration.Figure 15
Pu239out/Pu239in versus Th concentration.Figure 16
Pu240out/Pu240in versus Th concentration.Figure 17
Pu241out/Pu241in versus Th concentration.Figure 18
Pu242out/Pu242in versus Th concentration.Figure 19
Putotout/Putotin versus Th concentration.Starting from the previous figures, it is possible to highlight some interesting features.(i)
k-infinite at BOC increases with initial Th content (Figure 11).(ii)
Regardless of the initial fuel composition,k-infinite versus burn-up curves show the same qualitative trend (Figure 11).(iii)
Regardless of the initial fuel composition, the only Pu isotope increasing from BOC to EOC is Pu242 (Figure 18).(iv)
In order to achieve a reasonably high burn-up (i.e., 100 GWd/tHM at least), initial Th concentration should not be larger than 80% (Figure11).(v)
The presence of 30% Pu in fresh fuel is the threshold to obtain a discharge burn-up higher than 190 GWd/tHM (Figure11).(vi)
Increasing Th fraction entails an even smaller Pu consumption per unit energy, although this trend is not monotonous (Figure19). Up to a Th concentration not higher than 70%, the smallest fissile Pu consumption per unit energy corresponds to an initial Th percentage that is equal to 67.5% (Figures 15 and 17).(vii)
Up to a Pu concentration equal to 70%, we have the largest Pu mass at EOC, if the initial Th content is 67.5. Then, this composition corresponds to the best exploitation of Th itself. What is more, a smaller Pu consumption means a smaller production of Am and Cm, which generally build up in HTR cores. This behaviour can be considered positive or negative depending on what we are aiming at.(viii)
All the SNF compositions here obtained are strongly proliferation resistant, because of their high content in even Pu nuclides (Figures14, 16, and 18).(ix)
Except for Pu238 and Pu242, there are not significant differences changing the Th initial concentration from 0% to 50%. That means, up to 50% in Pu content, that Th contribution to energy production is substantially low (Figure 13). Conversely, Th content influences in a negative way the possibility of reaching higher burn-ups (Figure 11). Above the point corresponding to 50% initial Th percentage, the behaviour is much more determined by Th content. The highest Th exploitation per single cycle is reached with 67.5% initial Th content. That entails also a lower transmutation rate of Pu in Am and Cm (positive from the long-term radiotoxicity point of view).


## 5.1. Calculation Parameters and Models

In this research, on the basis of the advantages outlined in previous paragraphs, we investigated the main neutronic characteristics related to the utilization of Th as fertile material in PBMR loaded with a Pu/Th fuel, by means of MCNP5  [17] code.The model used to perform these analyses is described in the following lines.In our simulations Pu-RG (i.e., Plutonium produced during irradiation of UO2 fuel in a typical LWR) has been chosen as driver fuel. In principle,  Pu-RG is not suitable for proliferation purposes, due to its high content of even isotopes (Pu238, Pu240, Pu242), which are characterized by a high spontaneous fission rate and, at least in the case of Pu238, by a high decay heat as well. Th content in the  considered pebble has been varied from 0% to 100% of the HM constituting the fuel. Table 2 shows the main parameters of the analyzed PBMR fuel elements, while Table 3 shows the Pu-RG (1st   gen. Pu) isotopic vector.Table 2
(Pu-Th)O2 fuel general characteristics in accordance with PUMA Project [10, 11].
DescriptionValueFuel pebbleFuel pebble outer radius3.0 cmThickness of fuel free zone0.5 cmTotal heavy metal loading per fuel pebble2.0 gMatrix density1.74 g/cm3Packing fraction in pebble bed61%Coated particleKernel Coating MaterialC/C/SiC/CLayer thickness90/40/35/40μmLayer densities1.05/1.90/3.18/1.90 g/cm3Fuel kernelFuel kernel diameter200 micronKernel material type(Pu-Th)O2.0Th-oxide density9.55 g/cm3 (95% TD) [13]Pu-oxide density10.89 g/cm3 (95% TD) [12]Table 3
Composition of 1st  generation Pu [10, 11].
Pu isotopesFirst Generation (weight %)2382.5923953.8524023.6624113.132426.78Additionally, we fixed the following calculation parameters.(i)
An infinite lattice of stochastically arranged pebble has been considered for each fuel composition, by modelling a single sphere (with the He gap associated) with white boundary conditions.(ii)
Each of the considered pebbles contains 2 g HM (Pu/Th) [5].(iii)
Due to the lack of data and for the sake of simplicity, the density of (Pu,Th)O2 fuel (with Th content varying from 0% to 100% of HM) has been assumed as varying linearly from the density of PuO2 kernels (10.89 g/cm3  [12]) to the density of  ThO2 kernels (9.55 g/cm3  [13]), as a function of the Th fraction.(iv)
Geometric error reduction: the pebbles have been modelled as actual as possible (Figure15), also on the basis of the results obtained in [16].(v)
No temperature profile inside the pebble: the temperature is uniform everywhere and equal to 1200 K.(vi)
The cross section libraries used in this work are JEFF3.1 [18].(vii)
Helium Gap: according to the average packing fraction of the PBMR core (61% [8]) the helium gap has been explicitly modelled.Figure8 shows an x-y view of the pebble model adopted to perform these calculations. Please note that clipped CPs along the boundary of the fuelled zone have been eliminated.


## 5.2. Preliminary Results

The variation of thek-infinite at BOC as a function of Th percentage is shown in Figure 9. Analysing the trend of this curve, we can note that some trends can be recognised on the basis of Th content:(i)
0% ≤ Th% < 10%: the k-infinite value remains near 1.26 (flat trend);(ii)
10% ≤ Th% < 30%: a slightly positive gradient is shown (the k-infinite increases from 1.26 to 1.27);(iii)
30% ≤ Th% < 84.5%: the value of k-infinite increases linearly from ~1.27 and to ~1.39;(iv)
84.5% ≤ Th% < 100%: the k-infinite decreases from the peak at 84.5% (k∞=1.39274, σ=0.00154) to 1.Figure 9
Variation ofk-infinite at BOC as a function of Th content.In order to complete this analysis, the fuel behaviour as a function of burn-up has to be considered. All the following burn-up calculations has been  performed by means of the BGCore code [19].


## 5.3. BGCore Code

The burn-up calculations were performed with BGCore. BGCore is a software package for comprehensive computer simulation of nuclear reactor systems and their fuel cycles. It consists of a number of modules:(i)
3-dimensional coupled neutronic and thermal hydraulic module, which calculates the reactor core criticality, as well as the core power and temperature distributions; only steady state core analysis is considered at this stage of the BG-Core package development; no thermal feedback features were used in this analysis;(ii)
fuel depletion and decay module, which calculates the fuel isotopic composition;(iii)
auxiliary modules for basic data processing and management.


## 5.4.k-Effective versus k-Infinite Burn-Up Calculations Comparison

The “real’’k-effective in PBMR depends  not only from the material and the core geometry, as in LWR reactors, but also from a random contribution due to pebbles recirculation [20]. Then, a not-negligible approximation is introduced by calculating the  k-effective fixing the pebble arrangement in the core. In addition, even adopting this approximation, k-effective calculations are quite expensive form the computational point of view.Hence, in order to roughly estimate the differences betweenk-effective  model (i.e., fixed geometry) and k-infinite one, we made a preliminary evaluation, performing two burn-up calculations with a representative fuel composition. Particularly, a fuel composed of 50% Th and 50% Pu was adopted. This choice was also due to the fact that the results  obtained can be directly compared with an independent calculation performed by FZJ [15].Figure10 and Table 4 show our main results.Table 4
Puout/Puin for  k-infinite and  k-effective models.
Pu238out/Pu238inPu239out/Pu239inPu240out/Pu240inPu241out/Pu241inPu242out/Pu242inPutotout/Putotink-infinite43,76%5,46%23,25%48,63%189,23%28,64%k-effective42,57%9,05%53,65%55,51%185,46%38,38%Figure 10
k-infinite and k-effective trends as a function of burn-up (expressed in MWd/tHM).Thus, we can outline some considerations.(i)
If we adoptk-inf = 1 as exit conditions for the burn-up calculations, we obtain a reasonable discharge burnup (comparable to that one we can obtain in a real pebble-bed HTR for a discharged pebble with the same initial composition).(ii)
The shape ofk trends versus burn-up are quite similar.(iii)
The ratios (OUT/IN) among Pu isotopes have (in the two cases) the same order of magnitude.(iv)
The total ratioPuout/Puin in the two cases considered is comparable.(v)
Finally the discharge burn-up obtained in the FZJ calculation is quite similar (Δ~5%) to that we obtained assuming k=1 as exit condition for the k-infinite calculation (see Table 5).Table 5
Comparison betweenk-infinite model and FZJ [15] results.
Heavy metal load and composition per pebbleFZJ calculation—discharge burn-up  [MWd/kgHM]k-infinite calculation—discharge burn-up  [MWd/kgHM)]1 g Pu, 1.0 g Th-232332.38314.42For the reasons described above and because of the kind of study (preliminary parametric evaluation), we performed the calculations shown in the following paragraphs using a simplifiedk-infinite model instead of a k-effective one.


## 5.5. Burn-Up Results

As anticipated above, we analysed some Th/Pu compositions by means of BGCore, in order to assess the fuel behaviour as a function of burn-up. The basic assumptions for these calculations are the following:(i)
times steps of 25 days at full power,(ii)
power per pebble fixed at 1 kW,(iii)
uniform fuel temperature (1200 K) inside the pebble,(iv)
discharge burn-up is set at the time whenk-infinite becomes equal to 1, to give a preliminary estimation of the theoretical nuclear burn-up.(As already underlined, please note that the real PBMR will not be fuelled with pebbles at the same burn-up level, but with a mixture of pebbles at different burn-ups [21]: the real discharge burn-up in principle could be extremely high [21], i.e., up to twice the assumed discharge burnup. On the other hand, we calculated the k-infinite instead of k-effective (so we have not took into account leakages and effects related to the structures surrounding the core)).The results concerningk-inf are summarized in Figure 11. Please note that all the standard deviations are around 100 pcm, which is the value commonly accepted in this kind of calculation by means of MC codes.Figure 11
k-infinite variation as a function of burn-up (expressed in MWd/tHM).As anticipated in the introduction, we evaluated the fuel behaviour on the basis of the following parameters:(i)
the total energy produced,(ii)
the energy produced per initial Pu mass,(iii)
the ratio between final (discharged) and initial (loaded) masses for the Pu isotopes.The obtained results are shown in  Figures12–19.Figure 12
Produced energy (from BOC up to the time at whichk-infinite becomes equal to unity) versus Th concentration.Figure 13
Energy produced pergPuin  versus Th concentration.Figure 14
Pu238out/Pu238in versus Th concentration.Figure 15
Pu239out/Pu239in versus Th concentration.Figure 16
Pu240out/Pu240in versus Th concentration.Figure 17
Pu241out/Pu241in versus Th concentration.Figure 18
Pu242out/Pu242in versus Th concentration.Figure 19
Putotout/Putotin versus Th concentration.Starting from the previous figures, it is possible to highlight some interesting features.(i)
k-infinite at BOC increases with initial Th content (Figure 11).(ii)
Regardless of the initial fuel composition,k-infinite versus burn-up curves show the same qualitative trend (Figure 11).(iii)
Regardless of the initial fuel composition, the only Pu isotope increasing from BOC to EOC is Pu242 (Figure 18).(iv)
In order to achieve a reasonably high burn-up (i.e., 100 GWd/tHM at least), initial Th concentration should not be larger than 80% (Figure11).(v)
The presence of 30% Pu in fresh fuel is the threshold to obtain a discharge burn-up higher than 190 GWd/tHM (Figure11).(vi)
Increasing Th fraction entails an even smaller Pu consumption per unit energy, although this trend is not monotonous (Figure19). Up to a Th concentration not higher than 70%, the smallest fissile Pu consumption per unit energy corresponds to an initial Th percentage that is equal to 67.5% (Figures 15 and 17).(vii)
Up to a Pu concentration equal to 70%, we have the largest Pu mass at EOC, if the initial Th content is 67.5. Then, this composition corresponds to the best exploitation of Th itself. What is more, a smaller Pu consumption means a smaller production of Am and Cm, which generally build up in HTR cores. This behaviour can be considered positive or negative depending on what we are aiming at.(viii)
All the SNF compositions here obtained are strongly proliferation resistant, because of their high content in even Pu nuclides (Figures14, 16, and 18).(ix)
Except for Pu238 and Pu242, there are not significant differences changing the Th initial concentration from 0% to 50%. That means, up to 50% in Pu content, that Th contribution to energy production is substantially low (Figure 13). Conversely, Th content influences in a negative way the possibility of reaching higher burn-ups (Figure 11). Above the point corresponding to 50% initial Th percentage, the behaviour is much more determined by Th content. The highest Th exploitation per single cycle is reached with 67.5% initial Th content. That entails also a lower transmutation rate of Pu in Am and Cm (positive from the long-term radiotoxicity point of view).


## 6. Conclusions

Some aspects related to Thorium exploitation were outlined, with a particular emphasis on its suitability for working in pebble-bed HTRs as a fertile element in a Th-Pu fuel cycle. That is a promising kind of fuel, which was also considered in the framework of the EU PUMA project. Particularly, the influence of the Th/Pu weight fraction at BOC in a typical HTR pebble was analysed as far as the reactivity trend versus burn-up, the energy produced per Pu unit mass, and the Pu isotopic composition at EOC are concerned. On the basis of a preliminary comparative analysis (between ak-infinite and a k-effective model) and because of the kind of assessment presented here (i.e., a preliminary parametric evaluation), we chose a (simplified) k-infinite model instead of a k-effective one to perform our calculations.Some general trends depending on the initial Th content can be recognized, and particularly the following.(i)
k-infinite versus burn-up shows the same qualitative trend regardless of the initial Th content (Figure 11).(ii)
k-infinite at BOC increases with initial Th content in a proportional way, and its decrease is as sharp as the Th content itself is high (Figure 11). That may make the reactivity control somewhat difficult.(iii)
The final burn-up increases with initial Th percentage decrease (Figure11). If the content of Th is larger than 80%, the system becomes subcritical after less than 100 GWd/tHM burn-up (Figure 11).(iv)
Up to a Th concentration not higher than 70%, the smallest fissile Pu consumption per unit energy corresponds to an initial Th percentage that is equal to 67.5% (Figures15 and 17).(v)
All Pu nuclides except to Pu242 decrease from BOC to EOC regardless of the initial Th content (Figure 15).(vi)
The point corresponding of a minimum in Pu consumption per unit energy produced (i.e., 67.5% Th, Figure13) corresponds to the highest Th exploitation and to the lowest transmutation of Pu into heavier nuclides (Pu242, Am, and Cm).These results suggest us that there is not a “best” composition for the PBMR Pu/Th fuel, although deeper investigations need to be performed in order to draw final conclusions. Particularly, at least from the neutronic point of view, the influence of Th percentage on the reactivity coefficients as well as on the delayed neutron fraction should be determined. Thus, at the moment it is possible to state that some optimized Th percentage in the initial Pu/Th fuel could be suggested on the basis of the aim we are trying to reach.(i)
If we aim at Pu reducing, a Pu fertile-free fuel or a Th concentration~20% has to be adopted (Figure 19), because these fuels allow us to reach the highest burn-up and the highest Pu consumption with a single irradiation cycle.(ii)
If we aim at maximising the natural resource exploitation (i.e., of both Th and U), Th content has to be around 67.5%.(iii)
In any case, there does not seem to be any advantages in inserting more than 70% of Th in the fresh fuel (higher initial reactivity to be controlled, higher content of heavier elements, shorter cycle length, smaller consumption of the driver fuel).This work, although quite deep, has to be considered not completed; so it needs of some further work. Looking at the future developments of these research, in order to have a complete knowledge of the analyzed problem, it would be necessary to evaluate the following:(i)
reactivity coefficients,(ii)
short-  and long-term radiotoxicity evolution,(iii)
materials behaviour and compatibility in the analyzed cases (e.g., under high fluence, at high temperature, etc.),(iv)
fuel cycles in the frame of the energetic scenarios (taking into account also environmental and economical aspects),(v)
nonproliferation issues.



---

*Source: 749736-2009-11-30.xml*