Research Article  Open Access
KienCuong Nguyen, VinhVinh Le, TonNghiem Huynh, BaVien Luong, NhiDien Nguyen, "SteadyState ThermalHydraulic Analysis of the LEUFueled Dalat Nuclear Research Reactor", Science and Technology of Nuclear Installations, vol. 2021, Article ID 6673162, 10 pages, 2021. https://doi.org/10.1155/2021/6673162
SteadyState ThermalHydraulic Analysis of the LEUFueled Dalat Nuclear Research Reactor
Abstract
This paper presents results of steadystate thermalhydraulic analysis for the designed working core of the Dalat Nuclear Research Reactor (DNRR) using the PLTEMP/ANL code. The core was designed to be loaded with 92 lowenriched uranium (LEU) VVRM2 fuel bundles (FBs) and 12 beryllium rods surrounding a neutron trap at the core center, for replacement of the previous core with 104 highenriched uranium (HEU) VVRM2 FBs. Before using this code for thermohydraulic analysis of the designed LEU working core, it was validated by comparing calculation results with experimental data collected from the HEU working core of the DNRR. The discrepancy between calculated results and measured data was at the maximum about 0.8°C and 1.5°C of fuel cladding and outlet coolant temperatures, respectively. In the design calculation, thermohydraulic safety was confirmed through evaluation of the fuel cladding and coolant temperatures, as well as of other safety parameters such as Departure from Nucleate Boiling Ratio (DNBR) and Onset of Nucleate Boiling Ratio (ONBR). The calculation results showed that, in normal operation conditions at full nominal thermal power of 500 kW without uncertainty parameters, the maximum fuel cladding temperature of the hottest FB was about 90.4°C, which is lower than its limit value of 103°C, the minimum DNBR was 32.0, which is much higher than the recommended value of 1.5, and the minimum ONBR was 1.43, which is higher than the recommended value of 1.4 for VVRM2 LEU fuel type. When the global and local hot channel factors were taken into account, the maximum temperature of fuel cladding at the hottest FB was about 98.4 °C, for global only, and 114.3°C, for global together with local hot channel factors. The calculation results confirm the safety operation of the designed LEU core loaded with 92 fresh VVRM2 FBs.
1. Introduction
The DNRR is a 500kW, pooltype research reactor that uses light water as a moderator and coolant. The reactor was modified and upgraded from an original 250kW TRIGA Mark II reactor built in 1963. In November 1983, the initial core of the upgraded reactor was loaded with Russian VVRM2 HEU FBs of 36% U235 enrichment. The natural convection mechanism was reinforced by installing a 2 m high “chimney” above the reactor core, which permits the reactor to operate at a nominal thermal power of 500 kW [1]. The upgraded reactor retains some structures of the original TRIGA Mark II, including the graphite reflector, the horizontal beam tubes, the thermal column, and the biological shielding. From 2005 to 2012, the project on conversion of the DNRR core from HEU to LEU fuel was implemented. This was a twophase project. In the first phase, the reactor core was partially converted in September 2007 to a mixed core of 98 irradiated HEU and 6 fresh LEU VVRM2 FBs, and in the second phase, the reactor with an entire LEU core of 92 FBs was commissioned and put into operation in early 2012. At present, the reactor is operated continuously for about 150 hrs at full power, once every 1 or 2 weeks, for radioisotope production, neutron activation analysis, basic and applied researches, and education and training [2, 3].
For neutronics design calculation of the LEU working core, the MCNPREBUS linkage system code was used to estimate the shutdown margin, the excess reactivity, the control rod’s worth, the kinetics parameters, the burnup distribution, and the power density distribution [4–6]. The power peaking factors with different positions of four shim rods inside the reactor core were also determined. The maximum power peaking factor was used for thermal hydraulic analysis in a conservative method. Because of the complicated geometry of the reactor and its FB, the DNRR was modelled by the MCNP5 code to meet the requirements of high fidelity of calculation results with 3 dimensions in space and of the continuous energy calculation library ENDF/BVII.0 [7]. The PLTEMP/ANL code has been used for steadystate thermalhydraulic analysis of the DNRR’s LEU core with 92 FAs, because of its good correspondence to both the calculation model and numerical method [8]. After modification of the code, the calculation model of the DNRR with VVRM2 fuel type, using the “chimney” and natural convection mode for heat removal, was satisfied.
Practically, the PLTEMP/ANL code was applied for steadystate thermalhydraulic analysis for fuel conversion of research reactor cores using Russian tube fueltypes of VVRKN and IRT4M. A comparison of calculation results of the Kazakhstan WWRK reactor core loaded with LEU VVRKN fuel, using the PLTEMP/ANL and RELAP5 codes, was performed, and a good agreement was found between cladding surface temperatures calculated using PLTEMP/ANL and those calculated using RELAP53D [9, 10]. The PLTEMP/ANL was also used for steadystate analysis of the Uzbekistan WWRSM and IRTSofia reactors loaded with LEU IRT4M [11, 12]. Besides, the PLTEMP/ANL code was used to evaluate the safety parameters at a steady state of the hottest FB of the newly designed 10 MW multipurpose research reactor using VVRKN fuel [13]. In addition, the PLTEMP/ANL code was also applied for modeling and simulation of coupled nuclear heat energy deposition and transfer in the HEU pintype fuel assembly, as well as for steadystate safety analysis of LEU core of the GHARR1 reactor [14, 15]. The verification and validation of the PLTEMP/ANL code were also carried out by comparison with experimental data and calculation results for thermohydraulic analysis of research reactors using MTR fuel [16].
In addition to thermohydraulic analysis of the designed 92 LEU VVRM2 core for reactor core conversion and licencing its operation, the other purpose of this work was to validate the PLTEMP/ANL V3.8 code using the DNRR core loaded with 104 HEU VVRM2 FBs. The obtained calculation results were compared with experimental data collected using a VVRM2 instrumented fuel bundle (IFB) mounted with nine incorporated thermocouples on the fuel cladding [1]. This experimental method was also applied in validation of the code for steadystate thermalhydraulic analysis of nuclear reactors with plate type fuels using experimental data from the IEAR1 reactor in Brazil [17].
2. VVRM2 Fuel Types and Reactor Core Arrangement with LEU Fuel
The standard VVRM2 HEU and LEU fuel types [1, 18] have nearly the same characteristics in annulus geometry, with three coaxial tubes (fuel elements), including 2 cylindrical inner tubes and a hexagonally shaped outermost tube. Each fuel element is composed of a fuel meat, two aluminum cladding layers, and coolant channels on both sides of the fuel cladding. Both the HEU and LEU fuel types have the same thickness and total height of fuel element, and the same length of fuel meat. However, the thickness of the fuel meat and aluminum cladding differs between HEU and LEU fuel types. The U235 mass and density of the LEU fuel are greater than those of HEU fuel. Detailed parameters of HEU and LEU VVRM2 fuel types are described in Table 1 and illustrated in Figure 1.

The working core was loaded with 92 LEU FBs and 12 beryllium rods located around the neutron trap. The beryllium rods arrangement in the core center region can be avoided to get a high power density of FBs to be loaded near the neutron trap. In addition, the beryllium layers contribute to the thermalizing process of the neutron, creating high thermal neutron flux at the neutron trap that approximates the previous HEU core or mixed HEU and LEU core, with maximum thermal neutron flux of 2.2 × 10^{13} n cm^{−2}.s^{−1}.
The reactor core is immersed in demineralized water in the reactor tank of diameter 2 m and height 6.2 m. The extracting well (chimney) of length 2 m is located at the top of the reactor core. The main characteristics and detailed geometry of the DNRR are shown in Table 1 and depicted in Figures 2 and 3(a), respectively. Compared with the HEU core and the mixed core of HEU and LEU, the new LEU core has a larger total U235 loaded mass. Consequently, the operation time can be longer.
(a)
(b)
The two dry irradiation channels at cells 7–1 and 13–2 are mainly exploited for neutron activation analysis of shortlived isotopes. The wet channel at cell 1–4, the rotary specimen with 40 irradiation holes, and the neutron trap are used mainly for production of ^{131}I and ^{32}P radioisotopes. Four horizontal beam tubes, including 1 tangential and 3 radial beam tubes, are used for nuclear data measurement and neutron basic research based on silicon and bismuth filtering techniques for creating monoenergy neutrons [2].
To control the reactor, 2 safety rods and 4 shim rods made of boron carbide (B_{4}C) and an automatic regulating rod made of stainless steel are used; these rods are symmetrically arranged inside the reactor core. During the reactor operation, two safety rods are fully withdrawn and ready to drop freely into the core by natural gravity.
The instrumented fuel bundle (IFB) mounted with nine chromecopper thermocouples was used to measure the temperatures of the fuel cladding surfaces and the inlet and outlet coolants of the reactor core. The IFB is a HEU VVRM2 fuel type with 40.2 g of U235 and considered as fresh fuel because of its limited irradiation time inside the reactor core, which is suitable for shorttime experiments only. The 5 thermocouples from T_{1} to T_{5} are located in the aluminum cladding of the hexagonal outer tube, T_{6} and T_{7} are incorporated in the aluminum cladding of the two cylindrical inner tubes, and T_{8} and T_{9} are in direct contact with the coolant to measure the inlet and outlet coolant temperatures. A detailed illustration of the IFB is shown in Figure 3(b).
3. Calculation Method
3.1. PLTEMP/ANL Code
The PLTEMP code was originally developed by Mishima [19] for thermal hydraulic analysis of the Kyoto University Research Reactor loaded with curved plate MTR FBs and operated at a thermal power level of 5 MW. The code was then improved and extended to cover thermal hydraulic analysis for research reactors that use Russian concentric tubetype FBs, such as IRT4M, VVRKN, and VVRM2, at the Argonne National Laboratory, USA. The code has the ability to solve the problem of hottest single FB or few FBs, as well as of whole core calculation. The PLTEMP/ANL3.8 code was improved by adding some new features for thermal hydraulic analysis of research reactors at steadystate condition, in which one important feature was a revision of the coolant property library tables that were made identical to those used by the PARET/ANL code [20]. Many heat transfer correlations were also added for use if the relevant problems arise. The six hot channel factors used in these versions are also applicable to natural convection. In general, the code has three main solutions for temperature profiles, which include the Broyden method [21] and the analytical method for 3layer plates (fuel meat and claddings), and the analytical method for 5layer plates (fuel meat, gas gaps, and claddings).
The PLTEMP/ANL3.8, with a “chimney” or extracting well model and Collier’s heat transfer correlation [22], is quite adequate for the DNRR, which operates at low pressure and uses natural convection for heat removal. The laminar, turbulent, and transient modes of mixed convection are included in the code as equations (1) ÷ (3), respectively:where the recommended values are as follows: CL1 = 4.0, CL2 = 0.17, CL3 = 0.33, CL4 = 0.43, CL5 = 0.25, CL6 = 0.1, RE1 = 2000, and RE2 = 2500. The subscript b refers to the bulk coolant and refers to the coolant at the wall temperature. Re is the Reynolds number (ρVD_{e}/μ), Pr is the Prandtl number (μC_{p}/k), μ_{b} is the dynamics viscosity of the bulk liquid coolant (kg/(m.s)), is the dynamics viscosity of the coolant at the wall temperature (kg/(m.s)), k_{b} is the bulk coolant conductivity, D_{e} is the hydraulic diameter (m), β is the gap of the rectangular channel or annulus (m), s is the span of the channel (m), and is the acceleration due to gravity (9.80665 m/s^{2}).
Shah’s CHF correlation [23] was implemented to make the code suitable for the DNRR thermal hydraulic calculation and analysis. The ForsterGreif correlation was used to determine the ONBR in the case of the hot channel factors for safety analysis aswhere q is the heat flux in W/m^{2}, P is the pressure of the coolant in bar, and T or ΔT (w: water and sat: saturation) is the temperature in °C.
In the PLTEMP/ANL3.8 version, the hot channel factors can be applied in thermal hydraulic analysis to estimate safety margin parameters such as the ONBR, DNBR, and flow instability ratio (FIR) [24]. When using the hot channel factors, three steps of calculation are carried out. The first step is done as a normal calculation. The second step is a repetition of the first step with an increase in the reactor power and a decrease in the reactor flow to determine uncertainty in the Nusselt number correlation. The last step applies the hot channel factors to the bulk coolant, film temperature, and cladding surface heat fluxes obtained in the previous step. The hot channel factors of the code have two parts: global and local. In the global hot channel factors, the reactor power, flow rate, and heat transfer coefficient (correlation of the Nusselt number) are presented, while the local hot factors include bulk coolant temperature rise, local coolant film temperature rise, and heat flux from the cladding surface. The combined random uncertainties and the combined systematic errors are used to estimate the local and global hot channel factors, respectively.
In each node at an axial direction of the hottest or average channel, the output of the code includes the temperature, the ONBR, and the temperature of water saturation and the ONB temperature condition of the fuel cladding. Other safety parameters, including the ONBR, DNBR, and FIR, are also calculated and shown in the output file.
3.2. Modeling
The LEU VVRM2 fuel bundle has three coaxial annular tubes (fuel elements). The outermost fuel element has a hexagonal shape and is of 32 mm, with the 35 mm being of FB pitch in width across parallel sides. The other two inner fuel elements have circular shapes of outer diameters 22 mm and 11 mm, respectively (see Figure 1). The thickness of the fuel meat of the UO_{2}–Al dispersion is 0.94 mm and that of the aluminum clad is 0.78 mm on each side. There exists a gap of about 2.5–3 mm between adjacent fuel elements for coolant flow. In the PLTEMP code, the VVRM2 fuel was modelled with three coaxial cylinders as detailed in Figure 4(b). In the calculation model, the outermost cylinder preserves the conversion area of the original hexagonal tube. The working core of 92 LEU FBs was modelled with 2 hot channels (at cell 10–5 and cell 4–5 in Figure 4(d)) and 90 average channels, with a 2 m height “chimney” also taken into account. At the design calculation stage, the channel in cell 13–2 in Figure 4(d) was considered as a fullwater channel. The 12 beryllium rods were arranged in such a way as to avoid high power density of the FBs located near the neutron trap and to increase neutron thermalization for radioisotope production.
(a)
(b)
(c)
(d)
3.3. Code Validation
Before the PLTEMP/ANL3.8 code was used for thermal hydraulic analysis of the DNRR loaded fully with LEU FBs, it was validated by comparing the analytical results with the experimental data of the working core loaded with 104 HEU FBs at the full power of 500 kW. The power peaking factors of the IFB for the PLTEMP code input were calculated by the MCNP5 code. The hottest FB was determined at cell 9–6 (see Figure 3). The IFB was then put in this position to measure temperatures using seven thermocouples, five of which (T_{1} ÷ T_{5}) were used to measure the cladding surface temperature of the outermost tube and the two others (T_{8} and T_{9}) to measure the coolant temperature at the entrance and exit of the flow channel, respectively. The error of this measurement system is about ±3°C. The comparative results of the calculated and measured data in Figure 5 show that the discrepancy between calculated results and experimental data is not significant. Except for the difference of about 1.5°C at the outlet coolant channel, the difference in other temperatures was below 1°C and fully addressed by the error of the measurement system. A good agreement between calculation results and measurement data for DNRR’s HEU core allows the PLTEMP/ANL3.8 code to be used “adequately” for its LEU core due to the very similar geometry, power level, and coolant flow rate.
3.4. Power Density Distribution of the LEU Core
The MCNP5 code was used to calculate power distribution with 30 nodes in the axial direction of the hottest FB at cell 10–5 (see Figure 4(d)). The average power distribution in the radial direction of each FB was estimated, together with local power peaking of each cylindrical and hexagonal tube. Figure 6 shows the axial relative power distributions (axial power peaking factors) of the hottest FB by different insertions of the shim rods. The calculated results show that the maximum axial power peaking factor was 1.37, corresponding to the 25 cm insertion of 4 shim rods. The radial power peaking factors were 1.41 and 1.38 for the reactor core and the FB in the local hot channel, respectively. The most powerful FB was found to be next to a beryllium layer around the neutron trap (cell 4–5 and cell 10–5 in Figure 4(d)). The details of the power peaking factor of the hottest FB, according the control rod positions, were also obtained. Two of the six sides of the hottest FB conduct no heat, because this FB is located near the beryllium rod around the neutron trap.
(a)
(b)
3.5. Hot Channel Factors of the LEU Core
There are six hot channel factors of the PLTEMP code: the reactor power, the coolant flow, and the Nusselt number correlation in the global part and the bulk coolant temperature, the coolant film temperature, and the heat flux from the cladding surface in the local part. These factors are used to calculate the ONBR, the DNBR, and the maximum fuel cladding and coolant temperatures. The separate calculation of the three global factors is available. Table 2 presents the hot channel factors of the DNRR LEU core prepared for the PLTEMP/ANL code according to the method proposed by Feldman [24]. The first three factors were statistically combined from the random uncertainties of the fuel meat thickness, the U235 homogeneity, the U235 mass per fuel element, the power density, the cooling channel spacing, and the flow distribution, as given by the VVRM2 fuel catalogue [18]. This statistical combination was used for the random uncertainties because it was unlikely that all of these sources could adversely affect the limiting value of the core simultaneously. Three other factors are the systematic uncertainties of the reactor power measurement, coolant flow measurement, and heat transfer coefficient, which affect all regions of the reactor core equally. The combined factor of the random uncertainties [24] can be calculated according to and the statistic combination errors [24] can be determined by multiplying all factors with by

4. Results and Discussion
Steadystate thermalhydraulic calculation for the DNRR loaded with LEU fuel was implemented by the PLTEMP/ANL in three steps. At the beginning, a normal calculation without hot channel factors at power levels of 500 kW, 550 kW, and 600 kW together with the inlet coolant temperature of 32°C was carried out. Then, the calculation was performed with global hot channel factors, and the last calculation was carried out by using all global and local hot channel factors.
4.1. Without Hot Channel Factors
The bestestimate calculation results of the DNRR are shown in Figure 7. The maximum cladding surface temperature was determined to be 90.4°C, which was far below the estimated ONB temperature of 115.8°C using the ForsterGreif correlation. It met the requirement that the maximum fuel cladding temperature must be lower than 103°C at the nominal power as recommended by the fuel supplier [18]. Using the maximum permissible inlet temperature of 32 °C, the maximum coolant temperature at the core exit was calculated to be 57.7°C. This value was very low compared to the saturation temperature of 108°C. For a thin layer of the fuel element tube, the temperature of the fuel meat was only slightly higher than that of the cladding, with the maximum difference only 0.4°C. The calculated results showed that the minimum DNBR was 32.0, the minimum ONBR was 1.43, and the total coolant flow rate through the reactor core was about 26.3 m^{3}/h. The calculation results were in good confirmation for the safe operation of the DNRR when using LEU fuel.
If the reactor power increased up to 550 kW and 600 kW, we can see from the calculation results that the maximum fuel cladding surface temperature would rise to 94.2 and 97.7°C, corresponding with the maximum coolant temperature of 59.2 and 60.0°C, respectively. Figure 8 shows the calculation results of the fuel cladding and coolant temperatures at different reactor power levels.
The maximum fuel cladding surface temperature at the thermal power of 600 kW was still below its limitation temperature of 103°C. In practice, however, the maximum operational power is about 550 kW, corresponding to 110% of the nominal reactor power, so the fuel cladding temperature still remained less than 95°C. According to the recorded technological operation data in the operating logbook, the uncertainty in the power is only 5%, which meets the requirement of the maximum fuel cladding temperature of less than 103°C.
4.2. With Global Hot Channel Factors
Taking into account the systematic uncertainties in the global hot channel factors increased, the maximum cladding surface temperature became 98.4°C, which was well below the maximum fuel cladding temperature of 103°C and far below the estimated ONB temperature of 115.9°C as shown in Figure 9. All global hot channel factors of the DNRR can be determined through experiments or operational data. Based on the flow rate of the primary and secondary cooling systems and light water properties, the thermal power of the reactor was estimated and recorded during normal operation for each hour. The total data were analyzed to validate the PLTEMP/ANL thermalhydraulic analysis code. When the DNBR value was 17.79, the ONBR value was 1.26, the total coolant flow rate through the reactor core was 26.3 m^{3}/h, and the coolant temperature of 58.8°C at the outlet of the hottest FB was still extremely low compared with the saturation temperature of 107°C. Yet, the DNRR still met thermalhydraulic safety at the steadystate condition.
4.3. With all Global and Local Hot Channel Factors
Calculation results for nominal power with all systematic and random uncertainties taken into account are shown in Figure 10. In this calculation step, hot channel factors were applied to the bulk coolant, local coolant film temperature rises, and the cladding surface heat fluxes were obtained in the previous step. The maximum fuel cladding temperature was 114.3°C, which was only a few degrees less than the predicted ONB temperature of 116.4°C. Using Shah’s correlation, the obtained minimum DNBR value of 15.2 was far above the limit value of 1.5 specified for the VVRM2 LEU fuel of the DNRR [1, 18]. The maximum coolant temperature with all global and local hot channel factors considered was about 68°C, which was lower than the saturation value of 107°C. The ONBR value in this case was only 1.02 and obtained nearly to the boiling point. However, the integrity of the fuel cladding was still confirmed because the melting point of the aluminum material is about 650°C.
5. Conclusions
The validation of the PLTEMP/ANLV3.8 code was carried out by comparing the calculation results of the designed LEU core to the experimental data of the HEU core of the DNRR with extremely good consistency. The maximum discrepancy in the fuel cladding temperatures was about 0.8°C, as shown in Figure 5. This implies that this code can be applied to analyze steadystate thermalhydraulics for the DNRR’s LEU core with annulus fuel VVRM2 tubetype and natural convection, with an extracting well at the core top.
The designed working core with 92 LEU FBs and 12 beryllium rods was evaluated and subjected to steadystate thermal hydraulic analysis without hot channel factors at the nominal thermal power of 500 kW. The obtained results show that, at the hottest FB, the fuel cladding temperature was only 90.4°C, as shown in Figure 7, which was far below the limitation of the VVRM2 fuel cladding temperature of 103°C, the ONB temperature was about 115.8°C, according to the ForsterGreif correlation, the minimum DNBR value was about 32.0, and the minimum ONBR was about 1.43. Thus, the operation condition in this case met all the requirements for thermalhydraulic safety.
When the systematic uncertainties were taken into account, the maximum fuel cladding temperature was predicted to be 98.4°C, as shown in Figure 9, which again was well below the limit value of 103°C. The ONB temperature was about 115.9°C, the minimum DNBR value was 17.79, and the minimum ONBR was about 1.26, which was below recommended value of 1.4. However, the coolant temperature at the outlet of the hottest FB was about 58.8°C, which was extremely low compared with the saturation temperature of 107°C. This means that the DNRR still met the thermal hydraulic safety at steadystate condition with application of the global hot channel factors.
When all systematic as well as random uncertainties were applied to a limiting calculation, the maximum fuel cladding temperature obtained was 114.3°C, which was several degrees below the ONB point of 116.4°C. The minimum DNBR value, according to Shah’s correlation, was estimated as 15.2, which was much higher than the acceptable criterion of 1.5, but ONBR value was only 1.02 and obtained near to the boiling point for the DNRR.
From the above results we can conclude that the DNRR core loaded with VVRM2 LEU FBs is safe at the steadystate condition, from a thermalhydraulic point of view in cases without hot channel factors and with global hot channel factors, and the reactor can be safely operated within the thermal design limit and nominal power of 500 kW, with natural convection adopted for heat removal.
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 staff of DNRI for their kind support. The staff at Reactor Physics and Engineering Department of DNRI are acknowledged for their valuable technical discussions. Collaboration with the staff of Argonne National Laboratory is greatly appreciated. This research was supported by Ministry of Science and Technology of Vietnam under Grant No. DTCB.10/19/VNCHN.
References
 N. D. Nguyen, Safety Analysis Report for the Dalat Nuclear Research Reactor, Nuclear Research Institute, Vietnam Atomic Energy Institute, Hanoi, Vietnam, 2009.
 N. D. Nguyen, L. B. Vien, V. Le et al., “Results of operation and utilization of the Dalat nuclear research reactor,” Nuclear Science and Technology, vol. 4, no. 1, Vietnam Atomic Energy Association, 2014. View at: Google Scholar
 G. T. T. 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 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
 X5 Monte Carlo Team, MCNP  A General Monte Carlo NParticle Transport Code. Version 5. LAUR031987, Los Alamos National Security LLC, Los Alamos, NM, USA, 2003.
 J. G. Stevens, The RebusMCNP Linkage, Argonne National Laboratory, Lemont, IL, USA, 2008.
 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, no. 12, pp. 2931–3060, 2006. View at: Publisher Site  Google Scholar
 A. P. Olson and M. Kalimullah, “A user’s guide to the PLTEMP/ANL V3.8 code,” Tech. Rep., ANL/RERTR, Argonne National Laboratory, Lemont, IL, USA, 2009, TRN: US1600226. View at: 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: Google Scholar
 S. N. Koltochnik and A. A. Shaimerdenov, “Safety analysis for the WWRK research reactor converted to LEU fuel,” Eurasian Journal of Physics and Functional Materials, vol. 3, no. 3, pp. 204–218, 2019. View at: Publisher Site  Google Scholar
 P. L. Garner and N. K. Hanan, “Transient analyses for the Uzbekistan VVRSM reactor with IRT3M HEU fuel and IRT4M LEU fuel. ANL independent verification results,” Tech. Rep., Argonne National Laboratory, Lemont, IL, USA, 2007, TRN: US0803715. View at: Google Scholar
 S. Kadalev, “Thermal hydraulic calculations of the IRTSOFIA research reactor with LEU Fuel IRT4M, Sofia,” Comptes rendus de l’Académie bulgare des sciences: sciences mathématiques et naturelles, vol. 65, no. 9, pp. 1185–1192, 2012. View at: Google Scholar
 N.D. Nguyen, K.C. Nguyen, T.N. Huynh, D.H. Dang Vo, and H.N. 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. View at: Publisher Site  Google Scholar
 F. Ameyaw, A. Ayensu, and E. H. K. Akaho, “Modeling and simulation of coupled nuclear heat energy deposition and transfer in the fuel assembly of the Ghana Research Reactor1 (GHARR1),” Nuclear Engineering and Design, vol. 241, no. 12, pp. 5183–5188, 2011. View at: Publisher Site  Google Scholar
 E. Shitsi, P. Amoah, E. AmpomahAmoako, and H. C. Odoi, “Steadystate safety analysis of Ghana research reactor1 with lowenricheduranium core,” Journal of Thermal Science and Engineering Applications, vol. 12, no. 5, pp. 1–18, 2020. View at: Publisher Site  Google Scholar
 M. Kalimullah, A. P. Olson, E. E. Feldman et al., “Verification and validation of the PLTEMP/ANL code for thermalhydraulics analysis of experimental and test reactors,” Tech. Rep., Argonne National Laboratory, Lemont, IL, USA, 2015, TRN: 161160. View at: Google Scholar
 D. A. CastellanosGonzalez, J. M. Losada Moreira, J. Rubens Maiorino, and C. Pedro, “Validation of the cotenp code: a steadystate thermalhydraulic analysis code for nuclear reactors with plate type fuel assemblies,” Science and Technology of Nuclear Installations, vol. 2018, Article ID 9874196, 17 pages, 2018. View at: Publisher Site  Google Scholar
 VVRM2 and VVRM5 Fuel Bundles  Operation Manual”, 0001.04.00.000 PЭ, 2006.
 K. Mishima, K. Kanda, and T. Shibata, Thermalhydraulic Analysis for Core Conversion to the Use of LowEnrichment Uranium Fuels in KUR, KURRITR258, Research Reactor Institute, Kyoto University, Kyoto, Japan, 1984.
 W. L. Woodruff and R. S. Smith, A Users Guide to PARET/ANL, Version 7.5, Argonne National Laboratory, Lemont, IL, USA, 2006.
 H. William, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in FORTRAN, the Art of Scientific Computing, Cambridge University Press, Cambridge, England, 2nd edition, 1992.
 J. G. Collier and J. R. Thome, Convective Boiling and Condensation, Clarendon Press, Oxford, UK, 3rd edition, 1994.
 M. M. Shah, “Improved general correlation for critical heat flux during upflow in uniformly heated vertical tubes,” International Journal of Heat and Fluid Flow, vol. 8, no. 4, pp. 326–335, 1987. View at: Publisher Site  Google Scholar
 E. E. Feldman, Recommended Treatment of Hot Channel Factors in the PLTEMP Code, IntraLaboratory Memorandum, Argonne National Laboratory, Lemont, IL, USA, 2006.
Copyright
Copyright © 2021 KienCuong 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.