• Nuclear Engineering and Technology
• /
• 제55권1호
• /
• pp.310-323
• /
• 2023

The k-eigenvalue neutronics/thermal-hydraulics coupling calculation is a key issue for reactor design and analysis. Jacobian-free Newton-Krylov (JFNK) method, featured with super-linear convergence rate and high efficiency, has been attracting more and more attention to solve the multi-physics coupling problem. However, it may converge to the high-order eigenmode because of the multiple solutions nature of the k-eigenvalue form of multi-physics coupling issue. Based on our previous work, a modified JFNK with a line search method is proposed in this work, which can find the fundamental eigenmode together with thermal-hydraulics feedback in a wide range of initial values. In detail, the existing modified JFNK method is combined with the line search strategy, so that the intermediate iterative solution can avoid a sudden divergence and be adjusted into a convergence basin smoothly. Two simplified 2-D homogeneous reactor models, a PWR model, and an HTR model, are utilized to evaluate the performance of the newly proposed JFNK method. The results show that the performance of this proposed JFNK is more robust than the existing JFNK-based methods.

• Nuclear Engineering and Technology
• /
• 제42권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
• /
• 제55권12호
• /
• pp.4685-4694
• /
• 2023

The ability to calculate the material density sensitivity coefficients of power with respect to the material density has broad application prospects for accelerating Monte Carlo-Thermal Hydraulics iterations. The second-order material density sensitivity coefficients for the general Monte Carlo score have been derived based on the differential operator sampling method in this paper, and the calculation of the sensitivity coefficients of cell power scores with respect to the material density has been realized in continuous-energy Monte Carlo code RMC. Based on the power-density sensitivity coefficients, the sensitivity coefficients of power scores to some other physical quantities, such as power-boron concentration coefficients and power-temperature coefficients considering only the thermal expansion, were subsequently calculated. The effectiveness of the proposed method is demonstrated in the power-density coefficients problems of the pressurized water reactor (PWR) moderator and the heat pipe reactor (HPR) reflectors. The calculations were carried out using RMC and the ENDF/B-VII.1 neutron nuclear data. It is shown that the calculated sensitivity coefficients can be used to predict the power scores accurately over a wide range of boron concentration of the PWR moderator and a wide range of temperature of HPR reflectors.

• Nuclear Engineering and Technology
• /
• 제54권6호
• /
• pp.2094-2106
• /
• 2022

The solid-state core of a heat pipe cooled reactor operates at high temperatures over 1000 K with thermal and irradiation-induced expansion during burnup. The expansion changes the gap thickness between the solid components and the material properties, and may even cause the gap closure, which then significantly influences the thermal and mechanical characteristics of the reactor core. This study developed an irradiation behavior model for HPRTRAN, a heat pipe reactor system analysis code, to introduce the irradiation effects such as swelling and creep. The megawatt heat pipe reactor MegaPower was chosen as an application case. The coupled irradiation-thermal-mechanical model was developed to simulate the irradiation effects on the heat transfer and stresses of the whole reactor core. The results show that the irradiation deformation effect is significant, with the irradiation-induced strains up to 2.82% for fuel and 0.30% for monolith at the end of the reactor lifetime. The peak temperatures during the lifetime are 1027:3 K for the fuel and 956:2 K for monolith. The gap closure enhances the heat transfer but caused high stresses exceeding the yield strength in the monolith.

• Nuclear Engineering and Technology
• /
• 제51권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
• /
• 제50권7호
• /
• pp.1017-1023
• /
• 2018

The measurement and calculation of neutronic parameters in nuclear research reactors has an important influence on control and safety of the nuclear reactor. The power peaking factors, reactivity coefficients and kinetic parameters are the most important neutronic parameter for determining the state of the reactor. The position of the control shim safety rods in the core configuration affects these parameters. The main purpose of this work is to use the MTR_PC package to evaluate the effect of the partially insertion of the control rod on the neutronic parameters at the operating core of the Tehran Research Reactor. The simulation results show that by increasing the insertion of control rods (bank) in the core, the absolute values of power peaking factor, reactivity coefficients and effective delayed neutron fraction increased and only prompt neutron life time decreased. In addition, the results show that the changes of moderator temperature coefficients value versus the control rods positions are very significant. The average value of moderator temperature coefficients increase about 98% in the range of 0-70% insertion of control rods.

• Nuclear Engineering and Technology
• /
• 제53권2호
• /
• pp.423-429
• /
• 2021

The investigation of the utilization of enriched ²⁰⁸Pb as a coolant to enhance the performance of a long-life fast reactor with a Modified CANDLE (Constant Axial shape of Neutron flux, nuclide densities, and power shape During Life of Energy production) burnup scheme has performed. The analyzes were performed on a reactor with thermal power of 800 MegaWatt Thermal (MWTh) with a refueling process every 15 years. Uranium Nitride (enriched ¹⁵N), ²⁰⁸Pb, and High-Cr martensitic steel HT-9 were employed as fuel, coolant, and cladding materials, respectively. One of the Pb-nat isotopes, ²⁰⁸Pb, has the smallest neutron capture cross-section (0.23 mb) among other liquid metal coolants. Furthermore, the neutron-producing cross-section (n, 2n) of ²⁰⁸Pb is larger than sodium (Na). On the other hand, the inelastic scattering energy threshold of ²⁰⁸Pb is the highest among Na, ^natPb, and Bi. The small inelastic scattering cross-section of ²⁰⁸Pb can harden the neutron energy spectrum. Therefore, ²⁰⁸Pb is a better neutron multiplier than any other liquid metal coolant. The excess neutrons cause more production than consumption of ²³⁹Pu. Hence, it can reduce the initial fuel loading of the reactor. The selective photoreaction process was developing to obtain enriched ²⁰⁸Pb. The neutronic was calculated using SRAC and JENDL 4.0 as a nuclear data library. We obtained that the modified CANDLE reactor with enriched ²⁰⁸Pb as coolant and reflector has the highest k-eff among all reactors. Meanwhile, the ^natPb cooled reactor has the lowest k-eff. Thus, the utilization of the enriched ²⁰⁸Pb as the coolant can reduce reactor initial fuel loading. Moreover, the enriched ²⁰⁸Pb-cooled reactor has the smallest power peaking factor among all reactors. Therefore, the enriched ²⁰⁸Pb can enhance the performance of a long-life Modified CANDLE fast reactor.

• Nuclear Engineering and Technology
• /
• 제54권8호
• /
• pp.2877-2882
• /
• 2022

The paper presents the main results of computational and experimental optimization of the automatic power control system (AC) of the IBR-2M pulsed reactor in the presence of a high level of oscillatory instability. Optimization of the parameters of the AC made it possible to significantly reduce the influence of random and deterministic oscillations of reactivity on the noise of the pulse energy, as well as to sharply reduce the manifestation of the oscillatory instability of the reactor. As a result, the safety and reliability of operation of the reactor has increased substantially.

• Nuclear Engineering and Technology
• /
• 제56권3호
• /
• pp.980-992
• /
• 2024

Steam line break accident (SLB) in the nuclear reactor is one of the representative Non-LOCA accidents in which thermal-hydraulics and neutron kinetics are strongly coupled each other. Thus, the multi-scale and multi-physics approach is applied in this study in order to examine a realistic safety margin. An entire reactor coolant system is modelled by system scale node, whereas sub-channel scale resolution is applied for the region of interest such as the reactor core. Fuel performance code is extended to consider full core pin-wise fuel behaviour. The MARU platform is developed for easy integration of the codes to be coupled. An initial stage of the steam line break accident is simulated on the MARU platform. As cold coolant is injected from the cold leg into the reactor pressure vessel, the power increases due to the moderator feedback. Three-dimensional coolant and fuel behaviour are qualitatively visualized for easy comprehension. Moreover, quantitative investigation is added by focusing on the enhancement of safety margin by means of comparing the minimum departure from nucleate boiling ratio (MDNBR). Three factors contributing to the increase of the MDNBR are proposed: Various geometric parameters, realistic power distribution by neutron kinetics code, Radial coolant mixing including sub-channel physics model.

• Nuclear Engineering and Technology
• /
• 제53권9호
• /
• pp.2830-2838
• /
• 2021

This work deals with the assessment of the burnup capabilities of the Serpent Monte Carlo code to predict spent nuclear fuel (SNF) isotopic concentrations for low-enriched uranium (LEU) fuel at different burnup levels up to 47 MWd/kgU. The irradiation of six UO₂ experimental samples in three different VVER-1000 reactor units has been simulated and the predicted concentrations of actinides up to ²⁴⁴Cm have been compared with the corresponding measured values. The results show a global good agreement between calculated and experimental concentrations, in several cases within the margins of the nuclear data uncertainties and in a few cases even within the reported experimental uncertainties. The differences in the performances of the JEFF3.1.1, ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear data libraries (NDLs) have also been assessed and the use of the newly released ENDF/B-VIII.0 library has shown an increased accuracy in the prediction of the C/E's for some of the actinides considered, particularly for the plutonium isotopes. This work represents a step forward towards the validation of advanced simulation tools against post irradiation experimental data and the obtained results provide an evidence of the capabilities of the Serpent Monte-Carlo code with the associated modern NDLs to accurately compute SNF nuclide inventory concentrations for VVER-1000 type reactors.