# A Simplified Supercritical Fast Reactor with Thorium Fuel

Authors: Peng Zhang; Kan Wang; Ganglin Yu
Journal: Science and Technology of Nuclear Installations
 (2014)

---

## Abstract

Super-Critical water-cooled Fast Reactor (SCFR) is a feasible option for the Gen-IV SCWR designs, in which much less moderator and thus coolant are needed for transferring the fission heat from the core compared with the traditional LWRs. The fast spectrum of SCFR is useful for fuel breeding and thorium utilization, which is then beneficial for enhancing the sustainability of the nuclear fuel cycle. A SCFR core is constructed in this work, with the aim of simplifying the mechanical structure and keeping negative coolant void reactivity during the whole core life. A core burnup simulation scheme based on Monte Carlo lattice homogenization is adopted in this study, and the reactor physics analysis has been performed with DU-MOX and Th-MOX fuel. The main issues discussed include the fuel conversion ratio and the coolant void reactivity. The analysis shows that thorium-based fuel can provide inherent safety for SCFR without use of blanket, which is favorable for the mechanical design of SCFR.


---

## Body

## 1. SCFR Conceptual Design

In SCWR much less coolant water is needed for cooling the reactor, so a fast spectrum option (SCFR) is possible. With a harder spectrum, SCFR can have a higher conversion ratio. Thus, physically, SCFR is a kind of high conversion LWR design, for which the main approach to increase conversion ratio is to decrease the ratio of water to heavy metal, so as to decrease the moderation effect of water and make the spectrum hard [1–4].The coolant void reactivity (CVR) is a crucial safety aspect of fast reactor design. Main concern arising from coolant voiding is the hardening of the neutron spectrum, which increases fast fission in both seed and blanket fuel regions and also increases neutron leakage at the same time. In order to decrease the CVR, the flat core design is usually adopted to increase the neutron leakage. But it is not an economical option for SCFR, since it is operated at a very high pressure and so a bigger and thicker pressure vessel is needed to contain the core. Thus, increasing neutron absorption at void condition is the key for achieving negative CVR in SCFR, and so blankets are very important since they are the only regions in the core where neutron capture strongly prevails over fission. A typical SCFR blanket design is shown in Figure1, in which a ZrH layer is adopted to slow down the fast neutron leaked from seed assemblies at void condition and decrease the fast fission at blanket assemblies [5–8].Figure 1
A typical SCFR blanket assembly design.However, in SCFR the coolant in blanket assembly is very cold as very little heat generated, so the flow channels should be separated from those of the seed assemblies. This will make the flow path and the core mechanical design complicated. In order to simplify the mechanical design, we need to think out another way to keep CVR negative. Introducing thorium fuel is a possible option [4, 9–11], so this paper work mainly focuses on it.A SCFR core conceptual design is proposed, as shown in Figure2. No blanket assembly is adopted in the design. The core power density will increase in this case since a blanket assembly has much lower power than a seed assembly, and the mechanical design can be simplified as it is not needed to separate coolant flow in blanket and seed assemblies. In addition, it is possible to achieve negative CVR by changing the fuel design such as thorium fuel, which is a main purpose of this work.Figure 2
Conceptual design of SCFR core and assembly.The SCFR core can be divided into two radial regions, the inner core and the outer core. The numbers of assemblies for the two parts are generally the same. The inlet coolant flows downwards through outer core assemblies first. After mixing in the lower plenum, the coolant flows upwards through the inner core assemblies to the upper plenum and then to the outlet. The main design parameters of SCFR core and assembly are listed in Tables1 and 2.Table 1
SCFR core design parameters.
Parameters
Value
Number of assembly
241 (inner 121, outer 120)
Equivalent diameter/m
2.97
Active length/m
3.6
Thermal power/MWt
3498
Electric power/MWe
1532
Average linear heat generation rate/kW/m
16
Average power density/MW/m3
140Table 2
SCFR assembly design parameters.
Parameters
Value
Pellet OD/mm
8.2
Cladding OD/mm
9.62
Cladding thickness/mm
0.63
Cladding material
MA956, Fe-21.5Cr-5.75Al-0.6Ti-0.7Y2O3,7.25 g/cm3
Pitch/mm
10.62
P/D
1.1
Pin-to-channel-box gap/mm
0.5
Channel box thickness/mm
2.5
Gap between adjacent channel box/mm
1.0
Assembly width/mm
182.2


## 2. The “Two-Step” MC Core Burnup Analysis Method

Monte Carlo method has been widely used for the verification and validation of deterministic codes, and also for the analysis of many newly developed nuclear energy systems since it can deal with arbitrary geometry and spectrum configuration [12]. However, the efficiency of Monte Carlo calculation is too low for the ordinary large scale core designs, due to the long simulation time to get reliable results. With the rapid developing of computer technologies and parallel algorithms, the efficiency of Monte Carlo simulation has been greatly increased. Coupled neutronic thermal-hydraulic analysis of full core with Monte Carlo has been done [13], but core burnup analysis with continuous energy Monte Carlo is still difficult nowadays.A “two-step” Monte Carlo core burnup analysis method is proposed in this work. Similar to the deterministic code systems, the core calculation scheme is divided into two steps. First, the assembly simulation is done with continuous energy Monte Carlo, and the assembly group cross-sections are tallied at the same time. Second, with these group constants, the core simulation is done with multigroup Monte Carlo [14], which can greatly increase the core simulation efficiency. The core burnup analysis can be done if the assembly group constants can be prepared for different burnup states.The “two-step” core simulation scheme is illustrated in Figure3. MCBurn, a coupling code of MCNP_CE (continuous energy mode of MCNP) and ORIGEN, is adopted for assembly burnup simulations [15]. The assembly group constants are gathered and organized to prepare group cross-sections for core calculations. MCCBurn code is developed to manage core burnup simulation, which uses MCNP_MG (multigroup mode of MCNP) to do core criticality simulations.Figure 3
Two-step Monte Carlo core burnup analysis system.In order to justify the calculation tools, a BWR assembly with moderator density (DM) 0.3 g/cm3 has been modeled using MCBurn and MCCBurn, with KCODE card of MCNP as “KCODE 5000 1.0 30 90.” The burnup calculation results are shown in Figure 4. It can be seen that the differences of k
inf
⁡ between MCBurn and MCCBurn are well below 0.2%, which reveals that the group cross-sections generated using our method is right.Figure 4
A BWR assembly model and its burnup calculation results.


## 3. Reactor Physics Characteristics of SCFR with DU-MOX and Th-MOX

In order to concentrate on the reactor physics characteristics analysis of SCFR, the thermal-hydraulic parameters have been simplified. The coolant densities of inner and outer core are set to be the average values along core height, 0.13 g/cm3 and 0.6 g/cm3, respectively. Besides, a pin cell is modeled to represent the assembly, which means the radial heterogeneity of assembly is ignored. This may introduce some errors to the results, but we think this can be neglected for the preliminary feasibility study.We have chosen DU-MOX and Th-MOX for the comparing study, which are listed in Table3. The isotope weight vector of reactor grade plutonium (RGPu) is as follows: 2.8% 238Pu, 54.4% 239Pu, 22.8% 240Pu, 11.8% 241Pu, 7% 242Pu, and 1.2% 241Am. The depleted uranium (DU) is assumed to have 0.25 weight% 235U.Table 3
The different fuel forms chosen for study.
Fuel type
Fuel composition
Enrichment definition
DU-MOX
RGPuO2 + DUO2
(239Pu + 241Pu)/HM
Th-MOX
RGPuO2 + ThO2
(239Pu + 241Pu)/HM


### 3.1. Spectrum Comparison

Supposing the fuel is DU-MOX with an enrichment of 7.5%, the inner and outer core spectrums are shown in Figure5. The spectrums of PWR and sodium-cooled Fast Reactor (SFR) are also included for comparison. It can be seen that there is no PWR-like thermal peak for SCFR’s spectrum. Spectrum of inner core is harder than outer core since the coolant density is lower.Figure 5
Spectrum comparison of SCFR (inner and outer core), PWR, and SFR.


### 3.2. Cell Burnup Performance

#### 3.2.1. DU-MOX Fuel Cell

The main burnup performances of inner and outer cells with different enrichments are shown in Figures6 and 7. Where CVR is the coolant void reactivity; CR, conversion ratio, is defined as the ratio of capture rate of fertile nuclides (232Th, 238U, 238Pu, and 240Pu) to absorption rate of fissile nuclides (233U, 235U, 239Pu, and 241Pu); and FIR, fissile inventory ratio, is defined as the ratio of total amount of fissile material to the initial loading fissile material.Figure 6
Inner cell burnup performances with DU-MOX fuel.Figure 7
Outer cell burnup performances with DU-MOX fuel.From these figures we can see that the cell burnup performances are sensitive to fuel enrichment. With the increasing of fuel enrichment: cell initialk
inf
⁡ increases notably, and the differences decrease with burnup; CVR goes positive rapidly, and it increases slowly with burnup; CR and FIR decreases significantly.


#### 3.2.2. Th-MOX Fuel Cell

The Th-MOX fuel cell burnup performances for inner and outer core are shown in Figures8 and 9. The same trends can be observed for the parameters when fuel enrichment and burnup changed. A great advantage of adopting thorium fuel is that the CVR of cell can be kept negative for a wide range of enrichment and for the whole life cycle. Considering the negative effect of neutron leakage, the CVR of SCFR core will be negative. So in the following core analyses, CVR is not discussed in detail.Figure 8
Inner cell burnup performances with Th-MOX fuel.Figure 9
Outer cell burnup performances with Th-MOX fuel.Thus, thorium fuel is very suitable for SCFR, and Th-MOX is adopted in the following study. Considering the critical ability, the fuel enrichment for inner and outer core is chosen as 8.5% and 7.5%, respectively.


### 3.3. Core Burnup Performance

#### 3.3.1. Core Loading Pattern

In order to get the equivalent fuel cycle performances, we need to design the fuel loading patterns. The inner and outer core are considered independently, since they have different reactor physics characteristics.Assuming that the average linear heat generation rate is 16 kW/m, the specific power density should be about 34 kW/kgHM. The object burnup is supposed to be 60~70 MWd/kgHM, so an assembly should stay in core for 1800~2000 days, and loading pattern can be 3-, 4-, or 5-batch loading pattern. In this work, the cycle length is fixed as 550 days. Considering the higher power density of inner core, 3-batch loading pattern is suitable for inner core, while 4- or 5-batch loading pattern for outer core. The loading schemes for 1/6 symmetric core are shown in Figure10.Core loading scheme (m.n: m—the assembly group number; n—cycles of the assembly has been burned).
(a)
Inner: 3-batch; outer: 4-batch
(b)
Inner: 3-batch; outer: 5-batch


#### 3.3.2. Multicycle Simulation

The multicycle core simulation results are shown in Figure11. The in-3-out-4 batch loading pattern is adopted. PPF (peak power factor) is the ratio of the maximum assembly power to the average assembly power. The k
eff of core increases slowly with burnup. The average and maximum discharge burnups of inner core (57 MWd/kgHM and 61 MWd/kgHM) are lower than that of the outer core (71 MWd·kgHM and 74 MWd/kgHM).Figure 11
Multicycle burnup simulation results for Th-MOX core.Since the inner and outer core regions are not coupled tightly, the fuel forms can be different for the two regions. Considering the CVR, DU-MOX should not be loaded in the inner region, since the cell burnup results indicate that local CVR may be positive at a burnup higher than 40 GWd/tHM. Considering that many neutrons may leak from the core at void condition, it is possible to load DU-MOX fuel in the outer region while still keeping a negative CVR for the whole core. Th-MOX fuel is loaded in the inner region. The enrichments for Th-MOX and DU-MOX are 8.5% and 6%. In-3-out-5 batch loading pattern is adopted, and the multicycle core simulation results are shown in Figure12. The difference of discharge burnups between inner and outer core region is decreased. The main core parameters at equivalent cycle are listed in Table 4.Table 4
SCFR core parameters with different fuel options.
Parameters
Th-MOX
ThMOX + MOX
Fuel forms
Inner
(RGPu + Th)O2
(RGPu + Th)O2
Outer
(RGPu + Th)O2
(RGPu + DU)O2
Enrichment/%
Inner
8.5
8.5
Outer
7.5
6.0
Batch number
Inner
3
3
Outer
4
5
Corek
eff (0~550 EFPDs)
1.031~1.049
1.030~1.050
Maximum PPF
1.48
1.98
Average discharge burnup/MWd/kgHM
Inner
57.2
68.4
Outer
71.0
64.5
Whole core
62.9
67.0
Maximum discharge burnup/MWd/kgHM
Inner
61.1
76.0
Outer
74.2
71.8
FIR
Inner
1.01
1.01
Outer
0.81
1.01
Whole core
0.93
1.01Figure 12
Multicycle burnup simulation results for Th-MOX+DU-MOX core.


## 3.1. Spectrum Comparison

Supposing the fuel is DU-MOX with an enrichment of 7.5%, the inner and outer core spectrums are shown in Figure5. The spectrums of PWR and sodium-cooled Fast Reactor (SFR) are also included for comparison. It can be seen that there is no PWR-like thermal peak for SCFR’s spectrum. Spectrum of inner core is harder than outer core since the coolant density is lower.Figure 5
Spectrum comparison of SCFR (inner and outer core), PWR, and SFR.


## 3.2. Cell Burnup Performance

### 3.2.1. DU-MOX Fuel Cell

The main burnup performances of inner and outer cells with different enrichments are shown in Figures6 and 7. Where CVR is the coolant void reactivity; CR, conversion ratio, is defined as the ratio of capture rate of fertile nuclides (232Th, 238U, 238Pu, and 240Pu) to absorption rate of fissile nuclides (233U, 235U, 239Pu, and 241Pu); and FIR, fissile inventory ratio, is defined as the ratio of total amount of fissile material to the initial loading fissile material.Figure 6
Inner cell burnup performances with DU-MOX fuel.Figure 7
Outer cell burnup performances with DU-MOX fuel.From these figures we can see that the cell burnup performances are sensitive to fuel enrichment. With the increasing of fuel enrichment: cell initialk
inf
⁡ increases notably, and the differences decrease with burnup; CVR goes positive rapidly, and it increases slowly with burnup; CR and FIR decreases significantly.


### 3.2.2. Th-MOX Fuel Cell

The Th-MOX fuel cell burnup performances for inner and outer core are shown in Figures8 and 9. The same trends can be observed for the parameters when fuel enrichment and burnup changed. A great advantage of adopting thorium fuel is that the CVR of cell can be kept negative for a wide range of enrichment and for the whole life cycle. Considering the negative effect of neutron leakage, the CVR of SCFR core will be negative. So in the following core analyses, CVR is not discussed in detail.Figure 8
Inner cell burnup performances with Th-MOX fuel.Figure 9
Outer cell burnup performances with Th-MOX fuel.Thus, thorium fuel is very suitable for SCFR, and Th-MOX is adopted in the following study. Considering the critical ability, the fuel enrichment for inner and outer core is chosen as 8.5% and 7.5%, respectively.


## 3.2.1. DU-MOX Fuel Cell

The main burnup performances of inner and outer cells with different enrichments are shown in Figures6 and 7. Where CVR is the coolant void reactivity; CR, conversion ratio, is defined as the ratio of capture rate of fertile nuclides (232Th, 238U, 238Pu, and 240Pu) to absorption rate of fissile nuclides (233U, 235U, 239Pu, and 241Pu); and FIR, fissile inventory ratio, is defined as the ratio of total amount of fissile material to the initial loading fissile material.Figure 6
Inner cell burnup performances with DU-MOX fuel.Figure 7
Outer cell burnup performances with DU-MOX fuel.From these figures we can see that the cell burnup performances are sensitive to fuel enrichment. With the increasing of fuel enrichment: cell initialk
inf
⁡ increases notably, and the differences decrease with burnup; CVR goes positive rapidly, and it increases slowly with burnup; CR and FIR decreases significantly.


## 3.2.2. Th-MOX Fuel Cell

The Th-MOX fuel cell burnup performances for inner and outer core are shown in Figures8 and 9. The same trends can be observed for the parameters when fuel enrichment and burnup changed. A great advantage of adopting thorium fuel is that the CVR of cell can be kept negative for a wide range of enrichment and for the whole life cycle. Considering the negative effect of neutron leakage, the CVR of SCFR core will be negative. So in the following core analyses, CVR is not discussed in detail.Figure 8
Inner cell burnup performances with Th-MOX fuel.Figure 9
Outer cell burnup performances with Th-MOX fuel.Thus, thorium fuel is very suitable for SCFR, and Th-MOX is adopted in the following study. Considering the critical ability, the fuel enrichment for inner and outer core is chosen as 8.5% and 7.5%, respectively.


## 3.3. Core Burnup Performance

### 3.3.1. Core Loading Pattern

In order to get the equivalent fuel cycle performances, we need to design the fuel loading patterns. The inner and outer core are considered independently, since they have different reactor physics characteristics.Assuming that the average linear heat generation rate is 16 kW/m, the specific power density should be about 34 kW/kgHM. The object burnup is supposed to be 60~70 MWd/kgHM, so an assembly should stay in core for 1800~2000 days, and loading pattern can be 3-, 4-, or 5-batch loading pattern. In this work, the cycle length is fixed as 550 days. Considering the higher power density of inner core, 3-batch loading pattern is suitable for inner core, while 4- or 5-batch loading pattern for outer core. The loading schemes for 1/6 symmetric core are shown in Figure10.Core loading scheme (m.n: m—the assembly group number; n—cycles of the assembly has been burned).
(a)
Inner: 3-batch; outer: 4-batch
(b)
Inner: 3-batch; outer: 5-batch


### 3.3.2. Multicycle Simulation

The multicycle core simulation results are shown in Figure11. The in-3-out-4 batch loading pattern is adopted. PPF (peak power factor) is the ratio of the maximum assembly power to the average assembly power. The k
eff of core increases slowly with burnup. The average and maximum discharge burnups of inner core (57 MWd/kgHM and 61 MWd/kgHM) are lower than that of the outer core (71 MWd·kgHM and 74 MWd/kgHM).Figure 11
Multicycle burnup simulation results for Th-MOX core.Since the inner and outer core regions are not coupled tightly, the fuel forms can be different for the two regions. Considering the CVR, DU-MOX should not be loaded in the inner region, since the cell burnup results indicate that local CVR may be positive at a burnup higher than 40 GWd/tHM. Considering that many neutrons may leak from the core at void condition, it is possible to load DU-MOX fuel in the outer region while still keeping a negative CVR for the whole core. Th-MOX fuel is loaded in the inner region. The enrichments for Th-MOX and DU-MOX are 8.5% and 6%. In-3-out-5 batch loading pattern is adopted, and the multicycle core simulation results are shown in Figure12. The difference of discharge burnups between inner and outer core region is decreased. The main core parameters at equivalent cycle are listed in Table 4.Table 4
SCFR core parameters with different fuel options.
Parameters
Th-MOX
ThMOX + MOX
Fuel forms
Inner
(RGPu + Th)O2
(RGPu + Th)O2
Outer
(RGPu + Th)O2
(RGPu + DU)O2
Enrichment/%
Inner
8.5
8.5
Outer
7.5
6.0
Batch number
Inner
3
3
Outer
4
5
Corek
eff (0~550 EFPDs)
1.031~1.049
1.030~1.050
Maximum PPF
1.48
1.98
Average discharge burnup/MWd/kgHM
Inner
57.2
68.4
Outer
71.0
64.5
Whole core
62.9
67.0
Maximum discharge burnup/MWd/kgHM
Inner
61.1
76.0
Outer
74.2
71.8
FIR
Inner
1.01
1.01
Outer
0.81
1.01
Whole core
0.93
1.01Figure 12
Multicycle burnup simulation results for Th-MOX+DU-MOX core.


## 3.3.1. Core Loading Pattern

In order to get the equivalent fuel cycle performances, we need to design the fuel loading patterns. The inner and outer core are considered independently, since they have different reactor physics characteristics.Assuming that the average linear heat generation rate is 16 kW/m, the specific power density should be about 34 kW/kgHM. The object burnup is supposed to be 60~70 MWd/kgHM, so an assembly should stay in core for 1800~2000 days, and loading pattern can be 3-, 4-, or 5-batch loading pattern. In this work, the cycle length is fixed as 550 days. Considering the higher power density of inner core, 3-batch loading pattern is suitable for inner core, while 4- or 5-batch loading pattern for outer core. The loading schemes for 1/6 symmetric core are shown in Figure10.Core loading scheme (m.n: m—the assembly group number; n—cycles of the assembly has been burned).
(a)
Inner: 3-batch; outer: 4-batch
(b)
Inner: 3-batch; outer: 5-batch


## 3.3.2. Multicycle Simulation

The multicycle core simulation results are shown in Figure11. The in-3-out-4 batch loading pattern is adopted. PPF (peak power factor) is the ratio of the maximum assembly power to the average assembly power. The k
eff of core increases slowly with burnup. The average and maximum discharge burnups of inner core (57 MWd/kgHM and 61 MWd/kgHM) are lower than that of the outer core (71 MWd·kgHM and 74 MWd/kgHM).Figure 11
Multicycle burnup simulation results for Th-MOX core.Since the inner and outer core regions are not coupled tightly, the fuel forms can be different for the two regions. Considering the CVR, DU-MOX should not be loaded in the inner region, since the cell burnup results indicate that local CVR may be positive at a burnup higher than 40 GWd/tHM. Considering that many neutrons may leak from the core at void condition, it is possible to load DU-MOX fuel in the outer region while still keeping a negative CVR for the whole core. Th-MOX fuel is loaded in the inner region. The enrichments for Th-MOX and DU-MOX are 8.5% and 6%. In-3-out-5 batch loading pattern is adopted, and the multicycle core simulation results are shown in Figure12. The difference of discharge burnups between inner and outer core region is decreased. The main core parameters at equivalent cycle are listed in Table 4.Table 4
SCFR core parameters with different fuel options.
Parameters
Th-MOX
ThMOX + MOX
Fuel forms
Inner
(RGPu + Th)O2
(RGPu + Th)O2
Outer
(RGPu + Th)O2
(RGPu + DU)O2
Enrichment/%
Inner
8.5
8.5
Outer
7.5
6.0
Batch number
Inner
3
3
Outer
4
5
Corek
eff (0~550 EFPDs)
1.031~1.049
1.030~1.050
Maximum PPF
1.48
1.98
Average discharge burnup/MWd/kgHM
Inner
57.2
68.4
Outer
71.0
64.5
Whole core
62.9
67.0
Maximum discharge burnup/MWd/kgHM
Inner
61.1
76.0
Outer
74.2
71.8
FIR
Inner
1.01
1.01
Outer
0.81
1.01
Whole core
0.93
1.01Figure 12
Multicycle burnup simulation results for Th-MOX+DU-MOX core.


## 4. Conclusion

In this work, a SCFR core conceptual design is proposed. The core can be divided into two radial regions, the inner region and the outer region. The spectrums of the two regions are different, which results in their different reactor physics characteristics. The fuel forms and the loading patterns of the two regions can be designed independently. Unlike the significantk
eff  decrease in the traditional LWRs, the k
eff of the proposed SCFR core increases slowly during the whole life cycle. This reveals that SCFR has good fuel conversion ability. With the adoption of thorium fuel, the CVR can be kept negative during the whole life cycle without use of blanket, which is also beneficial for increasing the power density and simplifying the mechanical design.



---

*Source: 405654-2014-03-10.xml*