Conceptual Design of a 10MW Multipurpose Research Reactor Using VVR-KN Fuel

-e paper presents a conceptual design of a 10MWmultipurpose 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. -e core consists of 36 fuel assemblies: 27 standard fuel assemblies and 9 control fuel assemblies with the uranium density of 2.8 gU/cm and the U enrichment of 19.75 wt.%. -e 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. -e highest power rate and the highest discharged burnup of fuel assembly are 393.49 kW and 56.74% loss of U, respectively. -ermal hydraulics analysis has also been conducted using the PLTEMP4.2 code for evaluating the safety parameters at a steady state of the hottest channel. -e maximum temperatures of coolant and fuel cladding are 66.0°C and 83.0°C, respectively. -is value is lower than the design limit of 98°C for cladding temperature. -ermal fluxes at the vertical irradiation channels and the horizontal beam ports have been evaluated. -e maximum thermal fluxes of 2.5×10 and 8.9×10 n·cm·s are found at the neutron trap and the beryllium reflector, respectively.


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]. e 500 kW Dalat Nuclear Research Reactor (DNRR) (located in Dalat, Vietnam) is currently a unique research reactor in the country. e DNRR was reconstructed and upgraded from a 250 kW TRIGA Mark II reactor in the early 1980s. e core is loaded with Russian VVR-M2 fuel. e 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. e DNRR has been operational safely for about 35 years, and it is planned to continue the reactor operation for about a decade. en, 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 × 10 14 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 platetype material test reactor (MTR) fuel and heavy water as a reflector. e 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/cm 3 and burnable poison integrated fuel elements was mostly used in the designs of high-performance MPRRs. e KJRR reactor in Korea is designed using the MTR fuel with a high density of 8.0 gU/cm 3 and beryllium and graphite reflectors [8]. e plate-type fuel with the uranium density of 3 gU/cm 3 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 × 10 14 n·cm −2 ·s −1 [6]. e OPAL reactor achieves the maximum thermal flux of 3 × 10 14 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 × 10 14 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. e VVR-KN fuel assembly contains low enriched uranium fuel with the density of 2.8 gU/cm 3 . e 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.
us, 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 × 10 14 and 0.5-1.0 × 10 14 n·cm −2 ·s −1 , respectively. ermal hydraulics analysis has also been performed using the PLTEMP4.2 code to evaluate the safety parameters at a steady state of the MPRR core. e 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.
e 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 safetyrelated parameters in the hottest channel of the MPRR core. Finally, concluding remarks are derived in Section 6. 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. e 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. e fuel meat is ceramic-metal composition of UO 2 -Al with the uranium density of 2.8 gU/cm 3 . e 235 U enrichment is 19.75 wt.%. e total 235 U 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. e active fuel length is 600 mm with a standard deviation of 2 mm. e total length of neutron absorption (B 4 C or stainless steel material) in control rods (CRs) is 680 mm. e 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. e outer dimensions of the beryllium rod are the same as those of the fuel assembly [17,18].

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]. e reactor type is open swimming pool. e core is cooled and moderated by light water and reflected by beryllium. e 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 crosssectional view of the MPRR core. e 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. e peripheral ring of the core is loaded with 22 beryllium rods and 2 vertical channels. e 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. e core is controlled by nine control rods: two safety rods (SR1 and SR2), six shim rods (ShR1-ShR6), and an   automatic regulating rod (AR). e 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. e insertion of the shim rods is adjusted during the core operation to maintain criticality. e 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. e safety and shim rods use B 4 C as a neutron absorbing material with the density of 1.69 g/cm 3 . e AR rod is made of stainless steel with the density of 7.8 g/cm 3 for a smaller reactivity worth.

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]. e 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]. e 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 10 6 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 Cold neutron source Figure 2: Horizontal cross-sectional view of the MPRR core. HBT: horizontal beam tubes; SiD-1 to SiD4: neutron transmutation doping of silicon; S1-S9: small irradiation channels; B1-B2: big irradiation channels. 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. e 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. e Russian built-in correlation in the PLTEMP4.2 code for the VVR-KN fuel was used for thermal hydraulics analysis of the hottest FA. e heat removal of the reactor is carried out by forced convection with the direction from the core top to the core bottom.

Burnup Reactivity.
e 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. e 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 k eff curve of the MPRR core as a function of effective fullpower days (EFPDs). e operation time of the core with full power is 86 days. e 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. e xenon poisoning effect at the BOC is evaluated at about 2970 pcm. e excess reactivity at the EOC after a cooling time of about 7 days, when the xenon poisoning effect is ended, is about 3879 pcm. Figure 6 displays the radial relative power distribution and the radial power peaking factor (PPF) in the core at the beginning of cycle (BOC). e 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 (F xy ) is 1.415, which appears at the FA-1 assembly located at cell 8-4. is 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. e 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. e highest power rate of 106.21 kW and the average power density per fuel volume of 1148.46 W/cm 3 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/cm 3 . is value corresponds to a power peaking factor of 2.2. Figure 8 depicts the radial distribution of fuel burnup in the core at the end of cycle (EOC). e 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 235 U, respectively. e highest burnup of the FAs is 56.74% loss of 235 U, which appears at location 5-4. In the burnup calculations, the fuel assembly is axially divided into five nodes with equal volumes. e highest burnup level in the axial nodes is about 66% loss of 235 U. It can be seen from Figure 8 that the FA-2 assemblies have higher burnup levels of 49-56% loss of 235 U.

Burnup Distribution.
is is because the FA-2 assemblies were used together with the CRs during the reactor lifetime. e FAs with the burnup values less than 25% loss of 235 U 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.

Control Rod Reactivity Worths.
e 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. e 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. e 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. e lowest control rod worth of 1156 pcm is found at ShR1, while the highest value of 2507 pcm is found at ShR3. e 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 B 4 C. e shutdown margin is about −2330 pcm corresponding to the k eff value of 0.97722.

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. e 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. en, the FTC at the BOC is obtained as −2.467 pcm/K. e 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     Science and Technology of Nuclear Installations 7 found that the CTC values are about −11.594 and −12.327 pcm/K in the two ranges of coolant temperatures, respectively. e void reactivity coefficient was also evaluated when the coolant in the active core was voided from 0 to 5%. e void reactivity coefficient at the BOC is obtained as −259.9 pcm per one percent of coolant voided. e results show that the reactivity coefficients of the MPRR core are all negative. ese values satisfy the design targets and the safety requirements recommended in IAEA guidelines [24,25]. e maximum thermal flux at the neutron cold source is 1.57 × 10 13 n·cm −2 ·s −1 , while these values at the SiD positions are from 7.16 × 10 12 to 7.60 × 10 13 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. e values of thermal fluxes at the horizontal beam tubes close to the core are in the range of 1.15-7.91 × 10 13 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]. e results imply that the neutron fluxes at the vertical irradiation channels and the horizontal neutron beam tubes are sufficient for the application purposes.

Thermal Hydraulics Analysis
ermal 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. e hottest fuel assembly with the highest power rate of 393.49 kW is located at cell 8-4 as shown in Figure 3. e 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. e 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. e temperatures of fuel, cladding surface, and coolant and the ONBR and DNBR parameters at a steady state have been evaluated.     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. e values are lower than the safety limits of the design. e minimum values of the ONBR and the DNBR are 1.80 and 7.02, respectively, which satisfy the safety limits. e 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].

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. e core is loaded with 36 LEU fuel assemblies: 27 standard fuel assemblies and 9 control fuel assemblies. e cycle length is 86 effective full-power days with the excess reactivity of 9600 and 1039 pcm at the BOC and the EOC, respectively. e highest discharged burnup of the fuel assembly is 56.74% loss of 235 U. e highest power rate per assembly is 393.49 kW. e maximum power density per fuel volume is about 1442.59 W/cm 3 , 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. e reactivity worths of the nine control rods ensure the safe operation of the core and provide sufficient shutdown margin. e results show that the neutronics performance of the newly designed MPRR core meets the design targets. ermal hydraulics analysis at the steady state has been performed for the hottest fuel assembly using the PLTEMP4.2 code. e 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. e minimum ONBR and DNBR are about 1.80 and 7.02, respectively. ese values satisfy the safety limits. For the purposes of application, the neutron fluxes at the vertical and horizontal irradiation channels have been analyzed. e maximum thermal neutron flux of 2.5 × 10 14 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 × 10 13 n·cm −2 ·s −1 .

Data Availability
Data will be made available upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this paper.