Multi-Objective Optimization of Energy Storage Siting and Sizing Problem
The code for the PESGM 2026 Paper Multi-Objective Optimization of Energy Storage Siting and Sizing Problem
Aboumrad, Willie [IonQ, Inc]↗
SEARCH · Engineering Papers
Results for “energy code”
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.
Low Precision and Efficient Programming Languages for Sustainable AI: Final Report for the Summer Project of 2024
This document contains all relevant material generated during the authors' summer internship at NREL in 2024. This report shows how to improve energy efficiency of a few code samples by using low-precision data types combined with mixed-precision algorithms. The main applications considered here are (i) linear system solvers using mixed precision, and (ii) neural networks using mixed precision. This report also discusses how programming languages affect energy consumption of algorithms, energy metrics for a code and tools, and the available current software and hardware infrastructure.
97 MATHEMATICS AND COMPUTING↗
Leveraging MARVEL and SPHERE to Demonstrate NEAMS Thermal Hydraulics Codes
The Single Primary Heat Extraction and Removal Emulator (SPHERE) and Microreactor Applications Research Validation and Evaluation (MARVEL) were utilized to demonstrate thermal hydraulics codes of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. A recent long-duration test of a sodium heat pipe conducted at the SPHERE facility was utilized for the validation of the heat pipe code Sockeye. The first 10 hours of operation featured a startup of the heat pipe from room temperature; this portion was modeled with Sockeye and compared to the experimental data. Measurements included several thermocouple measurements across the length of the heat pipe surface. Two of Sockeye’s heat pipe models were applied to experiment, as well as several models of the cooling at the condenser end of the heat pipe. Sockeye results gave a decent match to the experimental data, with notable discrepancy in the condenser section, where experimental data indicated a significant inactive length of the heat pipe. The MARVEL microreactor was modeled through a coupled simulation of System Analysis Module (SAM) and SubChannel Module (SCM). This multiscale coupled simulation was performed as an exercise to demonstrate the compatibility and applicability of these codes for thermal-hydraulic analysis of microreactors. The coupling between the SCM and SAM codes for multiscale modeling of MARVEL, was achieved through a domain overlapping approach. In this example a transient SAM simulation provided boundary conditions to SCM, and SCM provided pressure drop to SAM.
97 - MATHEMATICS AND COMPUTING↗
Stochastic GW -GPU: Rapid Quasi-Particle Energies for Molecules beyond 10,000 Atoms
StochasticGW is a code for computing accurate quasi-particle (QP) energies of molecules and material systems in the GW approximation. StochasticGW utilizes the stochastic Resolution of the Identity (sROI) technique to enable a massively parallel implementation with computational costs that scale semilinearly with system size, allowing the method to access systems with tens of thousands of electrons. Here, we introduce a new implementation, StochasticGW-GPU, for which the main bottleneck steps have been ported to GPUs and give substantial performance improvements over previous versions of the code. We showcase the new code by computing band gaps of hydrogenated silicon clusters (Si x H y ) containing up to 10,001 atoms and 35,144 electrons, and we obtain individual QP energies with a statistical precision of better than ±0.03 eV with times-to-solution of less than 1 h.
Thomas, Phillip S. [Lawrence Berkeley National Lab↗
Particle-in-cell simulations of burning inertial confinement fusion capsule implosions
Anomalies observed in the neutron spectral shift of high-yield shots at the National Ignition Facility (NIF) suggest the presence of suprathermal ions [E. P. Hartouni et al., Nat. Phys. 19, 72 (2023)], implying that kinetic effects play a significant role in burning inertial confinement fusion (ICF) plasmas. Furthermore, recent measurements of reaction-in-flight (RIF) neutrons offer a direct probe of the stopping power in the burning fuel region of high-energy alpha particles and up-scattered fuel ions. We have developed the particle-in-cell code PICNIC, an exactly energy-conserving particle-in-cell Monte-Carlo collision (PIC-MCC) code to simulate the burn stage in ICF. We present results from 1D spherical simulations of NIF shot N210808, which was the first to exceed the Lawson criterion for laser fusion. We find that the suprathermal ions generated by large-angle Rutherford and nuclear elastic scattering (NES) with fusion alphas produce an alpha knock-on neutron (AKN) signal consistent with the extent and relative yield of the AKN spectrum identified in ignition experiments at the NIF. Furthermore, we also find that the inclusion of large-angle scattering physics does not explain the anomalously large spectral shift observed in the experiment.
High-energy-density plasmas↗
Studying the Power Sources behind Type Ic Supernovae
Astrophysical transients can be powered by a broad range of energy sources, including shock heating (internal and external shocks), decay of radioactive isotopes, and long-lived central engines (magnetar and fallback). The dominant energy source for astrophysical transients depends on the nature of the explosive engine and its progenitor. To model all transients, light-curve codes must include all of these energy sources. Here, we present a supernova light-curve code implementing analytic source models to compare the role of different energy sources in these transients. To demonstrate the utility of this code, we conduct an extensive study of Type Ic broad-line (Ic-BL) supernovae. A diverse set of energy sources have been linked to Ic-BL supernovae, making them an excellent candidate for this light-curve code. In this paper, we explore which features of the explosion (mass, velocity, etc.) affect the Type Ic supernovae light curves, focusing on shock-interaction and radioactive-decay energy sources. Although the explosion properties under both energy sources can be tuned to match the peak emission, matching the light-curve evolution in many Ic-BL supernovae requires fine-tuned conditions. We find that shock interactions in the stellar wind are likely to be the dominant energy source at peak for these supernovae.
79 ASTRONOMY AND ASTROPHYSICS↗
HERO WEC Belt Test Data
The following submission includes raw and processed data from the 2024 Hydraulic and Electric Reverse Osmosis Wave Energy Converter (HERO WEC) belt tests conducted using NREL's Large Amplitude Motion Platform (LAMP). A description of the motion profiles run during testing can be found in the run log document. Data was collected using NREL's Modular Ocean Data AcQuisition (MODAQ) system in the form of TDMS files. Data was then processed using Python and MATLAB and converted to MATLAB workspace, parquet, and csv file formats. During Data processing, a low pass filter was applied to each array and the arrays were then resampled to common 10Hz timestamps. A MATLAB data viewer script is provided to quickly visualize these data sets. The following arrays are contained in each test data file: - Time: Unix seconds timestamp - Test_Time: Time in seconds since beginning of test - POS_OS_1001: Encoder position in degrees (the encoder is located on the secondary shaft of the spring return and is driven by the winch after a 4.5:1 gear reduction) - LC_ST_1001: Anchor load cell data in lbf - PRESS_OS_2002: Air spring pressure in psi This data set has been developed by the National Renewable Energy Laboratory, operated by Alliance for Sustainable Energy, LLC, for the U.S. Department of Energy (DOE) under Contract No. DE-AC36-08GO28308. Funding provided by the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Water Power Technologies Office.
16 TIDAL AND WAVE POWER↗
Reported Energy and Cost Savings from the DOE ESPC IDIQ Program: FY 2024
Energy Savings Performance Contracts (ESPCs)are a contractual mechanism that allow a federal agency to procure energy savings and facility improvements without upfront capital costs to reduce costs and resiliency. ESPCs are covered under FAR Part 23.2, and 42 USC § 8287. Section 8287(a)(2)(A) of Title 42 of the U.S. Code requires that each energy savings performance contract (ESPC) undergo an annual energy audit, resulting in a separate audit report for every project. The objective of the present report is to compile and analyze all annual ESPC audit reports issued between October 1, 2023, and September 30, 2024, for projects awarded under Generations 1, 2, and 3 of DOE’s ESPC IDIQ contracts. During this period, 205 measurement and verification (M&V) reports were produced for 200 projects; the total number of reports exceeds the number of projects because some projects generated more than one report(for example, a few projects measure savings twice per year and produce two audit reports annually, each covering a different six-month period). By aggregating the results from these individual audits, the report determines the portfolio-wide realization rate of energy and cost savings for all active ESPC projects awarded under DOE’s IDIQ program. For all 205audit reports, sufficient information was available to compare project-level estimated, reported, and guaranteed cost savings. Reported cost savings accounted for ESCO verified savings per each project’s M&V plan. The total reported cost savings for the period addressed were $\$$647.8million,compared with the total guaranteed cost savings of $\$$601.6million. On average across the reported projects: •ESPC contractors guaranteed 92.8% of the estimated cost savings• projects reported achieving 100.0% of the estimated cost savings• projects reported achieving 107.7% of the guaranteed cost savings. The M&V performed for the period indicated adjustments for government operations and maintenance impacts to savings amount to$\$$43.9millionandcould be restored with the original operational parameters for impacted projects. Accounting for this potential cost savings impact, these projects still realized 100.4% of the guaranteed cost savings. The information on estimated and reported energy savings was collected and compared for all 205of the reports examined. Based on site energy, estimated savings totaled 14.88million MMBtu, and reported savings were 15.33million MMBtu; 3.1% greater than the estimated energy savings. All of the reports examined contained sufficient information to calculate source energy savings. Based on site-adjusted source energy, total estimated energy savings were 20.90 million MMBtu, and reported savings were 21.22million MMBtu, 101.5% of the estimated energy savings. For water savings, the estimated savings were 11,539,055 kGal and the reported savings were 13,315,930 kGal. This means 1,776,875 kGal more water was saved than estimated, which is about 15% higher than the estimate. These results indicate that, overall, the reported energy savings slightly exceeded the estimated values, while estimated water savings significantly exceeded estimated values, suggesting that the projects achieved greater cost savings than originally projected. The total annual expense for the ESCOs to perform annual M&V audits and reporting was $\$$10.02million. Through this effort, $\$$647.8 million in annual cost savings was verified. The M&V results indicated that $\$$43.9 million of these verified savings reflected adjustments due to government operations A-6and maintenance impacts, which could be restored under the original operational parameters for the affected projects. These findings show the value of M&V that only costs 1.7%of the guaranteed cost savings to ensure guarantees are met.
29 ENERGY PLANNING, POLICY, AND ECONOMY↗
Reported Energy and Cost Savings from the DOE ESPC IDIQ Program: FY 2024
Energy Savings Performance Contracts (ESPCs) are a contractual mechanism that allow a federal agency to procure energy savings and facility improvements without upfront capital costs to reduce costs and enhance mission resiliency. ESPCs are covered under FAR Part 23.2, and 42 USC § 8287. Section 8287(a)(2)(A) of Title 42 of the U.S. Code requires that each energy savings performance contract (ESPC) undergo an annual energy audit, resulting in a separate audit report for every project. The objective of the present report is to compile and analyze all annual ESPC audit reports issued between October 1, 2023, and September 30, 2024, for projects awarded under Generations 1, 2, and 3 of DOE’s ESPC IDIQ contracts. During this period, 205 measurement and verification (M&V) reports were produced for 200 projects; the total number of reports exceeds the number of projects because some projects generated more than one report (for example, a few projects measure savings twice per year and produce two audit reports annually, each covering a different six-month period). By aggregating the results from these individual audits, the report determines the portfolio-wide realization rate of energy and cost savings for all active ESPC projects awarded under DOE’s IDIQ program.
29 ENERGY PLANNING, POLICY, AND ECONOMY↗
Fusion Neutron Generator
The proposed code, named FROG (Fusion neutron Generator) is built upon the open-source particle transport Monte Carlo toolkit Geant4. Geant4 provides C++ classes that can be leveraged to build application-specific codes dealing with the transport of particles through matter. Geant4-based codes are applied in high-energy particle physics experiments, medical applications, shielding, and space applications for example. The FROG code allows the user to define the geometry of a neutron converter device shaped as a hollow cylinder, where a neutron breeding material such as lithium deuteride (LiD) is cladded by two concentric cylinders. Such neutron converter is then placed inside a regular nuclear fission reactor, where thermal neutrons will react with the neutron breeder material (typically, Lithium 6), and through a series of reactions, will generate high-energy neutrons – neutrons whose kinetic energy are around 14 MeV. The hollowed central portion can hold a specimen that will be bombarded by high-energy neutrons created inside the neutron breeding material. Figuratively speaking, this type of device transforms neutrons from thermal (~0.625 eV) to fusion (~14 MeV) energies and is sometimes termed “fusion-to-thermal neutron converters” in the literature. The code consists of C++ source file compiled and linked to generate an executable. The user can select the dimensions of the converter (radius, length, and thickness of the breeder material), the breeder material type, the cladding material, and the specimen material that will be activated or irradiated. As input, the neutron flux for a specific location inside a reactor, for instance, positions in ATR, is required. As output, the code predicts the number of high-energy neutrons produced, the total neutron flux and fluence as well as its detailed spectrum. The physics involved in such device is very complex, as it requires modeling neutron transport, light-ion (tritons) transport, as well as fusion reactions. The Geant4 toolkit provides the required physical models.
Martin, NicholasP. [Idaho National Laboratory (INL↗
Lessons Learned and Scalability Achieved When Porting Uintah to DOE Exascale Systems
A key challenge faced when preparing codes for Department of Energy (DOE) exascale systems was designing scalable applications for systems featuring hardware and software not yet available at leadership-class scale. With such systems now available, it is important to evaluate scalability of the resulting software solutions on these target systems. One such code designed with the exascale DOE Aurora and DOE Frontier systems in mind is the Uintah Computational Framework, an open-source asynchronous many-task (AMT) runtime system. To prepare for exascale, Uintah adopted a portable MPI+X hybrid parallelism approach using the Kokkos performance portability library (i.e., MPI+Kokkos). This paper complements recent work with additional details and an evaluation of the resulting approach on Aurora and Frontier. Results are shown for a challenging benchmark demonstrating interoperability of 3 portable codes essential to Uintah-related combustion research. These results demonstrate single-source portability across Aurora and Frontier with scaling characteristics shown to 3,072 Aurora nodes and 9,216 Frontier nodes. In addition to showing results run to new scales on new systems, this paper also discusses lessons learned through efforts preparing Uintah for exascale systems.
Holmen, John [ORNL] (ORCID:0000000259342641)↗
HERO WEC V1.0 - WEC-Sim Model (July 2024)
**This submission supersedes submission MHKDR-483** This submission file contains the files that are needed to simulate NREL's HERO WEC (hydraulic and electric reverse osmosis wave energy converter). This requires the user to have already installed WEC-Sim. In addition to the standard toolboxes that are required to run WEC-Sim the user will also need the Simscape Fluids and Simscape Driveline packages. The zip file (HERO_V1_WECSim_2024.zip) contains the following: - HERO_HPTO_2024.slx: Simulink-based WEC Sim model of the first gen (V1.0) Hydraulic PTO (power take-off) that was designed for the HERO WEC. This model has been updated since submission #483 based on in-laboratory experimental results. - wecSimInputFile.m: Input file needed to run the model - userDefinedFunctionsMCR.m: MCR (multi condition run) script that is needed if a use wants to simulate multiple wave conditions. - geometry (folder): Includes the geometry file that is needed for visualization - hydroData (folder): Includes the required WAMIT data to run WEC-Sim -HydVisualization.mlx: Visualization script to plot simulation results (not needed to run)
16 TIDAL AND WAVE POWER↗
Technoeconomic Design Optimization for Fast Reactors. Part II: Impact of Technoeconomic Constraints on Optimal Design
There is a current drive toward optimizing reactors, particularly small/micro reactors to minimize cost and maximize performance. Previous work has investigated the development of technoeconomic workflows for the design optimization of pool-type fast reactors that aim to deploy into district energy grids. Initial scoping studies verified that the workflow was capable of capturing design trends throughout a variety of design configurations and problem formulations while remaining sufficiently flexible. In this paper, this methodology is applied to understand how cost functions and technoeconomic constraints can drive optimal reactor design. Specifically, the UPu10Zr-fueled fast reactor model from Part I is adapted to include changes in the fissile content limits, control rod worth limits, control rod drive cost, and assumed fuel form. In the case of constraint relaxation at fixed power (fissile content and control rod worth limits), cost sensitivities of 5% to 10% were uncovered. Multi-objective optimization at varying reactor power levels with individualized control rod drives for each assembly (as opposed to one operational and one safety drive) increased cost by $\$10$ to $\$25$ million and substantially altered the optimal core geometry, favoring geometries with substantially fewer control rod placements relative to baseline optimization. Finally, a multi-objective optimization was performed at varying power levels with the fuel form overhauled to metallic, high-assay low-enriched uranium–based U10Zr with more refined fuel cost models. In the case of uranium fueling, the costs increased by at least $50 million relative to the baseline case. Furthermore, economic fuel zoning and lower reactivity swing cores were recovered. Each case serves to demonstrate the value of applying technoeconomic workflows to initial reactor design scoping studies to better understand the trade-off for a proposed concept between different design options.
Argonne Reactor Computation (ARC) codes↗
Comparison of tungsten versus molybdenum for double shell capsules using machine learning design optimization
Double shell targets are an alternative ignition platform for inertial confinement fusion. One design consideration for double shell targets is the choice of inner shell material to help trap radiation emitted by the hot fuel to aid ignition. Materials such as molybdenum and tungsten are of interest for the inner shell layer of the targets. While molybdenum has a lower density that could inhibit instability growth and allow for radiography and code benchmarking, tungsten has a higher density that could provide more compression and confinement. These tradeoffs have been explored using optimized designs for each material. Our previous work [Vazirani et al., “Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 28, 122709 (2021); Vazirani et al., “Coupling multi-fidelity xRAGE with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 30, 062704 (2023); and Vazirani et al., “Bayesian batch optimization for molybdenum versus tungsten inertial confinement fusion double shell target design,” Stat. Anal. Data Min. 17, e11698 (2024)] resulted in a multi-fidelity Bayesian optimization framework to find yield-optimized double shell target geometries. By leveraging simulations of varying fidelities (one-dimensional and two-dimensional) to inform one another, the multi-fidelity optimization was able to optimize a design in the highest fidelity with significantly fewer simulations than would be used in a systematic parameter scan. In this work, we apply the multi-fidelity Bayesian optimization to explore the optimized designs of double shell targets with molybdenum and tungsten inner shells as well as the physics producing the high performing implosions. A physics exploration of all the simulations used in this study shows trends in designs that contribute to high yields, ion temperatures, and fuel areal densities. Comparison of molybdenum and tungsten simulations shows that they can produce similar implosion conditions with different geometries, which would be important to study in experiments. Graded density layers produce varying performances with the two materials but continue to be of interest for future studies along with studies of doped inner shell materials and applied surface roughness.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Value-added catalog of physical properties for more than 1.3 million galaxies from the DESI survey
We present an extensive catalog of the physical properties of more than a million galaxies investigated with the Dark Energy Spectroscopic Instrument (DESI), one of the largest spectroscopic surveys to date. Spanning a full range of target types, including emission-line galaxies, luminous red galaxies, and quasars, our survey encompasses an unprecedented range of spectroscopic redshifts, all the way from 0 to 6. The physical properties, such as stellar masses and star formation rates, were derived via the CIGALE spectral energy distribution (SED) fitting code accounting for the contribution coming from active galactic nuclei (AGNs). Based on the modeling of the optical-mid-infrared (grz supplemented with WISE photometry) SEDs, we studied the galaxy properties with respect to their location on the main sequence. We have revised the dependence of stellar mass estimates on model choices and on the availability of WISE photometry. Indeed, the WISE data are required to minimize the misclassification of star-forming galaxies as AGNs. The lack of WISE bands in SED fits leads to elevated AGN fractions for 68% of star-forming galaxies identified using emission line diagnostic diagrams, but this does not significantly affect their stellar mass or star formation estimates.
79 ASTRONOMY AND ASTROPHYSICS↗
A future of inertial confinement fusion without laser-plasma instabilities
From the beginning of inertial confinement fusion (ICF) research, laser-plasma instabilities excited by narrowband lasers have limited the hydrodynamic design space of all laser-based approaches to inertial fusion energy (IFE). With advances in broadband laser technologies, the next generation of ICF drivers will likely have large bandwidth that is engineered to mitigate laser-plasma instabilities, thereby expanding the hydrodynamic design space to include both robust high yields (>200−MJ) with large-energy laser systems (>4 MJ) and high gains (>10) with moderate-energy laser facilities (<2 MJ). State-of-the-art simulations indicate that laser bandwidths of a few percent are required to mitigate instabilities for IFE-relevant conditions. To test these models and demonstrate that high-bandwidth lasers mitigate laser-plasma instabilities, the Fourth-generation Laser for Ultrabroadband eXperiments will be used with the OMEGA Laser System. The goal is to provide the community with the confidence to invest in a multiple-beam high-bandwidth laser facility that will demonstrate the necessary ablation pressures for robust direct-drive ignition without detrimental levels of hot electrons that degrade the implosion performance. It is important to recognize that the highest performing ICF implosions will likely never be completely LPI free because of the significant advantages to maximizing the laser intensity.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Application of a Physics-Informed Convolutional Neural Network for Monitoring the Temperature Fields in High-Temperature Gas Reactors
Here, this work presents current advances in applying a physics-informed convolutional neural network (CNN) to evaluate temperature distributions in advanced reactors. Our goal is to demonstrate that the CNN can reconstruct temperature fields within the solid region of a prismatic fuel assembly in a high-temperature gas reactor (HTGR) with sensor data available in only a few cooling channels. Before that, we showcase the superior performance of the physics-informed CNN in comparison to a purely data-driven multilayer perceptron (MLP), considering a canonical heated channel setup. This analysis shows the advantages of our approach and justifies its choice. The datasets employed here are obtained upon numerical simulations performed with codes under the Nuclear Energy Advanced Modeling and Simulation program. This work is important, as industry experience indicates that the assembly material in HTGR concepts is prone to large thermal-mechanical loads nearing operational limits. This makes it crucial to characterize peak temperatures and their distributions near hot spots. Modern thermocouples are unreliable in these types of harsh environments because of the high neutron fluxes and elevated temperatures involved. The CNN-based field reconstruction represents an attractive solution, enabling sensor arrays in less aggressive locations and augmenting indirect predictions for less accessible regions. The results show that the CNN reduces prediction errors by orders of magnitude in comparison to the MLP, considering the simple yet well-representative heated channel case. In the case of the HTGR fuel assembly, the CNN can successfully reconstruct temperature fields over various cooling regimes. Furthermore, we also explore the algorithm’s ability to detect abnormalities. Interestingly, the CNN proves it has the capacity to detect blockage in one of the noninstrumented cooling channels.
Machine learning↗
DISARM: Target Electronic Device Informed Mitigation of Software Runtime Side-Channel Vulnerabilities
Program runtime/timing attacks exploit variations in a program’s execution times to extract sensitive information from the program (e.g. encryption keys, sensitive variable data, intellectual property). State-of-the-art solutions to runtime side-channel attacks attempt to balance the execution time of the sensitive code for different control flow paths to eliminate the timing leakage. However, during the mitigation process, most techniques do not consider the underlying hardware/device on which the target program is supposed to run on. This can lead to over-fixing (unnecessary extra operations), under-fixing (not solving the imbalance properly), and even failures. Here, we propose DISARM, a joint hardware-software methodology (unlike any existing solution) for mitigating runtime side-channel vulnerabilities that utilizes timing values from real embedded devices to generate targeted software fixes. We implement DISARM to support C/C++/Java source codes and validate it across 22 standard benchmarks. DISARM outperforms state-of-the-art solutions such as PENDULUM and DifFuzzaR in terms of execution time overhead, code size overhead, and correctness on five different embedded/edge devices.
Timing/runtime side-channel↗