Abstract

The paper presents a conceptual design of a 10 MW multipurpose nuclear research reactor (MPRR) loaded with the low-enriched uranium (LEU) VVR-KN fuel type. Neutronics and burnup calculations have been performed using the REBUS-MCNP6 linkage system code and the ENDF/B-VII.0 data library. The core consists of 36 fuel assemblies: 27 standard fuel assemblies and 9 control fuel assemblies with the uranium density of 2.8 gU/cm3 and the 235U enrichment of 19.75 wt.%. The cycle length of the core is 86 effective full-power days with the excess reactivity of 9600 and 1039 pcm at the beginning of cycle and the end of cycle, respectively. The highest power rate and the highest discharged burnup of fuel assembly are 393.49 kW and 56.74% loss of 235U, respectively. Thermal hydraulics analysis has also been conducted using the PLTEMP4.2 code for evaluating the safety parameters at a steady state of the hottest channel. The maximum temperatures of coolant and fuel cladding are 66.0°C and 83.0°C, respectively. This value is lower than the design limit of 98°C for cladding temperature. Thermal fluxes at the vertical irradiation channels and the horizontal beam ports have been evaluated. The maximum thermal fluxes of 2.5 × 1014 and 8.9 ×1013 n·cm−2·s−1 are found at the neutron trap and the beryllium reflector, respectively.

1. Introduction

Research reactor is an important tool for supporting research and development in a wide range of scientific aspects such as neutron scattering, neutron activation analysis, material testing, medicine application, biological science, and education and training [1, 2]. The 500 kW Dalat Nuclear Research Reactor (DNRR) (located in Dalat, Vietnam) is currently a unique research reactor in the country. The DNRR was reconstructed and upgraded from a 250 kW TRIGA Mark II reactor in the early 1980s. The core is loaded with Russian VVR-M2 fuel. The operation of the DNRR is mainly for radioisotope production (RI), neutron activation analysis (NAA), basic and applied researches, and training [3, 4]. However, due to the limitation of power and experimental facilities, the DNRR cannot meet the increasing demands of radioisotope production for medical and industrial applications, facilities for nuclear physics experiments, and support to nuclear power development in the future. The DNRR has been operational safely for about 35 years, and it is planned to continue the reactor operation for about a decade. Then, it is highly desirable to replace the DNRR by a new research reactor with a higher power output. For that purpose, the country is planning to build a new multipurpose research reactor (MPRR) with the power of 10 MW. In addition to traditional applications, such as RI, NAA, and neutron imaging (NI), other advanced researches and applications, such as neutron transmutation doping of silicon (NTD), material science research and characterization by neutron scattering, and production of special isotopes using high neutron fluxes and so forth are being considered with the new MPRR. For most of the applications, high thermal neutron fluxes and high irradiation volumes are desirable, for example, thermal flux higher than 2.0 × 1014 n·cm−2·s−1, though there exist a number of constraints and difficulties to achieve this goal [5, 6].

A number of MPRRs with the power output greater than 10 MW currently in operation, under construction, and planned were reviewed [6, 7]. Most of the MPRRs use plate-type material test reactor (MTR) fuel and heavy water as a reflector. The IRT-1 reactor in Libya is the only one using tube-type fuel. In the last 30 years, the MTR fuel with the density of 4.8 gU/cm3 and burnable poison integrated fuel elements was mostly used in the designs of high-performance MPRRs. The KJRR reactor in Korea is designed using the MTR fuel with a high density of 8.0 gU/cm3 and beryllium and graphite reflectors [8]. The plate-type fuel with the uranium density of 3 gU/cm3 was also used in the MPRR design with the power of 20 MW for achieving high neutron fluxes at irradiation channels [5]. Several designs of MPRRs include the OPAL reactor (20 MW) in Australia [9, 10], the CARR reactor (60 MW) in China, the RA-10 reactor (30 MW) in Argentina [11], the HFIR reactor (85 MW) in USA, the FRM-II reactor (20 MW) in Germany [12], the HANARO reactor (30 MW) in Korea [13, 14], JRR-3M (20 MW) in Japan, the RMB reactor (30 MW) in Brazil [15], and the JRTR reactor (5 MW) in Jordan [16]. A comparison among several MPRRs shows that the FRM-II reactor achieves the highest thermal flux per unit power with a maximum thermal neutron flux of 8 × 1014 n·cm−2·s−1 [6]. The OPAL reactor achieves the maximum thermal flux of 3 × 1014 n·cm−2·s−1, allowing various applications of cold, thermal, and hot neutron beams, radioisotopes production, and silicon doping. Based on the capabilities of the OPAL reactor and the high thermal flux of the FRM-II reactor, an innovative design of a compact MPRR core with a power level of 10 MW and maximum thermal flux of 5 × 1014 n·cm−2·.s−1 was investigated [6].

In the present work, a conceptual design of 10 MW MPRR core loaded with the Russian tube-type VVR-KN fuel assemblies has been conducted. Neutronics and burnup calculations have been performed using the REBUS-MCNP6 linkage system code. The VVR-KN fuel assembly contains low enriched uranium fuel with the density of 2.8 gU/cm3. The active core is surrounded by a beryllium reflector to attain high thermal neutron fluxes and high performance of irradiation channels. In the design process, the numbers of fuel assemblies and control rods and their arrangement in the core are determined for obtaining enough excess reactivity during the core operation for compensating a number of irradiation applications, while ensuring safety requirements. Thus, neutronic characteristics of the MPRR core including burnup performance, power distribution and power peaking factor, control rod worths, reactivity coefficients, sizes and types of irradiation holes for a variety of applications, neutron fluxes at irradiation channels, and so forth have been investigated. One of the design objectives is to attain high thermal neutron fluxes in the core and irradiation channels with maximum thermal neutron fluxes at the core center and at the beryllium reflector of about 2.0–3.0 × 1014 and 0.5–1.0 × 1014 n·cm−2·s−1, respectively. Thermal hydraulics analysis has also been performed using the PLTEMP4.2 code to evaluate the safety parameters at a steady state of the MPRR core.

The paper is organized as follows. Section 2 presents a brief description of the VVR-KN fuel assembly and the design targets of the MPRR core. Section 3 describes the calculation methods using the REBUS-MCNP6 linkage system code and the PLTEMP4.2 code. The neutronics performance such as burnup reactivity, power distribution, burnup distribution, kinetic parameters, control rod worths, reactivity coefficients, and thermal neutron flux levels of the newly designed MPRR core is presented in Section 4. Section 5 describes the results of thermal hydraulic analysis for obtaining fuel and cladding temperatures and other safety-related parameters in the hottest channel of the MPRR core. Finally, concluding remarks are derived in Section 6.

2. VVR-KN Fuel and MPRR Core Design

2.1. VVR-KN Fuel

Figure 1 displays the cross-sectional views of two low-enriched uranium (LEU) VVR-KN fuel assemblies (FAs) denoted as FA-1 and FA-2, respectively. FA-1 is a standard fuel assembly consisting of eight concentric tubular fuel elements (FEs). Seven outer FEs of the FA-1 assembly have hexagonal shapes, while the inner FE of the FA-1 assembly is cylindrical. A cylindrical structural tube is located in the center of the FA-1 assembly. FA-2 is a control fuel assembly, which consists of five concentric tubular FEs similar to the outermost FEs of the FA-1 assembly. A cylindrical guide tube for control rod (CR) insertion is located at the center of the FA-2 assembly. The FE has the thickness of 1.6 mm, including the fuel meat with the thickness of 0.7 mm and two cladding sides with the thickness of 0.45 mm. The fuel meat is ceramic-metal composition of UO2-Al with the uranium density of 2.8 gU/cm3. The 235U enrichment is 19.75 wt.%. The total 235U amounts in the FA-1 and FA-2 assemblies are 248.2 g and 197.6 g, respectively. Cladding and structural materials are made of aluminum-alloy SAV-1. Ribs of 1.5 mm height on the outer surfaces of the FEs are for maintaining the water gap of 2 mm between the FEs and providing additional stiffening. The active fuel length is 600 mm with a standard deviation of 2 mm. The total length of neutron absorption (B4C or stainless steel material) in control rods (CRs) is 680 mm. The follower of the CR is an aluminum rod with the radius of 12.5 mm and the length of 700 mm. Beryllium rod is solid beryllium with the pitch of 66.3 mm and the height of 600 mm. The outer dimensions of the beryllium rod are the same as those of the fuel assembly [17, 18].

2.2. MPRR Core Design

In the present work, the WWR-K reactor with the power of 6 MW in Kazakhstan is taken as a reference for designing a new MPRR with the power of 10 MW [17]. The reactor type is open swimming pool. The core is cooled and moderated by light water and reflected by beryllium. The reactor core has the height of 180 cm and the radius of 113.5 cm. Table 1 presents the main design parameters and the design targets of the MPRR core using the VVR-KN FAs [17, 18]. Figure 2 shows the horizontal cross-sectional view of the MPRR core. The active core consists of 61 hexagonal cells for containing fuel assemblies, neutron trap, irradiation channels, and beryllium rods. Figure 3 shows the core map with the coordinates of hexagonal cells. In this design, the final core consists of 36 fuel assemblies: 27 assemblies of FA-1 type and 9 assemblies of FA-2 type. A neutron trap is located at the core center, and two vertical irradiation channels are located at the peripheral ring. The peripheral ring of the core is loaded with 22 beryllium rods and 2 vertical channels. The beryllium rods provide an additional reflector to the core. It is also possible to create additional irradiation channels by replacing the beryllium rods with vertical channels.

The core is controlled by nine control rods: two safety rods (SR1 and SR2), six shim rods (ShR1-ShR6), and an automatic regulating rod (AR). The safety rods are normally withdrawn from the core during the reactor operation and are inserted into the core only in the case of reactor shutdown. The insertion of the shim rods is adjusted during the core operation to maintain criticality. The AR rod is partially inserted into the core during the reactor operation. In the case of emergency shutdown, the control rods fall into the core by the gravity. The safety and shim rods use B4C as a neutron absorbing material with the density of 1.69 g/cm3. The AR rod is made of stainless steel with the density of 7.8 g/cm3 for a smaller reactivity worth.

3. Calculation Methods

Neutronic calculations for designing the 10 MW MPRR core have been conducted using the MCNP6 and the ENDF/B-VII.0 nuclear data library [19, 20]. The MCNP code is a general-purpose Monte Carlo radiation-particle transport code developed at Los Alamos National Laboratory [19]. MCNP6 has an ability to describe precisely the complex geometry of the research reactor. Burnup calculations have been performed using the MCNP6-REBUS linkage system code [21, 22]. The REBUS-PC code was extended to interface with the MCNP code to perform burnup calculations using neutron fluxes and cross sections from the MCNP code [22]. To determine a final core configuration, burnup calculations have been performed for six operating cycles with different core configurations and power levels. In the MCNP6 calculations, the number of neutron history per cycle was chosen as 106 with 10 skip cycles and 200 calculated cycles to obtain the statistic error of the effective multiplication factor of about 0.00006 and the error of neutron fluxes within 0.1%.

PLTEMP4.2 is a thermal hydraulics code developed by Los Alamos National Laboratory for analyzing a steady-state flow and temperature solution for a fuel assembly or a reactor core in the subcooled boiling regime [23]. A number of thermal-hydraulic correlations are available in the code to evaluate related safety margins such as onset of nucleate boiling ratio (ONBR), departure from nucleate boiling ratio (DNBR), and onset of flow instability (FI). In the present work, the PLTEMP4.2 code has been used to perform thermal hydraulics analysis at a steady state of the newly designed MPRR core for obtaining the safety parameters such as fuel and cladding temperatures, ONBR, and DNBR. The PLTEMP4.2 code models hexagonal FEs of the FA-1 and FA-2 assemblies as cylindrical tubes by keeping the equivalent areas of the FEs. The Russian built-in correlation in the PLTEMP4.2 code for the VVR-KN fuel was used for thermal hydraulics analysis of the hottest FA. The heat removal of the reactor is carried out by forced convection with the direction from the core top to the core bottom.

4. Neutronics Performance of the MPRR Core

4.1. Burnup Reactivity

The first core configuration was established with 26 FAs (17 FA-1 assemblies and 9 FA-2 assemblies) and the power of 6 MW. Burnup calculations were then conducted to establish the next cycle with the increase of power level. After five cycles with the increase of power levels and the number of FAs loaded into the core, a final cycle was obtained. Figure 4 shows the final core configuration of the MPRR with the power of 10 MW. Detailed neutronic and burnup calculations have been conducted to determine the neutronic and burnup performance of the MPRR core, such as burnup reactivity, power distribution, burnup distribution, and reactivity coefficients. In the burnup calculations, the beryllium poisoning effect of the beryllium rods and the reflector was also taken into account. The positions of the CRs, especially the shim rods, were adjusted to ensure nearly critical conditions as in the operation conditions of the reactor core. Figure 5 shows the keff curve of the MPRR core as a function of effective full-power days (EFPDs). The operation time of the core with full power is 86 days. The excess reactivity at the BOC and the EOC of the core are 9600 and 1039 pcm, respectively. It is worth noticing that the excess reactivity at the BOC is calculated without the xenon poisoning effect, whereas this effect is taken into account in the value at the EOC. The xenon poisoning effect at the BOC is evaluated at about 2970 pcm. The excess reactivity at the EOC after a cooling time of about 7 days, when the xenon poisoning effect is ended, is about 3879 pcm.

4.2. Power Distribution

Figure 6 displays the radial relative power distribution and the radial power peaking factor (PPF) in the core at the beginning of cycle (BOC). The relative power distribution was obtained by normalizing the power distribution over the core volume. Since the total power is 10 MW and the core consists of 36 FAs, the average power rate per assembly is about 277.78 kW. One can see that the radial PPF (Fxy) is 1.415, which appears at the FA-1 assembly located at cell 8–4. This PPF corresponds to the power rate of 393.49 kW of the hottest FA. Since the power distribution of the hottest FA is used as an input for thermal hydraulics and safety analysis, detailed tally option in the MCNP6 has been applied to calculate the axial power profiles of the FEs in the hottest FA. The FEs are numbered from the outer to the inner fuel elements in the assembly, that is, FE-1 is the outermost FE, and FE-8 is the element at the center of the assembly. Table 2 presents the total power rates of the FEs in the hottest FA. One can see that the power rates and power densities increase from the central tube to the outer tube in the assembly. The highest power rate of 106.21 kW and the average power density per fuel volume of 1148.46 W/cm3 are found at the outermost fuel element FE-1 as shown in Table 2. Figure 7 displays the axial power profiles in the FEs of the hottest assembly, that is, the FA-1 assembly located at cell 8–4. It can be seen that the power peak appears at the middle core, and the maximum power density per fuel volume is 1442.59 W/cm3. This value corresponds to a power peaking factor of 2.2.

4.3. Burnup Distribution

Figure 8 depicts the radial distribution of fuel burnup in the core at the end of cycle (EOC). The fuel burnup levels can be categorized into three classes with the values in the ranges of 20–25%, 35–45%, and 45–55% loss of 235U, respectively. The highest burnup of the FAs is 56.74% loss of 235U, which appears at location 5–4. In the burnup calculations, the fuel assembly is axially divided into five nodes with equal volumes. The highest burnup level in the axial nodes is about 66% loss of 235U. It can be seen from Figure 8 that the FA-2 assemblies have higher burnup levels of 49–56% loss of 235U. This is because the FA-2 assemblies were used together with the CRs during the reactor lifetime. The FAs with the burnup values less than 25% loss of 235U are the fresh FAs loaded at the BOC of the core configuration. Since the locations of the CRs in the core can be changed from cycle to cycle upon the configuration design, the loading pattern for the next cycle can be determined with more flexible options.

4.4. Control Rod Reactivity Worths

The reactivity worths of the control rods and the shutdown margin have been evaluated and met requirements with condition when all safety CRs were withdrawn from the reactor core. The reactivity worths were calculated with the full insertion of the control rods. Meanwhile, the shutdown margin is defined when all ShR rods are fully inserted, SR1 is completely withdrawn, SR2 is fully inserted, and AR is at the middle of the core. Table 3 presents the reactivity worths of the nine control rods in the MPRR at the BOC. The reactivity worths of the SR1 and SR2 rods are 2051 pcm and 2102 pcm, respectively. When the two SR rods are simultaneously fully inserted into the core, the total reactivity worth is 4629 pcm. Meanwhile, the reactivity worths of the ShR rods are within the range of 1156–2507 pcm. One can see that the values of control rod worths are varied in a wide range depending on the positions of the ShR rods. The lowest control rod worth of 1156 pcm is found at ShR1, while the highest value of 2507 pcm is found at ShR3. The AR rod has the reactivity worth of 395 pcm, which is the smallest value compared to the other rods. It is because the AR rod is made of stainless steel, while other safety and shim rods are made of B4C. The shutdown margin is about −2330 pcm corresponding to the keff value of 0.97722.

4.5. Reactivity Coefficients

In order to evaluate the safety parameters of the new MPRR core, the feedback reactivity coefficients associated with the change of fuel and coolant temperatures and coolant void fraction at the BOC have been calculated. The temperature of the reflector (water and beryllium) was assumed to be equal to room temperature. To calculate the fuel temperature coefficient (FTC), the fuel temperature was assumed to vary in the range of 294–600 K. Then, the FTC at the BOC is obtained as −2.467 pcm/K. The coolant temperature coefficients (CTC) were evaluated in two ranges of coolant temperature of 294–350 K and 350–400 K, where the linear change of the CTC is valid. It is found that the CTC values are about −11.594 and −12.327 pcm/K in the two ranges of coolant temperatures, respectively. The void reactivity coefficient was also evaluated when the coolant in the active core was voided from 0 to 5%. The void reactivity coefficient at the BOC is obtained as −259.9 pcm per one percent of coolant voided. The results show that the reactivity coefficients of the MPRR core are all negative. These values satisfy the design targets and the safety requirements recommended in IAEA guidelines [24, 25].

4.6. Neutron Fluxes at Irradiation Channels and Beam Tubes

Neutron fluxes in thermal energy range (E < 0.625 eV), epi-thermal range (0.625 eV–0.821 MeV), and fast energy range (0.821–10.0 MeV) at the vertical irradiation channels and the horizontal beam tubes have been evaluated at the BOC. Table 4 shows the average and maximum values of neutron fluxes in three groups at the vertical irradiation channels in the core and in the reflector. The highest thermal flux is found at the neutron trap with the average and maximum values of 2.02 × 1014 and 2.54 × 1014 n·cm−2·s−1, respectively. Thus, the neutron trap can be exploited for reactor fuel and material testing. The maximum thermal flux at the beryllium reflector is about 8.9 × 1013 n·cm−2·s−1. The maximum thermal flux at the neutron cold source is 1.57 × 1013 n·cm−2·s−1, while these values at the SiD positions are from 7.16 × 1012 to 7.60 × 1013 n·cm−2·s−1. Table 5 shows the neutron fluxes at the six horizontal beam tubes at the positions close to the core and near the outer reflector. The values of thermal fluxes at the horizontal beam tubes close to the core are in the range of 1.15–7.91 × 1013 n·cm−2·s−1. It can be seen that the neutron fluxes at the irradiation channels of the newly designed MPRR meet the design targets as shown in Table 1 and the utilization requirements shown in [1, 2]. The results imply that the neutron fluxes at the vertical irradiation channels and the horizontal neutron beam tubes are sufficient for the application purposes.

5. Thermal Hydraulics Analysis

Thermal hydraulics analysis at the steady state was performed for the hottest FA using the PLTEMP4.2 code to determine the safety parameters of the newly designed MPRR core. The hottest fuel assembly with the highest power rate of 393.49 kW is located at cell 8–4 as shown in Figure 3. The power rates and the axial power profiles of the FEs of the hottest assembly are displayed in Table 2 and Figure 7, respectively. For ensuring the safety of cladding temperature, the maximum cladding surface temperature must be lower than 98°C, and the ONBR must be greater than 1.3 as given in Table 1. In the thermal hydraulics analysis using the PLTEMP4.2 code, the built-in correlation for the VVR-KN fuel was used. The conditions of the inlet coolant are specified as follows: (1) the inlet coolant temperature is 45°C, which is the maximum inlet temperature allowed during operation, (2) the inlet pressure is 0.135 MPa at the core top, and (3) the coolant flow rate is 5.58 kg/s. The temperatures of fuel, cladding surface, and coolant and the ONBR and DNBR parameters at a steady state have been evaluated.

Figures 9 and 10 display the axial profiles of the coolant and cladding temperatures in the hottest fuel assembly in the case without the application of a hotspot factor. One can see that the highest coolant temperature and the maximum cladding temperature are of 66°C and 83°C, respectively. The values are lower than the safety limits of the design. The minimum values of the ONBR and the DNBR are 1.80 and 7.02, respectively, which satisfy the safety limits. The results of the thermal hydraulics analysis at the steady state have confirmed that the new design of the MPRR core meets the safety criteria in the IAEA guidelines [24, 25].

6. Conclusions

Conceptual design of a 10 MW MPRR core using the VVR-KN fuel has been conducted using the REBUS-MCNP6 linkage code and the ENDF/B-VII.0 data library. The core is loaded with 36 LEU fuel assemblies: 27 standard fuel assemblies and 9 control fuel assemblies. The cycle length is 86 effective full-power days with the excess reactivity of 9600 and 1039 pcm at the BOC and the EOC, respectively. The highest discharged burnup of the fuel assembly is 56.74% loss of 235U. The highest power rate per assembly is 393.49 kW. The maximum power density per fuel volume is about 1442.59 W/cm3, which corresponds to a power peaking factor of 2.2. Feedback reactivity coefficients regarding the change of fuel and coolant temperatures and coolant void fraction are all negative. The reactivity worths of the nine control rods ensure the safe operation of the core and provide sufficient shutdown margin. The results show that the neutronics performance of the newly designed MPRR core meets the design targets. Thermal hydraulics analysis at the steady state has been performed for the hottest fuel assembly using the PLTEMP4.2 code. The maximum coolant and cladding temperatures at the hottest channel are 66.0°C and 83.0°C, respectively, without the application of a hotspot factor. The minimum ONBR and DNBR are about 1.80 and 7.02, respectively. These values satisfy the safety limits. For the purposes of application, the neutron fluxes at the vertical and horizontal irradiation channels have been analyzed. The maximum thermal neutron flux of 2.5 × 1014 n·cm−2·s−1 is obtained at the neutron trap located in the core center, whereas this value at the beryllium reflector is about 8.9 × 1013 n·cm−2·s−1.

Data Availability

Data will be made available upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors are grateful to the administrative staffs of VINATOM and DNRI for their kind support. The staffs at Reactor Physics and Engineering Department of DNRI are acknowledged for their valuable technical discussions. Collaborations with the staffs of Argonne National Lab. and Institute of Nuclear Physics of Kazakhstan are greatly appreciated. This research was supported by Ministry of Science and Technology of Vietnam under Grant no. DTDL.CN-50/15.