Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code Assessment”

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 1,369 records · Page 76

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Dataset for "Machine Learning Ensembles Can Enhance Hydrologic Predictions and Uncertainty Quantification" Willard et al. (2025).

This data release provides all data and code used in the paper " "Machine Learning Ensembles Can Enhance Hydrologic Predictions and Uncertainty Quantifications" Willard et al. (2025)" to model stream temperature, evaluate, and assess results. The associated manuscript explores the effect of different ensemble construction techniques across different common machine learning (ML) architectures for predictions in unmonitored basins. Modeling was done using long short-term memory (LSTM), gated recurrent unit (GRU), temporal convolution network (TCN), and extreme gradient boosting (XGBoost) models, and stream site coverage spans 1362 locations across the conterminous United States. The ensemble construction techniques investigated include ensemble by random weight initialization, differing hyperparameters, different random subsets of training data, different subselections of input features, different architectures, and Monte Carlo Dropout. The data is organized into these items items:Code repository and data for the paper " "Machine Learning Ensembles Can Enhance Hydrologic Predictions and Uncertainty Quantifications" Willard et al. (2025).Code: stream_temp_ml_regionalization.zip contains the code repositoryData to run the code:- data_dir.zip -- contains all files that should be moved to the "DATA_DIR" variable defined in the "set_env_vars.sh" script in the code repository- metadata_dir.zip -- contains all files that should be moved to the "METADATA_DIR" variable defined in the "set_env_vars.sh" script in the code repositoryData produced by the code and used in the paper:- outputs_dir.zip - contains model output and results (outputs_dir/results), model weights (outputs_dir/models), and all other outputs used for the paper including feature importances.To cite this code, please use the following BibTeX or MLA entries:bibtex:@misc{willard2025streamensembles,author = {Jared Willard and Charuleka Varadharajan},title = {Dataset for "Machine Learning Ensembles Can Enhance Hydrologic Predictions and Uncertainty Quantification"},year = {2024},doi = {10.15485/2527393},publisher = {ESS-DIVE Repository},url = {https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2527393}}MLA: Willard, Jared, et al. Dataset for "Machine Learning Ensembles Can Enhance Hydrologic Predictions and Uncertainty Quantification". 2025. ESS-DIVE Repository, doi:10.15485/2448016.

54 ENVIRONMENTAL SCIENCES↗

Stochastic Analysis of Levee Stability Subject to Variable Seepage Conditions

© ASCE. Levee stability is highly influenced by seepage. Specifically, hydraulic conductivity distribution within a levee influences pore pressure distribution and controls the drained strength of the soil. In this study the influence of hydraulic conductivity and blanket layer thickness on failure probability is evaluated within the context of reliability analyses that also include soil strength and unit weight as random variables. First-Order Reliability Method (FORM) is used to evaluate reliability, rank random variables by importance and to obtain sensitivity of the solution to each random variable and its distribution parameters. Stability is computed using Spencer's method of slices coupled to a finite element seepage code to directly evaluate pore pressure. In addition, response surface method solutions are compared to the direct reliability solution to assess response surface accuracy. The results show that blanket layer thickness is more important than most strength parameters and that the uncertainty in hydraulic conductivity is less important to the variance in safety factor for the case of a low-permeability blanket layer. Numerical challenges caused by implementing a finite element limit-state function are discussed. Response surface methods are found to give a reasonable approximation to the direct reliability solution when the design point is between response surface fitting points.

Lanzafame, Robert↗

3-center and 4-center 2-particle Gaussian AO integrals on modern accelerated processors

We report an implementation of the McMurchie–Davidson (MD) algorithm for 3-center and 4-center 2-particle integrals over Gaussian atomic orbitals (AOs) with low and high angular momenta l and varying degrees of contraction for graphical processing units (GPUs). This work builds upon our recent implementation of a matrix form of the MD algorithm that is efficient for GPU evaluation of 4-center 2-particle integrals over Gaussian AOs of high angular momenta (l ≥ 4) [A. Asadchev and E. F. Valeev, J. Phys. Chem. A 127, 10889–10895 (2023)]. The use of unconventional data layouts and three variants of the MD algorithm allow for the evaluation of integrals with double precision and sustained performance between 25% and 70% of the theoretical hardware peak. Performance assessment includes integrals over AOs with l ≤ 6 (a higher l is supported). Preliminary implementation of the Hartree–Fock exchange operator is presented and assessed for computations with up to a quadruple-zeta basis and more than 20 000 AOs. The corresponding C++ code is part of the experimental open-source LibintX library available at https://github.com/ValeevGroup/libintx.

Chemistry↗

Fuel Fabrication Specification Impact Analysis for NBSR LEU Conversion

As part of a national initiative to enhance nuclear security and reduce proliferation risks, significant efforts have been undertaken by the National Nuclear Security Administration Material Management and Minimization Office of Reactor Conversion Program to convert U.S. high performance research reactors (USHPRRs) from the use of highly enriched uranium (HEU) to low-enriched uranium (LEU), including the National Bureau of Standards Reactor (NBSR). The current plan is to procure LEU fuel assemblies from commercial fabricators according to fuel specifications tailored for each USHPRR. The analysis conducted at Brookhaven National Laboratory was part of an effort to identify the sources of uncertainty in the fuel specifications that may impact the performance of the NBSR core after its conversion and, in particular, to assess the range of acceptable tolerance limits from the perspective of core safety and reactor performance. Using the stochastic neutronics code MCNP 6.2, the variations in important NBSR neutronics characteristics were analyzed as a function of the specification parameters independently and in combination. The important NBSR specification parameters analyzed were the fuel isotopic composition, the amount of impurity content in cladding, the fuel plate thickness, and the fuel element 235U mass loading. The range of variation of each specification parameter was based on the technical specification limit or available as-fabricated assay data and uncertainties. The NBSR neutronics characteristics selected for analysis were the reactor reactivity characteristics at equilibrium and the equilibrium fuel cycle length. Results show that with variations in the fabrication parameters of the as-fabricated U-10Mo fuel within the specification limitations, the excess reactivity of the NBSR LEU core remains well below the 15% Δk/k technical specification limit, and the shutdown margin is always significantly greater than the required 0.68% Δk/k. This ensures that the NBSR can be operated safely and reliably shut down for all analyzed cases within the specified fabrication limits after the LEU conversion. In the prototypic case, the fuel cycle length was 1.5 days longer than the targeted 38.5 days. In a credible worst-case scenario, where all low-reactivity parameters were combined, the fuel cycle length was reduced to 35.5 days, which is still considered manageable for reactor operations. Variations in cycle length are primarily driven by changes in 235U loading, with other parameters having secondary effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Overview of NSTX Upgrade initial results and modelling highlights

The National Spherical Torus Experiment (NSTX) has undergone a major upgrade, and the NSTX Upgrade (NSTX-U) Project was completed in the summer of 2015. NSTX-U first plasma was subsequently achieved, diagnostic and control systems have been commissioned, the H-mode accessed, magnetic error fields identified and mitigated, and the first physics research campaign carried out. During ten run weeks of operation, NSTX-U surpassed NSTX record pulse-durations and toroidal fields (TF), and high-performance ~1 MA H-mode plasmas comparable to the best of NSTX have been sustained near and slightly above the n = 1 no-wall stability limit and with H-mode confinement multiplier H 98y,2 above 1. Transport and turbulence studies in L-mode plasmas have identified the coexistence of at least two ion-gyro-scale turbulent micro-instabilities near the same radial location but propagating in opposite (i.e. ion and electron diamagnetic) directions. These modes have the characteristics of ion-temperature gradient and micro-tearing modes, respectively, and the role of these modes in contributing to thermal transport is under active investigation. The new second more tangential neutral beam injection was observed to significantly modify the stability of two types of Alfven eigenmodes. Improvements in offline disruption forecasting were made in the areas of identification of rotating MHD modes and other macroscopic instabilities using the disruption event characterization and forecasting code. Finally, the materials analysis and particle probe was utilized on NSTX-U for the first time and enabled assessments of the correlation between boronized wall conditions and plasma performance. These and other highlights from the first run campaign of NSTX-U are described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Verification of the prompt radiation shielding design of the ESS Accelerator (Stage 2)

The purpose of this CRADA is to perform additional analyses using MARS Monte-Carlo calculations code for verification of the prompt radiation shielding design of the ESS Accelerator. Stage 2 includes the following activities: 1.Activation assessment of the soil surrounding the wall of the ESS accelerator that maymigrate with the groundwater 2.Shielding assessment of the HVAC exhaust duct The existing ESS accelerator model (with possible minor modifications) shall be used in calculations. This work builds on the analysis originally performed under agreement CRADA FRA-2015-0302.

43 PARTICLE ACCELERATORS↗

Analysis of Space-Conditioning Loads in Commercial Buildings

Space conditioning end-uses, which include heating, cooling, and ventilation, represent a significant fraction of commercial building energy use, with a wide variety of heating and cooling technology options available in the market. In the interest of improving the overall efficiency of heating, ventilation and air-conditioning (HVAC) technologies, governments, utilities and private sector entities have implemented a variety of market transformation policies that aim to influence consumer purchase decisions. To evaluate the costs and benefits of such programs, analysts typically postulate a hypothetical default equipment choice, and compare it to one that provides comparable service with lower energy and/or power use. The corresponding reduced operating cost provides a benefit that offsets the potential higher cost of improved efficiency. Typically, life-cycle cost or cash-flow analyses are used to quantify the net economic benefit. These analyses require the capability to assess how a given equipment design would perform across a broad range of characteristics, both of the building and of the local weather. While these assessments can be performed using customized building simulations, it is generally not practical to develop and validate detailed building simulation code to cover all the potential variations of equipment design and installation. An alternative, and somewhat simpler, approach is to solely use detailed building simulations to generate time series of heating and cooling loads in commercial buildings. These loads can then be used as input to more detailed, stand-alone engineering models that simulate HVAC system performance under different equipment designs. This approach was used to evaluate a range of high-efficiency commercial packaged air conditioner design options for the Department of Energy’s Appliance and Equipment Standards Program (DOE-EERE 2015). While there may be some loss of precision relative to full simulation, the accuracy of this approach is sufficient for practical applications of cost-benefit analysis. This report describes the development of a database of commercial building heating and cooling loads, generated using the EnergyPlus software package, a whole building energy use model supported by the Department of Energy (DOE-EERE 2020a). EnergyPlus takes as input a set of configuration files that describe the building itself (size, zoning, envelope characteristics, etc.) and the various systems within it (HVAC, lighting, water heating, etc.). This analysis uses a publicly available collection of commercial reference buildings (CRB), comprised of sixteen building types and three vintages (DOE-EERE 2020b; Deru et al 2011). Each building is simulated in eighteen different locations, covering a wide range of climatic conditions. The prototype building description files assign the type of HVAC equipment used, and capacities across climate zones.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

MELCOR Code Change History (Revision 14959 to 18019)

This document summarily provides brief descriptions of the MELCOR code enhancement made between code revision number 14959and 18019. Revision 14959 represents the previous official code release; therefore, the modeling features described within this document are provided to assist users that update to the newest official MELCOR code release, 18019. Along with the newly updated MELCOR Users Guide and Reference Manual, users are aware and able to assess the new capabilities for their modeling and analysis applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Consistent Evaluation of the Prompt-fission Neutron Spectrum and Multiplicity for n+ 235,238 U and n+ 239 Pu

This report was written to satisfy a FY20 NCSP milestone on 235,238 U and 239 Pu. The milestone requires to “finalize a report assessing our methodology to evaluate prompt-fission neutron spectrum (PFNS) and multiplicity consistently”. More specifically, we study whether the code CGMF can reproduce ENDF/B-VIII.0 evaluated PFNS and average prompt-fission neutron multiplicities, $\bar{v}$, for 235,238 U and 239 Pu using one joint parameter set per isotope. If CGMF is shown to be able to reasonably reproduce ENDF/B-VIII.0 within its model-parameter space, this code could be used for future consistent evaluations of PFNS and $\bar{v}$. To answer this question, we explore here the parameter space of CGMF and its impact on calculated values and whether they are close to evaluated and experimental data. We also list experimental data that would enter a future evaluation and statistics method that could be used to obtain evaluated data and covariances. We conclude that values of $\bar{v}$ calculated by CGMF are reasonably close to ENDF/B-VIII.0 data, while more work on modeling the PFNS is needed (parameter optimization and model improvements) to reliably use it for evaluations.

07 ISOTOPE AND RADIATION SOURCES↗

Light Water Reactor Sustainability Program: Use of Time Distributions to Predict Operator Procedure Performance in Dynamic Human Reliability Analysis

The Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER) framework affords software capable of conducting human reliability analysis (HRA) using a dynamic approach built around operating procedures (OPs) from nuclear power plants (NPPs). Previous HUNTER reports document the development of this software tool, the coupling of HUNTER to the simulator code, the collection of operator performance data by using simulators to calibrate HUNTER models, and linking HUNTER to probabilistic risk assessment (PRA) software. The present report largely addresses two topics. The first is a new function in HUNTER called the HUNTER Procedure Performance Predictor (P3). HUNTER P3 uses HUNTER’s built in Monte Carlo tools featuring human performance variability to identify potential error traps in procedures. The second topic is time distribution analysis to generate time inputs for dynamic HRA. The current analysis was performed to investigate time distributions for task primitives, which are the minimum task unit of analysis used in dynamic HRA modeling. Using the time distribution data, the elapsed time for human actions in an extended loss of AC power (ELAP) scenario is then investigated. Time data and prediction are essential for modeling procedure performance.

99 GENERAL AND MISCELLANEOUS↗

Wing design by numerical optimization

A study was conducted to assess the feasibility of performing computerized wing design by numerical optimization. The design program combined a full potential, inviscid aerodynamics code with a conjugate gradient optimization algorithm. Three design problems were selected to demonstrate the design technique. The first involved modifying the upper surface of the inboard 50% of a swept wing to reduce the shock drag subject to a constraint on wing volume. The second involved modifying the entire upper surface of the same swept wing (except the tip section) to increase the lift-drag ratio subject to constraints on wing volume and lift coefficient. The final problem involved modifying the inboard 50% of a low-speed wing to achieve good stall progression. Results from the three cases indicate that the technique is sufficiently accurate to permit substantial improvement in the design objectives.

Hicks, R. M.↗

Comparison of analytical and experimental supersonic aerodynamic characteristics of a forward control missile

Techniques to predict the aerodynamic characteristics of slender cruciform missiles have been developed and are constantly being updated and improved. This paper presents comparisons between analytical and experimental supersonic aerodynamic data for a class of canard-controlled missile configurations similar to the sidewinder missile. Three aerodynamic prediction computer codes, including program MISSILE2, a recently improved version of program MISSILE, are evaluated by comparison with the test data to assess their accuracy. The major emphasis is placed on the roll control characteristics. In addition, tail span optimization, longitudinal and lateral control, induced roll, and missile roll orientation effects are addressed.

Allen, J. M.↗

Pressure measurements on a thick cambered and twisted 58 deg delta wing at high subsonic speeds

A pressure experiment at high subsonic speeds was conducted by a cambered and twisted thick delta wing at the design condition (Mach number 0.80), as well as at nearby Mach numbers (0.75 and 0.83) and over an angle-of-attack range. Effects of twin vertical tails on the wing pressure measurements were also assessed. Comparisons of detailed theoretical and experimental surface pressures and sectional characteristics for the wing alone are presented. The theoretical codes employed are FLO-57, FLO-28, PAN AIR, and the Vortex Lattice Method-Suction Analogy.

Chu, Julio↗

Probability of Failure and Risk Assessment of Propulsion Structural Components

Due to increasing need to account for the uncertainties in material properties, loading conditions, or geometries, a methodology was developed to determine structural reliability and the assess the risk associated with it. The methodology consists of a probabilistic structural analysis by a probabilistic finite element computer code Nonlinear Evaluation of Stochastic Structures Under Stress (NESSUS) and a generic probabilistic material properties model. The methodology is versatile and is equally applicable to high and cryogenic temperature structures. Results obtained demonstrate that the whole issue of structural reliability and risk can be formally evaluated using the methodology developed which is inclusive of uncertainties in material properties, structural parameters and loading conditions. The methodology is described in some detail with illustrative examples.

Shiao, Michael C.↗

System simulation on advanced earth-orbital satellite systems

A system simulation capability for advanced earth-orbital satellite systems (AEUSS) is presented. It is based on analytical assessment capabilities formulated for generic satellite systems in earth-orbital missions. The resulting AEOSS is a custom tailored software coded within the framework of the Macintosh version of the Acius relational database program, 4th Dimension. It projects the required power, weight, and cost for a satellite system. These key elements are computed on the component and subsystem levels, and then the system level. Selected performance analyses for essential components and subsystems are included.

Lee, Hwa-Ping↗

Transition to chaos in an open unforced 2D flow

The present numerical study of unsteady, low Reynolds number flow past a 2D airfoil attempts to ascertain the bifurcation sequence leading from simple periodic to complex aperiodic flow with rising Reynolds number, as well as to characterize the degree of chaos present in the aperiodic flow and assess the role of numerics in the modification and control of the observed bifurcation scenario. The ARC2D Navier-Stokes code is used in an unsteady time-accurate mode for most of these computations. The system undergoes a period-doubling bifurcation to chaos as the Reynolds number is increased from 800 to 1600; its chaotic attractors are characterized by estimates of the fractal dimension and partial Liapunov exponent spectra.

Pulliam, Thomas H.↗

On the assessment of visual communication by information theory

This assessment of visual communication integrates the optical design of the image-gathering device with the digital processing for image coding and restoration. Results show that informationally optimized image gathering ordinarily can be relied upon to maximize the information efficiency of decorrelated data and the visual quality of optimally restored images.

Huck, Friedrich O.↗