• Nuclear Engineering and Technology
• /
• v.54 no.8
• /
• pp.3073-3084
• /
• 2022

The propagation of radiation source uncertainties in spent nuclear fuel (SNF) cask shielding calculations is presented in this paper. The uncertainty propagation employs the depletion and source term outputs of the deterministic code STREAM as input to the transport simulation of the Monte Carlo (MC) codes MCS and MCNP6. The uncertainties of dose rate coming from two sources: nuclear data and modeling parameters, are quantified. The nuclear data uncertainties are obtained from the stochastic sampling of the cross-section covariance and perturbed fission product yields. Uncertainties induced by perturbed modeling parameters consider the design parameters and operating conditions. Uncertainties coming from the two sources result in perturbed depleted nuclide inventories and radiation source terms which are then propagated to the dose rate on the cask surface. The uncertainty analysis results show that the neutron and secondary photon dose have uncertainties which are dominated by the cross section and modeling parameters, while the fission yields have relatively insignificant effect. Besides, the primary photon dose is mostly influenced by the fission yield and modeling parameters, while the cross-section data have a relatively negligible effect. Moreover, the neutron, secondary photon, and primary photon dose can have uncertainties up to about 13%, 14%, and 6%, respectively.

• Nuclear Engineering and Technology
• /
• v.53 no.1
• /
• pp.44-60
• /
• 2021

In this paper, we validate the decay heat calculation capability via a two-step method to analyze spent nuclear fuel (SNF) discharged from pressurized water reactors (PWRs). The calculation method is implemented with a lattice code STREAM and a nodal diffusion code RAST-K. One of the features of this method is the direct consideration of three-dimensional (3D) core simulation conditions with the advantage of a short simulation time. Other features include the prediction of the isotope inventory by Lagrange non-linear interpolation and the use of power history correction factors. The validation is performed with 58 decay heat measurements of 48 fuel assemblies (FAs) discharged from five PWRs operated in Sweden and the United States. These realistic benchmarks cover the discharge burnup range up to 51 GWd/MTU, 23.2 years of cooling time, and spanning an initial uranium enrichment range of 2.100-4.005 wt percent. The SNF analysis capability of STREAM is also employed in the code-to-code comparison. Compared to the measurements, the validation results of the FA calculation with RAST-K are within ±4%, and the pin-wise results are within ±4.3%. This paper successfully demonstrates that the developed decay heat calculation method can perform SNF back-end cycle analyses.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1995.10a
• /
• pp.41-47
• /
• 1995

A parametric study bas been performed for the various refueling schemes of CANDU 6 reactor loaded with reference DUPIC fuel. The optimum discharge burnup was determined such that the peak bundle power is minimized for the equilibrium core. Based on the results of instantaneous core calculation using patterned random age distributions, it was decided to perform the refueling simulations only for 2-bundle and 4-bundle shift refueling schemes. The 600 FPD simulation has shown that the operational margins of the channel and bundle power to the license limits are 7.9% and 17.1%, respectively, for 2-bundle shift refueling scheme. The 4-bundle shift refueling scheme also satisfies the license limits and the operational margins of the channel and bundle power are 7.1% and 9.8%, respectively. The result of refueling simulation indicate the possibility of using reference DUPIC fuel in current CANDU 6 reactor.

• Proceedings of the Korean Nuclear Society Conference
• /
• 2004.10a
• /
• pp.1047-1048
• /
• 2004

Simulation Coated layers of a nuclear fuel particle were evaluated by field emission scanning electron microscopy and nano-indentation method to give basic data to estimate 'Amoeba effect' and give an optimum fabrication condition and high quality control. Coated layers on the fuel kernel are in the order of buffer pyrolytic carbon, inner pyrolytic carbon, silicon carbide and outer pyrolytic carbon layers, which average thicknesses are 95, 25, 30 and 28 ${\mu}m$, respectively. Their densities and hardnesses are 1.08, 1.15, 3.18, 1.82 $g/cm^3$ and 0.522, 0.874, 9.641, and 2.726 GPa, respectively. Comparing theoretical density of pyrolytic carbon of 2.22 $g/cm^3$, the relative amount of porosity in each layer is about 52% for buffer, 48% for inner PyC and 18% for outer PyC.

• Nuclear Engineering and Technology
• /
• v.53 no.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
• /
• v.50 no.2
• /
• pp.297-305
• /
• 2018

Fuel assembly (FA) bow in pressurized water reactor (PWR) cores is considered to be a complex process with a large number of influencing mechanisms and several unknowns. Uncertainty and sensitivity analyses are a common way to assess the predictability of such complex phenomena. To perform such analyses, a structural model of a row of 15 FAs in the reactor core is implemented with the finite-element code ANSYS Mechanical APDL. The distribution of lateral hydraulic forces within the core row is estimated based on a two-dimensional Computational Fluid Dynamics model with porous media, assuming symmetric or asymmetric core inlet and outlet flow profiles. The influence of the creep rate on the bow amplitude is tested based on different creep models for guide tubes and fuel rods. Different FA initial states are considered: fresh FAs or FAs with higher burnup, which may be initially straight or exhibit an initial bow from previous cycles. The simulation results over one reactor cycle demonstrate that changes in the creep rate and the hydraulic conditions may have a considerable impact on the bow amplitudes and the bow patterns. A good knowledge of the specific creep behavior and the hydraulic conditions is therefore crucial for making reliable predictions.

• Nuclear Engineering and Technology
• /
• v.51 no.2
• /
• pp.410-423
• /
• 2019

Fire suppression using a water mist nozzle directly above an n-Heptane pool in a fire compartment with a door opening was numerically investigated using the Fire Dynamics Simulator (FDS) for the purpose of application in nuclear power plants. Input parameters for the numerical simulation were determined by experimental measurements. Water mist was activated 10 s after the fire began. The sensitivity analysis was conducted for three input parameters: total number of cubic cells of 6032-2,926,400, droplets per second of 1000-500,000, and extinguishing coefficient of 0-100. In a new simple calibration method of this study, the extinguishing coefficient yielding the fire suppression time closest to that measured by experiments was found for use as the FDS simulation input value. When the water mist jet flow made contact with the developed fire, the heat release rate instantaneously increased, and then rapidly decreased. This phenomenon occurred with a displacement of the flame near the liquid fuel pool. Changing the configuration of the door opening with different aspect ratios and opening ratios had impact on the maximum value of the heat release rate due to the flame displacement.

• Nuclear Engineering and Technology
• /
• v.54 no.8
• /
• pp.3154-3165
• /
• 2022

During a loss-of-coolant accident (LOCA) in the pressurized water reactor (PWR), there is a possibility that high temperature and internal pressure of the fuel rods lead to ballooning of the cladding, which causes a partial blockage of flow area in a subchannel. Such flow blockage would influence the core coolant flow, thus affecting the core heat transfer during a reflooding phase and subsequent severe accident. However, most of the system analysis codes simulate the accident process based on the assumed channel blockage ratio, resulting in the fact that the simulation results are not consistent with the actual situation. This paper integrates the developed core Fuel Rod Thermal-Mechanical Behavior analysis (FRTMB) module into the self-developed severe accident analysis code ISAA. At the same time, the existing flow blockage model is improved to make it possible to simulate the change of flow distribution due to fuel rod deformation. Finally, the ISAA-FRTMB is used to simulate the QUENCH-LOCA-0 experiment to verify the correctness and effectiveness of the improved flow blockage model, and then the effect of clad ballooning on core heat transfer and subsequent parts of core degradation is analyzed.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1575-1588
• /
• 2019

Fracture near the U-10Mo/cladding material interface impacts fuel service life. In this work, a mesoscale stress model is developed with the fuel foil considered as a porous medium having gas bubbles and bearing bubble pressure and surface tension. The models for the evolution of bubble volume fraction, size and internal pressure are also obtained. For a U-10Mo/Al monolithic fuel plate under location-dependent irradiation, the finite element simulation of the thermo-mechanical coupling behavior is implemented to obtain the bubble distribution and evolution behavior together with their effects on the mesoscale stresses. The numerical simulation results indicate that higher macroscale tensile stresses appear close to the locations with the maximum increments of fuel foil thickness, which is intensively related to irradiation creep deformations. The maximum mesoscale tensile stress is more than 2 times of the macroscale one on the irradiation time of 98 days, which results from the contributions of considerable volume fraction and internal pressure of bubbles. This study lays a foundation for the fracture mechanism analysis and development of a fracture criterion for U-10Mo monolithic fuels.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3450-3463
• /
• 2023

In the two-step analysis of Pebble-Bed type High-Temperature Gas-Cooled Reactor (PB-HTGR), the lattice physics calculation for the generation of homogenized cross-sections is based on the fuel pebble. However, the randomly-dispersed fuel particles in the fuel pebble introduce double heterogeneity and randomness. Compared to the deterministic method, the Monte Carlo method which is flexible in geometry modeling provides a high-fidelity treatment. Therefore, the Monte Carlo code NECP-MCX is extended in this study to perform the lattice physics calculation of the PB-HTGR. Firstly, the capability for the simulation of randomly-dispersed media, using the explicit modeling approach, is developed in NECP-MCX. Secondly, the capability for the generation of the homogenized cross-section is also developed in NECP-MCX. Finally, simplified PB-HTGR problems are calculated by a two-step neutronics analysis tool based on Monte Carlo homogenization. For the pebble beds mixed by fuel pebble and graphite pebble, the bias is less than 100 pcm when compared to the high-fidelity model, and the bias is increased to 269 pcm for pebble bed mixed by depleted fuel pebble. Numerical results show that the Monte Carlo lattice physics calculation for the two-step analysis of PB-HTGR is feasible.