PRLO Code Interface & Workflow - 2025 Advances in Nuclear Fuel Management (ANFM) Workshop
This is a workshop presentation material, introducing the Plant ReLoad Optimization (PRLO) Interface and workflow.
Genetic Algorithm↗
SEARCH · Engineering Papers
Results for “M CODES”
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.
Type Ia Supernova Nucleosynthesis: Metallicity-dependent Yields
Abstract Type Ia supernova explosions (SN Ia) are fundamental sources of elements for the chemical evolution of galaxies. They efficiently produce intermediate-mass (with Z between 11 and 20) and iron group elements—for example, about 70% of the solar iron is expected to be made by SN Ia. In this work, we calculate complete abundance yields for 39 models of SN Ia explosions, based on three progenitors—a 1.4 M ⊙ deflagration detonation model, a 1.0 M ⊙ double detonation model, and a 0.8 M ⊙ double detonation model—and 13 metallicities, with 22 Ne mass fractions of 0, 1 × 10 −7 , 1 × 10 −6 , 1 × 10 −5 , 1 × 10 −4 , 1 × 10 −3 , 2 × 10 −3 , 5 × 10 −3 , 1 × 10 −2 , 1.4 × 10 −2 , 5 × 10 −2 , and 0.1, respectively. Nucleosynthesis calculations are done using the NuGrid suite of codes, using a consistent nuclear reaction network between the models. Complete tables with yields and production factors are provided online at Zenodo:Yields ( https://doi.org/10.5281/zenodo.8060323 ). We discuss the main properties of our yields in light of the present understanding of SN Ia nucleosynthesis, depending on different progenitor mass and composition. Finally, we compare our results with a number of relevant models from the literature.
79 ASTRONOMY AND ASTROPHYSICS↗
Foundations of plasma standards
Abstract The field of low-temperature plasmas (LTPs) excels by virtue of its broad intellectual diversity, interdisciplinarity and range of applications. This great diversity also challenges researchers in communicating the outcomes of their investigations, as common practices and expectations for reporting vary widely in the many disciplines that either fall under the LTP umbrella or interact closely with LTP topics. These challenges encompass comparing measurements made in different laboratories, exchanging and sharing computer models, enabling reproducibility in experiments and computations using traceable and transparent methods and data, establishing metrics for reliability, and in translating fundamental findings to practice. In this paper, we address these challenges from the perspective of LTP standards for measurements, diagnostics, computations, reporting and plasma sources. This discussion on standards, or recommended best practices, and in some cases suggestions for standards or best practices, has the goal of improving communication, reproducibility and transparency within the LTP field and fields allied with LTPs. This discussion also acknowledges that standards and best practices, either recommended or at some point enforced, are ultimately a matter of judgment. These standards and recommended practices should not limit innovation nor prevent research breakthroughs from having real-time impact. Ultimately, the goal of our research community is to advance the entire LTP field and the many applications it touches through a shared set of expectations.
data and reaction mechanisms↗
Generation and suppression of runaway electrons in MST tokamak plasmas
This paper discusses the behavior of runaway electrons in tokamak plasmas at low electron density, plasma current, and magnetic field using experimental data from the Madison Symmetric Torus and computational data from the NIMROD nonlinear resistive 3D MHD code. Density thresholds for the onset and suppression of runaway electrons are determined experimentally in steady tokamak plasmas, and in plasmas with a population of runaway electrons, resonant magnetic perturbations with different poloidal mode numbers are applied. Poloidal mode number m = 3 perturbations suppress the runaway electrons, while perturbations with m = 1 have little effect. This difference is consistent with the difference in computed magnetic topologies. The m = 1 RMP has minimal effect on the topology, while the m = 3 RMP produces a broad region of stochasticity, which can allow for rapid loss of the runaway electrons.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Real-Time Interactive 4D-STEM Phase-Contrast Imaging From Electron Event Representation Data: Less computation with the right representation
The arrival of direct electron detectors (DED) with high frame-rates in the field of scanning transmission electron microscopy has enabled many experimental techniques that require collection of a full diffraction pattern at each scan position, a field which is subsumed under the name four dimensional-scanning transmission electron microscopy (4D-STEM). DED frame rates approaching 100 kHz require data transmission rates and data storage capabilities that exceed commonly available computing infrastructure. Current commercial DEDs allow the user to make compromises in pixel bit depth, detector binning or windowing to reduce the per-frame file size and allow higher frame rates. This change in detector specifications requires decisions to be made before data acquisition that may reduce or lose information that could have been advantageous during data analysis. The 4D Camera, a DED with 87 kHz frame-rate developed at Lawrence Berkeley National Laboratory, reduces the raw data to a linear-index encoded electron event representation (EER). Here we show with experimental data from the 4D Camera that linear-index encoded EER and its direct use in 4D-STEM phase contrast imaging methods enables real-time, interactive phase-contrast from large-area 4D-STEM datasets. Furthermore, we detail the computational complexity advantages of the EER and the necessary computational steps to achieve real-time interactive ptychography and center-of-mass differential phase contrast using commonly available hardware accelerators.
4D-STEM↗
Hyper-fidelity depletion coupled with discrete pebble motion in pebble bed reactors
Pebble bed reactors have raised new interest during the past decade due to their attractive characteristics. Therefore, accurate simulations must be performed to better understand these systems and ensure optimal and safe designs. Most current methods use lower fidelity approaches with representative unit-cells or macro-zones with uniform fluxes, which have accuracy and flexibility limitations. A novel hyper-fidelity method for pebble bed reactors depletion is presented and internally couples Serpent 2 and a pseudo-motion routine. Pseudo-motion is applied handling vertical shifts of compositions in a static pebble bed, random reinsertion of used pebbles, insertion of fresh pebbles and used pebbles discarding. Associated with individual depletion to correctly determine the flux spectrum and composition in each pebble, this hyper-fidelity approach paves the way towards more accurate depletion calculation in pebble bed reactors. Using this method, a demonstration is completed on a small-scale reactor. In this application, the core reaches equilibrium, and the following data is extracted: core-wise parameters evolution, pebble-wise spatial and statistical distribution. Discarded pebbles are analyzed, and relevant information is shown. This work proves the feasibility of hyper fidelity depletion with Serpent 2, and the range of use for this method: reactor design and analysis for equilibrium and slow transients, lower fidelity methods validation and feeding fuel performance, thermal-hydraulics, or waste management models. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Subcritical Multiplication with a Fixed Source
In a subcritical, multiplying medium, the system multiplication describes the expected total number of neutrons created by a single source neutron. Subcriticality plays a large role in criticality safety and thus it is vital for the subcritical multiplication factor be accurate, especially as a system approaches criticality. This work examines the accuracy of calculating the system multiplication using the MCNP6.2 ® k-eigenvalue power iteration (KCODE) method when a fixed-point source is present in a multiplying medium, for near critical systems. This work compares the standard approach for calculating system multiplication, using the fixed-source calculational approach, to a new, single k-eigenvalue power iteration approach that incorporates a fixed-source component and a fission-source component into a single calculation. For the remainder of this paper, some theoretical background and numerical results for an approximate k eigenvalue approach, an accurate fixed-source approach and a new and more accurate k-eigenvalue approach to computing system multiplication are provided.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Advanced Simulation and Experiments of Strongly Coupled Geomechanics and Flow for Gas Hydrate Deposits: Validation and Field Application. Final Report
The objectives of the proposed research are (1) to investigate geomechanical responses induced by depressurization experimentally and numerically; (2) to enhance the current numerical simulation technology in order to simulate complex physically coupled processes by depressurization and (3) to perform in-depth numerical analyses of two selected potential production test sites: one based on the deposits observed at the Ulleung basin UBGH2-6 site; and the other based on well-characterized accumulations from the westend Prudhoe Bay. To these ends, the recipient will have the following specific objectives: 1). Information obtained from multi-scale experiments previously conducted at the recipient’s research partner (the Korean Institute of Geoscience and Mineral Resources (KIGAM)) that were designed to represent the most promising known Ulleung Basin gas hydrate deposit as drilled at site UBGH2-6 will be evaluated (Task 2). These findings will be further tested by new experimental studies at Lawrence Berkeley National Laboratory (LBNL) and Texas A&M (TAMU) (Task 3) that are designed capture complex coupled physical processes between flow and geomechanics, such as sand production, capillarity, and formation of secondary hydrates. The findings of Tasks 2 and 3 will be used to further improve numerical codes. 2) Develop (in Tasks 4 through 6) an advanced coupled geomechanics and non-isothermal flow simulator (T+MAM) to account for large deformation and strong capillarity. This new code will be validated using data from the literature, from previous work by the project team, and with the results of the proposed experimental studies. The developed simulator will be applied to both Ulleung Basin and Prudhoe Bay sites, effectively addressing complex geomechanical and petrophysical changes induced by depressurization (e.g., frost-heave, strong capillarity, cryo-suction, induced fracturing, and dynamic permeability).
02 PETROLEUM↗
Determination of the spent fuel decay heat with the VERA core simulator
VERA core simulator capabilities to predict spent fuel decay heat have been explored. The comparison with the Serpent2 and NEWT/TRITON sequence from the SCALE package, performed on the NPP Krško fuel test cases, has enabled independent verification of the successful implementation of the VERA depletion/decay calculation sequence applied on the extended nuclide set. Since the VERA is a 3-D core simulator, the focus in the rest of the paper was on the investigation of the 3-D effects and non-linearity of the soluble boron concentration averaging process, which are difficult to assess with most of available computational tools. A difference in the decay heat prediction between the 3-D and 2-D model approaches was analyzed. In addition, the impact of the Inconel grids on the decay heat was determined. Averaging the concentration of soluble boron over the fuel irradiation period has been found to be a viable approach, since it results in slight decay heat overprediction, which is considered conservative for most applications. Finally, demonstration on the Watts Bar Unit 1 cycle 1 depletion has shown formidable VERA capabilities to accurately predict spent fuel decay heat.
Kromar, Marjan↗
PBjam: A Python Package for Automating Asteroseismology of Solar-like Oscillators
Asteroseismology is an exceptional tool for studying stars using the properties of observed modes of oscillation. So far the process of performing an asteroseismic analysis of a star has remained somewhat esoteric and inaccessible to nonexperts. In this software paper we describe PBjam, an open-source Python package for analyzing the frequency spectra of solar-like oscillators in a simple but principled and automated way. The aim of PBjam is to provide a set of easy-to-use tools to extract information about the radial and quadropole oscillations in stars that oscillate like the Sun, which may then be used to infer bulk properties such as stellar mass, radius, age, or even structure. Asteroseismology and its data analysis methods are becoming increasingly important as space-based photometric observatories are producing a wealth of new data, allowing asteroseismology to be applied in a wide range of contexts such as exoplanet, stellar structure and evolution, and Galactic population studies.
79 ASTRONOMY AND ASTROPHYSICS↗
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↗
A quantal diffusion approach for multinucleon transfer in heavy-ion collisions
The stochastic mean-field (SMF) approach beyond the time-dependent-Hartree-Fock theory is used to explore the primary production cross sections in 64 Ni + 130 Te at the bombarding energy E c.m. = 184.3 MeV and 206 Pb + 118 Sn at E c.m. = 436.8 MeV. Secondary production cross-sections in the same systems are calculated using a statistical de-excitation model with GEMINI++ code. Furthermore, the obtained results are compared with available experimental data. Analysis employing SMF and GEMINI++ exhibit a good agreement with the experimental data.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
2.0 - MOOSE: Enabling massively parallel multiphysics simulation
The last 2 years have been a period of unprecedented growth for the MOOSE community and the software itself. The number of monthly visitors to the website has grown from just over 3,000 to now averaging 5,000. In addition, over 1,800 pull requests have been merged since the beginning of 2020, and the new discussions forum has averaged 600 unique visitors per month. The previous publication has been cited over 200 times since it was published 2 years ago. This paper serves as an update on some of the key additions and changes to the code and ecosystem over the last 2 years, as well as recognizing contributions from the community.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles
The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.
42 ENGINEERING↗
NWChem: Past, present, and future
Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Toward exascale whole-device modeling of fusion devices: Porting the GENE gyrokinetic microturbulence code to GPU
GENE solves the five-dimensional gyrokinetic equations to simulate the development and evolution of plasma microturbulence in magnetic fusion devices. The plasma model used is close to first principles and computationally very expensive to solve in the relevant physical regimes. In order to use the emerging computational capabilities to gain new physics insights, several new numerical and computational developments are required. Here, we focus on the fact that it is crucial to efficiently utilize GPUs (graphics processing units) that provide the vast majority of the computational power on such systems. In this paper, we describe the various porting approaches considered and given the constraints of the GENE code and its development model, justify the decisions made, and describe the path taken in porting GENE to GPUs. We introduce a novel library called gtensor that was developed along the way to support the process. Performance results are presented for the ported code, which in a single node of the Summit supercomputer achieves a speed-up of almost 15× compared to running on central processing unit (CPU) only. Typical GPU kernels are memory-bound, achieving about 90% of peak. Our analysis shows that there is still room for improvement if we can refactor/fuse kernels to achieve higher arithmetic intensity. We also performed a weak parallel scalability study, which shows that the code runs well on a massively parallel system, but communication costs start becoming a significant bottleneck.
Germaschewski, K. (ORCID:0000000284956354)↗
Nonlinear light-output calibration of the oxygenated xylene scintillators used in OMEGA neutron time-of-flight spectrometers
Neutron time-of-flight (nTOF) spectrometers are essential instruments for measuring and evaluating the performance of inertial confinement fusion implosions. The neutron spectrometers utilized for the OMEGA laser include two liquid-based scintillators, each consisting of a large volume filled with xylene that is coupled to four photomultiplier tubes. Analysis of the signal from these detectors requires detailed knowledge of the scintillator’s light output, which is needed to fit the nTOF spectrum, from which the neutron energy spectrum is informed. The light output is nonlinearly proportional to the neutron energy, which, in turn, affects the interpretation of the neutron energy spectrum from a TOF signal. A recent campaign on OMEGA was performed to calibrate the xylene detectors and infer the shape of the light-output curve. The campaign utilized materials with increasing Z placed in the OMEGA target chamber to initiate scattering events with the 14 MeV fusion neutrons. This process leads to the production of backscatter neutrons of varying energies that appear as peaks in the nTOF data. Simulations using a neutron transport code were combined with the measured deuterium–tritium neutron yields to calculate the expected backscattered neutron yields from the well-known scattering cross sections of each material. Furthermore, the neutron-energy dependent light output of the scintillator inferred from the experiment is compared to the light-output curve simulated with a neutron transport code for the following neutron energies: 1.5, 2.5, 6, and 14 MeV.
47 OTHER INSTRUMENTATION↗
MONKES: a fast neoclassical code for the evaluation of monoenergetic transport coefficients in stellarator plasmas
Abstract MONKES is a new neoclassical code for the evaluation of monoenergetic transport coefficients in stellarators. By means of a convergence study and benchmarks with other codes, it is shown that MONKES is accurate and efficient. The combination of spectral discretization in spatial and velocity coordinates with block sparsity allows MONKES to compute monoenergetic coefficients at low collisionality, in a single core, in approximately one minute. MONKES is sufficiently fast to be integrated into stellarator optimization codes for direct optimization of the bootstrap current and to be included in predictive transport suites. The code and data from this paper are available at https://github.com/JavierEscoto/MONKES/ .
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗