• Nuclear Engineering and Technology
• /
• v.42 no.5
• /
• pp.500-519
• /
• 2010

AEGIS is a lattice physics code incorporating the latest advances in lattice physics computation, innovative calculation models and efficient numerical algorithms and is mainly used for light water reactor analyses. Though the primary objective of the AEGIS code is the preparation of a cross section set for SCOPE2 that is a three-dimensional pin-by-pin core analysis code, the AEGIS code can handle not only a fuel assembly but also multi-assemblies and a whole core geometry in two-dimensional geometry. The present paper summarizes the major calculation models and part of the verification/validation efforts related to the AEGIS code.

• Nuclear Engineering and Technology
• /
• v.53 no.3
• /
• pp.741-751
• /
• 2021

There has been a deep interest in trying to find better-performing fuel clad motivated by the desire to decrease the likelihood of the reactor barrier failure like what happened in Fukushima in recent years. In this study, the effect of move towards accident tolerant fuel (ATF) cladding as the most attracting concept for improving reactor safety is investigated for SMART modular reactor. These reactors have less production cost, short construction time, better safety and higher power density. The SiC and FeCrAl materials are considered as the most potential candidate for ATF cladding, and the results are compared with Zircaloy cladding material from reactor physics point of view. In this paper, the calculations are performed by generating PMAX library by DRAGON lattice physics code to be used for further reactor core analysis by PARCS code. The differential and integral worth of control and safety rods, reactivity coefficient, power and temperature distributions, and boric acid concentration during the cycle are analyzed and compared from the conventional fuel cladding. The rod ejection accident (REA) is also performed to study how the power changed in response to presence of the ATF cladding in the reactor core. The key quantitative finding can be summarized as: 20 ℃ (3%) decrease in average fuel temperature, 33 pcm (3%) increase in integral rod worth and cycle length, 1.26 pcm/℃ (50%) and 1.05 pcm/℃ (16%) increase in reactivity coefficient of fuel and moderator, respectively.

• Nuclear Engineering and Technology
• /
• v.54 no.1
• /
• pp.30-35
• /
• 2022

In a nuclear reactor, gamma radiation is the dominant energy deposition in non-fuel regions. Heat is generated upon gamma deposition and consequently affects the mechanical and thermal structure of the material. Therefore, the safety of samples should be carefully considered so that their integrity and quality can be retained. To evaluate relevant parameters, an in-core gamma thermometer (GT) was used to measure gamma heating (GH) throughout the operation of the McMaster nuclear reactor (MNR) at four irradiation sites. Additionally, a Monte Carlo reactor physics code (Serpent-2) was utilized to model the MNR with the GT located in the same irradiation sites used in the measurement to verify its predictions against measured GH. This research aids in the development of modeling, calculation, and prediction of the GH utilizing Serpent-2 as well as implementing a new GH measurement at the MNR core. After all uncertainties were quantified for both approaches, comparable GH profiles were observed between the measurements and calculations. In addition, the GH values found in the four sites represent a strong level of radiation based on the distance of the sample from the core. In this study, the maximum and minimum GH values were found at 0.32 ± 0.05 W/g and 0.15 ± 0.02 W/g, respectively, corresponding to 320 Sv/s and 150 Sv/s. These values are crucial to be considered whenever sample is planned to be irradiated inside the MNR core.

• Nuclear Engineering and Technology
• /
• v.53 no.6
• /
• pp.1736-1746
• /
• 2021

The results of full-scale numerical experiments of a hybrid thorium-containing fuel cell facility operating in a close-to-critical state due to a controlled source of fusion neutrons are discussed in this work. The facility under study was a complex consisting of two blocks. The first block was based on the concept of a high-temperature gas-cooled thorium reactor core. The second block was an axially symmetrical extended plasma generator of additional neutrons that was placed in the near-axial zone of the facility blanket. The calculated models of the blanket and the plasma generator of D-T neutrons created within the work allowed for research of the neutronic parameters of the facility in stationary and pulse-periodic operation modes. This research will make it possible to construct a safe facility and investigate the properties of thorium fuel, which can be continuously used in the epithermal spectrum of the considered hybrid fusion-fission reactor.

• Nuclear Engineering and Technology
• /
• v.54 no.4
• /
• pp.1253-1260
• /
• 2022

This paper reports the control method for the core power of the China initiative Accelerator Driven System (CiADS) facility. In the CiADS facility, an intense external neutron source provided by a proton accelerator coupled to a spallation target is used to drive a sub-critical reactor. Without any control rod inside the sub-critical reactor, the core power is controlled by adjusting the proton beam intensity. In order to continuously change the beam intensity, an adjustable aperture is considered to be used at the Low Energy Beam Transport (LEBT) line of the accelerator. The aperture size is adjusted based on the Proportional Integral Derivative (PID) controllers, by comparing either the setting beam intensity or the setting core power with the measured value. To evaluate the proposed control method, a CiADS core model is built based on the point reactor kinetics model with six delayed neutron groups. The simulations based on the CiADS core model have indicated that the core power can be controlled stably by adjusting the aperture size. The response time in the adjustment of the core power depends mainly on the adjustment time of the beam intensity.

• Nuclear Engineering and Technology
• /
• v.51 no.4
• /
• pp.968-976
• /
• 2019

This paper presents a number of verification case studies for a recently developed sensitivity/uncertainty code package. The code package, ROMUSE (Reduced Order Modeling based Uncertainty/Sensitivity Estimator) is an effort to provide an analysis tool to be used in conjunction with reactor core simulators, in particular the Virtual Environment for Reactor Applications (VERA) core simulator. ROMUSE has been written in C++ and is currently capable of performing various types of parameter perturbations and associated sensitivity analysis, uncertainty quantification, surrogate model construction and subspace analysis. The current version 2.0 has the capability to interface with the Design Analysis Kit for Optimization and Terascale Applications (DAKOTA) code, which gives ROMUSE access to the various algorithms implemented within DAKOTA, most importantly model calibration. The verification study is performed via two basic problems and two reactor physics models. The first problem is used to verify the ROMUSE single physics gradient-based range finding algorithm capability using an abstract quadratic model. The second problem is the Brusselator problem, which is a coupled problem representative of multi-physics problems. This problem is used to test the capability of constructing surrogates via ROMUSE-DAKOTA. Finally, light water reactor pin cell and sodium-cooled fast reactor fuel assembly problems are simulated via SCALE 6.1 to test ROMUSE capability for uncertainty quantification and sensitivity analysis purposes.

• Nuclear Engineering and Technology
• /
• v.48 no.3
• /
• pp.635-641
• /
• 2016

Because of prohibitive data storage requirements in large-scale simulations, the memory problem is an obstacle for Monte Carlo (MC) codes in accomplishing pin-wise three-dimensional (3D) full-core calculations, particularly for whole-core depletion analyses. Various kinds of data are evaluated and quantificational total memory requirements are analyzed based on the Reactor Monte Carlo (RMC) code, showing that tally data, material data, and isotope densities in depletion are three major parts of memory storage. The domain decomposition method is investigated as a means of saving memory, by dividing spatial geometry into domains that are simulated separately by parallel processors. For the validity of particle tracking during transport simulations, particles need to be communicated between domains. In consideration of efficiency, an asynchronous particle communication algorithm is designed and implemented. Furthermore, we couple the domain decomposition method with MC burnup process, under a strategy of utilizing consistent domain partition in both transport and depletion modules. A numerical test of 3D full-core burnup calculations is carried out, indicating that the RMC code, with the domain decomposition method, is capable of pin-wise full-core burnup calculations with millions of depletion regions.

• Nuclear Engineering and Technology
• /
• v.53 no.10
• /
• pp.3133-3150
• /
• 2021

This work deals with the validation of a high-fidelity multiphysics system coupling the Serpent 2 Monte Carlo neutron transport code with SUBCHANFLOW, a subchannel thermalhydraulics code, and TRANSURANUS, a fuel-performance analysis code. The results for a full-core pin-by-pin burnup calculation for the ninth operating cycle of the Temelín II VVER-1000 plant, which starts from a fresh core, are presented and assessed using experimental data. A good agreement is found comparing the critical boron concentration and a set of pin-level neutron flux profiles against measurements. In addition, the calculated axial and radial power distributions match closely the values reported by the core monitoring system. To demonstrate the modeling capabilities of the three-code coupling, pin-level neutronic, thermalhydraulic and thermomechanic results are shown as well. These studies are encompassed in the final phase of the EU Horizon 2020 McSAFE project, during which the Serpent-SUBCHANFLOW-TRANSURANUS system was developed.

• Nuclear Engineering and Technology
• /
• v.51 no.2
• /
• pp.356-368
• /
• 2019

This paper presents the verification and validation (V&V) of the STREAM/RAST-K 2.0 code system for a pressurized water reactor (PWR) analysis. A lattice physics code STREAM and a nodal diffusion code RAST-K 2.0 have been developed by a computational reactor physics and experiment laboratory (CORE) of Ulsan National Institute of Science and Technology (UNIST) for an accurate two-step PWR analysis. The calculation modules of each code were already verified against various benchmark problems, whereas this paper focuses on the V&V of linked code system. Three PWR type reactor cores, OPR-1000, three-loop Westinghouse reactor core, and APR-1400, are selected as V&V target plants. This code system, for verification, is compared against the conventional code systems used for the calculations in nuclear design reports (NDRs) and validated against measured plant data. Compared parameters are as follows: critical boron concentration (CBC), axial shape index (ASI), assembly-wise power distribution, burnup distribution and peaking factors. STREAM/RAST-K 2.0 shows the RMS error of critical boron concentration within 20 ppm, and the RMS error of assembly power within 1.34% for all the cycles of all reactors.

• Nuclear Engineering and Technology
• /
• v.53 no.11
• /
• pp.3517-3527
• /
• 2021

New 80-MW (electric) ultra-long-life sodium cooled fast reactor core having inherent safety characteristics is designed with heterogeneous fuel assemblies comprised of driver and blanket fuel rods. Several options using upper sodium plenum and SSFZ (Special Sodium Flowing Zone) for reducing sodium void reactivity are neutronically analyzed in this core concept in order to improve the inherent safety of the core. The SSFZ allowing the coolant flow from the peripheral fuel assemblies increases the neutron leakage under coolant expansion or voiding. The Monte Carlo calculations were used to design the cores and analyze their physics characteristics with heterogeneous models. The results of the design and analyses show that the final core design option has a small burnup reactivity swing of 618 pcm over ~54 EFPYs cycle length and a very small sodium void worth of ~35pcm at EOC (End of Cycle), which leads to the satisfaction of all the conditions for inherent safety with large margin based on the quasi-static reactivity balance analysis under ATWS (Anticipated Transient Without Scram).