Research Article  Open Access
NhiDien Nguyen, KienCuong Nguyen, TonNghiem Huynh, DoanHaiDang Vo, HoaiNam Tran, "Conceptual Design of a 10 MW Multipurpose Research Reactor Using VVRKN Fuel", Science and Technology of Nuclear Installations, vol. 2020, Article ID 7972827, 11 pages, 2020. https://doi.org/10.1155/2020/7972827
Conceptual Design of a 10 MW Multipurpose Research Reactor Using VVRKN Fuel
Abstract
The paper presents a conceptual design of a 10 MW multipurpose nuclear research reactor (MPRR) loaded with the lowenriched uranium (LEU) VVRKN fuel type. Neutronics and burnup calculations have been performed using the REBUSMCNP6 linkage system code and the ENDF/BVII.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/cm^{3} and the ^{235}U enrichment of 19.75 wt.%. The cycle length of the core is 86 effective fullpower 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 ^{235}U, 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 × 10^{14} and 8.9 ×10^{13} 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 VVRM2 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 × 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. The IRT1 reactor in Libya is the only one using tubetype 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 highperformance MPRRs. The 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]. The platetype 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 RA10 reactor (30 MW) in Argentina [11], the HFIR reactor (85 MW) in USA, the FRMII reactor (20 MW) in Germany [12], the HANARO reactor (30 MW) in Korea [13, 14], JRR3M (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 FRMII 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]. The 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 FRMII 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 tubetype VVRKN fuel assemblies has been conducted. Neutronics and burnup calculations have been performed using the REBUSMCNP6 linkage system code. The VVRKN fuel assembly contains low enriched uranium fuel with the density of 2.8 gU/cm^{3}. 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 × 10^{14} and 0.5–1.0 × 10^{14} 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 VVRKN fuel assembly and the design targets of the MPRR core. Section 3 describes the calculation methods using the REBUSMCNP6 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 safetyrelated parameters in the hottest channel of the MPRR core. Finally, concluding remarks are derived in Section 6.
2. VVRKN Fuel and MPRR Core Design
2.1. VVRKN Fuel
Figure 1 displays the crosssectional views of two lowenriched uranium (LEU) VVRKN fuel assemblies (FAs) denoted as FA1 and FA2, respectively. FA1 is a standard fuel assembly consisting of eight concentric tubular fuel elements (FEs). Seven outer FEs of the FA1 assembly have hexagonal shapes, while the inner FE of the FA1 assembly is cylindrical. A cylindrical structural tube is located in the center of the FA1 assembly. FA2 is a control fuel assembly, which consists of five concentric tubular FEs similar to the outermost FEs of the FA1 assembly. A cylindrical guide tube for control rod (CR) insertion is located at the center of the FA2 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 ceramicmetal composition of UO_{2}Al with the uranium density of 2.8 gU/cm^{3}. The ^{235}U enrichment is 19.75 wt.%. The total ^{235}U amounts in the FA1 and FA2 assemblies are 248.2 g and 197.6 g, respectively. Cladding and structural materials are made of aluminumalloy SAV1. 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 (B_{4}C 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].
(a)
(b)
2.2. MPRR Core Design
In the present work, the WWRK 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 VVRKN FAs [17, 18]. Figure 2 shows the horizontal crosssectional 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 FA1 type and 9 assemblies of FA2 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 (ShR1ShR6), 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 B_{4}C as a neutron absorbing material with the density of 1.69 g/cm^{3}. The AR rod is made of stainless steel with the density of 7.8 g/cm^{3} 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/BVII.0 nuclear data library [19, 20]. The MCNP code is a generalpurpose Monte Carlo radiationparticle 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 MCNP6REBUS linkage system code [21, 22]. The REBUSPC 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 steadystate flow and temperature solution for a fuel assembly or a reactor core in the subcooled boiling regime [23]. A number of thermalhydraulic 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 FA1 and FA2 assemblies as cylindrical tubes by keeping the equivalent areas of the FEs. The Russian builtin correlation in the PLTEMP4.2 code for the VVRKN 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 FA1 assemblies and 9 FA2 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 k_{eff} curve of the MPRR core as a function of effective fullpower 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 (F_{xy}) is 1.415, which appears at the FA1 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, FE1 is the outermost FE, and FE8 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/cm^{3} are found at the outermost fuel element FE1 as shown in Table 2. Figure 7 displays the axial power profiles in the FEs of the hottest assembly, that is, the FA1 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}. 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 ^{235}U, respectively. The 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. The highest burnup level in the axial nodes is about 66% loss of ^{235}U. It can be seen from Figure 8 that the FA2 assemblies have higher burnup levels of 49–56% loss of ^{235}U. This is because the FA2 assemblies were used together with the CRs during the reactor lifetime. The 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.
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 B_{4}C. The shutdown margin is about −2330 pcm corresponding to the k_{eff} 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), epithermal 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 × 10^{14} and 2.54 × 10^{14} 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 × 10^{13} n·cm^{−2}·s^{−1}. The 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. The 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]. The results imply that the neutron fluxes at the vertical irradiation channels and the horizontal neutron beam tubes are sufficient for the application purposes.
 
^{a}Read as 2.02 × 10^{14}. 
 
^{a}Read as 7.91 × 10^{13}. 
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 builtin correlation for the VVRKN 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 VVRKN fuel has been conducted using the REBUSMCNP6 linkage code and the ENDF/BVII.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 fullpower 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 ^{235}U. The highest power rate per assembly is 393.49 kW. The 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. 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 × 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
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.CN50/15.
References
 IAEA (International Atomic Energy Agency), “The applications of research reactors,” Tech. Rep., IAEA, Vienna, Austria, 2001, IAEATECDOC1234. View at: Google Scholar
 IAEA (International Atomic Energy Agency), “Applications of research reactors,” Tech. Rep., IAEA, Vienna, Austria, 2014, IAEA Nuclear Energy Series No. NPT5.3. View at: Google Scholar
 G. Phan, H. N. Tran, K. C. Nguyen et al., “Comparative analysis of the Dalat nuclear research reactor with HEU fuel using SRAC and MCNP5,” Science and Technology of Nuclear Installations, vol. 2017, Article ID 2615409, 10 pages, 2017. View at: Publisher Site  Google Scholar
 Q. B. Do, G. T. T. Phan, K. C. Nguyen, Q. H. Ngo, and H. N. Tran, “Criticality and control rod worth analysis of the DNRR research reactor using the SRAC and MCNP5 codes,” Nuclear Engineering and Design, vol. 343, pp. 197–209, 2019. View at: Publisher Site  Google Scholar
 V. K. Raina, K. Sasidharan, S. Sengupta, and T. Singh, “Multipurpose research reactor,” Nuclear Engineering and Design, vol. 236, pp. 770–783, 2006. View at: Publisher Site  Google Scholar
 F. E. Teruel and Rizwanuddin, “An innovative research reactor design,” Nuclear Engineering and Design, vol. 239, pp. 395–407, 2009. View at: Publisher Site  Google Scholar
 C. G. Seo and N. Z. Cho, “A core design concept for multipurpose research reactors,” Nuclear Engineering and Design, vol. 252, pp. 34–41, 2012. View at: Publisher Site  Google Scholar
 C. G. Seo, H. T. Chae, B. C. Lee, B. J. Jun, and I. C. Lim, “Conceptual nuclear design of the kijang research reactor,” in Proceedings of the Transactions of the RRFM2013 European Research Reactor Conference, pp. 491–496, Saint Petersburg, Russia, April 2013. View at: Google Scholar
 S. Kim, “The OPAL (open pool Australian lightwater) reactor in Australia,” Nuclear Engineering and Technology, vol. 38, no. 5, pp. 443–448, 2006. View at: Google Scholar
 E. Villarino and A. Doval, “INVAP’s research reactor designs,” Science and Technology of Nuclear Installations, vol. 2011, Article ID 490391, 6 pages, 2011. View at: Publisher Site  Google Scholar
 P. Cantero, P. Ramirez, F. Brollo et al., “Key technical characteristics related to the design of the RA10 Multipurpose Research Reactor,” in Proceedings of European Research Reactor Conference (RRFM/IGORR2016), pp. 683–693, Berlin, Germany, March 2016. View at: Google Scholar
 X. Lin, R. Henkelmann, A. Turler, H. Gerstenberg, and F. D. Corte, “Neutron flux parameters at irradiation positions in the new research reactor FRMII,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 564, no. 2, pp. 641–644, 2006. View at: Publisher Site  Google Scholar
 H. Kim, H. R. Kim, K. H. Lee, and J. B. Lee, “Design characteristics and startup tests of HANARO,” Journal of Nuclear Science and Technology, vol. 33, no. 7, pp. 527–538, 1996. View at: Publisher Site  Google Scholar
 K. N. Choo, M. S. Cho, S. W. Yang, and S. J. Park, “Contribution of HANARO irradiation technologies to national nuclear R&D,” Nuclear Engineering and Technology, vol. 46, no. 4, pp. 501–512, 2014. View at: Publisher Site  Google Scholar
 D. Campolina, A. C. L. da Costa, E. P. Andrade, A. A. C. Santos, and V. Vasconcelos, “Neutronic analysis of the fuel loaded irradiation loop device of the RMB Multipurpose Brazilian Reactor,” Progress in Nuclear Energy, vol. 104, pp. 109–116, 2018. View at: Publisher Site  Google Scholar
 M. K. Jaradat, V. Radulovi, C. J. Park, L. Snoj, and S. M. Alkhafaji, “Verification of MCNP6 model of the Jordan Research and Training Reactor (JRTR) for calculations of neutronic parameters,” Annals of Nuclear Energy, vol. 96, pp. 96–103, 2016. View at: Publisher Site  Google Scholar
 N. A. Hanan and P. L. Garner, “Neutronic, steady state and transient analyses for the kazakhstan WWRK reactor with LEU fuel: ANL independent verification results,” Tech. Rep., Argonne National Laboratory, Lemont, IL, USA, 2015, ANL/RTR/TM15/7, 120766. View at: Publisher Site  Google Scholar
 I. T. Tretiyakov, S. A. Sokolov, V. I. Trushkin et al., “Project development for promising pooltype Research Reactors,” in Proceedings of the Transactions of the RRFM2013 European Research Reactor Conference, pp. 505–510, Saint Petersburg, Russia, April 2013. View at: Google Scholar
 ORNL, MCNP6.1/MCNP5/MCNPX, Monte Carlo NParticle Transport Code System Including MCNP6.1, MCNP51.60, MCNPX2.7.0 and Data Libraries, Contributed by Los Alamos National Laboratory, RSICC Computer Code Collection, Oak Ridge National Laboratory, Oak Ridge, TN, USA, 2013.
 M. B. Chadwick, P. Oblozinsky, M. Herman et al., “ENDF/BVII.0: next generation evaluated nuclear data library for nuclear science and technology,” Nuclear Data Sheets, vol. 107, pp. 2931–3060, 2006. View at: Publisher Site  Google Scholar
 N. A. Hanan, A. P. Olson, R. B. Pond, and J. E. Matos, “A Monte Carlo burnup code linking MCNP and REBUS,” in Proceedings of the 20th International Meeting on Reduced Enrichment for Research and Test Reactors, RERTR 1998, Sao Paulo, Brazil, October 1998. View at: Google Scholar
 J. G. Stevens, “The REBUSMCNP linkage,” Tech. Rep., Argonne National Laboratory, Lemont, IL, USA, 2008, ANL/RERTR/TM0804. View at: Google Scholar
 A. P. Olson and M. Kalimullah, “A users guide to the PLTEMP/ANL code,” Tech. Rep., Argonne National Laboratory, Lemont, IL, USA, 2015, ANL/RERTR/TM1122 Version 4.2. View at: Google Scholar
 IAEA (International Atomic Energy Agency), “Code on the safety of nuclear research reactors: design,” Tech. Rep., IAEA, Vienna, Austria, 1992, Safety Series No. 35S1. View at: Google Scholar
 IAEA (International Atomic Energy Agency), “Safety of Research Reactors. IAEA Safety Standards,” Tech. Rep., IAEA, Vienna, Austria, 2016, Specific Safety Requirements No. SSR3. View at: Google Scholar
Copyright
Copyright © 2020 NhiDien Nguyen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.