Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation and 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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Convergence and Quantum Advantage of Trotterized MERA for Strongly-Correlated Systems

Strongly-correlated quantum many-body systems are difficult to study and simulate classically. We recently proposed a variational quantum eigensolver (VQE) based on the multiscale entanglement renormalization ansatz (MERA) with tensors constrained to certain Trotter circuits. Here, we determine the scaling of computation costs for various critical spin chains which substantiates a polynomial quantum advantage in comparison to classical MERA simulations based on exact energy gradients or variational Monte Carlo. Algorithmic phase diagrams suggest an even greater separation for higher-dimensional systems. Hence, the Trotterized MERA VQE is a promising route for the efficient investigation of strongly-correlated quantum many-body systems on quantum computers. Furthermore, we show how the convergence can be substantially improved by building up the MERA layer by layer in the initialization stage and by scanning through the phase diagram during optimization. For the Trotter circuits being composed of single-qubit and two-qubit rotations, it is experimentally advantageous to have small rotation angles. We find that the average angle amplitude can be reduced considerably with negligible effect on the energy accuracy. Benchmark simulations suggest that the structure of the Trotter circuits for the TMERA tensors is not decisive; in particular, brick-wall circuits and parallel random-pair circuits yield very similar energy accuracies.

Miao, Qiang [Duke Quantum Center, Duke University,↗
DOE OSTI2024

Performance Results on CPU/GPU Exascale Architectures for OMEGA: The Ocean Model for E3SM Global Applications

The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.

58 GEOSCIENCES↗
DOE OSTI2025

Design and Study of Inductively Coupled Plasma Chamber Components Using the SupRISE Test Device at DIII-D

The DIII-D National Fusion Facility aims to increase the auxiliary heating power for the tokamak by upgrading the Neutral Beam Injection (NBI) system. In collaboration with North Carolina State University, the conventional ‘arc and filament’ NBI ion sources will be converted to inductively coupled plasma (ICP) sources which utilize radio frequency (RF) coupling to maximize reliability for high power operation. In support of this initiative, a full scale test device (SupRISE, Superior Radiofrequency Ion Source Experiment) is currently under construction at the DIII-D Facility. In preparation for the construction of a full scale prototype that can be installed on the DIII-D NBI system, experiments on SupRISE will be conducted to determine the optimal RF frequency for high power coupling, Faraday shield slit configuration, and ICP chamber length. SupRISE is comprised of an approximately 30 x 70 cm quartz dielectric vessel with an internal Faraday shield enclosed in a secondary vacuum chamber to ensure structural stability of the dielectric. Actively cooled front and back plates are designed to reduce the thermal stresses on the plasma facing components and mate with the existing accelerator used by the NBI system at DIII-D. 50 kW of RF power will be coupled to the plasma through the quartz over a variable frequency range of 4-8 MHz to sustain a plasma density of ∼ 10 18 m −3 for a 10 s ON, 210 s OFF duty cycle. Various modelling efforts have been employed to simulate the thermal and stress profiles over the primary components of the SupRISE device as well as the inductance behavior of the RF antenna. These simulation results and the final design for SupRISE will be presented. An additional reduced-scale predecessor ICP source (called RISE) has been used to develop a predictive match model that will be applied to frequency optimization studies on SupRISE. Furthermore, the outcomes of this research and complementary efforts at North Carolina State University are essential for the incorporation of ICP NBI positive ion sources at the DIII-D facility.

DIII-D↗
DOE OSTI2024

Comparative study of machine learning techniques for post-combustion carbon capture systems

Computational analysis of countercurrent flows in packed absorption columns, often used in solvent-based post-combustion carbon capture systems (CCSs), is challenging. Typically, computational fluid dynamics (CFD) approaches are used to simulate the interactions between a solvent, gas, and column's packing geometry while accounting for the thermodynamics, kinetics, heat, and mass transfer effects of the absorption process. These simulations can then be used explain a column's hydrodynamic characteristics and evaluate its CO 2 -capture efficiency. However, these approaches are computationally expensive, making it difficult to evaluate numerous designs and operating conditions to improve efficiency at industrial scales. In this work, we comprehensively explore the application of statistical ML methods, convolutional neural networks (CNNs), and graph neural networks (GNNs) to aid and accelerate the scale-up and design optimization of solvent-based post-combustion CCSs. We apply these methods to CFD datasets of countercurrent flows in absorption columns with structured packings characterized by several geometric parameters. We train models to use these parameters, inlet velocity conditions, and other model-specific representations of the column to estimate key determinants of CO 2 -capture efficiency without having to simulate additional CFD datasets. We also evaluate the impact of different input types on the accuracy and generalizability of each model. We discuss the strengths and limitations of each approach to further elucidate the role of CNNs, GNNs, and other machine learning approaches for CO 2 -capture property prediction and design optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Comparing Control Performance Between Simulation and Experiment using the Microreactor Automated Control System Testbed

In the advanced reactor domain, a flexible and scalable software/hardware infrastructure is crucial for integrating and validating various control technologies. This study used the Microreactor Automated Control System (MACS) hardware platform as a testbed. MACS was originally designed to mirror Idaho National Laboratory (INL)'s Microreactor Applications Research Validation and Evaluation (MARVEL), a 85-kW thermal fission microreactor. It features control drums for simulated reactivity control; lights that function as a surrogate reactor core, with the brightness being proportional to the reactor power; and light sensors that emulate neutron detectors. To transform MACS into a physical twin of MARVEL for evaluating control methods, the Control and Optimization Modular Modeling Application for Nuclear Deployment (COMMAND) software was employed. This software integrated the hardware with two models of the MARVEL core, based on Reactor Excursion and Leak Analysis Program (RELAP5-3D) and Monte Carlo N-Particle (MCNP) models. The study aimed to demonstrate the gap between control theory and actual practice—a gap that often necessitates empirical adjustments such as control gain retuning, filters, time discretization, and integrator anti-windup measures. Controllers were developed based on increasingly complex simulations without hardware, starting from the base MARVEL model and then introducing actuator saturation constraints and sensor noise. The final control strategy was then tested using MACS, and a comparative performance analysis was conducted.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2026

Data driven drift correction for complex optical systems

To exploit the thousand-fold increase in spectral brightness of modern light sources, increasingly intricate experiments are being conducted that demand extremely precise beam trajectory. Maintaining the optimal trajectory over several hours of an experiment with the needed precision necessitates active drift control. Here, we outline time varying Bayesian optimization (TVBO) as a data driven approach for robust drift correction, and illustrate its application for a split and delay optical system composed of six crystals and twelve input dimensions. Using numerical simulations, we exhibit the application of TVBO for linear drift, non-smooth temporal drift as well as constrained TVBO for multi-objective control settings, representing real-life operating conditions. This approach can be easily adapted to other X-ray beam conditioning and guidance systems, including multi-crystal monochromators and grazing-incidence mirrors, to maintain sub-micrometer and nanoradian beam stability over the course of an experiment spanning several hours.

Bayesian optimization↗
DOE OSTI2025

Extending JuTrack’s capabilities to the FRIB accelerator to enhance online modeling

JuTrack is a Julia-based accelerator modeling and tracking package that utilizes compiler-level automatic differentiation (AD) to enable fast and accurate derivative calculations. While JuTrack provides a solid foundation for beam dynamics simulations, its capabilities must be extended to support the Facility for Rare Isotopes (FRIB) linac. This includes modeling heavy-ion linac accelerator components such as the liquid-lithium charge stripper, which facilitates efficient acceleration by remove electrons from heavy isotopes, and incorporating multi-charge state acceleration tracking, which allows for charge-dependent beam dynamics. These extensions address challenges such as the beam matching and optimization of multi charge state through various accelerating structures and beam-material interaction modeling while maintaining the auto differentiation capability. This work focuses on adapting JuTrack to incorporate these elements, enhancing its online modeling abilities. We present modifications to JuTrack’s framework and demonstrate their performance in FRIB simulations.

Accelerator Physics↗
DOE OSTI2026

Non-Invasive Electrode Impedance Estimation for Optimized Charge Profile Parameterization of Lithium-Ion Batteries

This work presents a non-invasive method for parameterizing a physically motivated equivalent circuit model of lithium-ion batteries using operando electrochemical impedance spectroscopy and time-domain data. The proposed model consists exclusively of linear circuit elements, enabling computationally efficient simulation and real-time implementation on battery management system chips. By integrating frequency- and time-domain measurements, the model accurately estimates internal states such as the negative electrode potential, achieving a root mean square error of 12.3 mV during fast charging validation. Parameterization requires only rate tests with sinusoidal perturbations at three different ambient temperatures, making the approach experimentally accessible. The model reveals key insights into battery behavior, including rate-dependent overpotentials primarily governed by charge transfer kinetics at the positive electrode, and temperature-dependent impedance contributions from both charge transfer and solid-state diffusion processes. Validation using reference electrodes confirms the model’s ability to detect lithium plating onset and reproduce impedance behavior across a wide range of operating conditions. The approach enables in situ optimization of fast charging profiles and lays the foundation for future extensions incorporating aging effects and plating dynamics.

25 ENERGY STORAGE↗
DOE OSTI2026

Sensitivity analysis of design parameters in RowWise borehole layout for ground heat exchangers

Ground source heat pump systems offer a promising pathway toward energy-efficient building heating and cooling. The performance and cost-effectiveness of these systems heavily depend on the design of the ground heat exchanger (GHE), particularly the spatial placement of boreholes used in large commercial buildings. Among various borefield layout strategies, the RowWise approach generates and optimizes borehole configurations within irregular polygonal land boundaries, providing land-use efficiency and installation flexibility. While multiple geometric design parameters constrain the RowWise layout, their influence on the system thermal performance and total drilling requirements remains unclear. Thus, this study presents a sensitivity analysis of key design parameters influencing the RowWise layout of vertical borehole GHEs, including the perimeter spacing ratio, borehole spacing, borehole field rotation angle, and borehole length. GHEDesigner and Ray Tune are employed to generate and assess different RowWise configurations. A series of parametric simulations was conducted to quantify the impact of each parameter on geometric distribution, economic cost, and computational speed. The results provide critical insights into the sensitivity and relative importance of different design parameters of RowWise method, offering practical guidance for designers aiming to optimize GHE layouts by balancing thermal efficiency, land constraints, and economic feasibility.

Xu, Dikai [Purdue University]↗
DOE OSTI2024

PV Performance Modeling and Stakeholder Engagement (Final Technical Report)

This core capability project’s objective is to increase the value of photovoltaic (PV) performance models by improving their functionality, demonstrating, and quantifying their validity, and offering a wide range of stakeholder engagement opportunities. In FY22-24, we developed new and improved modeling algorithms and functions to represent PV performance more accurately in a variety of environments and conditions. The “Model parameter toolkit” was developed and includes functions to translate between different module temperature models, incidence angle modifier models, and single-diode models. A new modeling capability named “PV Atlas” was also developed leveraging Sandia’s High Performance Computing resources. This capability allows us to investigate several questions and provide climate-specific best practices and geographic data files; all these are hosted on an interactive website on Sandia’s GitHub and can be used for training, system optimization, or to provide best practices for uncertainty reduction. For model validation, we published high-quality PV performance, and weather data; these data are well documented, filtered, and processed for quality and include examples on how to run PV simulations. We also developed well documented, standardized methods for validating PV models and ran independent model validation and 2 blind modeling intercomparisons engaging with 49 organizations from 17 countries. We co-led and contributed to a growing, well documented and maintained suite of open-source functions for PV modeling (i.e., the pvlib-python) and we outreached to the PV modeling stakeholders via the PVPMC workshops and web resources. In addition, this project supported US representation and leadership for the International Energy Agency (IEA) PVPS Task 13; specifically, members of our team led and supported 3 subtasks on: 1) Best practices for the optimization of bifacial photovoltaic tracking, 2) Extreme weather events and their multiple impact on PV power plants: Risks, failure mechanisms and mitigation strategies, and 3) Best practice guidelines for the use of economic and technical Key Performance Indicators (KPIs). This project resulted in the publications of 14 peer reviewed journal papers, 37 conference presentations, 6 SAND reports, 5 public datasets and 6 new webpages on the PVPMC website. It supported the release of 13 pvlib-python versions where 28 enhancements were from this PV Performance Modeling project. We co-organized 5 PVPMC workshops in FY22-24 with the participation of 214 unique institutions and around 700 participants. The PVPMC website was redesigned, and its reliability was improved; it receives over 50,000 visitors/year from 202 unique countries.

14 SOLAR ENERGY↗
DOE OSTI2024

Water transport through monolayer fullerene membrane

Water transport through nanoporous materials is important in water treatment, desalination, and nanofiltration. Two-dimensional (2D) membranes such as porous graphene have been explored for high-permeance water transport. However, water transport through a new class of 2D membranes based on two-dimensional covalently linked fullerene monolayers has not been fully explored. Here we use classical molecular dynamics simulations to investigate both vapor and liquid water transport through a monolayer fullerene membrane. We find that a quasi-tetragonal phase fullerene membrane possesses the right pore size and geometry that allows fast water vapor transport (∼ 50 g m −2 day −1 Pa −1 ) and water liquid transport (∼ 2.0 g m −2 day −1 Pa −1 ). Furthermore, simulation of sea water transport through the fullerene membrane shows 100 % salt rejection. The much faster vapor transport rate is attributed to the funnel-shaped pore and the optimal size that allows free rotation of water molecules permeating through, while the slower liquid transport is due to the need to desolvate a water molecule to break its hydrogen-bond network across the hydrophobic pore. This work shows the great potential of using monolayer fullerene membranes as 2D membranes for fast and selective water transport.

Desalination↗
DOE OSTI2026

Surface Segregation of Liquid Metal Plasma-Facing Component Alloys: A ReaxFF Investigation

Engineering liquid metal alloys offers a transformative pathway for plasma-facing components by enabling chemically tailored surfaces that can simultaneously optimize plasma-material interactions, reduce divertor heat flux, and enhance core plasma confinement, thereby advancing the commercial viability of nuclear fusion power plants. This study, employing an atomistic simulation approach, provides direct evidence that incorporating nonmetal surface-active agents (such as O and H, or their combination) enables strong surface segregation. This capability makes tin−aluminum (Sn−Al) and tin−lithium (Sn−Li) alloys, with suitable compositions, good candidates for PFC applications. Specifically, the presence of low-Z solutes (Li, Al) leads to preferential surface enrichment, which imparts low-Z sputtering characteristics, while the Sn solvent maintains thermophysical stability. To systematically examine this behavior, we developed a ReaxFF force field spanning the full Sn/Al/Li/O/H chemistry, validated it against formation energies and elastic constants, and applied it in reactive molecular dynamics simulations at fusion-relevant temperatures. We also introduced an overlapbased segregation index that captures interfacial compositional separation directly from atomistic density distributions. Here, this metric reveals a clear hierarchy of segregation regimes and provides a unified view across all systems studied. Together, these findings establish a mechanistic link between nonmetal chemistry and interfacial structure, providing a predictive framework for designing self-adaptive, low-sputtering liquid metal alloys for fusion applications.

Alloys↗
DOE OSTI2025

Quantum Monte Carlo Approaches to Na Intercalation on Bilayer Graphene

We have performed Quantum Monte Carlo (QMC) simulations on Na-intercalated bilayer graphene to study the evolution of electronic and optical properties upon Na intercalation into hard carbon layers. The objective was to model the optimal configuration of Na intercalation into a hard carbon matrix containing graphene regions. Our study showed that Na intercalation can be energetically stabilized at large interlayer distances (over 6 Å) in both AA- and AB-stacked bilayer graphene. In the QMC results, we found a significant band gap opening at the equilibrium interlayer distance of Na-intercalated bilayer graphene, while corresponding density functional theory (DFT) results showed no gap. This difference between DFT and QMC results indicates that the gap opening induced by Na intercalation into a hard carbon is underestimated within the DFT framework. In addition, a zigzag configuration of Na atoms was found to be energetically stable at interlayer distances up to 10 Å, leading us to predict the existence of a local minimum of Na intercalation at large interlayer distance. These computation and modeling results can provide guidance on how to synthesize and optimize hard carbon with bilayer graphene regions that permit a zigzag intercalation configuration that will maximize and stabilize sodium hosting.

Binding energy↗
DOE OSTI2024

Axial confinement of wire array Z-pinch precursor plasmas by a pulsed magnetic mirror field

A novel wire-array Z-pinch platform has been developed to study the effect of a pulsed magnetic mirror field on a collisional high energy density plasma. The mirror is driven in series with the Z-pinch target and, therefore, scales with driver current. Axial confinement is observed based on visible imaging and axial flow velocity measurements. The presence of axial compressing J xB force is determined indirectly based on Thomson scattering and interferometry measurements and corroborated by three-dimensional extended-MHD simulations. Compared to non-magnetized wire array Z-pinch, a modified pulsed mirror configuration is observed to increase precursor plasma thermal energy density by about 30%. If optimized, such a configuration could potentially improve magnetized liner inertial fusion performance by reducing axial plasma end loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Optimal Power Management for Large-Scale Battery Energy Storage Systems via Bayesian Inference

Large-scale battery energy storage systems (BESS) have found ever-increasing use across industry and society to accelerate clean energy transition and improve energy supply reliability and resilience. However, their optimal power management poses significant challenges: the underlying high-dimensional nonlinear nonconvex optimization lacks computational tractability in real-world implementation, and the uncertainty of the exogenous power demand makes exact optimization difficult. This paper presents a new solution framework to address these bottlenecks. The solution pivots on introducing power-sharing ratios to specify each cell’s power quota from the output power demand. To find the optimal power-sharing ratios, we formulate a nonlinear model predictive control (NMPC) problem to achieve power-loss-minimizing BESS operation while complying with safety, cell balancing, and power supply-demand constraints. We then propose a parameterized control policy for the power-sharing ratios, which utilizes only three parameters, to reduce the computational demand in solving the NMPC problem. This policy parameterization allows us to translate the NMPC problem into a Bayesian inference problem for the sake of 1) computational tractability, and 2) overcoming the nonconvexity of the optimization problem. We leverage the ensemble Kalman inversion technique to solve the parameter estimation problem. Concurrently, a low-level control loop is developed to seamlessly integrate our proposed approach with the BESS to ensure practical implementation. This low-level controller receives the optimal power-sharing ratios, generates output power references for the cells, and maintains a balance between power supply and demand despite uncertainty in output power. We conduct extensive simulations and experiments on a 20-cell prototype to validate the proposed approach.

Battery energy storage systems (BESSs)↗
DOE OSTI2025

A continuous calibration of the ATLAS flavour-tagging classifiers via optimal transportation maps

A calibration of the ATLAS flavour-tagging algorithms using a new calibration procedure based on optimal transportation maps is presented. Simultaneous, continuous corrections to the b-jet, c-jet, and light-flavour jet classification probabilities from jet-tagging algorithms in simulation are derived for b-jets using $t\bar{t} \rightarrow e\mu \nu \nu bb$ data. After application of the derived calibration maps, closure between simulation and observation is achieved for jet flavour observables used in ATLAS analyses of Large Hadron Collider (LHC) Run 2 proton-proton collision data. This continuous calibration opens up new possibilities for the future use of jet flavour information in LHC analyses and also serves as a guide for deriving high-dimensional corrections to simulation via transportation maps, an important development for a broad range of inference tasks.

Aad, G. [Aix-Marseille Université] (ORCID:00000002↗
DOE OSTI2024

Training NuGraph2 for ICARUS

This presentation describes the process of training NuGraph2, a Graphical Neural Network for event reconstruction, on simulated ICARUS neutrino event data. This began with an investigation into filtering ICARUS spacepoint data. Then NuGraph2 was repeatedly trained on three event samples, which were used for finding optimized machine-learning parameters and to find and fix the causes of several crashes in NuGraph2 s preprocessing and training scripts.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Surrogate-driven Variance-based Sensitivity Analysis of Thermal Storage Tanks in Integrated Energy Systems

Sensitivity analysis and uncertainty quantification are essential steps for enhancing the accuracy of computational models by identifying and mitigating uncertainties. This study focuses on these steps for the Thermal Energy Delivery System at Idaho National Laboratory, specifically targeting the thermocline tank. Using a Modelica/Dymola simulation model, the study perturbed various design parameters and boundary conditions, including shape factor, porosity, outlet temperature, inlet mass flow rate, and system pressure, to predict and quantify uncertainty in the tank’s ax- ial temperature. A dataset of over 1,000 simulations was generated, and surrogate models were developed using the pyMAISE (Michigan Artificial Intelligence Standard Environment) library, which is an Automatic Machine Learning library for nuclear engineering applications. The optimal model, a feedforward neural network with two hidden layers, achieved an R2 score above 0.99 and a mean absolute error below 1 Kelvin. Sensitivity analyses using Sobol indices and Fourier amplitude sensitivity testing methods on this surrogate model revealed that the inlet mass flow rate at initial timestamps and porosity significantly impacts predicted temperatures across all sensors and time steps.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗