• Journal of computational fluids engineering
• /
• v.20 no.4
• /
• pp.14-20
• /
• 2015

Without electirc power to cool down the hot reactor core, passive systems utilizing natural circulation are becoming a big specialty of recent neculear systems after the severe accident in Fukusima. When we consider the natural circulation in a pool, thermal mixing phenomena may start from single phase circulation and can continue to two phase condition. Since the CUPID code, which has been developed for two-phase flow analysis, can deal with the phase transition phenomena, the CUPID would be pertinent to natural convection problems in single- and two-phase conditions. Thus, the CUPID should be validated against single- and two-phase natural circulation phenomena. For the first step of the validation process, this study is focused on the validation of single-phase natural circulation. Moreover, the CUPID code solves the fluid properties by the relationship to pressure and temperature from the steam table considering non-condensable gas effects, so that the effects from variable properties are included. Simple square thermal cavity problems are tested for laminar and turbulent conditions against numerical and experimental data. Throughout the investigation, it is found that the variable properties can affect the flow field in laminar condition, but the effect becomes weak in turbulence condition, and the CUPID code implementing steam table is capable of analyzing single phase natural circualtion phenomena.

• Proceedings of the Korean Nuclear Society Conference
• /
• 2001.10a
• /
• pp.143.2-143
• /
• 2001

• Proceedings of the Korean Nuclear Society Conference
• /
• 2002.10a
• /
• pp.186.2-186
• /
• 2002

• Nuclear Engineering and Technology
• /
• v.55 no.5
• /
• pp.1911-1923
• /
• 2023

The nuclear safety assessment involving large transient simulations is forcing the community to develop methods for coupling thermal-hydraulics and neutronic codes and three-dimensional (3D) Computational Fluid Dynamics (CFD) codes. In this paper a set of dynamic boundary conditions are implemented in OpenFOAM^Ⓡ in order to apply zero-dimensional (0D) approaches coupling with 3D thermal-hydraulic simulation in a single framework. This boundary conditions are applied to model pipelines, tanks, pumps, and heat exchangers. On a first stage, four tests are perform in order to assess the implementations. The results are compared with experimental data, full 3D CFD, and system code simulations, finding a general good agreement. The semi-implicit implementation nature of these boundary conditions has shown robustness and accuracy for large time steps. Finally, an application case, consisting of a simplified open pool with a cooling external circuit is solved to remark the capability of the tool to simulate thermal hydraulic systems commonly found in nuclear installations.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.258-270
• /
• 2020

A new task of using the Jacobian-Free-Newton-Krylov (JFNK) method for the PWR core transient simulations involving neutronics, thermal hydraulics and mechanics is conducted. For the transient scenario of PWR, normally the Picard iteration of the coupled coarse-mesh nodal equations and parallel channel TH equations is performed to get the transient solution. In order to solve the coupled equations faster and more stable, the Newton Krylov (NK) method based on the explicit matrix was studied. However, the NK method is hard to be extended to the cases with more physics phenomenon coupled, thus the JFNK based iteration scheme is developed for the nodal method and parallel-channel TH method. The local gap conductance is sensitive to the gap width and will influence the temperature distribution in the fuel rod significantly. To further consider the local gap conductance during the transient scenario, a 1D mechanics model is coupled into the JFNK scheme to account for the fuel thermal expansion effect. To improve the efficiency, the physics-based precondition and scaling technique are developed for the JFNK iteration. Numerical tests show good convergence behavior of the iterations and demonstrate the influence of the fuel thermal expansion effect during the rod ejection problems.

• Nuclear Engineering and Technology
• /
• v.50 no.3
• /
• pp.356-367
• /
• 2018

A series of tests dedicated to station blackout (SBO) accident scenarios have been recently performed at the $Prim{\ddot{a}}rkreislauf-Versuchsanlage$ (primary coolant loop test facility; PKL) facility in the framework of the OECD/NEA PKL-3 project. These investigations address current safety issues related to beyond design basis accident transients with significant core heat up. This work presents a detailed analysis using the best estimate thermal-hydraulic code TRACE (v5.0 Patch4) of different SBO scenarios conducted at the PKL facility; failures of high- and low-pressure safety injection systems together with steam generator (SG) feedwater supply are considered, thus calling for adequate accident management actions and timely implementation of alternative emergency cooling procedures to prevent core meltdown. The presented analysis evaluates the capability of the applied TRACE model of the PKL facility to correctly capture the sequences of events in the different SBO scenarios, namely the SBO tests H2.1, H2.2 run 1 and H2.2 run 2, including symmetric or asymmetric secondary side depressurization, primary side depressurization, accumulator (ACC) injection in the cold legs and secondary side feeding with mobile pump and/or primary side emergency core coolant injection from the fuel pool cooling pump. This study is focused specifically on the prediction of the core exit temperature, which drives the execution of the most relevant accident management actions. This work presents, in particular, the key improvements made to the TRACE model that helped to improve the code predictions, including the modeling of dynamical heat losses, the nodalization of SGs' heat exchanger tubes and the ACCs. Another relevant aspect of this work is to evaluate how well the model simulations of the three different scenarios qualitatively and quantitatively capture the trends and results exhibited by the actual experiments. For instance, how the number of SGs considered for secondary side depressurization affects the heat transfer from primary side; how the discharge capacity of the pressurizer relief valve affects the dynamics of the transient; how ACC initial pressure and nitrogen release affect the grace time between ACC injection and subsequent core heat up; and how well the alternative feeding modes of the secondary and/or primary side with mobile injection pumps affect core quenching and ensure stable long-term core cooling under controlled boiling conditions.

• Nuclear Engineering and Technology
• /
• v.41 no.4
• /
• pp.391-402
• /
• 2009

This paper introduces the development of safety assessment technology in Korea, focusing on the activities of the Korea Atomic Energy Research Institute in the areas of system thermal hydraulics, severe accidents and probabilistic safety assessment. In the 1970s and 1980s, safety analysis codes and methodologies were introduced from the United States, France, Canada and other developed countries along with technology related to the construction and operation of nuclear power plants. The main focus was on understanding and utilizing computer codes that were sourced from abroad up to the early 1990s, when efforts to develop domestic safety analysis codes and methodologies became active. Remarkable achievements have been made over the last 15 years in the development and application of safety analysis technologies. In addition, significant experimental work has been performed to verify the safety characteristics of reactors and fuels as well as to support the development and validation of analysis methods.

• Nuclear Engineering and Technology
• /
• v.42 no.3
• /
• pp.279-296
• /
• 2010

For the analysis of transient two-phase flows in nuclear reactor components, a three-dimensional thermal hydraulics code, named CUPID, has been developed. The CUPID code adopts a two-fluid, three-field model for two-phase flows, and the governing equations were solved over unstructured grids, which are very useful for the analysis of flows in complicated geometries. To obtain numerical solutions, the semi-implicit numerical method for the REALP5 code was modified for an application to unstructured grids, and it has been further improved for enhanced accuracy and fast running. For the verification of the CUPID code, a set of conceptual problems and experiments were simulated. This paper presents the flow model, the numerical solution method, and the results of the preliminary assessment.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1219-1225
• /
• 2017

In this work, a scalable algorithm for model calibration in nuclear engineering applications is presented and tested. The algorithm relies on the construction of surrogate models to replace the original model within the region of interest. These surrogate models can be constructed efficiently via reduced order modeling and subspace analysis. Once constructed, these surrogate models can be used to perform computationally expensive mathematical analyses. This work proposes a surrogate based model calibration algorithm. The proposed algorithm is used to calibrate various neutronics and thermal-hydraulics parameters. The virtual environment for reactor applications-core simulator (VERA-CS) is used to simulate a three-dimensional core depletion problem. The proposed algorithm is then used to construct a reduced order model (a surrogate) which is then used in a Bayesian approach to calibrate the neutronics and thermal-hydraulics parameters. The algorithm is tested and the benefits of data assimilation and calibration are highlighted in an uncertainty quantification study and requantification after the calibration process. Results showed that the proposed algorithm could help to reduce the uncertainty in key reactor attributes based on experimental and operational data.

• Nuclear Engineering and Technology
• /
• v.50 no.3
• /
• pp.379-388
• /
• 2018

Models for the rate of atomization and deposition of droplets for stratified and annular flow in horizontal pipes are presented. The entrained fraction is the result of a balance between the rate of atomization of the liquid layer that is in contact with air and the rate of deposition of droplets. The rate of deposition is strongly affected by gravity in horizontal pipes. The gravitational settling of droplets is influenced by droplet size: heavier droplets deposit more rapidly. Model calculation and simulation results are compared with experimental data from various diameter pipes. Validation for the suggested models was performed by comparing the Safety and Performance Analysis Code for Nuclear Power Plants calculation results with the droplet experimental data obtained in various diameter horizontal pipes.