Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code Analysis”

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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

SAS4A/SASSYS-1 Modeling Improvements for the Transition to Natural Circulation

SAS4A/SASSYS-1 (SAS) is a simulation tool used to perform deterministic analyses of anticipated events as well as design basis and beyond design basis accidents for advanced liquid-metal-cooled nuclear reactors. With its origin as SAS1A in the late 1960s, the SAS series of codes has been under continuous use and development for over fifty years and represents a critical investment in safety analysis capabilities for the U.S. Department of Energy. In recent years, SAS has undergone a number of improvements to enable improved safety analyses that meet end users’ modernized needs while complying with the current regulatory environment. Improvements made in versions 5.6 and 5.7 released within the last year include the development of anisotropic Reynolds number dependent loss coefficients throughout the core and heat transport systems, the ability to distinguish the transition friction factor from the fully developed laminar and turbulent friction factors, and timedependent direct coolant and wall heating for pipe-like elements in the heat transport systems. While it was possible to capture loss coefficients, friction factors, and heat transfer from an element to a heat sink within SAS in previous versions of the code, users were required to end the simulation and restart it to adjust the input to account for any significant changes to the values during the transient. With these improvements, users can better capture flow reversal, pump heating, and the transition from forced to natural circulation without being limited to constant orifice coefficients, constant heat sinks, or the need to restart the simulation and modify input. In order to demonstrate the application of these improvements, a loss of flow transient is simulated for the Advanced Burner Test Reactor (ABTR).

SAS4A/SASSYS-1↗
DOE OSTI2022

Density reconstruction in convergent high-energy-density systems using x-ray radiography and Bayesian inference

X-ray radiography is a technique frequently used to diagnose convergent high-energy-density (HED) systems, such as inertial confinement fusion implosion, and provides unique information that is not available through self-emission measurements. In this work, we investigate the scope and limits of that information using a radiography simulation combined with the Bayesian inference workflow. The accuracy of density reconstruction from simulated radiographs of spherical implosions driven with 27 kJ laser energy is assessed, including the increase or decrease in accuracy due to the addition of Lagrangian marker layers, Poisson noise, and improved prior information. This work is the first to present the full uncertainty distributions inferred from radiography analysis in HED systems and demonstrates the importance of constructing the full posterior probability density, as opposed to a point estimate, due to the modal structure of the likelihood surface introduced by typical experimental noise sources. This general methodology can be used both for robust analysis of radiographic data and for an improved design of radiography experiments by modeling the full experimental system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)↗
DOE OSTI2024

TDCOSMO - XVI. Measurement of the Hubble constant from the lensed quasar WGD 2038–4008

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant (H0). The TDCOSMO Collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038−4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain DΔ t = 1.68−0.38+0.40 Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat ΛCDM cosmology, we constrain the Hubble constant to be H0 = 65−14+23 km s−1 Mpc−1. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.Key words: gravitational lensing: strong / cosmological parameters / distance scale⋆ Corresponding author; kcwong19@gmail.com.⋆⋆ NHFP Einstein fellow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

COMPUTATION FLUID DYNAMICS ANALYSIS FOR GENERIC SMALL MODULAR REACTOR CONTAINMENT SEPARATE EFFECTS TEST

It is desirable for fourth-generation Small Modular Reactors to be passively cooled in standard and accident operations. Passive Containment Cooling Systems can reject heat from the containment structure, without using pumps or blowers. The targeted design containment structure is a large, domed, stainless steel, cylindrical vessel. In a postulated Design Basis Accident, steam will flash inside containment. Steam condensation occurs on the inner containment wall and transfers heat through the steel containment into a large body of water known as the annular reservoir (AR) surrounding the vessel serving as the ultimate heat sink. Natural circulation drives the flow in the AR and heat will be released to the environment by evaporation of water. Unique containment geometry requires a separate effects test (SET) facility for the verification and validation of the computer code and evaluation model development and assessment for reactor licensing efforts. In this study, STAR-CCM+, a computational fluid dynamics (CFD) code was used to inform the decision-making process on the design of the SET. The CFD simulation modeled, a two-phase turbulent flow with fluid film development and heat transfer for different containment geometries. The Reactor Excursion and Leak Analysis Program will also be used in a code-to-code verification against the CFD results.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2023

Niowave Neutron Source Converter: Lead-Bismuth-Eutectic (LBE) Windowless Target Design and Evaluation

Los Alamos National Laboratory (LANL) is working with Niowave on the design and evaluation of their lead-bismuth-eutectic (LBE) windowless target (i.e., neutron source converter). Niowave plans to use 200 kW electron beam at 40 MeV beam energy to produce neutrons by photonuclear reaction with LBE. Then, the neutrons undergo fission at the surrounding uranium target assembly (UTA) to produce Molybdenum 99 (Mo-99) as a fission product, which eventually decays to Technetium-99 (Tc-99m). Tc-99m is one of important radioisotopes that is used for medical diagnostics. LANL conducted 3D multiphysics analysis for the Niowave neutron converter design and provided design assessment in thermohydraulic aspects. LANL conducted radiation transport calculations using Monte-Carlo N-Particle (MCNP) code with unstructured meshing scheme. The 3D volumetric heating profiles in the LBE and Stainless-Steel (SS) housing were imported into multiphase computational fluid dynamics (CFD) to obtain 3D temperature profiles of LBE and SS through conjugate heat transfer (CHT) analysis. The key findings are: LBE film thickness at the center of the beam is approximately 1.6 cm with a maximum LBE velocity of approximately 1.8 m/s, which is below a 2 m/s limit to avoid erosion issues on supporting structures; Heat deposition in the LBE peaks at ~1 cm depth from the LBE free-surface because of the forward interactions of electron, photon, and neutron with LBE; LBE maximum temperature is ~360 °C, which is below LBE evaporation initiative temperature, ~450 °C; LBE-SS interface temperature is ~350 °C, which is below the safety thermal limit to prevent severe corrosion on SS. The results indicate that Niowave’s neutron converter design satisfies both hydraulic and thermal criteria for safe operation. By virtue of such computational analysis, Niowave can move toward establishing an experimental setup to experimentally test their LBE neutron converter. The following sections describe the detailed work done by LANL.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Opacity on NIF: Anchor 2 Iron (Level 2 Milestone 7118 Report)

In FY2020 we fulfilled the exit criteria of milestone 7118. We fielded iron opacity experiments at higher densities and temperatures approaching anchor 2 conditions and made good progress on overcoming the problems that we found. This report describes in detail the progress that has been made. We presented the results in a LANL seminar on August 10 and we gave two talks on our results to the National Opacity Workshop Series, one on July 9 and the other on September 14. Before discussing the anchor 2 measurements, we present an update of the NIF anchor 1 data. We did a new analysis of the anchor 1 data and compared the results to a series of calculations using the LANL ATOMIC code. Figures 1 and 2 compare the data with the best-fit to calculations. This result is significantly improved over the result published in Atoms. The differences are small for the iron quasi-continuum between 8-9.5 Å and for the Mg Heα and Lyα lines. The opacities near the centers of the iron bound-bound features also match fairly well in the 9.5-12.5 Å range. However, the opacity windows between the large iron bound-bound features are deeper in the Atomic predictions than the measurement. Furthermore, large discrepancies exist for the Mg He β, γ, and δ lines and for the short-wavelength quasicontinuum. The experiment and analysis refinements described below have set the stage for unraveling these discrepancies. Specifically, we can repeat the Anchor 1 experiments using improved backlighter, hohlraum, and spectrometers to obtain reduced backgrounds and better spectral resolution. Exploiting the AlMg calibration experiments and repeating the Anchor 1 measurements multiple times will enable formal and rigorous uncertainty determination. This well-defined path gives us confidence that high quality opacity data can be obtained on NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Total System Performance Ratio—A Systems Based Approach for Evaluating HVAC System Efficiency

The prescriptive path is the most widely used approach for commercial code compliance in the United States. Though easy to implement, prescriptive approaches do not typically discriminate between minimally compliant, high-performing and poorly performing HVAC system configurations. Hence, to meet aggressive energy and carbon reduction goals, it is clear that energy codes will need to transition from prescriptive to performance-based approaches, a transition that is riddled with several challenges. This paper discusses a new HVAC system-based performance approach (HVAC System Performance) which provides a simpler solution to HVAV system evaluation compared to whole building performance, while keeping tradeoffs limited to specific building systems. The Total System Performance Ratio (TSPR) is a metric for evaluation of overall system efficiency instead of individual component efficiency, a solution which could also eventually facilitate the transition to a 100% performance-based code structure. TSPR is a ratio that compares the annual heating and cooling load of a building to the annual energy consumed by the building’s HVAC system. A calculation software tool has been developed for determining a building’s TSPR. Already incorporated into the 2018 Washington State Energy Code, this approach is also being evaluated by ASHRAE Standard 90.l Project Committee and has the potential to provide a comprehensive performance-based approach for HVAC system evaluation and analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Navigating United States Standards and Regulation for Digital Energy Systems

This report provides an analysis of the U.S. standards and regulatory landscape for digital energy systems, focusing on cybersecurity, safety, and reliability requirements. It examines the interplay between federal mandates, state regulations, voluntary industry standards, and utility-specific policies, highlighting critical gaps between compliance and real-world risk mitigation. While NERC CIP standards enforce cybersecurity for Bulk Electric System assets, distribution-level infrastructure and emerging technologies often fall outside mandatory oversight, creating vulnerabilities. The report identifies systemic challenges such as reliance on self-attestation, uneven state adoption of safety codes, and lagging standards for advanced technologies like battery energy storage and inverter-based resources. Through a detailed gap analysis, it underscores the urgency of proactive risk-based approaches, independent verification, and strategic engagement with state and federal entities. Recommendations include adopting tiered security frameworks, strengthening procurement practices, and addressing emerging technology risks to ensure resilient and secure digital energy infrastructure. This guidance is intended for utilities, regulators, and stakeholders navigating compliance obligations and seeking to enhance cybersecurity and safety beyond minimum standards.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

MATBOX, an Open-Source Microstructure Analysis Toolbox for Meshing, Generation, Segmentation, and Characterization of 3D Heterogenous Volumes

Battery performance is strongly correlated with electrode microstructural properties. To account for its impact, lithium-ion battery (LIB) models either abstract the microstructural heterogeneity of composite electrodes using effective macroscopic properties (macro- or meso- scale models) or directly solve the system of equations on the microstructure geometry or mesh (microstructure-scale models). Therefore, to be adequate, both families of models require information from the microstructure geometry, which can be provided by the numerical tool presented in this work. MATBOX is a MATLAB open-source application [1] developed by NREL for performing various microstructure-related tasks including microstructure numerical generation, image filtering and microstructure segmentation, microstructure characterization and correlation, visualization, and microstructure meshing. MATBOX was originally developed for the analysis of LIB electrode microstructures; however, the algorithms provided by the toolbox are widely applicable to other heterogeneous materials. The toolbox provides a user-friendly experience thanks to a Graphical-User Interface, requires no coding by the user, and is well documented. This presentation will illustrate various MATBOX features for the characterization of a LIB electrode, including a fully automated Representative Volume Element (RVE) analysis, the numerical generation of complex 'virtual' microstructure, including dual-layer electrodes and carbon-binder additive phase, and the meshing of a complex NMC/graphite full cell microstructure suitable for 3D finite-element modeling. Other modules (segmentation, visualization, and correlation) will be briefly presented. Thanks to its modular, open-source approach, MATBOX can easily incorporate third-party algorithms to eventually build a standard in the field that will benefit the whole scientific community. Effective diffusion coefficient [2], additive phase numerical generation [3], and meshing [4] third-party algorithms have been already integrated in the toolbox with more to come.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dynamic responses of submerged components in advanced reactors: experimental and numerical studies

The seismic design of an advanced nuclear reactor must consider the interaction of vessel internal components with the surrounding coolant: fluid–structure interaction (FSI). Available analytical solutions for FSI of submerged components do not accommodate multiple-component, intense seismic inputs and complex reactor and internal geometries. Physical testing of reactor vessels and internals to inform seismic design is impractical and cost-prohibitive, leaving the use of verified and validated, robust numerical models as the only plausible option for analysis and design. Physical data that could be used for validating such numerical models for multi-component shaking are not available. This article describes an experimental program performed on a 6-degree-of-freedom earthquake simulator to generate data that could support validation of seismic FSI numerical models for submerged components in commercial finite element codes. A scaled model of a base-supported reactor vessel with simplified representations of submerged internals was tested to generate submerged-component response histories for a range of seismic inputs. The generated data were used to validate numerical models in the finite element code LS-DYNA. Numerical models were validated for calculation of hydrodynamic pressure responses on internals, in-water frequencies of internals, and acceleration responses of internals. The generated data and the analysis recommendations could aid engineering analysts designing submerged components and systems for seismic effects.

Engineering↗
DOE OSTI2024

FEM Analysis of Hybrid LTS/HTS Cos-Theta Dipole Magnet With Heterogeneous Cable Model

The Lawrence Berkeley National Laboratory (LBNL) and the National High Magnetic Field Laboratory (NHMFL) have published results on Bi-2212 superconductive magnets realized and tested in the canted cosine-theta and solenoid designs, respectively. Fermilab is now preparing for the assembly of the first Bi-2212 stress-managed cosine-theta insert magnet. The insert will be part of the first hybrid cosine-theta magnet made of Nb$_3$Sn outer layers within the US-MDP effort to reach a 20 T bore field. This paper presents the analytical analysis of the cosine-theta Nb$_3$Sn/Bi-2212 hybrid magnet. We report the parameters, logic, and implementation method of the 2D electromagnetic and mechanical FEM analysis of the LTS/HTS hybrid magnet. Results from a detailed heterogeneous model are compared to the homogeneous model implemented in the past. A Python code has been developed to simulate the current degradation due to stresses in the detail-modeled conductor areas. The current degradation has been introduced in the simulation dynamics for the HTS conductor as an iteration process, updating the input load of Lorentz forces of the energization at each step. The magnetic and mechanical analysis results of the 2D cosine-theta LTS/HTS dipole magnet have been described and analyzed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Crowd cluster data in the USA for analysis of human response to COVID-19 events and policies

We provide data on daily social contact intensity of clusters of people at different types of Points of Interest (POI) by zip code in Florida and California. This data is obtained by aggregating fine-scaled details of interactions of people at the spatial resolution of 10 m, which is then normalized as a social contact index. We also provide the distribution of cluster sizes and average time spent in a cluster by POI type. This data will help researchers perform fine-scaled, privacy-preserving analysis of human interaction patterns to understand the drivers of the COVID-19 epidemic spread and mitigation. Current mobility datasets either provide coarse-level metrics of social distancing, such as radius of gyration at the county or province level, or traffic at a finer scale, neither of which is a direct measure of contacts between people. We use anonymized, de-identified, and privacy-enhanced location-based services (LBS) data from opted-in cell phone apps, suitably reweighted to correct for geographic heterogeneities, and identify clusters of people at non-sensitive public areas to estimate fine-scaled contacts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Scoreboard

Emerging HPC machines have given rise to enhanced compute power that far outstrips the machine's ability to save large scale results for post-processing. To combat this, in situ data analysis techniques are slowly being adopted. With in situ data management favoring workflows composed of multiple simulations and analyses connected in transit on heterogeneous machines, scientists and engineers need a tool that enables them to create data extracts, visualizations, and interactively monitor and steer their simulations. Scoreboard Phase II is a next generation analysis software that supports composite in transit workflows on heterogeneous architectures and restores interactivity to in situ data analysis through simulation monitoring and computational steering. Scoreboard provides a simulation dashboard with graphs of metrics over time, controls for setting custom simulation steering parameters, controls for managing the set of data extracts being produced in the simulation, as well as the ability to explore data extracts, all from a web browser. Realizing the vision outlined in this project required research into making a system that integrates end to end from simulations all the way to the user. In situ tools generally suffer from complexity and excessive software dependencies. Scoreboard, by contrast, is easy to build and integrate into simulation codes and it provides first class FORTRAN support. The Scoreboard library is capable of in situ and in transit data analysis that can produce data extracts commonly needed for Computational Fluid Dynamics (CFD) analysis. Simulations can transparently stage data in transit to a Scoreboard Endpoint program, which can accept their data and produce the requested data extracts. This lets simulations return to their work while the Endpoint works on the analysis. Efficiently staging the data at scale was a topic of this research. Scoreboard provides the means to let the user manage data extracts and monitor/steer many simulations from a web browser. This area of the research focused on discovery of in transit network components to expose and control their steering parameters within an interactive browser-based user interface that includes: system topology, gathered metrics, notifications, dynamically-generated steering controls, and exploration of visualization data products.

Whitlock, BradJoseph [Intelligent Light] (00000001↗
DOE OSTI2025

HFBTHO-AD: Differentiation of a nuclear energy density functional code

The HFBTHO code implements a nuclear energy density functional solver to model the structure of atomic nuclei. HFBTHO has previously been used to calibrate energy functionals and perform sensitivity analysis by using derivative-free methods. To enable derivative-based optimization and uncertainty quantification approaches, we must compute the derivatives of HFBTHO outputs with respect to the parameters of the energy functional, which are a subset of all input parameters of the code. Here, we use the algorithmic/automatic differentiation (AD) tool Tapenade to differentiate HFBTHO. We compare the derivatives obtained using AD against finite-difference approximation and examine the performance of the derivative computation.

Algorithmic differentiation↗
DOE OSTI2020

Scalability of high-performance PDE solvers

Performance tests and analyses are critical to effective high-performance computing software development and are central components in the design and implementation of computational algorithms for achieving faster simulations on existing and future computing architectures for large-scale application problems. In this article, we explore performance and space-time trade-offs for important compute-intensive kernels of large-scale numerical solvers for partial differential equations (PDEs) that govern a wide range of physical applications. We consider a sequence of PDE-motivated bake-off problems designed to establish best practices for efficient high-order simulations across a variety of codes and platforms. We measure peak performance (degrees of freedom per second) on a fixed number of nodes and identify effective code optimization strategies for each architecture. In addition to peak performance, we identify the minimum time to solution at 80% parallel efficiency. The performance analysis is based on spectral and p-type finite elements but is equally applicable to a broad spectrum of numerical PDE discretizations, including finite difference, finite volume, and h-type finite elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Thermal Spectrum Molten Salt-Fueled Reactor Reference Plant Model

This report details a fully coupled neutronics thermal hydraulics reference plant model for a thermal spectrum molten salt-fueled reactor. The Multiphysics model is developed on the Nuclear Regulatory Commission's Comprehensive Reactor Analysis Bundle available on the Idaho National Laboratory high-performance computer, which natively and seamlessly couples Griffin and Sam codes of the Multiphysics Object-Oriented Simulation Environment (MOOSE) based applications. Griffin provides the reactor physics capabilities, k-eigenvalue, delayed neutron precursor concentrations, and transient solutions for stationary and flowing fuels. Sam provides the solution of mass, momentum, and energy conservation equations of the whole system with a set of system components for one-dimensional single-phase flow. The neutronics feedback model relies primarily on fuel and moderator temperatures, fuel salt density, and delayed neutron precursors distributions. This report presents our results for the coupled steady state core and an unprotected loss of flow event. Although this model is prototypical regarding capabilities in the Comprehensive Reactor Analysis Bundle, its results are consistent with experimental data of the pump transient tests that were performed at zero power of the Molten Salt Reactor Experiment (MSRE).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

FRAM Isotopic Analysis of High-Resolution CZT Data [Slides]

Abstract: H3D Inc. introduced a new type of CZT detectors with resolution several times better than that of a typical CZT; The current FRAM v.6.1 can analyze uranium spectra with minimal modification of the parameter sets and get acceptable results; The exceptionally large high-energy tail of the peaks makes it difficult to achieve good uranium results and impossible to successfully analyze plutonium spectra; The FRAM code needed to be modified to better fit the high-energy tail of the peaks to analyze plutonium data and to get better uranium analysis results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗