• Nuclear Engineering and Technology
• /
• 제54권12호
• /
• pp.4798-4808
• /
• 2022

The porous corrosion deposits (known as CRUD) adhered to the cladding have an important effect on the heat transfer from fuel rods to coolant in PWRs. The vapor film is the main constituent in the two-phase film boiling model. This paper presents a vapor film thickness correlation, associated with CRUD porosity, CRUD chimney density, CRUD particle size, CRUD thickness and heat flux. The dependences of the vapor film thickness on the various influential factors can be intuitively reflected from this vapor film thickness correlation. The temperature, pressure, and boric acid concentration distributions in CRUD can be well predicted using the two-phase film boiling model coupled with the vapor film thickness correlation. It suggests that the vapor thickness correlation can estimate the vapor film thickness more conveniently than the previously reported vapor thickness calculation methods.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• 제20권2호
• /
• pp.171-183
• /
• 2022

The parent and daughter nuclides in a radioactive decay chain arrive at secular equilibrium once they have a large half-life difference. The characteristics of this equilibrium state can be used to estimate the production time of nuclear materials. In this study, a mathematical model and algorithm that can be applied to radio-chronometry using the radioactive equilibrium relationship were investigated, reviewed, and implemented. A Bateman equation that can analyze the decay of radioactive materials over time was used for the mathematical model. To obtain a differential-based solution of the Bateman equation, an algebraic numerical solution approach and two different matrix exponential functions (Moral and Levy) were implemented. The obtained result was compared with those of commonly used algorithms, such as the Chebyshev rational approximation method and WISE Uranium. The experimental analysis confirmed the similarity of the results. However, the Moral method led to an increasing calculation uncertainty once there was a branching decay, so this aspect must be improved. The time period corresponding to the production of nuclear materials or nuclear activity can be estimated using the proposed algorithm when uranium or its daughter nuclides are included in the target materials for nuclear forensics.

• Nuclear Engineering and Technology
• /
• 제53권2호
• /
• pp.393-398
• /
• 2021

Current risk assessment of Spent Nuclear Fuel (SNF) transportation has the problem of the incomplete risk factors consideration and the general particle diffusion model utilization. In this paper, the accident frequency calculation and the detailed simulation of the accident consequences are coupled by the integrated risk assessment method. The "man-machine-environment" three-dimensional comprehensive risk indicator system is established and quantified to characterize the frequency of the transportation accidents. Consideration of vegetation, building and turbulence effect, the standard k-ε model is updated to simulate radioactive consequence of leakage accidents under complex terrain. The developed method is applied to assess the risk of the leakage accident in the scene of the typical domestic SNF Road Transportation (SNFRT). The critical risk factors and their impacts on the dispersion of the radionuclide are obtained.

• Nuclear Engineering and Technology
• /
• 제34권4호
• /
• pp.286-300
• /
• 2002

An algorithm for the sodium boiling model has been developed for calculation of the void reactivity feedback as well as the fuel and cladding temperatures in the KALIMER core after onset of sodium boiling. Modeling of sodium boiling in liquid metal reactors using sodium as a coolant is necessary because of phenomenon difference comparing with that observed generally in light water reactor systems. The applied model to the algorithm is the multiple-bubble slug ejection model. It allows a finite number of bubbles in a channel at any time. Voiding is assumed to result from formation of bubbies that (ill the whole cross section of the coolant channel except for the liquid film left on the cladding surface. The vapor pressure, currently, is assumed to be uniform within a bubble The present study is focused on not only demonstration of the vapor bubble behavior predicted by the developed model, but also confirmation of a qualitative acceptance for the model. As a result, the model can represent important phenomena in the sodium boiling, but it is found that further effort is also needed for its completition.

• Nuclear Engineering and Technology
• /
• 제54권3호
• /
• pp.1109-1114
• /
• 2022

The SMART (System-integrated Modular Advanced ReacTor) is an integral-type small modular reactor developed by KAERI (Korea Atomic Energy Research Institute). This paper discusses the feasibility and applicability of a 3D-based equivalent model using dynamic condensation method for seismic analysis of a SMART control rod drive mechanism. The equivalent model is utilized for complicated seismic analysis during the design of the SMART. While the 1D-based beam-mass equivalent model is widely used in the nuclear industry for its calculation efficiency, the 3D-based equivalent model is suggested for the seismic analysis of SMART to enhance the analysis accuracy of the 1D-based equivalent model while maintaining its analysis efficiency. To verify the suggested model, acceleration response spectra from seismic analysis based on the 3D-based equivalent model are compared to those from the 1D-based beam-mass equivalent model and experiments. The accuracy and efficiency of the dynamic condensation method are investigated by comparison to analysis results based on the conventional modeling methodology used for seismic analysis.

• Journal of Radiation Protection and Research
• /
• 제45권1호
• /
• pp.45-52
• /
• 2020

Background: Recently, the International Commission on Radiological Protection (ICRP) lowered the dose limit for the eye lens from 150 mSv to 20 mSv, highlighting the importance of accurate lens dose estimation. The ICRP reference computational phantoms used for lens dose calculation are mostly based on the data of Caucasian population, and thus might be inappropriate for Korean population. Materials and Methods: In the present study, a detailed Korean eye model was constructed by determining nine ocular dimensions using the data of Korean subjects. The developed eye model was then incorporated into the adult male and female mesh-type reference Korean phantoms (MRKPs), which were then used to calculate lens doses for photons and electrons in idealized irradiation geometries. The calculated lens doses were finally compared with those calculated with the ICRP mesh-type reference computational phantoms (MRCPs) to observe the effect of ethnic difference on lens dose. Results and Discussion: The lens doses calculated with the MRKPs and the MRCPs were not much different for photons for the entire energy range considered in the present study. For electrons, the differences were generally small, but exceptionally large differences were found at a specific energy range (0.5-1 MeV), the maximum differences being about 10 times at 0.6 MeV in the anteroposterior geometry; the differences are mainly due to the difference in the depth of the lens between the MRCPs and the MRKPs. Conclusion: The MRCPs are generally considered acceptable for lens dose calculations for Korean population, except for the electrons at the energy range of 0.5-1 MeV for which it is suggested to use the MRKPs incorporating the Korean eye model developed in the present study.

• Nuclear Engineering and Technology
• /
• 제54권3호
• /
• pp.790-802
• /
• 2022

Large-scale reactor simulation often requires the use of Monte Carlo calculation techniques to estimate important reactor parameters. One drawback of these Monte Carlo calculation techniques is they inevitably result in some uncertainty in calculated quantities. The present study includes parametric uncertainty quantification (UQ) and sensitivity analysis (SA) on the Advanced Test Reactor Critical (ATRC) facility housed at Idaho National Laboratory (INL) and addresses some complications due to Monte Carlo uncertainty when performing these analyses. This approach for UQ/SA includes consideration of Monte Carlo code uncertainty in computed sensitivities, consideration of uncertainty from directly measured parameters and a comparison of results obtained from brute-force Monte Carlo UQ versus UQ obtained from a surrogate model. These methodologies are applied to the uncertainty and sensitivity of k_eff for two sets of uncertain parameters involving fuel plate geometry and fuel plate composition. Results indicate that the less computationally-expensive method for uncertainty quantification involving a linear surrogate model provides accurate estimations for k_eff uncertainty and the Monte Carlo uncertainty in calculated k_eff values can have a large effect on computed linear model parameters for parameters with low influence on k_eff.

• Journal of the Korean Society of Safety
• /
• 제35권3호
• /
• pp.96-104
• /
• 2020

The numerical simulations were conducted to investigate the thermal-fluid phenomena occurred inside the experimental apparatus during a PCCS, used to remove heat released in accidents from a containment of light water nuclear power plant, operation. Numerical simulations of the flow and heat transfer caused by wall condensation inside the containment simulation vessel (CSV), which equipped with 18 vertical heat exchanger tubes, were conducted using the commercial computational fluid dynamics (CFD) software ANSYS-CFX. Shear stress transport (SST) and the wall condensation model were used for turbulence closure and wall condensation, respectively. The simulation using the actual size of the apparatus. However, rather than simulating the whole experimental apparatus in consideration of the experimental cases, calculation resources, and calculation time, the simulation model was prepared only in CSV. Selective simulation was conducted to verify the effects of non-condensable gas(NC gas) concentration, CSV internal pressure, and wall sub-cooling conditions. First, as a result of the internal flow of CSV, it was observed that downward flow due to condensation occurred surface of the vertical tube and upward flow occurred in the distant place. Natural convection occurred actively around the heat exchanger tube. Due to this rising and falling internal flow, natural circulation occurred actively around the heat exchanger tubes. Next, in order to check the performance of built-in condensation model using according to the non-condensable gas concentration, CSV internal flow and wall sub-cooling, the heat flux values were compared with the experimental results. On average, the results were underestimated with and error of about 25%. In addition, the influence of CSV internal pressure and wall sub-cooling was small, but when the condensate was highly generated due to the low non-condensable gas concentration, the error was large compared to the experimental values. This is considered to be due to the nature of the condensation model of the CFX code. However, in spite of the limitations of CFD, it is valid to use the built-in condensation model of CFD for PCCS performance prediction from a conservative perspective.

• Nuclear Engineering and Technology
• /
• 제34권2호
• /
• pp.154-164
• /
• 2002

Neutron cross section data on Dy-160, Dy-161, Dy-162, Dy-163 and Dy-164 were calculated and evaluated in the energy range of 1 keV to 20 MeV using a spherical optical model, statistical model and pre-equilibrium model. The energy dependent optical model potential parameters were obtained based on the recent experimental data. The width fluctuation correction in Hauser-Feshbach particle decay and the quantum mechanical approach in pre-equilibrium analysis were introduced and gave a better cross section calculation in EMPIRE-II. The total, elastic scattering and threshold reaction cross sections were evaluated and compared with the evaluated files. The model calculated (n, tot), (n, ${\gamma}$) and (n, p) cross sections were in good agreement with the experimental data in the measured energy range. The results will be applied to ENDF/B-VI for data improvement.

• Journal of the Korean Society of Safety
• /
• 제37권5호
• /
• pp.80-88
• /
• 2022

The Korean nuclear industry had developed the SPACE (Safety and Performance Analysis Code for nuclear power plants) code, which adopts a two-fluid, three-field model that is comprised of gas, continuous liquid and droplet fields and has the capability to simulate three-dimensional models. According to the revised law by the Nuclear Safety and Security Commission (NSSC) in Korea, the multiple failure accidents that must be considered for the accident management plan of a nuclear power plant was determined based on the lessons learned from the Fukushima accident. Generally, to improve the reliability of the calculation results of a safety analysis code, verification is required for the separate and integral effect experiments. Therefore, the goal of this work is to verify the calculation capability of the SPACE code for multiple failure accidents. For this purpose, an experiment was conducted to simulate a Control Element Drive Mechanism (CEDM) break with a safety injection failure using the ATLAS test facility, which is operated by Korea Atomic Energy Research Institute (KAERI). This experiment focused on the comparison between the experiment results and code calculation results to verify the performance of the SPACE code. The results of the overall system transient response using the SPACE code showed similar trends with the experimental results for parameters such as the system pressure, mass flow rate, and collapsed water level in component. In conclusion, it can be concluded that the SPACE code has sufficient capability to simulate a CEDM break with a safety injection failure accident.