Core and blanket thermal–hydraulic analysis of a molten salt fast reactor based on coupling of OpenMC and OpenFOAM

Bin Deng· Yong Cui· Jin-Gen Chen· Long He·

Shao-Peng Xia1,2· Cheng-Gang Yu1,2· Fan Zhu1,2,3 · Xiang-Zhou Cai1,2,3

Abstract In the core of a molten salt fast reactor(MSFR),heavy metal fuel and fission products can be dissolved in a molten fluoride salt to form a eutectic mixture that acts as both fuel and coolant. Fission energy is released from the fuel salt and transferred to the second loop by fuel salt circulation.Therefore,the MSFR is characterized by strong interaction between the neutronics and the thermal hydraulics. Moreover, recirculation flow occurs, and nuclear heat is accumulated near the fertile blanket, which significantly affects both the flow and the temperature fields in the core. In this work, to further optimize the conceptual geometric design of the MSFR, three geometries of the core and fertile blanket are proposed, and the thermal–hydraulic characteristics, including the three-dimensional flow and temperature fields of the fuel and fertile salts, are simulated and analyzed using a coupling scheme between the open source codes OpenMC and OpenFOAM. The numerical results indicate that a flatter core temperature distribution can be obtained and the hot spot and flow stagnation zones that appear in the upper and lower parts of the core center near the reflector can be eliminated by curving both the top and bottom walls of the core.Moreover,eight cooling loops with a total flow rate of 0.0555 m3 s-1 ensure an acceptable temperature distribution in the fertile blanket.

Keywords Molten salt fast reactor · Core and blanket thermal–hydraulic analysis · Neutronics and thermal hydraulics coupling

1 Introduction

The molten salt reactor(MSR)was selected as one of six advanced reactor concepts by the Generation IV International Forum in 2002 because of its considerable advantages of inherent safety, sustainability, excellent neutron economy, and nuclear nonproliferation [1]. Most MSR concepts were designed as thermal systems in the early stage of research and development. In 1954, the first MSR in the world,the Aircraft Reactor Experiment[2],was built at Oak Ridge National Laboratory; it had a power of 2.5 MW.In the 1960s,the Molten Salt Reactor Experiment[3], with a power of 8 MW, achieved criticality, further demonstrating the feasibility of operating a liquid-fueled reactor. After decades of MSR development, renewed interest in MSRs focuses on advanced goals by developing reactor and fuel cycle technologies [4–7]. The sixth and seventh EURATOM Framework Programs (the ALISIA and EVOL projects) proposed the concepts of the molten salt actinide recycler and transmuter[8]and the molten salt fast reactor [9]. The former was developed to burn transuranium elements (TRU) and minor actinides,whereas the latter is used mainly for breeding thorium fuel.Compared to traditional graphite-moderated MSRs, the MSFR not only has better232Th–233U breeding and TRU burning capability, but is also inherently safer, because it can avoid the positive reactivity feedback caused by a graphite moderator [10].

The innovative MSFR was initially designed as a singlecylinder core with a radial fertile blanket and absorber(the benchmark model of the MSFR) [9]. Thermal–hydraulic(TH) analyses of the MSFR benchmark (see Figs. S1 and S2) show that a large recirculation flow of fuel salt is formed near the fertile blanket,causing a hot spot to appear in this region and thus degrading the heat transfer performance of the reactor core [11, 12]. The corrosion rate of structural materials in the core is known to increase rapidly with increasing operating temperature. Thus, the core model should be further optimized to eliminate the hot spot and improve the core lifetime of the MSFR. For the channel-type MSR, the flow rate distribution is critical for optimizing the core temperature field,because the location and peak value of local hot spots are determined by the fuel flow in the graphite channels [13, 14]. Furthermore, TH studies of the MSFR indicate that it is essential to obtain a relatively flat flow field. Rubiolo et al. analyzed the core thermal hydraulics of the MSFR under nominal conditions[15],and Laureau et al.performed a steady-state neutronics and TH coupling simulation of an MSFR[16].Both groups adopted an optimized MSFR model in which curves were added to the radial walls of the core cavity in the benchmark model. Although the recirculation flow reported in these studies is slightly improved, a hot spot still exists on the wall near the fertile blanket.That is,the core geometry could be further optimized. Rouch et al. proposed a preliminary core TH design for an MSFR with the main objective of optimizing the core cavity walls, which can increase the overall flow mixing and reduce the temperature peak in the fuel salt and on the core walls [17]. The heat source used in the computational fluid dynamics(CFD) simulation was obtained from MCNP, whereas two nonidentical geometrical models were adopted in the CFD and MCNP simulation, which may make the coupling result unreliable. Further, temperature differences still appear even in a later study in which a geometrical homothetic transformation [17]was employed to adapt the heat source to the real geometry adopted for TH analysis.Laureau et al. further optimized the core geometry of the MSFR,and the hot spot on the curved wall near the fertile blanket was eliminated [18]. However, a new hot spot at 1080 K appeared at the top of the reactor core owing to the relatively low velocity in this region;the hot spot exceeded the upper temperature tolerance (1073 K) of the structural material [19]. For the safety assessment for MSFR, evaluating the THs behavior of the fertile blanket is also of great importance.The small fission fraction of233U and the neutron capture of232Th in the blanket heats the fertile salt,and the heat cannot be removed by heat convection between the fuel salt and the interface wall of the core and the fertile blanket [9, 19, 20]. Therefore, an external cooling system for the fertile blanket has to be considered in the reactor design.

In this work, a new coupling scheme between the open source codes OpenMC and OpenFOAM is developed to physically describe the strong coupling between the neutronics and thermal hydraulics in the MSFR. The distributions of the core power, fuel temperature, and fuel density are exchanged between the two codes to perform the coupling calculation, with the goal of eliminating the hot spot and flow stagnation zones in the core.In addition,a preliminary TH analysis of the fertile blanket is also performed. As mentioned above, the heat deposited in the fertile blanket cannot be evacuated by heat convection,which indicates that the heat transfer between the fuel salt and fertile blanket is relatively small. In addition, the influence of the core heat on the fertile blanket is analyzed in this work by using a conjugate heat transfer calculation for the core and fertile blanket. The results demonstrate that separate modeling of the fertile blanket is reasonable for MSFR TH analysis.

This paper is organized as follows. Section 2 describes the methodology and tools.Section 3 presents the core and blanket TH analyses. The conclusions of this work are given in Sect. 4.

2 Methodology and tools

2.1 Neutronics model

OpenMC is an open source Monte Carlo particle transport simulation code that can perform fixed source, keigenvalue, and subcritical multiplication calculations on three-dimensional solid geometry models [21]. It was originally developed by the Massachusetts Institute of Technology and is widely used for calculations of the criticality, depletion, multigroup cross section generation,etc.The Direct Accelerated Geometry technology has been integrated into OpenMC (DAG-OpenMC hereafter) to perform neutron transport simulation directly on a computer-aided design model [22]. In this work, this technology contributes to the solution of the problem of curved wall modeling for an MSFR owing to its advantages for precise modeling of an arbitrary complex geometry.

Based on the ENDF/B-VII.1 nuclear data,a microscopic cross section library between 500 and 1600 K with 100 K internals was produced by NJOY2016 to support the OpenMC continuous-energy transport calculation; a 20 K interval was used between 800 and 1400 K (the operating temperature range of MSFRs).The cross sections are given exactly for only a few points in the table; the temperature values for other materials can be obtained by statistical interpolation[23].Suppose that cross sections are available at temperatures T1, T2, …, Tn, and a material is assigned a temperature T, where Ti＜T ＜Ti+1. Statistical interpolation is applied as follows. A uniformly distributed random number (ξ) of unit intervals (0–1) is sampled. If ξ ＜(T - Ti)/(Ti+1- Ti),the cross sections at temperature Ti+1are used. Otherwise, the cross sections at Tiare used.To obtain the power distribution of the MSFR, the calculation model is divided into numerous regions, and each region is labeled with a unique material number.

2.2 Thermal–hydraulic model

The open source code OpenFOAM features automatic matrix construction and solution capabilities for scalar and vector equations [24]. Each fluid flow can be represented by a series of partial differential equations; the solvers for these fluid mechanics equations have been developed and implemented in OpenFOAM.

The fuel salt flow in the MSFR core without a moderator is considered to be turbulent [25], because the Reynolds number is estimated to be 2.5 × 106in the core region and 5.3 × 105in the out-of-core region [11]. The fuel and fertile salt flows can be assumed to be incompressible,because the effect of fuel density variation in the gravity force term is taken into account by using the Boussinesq approximation. Moreover, the ReNormalization Group(RNG) k–ε turbulence model can provide a better description of the flow on the curved wall than the standard k–ε turbulence model [26]. Therefore, considering that the fuel salt flow is incompressible, the Reynolds-averaged Navier–Stokes equations with the Boussinesq closure hypothesis and the RNG k–ε turbulence model are adopted,and buoyantBoussinesqSimpleFoam [24] is selected as the OpenFOAM solver. The steady-state continuity, momentum, and energy conservation equations for the fuel and fertile salts can be written as

2.3 Coupling scheme

The coupling calculation flowchart of DAG-OpenMC and OpenFOAM is presented in Fig. 1, illustrating the exchange of core power distribution, temperature, and density distribution of the fuel and fertile salt using a Python script (pyMF). As each region of the neutronics calculation model contains many meshes for the TH calculation by OpenFOAM, the average fuel temperature and density of each region are assigned according to the material number (where each region has a unique material number) used in DAG-OpenMC. Similarly, the Open-FOAM meshes located in the same region are assigned the same power density from the neutronics calculation.

The key points of the coupling code pyMF are summarized as follows:

(1) The temperature and density of the fuel salt and structural material are first assumed to be uniform in the DAG-OpenMC input file to start the neutronics calculation.

(4) The above iteration is performed until either keffand the power density converge or the maximum number of iterations is reached [27–29].

In our previous work[30],to verify the reliability of the coupled code, a benchmark model of the MSFR was established, and the effects of the region number in the neutronics calculation and the initial conditions on the keff,fuel salt velocity, and temperature distributions were studied. The TH results simulated using pyMF were also compared to the reference result for the MSFR [31], and good agreement was obtained for the velocity and temperature fields (Table S1 and Figs. S1 and S2). It can be concluded that the pyMF code can provide reliable results for the steady-state neutronics and TH coupling calculations for an MSFR.

3 Results and discussion

3.1 Description of the molten salt fast reactor

A schematic view of the MSFR system is shown in Fig. 2.The MSFR is a 3 GWth liquid-fueled reactor with a core fuel salt volume of 9 m3[9, 19, 20]. In the primary loop,a molten fluoride salt system functioning as both fuel and coolant works at approximately 923 and 1023 K at the inlet and outlet of the core, respectively. A Hastelloy[19,20]is selected as the axial reflector,and a radial fertile blanket is adopted to absorb the neutrons escaping from the core and thus improve the breeding capability of the reactor. The primary loop of the MSFR has 16 fuel salt loops, and each loop is equipped with a pump and a heat exchanger.The fuel salt circulates through the core and the 16 external loops,and the nuclear heat generated by fission is transferred to the secondary loop salt by the heat exchangers.The main parameters of the MSFR reactor core are listed in Table 1 [12, 20, 31].

3.2 Core thermal–hydraulic analysis

As illustrated in Fig. 3, three core geometries are proposed in this work to analyze the TH behavior of the MSFR. Geometry I (Fig. 3a) is described in the literature[17, 18]. Geometry II (Fig. 3b) is obtained by further curving the upper wall of the core, and Geometry III(Fig. 3c) is obtained by curving both the upper and lower walls of the core for comparison.All the arcs of the blanket and absorber in the three geometries are half-ellipses.

In the OpenFOAM simulation, only the core, fertile blanket, and primary loop are modeled, and the solid regions, including the reflector and absorber, are not considered,because the coupling calculation is performed only in the above three regions. To obtain a good balance between the calculation efficiency and accuracy, the maximum bulk cell size is set to 3 cm [16, 17]. The RNG k–ε turbulence model and a wall function are adopted to model the flow behavior near the wall. In fluid dynamics, a dimensionless distance to the wall (y+) is used as an evaluation criterion of near-wall mesh quality.Therefore,a suitable y+value is needed to obtain accurate results when the k–ε turbulence model and wall function are employed.To ensure that the y+values of all the volume elements adjacent to the wall are distributed within a range of 11.63–300 [32], the cell size is reduced to approximately 3 mm near the core wall. Considering that the MSFR core configuration has good symmetry, only 1/16 of the core geometry (Fig. 4) is simulated in the TH calculation to conserve computational resources. The pump and heat exchanger are treated as a momentum source in the momentum conservation equation and as a heat sink in the energy conservation equation. The mesh numbers of these three geometries are 124,540 (Fig. 4a), 132,784 (Fig. 4b),and 138,562 (Fig. 4c), respectively.

A continuous-energy neutronics simulation is carried out using the DAG-OpenMC code; the calculation models(Fig. 5) include the core, fertile blanket,and solid regions.All of the components are simulated during the neutron transport simulation; a reflective boundary condition is set for the symmetry planes,and a vacuum boundary condition is set for the walls. No heat exchange is modeled between the core and the fertile blanket or between the fertile blanket and the solid region, and the temperatures of the fertile blanket and solid region are assigned as the average temperature of the fuel salt in the core TH analysis. To obtain the power distribution for the TH calculation, the core is divided into numerous regions, and each region is labeled with a unique material number during the neutronics calculation. Considering the efficiency and accuracy of the neutronics calculation,the size of these regions should be approximately 20 cm, as demonstrated in our previous work [30]. Only a small fraction of the power is released in the external loops, and the strong coupling between the neutronics and thermal hydraulics is significant in the core region. Therefore, no detailed regional division of the external loops is performed.For the neutron transport simulation,300 generations(the first 50 iterations are skipped) with 50,000 neutrons per generation are adopted, and the obtained standard deviation of keffis approximately 0.00023.

The fuel velocity fields of the three geometries are shown in Fig. 6, which shows that a low flow velocity occurs near the outer wall (close to the fertile blanket),whereas the large recirculation region near the fertile blanket in the reference MSFR model (see Fig. S1) is eliminated, and a relatively flat velocity distribution is obtained for all three geometries. Although a small vortex or recirculation pattern still occurs at the core inlet, no hot spot forms owing to the low power density in the area of interest and the low fuel salt temperature at the core inlet.Therefore, further optimization to eliminate this small vortex and recirculation pattern is not needed. The connecting pipe is much smaller than the core; thus, the flow area is widened or narrowed when the fuel salt is flowing into or out of the core.Furthermore,the connecting pipe is curved. Therefore, a large local pressure drop is obtained,which maximizes the flow velocity in the connecting pipe in the three geometries. The maximum fuel salt velocity obtained for Geometry I (Fig. 6a) is approximately 5.05 m s-1, which is slightly higher than those for Geometry II (Fig. 6b) and Geometry III (Fig. 6c). The reason is that the connecting pipe in Geometry I is smaller than that in the other two geometries.As shown in Fig. 6a,a small stagnation zone appears near the upper and lower reflectors in the velocity field of Geometry I, whereas a stagnation zone appears in Geometry II only at the core bottom, as shown in Fig. 6b. For Geometry III, the stagnation zone is eliminated (Fig. 6c), which indicates that Geometry III allows a more uniform velocity distribution.

Table 1 Main parameters of MSFR

The fuel salt temperature field in each geometry is presented in Fig. 7. The fuel salt temperature increases along the axial direction from the bottom to the top for all three geometries,and a small hot spot appears at the top of the central core in Geometry I (Fig. 7a), which might be caused by the low velocity in this region. The temperature profile for Geometry II (Fig. 7b) does not differ significantly from that for Geometry I,although the peak value of the fuel temperature is slightly lower owing to the difference in the geometry. Therefore, the hot spot at the upper core cannot be effectively eliminated by curving only the upper wall of the core. Unlike these two geometries,Geometry III (Fig. 7c) exhibits a more uniform temperature distribution throughout the core and a relatively small radial gradient of the fuel temperature.Moreover,the peak temperature drops sharply, and the difference between the maximum core temperature and average outlet temperature is only approximately 9 K.

The turbulent dynamic viscosity field of each geometry is presented in Fig. 8. These results are also compared to the result for the reference MSFR model(see Fig.S3).The location of maximum turbulent dynamic viscosity is shifted from the upper core cavity to the lower core cavity. In addition, a zone of high turbulent dynamic viscosity is formed near the entire arc wall of the fertile blanket.Therefore,there is a strong turbulent transport effect in this region,which prevents the formation of the hot spot on the arc wall near the fertile blanket even though the core flow velocity is low in this region (see Figs. 6, 7). The distributions of the turbulent dynamic viscosity for Geometry I(Fig. 8a)and Geometry II(Fig. 8b)are both relatively flat;hence, the heat transfer in the core due to turbulent transport is almost uniform. By contrast, the turbulent dynamic viscosity in the core center in Geometry III(Fig. 8c)is low,and it increases gradually along the radial direction;therefore, the turbulent transport effect also increases gradually. As shown in Fig. 6, the mean flow velocity in the core center is higher than that near the fertile blanket;therefore,the cooling capacity resulting from fuel salt flow in the core center is better than that near the fertile blanket.Owing to the competing effects of fuel salt flow and turbulent transport, Geometry III has the smallest radial temperature gradient in the core and the most uniform radial temperature distribution among the three optimized geometries.

Table 2 lists the average core inlet and outlet temperatures, the temperature difference across the core, the core volumes, the total volume flow rate in the primary loop,and the maximum core temperature of each geometry.The heat exchanges are modeled as suitable heat sinks in the energy conservation equations to obtain the nominal value(923 K) of the core inlet temperature. A different momentum source is added to the momentum conservation equation for each geometry to maintain the temperature difference across the core at the nominal level (100 K),because the core volume and shape of each geometry is different. The temperature difference across the core and the average core inlet and outlet temperatures for Geometry I, II, and III are all approximately 100, 923, and 1023 K,respectively. However, the maximum temperatures differ significantly. The maximum fuel salt temperature in Geometry III is 1032.5 K, which is much lower than those of Geometry I (1070.3 K) and Geometry III (1062.1 K).Geometry III exhibits the best TH performance.

Table 2 Calculation results for the three geometries

3.3 Preliminary thermal–hydraulic analysis of the fertile blanket

The fertile blanket has a much smaller heat source than the core; hence, it is neglected in the TH studies of the core. However, the effect of the core heat on the fertile blanket is unknown, and the feasibility of modeling the fertile blanket separately should be analyzed. In this work,a conjugate heat transfer model of the core and fertile blanket is established; it is shown schematically in Fig. 9.The calculation model presented in Fig. 5c is adopted for the neutronics simulation, and the corresponding TH calculation mesh is shown in Fig. 10.A 3 GWth heat source is loaded in the core region,and three different power loading cases for the fertile blanket are analyzed. An adiabatic boundary condition is set for the blanket wall. In the OpenFOAM code,a heat flux balance is maintained on the interface wall between the core and the fertile blanket,which can be described as

where q is the heat flux; TBis the temperature of the interface wall;and λcand λbare the thermal conductivity of the fuel salt and fertile salt, respectively. Tcand Tbare the core and fertile blanket cell center temperature of the first control volume on the interface wall, respectively; δcand δbare the distance of the interface wall from the core and the fertile blanket cell center, respectively.

The temperature distribution including the conjugate heat transfer of the fuel and fertile salts is illustrated in Fig. 11. A total heat source of 0 MWth (Fig. 11a), 0.1 MWth (Fig. 11b), and 1 MWth (Fig. 11c) is loaded in the fertile blanket. Figure 11a shows that the temperature of the fertile blanket matches that of the core in the axial direction even when no heat source is loaded in the fertile blanket, which reveals that the nuclear heat of the core affects the temperature of the fertile blanket.Figure 11b,c indicates that the nuclear heat generated in the fertile blanket cannot be removed by the fuel salt flow through the interface wall between the core and the fertile blanket.Therefore,an external cooling system for the fertile blanket should be considered.233U will be generated in the fertile blanket, because the thorium loaded in this region is exposed to the core neutron flux, and the heat source (including fission and neutron capture by thorium) generated in the fertile salt is estimated to be approximately 37 MW[19, 20]. However, the temperature of the fertile blanket greatly exceeds that of the core when a heat source of only 1 MWth is loaded in the fertile blanket (Fig. 11c). Therefore,from Fig. 11a,c,one can deduce that the effect of the core heat on the fertile blanket is much smaller than that of the heat source in the fertile blanket itself, which demonstrates that modeling the fertile blanket separately is reasonable for the MSFR TH analysis.

Three cases are designed to study the required flow rate for obtaining an acceptable temperature distribution in the fertile blanket; the external circuit is ignored for convenience. A sketch map of the three cases is presented in Fig. 12 [case 1 (Fig. 12a), case 2 (Fig. 12b) and case 3(Fig. 12c)]. To reduce the computational demands, a minimal symmetrical unit is adopted for the TH calculation; the CFD meshes for the three cases are shown in Fig. 13. According to the grid independence verification,the mesh numbers for case 1 (Fig. 13a), case 2 (Fig. 13b),and case 3 (Fig. 13c) are 108,912, 56,576, and 29,772,respectively.

The core region of Geometry III (Fig. 3c) is selected to consider the effect of the core neutronics on the neutronics behavior of the fertile blanket. Consequently, the fertile blanket is also divided into several regions to obtain its power distribution,as illustrated in Fig. 14.All of the parts are simulated during the neutron transport simulation; a reflective boundary condition is set for the symmetry planes, and a vacuum boundary condition is set for the walls. The heat source of the fertile blanket is estimated to be approximately 37 MW, as mentioned above. Therefore,the proportion of233U loaded in the fertile blanket is set to match the power deposition.

Considering the temperature tolerance of the structural material,the temperature difference(temperature gradient)between the core and the fertile blanket should be as small as possible. Therefore, the inlet temperature of the fertile blanket is set to 923 K, which is consistent with the average core inlet temperature, by setting the inlet boundary condition to 923 K. In addition, the volume flow rate is adjusted so that the average outlet temperature of the fertile blanket is consistent with the average core outlet temperature (1023 K).

The velocity vector profile of the fertile blanket for the three cases is shown in Fig. 15. The inlet jet for case 1(Fig. 15a) is larger than those for case 2 (Fig. 15b) and case 3 (Fig. 15c). The main reason that the inlet velocity decreases as the number of external circuits increases is that when there are more external circuits, a lower inlet velocity for each pipe can maintain the same volume flow rate of the fertile blanket. A strong inlet jet causes a large reverse flow phenomenon in the entire fertile blanket, as illustrated in Fig. 15a. A small vortex appears in Fig. 15b,c,but no reverse flow occurs near the arc wall owing to the relatively low inlet flow velocity.

Figure 16 shows the temperature distribution of the fertile blanket for the three cases. The temperature distribution of the fertile blanket becomes more uniform as the number of external circuits increases. To maintain a constant outlet temperature of approximately 1023 K,the total volume flow rates of the fertile blanket for the three cases are almost identical to each other owing to their similar geometries. As illustrated above, the required inlet flow velocity for case 1 is highest among the three cases to match the equivalent volume flow rate.As a result,the inlet jet in case 1 causes a highly nonuniform temperature distribution, and a relatively high temperature appears at the bottom arc, which is higher than the outlet temperature of the fertile blanket. Furthermore, a local high temperature,which is much higher than the outlet temperature, appears in the top cavity of the fertile blanket.By contrast,because no reverse flow occurs near the arc wall, the temperature distribution of the top cavity for cases 2 and 3 is relatively flat, and no hot spot forms.

Table 3 shows the maximum temperature,average outlet temperature, and required total volume flow rate of the fertile blanket. Case 3 has the lowest maximum temperature and the most uniform temperature distribution.Therefore,it has the best TH performance among the three cases. However, the calculated maximum temperature and temperature distribution for case 2 are also acceptable;this result indicates that eight cooling circuits with a total volume flow rate of approximately 0.0555 m3s-1can also produce an acceptable temperature distribution. Considering the engineering requirements, the geometry model for case 2 is recommended.

4 Conclusion

In the core of the MSFR benchmark model,a large flow recirculation zone appears near the fertile blanket, and the fuel salt is almost stagnant at the core center,which is close to the axial reflectors. The recirculation in this region produces a local hot spot on the wall (near the fertile blanket),and thus results in unnecessary thermal stress and material problems [12]. A preliminary conjugate heat transfer analysis of the core and fertile blanket indicates the feasibility of using a separate model for the MSFR TH analysis and demonstrates that the nuclear heat generated in the fertile blanket cannot be removed by heat convection due to the fuel salt flow. Therefore, an external cooling system for the fertile blanket has to be considered to remove the heat generated by fission and neutron capture by thorium. In this work, a neutronics and TH coupling scheme between the DAG-OpenMC and OpenFOAM codes is adopted to further optimize the core design and conduct a preliminary TH analysis of the fertile blanket.The main conclusions are as follows:

(1) By optimizing the core geometry, the large recirculation region near the fertile blanket in the reference MSFR model is eliminated, and a relatively flat velocity distribution is obtained for all three studied geometries. Moreover, the stagnation zones disappear in the velocity field of Geometry III. The hot spot that appears in the upper core cannot be effectively eliminated by curving only the upper wall of the core. However, the peak temperature for Geometry III is much lower than that for the other two geometries. Furthermore, a more uniform distribution and a relatively small radial gradient of the fuel temperature are obtained for Geometry III.In terms of the thermal hydraulics, Geometry III,which has curved upper and lower core walls,has the best performance.

(2) For the fertile blanket, the strong inlet jet in case 1 produces a significant reverse flow in the entire fertile blanket, which results in a highly nonuniform temperature distribution and a relatively high temperature at the bottom arc. In addition, a local high temperature appears in the top cavity of the fertile blanket. For cases 2 and 3, no reverse flow occurs near the arc wall,and no hot spot forms owing to the relatively small inlet jet and relatively flat temperature distribution in the top cavity. Case 3 (with 16 circuits) has the most uniform temperature distribution and the lowest maximum temperature;the calculated result for case 2 (with 8 circuits) is also acceptable. Considering the engineering requirements, the geometry model of case 2 is recommended.

Table 3 Temperature results and required flow rate of fertile blanket