1. Introduction
Compared to large-scale commercial nuclear reactors, small modular reactors (SMRs) have captured widespread global interest as an energy source. This is because they have several advantages including lower total overnight expenses, less infrastructure and staff demands, suitability for smaller power grids, potential to add more capacity and other possible non-electric uses (potable water plants, hydrogen production and district heating) [
1]. A feasibility study on small modular reactors (SMR) carried out by the National Nuclear Laboratory found that SMRs have remarkable potential as an international industry because they meet market requirements that cannot be fulfilled by large nuclear power plants in all conditions. A global market assessment has predicted that the magnitude of the potential market for SMRs is about 65–85 GWe by the year 2035, provided the economics are competitive [
2]. The SMR design concepts that have been proposed can be classified into numerous types such as the pressurized water reactor (PWR) [
3,
4,
5], the boiling water reactor (BWR) [
6], the liquid metal cooled fast reactor (LMFR) [
7,
8,
9,
10,
11,
12,
13], the gas cooled fast reactor (GFR) [
14,
15] and the molten salt reactor (MSR) technologies [
16,
17]. Among these concepts, the liquid metal type lead or lead-bismuth eutectic (LBE) is generally chosen as the coolant for SMRs owing to its superior properties including a low melting point, high boiling point and chemical inertness with air and water. The LBE has better shielding properties against gamma rays and due to its high boiling point (1670 °C) there is less chance of a major void forming in the core due to overheating of the coolant. The LBE has fine thermal conductivity and due to its low kinetic viscosity and substantial expansion coefficient, it is more suitable for natural circulation. However, aside from these promising features, there are some limitations associated with the LBE coolant. One of the key disadvantages of using LBE is the erosion/corrosion it causes due to its high solubility with important structural and cladding materials at high temperature. This problem can be significantly mitigated by the selection of suitable structural materials (such as T91), oxygen control techniques and operating the reactor in a safe temperature range [
18,
19,
20,
21]. Besides the corrosion/erosion, the radiation-induced deterioration of mechanical properties and embrittlement of liquid LBE causes a significant decrease in the performance of the key structural materials. Another disadvantage of LBE is the formation of alpha emitter polonium-210 isotope by neutron capture of bismuth. It is possible that this radioactive polonium could be released into the air and inhaled by workers during maintenance and LBE leakage accidents. Protection measures need to be implemented to avoid personal irradiation when polonium-210 is above the tolerable limit during normal operation as well as during repair works, refueling, and emergency discharge of the coolant [
22]. Additional drawbacks of LBE coolant include the high material cost of bismuth, high toxicity of lead, huge static pressure on the bottom of pressure vessel because of its high density, high pumping demand for the coolant, and greater seismic isolation conditions due to high density. The LBE has a melting temperature of 125 °C and its liquid state is maintained with the help of a steam generator. The temperature of the water supplied to the steam generator is kept higher than the melting point of LBE. The steam supply system or electrical heating can be utilized to heat the reactor during startup or maintain the primary system in a hot condition at the lowest power level during normal operation [
23]. Another important problem related to the operation of LBE fast reactors is the chemical control of the coolant, which includes oxygen control and pollution source term analysis. The concentration of certain impurities is likely to activate and have an impact on erosion/corrosion and scale formation on the heat transfer surfaces. To control the chemistry, the following measures should be considered: quality control and purification of the coolant during operation and oxygen control, assessment of radiotoxicity of different components and removal, examination of core integrity, and maintenance.
Lead-bismuth eutectic coolant technology with advanced safety characteristics was first mastered in the former Soviet Union. The LBE coolant was chosen by the Russians for alpha class submarine reactors, which resulted in the widespread research and development of liquid metal coolant and related material technology [
24]. Since the 1990s, Russia has taken renewed interest in LBE coolant technology for civilian nuclear power reactors. Recently, the Russian Federation developed a modular-type innovative fast reactor design, SVBR-100 with LBE coolant [
25]. It employs a closed fuel cycle with mixed uranium and plutonium oxide fuel and the core design life is 60 years. In the early 21st century, conceptual designs for several fast spectrum secure transportable autonomous reactors (STAR) were presented by the Argonne National Laboratory. Three designs concepts were included: the encapsulated nuclear heat source (ENHS) [
26], STAR-LM [
27], and STAR-H
2 [
8]. The conceptual design of the ENHS reactor includes a LBE-cooled fast reactor of 125 MWth with 100% natural circulation, autonomous operation and a long core life of over 20 effective full power years (EFPY) of operation without refueling [
26]. STAR-LM is a 300–400 MWth LBE-cooled fast modular reactor that exhibits 100% natural circulation and uses nitride fuel cartridges, and generates electricity with almost autonomous operation for 15 years [
27]. STAR-H
2 is a 400 MWth, fast neutron spectrum, modular reactor concept with lead as a coolant and it features natural circulation and passive load following capabilities for hydrogen production [
11]. Recently, a conceptual design for a small STAR (SSTAR) was introduced, which integrates a supercritical CO
2 (S-CO
2) Brayton cycle energy conversion system with a small lead-cooled fast reactor (LFR) to attain a more compact reactor design. SSTAR is a pool-type LFR conception with 45 MWth power that can be operated for 30 years without refueling and it is capable of load following without control rod motion [
10]. In 2011, an upgraded 300 MWth fast spectrum, pool-type and lead-cooled SMR design named SUPERSTAR was presented at the Argonne National Laboratory (ANL). SUPERSTAR incorporates a 100% natural circulation primary coolant system, S-CO
2 Brayton cycle power converter and has a core lifetime of 15 years [
7]. In recent research, a 100 MWth lead-cooled small modular fast reactor SNCLFR-100 was proposed by the University of Science and Technology of China (USTC), which has a refueling interval of 10 years without core rearrangement and natural circulation during both normal and abnormal operations [
28]. Most modern nuclear power plant concepts incorporate natural circulation and load following features for design simplification and improved maneuverability, respectively.
Two conceptual designs of small modular portable reactors, SPARK and SPARK-LS, which are characterized by forced circulation of lead bismuth and chloride salt, respectively, were proposed in a previous study [
12,
17]. The study was mainly focused on developing a compact and small core design that can be conveniently transported by truck or barge. Due to the strong erosion/corrosion effect on materials and the high density of LBE, the long-term operation of these reactors installed in far-off areas demands long-run, safe operation of the primary coolant pumps. Up until now, this has proven to be a challenging task for the design and maintenance of the liquid metal pumps. Therefore, the preferred option for inherent safety and design simplification is to construct the primary coolant system based on natural circulation by removing the intermediate heat transport system.
In the present study, a LBE cooled fast reactor, named SPARK-NC, with advanced features of natural circulation and load following was proposed for remote power supply. The optimal rated power output was set as 10 MWe, and the power generated can be regulated from 5 MWe to 10 MWe during the reactor lifetime. The steam Rankine cycle
was selected as the power conversion system on account of its technology maturity. The scope of this paper does not cover an analysis of the steam Rankine cycle for the energy conversion system. The design constraint of transportation in the previous designs (SPARK and SPARK-LS) was alleviated to achieve natural circulation and improve the safety of the heat removal system. SPARK-NC can produce electricity for eight effective full power years (EFPYs) without refueling, and it can fulfil the power demands of a small and remote town with a population of approximately 20,000 [
29]. In off-grid and remote places, the diesel generators deployed for power supply have various disadvantages such as their use of extremely high-priced, inaccessible fuel, recurring maintenance, low efficiency, noise and air pollution. The SPARK-NC is a viable power source for remote locations compared with conventional diesel generators because of its better economics, simplified design, improved maneuverability, and lower operational and maintenance costs. The motivation to utilize nuclear reactors with the capability of load following arises from the extensive deployment of alternate sources of electricity such as wind and solar power. Due to various factors, alternate sources of electricity may not be constantly available, for example, solar power does not produce electricity at night or during bad weather conditions and differs from summer to winter. A sophisticated system is needed to ensure that electricity generation corresponds to varying power demands. Although most power reactors are designed to operate with base loads nowadays, in order to have the means to adjust power supply to daily or occasional fluctuations in electricity demand, changes are needed to enhance the maneuverability of power reactors.
The neutronics and thermal-hydraulics behavior of the optimized core were studied. The kinetics parameters and reactivity feedback coefficients were evaluated at different time stages in the fuel cycle, and the reactivity worth of the control and scram systems was calculated. The natural circulation and load following capabilities of SPARK-NC were considered systematically using energy conservation, the pressure drop and the quasi-static reactivity balance (QSRB) equations [
30]. Three possible methods for the regulation of the reactor power output were examined based on the systematical analysis of the load following and natural circulation capabilities.
The remainder of this paper is ordered as follows.
Section 2 presents the design goals and constraints.
Section 3 introduces the computational approaches and models employed in the design study.
Section 4 describes the optimized core design and performance.
Section 5 provides the natural circulation analysis,
Section 6 includes load following analysis, and the conclusions are presented in
Section 7.
5. Natural Circulation Analysis
A simple natural circulation study was carried out for the SPARK-NC by finding the height difference between the thermal and cold centers of the primary LBE coolant system that is necessary to remove the heat produced in the core and sustain natural circulation under a nominal power of 30 MWth. The thermal center is the center of the active reactor core while the cold center represents the center of the steam generator, which works as the heat sink in the primary system. For reactor power
and coolant density
, if
and
characterize the fuel assembly number and number of fuel pins per fuel assembly, respectively, we can write the energy conservation equation for the core as
where
represents the coolant velocity,
is the coolant temperature rise and
is the specific heat capacity of the coolant. The total pressure loss for LBE in the primary coolant system can be estimated as
where
is the total length of the fuel pin,
denotes the equivalent hydraulic diameter,
represents the Darcy–Weisbach friction factor, and M is the multiplication factor that expresses the pressure losses apart from the frictional losses of the fuel pins. In this case,
was applied according to previous research produced by the Massachusetts Institute of Technology [
41]. The expression for buoyancy pressure due to the gradient of coolant density can be written as
where
is the volumetric expansion coefficient,
and
is the height difference between the thermal and cold centers as shown in
Figure 9.
In natural circulation, only one of three variables
,
and
is independent. For a given reactor power, the coolant velocity
and rise in temperature
can be estimated by using the energy and pressure balance equations. The coolant velocity can be eliminated from the pressure drop equation by employing the energy conservation equation as follows:
By setting the total pressure drop equal to the buoyancy pressure, the velocity due to the natural circulation
and coolant temperature rise
are obtained as
and
On the other hand, with a predetermined coolant temperature rise, the height difference between the thermal and cold centers for natural circulation can be obtained as follows:
Assuming that SPARK-NC has a 10 MWe power output and a predetermined coolant temperature rise of 120 °C, the corresponding height difference between the thermal and cold centers that is essential for heat removal from the core was estimated to be 1.42 m.
6. Load Following Analysis
There are three possible ways to regulate the reactor power to the target level: (1) by keeping the control assemblies fixed and changing the coolant inlet temperature, (2) moving the control assemblies and keeping the coolant inlet temperature constant, and (3) adjusting both the coolant inlet temperature and the position of the control assemblies. Since method (3) is complicated and makes it difficult to control the reactor, methods (1) and (2) were examined in this study. Both methods are discussed in detail in the subsequent sections.
6.1. Adjusting the Coolant Inlet Temperature
In this method for load following analysis of SPARK-NC, the quasi-static reactivity balance (QSRB) method established by Wade et al. at ANL [
30] was implemented to calculate the required externally-imposed reactivity (by moving the control rods or adjustment of coolant inlet temperature) to change the reactor power to the desired level. The reactivity balance equation in the QSRB method is written as:
where
and
represents the normalized power and normalized flow, respectively,
is the externally-imposed reactivity and
is the change from typical coolant inlet temperature.
,
, and
are integral reactivity parameters, which are calculated from the reactivity feedback coefficients in Equations (9)–(11),
is the rise coolant temperature, and
is the increase in the average fuel temperature relative to the average coolant temperature.
In the case of slow transients, in which the non-equilibrium energy in the fuel and delayed neutron imbalance can be prevented, Equation (8) illustrates the relationship among the externally controlled flow and temperature parameters of the coolant, the externally imposed reactivity, and the new adjusted power level of the reactor core. The power level ascends or descends through the core power coefficients in order to compensate for any change in the reactivity. The reactivity balance equation can be correlated with the energy conservation and pressure drop equation to estimate externally-imposed reactivity (through the adjustment of the coolant inlet temperature or by movement of the control assemblies) required to attain the desired power level. The nominal power output is set to 10 MWe as the basis of the operation and the essential externally-imposed reactivity for different power outputs at various time points of the fuel cycle is estimated using method (1), as shown in
Figure 10. The inlet coolant temperature and the rise in coolant temperature at various power outputs and at different time points in the fuel cycle are plotted in
Figure 11.
It is clear that the power output can be changed from 10 MWe to 5 MWe at BOC, MOC or EOC by raising the coolant inlet temperature by around 37 °C, while the corresponding coolant outlet temperature is approximately 445 °C. When the power output is changed from 5 MWe to 9 MWe, the range in the variation of the coolant inlet and outlet temperatures overlaps with the range of the nominal coolant inlet and outlet temperatures in the basic operational case (10 MWe) throughout the entire fuel cycle, which indicates the viability of load following operation by adjusting the coolant inlet temperature.
6.2. Adjusting the Position of Control Assemblies
In the second method, manipulation of the power was achieved by moving the secondary control assemblies, which act as power regulating assemblies while keeping the inlet temperature constant. The reactivity insertion during load following and the resulting effect on the power distribution and thermal hydraulic parameters were investigated using DAISY code. Control assemblies with relatively low rod worth were chosen to avoid deformation in the power shape and degrading the neutron economy. The differential reactivity worth of secondary control assemblies at BOC is given in
Figure 12.
The change in the power distribution with the insertion of power regulating rods was studied at BOC. The radial power distribution was not deformed significantly and the peaking factor of the hottest assembly changed from 1.33 to 1.34 at BOC as shown in
Figure 13. The axial power profile in central and peripheral regions changed slightly when the low worth control rods were moved. The axial peaking factors of the lower half of the fuel assemblies bulged out slightly as shown in
Figure 14.
To lower the reactor power level by 50%,
$0.08 negative reactivity was inserted for 15 sec as shown in
Figure 15. The negative reactivity coefficients start working and begin to raise the reactivity until it reaches a steady state. while the power level reaches its desired level in 5 min and remains unchanged due to the reactivity coefficients and neutronics parameters. The power level was also reduced in 10% steps between the rated power of 10 MWe and 5 MWe by inserting different reactivities with the regulating assemblies as shown in
Figure 16. The temperatures are proportional to the power, when the power is decreased the temperature of the fuel and coolant also decreases. The change in temperature of fuel and coolant corresponding to a 50% drop in reactor power is shown in
Figure 17.
The initial values of the coolant outlet temperature, maximum cladding surface temperature, fuel maximum temperature and fuel average temperature were 450 °C, 484.4 °C, 556.5 °C and 441.7 °C, respectively. After the regulating rods were inserted, these values started to decrease and reached constant values of 390.7 °C, 407.8 °C, 441.1 °C, 385.7 °C, respectively. In fast spectrum reactors, the reactivity feedback coefficients are very small. Therefore, for the desired power output range, the adjustment range of the externally-imposed reactivity is quite narrow. This presents a challenge in fine-tuning the position of the control assemblies to provide small-scale reactivity. The adjustment of coolant inlet temperature is a reasonably easy task that involves simply altering the mass flow rate of water into the steam generator. Hence, method (1), which attains the desired power level through regulating the coolant inlet temperature is preferred.
7. Conclusions
In the present study, a conceptual design for a small lead-bismuth-cooled core fast reactor, SPARK-NC with natural circulation and load following was proposed for the purpose of providing remote power supply. The nominal power output was set to 10 MWe, and can be adjusted from 5 MWe to 10 MWe throughout the entire core lifetime. The core neutronics and thermal-hydraulics design calculations were performed using the SARAX code and the natural circulation capability of SPARK-NC was investigated by employing the energy conservation equation, pressure drop equation and quasi-static reactivity balance equation. In order to flatten the radial power distribution, three radial zones were constructed by using different fuel enrichments and fuel pin diameters. To ensure an adequate shutdown margin, two independent systems, i.e., a control system and a scram system, were introduced in the core. The control assemblies are further classified into two types: primary control assemblies used for reactivity control and power flattening and secondary control assemblies (with relatively smaller reactivity worth) used for power regulation. The tungsten ballast coated with B4C was used in the scram assemblies to counteract the buoyancy of the high-density coolant. The targeted fuel cycle length of 8 EFPYs was fulfilled with a compact core design, and the neutronics and thermal-hydraulics parameters were below the design constraints throughout the entire fuel cycle. Two methods for the load following capability (adjusting the coolant inlet temperature and the movement of secondary control assemblies) of SPARK-NC were assessed. By comparing various methods for adjusting the reactor power output, we found that adjusting the coolant inlet temperature was the most viable, easy to implement and favored for the load following operation.
The reliability of fuel and structural components in load following reactors is very important for maneuverability. The alteration in power with fluctuating demand produces variations in the coolant temperature, which consequently leads to periodic changes in the mechanical load of some reactor components. For large temperature gradient, this may create localized structural impairment (fatigue) in the reactor components. Moreover, the fuel performance could be constrained by the deterioration in the fuel cladding due to the pellet-cladding interaction (PCI), stress corrosion cracking and various other processes.
In future work, the mechanical characteristics and integrity of fuel and structural components need to be assessed during the load following operation of SPARK-NC. This examination should include the fatigue, tensile and fractural toughness of fuel and structural parts at high temperature, and also the study of long-term corrosion and irradiation effects.