Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy analysis at scale”

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 469 records · Page 26

AutoBEM: A scalable framework for nationwide building energy simulation and retrofit evaluation in the United States

This paper presents AutoBEM, an integrated, automated framework for nationwide building energy modeling and retrofit evaluation in the United States. Unlike prior UBEM platforms that either rely primarily on representative stock sampling or operate at city scale, AutoBEM automates the generation of building-resolved, physics-based EnergyPlus/OpenStudio simulation models at national scale using GIS-derived geometry, prototype-based assumptions, and standardized scalable workflows. Leveraging the Model America dataset and high-performance computing, AutoBEM generates and simulates energy models for 122.9 million buildings, representing 97.8% of the U.S. building stock. These models are being made publicly and freely available as the Model America v1.0 (MAv1) dataset. AutoBEM supports detailed, building-level assessments of energy consumption, CO2 emissions, and post-processed anthropogenic heat emissions (AHE), and evaluates 151 energy conservation measures (ECMs) using localized utility pricing and building characteristics. In addition, AutoBEM incorporates both typical and future climate conditions through integration with Typical Meteorological Year (TMY) and Future TMY (fTMY) weather data derived from IPCC scenarios. In a case study of Phoenix, Arizona, AutoBEM identified several high-efficiency HVAC upgrades and selected envelope measures with short modeled payback periods (1.5 years) for certain building types and standards. Simulations under future climate scenarios (SSP5–RCP8.5) project an 11.3% increase in electricity use and a 32% reduction in natural gas demand by 2100, underscoring the need for climate-adaptive retrofit planning. By enabling reproducible, bottom-up, and location-specific analysis at scale, AutoBEM provides a step toward a national digital twin of the built environment and supports data-driven screening and planning for decarbonization, resilience, and energy equity.

Li, Hang [ORNL] (ORCID:0000000306001920)↗

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗

Generation of lognormal synthetic Lyman-α forest spectra for P 1 D analysis

The one-dimensional flux power spectrum (P 1 D ) of the Lyman-α forest probes small-scale structure in the intergalactic medium (IGM) and is therefore sensitive to a variety of cosmological and astrophysical parameters. These include the amplitude and shape of the matter power spectrum, the thermal history of the IGM, the sum of neutrino masses, and potential small-scale fluctuations due to the nature of dark matter. However, P 1 D is also highly sensitive to observational and instrumental systematics, making accurate synthetic spectra essential for validating analyses and quantifying these effects, especially in high-volume surveys like the Dark Energy Spectroscopic Instrument (DESI). We present an efficient lognormal mock framework for generating one-dimensional Lyman-α forest spectra tailored for P 1 D analysis. Our method captures the redshift evolution of the mean transmitted flux and the scale-dependent shape and amplitude of the one-dimensional flux power spectrum by tuning Gaussian field correlations and transformation parameters. Across the DESI Early Data Release (EDR) redshift range (2.0 ≤ z ≤ 3.8), and a wide range of scales (10 -4 s km -1 ≤ k ≤ 1.0 s km -1 ), our mocks recover the mean flux evolution with redshift to sub-percent accuracy, and the P 1 D at the percent level. Additionally, we discuss potential extensions of this framework, such as the incorporation of astrophysical contaminants, continuum uncertainties, and instrumental effects. Such improvements would expand its utility in ongoing and upcoming surveys and enable a broader range of validation efforts and systematics studies for P 1 D inference and precision cosmology.

Lyman alpha forest↗

Anisotropic Energy Transfer and Conversion in Magnetized Compressible Turbulence

We present a spatial filtering (or coarse-graining) analysis on 3D magnetized magnetohydrodynamic (MHD) turbulence simulations. The filtered compressible MHD formulae show transfer of kinetic and magnetic energies from large to small scales, as well as energy conversion between kinetic, magnetic, and thermal energies. The anisotropic filtering enables separate analyses of the energy flows perpendicular and parallel to the global mean magnetic field. Anisotropy in energy cascade is demonstrated by the larger perpendicular energy cascade rate and also the larger perpendicular wavenumbers associated with the peak energy transfer rate. We also find that the “inertial range” along the parallel (perpendicular) direction in the anisotropic energy cascade formulation is no longer strictly dissipation-free, because it includes the dissipation in the perpendicular (parallel) direction. A change in the driving force (kinetic only versus kinetic and magnetic) affects the energy conversion between kinetic and magnetic energies. While the compressibility of the driving force changes the partition of different channels of energy transfer and conversion, and also increases the total energy transfer rate, the global energy flow remains unaffected by compressibility qualitatively. Our analysis can be applied to multispacecraft observations of turbulence in the solar wind or a planetary magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Ion beam probing of electrostatic fields

The determination of a cylindrically symmetric, time-independent electrostatic potential V in a magnetic field B with the same symmetry by measurements of the deflection of a primary beam of ions is analyzed and substantiated by examples. Special attention is given to the requirements on canonical angular momentum and total energy set by an arbitrary, nonmonotone V, to scaling laws obtained by normalization, and to the analogy with ionospheric sounding. The inversion procedure with the Abel analysis of an equivalent problem with a one-dimensional fictitious potential is used in a numerical experiment with application to the NASA Lewis Modified Penning Discharge. The determination of V from a study of secondary beams of ions with increased charge produced by hot plasma electrons is also analyzed, both from a general point of view and with application to the NASA Lewis SUMMA experiment. Simple formulas and geometrical constructions are given for the minimum energy necessary to reach the axis, the whole plasma, and any point in the magnetic field. The common, simplifying assumption that V is a small perturbation is critically and constructively analyzed; an iteration scheme for successively correcting the orbits and points of ionization for the electrostatic potential is suggested.

Persson, H.↗

A Cross-correlation Study of High-energy Neutrinos and Tracers of Large-scale Structure

The origin of the bulk of the astrophysical neutrinos detected by the IceCube Observatory remains a mystery. Previous source-finding analyses compared the directions of IceCube events and individual sources in astrophysical catalogs. The source association method is technically challenging when the number of source candidates is much larger than the number of the observed astrophysical neutrinos. In this paper, we show that in this large source number regime, a two-point cross-correlation analysis of neutrino data and source catalog can instead be used to constrain potential source populations for the high-energy astrophysical neutrinos, and provide spatial evidence for the existence of astrophysical neutrinos. We present an analysis of the cross-correlation of the IceCube 2010–2012 point-source data and a WISE–2MASS galaxy sample. While we find no significant detection of cross-correlation with the publicly available neutrino data set, we show that, when applied to the full IceCube data, which has a longer observation time and higher astrophysical neutrino purity, our method has sufficient statistical power to detect a cross-correlation signal if the neutrino sources trace the large-scale structure of the universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Scaling patch analysis of turbulent planar plume

Proper scaling in turbulent planar plumes is investigated in this study using a scaling patch approach. Based on the scaled boundary conditions, a proper velocity scale for the mean axial flow is the plume centerline velocity U ref = U ctr , and a proper temperature scale for the temperature excess is Θ ref = T ctr – T ∞ , where T ctr is the plume centerline temperature and T ∞ is the ambient fluid temperature. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity, mean momentum, and mean energy equations, respectively, the following is found: (1) a proper scale for the mean transverse flow is V ref = ( d δ / d x ) U ctr , where d δ / d x is the growth rate of the plume width. (2) A proper scale for the Reynolds shear stress is R vu , ref = U ctr V ref = ( d δ / d x ) U ctr 2 , a mix of the scales for the mean axial and transverse flows. (3) A proper scale for the turbulent heat flux is R v θ , ref = V ref Θ ctr , a mix of the scales for the mean transverse flow and mean temperature excess. The mean transverse flow thus plays a critical role in the scaling of turbulent planar plumes. Approximate functions are developed for the scaled mean transverse flow, Reynolds shear stress, and turbulent temperature flux, and are found to agree favorably with experimental and numerical simulation data. The integral analysis of the mean momentum equation yields a Richardson number Ri, which remains invariant in the axial direction. The Richardson number is defined as Ri = def g β Θ ctr δ t / ( U ctr V ref ) ≈ 1 / 2 , where g is the gravitational acceleration, β is the thermal expansion coefficient, and δ t is the plume half-width based on the mean temperature profile. This Richardson number arises directly from the scaling patch analysis of the mean momentum equation, including both the streamwise and transverse velocity scales.

42 ENGINEERING↗

Three-Dimensional Grid Visualization for Planning Activities: A Dubai Case Study

National Laboratory of the Rockies (NLR), in collaboration with the Dubai Electricity and Water Authority (DEWA) and Infra-X, has undertaken the Energy Visualization Analysis Project. The aim of this project is to enhance analytical and 3D visualization capabilities for distribution network planning and renewable energy integration. As modern grid continues to evolve with large-scale solar PV deployment and emerging distributed energy resources (DERs), the ability to effectively analyze, visualize, and communicate complex grid behaviors has become increasingly critical. The project focuses on developing empirical use cases based on real distribution feeder data and engineering workflows, ensuring the outcomes are directly aligned with operational environment. Through time-series power flow simulations and nodal hosting capacity analysis, the study quantifies the impacts of high PV penetration on voltage and thermal limits within representative 11 kV feeders. These analyses identify specific nodes and conditions where DER integration challenges arise. Furthermore, a Battery Energy Storage System (BESS) optimization algorithm was applied to determine the optimal size and placement of storage systems that can mitigate network constraints and enhance hosting capacity. The comparative results between base-case and BESS-augmented scenarios clearly demonstrate improvements in network stability and load management efficiency. In parallel, the NLR team developed an immersive 3D visualization framework, enabling interactive exploration of grid simulations using commodity head-mounted display (HMD) systems. This framework transforms conventional 2D simulation data into spatially intuitive visual environments - allowing engineers to analyze feeder conditions, PV hosting potential, and BESS effects in real time. This report represents the first foundational phase in establishing a visualization-driven analytical ecosystem. It provides a methodological foundation for data integration, visualization architecture, and simulation-based decision support, paving the way for large-scale adoption of immersive visualization across DEWA's Smart Grid Initiative, R&D activities, and future network resilience studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Studying the Peak Force of Magnetic Linear Motion Devices Using Dimensional Neutrality

The merits of different sizing equations for linear motion magnetic force devices are assessed. A peak force scaling analysis for a linear magnetic coupling with respect to the volumetric force density, mass force density, magnetic shear stress and energy density is compared. By utilizing a 3-D analytic based model of a Halbach cylinder magnetic coupling it is shown that the maximizing of volumetric and mass force density is not an appropriate sizing metric for a linear magnetic actuator. Only energy density and magnetic shear stress are dimensionally neutral sizing metrics for linear actuators.

16 TIDAL AND WAVE POWER↗

SCALE Code System

The SCALE Code System is a widely used modeling and simulation suite for nuclear safety analysis and design that is developed, maintained, tested, and managed by the Reactor and Nuclear Systems Division (RNSD) of Oak Ridge National Laboratory (ORNL). SCALE provides a comprehensive, verified and validated, user-friendly tool set for criticality safety, reactor and lattice physics, radiation shielding, spent fuel and radioactive source term characterization, and sensitivity and uncertainty analysis. Since 1980, regulators, licensees, and research institutions around the world have used SCALE for safety analysis and design. SCALE provides an integrated framework with dozens of computational modules, including three deterministic and three Monte Carlo radiation transport solvers that are selected based on the desired solution strategy. SCALE includes current nuclear data libraries and problem-dependent processing tools for continuous-energy (CE) and multigroup (MG) neutronics and coupled neutron-gamma calculations, as well as activation, depletion, and decay calculations. SCALE includes unique capabilities for automated variance reduction for shielding calculations, as well as sensitivity and uncertainty analysis. SCALE’s graphical user interfaces assist with accurate system modeling, visualization of nuclear data, and convenient access to desired results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Solar Futures Study Databook

The Solar Futures Study explores pathways for solar energy to drive deep decarbonization of the U.S. electric grid and considers how further electrification could decarbonize the broader energy system. This workbook contains data behind most figures in the U.S. Department of Energy's Solar Futures Study. Data are not provided for figures with data labels and certain figures generated with larger data sets. These data will be made available upon request.

14 SOLAR ENERGY↗

A Method for Calculating Strain Energy Release Rates in Preliminary Design of Composite Skin/Stringer Debonding Under Multi-Axial Loading

Three simple procedures were developed to determine strain energy release rates, G, in composite skin/stringer specimens for various combinations of unaxial and biaxial (in-plane/out-of-plane) loading conditions. These procedures may be used for parametric design studies in such a way that only a few finite element computations will be necessary for a study of many load combinations. The results were compared with mixed mode strain energy release rates calculated directly from nonlinear two-dimensional plane-strain finite element analyses using the virtual crack closure technique. The first procedure involved solving three unknown parameters needed to determine the energy release rates. Good agreement was obtained when the external loads were used in the expression derived. This superposition technique was only applicable if the structure exhibits a linear load/deflection behavior. Consequently, a second technique was derived which was applicable in the case of nonlinear load/deformation behavior. The technique involved calculating six unknown parameters from a set of six simultaneous linear equations with data from six nonlinear analyses to determine the energy release rates. This procedure was not time efficient, and hence, less appealing. A third procedure was developed to calculate mixed mode energy release rates as a function of delamination lengths. This procedure required only one nonlinear finite element analysis of the specimen with a single delamination length to obtain a reference solution for the energy release rates and the scale factors. The delamination was extended in three separate linear models of the local area in the vicinity of the delamination subjected to unit loads to obtain the distribution of G with delamination lengths. This set of sub-problems was Although additional modeling effort is required to create the sub- models, this local technique is efficient for parametric studies.

Krueger, Ronald↗

Arrhenius Analysis of the Degradation Modes in BackFLIP Study of Emerging Photovoltaic Backsheets: Preprint

With increasingly large-scale deployments of photovoltaic (PV) modules to meet global energy demands, interest in better establishing module lifespans is growing. We report on the performance of co-extruded polyolefin- (PO) based backsheets as environmentally friendly alternatives to fluoropolymer-reinforced polyethylene-terephthalate- (PET)-based backsheets using three hygrometric accelerated test conditions. After completing 4000 cumulative hours of aging, we analyzed data from electrical performance (I-V), surface roughness (gloss), and appearance (L, a*, b* color) characterizations are analyzed to quantify degradation rates, quantify the corresponding activation energy, and cross-correlate between the characteristics examined, thereby providing insights into the relationship between physical characteristics and operating performance.

BackFLIP↗

H2@Scale CRADA: CA Research Consort. (Ref. Station, Fueling Perf. Test Device, Station Cap Model)

Heavy-duty (HD) hydrogen stations are an essential part of building a clean energy infrastructure by lowering greenhouse gas emissions and criteria pollutants from the trucking industry. Clean fuel at HD stations also supports energy, environmental, and social justice by preventing pollution near the stations and on truck routes while strengthening U.S. manufacturing and providing pathways to private sector uptake. The goal of the CA Research Consortium is to continue the research partnership between national laboratories, DOE, and California public agencies for the advancement of hydrogen fueling infrastructure for HD vehicles. This cooperative work directly impacts California's transition to Zero Emission Trucks and will provide tools and information that lead to more efficient design, acceptance, and commissioning of these larger capacity, higher flow rate stations serving HD hydrogen trucks such as class 8 trucks in long haul applications.

AMR↗

A Business Case Evaluation of Gas Switching Reforming (GSR) Technology: A Promising Technology for Natural Gas Reforming with Integrated CO2 Capture

Hydrogen is essential in the transition to sustainable energy, and developing low-carbon production methods is a key research focus. Traditional steam methane reforming (SMR) dominates the hydrogen industry but contributes substantially to CO2 emissions. In response, Gas Switching Reforming (GSR) has emerged as a novel process that integrates carbon capture and utilizes process heat more efficiently. Unlike other reforming methods, GSR consolidates oxidation and reduction reactions within a single reactor, which minimizes external energy inputs and simplifies scaling. Like conventional steam methane reforming (SMR), GSR can be integrated with water-gas shift and pressure swing adsorption units for pure hydrogen production. This work presents a comprehensive business case analysis of GSR technology based on experimental results in Technology Readiness Level 3, Life Cycle Assessment (LCA) and Techno-Economic (TEA) evaluation incorporating ASPEN Plus process modeling considering different configurations and energy scenarios. The TEA incorporates data from kinetic experiments from various catalysts to evaluate the GSR process under various conditions. The goal of this work is to evaluate GSR’s potential to serve as a low-carbon alternative to SMR, focusing on global warming potential and additional impact categories to evaluate a wide spectrum of environmental impacts. Comparative assessments were conducted with SMR, chemical loop reforming (CLR), and proton exchange membrane (PEM) electrolysis to explore trade-offs across environmental metrics. The environmental impact assessment of this work encompasses the entire hydrogen production lifecycle from raw material extraction to plant decommissioning, using a cradle-to-gate boundary. Preliminary findings highlight that GSR, when integrated with low-carbon energy sources, could significantly reduce environmental impacts, making it a promising candidate for low-carbon hydrogen infrastructure. The insights from this business case evaluation aim to guide industry in scale-up and commercialization of this promising clean energy technology.

03 NATURAL GAS↗

Energy Finite Element Analysis for Computing the High Frequency Vibration of the Aluminum Testbed Cylinder and Correlating the Results to Test Data

The Energy Finite Element Analysis (EFEA) is a finite element based computational method for high frequency vibration and acoustic analysis. The EFEA solves with finite elements governing differential equations for energy variables. These equations are developed from wave equations. Recently, an EFEA method for computing high frequency vibration of structures either in vacuum or in contact with a dense fluid has been presented. The presence of fluid loading has been considered through added mass and radiation damping. The EFEA developments were validated by comparing EFEA results to solutions obtained by very dense conventional finite element models and solutions from classical techniques such as statistical energy analysis (SEA) and the modal decomposition method for bodies of revolution. EFEA results have also been compared favorably with test data for the vibration and the radiated noise generated by a large scale submersible vehicle. The primary variable in EFEA is defined as the time averaged over a period and space averaged over a wavelength energy density. A joint matrix computed from the power transmission coefficients is utilized for coupling the energy density variables across any discontinuities, such as change of plate thickness, plate/stiffener junctions etc. When considering the high frequency vibration of a periodically stiffened plate or cylinder, the flexural wavelength is smaller than the interval length between two periodic stiffeners, therefore the stiffener stiffness can not be smeared by computing an equivalent rigidity for the plate or cylinder. The periodic stiffeners must be regarded as coupling components between periodic units. In this paper, Periodic Structure (PS) theory is utilized for computing the coupling joint matrix and for accounting for the periodicity characteristics.

Vlahopoulos, Nickolas↗

A note on the interaction of unsteady flow with an acoustic liner

The effects of a mean grazing flow on the energy exchanges involved in the interaction of a bias-flow acoustic liner with, respectively, incident sound and boundary layer turbulence are contrasted. The analysis of model problems which make use of a line vortex to simulate large scale, unsteady boundary layer structures indicates that, whereas acoustic waves may be effectively attenuated, dissipation caused by 'jetting' in the apertures of the liner can result in a net transfer of energy from the mean flow to the turbulence in the boundary layer.

Howe, M. S.↗

Efficient and Reliable Power Takeoff for Ocean Wave Energy Harvesting

The project goal is to significantly improve the current ocean wave energy harvesting through innovative Power Take-off (PTO) design, advanced power electronics, and novel wave capture structures. The objective of the project is to design and demonstrate system-agnostic components for application across multiple MHK systems, and complete component designs, build scaled prototypes, and perform testing and analysis for metric validation of 25% increase in component rating/per unit cost and 50% reduction in failure rate. The major innovation of the PTO is the Mechanical Motion Rectifier (MMR) mechanism that rectifies the bi-directional oscillatory motion of the input from waves into a steady unidirectional rotation output to directly drive the electrical generator. Through this mechanism, the efficiency and the fatigue life of the PTO can be significantly improved to benefit the energy absorption and lifespan of the wave energy converters (WEC). During the period of performance, the component and system design are completed, the scaled prototypes are developed and performed testing. It is validated that the 25% increase in a component rating/per unit cost. The 50% reduction in failure rate is not directly validated by experiments, however, it can be explained qualitatively with analysis. Besides, 8 journal articles, 13 conference proceedings, 1 patent, 3 Master thesis and two Ph.D. dissertations are published based on the work related to this project. The list of all the publications can be found at the end of the project as an appendix. Over 30 students and postdocs were trained through this project. Three prototypes of 100W and 500W WECs and 10KW PTO were designed, built, and tested in ocean wave tank and using the NREL dynamometer. This project demonstrated 50-80% PTO efficiency, 90-98% power electronics efficiency, up to 66% capture width ratio in irregular waves, and 34% overall efficiency in regular waves.

16 TIDAL AND WAVE POWER↗