• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.71-77
• /
• 2016

The CUPID code is a transient, three-dimensional, two-fluid, thermal-hydraulic code designed for a component-scale analysis of nuclear reactor components. The primary objective of this study is to assess the applicability of CUPID to single-phase turbulent flow analyses of $2{\times}2$ rod bundle subchannel. The bulk velocity at the inlet varies from 1.0 m/s up to 2.0 m/s which is equivalent to the fully turbulent flow with the range of Re=12,500 to 25,000. Adiabatic single-phase flow is assumed. The velocity profile at the exit region is quantitatively compared with both experimental measurement and commercial CFD tool. Three different boundary conditions are simulated and quantitatively compared each other. The calculation results of CUPID code shows a good agreement with the experimental data. It is concluded that the CUPID code has capability to reproduce the turbulent flow behavior for the $2{\times}2$ rod bundle geometry.

• Nuclear Engineering and Technology
• /
• v.52 no.12
• /
• pp.2743-2759
• /
• 2020

A detailed computational fluid dynamics (CFD) simulation analysis model was developed using ANSYS CFX 16.1 and analyzed to simulate the basic design and internal flow characteristics of a 180 MW small modular reactor (SMR) with a natural circulation flow system. To analyze the natural circulation phenomena without a pump for the initial flow generation inside the reactor, the flow characteristics were evaluated for each output assuming various initial powers relative to the critical condition. The eddy phenomenon and the flow imbalance phenomenon at each output were confirmed, and a flow leveling structure under the core was proposed for an optimization of the internal natural circulation flow. In the steady-state analysis, the temperature distribution and heat transfer speed at each position considering an increase in the output power of the core were calculated, and the conceptual design of the SMR had a sufficient thermal margin (31.4 K). A transient model with the output ranging from 0% to 100% was analyzed, and the obtained values were close to the T_hot and T_cold temperature difference value estimated in the conceptual design of the SMR. The K-factor was calculated from the flow analysis data of the CFX model and applied to an analysis model in RELAP5/MOD3.3, the optimal analysis system code for nuclear power plants. The CFX analysis results and RELAP analysis results were evaluated in terms of the internal flow characteristics per core output. The two codes, which model the same nuclear power plant, have different flow analysis schemes but can be used complementarily. In particular, it will be useful to carry out detailed studies of the timing of the steam generator intervention when an SMR is activated. The thermal and hydraulic characteristics of the models that applied porous media to the core & steam generators and the models that embodied the entire detail shape were compared and analyzed. Although there were differences in the ability to analyze detailed flow characteristics at some low powers, it was confirmed that there was no significant difference in the thermal hydraulic characteristics' analysis of the SMR system's conceptual design.

• Nuclear Engineering and Technology
• /
• v.35 no.1
• /
• pp.45-54
• /
• 2003

A hyperbolic two-phase, two-fluid equation system developed in the previous work has been implemented in an existing nuclear safety analysis code, MARS. Although the implicit treatment of interfacial pressure force term introduced in momentum equation of the hyperbolic equation system is required to enhance the numerical stability, it is very difficult to implement in the code because it is not possible to maintain the existing numerical solution structure. As an alternative, two-step approach with stabilizer momentum equations has been selected. The results of a linear stability analysis by Von-Neumann method show the equivalent stability improvement with fully-implicit solution method. To illustrate the applicability, the new solution scheme has been implemented into the best-estimate thermal-hydraulic analysis code, MARS. This paper also includes the comparisons of the simulation results for the perturbation propagation and water faucet problems using both two-step method and the original solution scheme.

• Nuclear Engineering and Technology
• /
• v.41 no.10
• /
• pp.1347-1360
• /
• 2009

A multi-dimensional component for the thermal-hydraulic system analysis code, MARS, was developed for a more realistic three-dimensional analysis of nuclear systems. A three-dimensional and two-fluid model for a two-phase flow in Cartesian and cylindrical coordinates was employed. The governing equations and physical constitutive relationships were extended from those of a one-dimensional version. The numerical solution method adopted a semi-implicit and finite-difference method based on a staggered-grid mesh and a donor-cell scheme. The relevant length scale was very coarse compared to commercial computational fluid dynamics tools. Thus a simple Prandtl's mixing length turbulence model was applied to interpret the turbulent induced momentum and energy diffusivity. Non drag interfacial forces were not considered as in the general nuclear system codes. Several conceptual cases with analytic solutions were chosen and analyzed to assess the fundamental terms. RPI air-water and UPTF 7 tests were simulated and compared to the experimental data. The simulation results for the RPI air-water two-phase flow experiment showed good agreement with the measured void fraction. The simulation results for the UPTF downcomer test 7 were compared to the experiment data and the results from other multi-dimensional system codes for the ECC delivery flow.

• Journal of the Korean Society of Safety
• /
• v.27 no.5
• /
• pp.229-234
• /
• 2012

Operator's action time is evaluated from MAAP4 analysis used in conventional probabilistic safety assessment(PSA) of a nuclear power plant. MAAP4 code which was developed for severe accident analysis is too conservative to perform a realistic PSA. A best-estimate code such as RELAP5/MOD3, MARS has been used to reduce the conservatism of thermal hydraulic analysis. In this study, operator's action time of core cooling recovery operation is evaluated by using the MARS code, which its Fussell-Vessely(F-V) value was evaluated as highly important in a small break loss of coolant(SBLOCA) event and loss of component cooling water(LOCCW) event in previous PSA. The main conclusions were elicited : (1) MARS analysis provides larger time window for operator's action time than MAAP4 analysis and gives the more realistic time window in PSA (2) Sufficient operator's action time can reduce human error probability and core damage frequency in PSA.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05b
• /
• pp.379-384
• /
• 1996

Thermal-hydraulic performance has been analyzed for an advanced reactor core loaded with hexagonal fuel assemblies. Currently available CHF prediction models and data base for triangular lattice bundles have been thoroughly reviewed, and as a result the KfK-3 CHF correlation with limit CHFR of 1.235 has been determined to be most appropriate. The pressure drop model in COBRA-IV-I code has been modified for the analysis of triangular lattice rod bundles. In view of maximizing the thermal margin, the geometry of a hexagonal fuel assembly, such as rod diameter and rod pitch, has been optimized with a fixed fuel assembly cross sectional area The optimum value of the moderator-to-fuel volume ratio is estimated to lie between 0.65 to 1 with 9.5 mm rod diameter. The thermal margin of these hexagonal fuel assemblies in the AP600 core has been evaluated and compared with that of square lattice fuel assemblies such as VANTAGE-5H and KOFA. The analysis result shows that the performances of hexagonal fuel assemblies are more favorable than the square fuel assemblies in the aspect of steady-state overpower margin.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3388-3400
• /
• 2023

Detailed analysis of the neutron pathway through matter inside the nuclear reactor core is exceedingly needed for safety and economic considerations. Due to the constant development of high-performance computing technologies, neutronics analysis using computer codes became more effective and efficient to perform sophisticated neutronics calculations. In this work, a commercial pressurized water reactor (PWR) presented by Virtual Environment for Reactor Applications (VERA) Core Physics Benchmark are modeled and simulated using a high-fidelity simulation of OpenMC code in terms of criticality and fuel pin power distribution. Various problems have been selected from VERA benchmark ranging from a simple two-dimension (2D) pin cell problem to a complex three dimension (3D) full core problem. The development of the code capabilities for reactor physics methods has been implemented to investigate the accuracy and performance of the OpenMC code against VERA SCALE codes. The results of OpenMC code exhibit excellent agreement with VERA results with maximum Root Mean Square Error (RMSE) values of less than 0.04% and 1.3% for the criticality eigenvalues and pin power distributions, respectively. This demonstrates the successful utilization of the OpenMC code as a simulation tool for a whole core analysis. Further works are undergoing on the accuracy of OpenMC simulations for the impact of different fuel types and burnup levels and the analysis of the transient behavior and coupled thermal hydraulic feedback.

• Journal of Mechanical Science and Technology
• /
• v.17 no.9
• /
• pp.1380-1387
• /
• 2003

When a cold HPSI (High pressure Safety Injection) fluid associated with an overcooling transient, such as SGTR (Steam Generator Tube Rupture), MSLB (Main Steam Line Break) etc., enters the cold legs of a stagnated primary coolant loop, thermal stratification phenomena will arise due to incomplete mixing. If the stratified flow enters the downcomer of the reactor pressure vessel, severe thermal stresses are created in a radiation embrittled vessel wall by local overcooling. As general thermal-hydraulic system analysis codes cannot properly predict the thermal stratification phenomena, RG 1.154 requires that a detailed thermal-mixing analysis of PTS (pressurized Thermal Shock) evaluation be performed. Also. previous PTS studies have assumed that the thermal stratification phenomena generated in the stagnated loop side of a partially stagnated primary coolant loop are neutralized in the vessel downcomer by the strong flow from the unstagnated loop. On the basis of these reasons, this paper focuses on the development of a 3-dimensional thermal-mixing analysis model using PHOENICS code which can be applied to both partial and total loop stagnated cases. In addition, this paper verifies the fact that, for partial loop stagnated cases, the cold plume generated in the vessel downcomer due to the thermal stratification phenomena of the stagnated loop is almost neutralized by the strong flow of the unstagnated loop but is not fully eliminated.

• Nuclear Engineering and Technology
• /
• v.56 no.2
• /
• pp.762-769
• /
• 2024

The TMSR-SF0 simulator is an integral effect thermal-hydraulic experimental system for the development of thorium molten salt reactor (TMSR) program in China. The simulator has two heat transport loops with liquid FLiNaK. In literature, the 95% level confidence uncertainties of the thermophysical properties of FLiNaK are recommended, and the uncertainties of density, heat capacity, thermal conductivity and viscosity are ±2%, ±10, ±10% and ±10% respectively. In order to investigate the effects of thermophysical properties uncertainties on the molten salt heat transport system, the uncertainty and sensitivity analysis of the heat transfer characteristics of the simulator system are carried out on a RELAP5 model. The uncertainties of thermophysical properties are incorporated in simulation model and the Monte Carlo sampling method is used to propagate the input uncertainties through the model. The simulation results indicate that the uncertainty propagated to core outlet temperature is about ±10 ℃ with a confidence level of 95% in a steady-state operation condition. The result should be noted in the design, operation and code validation of molten salt reactor. In addition, more experimental data is necessary for quantifying the uncertainty of thermophysical properties of molten salts.

• Journal of the Korea Society for Simulation
• /
• v.13 no.2
• /
• pp.13-21
• /
• 2004

The thermal-hydraulic model ARTS which was based on the RETRAN-3D code adopted in the domestic full-scope power plant simulator which was provided in 1998 by KEPRI. Since ARTS is a generalized code to model the components with control volumes, the smaller time-step size should be used even if converged solution could not get in a single volume. Therefore, dedicated models which do not force to reduce the time-step size are sometimes more suitable in terms of a real-time calculation and robustness. In the case of PRT(Pressurizer Relief Tank) model, it is consist of subcooled water in bottom and non-condensable gas in top. The sparger merged under subcooled water enhances condensation. The complicated thermal-hydraulic phenomena such as condensation, phase separation with existence of non-condensable gas makes difficult to simulate. Therefore, the PRT volume can limit the time-step size if we model it with a general control volume. To prevent the time-step size reduction due to convergence failure for simulating this component, we developed a dedicated model for PRT. The dedicated model was expected to provide substantially more accurate predictions in the analysis of the system transients. The results were resonable in terms of accuracy, real-time simulation, robustness and education of operators, complying with the ANSI/ANS-3.5-1998 simulator software performance criteria and RETRAN-3D results.