Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical 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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Virtual Tokamak for Test and Development of Plasma Control Applied to NSTX-U

Safe and efficient operation of tokamak experiments depends on shot preparation prior to the experiment to maximize performance and avoid issues that are predictable consequences of known physics. The General Atomics TokSys toolbox is used to test and develop the plasma control system (PCS) on numerous tokamaks around the world. Tokamaks that use a version of the DIII-D PCS can connect it to a TokSys simulation and control a virtual version of the tokamak with the real PCS. Recent upgrades to the simulation include more detailed profile and transport modeling to support design and implementation of optimized scenarios. The simulation combines a module called “Profiles” that simulates the 1D profiles of density, pressure, and current and a module called “GSevolve'' that evolves the 2D Grad-Shafranov equilibrium. The “Profiles'' module primarily uses simplified models for fueling, heating, and current. The TRANSP code provides precomputed diffusion coefficients and heating/current drive profiles from radio-frequency and neutral beam sources for the scenarios being simulated. The fidelity of a simulation can be enhanced by analyzing it with TRANSP and making corrections to the precomputed values and iterating until convergence. Simulations of the National Spherical Torus Upgrade (NSTX-U) are presented. It is shown that the simulation reproduces the experimental profile evolution. Here, the simulation has also been upgraded to reproduce the experimental outcome at points of bifurcation. This is demonstrated with a simulation of vertical displacement in NSTX-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Convergence Analysis of the Alternating Anderson–Picard Method for Nonlinear Fixed-Point Problems

Anderson acceleration (AA) has been widely used to solve nonlinear fixed-point problems due to its rapid convergence. This work focuses on a variant of AA in which multiple Picard iterations are performed between each AA step, referred to as the Alternating Anderson–Picard (AAP) method. Furthermore, despite introducing more “slow” Picard iterations, this method has been shown to be efficient and even more robust in both linear and nonlinear cases. However, there is a lack of theoretical analysis for AAP in the nonlinear case. In this paper, we address this gap by establishing the equivalence between AAP and a multisecant-GMRES method that uses GMRES to solve a multisecant linear system at each iteration. From this perspective, we show that AAP “converges” to the Newton-GMRES method. Specifically, as the residual approaches zero, the multisecant matrix, the approximate Jacobian inverse, the search direction, and the optimization gain of AAP converge to their counterparts in the Newton-GMRES method. These connections provide insights for analyzing the asymptotic convergence properties of AAP. Consequently, we show that AAP is locally 𝑞-linear convergent and provide an upper bound for the convergence factor of AAP. To validate the theoretical results, numerical examples are provided.

Anderson acceleration↗
DOE OSTI2025

A mixed-integer PDE-constrained optimization formulation for constructing electromagnetic cloaks with multiple materials

We study the design of an electromagnetic cloak from multiple materials with an additional constraint on the mass of the cloak. Our problem is an example of a topology optimization problem, and we formulate this problem as a mixed-integer partial-differential equation constrained optimization (MIPDECO) problem, where Maxwell’s equation models the propagation of the wave through the cloak and surrounding medium. We use binary variables to model the assignment of the different materials, and their relevant properties (permittivity and density). The mass constraint adds a nontrivial constraint to this problem. We propose a two-phase strategy to solve this problem. In the first phase, we solve a continuous relaxation, and then propose a new variant of the feasibility pump that exploits the structure of the PDE to obtain an initial integral solution candidate. In the second phase, we use a trust-region approach to improve this incumbent. We also consider a continuation or mesh-sequencing approach to find better solutions faster on consecutively finer meshes. We present detailed numerical results to illustrate the effectiveness of our approaches for constructing multi-material cloaks with a mass constraint.

Calculus of Variations and Optimization↗
DOE OSTI2025

Chapter 4 - Recent Advances in Identification of Differential Equations from Noisy Data: IDENT Review

Differential equations and numerical methods are extensively used to model various real-world phenomena in science and engineering. With modern developments, we aim to find the underlying differential equation from a single observation of time-dependent data. If we assume that the differential equation is a linear combination of various linear and nonlinear differential terms, then the identification problem can be formulated as solving a linear system. The goal then reduces to finding the optimal coefficient vector that best represents the time derivative of the given data. We review some recent works on the identification of differential equations. We find some common themes for the improved accuracy: (i) The formulation of linear system with proper denoising is important, (ii) how to utilize sparsity and model selection to find the correct coefficient support needs careful attention, and (iii) there are ways to improve the coefficient recovery. We present an overview and analysis of recent developments on the topic.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Robust Optimal Control of Inverter-Based Resources Under Grid-Forming Operation

In this paper, we propose and solve a robust control problem for inverter-based resources under grid-forming operation to regulate the voltage and frequency. One major challenge is to mitigate the effect of unmeasurable load current disturbance, grid and load parametric uncertainties. Moreover, strong coupling between the state variables on both the AC and DC sides, as well as between the modulating control input and the frequency impose additional challenges. To address these challenges, first, a robust control problem is solved at the high level via transformation into an equivalent, but more tractable, optimal control problem. Then, in the middle layer a voltage control law is designed on the one side, and a frequency control law on the other side. Finally, an inverter filter current controller is designed to complete the controller design. Theoretical results are derived to provide stability guarantees for the resulting closed-loop system. Specifically, we show that the inverter current injection error is dissipative, the frequency error is semi-globally asymptotically stable, and the inverter terminal voltage error is globally asymptotically stable, all with provided sufficient conditions. Here, numerical simulation experiments are used to validate the theoretical claims. Furthermore, the developed controller is compared with existing work in literature to show the efficacy of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Review of Particle Physics - Scalar Mesons below 1 GeV

The summarizes much of particle physics and cosmology. Using data from previous editions, plus 2,717 new measurements from 869 papers, we list, evaluate, and average measured properties of gauge bosons and the recently discovered Higgs boson, leptons, quarks, mesons, and baryons. We summarize searches for hypothetical particles such as supersymmetric particles, heavy bosons, axions, dark photons, etc. Particle properties and search limits are listed in Summary Tables. We give numerous tables, figures, formulae, and reviews of topics such as Higgs Boson Physics, Supersymmetry, Grand Unified Theories, Neutrino Mixing, Dark Energy, Dark Matter, Cosmology, Particle Detectors, Colliders, Probability and Statistics. Most of the 120 reviews are updated, including many that are heavily revised. The is divided into two volumes. Volume 1 includes the Summary Tables and 97 review articles. Volume 2 consists of the Particle Listings and contains also 23 reviews that address specific aspects of the data presented in the Listings. The complete (both volumes) is published online on the website of the Particle Data Group () and in a journal. Volume 1 is available in print as the . A with the Summary Tables and essential tables, figures, and equations from selected review articles is available in print, as a web version optimized for use on phones, and as an Android app. The 2024 edition of the Review of Particle Physics should be cited as: S. Navas et al. (Particle Data Group), Phys. Rev. D 110, 030001 (2024)© 20242024

Navas, S. [Universidad de Granada]↗
DOE OSTI2025

The Snow Albedo Evolution (SALVO) Campaign at the North Slope of Alaska Field Campaign Report

The springtime surface-albedo transition in the Alaskan Arctic, along with the forces that determine its duration and nature, was the focus of the Snow ALbedo eVOlution (SALVO) campaign. The SALVO team investigated the reasons and durations of the stages of spring melt, during which albedo values decrease from 0.8 to 0.1, signifying the year’s largest and most significant radiative energy change. SALVO II (2022, 2024) built on findings from the successful SALVO I (2019) melt season campaign and previous research on the melt transition period conducted by Grenfell and Perovich (2004). SALVO used the ARM NSA observatory, thus aligning with a core principle of the ARM Decadal Vision to “provide comprehensive and impactful field measurements to support scientific advancement of atmospheric process understanding.” SALVO fieldwork was conducted in the spring near the ARM NSA Central Facility in Utqiagvik, Alaska (Figure 1). Comprehensive data sets were collected in 2019, 2022, and 2024, which included spectral and broadband albedos, snow and meltwater depths, snow stratigraphy, snow grain size data, and aerial imagery. In April of each project year, before the snow began to melt, we established survey lines at three or four locations: inland tundra at NSA E12, coastal tundra at the NSA C1 site, Elson Lagoon, and offshore on the Chukchi Sea sea ice (2022 only). We set up a 200-meter-long line, marked every five meters, to enable repeated measurements of surface albedo and snow depth at each location. During most site visits, the SALVO II team dug at least one snow pit to assess the characteristics of the vertical snow layers, including the types and sizes of snow grains, snow density, and the distribution of liquid water. We evaluated the characteristics of flowing or pooled water at the bottom of the snowpack. Watching snow melt in the Arctic has been compared by some to being as exciting as watching paint dry, but it is far more exciting and dynamic than that. Initially, changes happen slowly while the snow cover remains above 80%, but then reductions in the albedo of the landscape, from about 0.8 to 0.4, combined with rising spring temperatures, begin to accelerate the melt. It happens so quickly that no matter how hard the field team works, they cannot keep up with documenting the changes. At first, liquid water is scarce—only a few wet layers of snow in the snowpack. But then, water becomes ubiquitous, transporting melt energy and creating an extremely heterogeneous albedo and melt landscape. At the end of each field season, the SALVO team is exhausted, relieved to see the snow gone, and joyful to play in the melt ponds on the sea ice under the midnight sun. Overall, the 2024 SALVO field campaign was a highlight, reflecting numerous lessons learned from 2019 and 2022. Daily measurements were taken more than 10 times at each site, resulting in a comprehensive time series of snow conditions and albedo. The team deployed using snowmachines and then on foot when snowmachine travel was no longer permitted on the melting tundra. Instrument mounts and packaging were optimized for quick deployment, regardless of the mode of transportation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A two-dimensional numerical study of the magneto-Rayleigh–Taylor instability with FLASH: Application to the staged Z-pinch concept

Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.

Adaptive mesh refinement↗
DOE OSTI2022

Quantum Circuits for the Preparation of Spin Eigenfunctions on Quantum Computers

The application of quantum algorithms to the study of many-particle quantum systems requires the ability to prepare wave functions that are relevant in the behavior of the system under study. Hamiltonian symmetries are important instruments used to classify relevant many-particle wave functions and to improve the efficiency of numerical simulations. In this work, quantum circuits for the exact and approximate preparation of total spin eigenfunctions on quantum computers are presented. Two different strategies are discussed and compared: exact recursive construction of total spin eigenfunctions based on the addition theorem of angular momentum, and heuristic approximation of total spin eigenfunctions based on the variational optimization of a suitable cost function. The construction of these quantum circuits is illustrated in detail, and the preparation of total spin eigenfunctions is demonstrated on IBM quantum devices, focusing on three- and five-spin systems on graphs with triangle connectivity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

mzPeak: Designing a Scalable, Interoperable, and Future-Ready Mass Spectrometry Data Format

Advances in mass spectrometry (MS) instrumentation, such as higher resolution, faster scan speeds, and improved sensitivity, have significantly increased the volume and complexity of data. The growing adoption of imaging and ion mobility further amplifies these challenges across MS-based omics fields, including proteomics, metabolomics, and lipidomics. While these technologies unlock new possibilities, they also present significant challenges in data management, storage, and accessibility. Existing open formats, such as the XML-based community standards mzML and imzML, struggle to meet the demands of modern MS workflows due to their large file sizes, slow data access, and limited metadata support. Vendor-specific formats, while optimized for proprietary instruments, lack interoperability, comprehensive metadata support and long-term archival reliability. This white paper lays the groundwork for mzPeak, a next-generation community data format designed to address these challenges and support high-throughput, multi-dimensional MS workflows. By adopting a hybrid model that combines efficient binary storage for numerical data and both human and machine-readable metadata storage, mzPeak will reduce file sizes, accelerate data access, and offer a scalable, adaptable solution for evolving MS technologies. For researchers, mzPeak will enable enhanced interoperability across platforms, seamless support for complex workflows including ion mobility and MS imaging, and faster data access compared to existing community formats such as mzML. Its design will ensure data is managed in compliance with regulatory standards, essential for applications such as precision medicine and chemical safety, where long-term data integrity and accessibility are critical. For vendors, mzPeak provides a streamlined, open alternative to proprietary formats, reducing the burden of regulatory compliance while aligning with the industry's push for transparency and standardization. By offering a high-performance, interoperable solution, mzPeak positions vendors to meet customer demands for sustainable data management tools which will be able to handle emerging and future data types and workflows. mzPeak aspires to become the cornerstone of MS data management, empowering researchers, vendors, and developers to innovate and collaborate more effectively.

data formats↗
DOE OSTI2025

The effects of photosynthetic rate on respiration in light, starch/sucrose partitioning, and other metabolic fluxes within photosynthesis

In the future, plants may encounter increased light and elevated CO 2 levels. How consequent alterations in photosynthetic rates will impact fluxes in photosynthetic carbon metabolism remains uncertain. Respiration in light ( R L ) is pivotal in plant carbon balance and a key parameter in photosynthesis models. Understanding the dynamics of photosynthetic metabolism and R L under varying environmental conditions is essential for optimizing plant growth and agricultural productivity. However, measuring R L under high light and high CO 2 (HLHC) conditions poses challenges using traditional gas exchange methods. In this study, we employed isotopically nonstationary metabolic flux analysis (INST-MFA) to estimate RL and investigate photosynthetic carbon flux, unveiling nuanced adjustments in Camelina sativa under HLHC. Despite numerous flux alterations in HLHC, RL remained stable. HLHC affects several factors influencing RL, such as starch and sucrose partitioning, v o /v c ratio, triose phosphate partitioning, and hexose kinase activity. Analysis of A/C i curve operational points reveals that HLHC’s major changes primarily stem from CO 2 suppressing photorespiration. Integration of these fluxes into a simplified model predicts changes in CBC labeling under HLHC. This study extends our prior discovery that incomplete CBC labeling is due to unlabeled carbon reimported during R L , offering insights into manipulating labeling through adjustments in photosynthetic rates.

Elevated CO2↗
DOE OSTI2025

Ionization effects in single-shot carrier-envelope phase detection with gas-gap devices

The carrier-envelope phase (CEP) is a key parameter for attosecond waveform control of ultrashort laser pulses. For laser systems with high repetition rates, however, single-shot CEP detection is still challenging. Building on recent findings on phase detection with electric current generation in gases, we show that with an optimized detection method, a long-term stable single-shot phase detection is feasible. We investigate the achievable performance depending on various parameters such as pulse duration, gas pressure, or incident intensity. The latter exhibits a regime where the saturation of the ionization process boosts the phase sensitivity of the signal, which is elucidated by numerical calculations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Beam breakup instability studies of powerful energy recovery linac for experiments

The maximum achievable beam current in an energy recovery linac (ERL) is often constrained by beam breakup (BBU) instability. Our previous research highlighted that filling patterns have a substantial impact on BBU instabilities in multipass ERLs. In this study, we extend our investigation to the eight-cavity model of the Powerful ERL for Experiment (PERLE). We evaluate its requirements for damping cavity higher order modes (HOMs) and propose optimal filling patterns and bunch timing strategies. Our findings reveal a significant new insight: while filling patterns are crucial, the timing of bunches also plays a critical role in mitigating HOM beam loading and BBU instability. This previously underestimated factor is essential for effective BBU control. We estimated the PERLE threshold current using both analytical and numerical models, incorporating the designed PERLE HOM dampers. During manufacturing, HOM frequencies are expected to vary slightly. Our study found no significant difference in BBU suppression for relative rms frequency jitters of 0.001, 0.002, and 0.005 for the same HOM. Introducing a jitter of 0.001 into our models, we found that the dampers effectively suppressed BBU instability, achieving a threshold current an order of magnitude higher than the design requirement. Our results offer new insights into ERL BBU beam dynamics and have important implications for the design of future ERLs. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Qutrit Circuits and Algebraic Relations: A Pathway to Efficient Spin-1 Hamiltonian Simulation

Quantum information processing has witnessed significant advancements through the application of qubit- based techniques within universal gate sets. Recently, exploration beyond the qubit paradigm to d-dimensional quantum units or qudits has opened new avenues for improving computational efficiency. This paper delves into the qudit-based approach, particularly addressing the challenges presented in the high-fidelity implementation of qudit-based circuits due to increased complexity. As an innovative approach towards enhancing qudit circuit fidelity, we explore algebraic relations, such as the Yang-Baxter-like turnover equation, which may enable circuit compression and optimization. The paper introduces the turnover relation for the three-qutrit time propagator and its potential use in reducing circuit depth. We further investigate whether this relation can be generalized for higher-dimensional quantum circuits, including a focused study on the one-dimensional spin-1 Heisenberg model. Our paper outlines both rigorous and numerically efficient approaches to potentially achieve this generalization, providing a foundation for further explorations in the field of qudit-based quantum computing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modeling Strong Light-Matter Coupling in Correlated Systems: State-Averaged Cavity Quantum Electrodynamics Complete Active Space Self-Consistent Field Theory

The description of strongly correlated systems interacting with quantized cavity modes poses significant theoretical challenges due to the combinatorial scaling of electronic and photonic degrees of freedom. Recent advances addressing this complexity include cavity quantum electrodynamics (QED) generalizations of complete active space configuration interaction and density matrix renormalization group methods. In this work, we introduce a QED extension of state-averaged complete active space self-consistent field theory, which incorporates cavity-induced correlations through a second-order orbital optimization framework with robust convergence properties. The method is implemented using both photon number state and coherent state representations, with the latter showing robust origin invariance in the energies regardless of the completeness of the photonic Fock space. The implementation enables symmetry-free orbital relaxations to account for photon-mediated symmetry breaking in polaritonic systems. Numerical validation on lithium hydride, hydroxide anion, and magnesium hydride cation demonstrates that this method achieves significantly improved accuracy in modeling ground-state and polariton potential energy surfaces compared to QED-CASCI in a fixed orbital basis. In these studies, we reach sub-kcal/mol accuracy in potential energy surface in much smaller active spaces than are required for QED-CASCI. This advancement provides a more robust approach for studying cavity-altered chemical landscapes for ground and exited strongly coupled systems.

CASSCF↗
DOE OSTI2025

Learning error distribution kernel‐enhanced neural network methodology for multi‐intersection signal control optimization

Traffic congestion has substantially induced significant mobility and energy inefficiency. Many research challenges are identified in traffic signal control and management associated with artificial intelligence (AI)-based models. For example, developing AI-driven dynamic traffic system models that accurately capture high-resolution traffic attributes and formulate robust control algorithms for traffic signal optimization is difficult. Additionally, uncertainties in traffic system modeling and control processes can further complicate traffic signal system controllability. To partially address these challenges, this study presents a novel, hybrid neural network model enhanced with a probability density function kernel shaping technique to formulate traffic system dynamics better and improve comprehensive traffic network modeling and control. The numerical experimental tests were conducted, and the results demonstrate that the proposed control approach outperforms the baseline control strategies and reduces overall average delays by 11.64% on average. By leveraging the capabilities of this innovative model, this study aims to address major challenges related to traffic congestion and energy inefficiency toward more effective and adaptable AI-based traffic control systems.

Wang, Hong [Oak Ridge National Laboratory (ORNL), ↗
DOE OSTI2024

TEAMER: Crossflow Turbine Fairing Geometry Optimization - Report and CFD Modeling Files

The dataset includes computational fluid dynamics (CFD) models and simulation files for crossflow turbines as well as a detailed project report. The report documents the project undertaken by the Ocean Renewable Power Company (ORPC) to design and optimize a modular fairing for the Modular RivGen Marine Hydrokinetic (MHK) turbine, which enhances the efficient deployment and operation of turbine arrays. The project focused on optimizing the hydrodynamic performance of the fairing using CFD, with an emphasis on two key geometric parameters: the fairing's cross-sectional shape and the spacing between the rotor and the fairing. The analysis aimed to maximize net power output while also assessing discretized loading to evaluate ultimate and fatigue loads on the turbine components. The numerical modeling was conducted using both the commercial CFD software Star-CCM+ and the open-source code openFOAM, with the latter utilizing the actuator line library, turbinesFOAM. This dual-code approach was intended to increase confidence in the results and demonstrate the viability of using open-source tools for high-fidelity marine energy modeling. This dataset includes all necessary files for actuator line simulations in openFOAM, as well as 2D blade-resolved CFD results, along with Python and Java scripts for setting up and post-processing simulations.

16 TIDAL AND WAVE POWER↗