• Nuclear Engineering and Technology
• /
• v.55 no.5
• /
• pp.1775-1782
• /
• 2023

Several reactor physics commissioning experiments were conducted to obtain the neutronic parameters at the beginning of the G.A. Siwabessy Multi-purpose Reactor (RSG-GAS) operation. These parameters are essential for the reactor to safety operate. Leveraging the experimental data, this study evaluated the calculated core reactivity, control rod reactivity worth, integral control rod reactivity curve, and fuel reactivity. Calculations were carried out with Serpent 2 code using the latest neutron cross-section data ENDF/B-VIII.0. The criticality calculations were carried out for the RSG-GAS first core up to the third core configuration, which has been done experimentally during these commissioning periods. The excess reactivity for the second and third cores showed a difference of 510.97 pcm and 253.23 pcm to the experiment data. The calculated integral reactivity of the control rod has an error of less than 1.0% compared to the experimental data. The calculated fuel reactivity value is consistent with the measured data, with a maximum error of 2.12%. Therefore, it can be concluded that the RSG-GAS reactor core model is in good agreement to reproduce excess reactivity, control rod worth, and fuel element reactivity.

• Nuclear Engineering and Technology
• /
• v.53 no.11
• /
• pp.3665-3674
• /
• 2021

One of the important parts of the in-core fuel management is loading pattern optimization (LPO). The loading pattern optimization as a reasonable design of the in-core fuel management can improve both economic and safe aspects of the nuclear reactor. This work proposes the hybrid of fuzzy logic controller with harmony search algorithm (HS) for loading pattern optimization in a pressurized water reactor. The music improvisation process to find a pleasing harmony is inspiring the harmony search algorithm. In this work, the adjustment of the harmony search algorithm parameters such as the bandwidth and the pitch adjustment rate are increasing performance of the proposed algorithm which is done through a fuzzy logic controller. Hence, membership functions and fuzzy rules are designed to improve the performance of the HS algorithm and achieve optimal results. The objective of the method is finding an optimum core arrangement according to safety and economic aspects such as reduction of power peaking factor (PPF) and increase of effective multiplication factor (Keff). The proposed approach effectiveness has been tried in two cases, Michalewicz's bivariate function problem and NEACRP LWR core. The results show that by using fuzzy harmony search algorithm the value of the fitness function is improved by 15.35%. Finally, with regard to the new solutions proposed in this research it could be used as a trustworthy method for other optimization issues of engineering field.

• Nuclear Engineering and Technology
• /
• v.53 no.10
• /
• pp.3127-3132
• /
• 2021

Feasibility study of burn-up analysis and monitoring using delayed fast neutrons was investigated at Missouri University of Science and Technology Reactor (MSTR). Burnt and fresh fuel elements were used to collect delayed fast neutron data for different power levels. Total reactivity varied depending on the burn-up rate of fuel elements for each core configuration. The regulating rod worth was 2.07E-04 𝚫k/k/in and 1.95E-04 𝚫k/k/in for T121 and T122 core configurations at 11 inch, respectively. Delayed fast neutron spectrum of F1 (burnt) and F16 (fresh) fuel elements were analyzed further, and a strong correlation was observed between delayed fast neutron emission and burn-up. According to the analyzed peaks in burnt and fresh fuels, reactor power dependency was observed and it was determined that delayed neutron provided more reliable results at reactor powers of 50 kW and above.

• Nuclear Engineering and Technology
• /
• v.50 no.1
• /
• pp.35-42
• /
• 2018

This article presents clean core criticality calculations and control rod worth calculations for TRIGA (Training, Research, Isotope production-General Atomics) Mark II research reactor benchmark cores using Winfrith Improved Multi-group Scheme-D/4 (WIMS-D/4) and Program for Reactor In-core Analysis using Diffusion Equation (PRIDE) codes. Cores 133 and 134 were analyzed in 2-D (r, ${\theta}$) and 3-D (r, ${\theta}$, z), using WIMS-D/4 and PRIDE codes. Moreover, the influence of cross-section data was also studied using various libraries based on Evaluated Nuclear Data File (ENDF/B-VI.8 and VII.0), Joint Evaluated Fission and Fusion File (JEFF-3.1), Japanese Evaluated Nuclear Data Library (JENDL-3.2), and Joint Evaluated File (JEF-2.2) nuclear data. The simulation results showed that the multiplication factor calculated for all these data libraries is within 1% of the experimental results. The reactivity worth of the control rods of core 134 was also calculated with different homogenization approaches. A comparison was made with experimental and reported Monte Carlo results, and it was found that, using proper homogenization of absorber regions and surrounding fuel regions, the results obtained with PRIDE code are significantly improved.

• Nuclear Engineering and Technology
• /
• v.53 no.2
• /
• pp.455-465
• /
• 2021

The Molten Salt Reactor (MSR), one of the six advanced reactor types of the 4th generation nuclear energy systems, has many impressive features including economic advantages, inherent safety and nuclear non-proliferation. This paper introduces a system analysis code named TREND, which is developed and used for the steady and transient simulation of MSRs. The TREND code calculates the distributions of pressure, velocity and temperature of single-phase flows by solving the conservation equations of mass, momentum and energy, along with a fluid state equation. Heat structures coupled with the fluid dynamics model is sufficient to meet the demands of modeling MSR system-level thermal-hydraulics. The core power is based on the point reactor neutron kinetics model calculated by the typical Runge-Kutta method. An incremental PID controller is inserted to adjust the operation behaviors. The verification and validation of the TREND code have been carried out in two aspects: detailed code-to-code comparison with established thermal-hydraulic system codes such as RELAP5, and validation with the experimental data from MSRE and the CIET facility (the University of California, Berkeley's Compact Integral Effects Test facility).The results indicate that TREND can be used in analyzing the transient behaviors of MSRs and will be improved by validating with more experimental results with the support of SINAP.

• Journal of Radiation Protection and Research
• /
• v.41 no.2
• /
• pp.155-159
• /
• 2016

Background: Nuclear reactors produce a great number of antielectron neutrinos mainly from beta-decay chains of fission products. Such neutrinos have energies mostly in MeV range. We are interested in neutrinos in a region of keV, since they may take part in special weak interactions. We calculate reactor antineutrino spectra especially in the low energy region. In this work we present neutrino spectrum from a typical pressurized water reactor (PWR) reactor core. Materials and Methods: To calculate neutrino spectra, we need information about all generated nuclides that emit neutrinos. They are mainly fission fragments, reaction products and trans-uranium nuclides that undergo negative beta decay. Information in relation to trans-uranium nuclide compositions and its evolution in time (burn-up process) were provided by a reactor code MVP-BURN. We used typical PWR parameter input for MVP-BURN code and assumed the reactor to be operated continuously for 1 year (12 months) in a steady thermal power (3.4 GWth). The PWR has three fuel compositions of 2.0, 3.5 and 4.1 wt% $^{235}U$ contents. For preliminary calculation we adopted a standard burn-up chain model provided by MVP-BURN. The chain model treated 21 heavy nuclides and 50 fission products. The MVB-BURN code utilized JENDL 3.3 as nuclear data library. Results and Discussion: We confirm that the antielectron neutrino flux in the low energy region increases with burn-up of nuclear fuel. The antielectron-neutrino spectrum in low energy region is influenced by beta emitter nuclides with low Q value in beta decay (e.g. $^{241}Pu$) which is influenced by burp-up level: Low energy antielectron-neutrino spectra or emission rates increase when beta emitters with low Q value in beta decay accumulate Conclusion: Our result shows the flux of low energy reactor neutrinos increases with burn-up of nuclear fuel.

• Nuclear Engineering and Technology
• /
• v.51 no.1
• /
• pp.116-124
• /
• 2019

$M^2LFR-1000$ is a medium-power modular lead-cooled fast reactor, developed by University of Science and Technology of China (USTC), aiming at achieving a reactor design fulfilling the Gen IV nuclear system requirements and meanwhile emphasizing the optimum safety and economics. In order to evaluate the safety performance of $M^2LFR-1000$ reactor core, three typical transients are selected from initiating events, which are unprotected transient overpower (UTOP), unprotected loss of offsite power (ULOHS+ULOF) and increase of feedwater flowrate with primary pumps trip (IFW+PLOF). These three transients presented and discussed in this paper are performed with the code Analysis of THermal-hydraulics of LEaks and Transients (ATHLET), which is developed by Gesellschaft $f{\ddot{u}}r$ Anlagen-und Reaktorsicherheit gGmbH (GRS). The results indicate that the $M^2LFR$ is safe enough under these three transients due to the good inherent safety features of the reactor, without human intervention, the reactor will reach a new steady state under UTOP condition.

• Nuclear Engineering and Technology
• /
• v.37 no.1
• /
• pp.25-78
• /
• 2005

As a key and core knowledge for the design of various types of nuclear reactors, the discipline of reactor physics has been advanced continually in the past six decades and has led to a very sophisticated fabric of analysis methods and computer codes in use today. Notwithstanding, the discipline faces interesting challenges from next-generation nuclear reactors and innovative new fuel designs in the coming. After presenting a brief overview of important tasks and steps involved in the nuclear design and analysis of a reactor, this article focuses on the currently-used design and analysis methods, issues and limitations, and current activities to resolve them as follows: (1) Derivation of the multi group transport equations and the multi group diffusion equations, with representative solution methods thereof. (2) Elements of modem (now almost three decades old) diffusion nodal methods. (3) Limitations of nodal methods such as transverse integration, flux reconstruction, and analysis of UO2-MOX mixed cores. Homogenization and related issues. (4) Description of the analytic function expansion nodal (AFEN) method. (5) Ongoing efforts for three-dimensional whole-core heterogeneous transport calculations and acceleration methods. (6) Elements of spatial kinetics calculation methods and coupled neutronics and thermal-hydraulics transient analysis. (7) Identification of future research and development areas in advanced reactors and Generation-IV reactors, in particular, in very high temperature gas reactor (VHTR) cores.

• Nuclear Engineering and Technology
• /
• v.55 no.10
• /
• pp.3732-3753
• /
• 2023

The Real-time Analysis for Particle-transport and In-situ Detection (RAPID) Code System, developed based on the Multi-stage Response-function Transport (MRT) methodology, enables real-time simulation of nuclear systems such as reactor cores, spent nuclear fuel pools and casks, and sub-critical facilities. This paper presents the application of a novel fission matrix-based burnup methodology to the well-characterized JSI TRIGA Mark II research reactor. This methodology allows for calculation of nuclear fuel depletion by combination and interpolation of RAPID's burnup dependent fission matrix (FM) coefficients to take into account core changes due to burnup. The methodology is compared to experimentally validated Serpent-2 Monte Carlo depletion calculations. The results show that the burnup methodology for RAPID (bRAPID) implemented into RAPID is capable of accurately calculating the k_eff burnup changes of the reactor core as the average discrepancies throughout the whole burnup interval are 37 pcm. Furthermore, capability of accurately describing 3D fission source distribution changes with burnup is demonstrated by having less than 1% relative discrepancies compared to Serpent-2. Good agreement is observed for axially and pin-wise dependent fuel burnup and nuclear fuel nuclide composition as a function of burnup. It is demonstrated that bRAPID accurately describes burnup in areas with high gradients of neutron flux (e.g. vicinity of control rods). Observed discrepancies for some isotopes are explained by analyzing the neutron spectrum. This paper presents a powerful depletion calculation tool that is capable of characterization of spent nuclear fuel on the fly while the reactor is in operation.

• Nuclear Engineering and Technology
• /
• v.50 no.3
• /
• pp.356-367
• /
• 2018

A series of tests dedicated to station blackout (SBO) accident scenarios have been recently performed at the $Prim{\ddot{a}}rkreislauf-Versuchsanlage$ (primary coolant loop test facility; PKL) facility in the framework of the OECD/NEA PKL-3 project. These investigations address current safety issues related to beyond design basis accident transients with significant core heat up. This work presents a detailed analysis using the best estimate thermal-hydraulic code TRACE (v5.0 Patch4) of different SBO scenarios conducted at the PKL facility; failures of high- and low-pressure safety injection systems together with steam generator (SG) feedwater supply are considered, thus calling for adequate accident management actions and timely implementation of alternative emergency cooling procedures to prevent core meltdown. The presented analysis evaluates the capability of the applied TRACE model of the PKL facility to correctly capture the sequences of events in the different SBO scenarios, namely the SBO tests H2.1, H2.2 run 1 and H2.2 run 2, including symmetric or asymmetric secondary side depressurization, primary side depressurization, accumulator (ACC) injection in the cold legs and secondary side feeding with mobile pump and/or primary side emergency core coolant injection from the fuel pool cooling pump. This study is focused specifically on the prediction of the core exit temperature, which drives the execution of the most relevant accident management actions. This work presents, in particular, the key improvements made to the TRACE model that helped to improve the code predictions, including the modeling of dynamical heat losses, the nodalization of SGs' heat exchanger tubes and the ACCs. Another relevant aspect of this work is to evaluate how well the model simulations of the three different scenarios qualitatively and quantitatively capture the trends and results exhibited by the actual experiments. For instance, how the number of SGs considered for secondary side depressurization affects the heat transfer from primary side; how the discharge capacity of the pressurizer relief valve affects the dynamics of the transient; how ACC initial pressure and nitrogen release affect the grace time between ACC injection and subsequent core heat up; and how well the alternative feeding modes of the secondary and/or primary side with mobile injection pumps affect core quenching and ensure stable long-term core cooling under controlled boiling conditions.