• Nuclear Engineering and Technology
• /
• 제51권5호
• /
• pp.1209-1217
• /
• 2019

Due to the limitation of the computers, the conventional homogenization method is based on many assumptions and approximations, and some tough problems such as energy spectrum and boundary condition are faced. To deal with those problems, the Monte Carlo global homogenization is adopted. The Reactor Monte Carlo code RMC is used to study the global homogenization method, and the one-step global homogenization method is proposed. The superimposed mesh geometry is also used to divide the physical models, leading to better geometric flexibility. A set of multigroup homogenization cross sections is online generated for each mesh under the real neutron energy spectrum and boundary condition, the cross sections are adjusted by the superhomogenization method, and no leakage correction is required. During the process of superhomogenization, the author-developed reactor core program NLSP3 is used for global calculation, so the global flux distribution and equivalent homogenization cross sections could be solved simultaneously. Meanwhile, the calculated homogenization cross section could accurately reconstruct the non-homogenization flux distribution and could also be used for fine calculation. This one-step global homogenization method was tested by a PWR assembly and a small reactor model, and the results show the validity.

• Nuclear Engineering and Technology
• /
• 제54권7호
• /
• pp.2650-2659
• /
• 2022

The dense granular flow spallation target is a new target concept proposed for an Accelerator-Driven Sub-critical (ADS) system. In this paper, the beam-target configurations of a tungsten granular flow target for the ADS with a thermal power of 1 GW is explored. The beam profile options using different scanning methods are discussed. The critical geometry parameters are adjusted to investigate the performance of the granular target from the aspects of neutron efficiency, stability and temperature distribution in target medium. To figure out how the target under accident conditions would behave, different clogging conditions are induced in the simulation. The dynamic processes are analyzed and some important parameters such as abnormal temperature rise and beam cutoff time window are obtained. The response of the sub-critical reactor to a clogging accident is also investigated. It is indicated that the monitoring of the granular flow by the neutron detectors in the sub-critical core will be effective.

• Nuclear Engineering and Technology
• /
• 제54권7호
• /
• pp.2395-2407
• /
• 2022

The sodium-cooled fast reactor is among the innovative nuclear technologies selected in the framework of the development of Generation IV concepts, allowing the irradiation of uranium-plutonium mixed oxide fuels (MOX). A fundamental step for the safety assessment of MOX-fuelled pins for fast reactor applications is the evaluation, by means of fuel performance codes, of the integral thermal-mechanical behaviour under irradiation, involving the fission gas behaviour and release in the fuel-cladding gap. This work is dedicated to the performance analysis of an inner-core fuel pin representative of the ASTRID sodium-cooled concept design, selected as case study for the benchmark between the GERMINAL and TRANSURANUS fuel performance codes. The focus is on fission gas-related mechanisms and integral outcomes as predicted by means of the SCIANTIX module (allowing the physics-based treatment of inert gas behaviour and release) coupled to both fuel performance codes. The benchmark activity involves the application of both GERMINAL and TRANSURANUS in their "pre-INSPYRE" versions, i.e., adopting the state-of-the-art recommended correlations available in the codes, compared with the "post-INSPYRE" code results, obtained by implementing novel models for MOX fuel properties and phenomena (SCIANTIX included) developed in the framework of the INSPYRE H2020 Project. The SCIANTIX modelling includes the consideration of burst releases of the fission gas stored at the grain boundaries occurring during power transients of shutdown and start-up, whose effect on a fast reactor fuel concept is analysed. A clear need to further extend and validate the SCIANTIX module for application to fast reactor MOX emerges from this work; nevertheless, the GERMINAL-TRANSURANUS benchmark on the ASTRID case study highlights the achieved code capabilities for fast reactor conditions and paves the way towards the proper application of fuel performance codes to safety evaluations on Generation IV reactor concepts.

• Nuclear Engineering and Technology
• /
• 제41권7호
• /
• pp.885-892
• /
• 2009

This manuscript will discuss a numerical method where the six equations of two-phase flow, the solid heat conduction equations, and the two equations that describe neutron diffusion and precursor concentration are solved together in a tightly coupled, nonlinear fashion for a simplified model of a nuclear reactor core. This approach has two important advantages. The first advantage is a higher level of accuracy. Because the equations are solved together in a single nonlinear system, the solution is more accurate than the traditional "operator split" approach where the two-phase flow equations are solved first, the heat conduction is solved second and the neutron diffusion is solved third, limiting the temporal accuracy to $1^{st}$ order because the nonlinear coupling between the physics is handled explicitly. The second advantage of the method described in this manuscript is that the time step control in the fully implicit system can be based on the timescale of the solution rather than a stability-based time step restriction like the material Courant limit required of operator-split methods. In this work, a pilot code was used which employs this tightly coupled, fully implicit method to simulate a reactor core. Results are presented from a simulated control rod movement which show $2^{nd}$ order accuracy in time. Also described in this paper is a simulated rod ejection demonstrating how the fastest timescale of the problem can change between the state variables of neutronics, conduction and two-phase flow during the course of a transient.

• Nuclear Engineering and Technology
• /
• 제52권5호
• /
• pp.869-877
• /
• 2020

A 24-month Advanced Power Reactor 1400 (APR1400) core with a very-low-boron (VLB) concentration has been investigated for an inherently safe and high-performance PWR in this work. To develop a high-performance APR1400 which is able to do the passive frequency control operation, VLB feature is essential. In this paper, the centrally-shielded burnable absorber (CSBA) is utilized for an efficient VLB operation in the 24-month cycle APR1400 core. This innovative design of the VLB APR1400 core includes the optimization of burnable absorber and loading pattern as well as axial cutback for a 24-month cycle operation. In addition to CSBA, an Er-doped guide thimble is also introduced for partial management of the excess reactivity and local peaking factor. To improve the neutron economy of the core, two alternative radial reflectors are adopted in this study, which are SS-304 and ZrO₂. The core reactivity and power distributions for a 2-batch equilibrium cycle are analyzed and compared for each reflector design. Numerical results show that a VLB core can be successfully designed with 24-month cycle and the cycle length is improved significantly with the alternative reflectors. The neutronic analyses are performed using the Monte Carlo Serpent code and 3-D diffusion code COREDAX-2 with the ENDF/B-VII.1.

• Nuclear Engineering and Technology
• /
• 제54권8호
• /
• pp.3059-3072
• /
• 2022

A Jacobian-Free Newton Krylov Two-Nodal Coarse Mesh Finite Difference algorithm based on Nodal Expansion Method (NEM_TNCMFD_JFNK) is successfully developed and proposed to solve the three-dimensional (3D) and multi-group reactor physics models. In the NEM_TNCMFD_JFNK method, the efficient JFNK method with the Modified Incomplete LU (MILU) preconditioner is integrated and applied into the discrete systems of the NEM-based two-node CMFD method by constructing the residual functions of only the nodal average fluxes and the eigenvalue. All the nonlinear corrective nodal coupling coefficients are updated on the basis of two-nodal NEM formulation including the discontinuity factor in every few newton steps. All the expansion coefficients and interface currents of the two-node NEM need not be chosen as the solution variables to evaluate the residual functions of the NEM_TNCMFD_JFNK method, therefore, the NEM_TNCMFD_JFNK method can greatly reduce the number of solution variables and the computational cost compared with the JFNK based on the conventional NEM. Finally the NEM_TNCMFD_JFNK code is developed and then analyzed by simulating the representative PWR MOX/UO₂ core benchmark, the popular NEACRP 3D core benchmark and the complicated full-core pin-by-pin homogenous core model. Numerical solutions show that the proposed NEM_TNCMFD_JFNK method with the MILU preconditioner has the good numerical accuracy and can obtain higher computational efficiency than the NEM-based two-node CMFD algorithm with the power method in the outer iteration and the Krylov method using the MILU preconditioner in the inner iteration, which indicates the NEM_TNCMFD_JFNK method can serve as a potential and efficient numerical tool for reactor neutron diffusion analysis module in the JFNK-based multiphysics coupling application.

• Nuclear Engineering and Technology
• /
• 제55권2호
• /
• pp.452-459
• /
• 2023

The core power control is an important issue for the study of dynamic characteristics in China initiative accelerator driven subcritical system (CiADS), which has direct impact on the control strategy and safety analysis process. The CiADS is an experimental facility that is only controlled by the proton beam intensity without considering the control rods in the current engineering design stage. In order to get the optimized operation scheme with the stable and reliable features, the variation of beam intensity using the continuous and periodic control approaches has been adopted, and the change of collimator and the adjusting of duty ratio have been proposed in the power control process. Considering the neutronics and the thermal-hydraulics characteristics in CiADS, the physical model for the core power control has been established by means of the point reactor kinetics method and the lumped parameter method. Moreover, the multi-inputs single-output (MISO) logical structure for the power control process has been constructed using proportional integral derivative (PID) controller, and the meta-heuristic algorithm has been employed to obtain the global optimized parameters for the stable running mode without producing large perturbations. Finally, the verification and validation of the control method have been tested based on the reference scenarios in considering the disturbances of spallation neutron source and inlet temperature respectively, where all the numerical results reveal that the optimization method has satisfactory performance in the CiADS core power control scenarios.

• Nuclear Engineering and Technology
• /
• 제52권9호
• /
• pp.1896-1906
• /
• 2020

The BEAVRS (Benchmark for Evaluation and Validation of Reactor Simulation) benchmark calculations were performed by DeCART stand-alone and DeCART/MATRA multi-physics coupled code system to verify their accuracy. The solutions of DeCART stand-alone calculations for the control rod bank worth, detector signal, isothermal temperature coefficient, and critical boron concentration agreed very well with the measurements. The root-mean-square errors of the boron letdown curves for two-cycles were less than about 20 ppm, while the individual and total control rod bank worth agreed well within 7.3% and 2.4%, respectively. For the BEAVRS benchmark calculations at the beginning of burnup, the difference between DeCART simplified thermal-hydraulic stand-alone and DeCART/MATRA coupled calculations were not significantly large. Therefore, it is concluded that both the DeCART stand-alone code and the DeCART/MATRA multi-physics coupled code system have the capabilities to generate high fidelity transport solutions at core follow calculations.

• Nuclear Engineering and Technology
• /
• 제41권5호
• /
• pp.617-648
• /
• 2009

Under the government supported long-term nuclear R&D program, the severe accident research program at KAERI is directed to investigate unresolved severe accident issues such as core debris coolability, steam explosions, and hydrogen combustion both experimentally and numerically. Extensive studies have been performed to evaluate the in-vessel retention of core debris through external reactor vessel cooling concept for APR1400 as a severe accident management strategy. Additionally, an improvement of the insulator design outside the vessel was investigated. To address steam explosions, a series of experiments using a prototypic material was performed in the TROI facility. Major parameters such as material composition and void fraction as well as the relevant physics affecting the energetics of steam explosions were investigated. For hydrogen control in Korean nuclear power plants, evaluation of the hydrogen concentration and the possibility of deflagration-to-detonation transition occurrence in the containment using three-dimensional analysis code, GASFLOW, were performed. Finally, the integrated severe accident analysis code, MIDAS, has been developed for domestication based on MELCOR. The data transfer scheme using pointers was restructured with the modules and the derived-type direct variables using FORTRAN90. New models were implemented to extend the capability of MIDAS.

• Nuclear Engineering and Technology
• /
• 제48권3호
• /
• pp.660-672
• /
• 2016

In Part I of this paper, the two-temperature homogenized model for the fully ceramic microencapsulated fuel, in which tristructural isotropic particles are randomly dispersed in a fine lattice stochastic structure, was discussed. In this model, the fuel-kernel and silicon carbide matrix temperatures are distinguished. Moreover, the obtained temperature profiles are more realistic than those obtained using other models. Using the temperature-dependent thermal conductivities of uranium nitride and the silicon carbide matrix, temperature-dependent homogenized parameters were obtained. In Part II of the paper, coupled with the COREDAX code, a reactor core loaded by fully ceramic microencapsulated fuel in which tristructural isotropic particles are randomly dispersed in the fine lattice stochastic structure is analyzed via a two-temperature homogenized model at steady and transient states. The results are compared with those from harmonic- and volumetric-average thermal conductivity models; i.e., we compare $k_{eff}$ eigenvalues, power distributions, and temperature profiles in the hottest single channel at a steady state. At transient states, we compare total power, average energy deposition, and maximum temperatures in the hottest single channel obtained by the different thermal analysis models. The different thermal analysis models and the availability of fuel-kernel temperatures in the two-temperature homogenized model for Doppler temperature feedback lead to significant differences.