SEARCH · Engineering Papers
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.
EPRDATA in MCNP comes in 2 libraries, EPRDATA12 and EPRDATA14. The acronym EPR stands for Electron Photon, and atomic Relaxation. The libraries are in an extended form of the ACE format. This extended data includes more physics than the more common MCPLIB photo-atomic libraries, such as MCPLIB63 and MCPLIB84.
Quantifying the performance of a new evaluation is difficult and doing so for a new library is significantly more challenging. Cumulative χ 2 metric is an established approach to comparing distributions. In this case, C/E results from different libraries. Andrei Trkov and Roberto Capote have published cumulative χ 2 results for ENDF/B-VIII.0. Cumulative χ 2 metric was recently calculated for VALID library maintained at ORNL.
Researchers noted an unusually large bias for 8-in. nickel reflected HEU sphere in HMF-003 between the 252- group library and the CE library in SCALE 6.2.4. The bias was investigated by reviewing reactions that $k_{eff}$ is sensitive to using TSUNAMI and collecting reaction rate data using tools within SCALE. For this system, bias is primarily due to the elastic scattering in nickel. Researchers compared new multigroup structures in SCALE 6.3. For this system, reduction in CE-to-MG bias seen in the 302-group and further improved in the 1597-group structure. When researchers compared ENDF/B-VIII.0 in SCALE 6.3, library showed improved results with new nickel evaluation. However, the MG-to-CE bias can still be large as large biases can occur in any system. This example highlights the importance of validating results with measured systems.
The computational bias of criticality safety computer codes must be established through the validation of the codes to critical experiments. A large collection of suitable experiments has been vetted by the International Criticality Safety Benchmark Evaluation Project (ICSBEP) and made available in the International Handbook of Evaluated Criticality Safety Benchmark Experiments (ICSBEP Handbook). More than 600 cases from this handbook have been prepared and reviewed within the Verified, Archived Library of Inputs and Data (VALID), which is maintained by the Reactor and Nuclear Systems Division at Oak Ridge National Laboratory. The performance of the KENO V.a and KENO-VI Monte Carlo codes within the SCALE 6.2.4 code system is assessed using the VALID models of benchmark experiments. A range of nuclear cross section libraries based on Evaluated Nuclear Data File (ENDF)/B-VII.1 in both multigroup (MG) and continuous energy (CE) formats is considered. The critical experiments available to validate the KENO V.a code cover 15 broad categories of systems. These systems use a range of fissile materials, including a range of uranium enrichments, various plutonium isotopic vectors, and some mixed uranium/plutonium oxides. The physical forms of the fissile material also vary and are represented as metal, solutions, or arrays of rods or plates in a water moderator. The neutron energy spectra of the systems also vary and cover fast, intermediate, mixed, and thermal spectra. Over 550 of the total cases use the KENO V.a code for the four nuclear data libraries considered in this report.
This document pertains to the reporting requirements of DOE contract FE-0031767. The document covers the development of a reduced order model library that can be used to represent coal-fired power plants. The models are transient and physics-based. Model classes in the library that correspond to different components in CFPPs are presented with physical explanations and representative results. Also, the use of the library to build an overall model for a representative CFPP is presented along with simulation results.
This presentation highlights the process of developing a thermal data library for MCNP6 to support on-the-fly S(alpha, beta) sampling for temperature ranges applicable to nuclear criticality safety. Advanced on-the-fly (OTF) strategy driven libraries have been developed for six materials based on ENDF/B-VIII.0. Validation of the OTF libraries is being conducted currently.
In general, excitonic and plasmonic nanoscale materials in close proximity show high potential for significant breakthroughs in energy related materials research. The interactions between these two kinds of materials result in coupled optical transitions (plexcitons), distinct from those of both the individual exciton and plasmon as well as from those of the sum of their constituents (synergistic effects). By linking together materials-research and physical-research approaches, this project contributes to a concerted approach on nanomaterials energy research. The project’s overall goal is to accelerate the development of well-defined plexcitonic model systems consisting of carefully engineered plasmonic and excitonic nanomaterial— essential for both gaining a fundamental understanding of plexcitonic nanomaterials and the development of novel design principles for biomimetic solar energy harvesting. During the 3-year project period and the terminal renewal with limited support for a 12-month period, we successfully synthesized and characterized (1) a robust excitonic nanomaterial and (2) a library of plasmonic nanoparticles as well as developed (3) a microfluidic platform for homogenous nanosynthesis as summarized below: (1) Robust Excitonic Nanomaterial. Supramolecular assemblies are Nature’s most successful material system for solar energy harvesting. However, photovoltaic devices based on artificial supramolecular assemblies continue to be stymied by disappointing efficiencies and poor stability. The conceptual failure may lie in current solar cell architectures, which rely on solidifying supramolecular assemblies as an ensemble into a solid matrix, neglecting the intrinsic fragility of the assemblies’ internal structure, thus disrupting or even destroying their delicate optoelectronic properties, that is, delicate Frenkel excitonic properties. Supramolecular assemblies may finally serve as usable light harvesting material systems for solar energy conversion technologies, only if they meet the following criteria: (a) Stability, that is, the fragile structure including its delicate Frenkel excitonic character needs to be stable, (b) Robustness, that is, resistant against elevated and fluctuating temperatures, and (c) Viability for device integration, that is, capable of being immobilized onto solid substrates. Here, by developing a nanocomposite via a tunable, cage-like scaffold design, we successfully provided stable supramolecular nanocomposites, that inhabit robust Frenkel excitons despite harming environmental conditions such as extreme heat stress. (2) Library of Plasmonic Nanoparticles. Naturally, current models describing plasmonic hybrid quantum states—plasmonic hybridizations—parallel those developed for molecular orbitals, equating individual plasmonic nanostructures with “atoms” and the plasmonic nanoassemblies with “molecules.” In analogy to organic synthesis, a suitably robust fabrication method would allow for “atom-like” manipulation of “molecule-like” plasmonic nanoassemblies; of high value for next-generation energy nanotechnologies. Despite this frequent comparison, current plasmonic nanoassembly fabrication methods favor top-down templating over wet-chemical synthesis, however, achieving precise control over nanostructure’s geometry and surface characteristics remain an art and a scientific challenge. The conceptual failure may lie in the current wet-chemical synthesis paradigm, as it relies on the accessibility of a multi-dimensional synthesis parameter space through limited, rather one-dimensional synthesis procedures by employing step-by-step approaches. Solution-based nanoarchitectonics for rational design of precisely built plasmonic nanoassemblies via solution-based fabrication may finally be possible only if multi-dimensional syntheses approaches are available that allow for comprehensive control over the plasmonic nanomaterials’ (a) Structural Properties and (b) Surface Properties. Here, by developing an innovative multidimensional 1,3-propanediol based polyol synthesis, we successfully provided control over the plasmonic building-block’s geometry (size and shape) together with its surface characteristics. Our results present a critical step toward the vision of a “periodic table-like” system for plasmonic materials based on straightforward wet-chemical syntheses for energy nanotechnologies. Developing deliberate modifications on this synthesis, we generated a library of plasmonic nanostructures covering the vast parameter space—opening the door for fundamental investigation of plexcitonic model systems. (3) Microfluidic Platform for Homogenous Nanosynthesis. Control over structural properties of plexcitonic nanocomposites remains a challenge due to current limitations in nanosynthesis techniques. Slight variations in nanostructure’s geometry impact their optoelectronic properties, demanding precise synthesis beyond the capabilities of solution-based (batch) synthesis processes. In contrast, the small, confined liquid volumes used in microfluidics—a reaction technique where the manipulation of fluids takes place in channels with dimensions of tens of micrometers—allows for homogenous synthesis conditions, providing excellent control of the reaction and, as a result, of the materials’ geopmetry and composition. However, thus far, the majority of plexcitonic systems has been developed via batch synthesis. Here, by successfully developing a two-channel microreactor, our microfluidic-supported synthesis approach combines the advantages of both microfluidics and batch platforms, allowing for precise spatio-temporal control over all synthesis parameters opening the possibility for homogenous nanosythnesis of well-defined plexcitonic model systems.
The EOSPAC utility package is a collection of interface routines, which can be used to access the SESAME data library and perform various data adjustments and interpolations on the SESAME data. The SESAME data library contains both thermodynamic (e.g., equation of state) and transport coefficients (e.g., opacity and conductivity). Note, for simplicity, the term EOS (equation of state) used herein includes both thermodynamic variables and transport coefficients. The EOSPAC utility package is designed to be used by physics codes (henceforth ”host codes”) written in multiple languages and on multiple platforms. The remainder of this manual is organized into several sections. Chapter 2 discusses conventions such as data organization and routine names. Chapter 3 provides a general overview of basic theory and models implemented within EOSPAC. Chapter 4 provides a general overview of how to use the EOSPAC interface library. Chapters 5 to 7 describe the public interfaces of EOSPAC in detail. Chapter 8 provides a brief introduction to some related tools, which may be of use to the user. Chapter 9 provides details related to some selected numerical features of EOSPAC. Chapter 10 gives examples for using the interface routines described in chapters 5 to 7. Chapter 11 provides technical support contact information. Chapter 12 contains a brief set of acknowledgments. Chapter 13 contains a list of referenced documents. Finally, chapter 14 lists the “table types: mnemonic conventions”, “table types: grouped by category, sorted by name”, “table types: eospac version 5 cross reference”, “options: setup phase”, “data information parameters”, “meta-data information parameters”, “options: interpolation phase”, and the “error codes”.
The Intentional Forensics (IF) project seeks to design a suite of tagging materials that can be used to unobtrusively “barcode” nuclear fuel. The idea being that, if the fuel leaves institutional controls and is later recovered, the “barcode” can be used to determine the provenance of the fuel. The IF project is investigating several schemes for this “barcode” and taggants with tailored isotopic compositions are once leading concept. Recently 18 O was suggested as a potential taggant [IF2023]. This report assesses the quality of 16,17,18 O evaluations in the ENDF/B-VIII.0 library [ENDF8.0]. In a neutron-rich environment such as a nuclear reactor, we are interested in two main processes – neutron scattering and neutron absorption (or capture). For the typical neutron energies encountered in a reactor, absorption is determined exclusively by the neutron radiative capture cross section (in which a neutron is absorbed, and the residual nucleus emits one or more gamma rays). The capture cross sections are described by and can be completely determined from R-matrix parameters in the ENDF evaluations. Therefore, it is sufficient to consider the quality of the neutron resonance part of these evaluations. Neutron scattering is determined mainly by the elastic scattering cross section which is very similar for all oxygen isotopes and is essentially constant until around 100 keV. The cross section varies from 3-4 barns, depending on the isotope, and 16 O has the largest elastic cross section of all isotopes. Therefore, we expect that significant amounts of 17 O and 18 O could diminish the moderating role of oxygen in fuel. We note that, in traditional power reactor, the hydrogen in the coolant water provides the vast majority of the moderation in the reactor, so the impact of 17 O and 18 O in any taggant would be minimal. At higher energies, neutron resonances start to play a role and the R-matrix parameters determine the cross sections. Fig. 1 shows the mass region near stable oxygen isotopes. As nuclei capture neutrons, they transmute into a different isotope, one unit to the right. In particular, 16 O transmutes to 17 O, 17 O transmutes to 18 O and 18 O transmutes to 19 O. 19 O, being unstable, beta decays in 27 s to 19 F. So, even if the 18 O neutron capture data is good, the fact that it transmutes into 19 F may be problematic and needs investigation. Fig. 2 from [NuDat] shows the thermal neutron capture cross sections for all nuclei in the ENDF/B-VII.1 library [ENDF7.1]. Clearly oxygen thermal capture values are quite small compared to other nuclei. Note, there is no 18 O evaluation in the ENDF/B-VII.1 library. In any event, given that neutron capture events lower the reactor reactivity and potentially induce chemical changes in the fuel, we will focus our attention on the capture part of the neutron resonances.
This presentation touches on motivation and applications of Zirconium. Performance of evaluated libraries on double and quasi differential experiments is considered. Additionally, the performance of evaluated libraries in select ICSBEP benchmarks is considered. Provided is a detailed look in the performance of evaluated libraries in preliminary benchmark configurations for IER-516: Zirconium Test Assembly (ZTA). Conclusions and recommended experiments to satisfy Zirconium scattering nuclear data needs are presented.
The Multiphysics Object Oriented Simulation Environment (MOOSE) is a general finite element package meant for high performance solution of multiphysics problems in science and engineering. In this report, we present additions and enhancements to MOOSE funded by the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program. NEAMS-funded MOOSE framework library improvements include: expansion of support for multi-level multi-application restart; enhancement of coupling between native MOOSE applications and external libraries; addition of a sparse automatic differentiation (AD) container enabling non-local degree of freedom coupling; block-specific quadrature rules; further development of an eigenvalue executioner; faster setup of periodic boundary conditions; and creation of a mesh meta-data system streamlining simulation startup. Besides developing the framework library, NEAMS funds were used to overhaul a swath of important MOOSE infrastructure in order to substantially improve user experience. These critical infrastructure changes included adapting MOOSE to python 3, improving the test harness to catch race conditions, and most importantly transitioning MOOSE to use Conda, a globally known packaging system, that greatly eases software adoption.
Batched sparse linear algebra operations in general, and solvers in particular, have become the major algorithmic development activity and foremost performance engineering effort in the numerical software libraries work on modern hardware with accelerators such as GPUs. Many applications, ECP and non-ECP alike, require simultaneous solutions of many small linear systems of equations that are structurally sparse in one form or another. In order to move towards high hardware utilization levels, it is important to provide these applications with appropriate interface designs to be both functionally efficient and performance portable and give full access to the appropriate batched sparse solvers running on modern hardware accelerators prevalent across DOE supercomputing sites since the inception of ECP. To this end, we present here a summary of recent advances on the interface designs in use by HPC software libraries supporting batched sparse linear algebra and the development of sparse batched kernel codes for solvers and preconditioners. We also address the potential interoperability opportunities to keep the corresponding software portable between the major hardware accelerators from AMD, Intel, and NVIDIA, while maintaining the appropriate disclosure levels conforming to the active NDA agreements. The presented interface specifications include a mix of batched band, sparse iterative, and sparse direct solvers with their accompanying functionality that is already required by the application codes or we anticipated to be needed in the near future. This report summarizes progress in Kokkos Kernels and the xSDK libraries MAGMA, Ginkgo, hypre, PETSc, and SuperLU.
SOURCES-4C has some significant deficiencies that are mostly due to problems with the default data libraries that it uses. Some research groups have modified these libraries to fix most of the problems, but not all of them. Adding modern data libraries to SOURCES-4C is one of the most promising approaches to resolving these issues. Other promising options include using either NEDIS or a GEANT4 tool. Both of these codes approximately match the accuracy of the modified SOURCES-4C code. The NEDIS code is Russian. The GEANT4 tool might be computationally slow and it might also no longer be actively maintained. There are a handful of other (α,n) codes, but they are not accurate.
The Framework for Optimization of Resources and Economics (FORCE) tool suite is the U.S. Department of Energy’s Nuclear Integrated Energy Systems (IES) Program flagship tool suite for technoeconomic IES analysis of IES. This tool suite is useful for analysis designed to evaluate and improve the technoeconomics of energy production systems, particularly for systems including nuclear technology. In this report, we document the development activity for the FORCE tool suite to extend its capabilities as performed during fiscal year 2024. In addition to reliability and accessibility, capability is one of the three standards guiding the development of the FORCE tool suite and the software codes that are its constituent parts. Extending the capabilities of the FORCE tool suite allows analysis both within the IES program as well as industry, university, and laboratory partners to perform analysis with more accuracy, insight, and impactful narrative. Four areas of capability development were the focus of activity this year: economic parameter uncertainty quantification, multiresolution analysis, components-to-optimization workflow automation, and statespace construction workflows for real-time optimal control. In economic parameter uncertainty quantification, the ability of HERON to capture risk due to scenarios (weather and energy demand uncertainty) was expanded to also include uncertainties in financial parameters such as capital cost or operation and maintenance costs. By including these sources of uncertainty, which are sometimes very large compared with scenario uncertainty, HERON is better able to capture the risk posed by investment in various IES technology. Because of this, analysts can also consider the reduction in risks that can be realized by choice of some technologies. In multiresolution analysis, development activity extended on work completed previously. In fiscal year 2023, methods for decomposing time series signals, such as demand, solar and wind availability, and price profiles, were analyzed and down-selected to those most effective at splitting signals into different resolutions. These resolutions allow considering the influence of different energy demand and supply behaviors across different time scales. For example, energy demand might be divided into seasonal, weekly, and hourly profiles. In fiscal year 2024, this preliminary work was extended and implemented within the Risk Analysis Virtual Environment (RAVEN) risk and uncertainty analysis platform, which is used throughout the FORCE framework. This development of the “multi-resolution time series analysis” (MR-TSA) module in RAVEN allows training synthetic history generators on complex time series. These synthetic history generators can then be used in HERON for generating scenarios that represent possible market and weather scenarios that can be analyzed on different time scales. We envision completing this work in the future, implementing multiresolution dispatch optimization strategies that can make the most beneficial use of these stratified time histories. In components-to-optimization workflow development, workflows for translating user inputs of components into algorithms for algebraic optimization were selected and implemented. Similar algorithms within the Holistic Energy Resource Optimization Network (HERON) were separated from the main code base of HERON and gathered with the components-to-optimization workflows in the new Dispatch Optimization Variable Engine (DOVE) software library. This modularization allows FORCE users to analyze dispatch optimization and energy system duty cycles independently of HERON, which previously was a burdensome task. Additionally, these dispatch optimization algorithms, set up in an independent library, can now be used across all software applications within FORCE, especially including the real-time optimal control software Optimization of Real-time Capacity Allocation (ORCA). Allowing FORCE software to share dispatch optimization algorithms within a single library allows for improved software maintenance and reliability. In statespace characterization workflow development, alternative workflows for optimizing dispatch with additional technical accuracy was the focus, particularly to improve the real-time optimization decision making in ORCA. Using algorithms and workflows initially developed for the Feasible Actuator Range Modifier (FARM), workflows for determining the statespace representation of IES were identified and demonstrated. The resulting dispatch optimization required a more robust optimization algorithm than that originally used in HERON (and moved to DOVE), which required adding an alternate workflow to DOVE that can more accurately match the behavior of physical systems using a partial differential equation representation. In conclusion, capability developments in the FORCE tool suite in fiscal year 2024 have improved the ability of the FORCE tool suite to perform
There were many improvements and additions to ENDF/B-VIII.1. We tested four key moderator materials: Be, HF, FLiBe, and graphite. Be-metal moderated critical assemblies of HEU only see 26 pcm difference. Sd library adds new physics to the TSL which improves the agreement with the benchmark. Other changes besides the TSL impact the benchmark: 235U and 9Be. Free gas change affects library heavily, hundreds of pcm in wrong way – need to look at FG libraries. F and Cl TSLs added to 8.1. A 3k pcm disagreement was seen without new TSL. Adding 8.1 + TSL reduces pcm disagreement A LOT.
The R -matrix analysis of neutron-induced reactions for two stable chlorine isotopes ( 35 Cl and 37 Cl) was performed in the energy range of thermal up to 1.2 MeV. Starting from the repository of the ENDF/B-VIII.0 library and following recent measurement series, this work represents a significant improvement, particularly in the evaluation of the ( n , p ) reaction channel. The evaluation methodology used the R -matrix code SAMMY to generate a set of Reich–Moore resonance parameters for both stable chlorine isotopes. Consistent with recently measured 35 Cl( n , p ) data, the presented evaluation features a dramatic increase in the magnitude of the ( n , p ) reaction channel over the ENDF/B-VIII.1 nuclear data library and previous ENDF/B libraries.
The priority for AMPX development for FY2024 has been to support processing the ENDF/B-VIII.1 library. This includes the development to support Thermal Neutron Scattering Library, GNDS-2.1 support and AMPX maintenance, and transport libraries for SCALE 7.0.0.
Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.