Search

Energy discretization of the transport equation is difficult due to numerous strong, narrow cross-section (XS) resonances. The standard traditional multigroup (MG) method can be sensitive to approximations in the weighting spectrum chosen for XS averaging, which can lead to inaccurate treatment of important phenomena such as self-shielding. We generalize the concept of a group to a discontiguous range of energies to create the Finite-Element with Discontiguous-Support (FEDS) method. FEDS uses clustering algorithms from machine learning to determine optimal definitions of discontiguous groups. By combining parts of multiple resonances into the same group, FEDS can accurately treat resonance behavior even when the number of groups is orders of magnitude smaller than the number of resonances. In this paper, we introduce the theory of the FEDS method and describe the workflow needed to use FEDS, noting that ordinary MG codes can use FEDS XSs without modification, provided these codes can handle upscattering. This allows existing MG codes to produce FEDS solutions. In the context of light water reactors, we investigate properties of FEDS XSs compared to MG XSs and compare k -eigenvalue and reaction rate quantities of interest to continuous-energy Monte Carlo, showing that FEDS provides higher accuracy and less cancellation of error than MG with expert-chosen group structures. [ABSTRACT FROM AUTHOR]

AbstractArbitrary Lagrangian-Eulerian methods are a popular choice for hydrodynamic modeling in radiation (rad-hydro) simulations. Because these methods involve a relaxation step that moves the mesh relative to material boundaries, multimaterial spatial zones are generally present. Accurate treatments of these zones are needed to resolve various physical phenomena of interest for inertial confinement fusion applications. However, these codes are often paired with single-material, deterministic thermal radiative transfer (TRT) codes that are oblivious to the material compositions of each zone. These single-material TRT codes can only accept homogenized material properties (opacities, specific heats, etc.) from the hydrodynamic code and output homogenized solutions. After each TRT time step, the multimaterial hydrodynamic code must dehomogenize the quantities computed by the TRT package in order to update subzonal material temperatures.The process by which hydrodynamic codes perform this dehomogenization has not been well documented in previous literature, and the methods can vary significantly from code to code. The purpose of this paper is to document, study, and compare existing techniques used for rad-hydro simulations as well as present a new method with potentially promising results. We summarize several methods and give comparisons on infinite-medium problems as well a finite-medium problem for two of the methods. [ABSTRACT FROM AUTHOR]

It is challenging to select an appropriate group structure for any given multigroup neutron transport problem. Many group structures were designed long ago, and the reasoning behind the creator's choices may be unknown. In this work, we apply the simulated annealing optimization method to develop improved group structures for a set of test problems. We then use a random forest (a machine learning method) to identify which group structure will be the best for any new problem based on input characteristics, such as geometry and isotopics. Simulated annealing spans a large solution space before narrowing in on an optimal solution, avoiding local minima by jumping around. Our solution space, however, is large and inconsistent, making finding the optimal group structure infeasible. Instead, we find potentially optimal group structures, ones that yield more accurate solutions than our standard group structures, but are probably not the "best" possible. Group structures are obtained for six classes of problems, ranging from a fast 233U system to a thermal 239Pu system. These were chosen to encompass a series of critical assemblies from the International Criticality Safety Benchmark Evaluation Project (ICSBEP) handbook. These optimized group structures were used in PARTISN for a large range of ICSBEP critical assemblies and compared to the traditional Los Alamos National Laboratory group structures. Our reference solution was from 618-group PARTISN runs. The results were used to train a random forest regressor model with bagging, which was then tested on similar benchmarks. The bagging regressor model chose the best group structure from 52% to 65% of the time, and a subjectively "good" group structure up to 91% of the time. [ABSTRACT FROM AUTHOR]