• Nuclear Engineering and Technology
• /
• v.44 no.2
• /
• pp.151-160
• /
• 2012

The advantages for using Monte Carlo methods to analyze full-core reactor configurations include essentially exact representation of geometry and physical phenomena that are important for reactor analysis. But this substantial advantage comes at a substantial cost because of the computational burden, both in terms of memory demand and computational time. This paper focuses on the challenges facing full-core Monte Carlo for keff calculations and the prospects for Monte Carlo becoming a routine tool for reactor analysis.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.8-14
• /
• 2016

A numerical analysis is conducted to estimate the core expansion and the energy behaviors induced by a core disruptive accident in a sodium-cooled fast reactor. The numerical formulation based on underwater explosion theory is carried out to simulate the core explosion inside the reactor vessel. The transient pressure, temperature and expansion of the core are examined by solving the equation of state and nonlinear governing equation of momentum conservation in one-dimensional spherical coordinates. The energy balance inside the computation domain is examined during the core expansion process. Heat transfer between the core and the sodium coolant, and the bubble rise during the expansion process are briefly investigated.

• Nuclear Engineering and Technology
• /
• v.34 no.2
• /
• pp.117-131
• /
• 2002

Core disruptive accidents have been investigated at Korea Atomic Energy Research Institute(KAERI) as part of the work to demonstrate the inherent and ultimate safety of conceptual design of the Korea Advanced Liquid Metal Reactor(KALIMER), a 150 MWe pool- type sodium cooled prototype fast reactor that uses U-Pu-Zr metallic fuel. In this study, a simple method and associated computer program, SCHAMBETA, was developed using a modified Bethe-Tait method to simulate the kinetics and thermodynamic behavior of a homogeneous spherical core over the period of the super-prompt critical power excursion induced by the ramp reactivity insertion. Calculations of the energy release during excursions in the sodium-voided core of the KALIMER were subsequently performed using the SCHAMBETA code for various reactivity insertion rates up to 100 S/s, which has been widely considered to be the upper limit of ramp rates due to fuel compaction. Benchmark calculations were made to compare with the results of more detailed analysis for core meltdown energetics of the oxide fuelled fast reactor. A set of parametric studies were also performed to investigate the sensitivity of the results on the various thermodynamics and reactor parameters.

• Nuclear Engineering and Technology
• /
• v.53 no.6
• /
• pp.1731-1735
• /
• 2021

The operation of the nuclear reactor requires accurate and fast methods and techniques for analysing its kinetics. These techniques become even more important when the MOX-fuel is used due to the lower value of delayed neutron fraction 𝛽 for ²³⁹Pu. Based on a Birth-and-Death process review, the mathematical model of thermal reactor core has been proposed different from existing ones. The analytical method for thermal point-reactor parameters evaluation is described within this work. The proposed method is applied for analysis of the unsteady transient processes taking place in a thermal reactor at its start-up or shutdown power change, as well as during small accidental power variation from the rated value. Theoretical determination of MASURCA reactor core reactivity through the analysis of experimental data on neutron time spectra was made.

• 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.

• Transactions of the Korean Society of Pressure Vessels and Piping
• /
• v.15 no.2
• /
• pp.58-63
• /
• 2019

This paper presents preliminary analysis and results on IASCC sensitivity of a core shroud in the reactor pressure vessel. First, neutron irradiation flux distribution of the reactor internals was calculated by using the Monte Carlo simulation code, MCNP6.1 and the nuclear data library, ENDF/B-VII.1. Second, based on the neutron irradiation flux distribution, temperature and stress distributions of the core shroud during normal operation were determined by performing finite element analysis using the commercial finite element analysis program, ABAQUS, considering irradiation aging-related degradation mechanisms. Last, IASCC sensitivity of the core shroud was assessed by using the IASCC sensitivity definition of EPRI MRP-211 and the finite element analysis results. As a result of the preliminary analysis, it was found that the point at which the maximum IASCC sensitivity is derived varies over operating time, initially moving from the shroud plate located in the center of the core to the top shroud plate-ring connection brace over operating time. In addition, it was concluded that IASCC will not occur on the core shroud even after 60 years of operation (40EFPYs) because the maximum IASCC sensitivity is less than 0.5.

• Nuclear Engineering and Technology
• /
• v.50 no.1
• /
• pp.54-67
• /
• 2018

There have been recent efforts to establish methods for high-fidelity and multi-physics simulation with coupled thermal-hydraulic (T/H) and neutronics codes for the entire core of a light water reactor under accident conditions. Considering the computing power necessary for a pin-by-pin analysis of the entire core, subchannel-scale T/H analysis is considered appropriate to achieve acceptable accuracy in an optimal computational time. In the present study, the applicability of in-house code CUPID of the Korea Atomic Energy Research Institute was extended to the subchannel-scale T/H analysis. CUPID is a component-scale T/H analysis code, which uses three-dimensional two-fluid models with various closure models and incorporates a highly parallelized numerical solver. In this study, key models required for a subchannel-scale T/H analysis were implemented in CUPID. Afterward, the code was validated against four subchannel experiments under unheated and heated single-phase incompressible flow conditions. Thereafter, a subchannel-scale T/H analysis of the entire core for an Advanced Power Reactor 1400 reactor core was carried out. For the high-fidelity simulation, detailed geometrical features and individual rod power distributions were considered in this demonstration. In this study, CUPID shows its capability of reproducing key phenomena in a subchannel and dealing with the subchannel-scale whole core T/H analysis.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.214-219
• /
• 2004

The HANARO, a multi-purpose research reactor of 30 MWth, open-tank-in-pool type, has been under normal operation since its initial criticality in February, 1995. The HANARO is composed of inlet plenum, grid plate, core channel with flow tubes and chimney. The reactor core channel is located at about twelve m (12 m) depth of the reactor pool and cold by the upward flow that the coolant enters the lower inlet of the plenum, rises up through the grid plate and the core channel and exit through the outlet of chimney. A guide tube is extended from the reactor core to the top of the reactor chimney for easily un/loading a target under the reactor normal operation. But active coolant through the core can be Quickly raised up to the top of the chimney through the guide tube by jet flow. This paper is described an analytical analysis to study the flow behavior through the guide tube under reactor normal operation and unloading the target. As results, it was conformed through the analysis results that the flow rate, about fourteen kilogram per second (14 kg/s) suppressed the guide tube jet and met the design cooling flow rate in a circular flow tube, and that the fission moly target cooling flow rate met the minimum flow rate to cool the target.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.1-6
• /
• 2005

The HANARO, a multi-purpose research reactor of 30 MWth, open-tank-in-pool type, has been under normal operation since its initial criticality in February, 1995. The HANARO is composed af inlet plenum, grid plate, core channel with flow tubes and chimney. The reactor core channel is located at about twelve meters (12 m) depth of the reactor pool and cooled by the upward flow that the coolant enters the lower inlet of the plenum, rises up through the grid plate and the core channel and comes out from the outlet of chimney. A fission moly guide tube is extended from the reactor core to the top of the reactor chimney for easily loading a fission moly target under the reactor normal operation. But active coolant through the core can be quickly raised up to the top of the chimney through the guide tube by jet flow. This paper describes an analytical analysis that is the study of the flow behavior through the guide tube under reactor normal operation and unloading the target. As results, it was conformed through the analysis results that the flow rate, reduced to about fourteen kilogram per second (14 kg/s) from the original flow rate of sixteen point three kilogram per second (16.3 kg/s) did not show the guide tube jet.

• Nuclear Engineering and Technology
• /
• v.53 no.1
• /
• pp.93-102
• /
• 2021

The analysis of the fluid flow characteristics in reactor pressure vessel is an important part of the hydraulic design of nuclear power plant, which is related to the structure design of reactor internals, the flow distribution at core inlet and the safety of nuclear power plant. The flow distribution and mixing characteristics in the pressurized reactor vessel for the 1000MWe advanced pressurized water reactor is analyzed by using Computational Fluid Dynamics (CFD) method in this study. The geometry model of the full-scaled reactor vessel is built, which includes the cold and hot legs, downcomer, lower plenum, core, upper plenum, top plenum, and is verified with some parameters in DCD. Under normal condition, it is found that the flow skirt, core plate holes and outlet pipe cause pressure loss. The maximum and minimum flow coefficient is 1.028 and 0.961 respectively, and the standard deviation is 0.019. Compared with other reactor type, it shows relatively uniform of the flow distribution at the core inlet. The coolant mixing coefficient is investigated with adding additional variables, showing that mass transfer of coolant occurs near the interface. The coolant mainly distributes in the 90° area of the corresponding core inlet, and mixes at the interface with the coolant from the adjacent cold leg. 0.1% of corresponding coolant is still distributed at the inlet of the outer-ring components, indicating wide range of mixing coefficient distribution.