• Nuclear Engineering and Technology
• /
• v.55 no.6
• /
• pp.2112-2124
• /
• 2023

With the rise of economic and safety standards for nuclear reactors, new concepts of Gen-IV reactors and modular reactors showed more complex designs that challenge current tools for reactor physics analysis. A Monte Carlo (MC) two-step method was proposed in this work. This calculation scheme uses the continuous-energy MC method to generate multi-group cross-sections from heterogeneous models. The multi-group MC method, which can adapt locally-heterogeneous models, is used in the core calculation step. This calculation scheme is verified using a Gen-IV modular lead-based fast reactor (LFR) benchmark case. The influence of homogenized patterns, scatter approximations, flux separable approximation, and local heterogeneity in core calculation on simulation results are investigated. Results showed that the cross-sections generated using the 3D assembly model with a locally heterogeneous representation of control rods lead to an accurate estimation with less than 270 pcm bias in core reactivity, 0.5% bias in control rod worth, and 1.5% bias on power distribution. The study verified the applicability of multi-group cross-sections generated with the MC method for LFR analysis. The study also proved the feasibility of multi-group MC in core calculation with local heterogeneity, which saves 85% time compared to the continuous-energy MC.

• Nuclear Engineering and Technology
• /
• v.49 no.1
• /
• pp.267-276
• /
• 2017

In any reactor physics analysis, the instantaneous power distribution in the core can be calculated when the actual bundle-wise burnup distribution is known. Considering the fact that CANDU (Canada Deuterium Uranium) utilizes on-power refueling to compensate for the reduction of reactivity due to fuel burnup, in the CANDU fuel management analysis, snapshots of power and burnup distributions can be obtained by simulating and tracking the reactor operation over an extended period using various tools such as the $^*SIMULATE$ module of the Reactor Fueling Simulation Program (RFSP) code. However, for some studies, such as an evaluation of a conceptual design of a next-generation CANDU reactor, the preferred approach to obtain a snapshot of the power distribution in the core is based on the patterned-channel-age model implemented in the $^*INSTANTAN$ module of the RFSP code. The objective of this approach is to obtain a representative snapshot of core conditions quickly. At present, such patterns could be generated by using a program called RANDIS, which is implemented within the $^*INSTANTAN$ module. In this work, we present an alternative approach to derive the patterned-channel-age model where a simulated-annealing-based algorithm is used to find such patterns, which produce reasonable power distributions.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1998.05b
• /
• pp.771-776
• /
• 1998

TORT is employed for fast neutron fluence calculation at the reactor pressure vessel. KORI Unit 1 reactor at cycle 1 is modeled for this calculation. Three-dimensional cycle averaged assembly power distributions for KORI Vnit 1 at cycle 1 are calculated by using the core physics code, NESTLE 5.0. The root mean square error is within 4.3% compared with NDR (Nuclear Design Report) far all burnup steps. The C/E (Calculated/Experimental) values for the in-vessel dosimeters distribute between 0.98 and 1.36. The most updated cross-section library. BUGLE-96 based on ENDF/B-VI is used for the neutron fluence calculation. The makimum fast neutron nun calculated on reactor pressure vessel for KORI Unit 1 operated for 411.41 effgctive full power days is 1.784x10$^{18}$ n/$\textrm{cm}^2$. The position of the maximum neutron fluence in RPV wall 1/4 T is nearby 60cm below the midplane at zero degree.

• Nuclear Engineering and Technology
• /
• v.55 no.4
• /
• pp.1439-1447
• /
• 2023

The fuel cycle characteristics of the IVG.1M reactor were studied within the framework of the research reactor conversion program to modernize the IVG.1M reactor. Optimum use of the nuclear fuel and reactor was achieved through routine methods which included partial fuel reloading combined with scheduled maintenance operations. Since, the additional problem in planning the fuel cycle of the IVG.1M reactor was the poisoning of the beryllium parts of the core, reflector, and control system. An assessment of the residual power and composition of spent fuel is necessary for the selection and justification of the technology for its subsequent management. Computational studies were performed using the MCNP6.1 program and the neutronics model of the IVG.1M reactor. The proposed scheme of annual partial fuel reloading allows for maintaining a high reactor reactivity margin, stabilizing it within 2-4 β_eff for 20 years, and achieving a burnup of 9.9-10.8 MW × day/kg U in the steady state mode of fuel reloading. Spent fuel immediately after unloading from the reactor can be placed in a transport packaging cask for shipping or safely stored in dry storage at the research reactor site.

• Nuclear Engineering and Technology
• /
• v.54 no.8
• /
• pp.2792-2800
• /
• 2022

In the core of the WWR-K reactor, a long-term irradiation of tristructural isotopic (TRISO)-coated fuel particles (CFPs) with a UO₂ kernel was carried out under high-temperature gas-cooled reactor (HTGR)-like operating conditions. The temperature of this TRISO fuel during irradiation varied in the range of 950-1100 ℃. A fission per initial metal atom (FIMA) of uranium burnup of 9.9% was reached. The release of gaseous fission products was measured in-pile. The release-to-birth ratio (R/B) for the fission product isotopes was calculated. Aspects of fuel safety while achieving deep fuel burnup are important and relevant, including maintaining the integrity of the fuel coatings. The main mechanisms of fuel failure are kernel migration, silicon carbide corrosion by palladium, and gas pressure increase inside the CFP. The formation of gaseous fission products and carbon monoxide leads to an increase in the internal pressure in the CFP, which is a dominant failure mechanism of the coatings under this level of burnup. Irradiated fuel compacts were subjected to electric dissociation to isolate the CFPs from the fuel compacts. In addition, nondestructive methods, such as X-ray radiography and gamma spectrometry, were used. The predicted R/B ratio was evaluated using the fission gas release model developed in the high-temperature test reactor (HTTR) project. In the model, both the through-coatings of failed CFPs and as-fabricated uranium contamination were assumed to be sources of the fission gas. The obtained R/B ratio for gaseous fission products allows the finalization and validation of the model for the release of fission products from the CFPs and fuel compacts. The success of the integrity of TRISO fuel irradiated at approximately 9.9% FIMA was demonstrated. A low fuel failure fraction and R/B ratios indicated good performance and reliability of the studied TRISO fuel.

• Nuclear Engineering and Technology
• /
• v.54 no.3
• /
• pp.940-947
• /
• 2022

Meta-heuristic algorithms have found their place in optimization problems. Henry gas solubility optimization (HGSO) is one of the newest population-based algorithms. This algorithm is inspired by Henry's law of physics. To evaluate the performance of a new algorithm, it must be used in various problems. On the other hand, the optimization of the proportional-integral-derivative (PID) gains for load-following of a nuclear power plant (NPP) is a good challenge to assess the performance of HGSO. Accordingly, the power control of a pressurized water reactor (PWR) is targeted, based on the point kinetics model with six groups of delayed-neutron precursors. In any optimization problem based on meta-heuristic algorithms, an efficient objective function is required. Therefore, the integral of the time-weighted square error (ITSE) performance index is utilized as the objective (cost) function of HGSO, which is constrained by a stability criterion in steady-state operations. A Lyapunov approach guarantees this stability. The results show that this method provides superior results compared to an empirically tuned PID controller with the least error. It also achieves good accuracy compared to an established GA-tuned PID controller.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1540-1553
• /
• 2019

A nuclear power plant (NPP) is a highly complex system-of-systems as manifested through its internal systems interdependence. The negative impact of such interdependence was demonstrated through the 2011 Fukushima Daiichi nuclear disaster. As such, there is a critical need for new strategies to overcome the limitations of current risk assessment techniques (e.g. the use of static event and fault tree schemes), particularly through simulation of the nonlinear dynamic feedback mechanisms between the different NPP systems/components. As the first and key step towards developing an integrated NPP dynamic probabilistic risk assessment platform that can account for such feedback mechanisms, the current study adopts a system dynamics simulation approach to model the thermal dynamic processes in: the reactor core; the secondary coolant system; and the pressurized water reactor. The reactor core and secondary coolant system parameters used to develop system dynamics models are based on those of the Palo Verde Nuclear Generating Station. These three system dynamics models are subsequently validated, using results from published work, under different system perturbations including the change in reactivity, the steam valve coefficient, the primary coolant flow, and others. Moving forward, the developed system dynamics models can be integrated with other interacting processes within a NPP to form the basis of a dynamic system-level (systemic) risk assessment tool.

• Nuclear Engineering and Technology
• /
• v.52 no.3
• /
• pp.469-476
• /
• 2020

Research reactors in-core experimental facilities are designed to provide the highest steady state flux for user's irradiation requirements. However, fuel conversion from highly enriched uranium (HEU) to low enriched uranium (LEU) driven by the ongoing effort to diminish proliferation risk, will impact reactor physics parameters. Preserving the reactor capability to produce the needed flux to perform its intended research functions, determines the conversion feasibility. This study investigates the neutron flux in the central experimental facility of two material test reactors (MTR), the IAEA generic10 MW benchmark reactor and the 22 MW s Egyptian Test and Research Reactor (ETRR-2). A 3D full core model with three uranium enrichment of 93%, 45%, and 20% was constructed utilizing the OpenMC particle transport Monte Carlo code. Neutronics calculations were performed for fresh fuel, the beginning of life cycle (BOL) and end of life cycle (EOL) for each of the three enrichments for both the IAEA 10 MW generic reactor and core 1/98 of the ETRR-2 reactor. Criticality calculations of the effective multiplication factor (K_eff) were executed for each of the twelve cases; results show a reasonable agreement with published benchmark values for both reactors. The thermal, epithermal and fast neutron fluxes were tallied across the core, utilizing the mesh tally capability of the code and are presented here. The axial flux in the central experimental facility was tallied at 1 cm intervals, for each of the cases; results for IAEA 10 MW show a maximum reduction of 14.32% in the thermal flux of LEU to that of the HEU, at EOL. The reduction of the thermal flux for fresh fuel was between 5.81% and 9.62%, with an average drop of 8.1%. At the BOL the thermal flux showed a larger reduction range of 6.92%-13.58% with an average drop of 10.73%. Furthermore, the fission reaction rate was calculated, results showed an increase in the peak fission rate of the LEU case compared to the HEU case. Results for the ETRR-2 reactor show an average increase of 62.31% in the thermal flux of LEU to that of the HEU due to the effect of spectrum hardening. The fission rate density increased with enrichment, resulting in 34% maximum increase in the HEU case compared to the LEU case at the assemblies surrounding the flux trap.

• Nuclear Engineering and Technology
• /
• v.54 no.7
• /
• pp.2670-2689
• /
• 2022

STREAM - a lattice transport calculation code with method of characteristics for the purpose of light water reactor analysis - has been developed by the Computational Reactor Physics and Experiment laboratory (CORE) of the Ulsan National Institute of Science and Technology (UNIST). Recently, efforts have been taken to develop a photon module in STREAM to assess photon heating and the influence of gamma photon transport on power distributions, as only neutron transport was considered in previous STREAM versions. A multi-group photon library is produced for STREAM based on the ENDF/B-VII.1 library with the use of the library-processing code NJOY. The developed photon solver for the computation of 2D and 3D distributions of photon flux and energy deposition is based on the method of characteristics like the neutron solver. The photon library and photon module produced and implemented for STREAM are verified on VERA pin and assembly problems by comparison with the Monte Carlo code MCS - also developed at UNIST. A short analysis of the impact of photon transport during depletion and thermal hydraulics feedback is presented for a 2D core also from the VERA benchmark.

• Nuclear Engineering and Technology
• /
• v.52 no.8
• /
• pp.1603-1610
• /
• 2020

Due to ever-growing advancements in computers and relatively easy access to them, many efforts have been made to develop high-fidelity, high-performance, multi-physics tools, which play a crucial role in the design and operation of nuclear reactors. For this purpose in this study, the neutronic Monte Carlo and thermal-hydraulic sub-channel codes entitled MCNP and COBRA-EN, respectively, were applied for external coupling with each other. The coupled code was validated by code-to-code comparison with the internal couplings between MCNP5 and SUBCHANFLOW as well as MCNP6 and CTF. The simulation results of all code systems were in good agreement with each other. Then, as the second problem, the core of the VVER-1000 v446 reactor was simulated by the MCNP4C/COBRA-EN coupled code to measure the capability of the developed code to calculate the neutronic and thermohydraulic parameters of real and industrial cases. The simulation results of VVER-1000 core were compared with FSAR and another numerical solution of this benchmark. The obtained results showed that the ability of the MCNP4C/COBRA-EN code for estimating the neutronic and thermohydraulic parameters was very satisfactory.