# GCFR Coupled Neutronic and Thermal-Fluid-Dynamics Analyses for a Core Containing Minor Actinides

Authors: Diego Castelliti; Eleonora Bomboni; Nicola Cerullo; Guglielmo Lomonaco; Carlo Parisi
Journal: Science and Technology of Nuclear Installations
 (2009)

---

## Abstract

Problems about future energy availability, climate changes, and air quality seem to play an important role in energy production. While current reactor generations provide a guaranteed and economical energy production, new nuclear power plant generation would increase the ways and purposes in which nuclear energy can be used. To explore these new technological applications, several governments, industries, and research communities decided to contribute to the next reactor generation, called “Generation IV.” Among the six Gen-IV reactor designs, the Gas Cooled Fast Reactor (GCFR) uses a direct-cycle helium turbine for electricity generation and for aCO2-free thermochemical production of hydrogen. Additionally, the use of a fast spectrum allows actinides transmutation, minimizing the production of long-lived radioactive waste in an integrated fuel cycle. This paper presents an analysis of GCFR fuel cycle optimization and of a thermal-hydraulic of a GCFR-prototype under steady-state and transient conditions. The fuel cycle optimization was performed to assess the capability of the GCFR to transmute MAs, while the thermal-hydraulic analysis was performed to investigate the reactor and the safety systems behavior during a LOFA. Preliminary results show that limited quantities of MA are not affecting significantly the thermal-fluid-dynamics behavior of a GCFR core.


---

## Body

## 1. Introduction: The Generation IV Project

Generation IV proposals have been specifically advanced in order to meet the new challenges for nuclear energy production, increasing also the ways and purposes in  which nuclear energy can be used. Six new reactor designs have been proposed, all of them aiming to satisfy requirements of reliability, sustainability, economics, and resistance to proliferation that characterize the Gen IV [1].Among the six Gen IV reactor projects, there are two with a thermal  neutronic spectrum (High, Temperature Gas Reactor (HTGR), Supercritical Water Reactor (SCWR)), three with a fast spectrum (Gas-Cooled Fast Reactor  (GCFR), Sodium Fast Reactor (SFR), Lead Fast Reactor (LFR)), and one with epithermal spectrum (Molten Salt Reactor (MSR)).An important Gen IV target is Minor Actinides (MAs) transmutation. This paper presents an analysis of the GCFR fuel cycle with MAs and a thermal-hydraulic analysis of a GCFR-prototype, charged with such composed fuel, under steady-state and transient conditions. The fuel cycle has been designed in order to asses the capability of the GCFR to transmute MAs, while the thermal-hydraulic investigation was performed in order to investigate the reactor and the safety system behaviors during one of the Design Basis Accidents (DBAs), the Loss-of-Flow Accident (LOFA). Sensitivity analyses are also presented for assessing the effects of the presence of different MAs concentration on the fuel performance.


## 2. The Gas-Cooled Fast Reactor

Generation IV reactor projects are finalized to answer to new challenges that nuclear energy must face in this new century. Those mandatory projects have to meet Gen IV criteria [1], namely:(1)
sustainability;(2)
economics;(3)
proliferation-resistance;(4)
reliability.Criteria (1), (2), and (3)  include, as important target, the MAs transmutation.The capability of removing the heat produced in the core (both in nominal and accidental conditions) is a key point (particularly for this reactor) of criterion (4); this work (in conjunction with what already published in [2]) gives a positive (even if still partial) answer to this request.The GCFR proposed in the frame of Gen IV initiative (Figure1) uses a direct-cycle helium turbine for electricity generation. This reactor can be employed both for electricity generation and in CO2-free thermochemical production of hydrogen. Besides, the use of a fast spectrum allows the transmutation of the actinides, minimizing the long-lived radioactive waste production in an integrated fuel cycle (one of main targets of Gen IV).Figure 1
GCFR scheme [3].Furthermore, the GCFR fast spectrum makes the use of available fissile and fertile materials (including depleted uranium) possible in a more efficient way compared to the thermal reactors. GCFR has better safety features than liquid metals reactor, and it can be quite flexible in Pu management [3].The GCFR has a high power density, about 100 KW/l, similar to a PWR, but in this case the power has to be removed by a noncondensable gas (helium) instead of water.The proposed GCFR initial design shares several components with HTGR technology, and it uses a direct Brayton cycle (without any Intermediate Heat Exchanger, or IHX). This allows to combine a high thermal efficiency (thanks to the high coolant temperatures) with a high MAs transmutation and fertilization rate (thanks to the fast neutron spectrum). The high coolant output temperatures make the reactor suitable for industrial applications, too.However, the direct Brayton cycle design is considered very challenging; so a “more conventional”  version with a He,N2, or CO2 secondary loop could be preferred [4]. The inevitable performance loss is accepted, in order to achieve better plant safety features. The ETDR  experimental facility moves in this direction  (see Section 3).Preliminary analysis results confirmed that it is possible to have a self-breeding core without any fertile blanket (i.e., Breeding Gain slightly greater than zero in order to meet nonproliferation constraints) from a  starting fuel based on Pu and depleted U. More details on proposed technical solutions for GCFR can be found in  [3, 5–11].


## 3. ETDR: Core and Plant Characteristics

Actual realization of a commercial GCFR is scheduled not before the year 2040. An experimental and demonstration facility is scheduled for the year 2020, with the purpose to test the technical solutions proposed for the GCFR concept. The name of this demonstration facility is Experimental and Testing Demonstration Reactor (ETDR). The ETDR will represent a milestone on the way approaching the final GCFR configuration, thanks to the numerous technical solutions that will be tested despite a quite different power size. While the GCFR is projected to have a nominal thermal power of 2400 (or 3000)MWth (600 MWth in the modular version), ETDR nominal power will be limited at 50 MWth, with no electricity production and a secondary water-cooled loop. It will be the first gas cooled, fast neutron spectrumnuclear plant ever built.ETDR plant scheme is illustrated in Figure2.Figure 2
: ETDR plant overview [12].Figure 3
ETDR flux diagram [12].In Figures 1, 2 and  3 we can recognize.(1)
core and vessel;(2)
main heat exchanger;(3)
main blower;(4)
main cross duct;(5)
decay heat removal loop;(6)
decay heat removal blower (supposed inactive in the LOFA transient);(7)
tertiary air cooling system (not modeled).A more detailed ETDR plant description can be found in [10].Decay Heat Removal (DHR) system is composed by three loops (3×100% redundancy) connected to vessel and should have the same design that will be installed on the commercial 600 MWth GCFR. One of those three loops has no blowers and it relies on natural circulation (NC) only (full passive system), while the second loop is equipped with a blower that is supposed to work only for 24 hours after the DHR valve is opened (few seconds later incidental event) and before NC takes place. Finally, the third DHR loop is totally based on active safety systems (blower). The NC loop is particularly efficient in transient events where there is not a fast depressurization of primary circuit (structural integrity preserved), while active devices are especially suitable for fast depressurization events (LOCA) (Blower in second loop is meant to work with a system residual pressure of about 7÷10 bars (SBLOCA), while blower in third loop can work with He at nearly 1 bar (LBLOCA), requiring a higher pumping power.)As previously said, DHR system is connected to the reactor pressure vessel (RPV). Hot leg is connected to the upper plenum, while helium from cold leg flows into downcomer. Therefore the core is cooled by bottom-to-top coolant flow. During normal operation, reactor power is removed via a counter-current U-tubes helium-water HX. The secondary side is then cooled by forced air circulation.Concerning the DHR HXs, the primary to secondary heat removal is achieved by the use of counter-current straight vertical tubes helium-water HX. NC is established thanks to difference of 15 m between the primary and secondary side thermal barycenters.DHR secondary circuit, containing water pressurized at 10 bars, is exchanging the removed power from the primary side with the final heat sink, a water pool at atmospheric pressure. The height difference between secondary-side thermal barycenters is about 4 m.


## 4.  ETDR: Core Fuelled with Minor Actinides

Several studies and calculations have been performed on thermal-hydraulics of ETDR systems, both in steady-state and transients conditions, but most of these analyses have not considered any neutronic feedback in the core.We tried to obtain some results by analyzing a core charged with a small quantity of  MA. (In the present design it is not possible to exceed10% of MAs content because of reactor control problems.) In fact, one of most important features of this system is the capability to burn MAs, greatly reducing wastes and producing energy by fast fissions as well.In order to make some preliminary evaluations of the dynamic effects on a GCFR core due to MAs presence in the fuel, we performed some point neutron kinetics-(NK-)TH calculations on the ETDR design. For comparison purposes, we chose to analyze Steady State and a LOFA transient because there are validated data already available for these cases [2]. Neutronic data (Doppler coefficient, void coefficient, etc.) implemented in our calculations are taken from the MCNP calculations performed for the GCFR core (600 MHth) see [3, 5–11]. As already anticipated, to perform a complete ETDR LOFA analysis [2], a point NK analysis is necessary, taking into account the variations of reactivity coefficients (e.g., Doppler and void coefficients) due to the fuel isotopic changes.Calculations were performed with and without the presence of MAs in the fuel. The code used for the point NK-TH calculations was the RELAP5-3D  [13] .We performed, at first, MCNP5 [14] calculations, assessing Doppler and void coefficients, delayed neutron fractions, and prompt neutron average life, in order to run a point NK-TH  calculation. We considered, at first, isotope mass fractions shown in Table 1 (first column, LWR discharge). This composition appears to be quite performing, according to criticality, as previously shown in [15].Table 1
Core isotopic compositions considered [15]
IsotopeMass fraction LWR discharge (%)Mass fraction5% MAs (%)U2350.20.571U23879.81180.991Np2370.920.87Pu2380.460.446Pu2399.8049.251Pu2404.3024.278Pu2412.2811.222Pu2421.2411.206AM2410.56—AM242m—0.017AM2430.320.81Cm242—0.001Cm243—0.002Cm2440.10.265Cm245—0.065Cm246—0.005In conclusion, we assumed, from MCNP calculations.(i)
Delayed neutron fractionβ=332pcm.(ii)
Prompt neutron average life=0.364μs.(iii)
Reactivity insertion from 1200 K to 500 K=+98 pcm.(iv)
Reactivity insertion from full depressurization (7.0 MPa→ 0.1 MPa) =+68 pcm.Then, these data were used as input for the point NK model.


## 5. RELAP-3D Model Used

For the thermal-hydraulic calculations, a RELAP-3D  model has been built, based on the start-up core design [16]. The model is illustrated in Figure 4.Figure 4
RELAP5-3D ETDR nodalization [10].The RELAP5-3D  nodalization of the ETDR reactor can be divided in the following blocks:(i)
vessel and core (components from 100 to 240);(ii)
main cross duct (component 250, 450, 480);(iii)
main heat exchanger (primary side and secondary side components from 290 to 400 and from 700 to 760);(iv)
blower (component 350);(v)
decay heat removal system (components from 800 to 880 and from 900 to 999).It is important to note that secondary system has been modeled as boundary condition.


## 6. Point-NK TH Steady State Calculations with and without MA

Three calculations have been performed:(1)
“imposed power” calculation (performed with no NK data. The 50 MW reactor power is provided with no neutronic feedbacks; so power is constant, and system does not “feel” any changes in reactivity due to neutronic feedback);(2)
reference core composition with NK feedbacks added;(3)
reference core with MA and NK feedbacks added (different fuel composition; see Section4).A comparison between those calculations is shown in Table2.Table 2
Point Neutron Kinetics (NK) core parameters (with and without MA)
ParametersImposed powerReference core (point  NK values)Reference core + MA  (point NK values)Reactor power50 MW50.02 MW50.02 MWTotal system flow rate32.10 Kg/s32.13 Kg/s32.11 Kg/sFlow rate in hot channel0.581 Kg/s0.581 Kg/s0.581 Kg/sFlow rate in average channel0.595 Kg/s0.595 Kg/s0.595 Kg/sReactor pressure at top of vessel6.92 MPa6.924 MPa6.92 MPaReactor mass of coolant without DHR389.75 Kg389.75 Kg389.75 KgReactor mass of coolant with DHR473.54 Kg473.54 Kg473.54 KgCore inlet temperature263.9℃263.1℃262.8℃Core average exit temperature563.7℃564.5℃564.6℃Coolant temp rise in the core299.8℃300.5℃299.3℃Exit coolant temp in hot channel615.6℃616.5℃616.6℃Exit coolant temp in avg. Channel562.7℃563.5℃563.6℃Max clad temp in hot channel687.6℃693.5℃691.7℃Max clad temp in avg. channel621.2℃626.5℃624.9℃Max fuel temp in hot channel1069.5℃1073.9℃1070.4℃Max fuel temp in avg. channel946.0℃949.9℃946.8℃Main blower speed523.6 rad/s523.6 rad/s523.6 rad/sMain blower rated head8007.8m2/s28321m2/s28319m2/s2Main blower rated torque2009.1Pa·m32024.2Pa·m32024.5Pa·m3MHX heat transfer area157.3m2157.3m2157.3m2As it can be seen above, the differences between all the relevant TH parameters are negligible, even with different fuel vectors (→ different neutronic feedbacks, different reactivity insertions).It is interesting to notice that adding the MA into the fuel is not affecting the safety relevant parameters (e.g., cladding and fuel temperatures). Therefore reactor safety-related steady state TH parameters are roughly independent in the range of considered fuel compositions.


## 7.  Point-NK TH LOFA Calculations

In the previous paragraph we showed that the introduction of small quantities of MAs in the fuel will have substantially no effect on the steady state calculations. No particular differences were also found comparing two scram curves, one with standard imposed decay heat [10] and the other one considering the effects of the small fractions of MAs. The difference of relative power (relative error) is never going beyond 1.4%  [10].We present now, in the following figures, LOFA transient analyses performed with the MAs concentration illustrated in Table1, compared with results from the standard imposed scram curve.LOFA analyses were performed using the run-down blower curve provided by CEA [2].The sequence of events leading to the LOFA transient is hereafter illustrated (the steady state condition is achieved after about 50 seconds. Transient sequence is started at second 378 of run (as suggested by benchmark data [16])):(1)
time 0 second: main blower loss of power;(2)
time 0.75 second: reactor scram (relative blower speed going below 90%);(3)
time 173.66 second: closure of main loop valves and opening of DHR loop valves.It should be noted that we perform the MAs analyses referring to a concentration deriving from a typical LWR discharge (see Table1) and to an imposed concentration of  5% of MAs in the fuel composition.Some differences are expected because of the different isotopic composition of compared fuels, leading to different decay head profiles: in detail, higher pressures and temperatures are expected from the presence of higher quantity of MAs in the fuel, especially in the third transient phase (long term cooling), because of the higher activity in the fuel.In the next figures, the LOFA transient evolution is shown. After the scram event, pressure and temperatures in the primary system start to drop, due to the stopped power generation  in fuel. Before the closure of main circuit valves and opening of DHR valves, we can notice (~100 seconds after scram) pressures and temperatures to start rising again, because of mass flow becoming too low.After 173 seconds after scram, blower relative rotational velocity falls below10%, allowing the main valves to close and the DHR circuit to be opened, thus starting the insetting of Natural Circulation (NC).Helium in DHR tubes starts flowing downwards because of lower temperature and higher pressure than core, refraining the NC to incur and causing the raise of pressure and temperature (see Figures 4, 5 and 6).Eventually, after~470 seconds from start of transient, conditions are met for the NC to incur, starting the long term cooling phase.In Figure5, the behavior of pressure at top of the RPV for different types of fuels is shown. The increasing pressure in the RPV is due to greater decay heat of MA isotopes.Figure 5
Pressure at the top of the vessel.Also the rising of core average exit temperature (Figure6) is the consequence of the same cause.Figure 6
Core average exit temperature.The curve in Figure7 shows that the most important parameter of the thermal-fluid-dynamics core design (maximum clading temperature in hot channel) is basically the same for all the analyzed situations.Figure 7
Maximum fuel temperature in hot channel.Figures8 and 9 show the power exchanged in both DHR HXs (primary-secondary side helium-water HX and secondary-final heat sink water-water HX), showing the convergence of the exchanged power.Figure 8
Exchanged power in 1st DHR HX.Figure 9
Exchanged power in 2nd DHR HX.


## 8. Conclusions

From the calculations we can note that there are no relevant differences in both steady state and transient behavior between cases with and without MAs inside the fresh fuel.The transient analyses showed that a limited insertion of MAs (~5%) does not significantly affect transient behavior of a GCFR system. Small differences are noticed after about 400 seconds from the LOFA starting point, when NC occurs. In addition very small pressure differences were identified (0.04 MPa), and they are mostly due to a higher residual pressure in the system and to higher coolant, cladding, and fuel temperatures (about 20 ℃ higher than the reference case).DHR system performed well under all analyzed cases. Generally speaking the capability of removing the relevant amount of heat produced in the core (both in nominal and LOFA conditions) is demonstrated even in presence of (relatively) limited amount of MAs.In conclusion, we can state that the influence of adding limited quantities of MAs is negligible considering the proposed GCFR design.The MA quantity that can be added is currently limited by safety and control problems related to lowerβ values for MA. Therefore further investigations about this topic could clarify the capabilities of GCFR system in increasing the MA input quantities. In that case the (foreseen) introduction of greater quantities of MA into the fuel for maximizing transmutation should be supported by a larger set of neutronic and thermal-hydraulic analyses.



---

*Source: 573481-2009-05-27.xml*