Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MULTIGROUP THEORY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Investigations about iso-geometric analysis for self-shielding calculations with the subgroup method

The implementation of a self-shielding method for a neutron transport calculation code based on the iso-geometric analysis (IGA) method that can resolve the multi-group neutron transport equation for arbitrary spatial domain, is presented. The self-shielding model based on the subgroup theory is adopted because the subgroup method can be used to perform calculations for any arbitrary geometrical domain which is the main purpose of our IGA code. Some basic theory of the subgroup method is given. A self-shielding calculation based on the PWR fuel pin composed of MOX fuel is presented and compared with a Monte Carlo calculation. The result is that the combination of SN transport theory, IGA and the subgroup method gives correctly shielded cross sections. Validation of the combination of the IGA method and SN neutron transport for subgroup calculations in two steps are presented: the first step is to ascertain the correctness of the IGA solutions compared to a known analytical solution in diffusion theory; the second step is a validation of the IGA solutions compared to a reference calculation in SN transport theory. The relation between the settings for the IGA solver and the accuracy of the results are elucidated. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

A new approach for resonance treatment of doubly heterogeneous fuel using the RSE method

A new resonance calculation method for the doubly-heterogeneous (DH) fuels such as high- temperature gas-cooled reactor fuel is proposed based on the resonance calculation based on Spectral Expansion (RSE) method. The concept of pointwise disadvantage factor for fuel grain is taken into account to treat the DH fuels. The verification calculation is carried out for simplified single fuel cell and fuel compact consisting of five fuel cells and graphite moderator. The calculation results indicate that the present method can appropriately handle the space- dependent self-shielding effect for DH fuels. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Group structure selection with random forests

Choosing an appropriate group structure for multigroup transport is far from an exact science. For some applications, one blindly uses a group structure developed years ago by forgotten methods. Furthermore, one sometimes uses the same group structure for a variety of problems, even if the group structure was originally developed with a certain application in mind. In this work, we create optimized group structures with simulated annealing for critical assembly test problems and apply a random forest regressor with bagging to choose the best group structure based on parameters of the different test problems. The optimized group structures were generated using a simulated annealing optimizer for several simple, spherical, and unreflected problems. The optimization was performed to minimize a cost function that included fission rate, absorption rate, leakage, and k{sub eff}. A random forest regressor was then trained on a set of International Criticality Safety Benchmark Evaluation Project inputs and used to select one of these six group structures. The trained machine learning model chose the best group structure 65% of the time, and one of the three best 89% of the time. Furthermore, it decreased the L2 error over all test problems by a factor of 25 when compared to the standard Los Alamos 70-group structure. In other words, the model chose group structure that were far more appropriate for the test problems than the traditional LANL group structure. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Setup and verification of a SCALE/KENO platform for generic FHR benchmark calculations

The work presented in this article is preliminary to downstream analysis of a generic fluoride salt-cooled high-temperature reactor (gFHR) core performed by the University of Tennessee in collaboration with Kairos Power (KP). A Monte Carlo transport model of the publicly available gFHR equilibrium core is developed in SCALE/KENO with multigroup energy treatment. Several output quantities of interest are used to verify the simulations against a benchmark model developed by KP using the continuous energy Monte Carlo code Serpent 2. Good agreement is seen in flux and fission rate profiles with a maximum relative difference of 1.4% and 2.8% respectively. Furthermore, an effective multiplication factor bias of 44 pcm was observed between the two simulations. The fuel temperature reactivity coefficient calculated with SCALE is within uncertainty to the reference model. This verification acts as a publicly reproducible benchmark for the gFHR in SCALE/KENO. A simplified depletion model is also presented where a single fuel pebble is depleted to discharge burnup through the equilibrium core while the equilibrium core is assumed to be invariant. This method produces results that intercept the equilibrium core concentrations in every case, however, an interesting artifact of this particular depletion model is uncovered. The phenomenon is shown to be a fundamental feature of the differential rate equations and inspires questions about how this system behaves when the time evolution of the equilibrium core is considered. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Evaluating the X2 initial core zero power physics tests with Serpent-Ants

The validation of the Ants nodal neutronics code for VVER applications is started by modelling the zero power physics tests for the initial core of the Khmelnitsky 2 nuclear power plant as described in the X2 benchmark using the Serpent-Ants two step neutronics calculation chain. The Ants prediction compare favorably against an earlier continuous energy Monte-Carlo reference solution as well as the measured data, except for discrepancies in the SCRAM worth between predicted and measured values. Similar discrepancies have been previously reported for VVER-1000 reactors, and the Ants predicted SCRAM worths match well with the Serpent results. Additionally, the assembly power distribution and axial power distribution predicted by Ants for the critical HZP state of the reactor are also compared to a reference Serpent solution with very good results. Lastly, an investigation into the effects of the few-group structure used in the nodal calculations and the use or lack of leakage correction for group constants on the results shows that the best accuracy is reached with seven or eight energy groups without leakage correction, although a reasonable accuracy can also be obtained with two or three energy groups using fundamental mode leakage correction. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Adjoint sensitivity analysis and data assimilation for verification of dry storage cask contents

Dry cask storage is a method for interim storage of spent fuel assemblies which contain fissile isotopes of uranium and plutonium. These can present a proliferation concern and consequently there is a need for non-destructive testing methods to verify a dry cask's contents for proliferation protection. We present an application of adjoint sensitivity analysis and data assimilation to a multigroup diffusion model of dry cask storage. Adjoint sensitivity analysis allows the efficient calculation of sensitivities for use in data assimilation to calibrate imprecisely known parameter values and data consistency tests to detect diversion scenarios. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Neutronics and thermal-hydraulics simulation of generic pebble-bed fluoride-salt-cooled high-temperature reactor (gFHR)

The fluoride-salt-cooled high-temperature reactor (FHR) is one type of the advanced reactors and has been attracting great interest from the research institutes and commercial companies in the recent years. However, currently the technology is relatively immature. To facilitate the design and safety analysis of FHRs, Kairos Power published the generic FHR (gFHR) benchmark. In this paper, a hybrid method combining the stochastic code and deterministic code is developed to simulate the gFHR benchmark. Serpent 2 is employed as a few-group cross section generator and the cross sections are applied to the finite difference neutron diffusion code AGREE. A good consistency with the gFHR benchmark is achieved. The agreement between the Serpent 2 results and the AGREE results shows that the hybrid method is applicable to FHRs. The thermal-hydraulics are coupled to neutronics for the steady-state calculations and good agreement between AGREE and SAM is achieved. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

POLCA8 - modelling of cross section variations inside hexagonal assemblies

This paper presents the POLCA8 approach for modelling non-constant cross section distributions inside hexagonal fuel assemblies. The multigroup diffusion equation is modified to account for intranodal cross section variations. The obtained equation is solved in a node-wise manner based on the Fourier expansion method. As a result of varying cross sections, the solution includes a particular part additionally to the homogeneous one. A method for obtaining the particular solution is derived. Numerical tests on a VVER-1000 core are presented showing the impact of cross-section variations to some key parameters for reactor operation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Genetic algorithm-based optimisation of the few-group structure for lead fast reactors analysis

The optimal choice of the few-group structure for full-core transient analyses is still an open issue in reactor physics, especially for fast system like the lead fast reactor. One possible approach to select the group boundaries is represented by heuristic search algorithms, such as evolutionary ones. In this paper, a genetic algorithm coupled with the SIMMER code is employed to determine optimized six-group boundaries for the analysis of the ALFRED reactor. The Serpent Monte Carlo code is adopted to produce both the fine-group cross section library and the fine-group flux, used as a figure of merit to drive the genetic optimisation. The results show that the algorithm is indeed able to find satisfactory solutions that comply with the set objectives and can be reasonably interpreted in light of the underlying physics of the considered core. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

A flexible linear diffusion acceleration to k-eigenvalue neutron transport with SN discontinuous finite element method

In this paper, we derive a flexible linear diffusion acceleration (LDA) for k-eigenvalue neutron transport discretized with discontinuous finite element method (DFEM) and discrete ordinates(SN). This LDA is based on our two pieces of previous works: the flexible non linear diffusion acceleration (NDA) for DFEM-SN and LDA for k-eigenvalue neutron transport using pre-conditioned Jacobian-free Newton-Krylov with self-adjoint angular flux (SAAF), continuous finite element method(CFEM), and SN. We point out the differences between LDA and NDA for DFEM-SN and the difference between DFEM-SN and SAAF-CFEM-SN for LDA. Numerical tests are presented to compare the convergence behaviour of NDA and LDA. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Investigation of the impact of difference between FRENDY and NJOY2016 on neutronics calculations

A nuclear data library is used as a starting input for all subsequent neutronics calculations. NJOY has been used worldwide as a nuclear data processing code to create cross section libraries for a long time. For the verification of NJOY method and for providing an alternative nuclear data processing tool, JAEA has been developed the new nuclear library processing code FRENDY. In this paper, nuclear calculations were performed using the ACE files and the multigroup libraries created by both FRENDY and NJOY, and the impacts on the neutronics characteristics due to nuclear data processing were investigated using those libraries. MCNP was used to compare the ACE files by calculating many benchmark problems including ICSBEP and it was confirmed that the k-eff values generally agree with each other within the range of statistical errors. The multigroup cross sections are verified by the BWR design codes LANCR/AETNA through calculation of a commercial BWR5 equilibrium core loaded with 9*9 fuels. It was confirmed that fuel assembly and core characteristics are consistent with each other. From the above investigations, it was confirmed that FRENDY can provide comparable continuous/multi-group neutron cross sections with NJOY. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

On the effect of scalar flux weighing of linearly anisotropic scattering matrices in few-group transport calculations

The majority of the codes available for homogenized group constant generation for deterministic transport calculations apply the approximation of scalar flux weighting during energy group condensation of higher-order anisotropic scattering matrices. In this paper, we discuss the effect of scalar flux weighting of the linearly anisotropic scattering matrices in the frame of S P{sub 3} and S{sub 12} calculations performed for a two-dimensional VVER-440 reactor benchmark. To compare group constants generated for two neutron energy groups, an infinite pin cell was homogenized with Serpent 2 and ERANOS ECCO. Serpent 2 applies scalar flux, while ERANOS ECCO performs current weighting of the linearly anisotropic scattering matrices during energy group condensation. For analyzing the effect of the various weighting options, three simple reactor models were built assuming different core sizes using standard rectangular assemblies with 15*15 fuel pins. Diffusion, SP{sub 3} and S{sub 12} calculations were performed for the 3 models using group constants generated with Serpent 2 and ERANOS ECCO. The effect of scalar flux weighting of linearly anisotropic scattering matrices is shown by comparing the decrease in reactivity due to the decreased reactor size, as well as the assembly power distribution to reference results obtained with Serpent 2 Monte Carlo calculations. Neglecting higher than linearly anisotropic scattering and indirect application of diffusion coefficients in higher-order transport calculations is advised if angular flux-moment spectra weighted higher-order scattering matrices cannot be generated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Development and benchmarking of transient nodal code SIMULATE5-K neutron kinetics solver

SIMULATE5-K is Studsvik's next generation best estimate transient code. The time dependent diffusion equation is solved with a nodal method consistent with that implemented in the licensed core design code SIMULATE5. Arbitrary number of neutron and delayed neutron precursor groups can be used. For the solution of the spatial problem, the coupling coefficients used to relate the node leakages are found by first converting the time dependent diffusion equation to a static diffusion equation with the use of flux and delayed neutron precursor dynamic frequencies. Once the static-like equations are obtained, the multi-group analytical nodal model is used to obtain the coupling coefficients, expressing the node leakage in terms of adjacent node average fluxes. The coupling coefficients are then inserted into the time dependent nodal balance equation. For the time integration, the time dependent neutron balance equation is solved with the frequency transformation method. The treatment of the temporal dependence yields a fixed source problem which can be solved utilizing the existing fixed-source methodology. The primary purpose of this paper is to describe the neutron kinetics methodology implemented in SIMULATE5-K. The accuracy of the method is demonstrated for a series of well-known, neutronic-only benchmark problems. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Thermal feedback coupling in a transient Monte Carlo high-order/low-order scheme

Transient simulations of nuclear systems face the computational challenge of resolving both space and time during reactivity changes. A common strategy for tackling this issue is to split the neutron flux into shape and amplitude functions. This split can be solved with high- order/low-order methods. While this multi-fidelity approach has traditionally been reserved for deterministic methods, it is also possible to implement in Monte Carlo as an efficient alternative to Dynamic Monte Carlo. This work implements the frequency transform method with thermal feedback in high-order/low-order Monte Carlo by blending static coupling methods such as single-batch Monte Carlo, with a simple thermal-fluids calculation. While previous work focused solely on prescribed transients, the addition of time-dependent thermal-fluids allows transients to be self-propagating. Tests were run in a fluids-initiated transient to showcase the basic functionalities of this methodology. Preliminary results behave as expected, paving the way for studying more sophisticated thermal-fluids coupling methods. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Transient MOC with frequency transform and DSA on unstructured mesh

We present an implementation of the transient method of characteristics (MOC) with isotropic time derivatives, accelerated by diffusion synthetic acceleration (DSA). The fully implicit frequency transform method is used to solve the transient problem with analytic precursor integration. The code works on meshes composed of almost any of the commonly used non-curvilinear finite element types, and can handle the deformation of geometry in time-dependent transport calculations. We present results of a continuous Fourier analysis for the transient multigroup DSA problem, and representative benchmarking results are presented for the C5G7-TD benchmark in 2D showing reasonable performance and agreement compared to other codes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Reduced-order modeling of neutron transport separated in energy by Minimax Proper Generalized Decomposition

In this article, we demonstrate a Petrov-Galerkin Proper Generalized Decomposition (PGD) known as Minimax PGD for modeling neutron transport separated in energy. To compare the Minimax with the classical Galerkin PGD, we assess both on a model problem of UO{sub 2} or Mixed Oxide (MOX) fuel pins in an infinite lattice with 3 industry-standard energy meshes. We find the Minimax PGD achieves a superior decomposition to Galerkin PGD, both with and without update of the energy modes. This suggests Minimax PGD may be more computationally efficient, provided this reduction in modes (to achieve a given accuracy) outweighs the cost of solving the necessary adjoint problems. In either case, we note that PGD offers an a priori Reduced-Order Model (ROM) which may be dramatically cheaper to solve than the full-order model, especially in problems with fine to ultrafine energy meshes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Low level coupling scheme between neutronics and thermal-hydraulics based on Anderson acceleration

The simulation of nuclear reactors is a multiphysics problem mixing, amongst other fields, neutron transport and thermal-hydraulics. The simplest and most used approach in multiphysics simulation is based on the coupling of single-physics codes in a black-box fashion. However, in order to reduce the computational time needed for such simulations, case-dependent optimizations are often required. In this paper, we aim at reducing the computational time required to solve a coupled neutronic/thermal-hydraulic steady-state problem on a simplified Pressurized Water Reactor (PWR) core. The idea is to deal simultaneously with the coupling of the energy groups of the deterministic neutronic description of the core and its thermal-hydraulic description with the Anderson acceleration. By doing so, the fission source terms are directly accelerated instead of the power map as done in most cases. The power method used to solve the k-eigenvalue problem inside the neutronic solver is thus accelerated with the Anderson acceleration. The numerical experimentations conducted in this work are performed using APOLLO3 and THEDI, and indicate that such coupling strategy improves the convergence rates in terms of number of iterations required and the total computational time. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Optimizing group structures using hierarchical division

Creating group structures with few groups that give low errors is a difficult problem in reactor analysis. In recent years, automated optimization techniques have been applied to this task. We continue this trend by applying the hierarchical division algorithm to generate optimized group structures that minimize a cost function. At each stage, the algorithm adds a single group boundary to an existing group structure, dividing one group into two to increase the resolution of the group structure. The location of the added boundary is the one that gives the lowest error over all possible new boundary locations. Our implementation requires a beginning group structure, a set of candidate new boundary locations, and a set of reference reaction rates. As a proof of concept, we used WIMS-69 as the initial group structure, XMAS-172 as the ending group structures, and a 344-group reference group structure. Testing on two simple, homogenized reactor problems, we found that hierarchical division was able to reduce the error by a factor of around 5 with an increase of only 15% in the number of groups. Because hierarchical division can get stuck in local minima, it often reaches a plateau in its error reduction capability as many groups are added. Nevertheless, we find hierarchical division has strong potential to make good group structures into great group structures at a modest increase in computational cost. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗