• Nuclear Engineering and Technology
• /
• v.32 no.1
• /
• pp.46-56
• /
• 2000

Numerical results are presented for the 2-dimensional turbulent transient heat transfer of the shell/tube heat exchanger with a step change of the inlet temperature in the primary side. Heat transfer boundary conditions outside the pipe are given partially by the convection heat transfer conditions and partially by insulated conditions. Calculation results were obtained by solving the unsteady two-dimensional elliptic forms for the Reynolds-averaged governing equations for the mass, momentum and energy. Finite-difference method was used to obtain discretization equations, and the SIMPLER solution algorithm was employed for the calculation procedure. Turbulent model used is the algebraic model proposed by Cebeci-Smith. Results presented include the time variant Nusselt number distribution, average temperature distribution and outlet temperatures for the various inlet temperatures and flow rates.

• Nuclear Engineering and Technology
• /
• v.50 no.7
• /
• pp.1190-1197
• /
• 2018

Analysis of the GRSW-A model coupled into the FALCON code is extended by simulation of central void formation in fuel pellets due to high-temperature fuel restructuring. The extended calculation is verified against published, well-known experimental data. Good agreement with the data for a central void diameter in pellets of the rod irradiated in an Experimental Breeder Reactor is shown. The new calculation methodology is employed in comparative analysis of modern BWR fuel behavior under assumed high-power operation. The initial fuel porosity is shown to have a major effect on the predicted central void diameter during the operation in question. Discernible effects of a central void on peak fuel temperature and Pellet-Cladding Mechanical Interaction (PCMI) during a simulated power ramp are shown. A mitigating effect on PCMI is largely attributed to the additional free volume in the pellets into which the fuel can creep due to internal compressive stresses during a power ramp.

• Journal of Radiation Protection and Research
• /
• v.41 no.2
• /
• pp.155-159
• /
• 2016

Background: Nuclear reactors produce a great number of antielectron neutrinos mainly from beta-decay chains of fission products. Such neutrinos have energies mostly in MeV range. We are interested in neutrinos in a region of keV, since they may take part in special weak interactions. We calculate reactor antineutrino spectra especially in the low energy region. In this work we present neutrino spectrum from a typical pressurized water reactor (PWR) reactor core. Materials and Methods: To calculate neutrino spectra, we need information about all generated nuclides that emit neutrinos. They are mainly fission fragments, reaction products and trans-uranium nuclides that undergo negative beta decay. Information in relation to trans-uranium nuclide compositions and its evolution in time (burn-up process) were provided by a reactor code MVP-BURN. We used typical PWR parameter input for MVP-BURN code and assumed the reactor to be operated continuously for 1 year (12 months) in a steady thermal power (3.4 GWth). The PWR has three fuel compositions of 2.0, 3.5 and 4.1 wt% $^{235}U$ contents. For preliminary calculation we adopted a standard burn-up chain model provided by MVP-BURN. The chain model treated 21 heavy nuclides and 50 fission products. The MVB-BURN code utilized JENDL 3.3 as nuclear data library. Results and Discussion: We confirm that the antielectron neutrino flux in the low energy region increases with burn-up of nuclear fuel. The antielectron-neutrino spectrum in low energy region is influenced by beta emitter nuclides with low Q value in beta decay (e.g. $^{241}Pu$) which is influenced by burp-up level: Low energy antielectron-neutrino spectra or emission rates increase when beta emitters with low Q value in beta decay accumulate Conclusion: Our result shows the flux of low energy reactor neutrinos increases with burn-up of nuclear fuel.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1998.05a
• /
• pp.635-640
• /
• 1998

Thermal fluid mixing behavior during a postulated medium-size hot leg break loss of coolant accident is analyzed for the international comparative assessment study on pressurized thermal shock (PTS-ICAS) proposed by OECD-NEA. The applicability of RELAP5 code to analyze thermal fluid mixing behavior is evaluated through a simple modeling relevant to the problem constraints. Based on the calculation result, the onset of Thermal stratification is investigated using Theofanous's empirical correlation. Sensitivity calculations using a fine node model and crossflow model are also performed to evaluate the modeling capability on multi-dimensional characteristics related to thermal fluid mixing.

• Nuclear Engineering and Technology
• /
• v.53 no.11
• /
• pp.3528-3542
• /
• 2021

A detailed geometrical model of a Krško reactor core was developed using a Monte Carlo neutron transport code MCNP. The main goal of developing an MCNP core model is for it to be used in future research focused on ex-core calculations. A script called McCord was developed to generate MCNP input for an arbitrary fuel cycle configuration from the diffusion based core design package CORD-2, taking advantage of already available material and temperature data obtained in the nuclear core design process. The core model was used to calculate 3D power density profile inside the core. The applicability of the calculated power density distributions was tested by comparison to the CORD-2 calculations, which is regularly used for the nuclear core design calculation verification of the Krško core. For the hot zero power and hot full power states differences between MCNP and CORD-2 in the radial power density profile were <3%. When studying axial power density profiles the differences in axial offset were less than 2.3% for hot full power condition. To further confirm the applicability of the developed model, the measurements with in-core neutron detectors were compared to the calculations, where differences of 5% were observed.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.11a
• /
• pp.301-306
• /
• 1996

To predict the average liquid level under the condition of the countercurrent stratified two-phase flow in a pipe, an analytical model has been suggested. This is made by introducing the interfacial level gradient into the liquid-phase and the gas-phase momentum equations. The analytical method for the gas-phase pressure drop calculation with f$_i$ $\neq$ f$_G$ has also been described using the liquid level prediction model developed in the present study.

• Journal of the Korean Society of Systems Engineering
• /
• v.14 no.2
• /
• pp.24-32
• /
• 2018

System Complexity one of the most common cause failure of the projects, it leads to a lack of understanding about the functions of the system. Hence, the model is developed for communication and furthermore modeling help analysis, design, and understanding of the system. On the other hand, the text-based specification is useful and easy to develop but is difficult to visualize the physical composition, structure, and behaviour or data exchange of the system. Therefore, it is necessary to transform system description into a diagram which clearly depicts the behaviour of the system as well as the interaction between components. According to the International Atomic Energy Agency (IAEA) Safety Glossary, The safety system is a system important to safety, provided to ensure the safe shutdown of the reactor or the residual heat removal from the reactor core, or to limit the consequences of anticipated operational occurrences and design basis accidents. Core Protection Calculator System (CPCS) in Advanced Power Reactor 1400 (APR 1400) Nuclear Power Plant is a safety critical system. CPCS was developed using systems engineering method focusing on Departure from Nuclear Boiling Ratio (DNBR) calculation. Due to the complexity of the system, many diagrams are needed to minimize the risk of ambiguities and lack of understanding. Using Model-Based Systems Engineering (MBSE) software for modeling the DNBR algorithm were used. These diagrams then serve as the baseline of the reverse engineering process and speeding up the development process. In addition, the use of MBSE ensures that any additional information obtained from auxiliary sources can then be input into the system model, ensuring data consistency.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05b
• /
• pp.637-642
• /
• 1996

This paper presents the analysis results for the core degradation processes and the fission product release of the PHEBUS FPT0 experiment using MELCOR1.8.3. The objective of this study is to assess models associated with the core damage and fission product behavior in MELCOR. The calculation results were much improved through sensitivity studies. Thermal/hydraulic behavior in the core and the circuit was well predicted under the intact core geometry. In non-eutectic model case. the UO$_2$ dissolution model in the MELCOR always showed such a tendency that the resulting dissolved UO$_2$ mass was small at the highly oxidized condition due to the model logic. Total H$_2$ generation mass was underpredicted because the stiffner was not modeled and the liner in the shroud was not allowed to be oxidized in MELCOR. Some difficulties were found in modeling the activation product were solved by manipulating the RN input associated with the initial fission product inventory. These problem were occurred because there are no control rod model in MELCOR. Generally the fission product release ratio showed a similar trend compared with the measured data except the activation product. which have no model to simulate in MELCOR.

• Journal of the Korean Society of Systems Engineering
• /
• v.17 no.2
• /
• pp.37-48
• /
• 2021

Improving the plant protection system against unforeseen changes/transients during operation is essential to maintain plant safety. Under this condition, it requires rapid and accurate signal processing. The use of an Intellectual Property (IP) core for floating point calculations for Safety Critical MMIS can make numerical computations easier and more precise, improving system accuracy. It can represent and manipulate rational numbers as well as a much broader range of values with dynamic range in nuclear power plant. Systems engineering approach (SE) is used through the development process, it helps to reduce complexity and avoid omissions and invalid assumptions as delivers a better understanding of the stakeholders needs. For the implementation on the FPGA target board, the 32-bit floating-point arithmetic with IEEE-754 standards has designed using Simulink model in Matlab for all operations of addition, subtraction, multiplication and division and VHDL code generated.

• Nuclear Engineering and Technology
• /
• v.51 no.8
• /
• pp.1871-1885
• /
• 2019

The deterministic MOC code STREAM of the Computational Reactor Physics and Experiment (CORE) laboratory of Ulsan National Institute of Science and Technology (UNIST), was initially designed for the calculation of pressurized water reactor two- and three-dimensional assemblies and cores. Since fast reactors play an important role in the generation-IV concept, it was decided that the code should be upgraded for the analysis of fast neutron spectrum reactors. This paper presents a coupled code - TULIP/STREAM, developed for the fast reactor assembly and core calculations. The TULIP code produces self-shielded multi-group cross-sections using a one-dimensional cylindrical model. The generated cross-section library is used in the STREAM code which solves eigenvalue problems for a two-dimensional assembly and a multi-assembly whole reactor core. Multiplication factors and steady-state power distributions were compared with the reference solutions obtained by the continuous energy Monte-Carlo code MCS. With the developed code, a sensitivity study of the number of energy groups, the order of anisotropic P_N scattering, and the multi-group cross-section generation model was performed on the k_eff and power distribution. The 2D core simulation calculations show that the TULIP/STREAM code gives a k_eff error smaller than 200 pcm and the root mean square errors of the pin-wise power distributions within 2%.