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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Scalable freeform optimization of wide-aperture 3D metalenses by zoned discrete axisymmetry

We introduce a novel framework for design and optimization of 3D freeform metalenses that attains nearly linear scaling of computational cost with diameter, by breaking the lens into a sequence of radial “zones” with 𝑛-fold discrete axisymmetry, where 𝑛 increases with radius. This allows vastly more design freedom than imposing continuous axisymmetry, while avoiding the compromises of the locally periodic approximation (LPA) or scalar diffraction theory. Using a GPU-accelerated finite-difference time-domain (FDTD) solver in cylindrical coordinates, we perform full-wave simulation and topology optimization within each supra-wavelength zone. We validate our approach by designing millimeter and centimeter-scale, poly-achromatic, 3D freeform metalenses which outperform the state of the art. By demonstrating the scalability and resulting optical performance enabled by our “zoned discrete axisymmetry” (ZDA) and supra-wavelength domain decomposition, we highlight the potential of our framework to advance large-scale meta-optics and next-generation photonic technologies.

Sun, Mengdi [Wesleyan University]↗

DOE OSTI2024

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Engagement: Hyperparameter Optimization of Generative Adversarial Network Models for High-Energy Physics Simulations

We present our SciDAC FASTMath-HEP partnership results for tuning generative adversarial models (GANs) for high energy physics applications. The GANs are used in hybrid simulations to accelerate otherwise time-consuming computations. We optimize for both, prediction accuracy and variability with the goal to find GAN architectures that are reliable and robust.

high energy physics↗

DOE OSTI2025

Next-generation perovskite photovoltaics: improving, stabilizing, and lead-sealing of record-setting laboratory solar cells towards commercialization

Summary: In the proposed program we plan to improve perovskite photovoltaic performance by developing (1) orientational control of 3D/2D perovskite heterostructures to simplify device architectures, thus improving device efficiencies and stability; (2) high-throughput optical measurements and real-time device simulations for device optimization; (3) robust, dual-pronged lead-sealing and oxygen/moisture/UV barrier films for long-term stability. Specifically, we seek to develop an in-depth understanding of the perovskite film formation, interface passivation, device stability, and environmentally friendly encapsulation, which together will lead to perovskite devices with PCE of over 28%, stability of T80 at 85/85 (85% humidity at 85 degrees C) for 10,000, expected to be equivalent to T90 of 100 hours), and architectures that would satisfy the U.S. EPA and RoHS limits of lead-leaching. The knowledge generated with this project will be applicable to tandem devices with wider-bandgap perovskites.

14 SOLAR ENERGY↗

DOE OSTI2025

Evaluation of the response to electrons and pions in the scintillating fiber and lead calorimeter for the future electron-ion collider

The performance of the Baby Barrel Electromagnetic Calorimeter (Baby BCAL) — a small-scale lead-scintillating-fiber (Pb/ScFi) prototype of the GlueX Barrel Electromagnetic Calorimeter (BCAL) — was tested in a dedicated beam campaign at the Fermilab Test Beam Facility (FTBF). This study provides a benchmark for the Pb/ScFi component of the future Barrel Imaging Calorimeter (BIC) in the ePIC detector at the Electron-Ion Collider (EIC). The detector response to electrons and pions was studied at beam energies between 4 and 10 GeV, extending previous GlueX tests to a higher energy regime. The calibrated detector exhibits good linearity within uncertainties, and its electron energy resolution meets EIC requirements. The data further constrain the constant term in the energy resolution to below 1.9%, improving upon previous constraints at lower energies. Simulations reproduce key features of the electron and pion data within the limitations of the collected dataset and the FTBF test environment. Electron-pion separation in the test beam setup was analyzed using multiple methods, incorporating varying degrees of beam-related effects. The inclusion of longitudinal shower profile information enhanced the separation performance, underscoring its relevance for the full-scale BIC in ePIC. These results provide essential benchmarks for the Pb/ScFi section of the future BIC, validating detector simulations and guiding optimization strategies for electron-pion discrimination.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Surrogate optimization of variational quantum circuits

Variational quantum eigensolvers are touted as a near-term algorithm capable of impacting many applications. However, the potential has not yet been realized, with few claims of quantum advantage and high resource estimates, especially due to the need for optimization in the presence of noise. Finding algorithms and methods to improve convergence is important to accelerate the capabilities of near-term hardware for VQE or more broad applications of hybrid methods in which optimization is required. To this goal, we look to use modern approaches developed in circuit simulations and stochastic classical optimization, which can be combined to form a surrogate optimization approach to quantum circuits. Using an approximate (classical CPU/GPU) state vector simulator as a surrogate model, we efficiently calculate an approximate Hessian, passed as an input for a quantum processing unit or exact circuit simulator. This method will lend itself well to parallelization across quantum processing units. We demonstrate the capabilities of such an approach with and without sampling noise and a proof-of-principle demonstration on a quantum processing unit utilizing 40 qubits.

Gustafson, Erik J. [RIACS, Mtn. View] (ORCID:00000↗

DOE OSTI2024

A compact and portable gamma-ray spectrometer (GRASP) for inertial confinement fusion and basic science experiments

A compact and portable gamma-ray spectrometer has been designed to diagnose different components of the inertial confinement fusion-relevant γ-ray spectrum with energies between ∼3.7–17.9 MeV. The system is designed to be as compact as possible for convenient transportation and fielding in diagnostic ports on the OMEGA laser, the National Ignition Facility, and other photon-source facilities. The system consists of a conversion foil for Compton scattering in front of four magnetic spectrometer “arms,” each covering a different energy range and constructed out of cylindrical permanent magnet Halbach arrays. Monte Carlo simulations have been used to optimize and assess the performance of the conversion foil, and COSY INFINITY ion-optical simulations have been used to optimize the spectrometer magnets. The performance of the design is assessed for a simulated direct-drive γ-ray spectrum. Spanning its total γ-ray energy bandwidth and using a 1.7 mm thick boron conversion foil, the system’s total energy resolution and efficiency are ∼15.8%–4.5% and 5.4 × 10−7–3.7 × 10−7e−/γ, respectively, with room for improvement. Spectral γ-ray measurements will provide guidance to the inertial confinement fusion program toward achieving high-energy gain relevant to inertial fusion energy and enable new measurement capabilities for basic discovery science.

Instruments & Instrumentation↗

DOE OSTI2025

Process intensification approach to enhancing heat and mass transfer during drying: Ultrasonic (US) assisted drying of paper and board

Drying of paper and board is conventionally achieved through alternating conduction (steam-heated cylinders) and pocket convection (heated air over the paper web surface). These conventional drying systems rely heavily on steam from fossil fuels, resulting in inefficiencies, high energy usage, and thermal losses due to surface-driven mechanisms. Here, to address these challenges, an experimental system, with in-situ drying characteristics measurements, was developed to investigate process intensification using ultrasonic-based dewatering—a volumetric, pressure-driven acoustic energy system—integrated with conventional drying. The objectives of this study are to assess the impact of ultrasonics (US) on dewatering; compare performances to conventional drying systems; identify improvements in drying rate and energy use as a function of moisture content; and gain potential insights on heat and mass transfer mechanisms during US-assisted drying. US performance was evaluated across frequencies, power levels, pulp types, and basis weights. Results show that improvements to ultrasonic applications in conjunction with convection were 30-43% in drying rate and 20-35% in drying time over continuous and intermittent applications. When combined with conduction and convection, ultrasonics yielded up to 20% improvement in both rate and time and up to 20% reduction in energy consumption. Observations support a hypothesis of extension of the constant rate period due to improved capillary flow at higher moisture content and enhancing vapor diffusion and boundary layer disruption at lower moisture contents during falling rate period. These findings will inform future modeling, simulation, design and optimization of advanced drying systems.

42 ENGINEERING↗

DOE OSTI2024

Adsorption Thermodynamics for Process Simulation

Adsorption has rapidly evolved in recent decades and is an established separation technology extensively practiced in gas separation industries and others. However, rigorous thermodynamic modeling of multicomponent adsorption equilibrium remains elusive, and industrial practitioners rely heavily on expensive and time-consuming trial-and-error pilot studies to develop adsorption units. Here, this article highlights the need for rigorous adsorption thermodynamic models and the limitations and deficiencies of existing models such as the extended Langmuir isotherm, dual-process Langmuir isotherm, and adsorbed solution theory. It further presents a series of recent advances in the generalization of the classical Langmuir isotherm of single-component adsorption by deriving an activity coefficient model to account for the adsorbed phase adsorbate–adsorbent interactions, substituting adsorbed phase adsorbate and vacant site concentrations with activities, and extending to multicomponent competitive adsorption equilibrium, both monolayer and multilayer. Requiring a minimum set of physically meaningful model parameters, the generalized Langmuir isotherm for monolayer adsorption and the generalized Brunauer–Emmett–Teller isotherm for multilayer adsorption address various thermodynamic modeling challenges including adsorbent surface heterogeneity, isosteric enthalpies of adsorption, BET surface areas, adsorbed phase nonideality, adsorption azeotrope formation, and multilayer adsorption. Also discussed is the importance of quality adsorption data that cover sufficient temperature, pressure, and composition ranges for reliable determination of the model parameters to support adsorption process simulation, design, and optimization.

09 BIOMASS FUELS↗

DOE OSTI2024

Multivariate Testing of Sampling Techniques to Address Class Imbalance in Building Use Type Classification

This study addresses the challenges inherent in building use type classification, particularly focusing on the issue of class imbalance in the training datasets for machine learning classifiers. We comprehensively analyze the efficacy of various class-balancing sampling techniques. Employing Monte Carlo simulations and Bayesian optimization, we evaluated the performance of multiple sampling methods, including Random Oversampling, Random Undersampling, SMOTE, Borderline-SMOTE, and ADASYN, across a dataset encompassing nine southeastern coastal states of the United States. Our findings reveal that simple random over- and undersampling techniques outperform more sophisticated methods. Additionally, we show inherent value in creating an imbalance in training data to effectively train a machine learning classifier for distinguishing between residential and nonresidential buildings. This study provides valuable guidance for future research on building use type classification research and lays essential groundwork for developing attribute-rich building stock datasets.

Adams, Daniel↗

DOE OSTI2024

Samoa Updater: An Application of the Levenberg-Marquardt Method to Update DELFIC Predictions Using Field Measurements

The US Department of Energy (DOE) Forensics Operations (DFO) is a member of the Ground Collections Task Force (GCTF), which is responsible for sample collection of radiological debris for attribution should a nuclear detonation ever occur in the United States. The DFO runs the Defense Land Fallout Interpretive Code (DELFIC) Fallout Planning Tool to predict the deposition of fallout from a nuclear detonation. This prediction is refined using the DELFIC Updater tool, which takes ground measurements and adjusts DELFIC inputs to minimize the difference between prediction and observation, yielding improved predictions of fallout in locations both measured and not yet measured. Samoa, a framework for uncertainty analysis and optimization, is used to improve DELFIC predictive fallout modeling. This new capability using Samoa, dubbed “Samoa Updater,” is compared with the current DELFIC Updater, a brute-force sampling approach. Samoa Updater uses the Levenberg– Marquardt (LM) method, a gradient-based nonlinear least squares approach that uses the functional shape of the input space to increase optimization speed. In simulated test cases Samoa Updater yields faster and more accurate solutions than the current Updater.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2024

MULTI-LEADER: MULTI-source LEarning-Accelerated Design of high-Efficiency multi-stage compRessor (Final Technical Report)

The objective of MULTI-LEADER is to cut design costs by 80% while generating more energy-efficient designs of multi-stage compressors by developing and implementing novel machine learning (ML) techniques, which enable faster and fewer design iterations, improved solver performance, and concurrent multi-disciplinary design. Current industrial practices for the design of multi-stage compressors involve simulation-based design optimization with successive levels of model fidelity, iteratively evaluated between distinct disciplines, one stage at a time to tackle the high dimensional design variations. This project addresses these key design challenges: (1) concurrent optimization of multiple stages under many non-linear constraints; (2) multitude of evaluation of high-fidelity and expensive solvers and their gradients during optimization convergence in high-dimensional design; (3) multi-disciplinary design to maximize aerodynamic performance while guaranteeing structural integrity and additive manufacturability; (4) utilization of multiple fidelity of solvers with disparate parameterization and modeling assumptions. MULTI-LEADER achieved more than 5x speed up in detailed design of more energy-efficient compressors via these machine learning (ML) innovations: (i) rapid design surrogates by multi-source learning from diverse fidelities across multiple disciplines, (ii) physics-constrained data-augmented modeling for improved empiricism, (iii) generative manifold embedding for high dimensional concurrent design without gradient information; (iv) budget-constrained fidelity-adaptive sampling towards fewer design iterations.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Region of Interest Filter in DUNE DAQ

The Deep Underground Neutrino Experiment (DUNE) is the next generation neutrino experiment currently under construction. It consists of a broadband neutrino beam at Fermilab, a high precision near detector, and the largest liquid argon time projection chamber far detector ever designed at the Sanford Underground Research Facility (SURF). The Region of Interest (ROI) filter is designed for DUNE’s online Data Acquisition (DAQ) system to address data rate constraints and enable low energy physics in the <10 MeV range. The filter employs zero suppression on the detector signal, and by tuning the readout window and threshold the data rates can be reduced by >90%. Performance of the ROI is analyzed in LArsoft on MARLEY generated low energy MC events propagated through the detector simulation. Notably, the optimized ROI filter enables a lower trigger threshold for readout at ~5-10 MeV, allowing DUNE to explore low-energy physics, specifically focusing on solar boron 8 neutrinos which are relevant in this energy range. This advancement enhances DUNE's scientific capabilities, opening avenues for detailed analyses of previously inaccessible low-energy neutrino interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Harnessing Quantum Information Science for Enhancing Sensors in Harsh Fossil Energy Environments

The main goals of this project are to utilize real-time quantum dynamics simulations and quantum optimal control algorithms to (1) harness near-surface nitrogen vacancy (NV) centers to detect chemical analytes in harsh fossil energy environments, and (2) design optimally constructed electromagnetic fields for initializing these near-surface NV center spins for efficient sensor performance and detectivity.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2025

Radio Frequency Algorithms

Final Technical Report for DOE SC0019287 award titled Radio Frequency Algorithms. This project supported the design, simulation, analysis, and optimization of robust algorithms to configure, optimize, and control Radio Frequency and Low-Level RF systems for several modern accelerator systems.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Block encoding bosons by signal processing

Block Encoding (BE) is a crucial subroutine in many modern quantum algorithms, including those with near-optimal scaling for simulating quantum many-body systems, which often rely on Quantum Signal Processing (QSP). Currently, the primary methods for constructing BEs are the Linear Combination of Unitaries (LCU) and the sparse oracle approach. In this work, we demonstrate that QSP-based techniques, such as Quantum Singular Value Transformation (QSVT) and Quantum Eigenvalue Transformation for Unitary Matrices (QETU), can themselves be efficiently utilized for BE implementation. Specifically, we present several examples of using QSVT and QETU algorithms, along with their combinations, to block encode Hamiltonians for lattice bosons, an essential ingredient in simulations of high-energy physics. We also introduce a straightforward approach to BE based on the exact implementation of Linear Operators Via Exponentiation and LCU (LOVE-LCU). We find that, while using QSVT for BE results in the best asymptotic gate count scaling with the number of qubits per site, LOVE-LCU outperforms all other methods for operators acting on up to qubits, highlighting the importance of concrete circuit constructions over mere comparisons of asymptotic scalings. Using LOVE-LCU to implement the BE, we simulate the time evolution of single-site and two-site systems in the lattice theory using the Generalized QSP algorithm and compare the gate counts to those required for Trotter simulation.

Kane, Christopher F↗

DOE OSTI2025

Modeling and Analysis of Wind Turbine Wake Vortex Evolution Due to Time-Constant Spatial Variations in Atmospheric Flow

Modern utility-scale wind turbines are evolving toward larger, lighter, and more flexible designs to meet the growing demand for renewable energy while minimizing logistical costs. However, these advancements in lightweight design result in heightened aeroelastic sensitivity, leading to complex interactions which affect the rotor’s capacity to withstand aerodynamic loading and the cascading effects that manifest in the wake’s vortex-structure evolution under variable atmospheric conditions. In this paper, we analyze the influence of stream-wise fluctuating atmospheric flow conditions on wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex–wake interactions for a range of conditions with spatially variant atmospheric flow. CODEF solutions for turbine operation in wind conditions featuring only one parameter fluctuation are compared to wind conditions with several wind parameter variations in combination. By isolating individual inflow variations and comparing them to multi-parameter scenarios, we determine the contributions of each atmospheric factor to rotor dynamics, wake evolution, and downstream wind farm interactions. The purpose of this paper is to analyze the effects of spatial variations in atmospheric flow on the topological evolution of wind turbine vortex wakes, which constitutes a gap in the current understanding of wind turbine wake dynamics. The insights gained from this study are particularly valuable for the development of wind farm control strategies aimed at mitigating the adverse effects of wake interactions, enhancing energy capture, and improving the overall stability of wind farm operations. With these insights, we aim to contribute to the development of modeling and simulation tools to optimize utility-scale wind power plants operating in diverse atmospheric environments.

Farrell, Alayna (ORCID:000000023555720X)↗

DOE OSTI2024

An integrated approach to optimizing concentration shock wave electrodialysis using 2D multicell simulation and response surface models

Shock wave electrodialysis (SWED) is a highly promising technique for energy-efficient ion separation in the context of a circular economy. This paper presents a approach way of modeling and improving SWED using a two-dimensional multicell model combined with the COMSOL program and response surface methodology. The model integrates the Nernst-Planck equation, Darcy's law, and first-order electroosmosis to examine the local concentration, flux of ionic species, distribution of current, and velocity of flow in SWED cells under various operating conditions. We first illustrate the clear depiction of concentration, velocity, and electric potential distribution through contours which aids in identifying optimal operating conditions and designing scalable SWED systems. The results emphasize the significance of surface charge density and voltage in influencing the features of shock waves for obtaining effective ion separation while optimizing energy consumption and improving current efficiency by controlling the retention time of feed flow. Here, this study defines two crucial characteristics of shock waves, namely the length of the flat depletion zone of a fully developed shock wave (shock wave height) and the distance of shock wave propagation (shock wave length). These properties significantly impact separation performance, as determined by the simulation results. Additionally, the response surface methodology is incorporated with the COMSOL models to develop predictive models and graph responses, enabling a more comprehensive understanding of the interactions between parameters and performance indicators, such as removal ratio, energy consumption, and water recovery. Finally, this work suggests design tactics for expanding SWED processes and outlines potential areas for further research. This research provides valuable insights into the prospective applications, design optimization, and scalability of SWED in the field of electrokinetic separation technologies for green chemistry and a circular economy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗