Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation optimization”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Variational quantum state preparation for quantum-enhanced metrology in noisy systems

Here, we investigate optimized quantum state preparation for quantum metrology applications in noisy environments. Using the QFI-OPT package, we simulate a low-depth variational quantum circuit (VQC) composed of a sequence of global rotations and entangling operations applied to a chain of qubits that are subject to dephasing noise. The parameters controlling the VQC are numerically optimized to maximize the quantum Fisher information, which characterizes the ultimate metrological sensitivity of a quantum state with respect to a global rotation. We find that, regardless of the details of the entangling operation implemented in the VQC, the optimal quantum states can be broadly classified into a trio of qualitative regimes, i.e., catlike, squeezed like, and product states, associated with different dephasing rates. Our findings are relevant for designing optimal state-preparation strategies for next-generation quantum sensors exploiting entanglement, such as time and frequency standards and magnetometers, aimed at achieving state-of-the-art performance in the presence of noise and decoherence.

quantum Fisher information↗
DOE OSTI2024

Dynamically Learning Incentives for Load Control

As electrical generation becomes more distributed and volatile, and loads become more uncertain, controllability of distributed energy resources (DERs), regardless of their ownership status, will be necessary for grid reliability. Grid operators lack direct control over end-users' grid interactions, such as energy usage, but incentives can influence behavior -- for example, an end-user that receives a grid-driven incentive may adjust their consumption or expose relevant control variables in response. A key challenge in studying such incentives is the lack of data about human behavior, which usually motivates strong assumptions, such as distributional assumptions on compliance or rational utility-maximization. In this paper, we propose a general incentive mechanism in the form of a constrained optimization problem -- our approach is distinguished from prior work by modeling human behavior (e.g., reactions to an incentive) as an arbitrary unknown function. We propose feedback-based optimization algorithms to solve this problem that each leverage different amounts of information and/or measurements. We show that each converges to an asymptotically stable incentive with (near)-optimality guarantees given mild assumptions on the problem. Finally, we evaluate our proposed techniques in voltage regulation simulations on standard test beds. We test a variety of settings, including those that break assumptions required for theoretical convergence (e.g., convexity, smoothness) to capture realistic settings. In this evaluation, our proposed algorithms are able to find near-optimal incentives even when the reaction to an incentive is modeled by a theoretically difficult (yet realistic) function.

demand response↗
DOE OSTI2025

Enhanced Monte Carlo Simulations for Electron Energy Loss Mitigation in Real-Space Nanoimaging of Thick Biological Samples and Microchips

High-resolution imaging using Transmission Electron Microscopy (TEM) is essential for applications such as grain boundary analysis, microchip defect characterization, and biological imaging. However, TEM images are often compromised by electron energy spread and other factors. In TEM mode, where the objective and projector lenses are positioned downstream of the sample, electron–sample interactions cause energy loss, which adversely impacts image quality and resolution. This study introduces a simulation tool to estimate the electron energy loss spectrum (EELS) as a function of sample thickness, covering electron beam energies from 300 keV to 3 MeV. Leveraging recent advances in MeV-TEM/STEM technology, which includes a state-of-the-art electron source with 2-picometer emittance, an energy spread of 3 × 10 -5 , and optimized beam characteristics, we aim to minimize energy spread. By integrating EELS capabilities into the BNL Monte Carlo (MC) simulation code for thicker samples, we evaluate electron beam parameters to mitigate energy spread resulting from electron–sample interactions. Based on our simulations, we propose an experimental procedure for quantitively distinguishing between elastic and inelastic scattering. The findings will guide the selection of optimal beam settings, thereby enhancing resolution for nanoimaging of thick biological samples and microchips.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Additive Manufacturing with Cellulose-Based Composites: Materials, Modeling, and Applications

Recent advances in large-scale additive manufacturing (AM) with polymer-based composites have enabled efficient production of high-performance materials. Cellulose nanomaterials (CNMs) have emerged as bio-based feedstocks due to their exceptional strength and sustainability. However, challenges such as hornification and poor dispersion in polymer matrices still limit large-scale CNM–polymer composite manufacturing, requiring novel strategies. Here, this review outlines an approach starting with atomic-level simulations to link molecular composition to key parameters like bulk density, viscosity, and modulus. These simulations provide data for finite element analysis (FEA), which informs large-scale experiments and reduces the need for extensive trials. The strategy explores how atomic interactions impact the morphology, adhesion, and mechanical properties of CNM-based composites in AM processes. The review also discusses current developments in AM, along with predictions of mechanical and thermal properties for structural applications, packaging, flexible electronics, and hydrogel scaffolds. By integrating experimental findings with molecular dynamics (MD) simulations and finite element modeling (FEM), valuable insights for material design, process optimization, and performance enhancement in CNM-based AM are provided to address ongoing challenges.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Local Electric Field Effects on Water Dissociation in Bipolar Membranes Studied Using Core–Shell Catalysts

The local electric field strength is thought to affect the rate of water dissociation (WD) in bipolar membranes (BPMs) at the catalyst–nanoparticle surfaces. Here, we study core–shell nanoparticles, where the core is metallic, semiconducting, or insulating, to understand this effect. The nanoparticle cores were coated with a WD catalyst layer (TiO2 or HfO2) via atomic layer deposition (ALD), and the morphology was imaged with transmission electron microscopy. Irrespective of the core material, these core–shell catalysts displayed comparable WD overpotentials at optimal mass loading, despite the hypothesized differences in the electric field strength across the catalyst particle suggested by continuum electrostatic simulations. Substantial atomic interdiffusion between the core and shell was ruled out by X-ray absorption spectroscopy, X-ray photoelectron spectroscopy, and diffuse reflectance optical measurements. However, the optimal mass loading of catalyst was roughly 1 order of magnitude higher for the conductive and high dielectric core materials than for the low dielectric insulating cores. These findings are consistent with the hypothesis that electric field screening within the core material focuses the electric field drop between particles such that larger film thicknesses can be tolerated. Collectively, these data support the idea that it is the local electric field at the molecular level that controls proton-transfer rates and that the metal core/dielectric-shell constructs introduced here modulate that field. Further materials and synthetic design may enable optimization of the electric field strength across the proton-transfer trajectory at the material surface.

Sarma, Prasad V↗
DOE OSTI2026

Exciton-Defect Interaction and Optical Properties from a First-Principles T-Matrix Approach

Understanding exciton-defect interactions is critical for optimizing optoelectronic and quantum information applications in many materials. However, ab initio simulations of material properties with defects are often limited to high defect density. Here, we study effects of exciton-defect interactions on optical absorption and photoluminescence spectra in monolayer MoS 2 using a first-principles T-matrix approach. We demonstrate that exciton-defect bound states can be captured by the disorderaveraged Green’s function with the T-matrix approximation and further analyze their optical properties. Our approach yields photoluminescence spectra in good agreement with experiments and provides a new, computationally efficient framework for simulating optical properties of disordered 2D materials from firstprinciples.

T-matrix↗
DOE OSTI2024

Accelerating the design of lattice structures using machine learning

Lattices remain an attractive class of structures due to their design versatility; however, rapidly designing lattice structures with tailored or optimal mechanical properties remains a significant challenge. With each added design variable, the design space quickly becomes intractable. To address this challenge, research efforts have sought to combine computational approaches with machine learning (ML)-based approaches to reduce the computational cost of the design process and accelerate mechanical design. While these efforts have made substantial progress, significant challenges remain in (1) building and interpreting the ML-based surrogate models and (2) iteratively and efficiently curating training datasets for optimization tasks. Here, we address the first challenge by combining ML-based surrogate modeling and Shapley additive explanation (SHAP) analysis to interpret the impact of each design variable. We find that our ML-based surrogate models achieve excellent prediction capabilities (R 2 > 0.95) and SHAP values aid in uncovering design variables influencing performance. We address the second challenge by utilizing active learning-based methods, such as Bayesian optimization, to explore the design space and report a 5 × reduction in simulations relative to grid-based search. Collectively, these results underscore the value of building intelligent design systems that leverage ML-based methods for uncovering key design variables and accelerating design.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Three-Dimensional Heat Flux and Thermal Analysis of Angled Tungsten Samples on DIII-D

ITER-grade tungsten and dispersoid-strengthened tungsten samples with the top surface angled at ~15° towards the incident plasma flux were exposed to 9 H-mode discharges with edge-localized modes (ELMs) in the lower divertor of DIII-D tokamak using the Divertor Material Evaluation System (DiMES). Surface damage included cracking and flaking of material on the two samples farthest away from the plasma strike point, and significant melting of the two samples closest to the strike point. Heat flux and thermal analysis tools new to DIII-D have been applied to better understand this material response and to help optimize the exposure conditions for future experiments. SMITER field-line tracing simulations based on IRTV data and EFIT equilibria estimate an average inter-ELM perpendicular heat flux, 𝑞⊥,𝑖nter−𝐸LM , on the angled surfaces of 10.1 – 19.6 MW/m² for a majority of the 9 discharges, increasing to 15.6 – 24.5 MW/m² for the single, higher-power shot where samples melted. Fast camera data showed shallow intra-ELM melting and re-solidification, which transitioned to bulk inter-ELM melting with melt motion in the 𝐽⃗ 𝑥 𝐵⃗ direction. About 50% of the protruding volume of the most affected sample was displaced via melt-motion. SIERRA thermal modeling software was able to reproduce an onset time of melting consistent with fast camera data and final sample conditions, within < 200 ms. Maximum surface temperatures of 3122 K and 2787 K are estimated for the samples farthest away from the strike point, while the closest samples achieve melting at 4067 ms and 4750 ms into the ~5000 ms plasma exposure. A +10% increase in both the SMITER 𝑞⊥,𝑖nter−𝐸LM calculations and the estimated ELM heat loads 𝑞⊥, 𝐸LM was required to achieve this result, which is within the uncertainty of the diagnostic data but likely accounts for non-ideal geometry effects plus other physics uncertainties not included in this first iteration of modeling. This work provided valuable estimates of the 3D temperature evolution to help better understand the observed surface morphology and internal recrystallization of samples, which are discussed in detail in a complementary manuscript [1]. Benchmarking efforts with more diagnosed DIII-D experiments are underway to further refine the SMITER and SIERRA models for DiMES. Future use of these tools will enable researchers to precisely target heat flux exposure conditions in DIII-D to test, but not exceed, the thermomechanical limitations of novel plasma-facing materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Existence of a robust optimal control process for efficient measurements in a two-qubit system

The verification of quantum entanglement is essential for quality control in quantum communication. In this work we propose an efficient protocol to directly verify the two-qubit entanglement of a known target state through a single-expectation-value measurement. Our method provides exact entanglement quantification using the concurrence measure without performing quantum state tomography. We prove the existence of a unitary transformation that drives a known initial state of a two-qubit system to a designated final state, where the trace over a chosen observable directly yields the concurrence of the initial state. Furthermore, we implement an optimal control process of that transformation and demonstrate its effectiveness through numerical simulations. We also show that this process is robust to environmental noise. Our approach offers advantages in directly verifying entanglement with low circuit depth, making it suitable for industrial-scale quality control of entanglement generation. Our results presented here provide mathematical justification for our earlier computational experiments.

Rodriguez, Ricardo [McPherson Coll.] (ORCID:000000↗
DOE OSTI2025

State Dependent Optimization with Quantum Circuit Cutting

Quantum circuits can be reduced through optimization to better fit the constraints of quantum hardware. One such method, initial-state dependent optimization (ISDO), reduces gate count by leveraging knowledge of the input quantum states. Surprisingly, we found that ISDO is broadly applicable to the downstream circuits produced by circuit cutting. Circuit cutting also requires measuring upstream qubits and has some flexibility of selection observables to do reconstruction. Therefore, we propose a state-dependent optimization (SDO) framework that incorporates ISDO, our newly proposed measure-state dependent optimization (MSDO), and a biased observable selection strategy. Building on the strengths of the SDO framework and recognizing the scalability challenges of circuit cutting, we propose nonseparate circuit cutting-a more flexible approach that enables optimizing gates without fully separating them. We validate our methods on noisy simulations of QAOA, QFT, and BV circuits. Results show that our approach consistently mitigates noise and improves overall circuit performance, demonstrating its promise for enhancing quantum algorithm execution on near-term hardware.

Li, Xinpeng↗
DOE OSTI2026

Digital Twin Based Condition Monitoring of LCC-LCC Inductive Power Transfer Systems

Inductive power transfer (IPT) systems provide a flexible, hands-free charging opportunity to electric vehicles (EV). The resonant network components and the transmitter and receiver coils are often subjected to high voltages or currents. Component aging in the compensation network and coils of resonant IPT systems is detrimental to the reliability and power transfer efficiency of the IPT system. Monitoring the component health of such multi-element complex systems requires robust optimization algorithms. This paper discusses condition monitoring of a resonant IPT system for an EV charger using a digital twin model. A hybrid estimation algorithm based on genetic algorithms and adaptive particle swarm optimization is developed to estimate the parameters of the digital twin model. Simulation results are used to verify the monitoring capabilities of the developed algorithm under various operating conditions of the IPT system.

Weldehawaryat, Lidya Mussie [graduate research ass↗
DOE OSTI2024

Battery State of Health Estimator: Cooperative Research and Development Final Report

NREL has developed a software tool to enable Renewance to estimate the degradation of batteries from basic information such as the type of battery and the application of that battery during its first life, so that used batteries may be evaluated for potential repurposing at low cost. This software tool utilizes NREL's BLAST-Lite battery degradation modeling code, which was updated with additional models for commercially produced lithium-ion batteries as a part of this CRADA. The software tool enables users to input details such as battery type and application so that lifetime estimates can be made without any programming or expert battery knowledge. The application input loads in saved values for parameters such as cycles per year, depth-of-discharge, and other battery operating parameters from a file defined by Renewance. These parameters may be modified to refine simulations for specific batteries. The software tool also incorporates a degradation model optimization tool, whereby existing battery degradation models may be tuned according to measured battery health. This ensures that new models still predict degradation behaviors expected from a certain battery chemistry, but with the overall degradation rate tuned to a specific battery make and model. The new model can then be saved for estimating the degradation of other similar batteries. An additional task was planned to utilize machine-learning to enable battery health diagnosis from rapid EIS measurements to accelerate the screening of used batteries. This task was not completed due to lack of available data for training a machine-learning model. CRADA benefit to DOE, Participant, and US Taxpayer: Further development of open-source software tool BLAST-Lite for predicting the lifetime of commercially produced Lithium-ion batteries (NREL SWR-22-69).

25 ENERGY STORAGE↗
DOE OSTI2026

Effects of Media Nutrient Variation on Microalgae Productivity and Economics During Semi-Continuous Cultivation

The development of large-scale microalgae growth for biofuel production is currently limited by the cost of biomass production. However, new approaches to infrastructure and cultivation practices are bringing the field closer to realization. Macronutrients in the cultivation media contribute significant costs, especially since their concentrations have not been optimized for specific strains and conditions. Environmental photobioreactors (ePBRs) were used to simulate cultivation under outdoor conditions, during which the nitrogen and phosphorus levels in the media were varied. The growth of two potential biofuel production strains, Picochlorum celeri and Tetraselmis striata, with varying nutrient inputs during summer and winter scripts, respectively, was studied. This study demonstrated that nitrogen and phosphorus in f/2 media could be reduced by more than 60% from the standard formulation, while maintaining growth rates in a semi-continuous harvesting approach. Experiments comparing the standard and reduced nutrient input concentrations were also conducted for both species in 820 L outdoor raceway ponds, in Mesa, AZ. P. celeri grown in these ponds in October had a growth rate of 10.6 ± 0.7 g/m 2 /day and 10.6 ± 0.3 g/m 2 /day for the standard and low-nutrient P. celeri ponds, respectively. T. striata grown in April–May had a growth rate of 16.6 ± 1.4 g/m 2 /day for the standard nutrient input ponds and 17.4 ± 1.1 g/m 2 /day for the low-nutrient input ponds, and in October 14.5 ± 0.6 g/m 2 /day for standard nutrient ponds and 14.4 ± 0.6 g/m 2 /day for low-nutrient ponds. These outdoor data therefore confirmed the indoor ePBR data. Techno-economic analysis shows that, if high growth rates can be attained at lower nutrient concentrations, a reduction of at least 60% in nutrient costs can be achieved. Such results highlight the importance of managing macronutrient media inputs, as these have a considerable contribution to biomass production costs in large-scale facilities. The analysis also points to the importance of maintaining high spent medium recycling rates in an industrial deployment, so as to minimize the losses of nitrogen and phosphorus compounds.

09 BIOMASS FUELS↗
DOE OSTI2024

Automated Calibration for Rapid Optical Spectroscopy Sensor Development for Online Monitoring

An automated platform has been developed to assist researchers in the rapid development of optical spectroscopy sensors to quantify species from spectral data. This platform performs calibration and validation measurements simultaneously. Real-time, in situ monitoring of complex systems through optical spectroscopy has been shown to be a useful tool; however, building calibration models requires development time, which can be a limiting factor in the case of radiological or otherwise hazardous systems. While calibration time can be reduced through optimized design of experiments, this study approached the challenge differently through automation. The ATLAS (Automated Transient Learning for Applied Sensors) platform used pneumatic control of stock solutions to cycle flow profiles through desired calibration concentrations for multivariate model construction. Additionally, the transients between desired concentrations based on flow calculations were used as validation measurements to understand model predictive capabilities. This automated approach yielded an incredible 76% reduction in model development time and a 60% reduction in sample volume versus estimated manual sample preparation and static measurements. The ATLAS system was demonstrated on two systems: a three-lanthanide system with Pr/Nd/Ho representing a use case with significant overlap or interference between analyte signatures and an alternate system containing Pr/Nd/Ni to demonstrate a use case in which broad-band corrosion species signatures interfered with more distinct lanthanide absorbance profiles. Both systems resulted in strong model prediction performance (RMSEP < 9%). Lastly, ATLAS was demonstrated as a tool to simulate process monitoring scenarios (e.g., column separation) in which models can be further optimized to account for day-to-day changes as necessary (e.g., baseline correction). Ultimately, ATLAS offers a vital tool to rapidly screen monitoring methods, investigate sensor fusion, and explore more complex systems (i.e., larger numbers of species).

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Ice storage model-predictive control in an office building with PV: scenario, error and sensitivity analysis

Thermal energy storage (TES) can enable more building-sited renewable electricity generation and lower utility bill costs for buildings owners and occupants, especially when there are high demand and variable time-of-use (TOU) charges. A model predictive control (MPC) strategy can offer additional savings over a schedule-based control with added complexity and reliance on forecasts. Here, this study examines savings for medium office buildings with chiller plants in three locations with building-installed solar photovoltaics (PV) to understand the impact of MPC. Control setpoints are fixed by a schedule-based control or optimized by nonlinear MPC. These control setpoints are actuated within EnergyPlus building models to simulate the utility cost of the chiller plant. NLP solutions can be unstable or unrealistic, but our results show that by regularizing the NLP, the solutions can be reasonably followed by the building model. MPC models make simplifications that lead to errors once the controller is participating in and changing the operation of the building. These errors average 9 % across the cases, showing that the most important parts of the system are represented. The no-thermal load costs are computed to show that the optimization can in some cases achieve both the minimum TOU and minimum monthly demand costs by demand management while reducing TOU energy costs by energy arbitrage. The MPC saves 35–66 % in the annual chiller plant operating costs, which is an additional savings above the schedule by 1–33 %. PV and TES are complementary and mostly independent, but a load with PV often results in better performance for the schedule. Our case study and sensitivity analysis show the importance of modeling and optimization for complex rates, but also the circumstances wherein a simpler strategy achieves the same performance with less potential for error.

14 SOLAR ENERGY↗
DOE OSTI2024

eReaxFF force field development for BaZr 0.8 Y 0.2 O 3-δ solid oxide electrolysis cells applications

The use of solid-oxide materials in electrocatalysis applications, especially in hydrogen-evolution reactions, is promising. However, further improvements are warranted to overcome the fundamental bottlenecks to enhancing the performance of solid-oxide electrolysis cells (SOECs), which is directly linked to the more-refined fundamental understanding of complex physical and chemical phenomena and mass exchanges that take place at the surfaces and in the bulk of electrocatalysis materials. Here, we developed an eReaxFF force field for barium zirconate doped with 20 mol% of yttrium, BaZr 0.8 Y 0.2 O 3-δ (BZY20) to enable a systematic, large-length-scale, and longer-timescale atomistic simulation of solid-oxide electrocatalysis for hydrogen generation. All parameters for the eReaxFF were optimized to reproduce quantum-mechanical (QM) calculations on relevant condensed phase and cluster systems describing oxygen vacancies, vacancy migrations, electron localization, water adsorption, water splitting, and hydrogen generation on the surfaces of the BZY20 solid oxide. Using the developed force field, we performed both zero-voltage (excess electrons absent) and non-zero-voltage (excess electrons present) molecular dynamics simulations to observe water adsorption, water splitting, proton migration, oxygen-vacancy migrations, and eventual hydrogen-production reactions. Based on investigations offered in the present study, we conclude that the eReaxFF force field-based approach can enable computationally efficient simulations for electron conductivity, electron leakage, and other non-zero-voltage effects on the solid oxide materials using the explicit-electron concept. Moreover, we demonstrate how the eReaxFF force field-based atomistic-simulation approach can enhance our understanding of processes in SOEC applications and potentially other renewable-energy applications.

08 HYDROGEN↗
DOE OSTI2024

Geometry Optimization of Cable-Based Actuation for Small-Scale Model Testing of a Floating Marine Turbine: Preprint

This research aims to apply combined wave and tidal current loads to a small-scale floating marine current turbine in a wave tank, where an actuation system applies hydrodynamic and mooring forces on the hardware based on results from a simulation. We use a real-time hybrid test setup with physical wave forcing from the wave tank and simulated current and mooring forces implemented through a tensioned cable array. For this paper, we developed an optimization algorithm that adjusts the hardware geometry of the cable array to achieve more efficient tension allocation across the cables for the loads that need to be actuated. By adjusting the points where the cables are attached to the floating platform and the angles between the platform and the winches, an optimal cable geometry can be found to minimize tension variations in the lines, maintain the desired pretension, and prevent excessive tensions. We present the optimization problem formulation, the actuation system evaluation approach, and optimization results that show effective cable actuation setups that are being considered for implementation in the wave tank tests.

cable-based actuation↗
DOE OSTI2026

Polarization transmission in the Hadron Storage Ring of the Electron-Ion Collider

The successful operation of the future Electron-Ion Collider is contingent on maintaining high hadron beam polarization up to 275 GeV. The Hadron Storage Ring lattice, however, features a symmetry-breaking interaction region that excites strong, nonsystematic spin resonances, posing a significant threat to polarization preservation. This paper systematically investigates two complementary strategies to ensure high polarization transmission. The first method involves optimizing the vertical betatron phase advance between Siberian snakes to orchestrate a cancellation of depolarizing kicks across the ring. We demonstrate through simulations that both dynamic and fixed-optics solutions based on this principle can successfully preserve polarization through the strongest resonances. The second, more powerful approach involves optimizing the snake rotation axes to suppress resonance driving terms at their source. We revisit established symmetric configurations, such as the Lee-Courant schemes, and introduce a novel, highly symmetric “Doubly Lee-Courant” (DLC) scheme, which enforces a local 𝜋 spin phase advance across every consecutive pair of snakes. Our analysis reveals a clear performance hierarchy, with the DLC configuration providing an exceptionally robust and energy-insensitive baseline for polarization preservation. We conclude that a hybrid strategy, using a DLC snake scheme as a symmetric foundation and betatron phase tuning for fine corrections, offers the most effective path forward for the EIC and future high-energy polarized-beam facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗