Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Higher-order 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Higher-order factorization machine for accurate surrogate modeling in material design

Efficient and robust optimization is important in material science for identifying optimal structural parameters and enhancing material performance. Surrogate-based active learning algorithms have recently gained great attention for their ability to efficiently navigate large, high-dimensional design spaces. Among surrogate models, 2 nd -order factorization machine (FM) models are widely employed as the surrogate model in active learning algorithms due to their balance between simplicity and effectiveness. However, their quadratic nature limits their capacity to capture complex, higher-order interactions among variables, often leading to suboptimal solutions. To overcome this limitation, we propose an active learning scheme integrating a 3 rd -order FM model, capable of modeling three-variable interactions and more intricate relationships in material systems. We comprehensively evaluate the surrogate modeling performance of the 3 rd -order FM case using various objective functions. Furthermore, we examine the optimization reliability and efficiency of the 3 rd -order FM-based active learning in a real-world material design task (e.g., nanophotonic structures for transparent radiative cooling). Our study shows that the 3 rd -order FM outperforms the 2 nd -order model in both surrogate accuracy and optimization performance, highlighting higher-order models’ promises for material design and optimization problems.

Factorization machine↗
DOE OSTI2024

Scaling whole-chip QAOA for higher-order ising spin glass models on heavy-hex graphs

Abstract We show that the quantum approximate optimization algorithm (QAOA) for higher-order, random coefficient, heavy-hex compatible spin glass Ising models has strong parameter concentration across problem sizes from 16 up to 127 qubits for p = 1 up to p = 5, which allows for computationally efficient parameter transfer of QAOA angles. Matrix product state (MPS) simulation is used to compute noise-free QAOA performance. Hardware-compatible short-depth QAOA circuits are executed on ensembles of 100 higher-order Ising models on noisy IBM quantum superconducting processors with 16, 27, and 127 qubits using QAOA angles learned from a single 16-qubit instance using the JuliQAOA tool. We show that the best quantum processors find lower energy solutions up to p = 2 or p = 3, and find mean energies that are about a factor of two off from the noise-free distribution. We show that p = 1 QAOA energy landscapes remain very similar as the problem size increases using NISQ hardware gridsearches with up to a 414 qubit processor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Robust Dark Energy Constraints with the Dark Energy Spectroscopic Survey (Final Technical Report)

This project developed and applied advanced theoretical, computational, and data-analysis methodologies to extract robust and precise cosmological constraints from the Dark Energy Spectroscopic Instrument (DESI). The work focused on maximizing the scientific return of DESI through optimized survey strategy, novel higher-order clustering statistics, improved modeling of small-scale structure, and rigorous mitigation of observational systematics. Over the award period, the project made substantial contributions to DESI science planning, produced new methods for bispectrum and three-point correlation function analyses, advanced constraints on primordial non-Gaussianity, and delivered widely used software tools. The project also played a major role in training graduate students and a postdoctoral researcher who contributed directly to DESI key projects. The results have significantly enhanced the cosmological reach of DESI and provide a strong foundation for future surveys such as DESI-II and Stage-V experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Lyapunov controlled counterdiabatic quantum optimization

We introduce a quantum algorithm that integrates counterdiabatic (CD) protocols with quantum Lyapunov control (QLC) to address combinatorial optimization problems. This approach offers versatility, allowing implementation as either a digital-analog or purely digital algorithm based on selected control strategies. By examining spin-glass Hamiltonians, we illustrate how the algorithm can explore alternative paths to enhance solution outcomes compared to conventional CD techniques. This method reduces dependence on extensive higher-order CD terms and on classical optimization techniques, making it more suitable for existing quantum computing platforms. The combination of digital compression via CD protocols and the adaptable nature of QLC methods positions this approach as a promising candidate for near-term quantum devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Gradient-informed Hamiltonian Monte Carlo for multicomponent CALPHAD model optimization and uncertainty quantification

CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Higher-order interaction effects among operating conditions and feedstocks shape reactor microbiomes and fatty acid production profiles

Arrested anaerobic digestion (AAD) offers a promising route for producing fatty acids (FAs) from organic residues, yet optimal conditions for selectively generating medium-chain fatty acids (MCFAs) remain poorly defined. Here, we systematically evaluated the main and interaction effects of pH (5, 7, 9), feedstock (food waste, manure), temperature (35 and 45 °C), and inoculum source on microbiome composition and FA production. Anaerobic digester sludge and a novel bison rumen inoculum were compared. Significant higher-order interactions among operating parameters governed FA profiles and microbiome structure. Butyric acid production was driven by a three-way interaction among pH, feedstock, and temperature (p < 0.001), with maximum concentrations achieved in food waste reactors at pH 5.0 and 35 °C (1.2 ± 0.1 g L −1 with sludge and 1.1 ± 0.3 g L −1 with rumen). MCFA production exhibited significant four-way interactions (p < 0.1 to p < 0.001). At 45 °C and pH 5.0, inoculum source tuned MCFA selectivity: sludge favored pentanoic acid (0.4 ± 0.1 g L −1 ), whereas rumen favored hexanoic and heptanoic acids (up to 0.4 ± 0.2 g L −1 ). Manure reactors produced < 0.2 g L −1 MCFAs under all conditions. Genera, including Megasphaera, Prevotella, and Lactobacillus, were associated with production of specific MCFAs. PICRUSt2-based pathway predictions were consistent with MCFA production patterns and suggested a potential role for lactic acid–driven chain elongation pathways. This study provides insights into how interacting operating conditions shape AAD microbiomes, their FA profiles, and advances the trajectory of research aimed at engineering robust and controllable microbiomes for waste valorization.

09 BIOMASS FUELS↗
DOE OSTI2026

Faster Randomized Dynamical Decoupling

We present a randomized dynamical decoupling (DD) protocol that can substantially improve the performance of any given deterministic DD scheme for suppressing coherent noise by using no more than two additional pulses. Our construction is implemented by probabilistically applying sequences of pulses, which, when combined, effectively eliminate the error terms that scale linearly with the system-environment coupling strength. As a result, we show that a randomized protocol using a few pulses can outperform deterministic DD protocols that require considerably more pulses. Furthermore, we prove that the randomized protocol provides an improvement compared to deterministic DD sequences that aim to reduce the error in the system’s Hilbert space, such as Uhrig DD, which had been previously regarded to be optimal. To rigorously evaluate the performance, we introduce new analytical methods suitable for analyzing higher-order DD protocols that might be of independent interest. Here, we also present numerical simulations confirming the significant advantage of using randomized protocols compared to widely used deterministic protocols.

Quantum algorithms & computation↗
DOE OSTI2025

SymProp: Scaling Sparse Symmetric Tucker Decomposition via Symmetry Propagation

Sparse symmetric tensors are an important class of tensors, and their decompositions serve as powerful tools for revealing low-rank structures. This paper introduces SymProp, a novel approach for scaling sparse symmetric Tucker decomposition by propagating symmetry through intermediate computations. SymProp optimizes two key computational kernels: Sparse Symmetric Tensor Times Same Matrix chain (S3 TTMc) for Higher-Order Orthogonal Iteration (HOOI) and Sparse Symmetric Tensor Times Same Matrix chain Times Core (S3 TTMcTC) for Higher-Order QR Iteration (HOQRI). Our method employs a metaprogramming-based index iteration approach to efficiently handle the upper triangular parts of intermediate dense symmetric tensors. SymProp achieves up to 50.9× speedup over SPLATT and up to 360.8× over Compressed Sparse Symmetric (CSS) format on the S3 TTMc operation. Moreover, our S3 TTMc and S3 TTMcTC implementations support tensor orders four levels higher than state-of-the-art methods. Our HOQRI demonstrates superior scalability and up to a 33.6× speedup over optimized HOOI. By enabling more scalable Tucker decompositions for higher orders, decomposition ranks, and dimension sizes, SymProp opens new possibilities for analyzing complex hypergraph structures in fields such as network science, data mining, and machine learning.

Li, Zecheng [North Carolina State University]↗
DOE OSTI2024

Improving ADAM through an implicit-explicit (IMEX) time-stepping approach

The ADAM optimizer, often used in machine learning for neural network training, corresponds to an underlying ordinary differential equation (ODE) in the limit of very small learning rates. Here, this work shows that the classical ADAM algorithm is a first-order implicit-explicit (IMEX) Euler discretization of the underlying ODE. Employing the time discretization point of view, we propose new extensions of the ADAM scheme obtained by using higher-order IMEX methods to solve the ODE. Based on this approach, we derive a new optimization algorithm for neural network training that performs better than classical ADAM on several regression and classification problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Rapid Tuning of Synchrotron Surrogate Model at the Recycler Ring

The 8 GeV proton-storage Recycler Ring (RR) is essential for reaching megawatt beam intensity goals for the DUNE neutrino beam at Fermilab. Custom shims on each RR permanent magnet were designed to cancel manufacturing defects and bring magnetic fields to the design values. Remaining imperfections cause the observed tune variation vs energy to deviate from what is calculated using the design fields. Using the POUNDERS (“Practical Optimization Using No Derivatives for sums of Squares”) optimization method with Synergia in the loop, we demonstrate rapid convergence to a set of additive, higher-order multipole moments of these magnetic shims which reproduce that observed variation, and show that the convergence advantage grows with the parameter-space dimensionality.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Plasma processing of SRF cavities at Jefferson Lab: Experiment results and simulation insight

Plasma processing of superconducting radio frequency (SRF) cavities has been an active research effort at Jefferson Lab (JLab) since 2019, aimed at enhancing cavity performance by removing hydrocarbon contaminants and reducing field emission. In this experiment, processing using argon-oxygen and helium-oxygen gas mixtures to find minimum ignition power at different cavity pressure was investigated. Ongoing simulations are contributing to a better understanding of the plasma surface interactions and the fundamental physics behind the process. These simulations, combined with experimental studies, guide the optimization of key parameters such as gas type, RF power, and pressure to ignite plasma using selected higher-order mode (HOM) frequencies. This paper presents experimental data from argon-oxygen and helium-oxygen gas mixture C75 and C100 cavity plasma ignition studies, as well as simulation results for the C100-type cavity based on the COMSOL model previously applied to the C75 cavity.

Accelerator Physics↗
DOE OSTI2026

In-Band Scattering and Absorption of Infrared Blocking Foam Filters for Millimeter-wave Cameras

Expanded closed-cell polymer foams are widely used as thermal infrared (IR) blocking filters in millimeter-wave cameras, particularly for Cosmic Microwave Background observations. Precise knowledge of their millimeter-wave properties is essential for optimizing sensitivity. We present broadband (150 GHz - 2 THz) transmittance spectroscopy of Styroace-II and several Zotefoam filters, fitting their spectra with a radiative transfer model incorporating dielectric absorption and Rayleigh, Mie, and higher-order scattering. For a typical 5~cm thick filter stack at 280~GHz, Styroace-II exhibits ${\sim}10\%$ scattering with absorption estimated as ${\lesssim}5\%$ by effective-medium theory, while Zotefoam HD30 offers superior performance at ${\sim}3\%$ scattering and absorption likewise bounded to ${{\lesssim}0.3\%}$. Each model component is constrained at the ${\sim}0.1\%$ transmittance level for millimeter wavelengths. We observe batch-to-batch scattering variability of up to 2 percentage points in foams with multiple tested batches. Less commonly used Zotefoam formulations (LD15 and LD24) can further reduce in-band scattering to ${<}1\%$ while maintaining negligible in-band absorption and likely comparable IR blocking due to shared polyethylene absorption features and similar cell sizes. Based on this work, a filter constructed from the best measured LD24 batch has replaced the Styroace-II filter in a Simons Observatory 220/280 GHz Small Aperture Telescope.

Thomas, Alex [Chicago U., Astron. Astrophys. Ctr.;↗
DOE OSTI2024

Probing inhomogeneous cuprate superconductivity by terahertz Josephson echo spectroscopy

Abstract Inhomogeneities crucially influence the properties of quantum materials, yet methods that can measure them remain limited and can access only a fraction of relevant observables. For example, local probes such as scanning tunnelling microscopy have documented that the electronic properties of cuprate superconductors are inhomogeneous over nanometre length scales. However, complementary techniques that can resolve higher-order correlations are needed to elucidate the nature of these inhomogeneities. Furthermore, local tunnelling probes are often effective only far below the critical temperature. Here we develop a two-dimensional terahertz spectroscopy method to measure Josephson plasmon echoes from an interlayer superconducting tunnelling resonance in a near-optimally doped cuprate. The technique allows us to study the multidimensional optical response of the interlayer Josephson coupling in the material and disentangle intrinsic lifetime broadening from extrinsic inhomogeneous broadening for interlayer superconducting tunnelling. We find that inhomogeneous broadening persists up to a substantial fraction of the critical temperature, above which this is overcome by the thermally increased lifetime broadening.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Half-Life and Precision Shape Measurement of the 2⁢𝜈⁢𝛽⁢𝛽 Decay of 130 Te

Here, we present a new measurement of the 2⁢𝜈⁢𝛽⁢𝛽 half-life of 130 Te (𝑇$^{2⁢𝜈}_{1/2}$) using the first complete model of the CUORE data, based on 1038 kg yr of collected exposure. Thanks to optimized data selection, we achieve a factor of two improvement in precision, obtaining 𝑇$^{2⁢𝜈}_{1/2}$ = (9.32⁢$^{+0.05}_{−0.04}$⁢stat⁢ $^{+0.07}_{−0.07}$⁢syst)×10 20 yr. The signal-to-background ratio is increased by 70% compared to our previous results, enabling the first application of the improved 2⁢𝜈⁢𝛽⁢𝛽 formalism to 130 Te . Within this framework, we determine a credibility interval for the effective axial coupling in the nuclear medium as a function of nuclear matrix elements. We also extract values for the higher-order nuclear matrix element ratios: second-to-first and third-to-first. The second-to-first ratio agrees with nuclear model predictions, while the third-to-first ratio deviates from theoretical expectations. These findings provide essential tests of nuclear models and key inputs for future 0⁢𝜈⁢𝛽⁢𝛽 searches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Deriving cloud droplet number concentration from surface-based remote sensors with an emphasis on lidar measurements

Abstract. Given the importance of constraining cloud droplet number concentrations (Nd) in low-level clouds, we explore two methods for retrieving Nd from surface-based remote sensing that emphasize the information content in lidar measurements. Because Nd is the zeroth moment of the droplet size distribution (DSD), and all remote sensing approaches respond to DSD moments that are at least 2 orders of magnitude greater than the zeroth moment, deriving Nd from remote sensing measurements has significant uncertainty. At minimum, such algorithms require the extrapolation of information from two other measurements that respond to different moments of the DSD. Lidar, for instance, is sensitive to the second moment (cross-sectional area) of the DSD, while other measures from microwave sensors respond to higher-order moments. We develop methods using a simple lidar forward model that demonstrates that the depth to the maximum in lidar-attenuated backscatter (Rmax⁡) is strongly sensitive to Nd when some measure of the liquid water content vertical profile is given or assumed. Knowledge of Rmax⁡ to within 5 m can constrain Nd to within several tens of percent. However, operational lidar networks provide vertical resolutions of > 15 m, making a direct calculation of Nd from Rmax⁡ very uncertain. Therefore, we develop a Bayesian optimal estimation algorithm that brings additional information to the inversion such as lidar-derived extinction and radar reflectivity near the cloud top. This statistical approach provides reasonable characterizations of Nd and effective radius (re) to within approximately a factor of 2 and 30 %, respectively. By comparing surface-derived cloud properties with MODIS satellite and aircraft data collected during the MARCUS and CAPRICORN II campaigns, we demonstrate the utility of the methodology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

SDA: a symbolic differential algebra package in C++

Truncated Power Series Algebra (TPSA), or Differential Algebra (DA), is a well-established tool in accelerator physics, commonly used for generating high-order maps of dynamic systems, as well as in symplectic tracking, normal form analysis, verified integration, optimization, and fast multipole methods. This package is the first to perform symbolic DA computations, enabling traceability of initial condition contributions and runtime reduction for repeated DA calculations, potentially expanding DA’s applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Efficient near-field ptychography reconstruction using the Hessian operator

X-ray ptychography is a powerful and robust coherent imaging method providing access to the complex object and probe (illumination). Ptychography reconstruction is typically performed using first-order methods due to their computational efficiency. Higher-order methods, while potentially more accurate, are often prohibitively expensive in terms of computation. In this study, we present a mathematical framework for reconstruction using second-order information derived from an efficient computation of the bilinear Hessian and Hessian operator. The formulation is provided for Gaussian-based models, enabling the simultaneous reconstruction of the object, probe, and object positions. Synthetic data tests, along with experimental near-field ptychography data processing, demonstrate a ten-fold reduction in computation time compared to first-order methods. The derived formulas for computing the Hessians, along with the strategies for incorporating them into optimization schemes, are well-structured and easily adaptable to various ptychography problem formulations.

Carlsson, Marcus [Lund Univ. (Sweden)] (ORCID:0000↗