Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Combining Generative Modeling and Advanced Control for Building Scenario Generation

Buildings make up a large portion of energy consumption in the U.S. today. Understanding their energy consumption patterns can improve their efficiency, but requires detailed models that rely on incomplete or unknown information. Previous work has shown that artificial intelligence (AI) can be used to predict missing information and even suggest upgrades to improve building efficiency. However, building upgrades may require undesirable upfront costs. Oppositely, advanced control could improve building efficiency with negligible upfront cost. To explore the tradeoffs between these two approaches, in this work we propose a workflow to compute optimal temperature setpoint schedules to minimize energy consumption and operational cost. Results show that modifying the temperature setpoints in a building using model predictive control (MPC) can effectively reduce its energy consumption and operational cost. This optimal operation cannot fully meet a desired goal. However, we show that by considering MPC in addition to component upgrades, a desired goal can be met with significantly less upfront costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Towards Scaling Law Analysis For Spatiotemporal Weather Data

Compute-optimal scaling laws are relatively well studied for NLP and CV, where objectives are typically single-step and targets are comparatively homogeneous. Weather forecasting is harder to characterize in the same framework: autoregressive rollouts compound errors over long horizons, outputs couple many physical channels with disparate scales and predictability, and globally pooled test metrics can disagree sharply with per-channel, late-lead behavior implied by short-horizon training. We extend neural scaling analysis for autoregressive weather forecasting from single-step training loss to long rollouts and per-channel metrics. We quantify (1) how prediction error is distributed across channels and how its growth rate evolves with forecast horizon, (2) if power law scaling holds for test error, relative to rollout length when error is pooled globally, and (3) how that fit varies jointly with horizon and channel for parameter, data, and compute-based scaling axes. We find strong cross-channel and cross-horizon heterogeneity: pooled scaling can look favorable while many channels degrade at late leads. We discuss implications for weighted objectives, horizon-aware curricula, and resource allocation across outputs.

Kiefer Jr, Alexander [ORNL] (ORCID:000000025398874↗
DOE OSTI2023

A new efficient grain growth model using a random Gaussian-sampled mode filter

This paper presents the use of a Gaussian neighborhood mode filter for predicting grain growth in a manner similar to the solutions obtained by a Monte Carlo Potts model. This flexible grain growth model can quickly utilize modern, computationally optimized data science strategies on graphics processing units to simulate grain growth up to 100 times faster than the state-of-the-art, publicly available Monte Carlo Potts model. We show that, given the correct neighborhood, the mode filter can replicate normal grain growth in two or three dimensions. In addition, the paper briefly demonstrates the ability to model limited anisotropic in grain boundary energy and mobility. Anisotropic grain boundary energy is modeled by defining a weighted mode filter operation. Anisotropic grain boundary mobility is modeled by scaling and orienting the Gaussian neighborhood in a particular direction.

Anisotropy↗
DOE OSTI2021

Analysis of Hydrocyclone Geometry via Rapid Optimization Based on Computational Fluid Dynamics

Abstract Hydrocyclones exploit density gradients for the centrifugal separation of dispersions in a continuous liquid. Selection of the geometrics for optimal separation is case specific, like the media characteristics. The existing optimization method based on computational fluid dynamics (CFD) provides a powerful analytical tool but requires long computational times. The most common praxis for CFD optimization is via the single‐factor optimization method (SFOM). In this study, a novel approach is presented as an improved rapid optimization method that implements a dynamic‐mesh and user‐defined function optimization method (DUOM). The DUOM adapts the dynamic‐mesh approach from other applications to the optimization analysis of hydrocyclones. The DUOM reduced the computational time by 31.1 %, compared to the SFOM.

Liu, Lin↗
DOE OSTI2022

Pyomo.DOE: An open-source package for model-based design of experiments in Python

Predictive mathematical models are a cornerstone of science and engineering. Yet selecting, calibrating, and validating said science-based models often remains an art in practice. Model-based design of experiments (MBDoE) provides a systematic framework to maximize information gain from experiments while minimizing time and resource costs. But MBDoE remains limited to niche application areas, in part because practitioners must integrate expertise in statistics, computational optimization, and modeling. To help reduce this barrier, we introduce Pyomo.DOE, an open-source package for MBDoE. Pyomo.DOE uses a nonlinear sensitivity analysis code k_aug to quickly approximate the Fisher information matrix and leverages a new stochastic programming abstraction. We demonstrate Pyomo.DOE with the first application of MBDoE to fixed-bed breakthrough experiments, which highlights the power of Pyomo.DOE to quantify the value of experimental modifications a priori for large-scale partial differential-algebraic equation (PDAE) models. Here we also provide a mathematical primer on MBDoE targeted at general chemical engineers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Coaxial color channel focus evaluation to estimate standoff height in directed energy deposition additive manufacturing

Directed energy deposition (DED) is an additive manufacturing process that is being rapidly adopted by industry and is well suited for the fabrication of complex components in a variety of metal alloys. In laser cladding systems such as DED, powder is blown in a stream to a metal substrate coincident with a laser necessary to deposit molten metal with 3D spatial control. The focus of both the laser and the powder stream are crucial, and best deposition occurs at a predetermined standoff height between the build surface and the print head. Generally, no monitoring of this distance is implemented in commercial DED systems. Due to potential over or under building, the standoff height often changes over time but tends to self-correct. However, inexpensive and minimally intrusive methods to identify optimal standoff are required to provide real-time control to maintain the optimal distance. The present work explores the quantification of the focus of the three-color channels of a coaxial camera to determine the standoff height. An experiment was performed in which a 254 mm wall is built and the standoff height, initially 5.0 mm below the optimal position, was then intentionally increased every 25.4 mm of wall length by an amount of 1.0 mm to a final position 7.0 mm above optimal. Computer vision is demonstrated to monitor the amount of focus in each color band and estimate standoff distance. Finally, a response can be calculated in under 40 ms using simple hardware and can work in most laser-based DED systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Computational Framework to design 3D stiffness gradient acoustic metamaterials for impedance matching

Acoustic waves play a crucial role in various applications, including medical imaging, non-destructive testing, and sonar systems. One of the significant challenges in these applications is impedance matching, which is essential for minimizing reflections and maximizing the transfer of acoustic energy between different media. Acoustic metamaterials offer a promising solution to this challenge. In addition to impedance control, gradient stiffness can enhance structural efficiency and enable spatial control of wave propagation, making it a valuable feature in acoustic metamaterial design. In this pa- per, we present our developed computational method to design 3D stiffness gradient acoustic metamaterials for impedance matching. The key steps in our approach include generating initial designs using a periodic covariance function to provide unit cells that are both periodic on the boundaries and randomly formed inside the unit cell. Furthermore, we integrated manufacturing constraints into the design process, ensuring that the structures are interconnected for fabrication. We propose two computational optimization algorithms: GenUnit, based on a non-dominated sorting genetic algorithm (NSGA-II), and MLMatch, which leverages differentiable machine learning. The two approaches are not separate contributions but complementary com- ponents of a unified framework. GenUnit requires no training data and directly interfaces with physics-based simulations, making it highly accurate but slower for large-scale exploration. In contrast, MLMatch is data-hungry during training but, once trained, enables near-instantaneous inference and broad design-space coverage. Together, they form a hybrid strategy: ML- Match rapidly explores the global design space, and GenUnit provides local refinement with high-fidelity accuracy. This balance between training cost, inference time, and precision is the motivation for including both methods in the same study. We applied this dual-algorithm framework to generate two metallic-based metamaterial designs that match the acoustic impedance of water while exhibiting a controlled gradient in stiffness (from stiff to soft). The stiffness gradient is particularly advantageous in applications where one side of the structure must interface with soft or sensitive surfaces, such as human tissue or delicate components. Here, this work paves the way for improved materials in various acoustic applications, particularly in ultrasound devices, by providing better impedance.

Metamaterial↗
DOE OSTI2024

Dyn$\mathrm{AMO}$: Multi-agent reinforcement learning for dynamic anticipatory mesh optimization with applications to hyperbolic conservation laws

Here we introduce DynAMO, a reinforcement learning paradigm for Dynamic Anticipatory Mesh Optimization. Adaptive mesh refinement is an effective tool for optimizing computational cost and solution accuracy in numerical methods for partial differential equations. However, traditional adaptive mesh refinement approaches for time-dependent problems typically rely only on instantaneous error indicators to guide adaptivity. As a result, standard strategies often require frequent remeshing to maintain accuracy. In the DynAMO approach, multi-agent reinforcement learning is used to discover new local refinement policies that can anticipate and respond to future solution states by producing meshes that deliver more accurate solutions for longer time intervals. By applying DynAMO to discontinuous Galerkin methods for the linear advection and compressible Euler equations in two dimensions, we demonstrate that this new mesh refinement paradigm can outperform conventional threshold-based strategies while also generalizing to different mesh sizes, remeshing and simulation times, and initial conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

First-Principles Studies on Sc 2 RuZ (Z = Si, Ge, Sn) Inverse Heusler Alloys: Structural, Electronic, and Transport Properties

The continuous demand for efficient, nontoxic, and thermally stable materials for room-temperature energy conversion motivates the exploration of novel thermoelectric systems beyond the traditional magnetic Heusler alloys. While full and half-Heusler compounds, especially Co-, Ni-, and Mn-based systems, have demonstrated promising thermoelectric properties, their typically high operating temperatures and magnetic complexities limit their applicability in ambient thermal management. In this context, we investigate whether Sc-based inverse Heusler alloys can offer a viable nonmagnetic alternative with competitive thermoelectric performance. In this work, we perform a systematic first-principles study of the inverse Heusler compounds Sc 2 RuZ (Z = Si, Ge, Sn), focusing on their structural, electronic, mechanical, and thermodynamic-thermoelectric properties. Density Functional Theory (DFT) was employed to compute optimized lattice structures and band dispersion, while dynamical stability was assessed via phonon calculations. Thermoelectric transport coefficients, including Seebeck coefficient, electrical conductivity, and thermal conductivity, were estimated using the semiclassical Boltzmann transport theory within the constant relaxation time approximation. Our results show that all Sc 2 RuZ compounds are thermodynamically stable semiconductors with indirect band gaps of 0.12–0.16 eV and exhibit high elastic moduli, especially Sc 2 RuSn, which demonstrates superior stiffness and incompressibility. Importantly, all compounds display promising room-temperature thermoelectric characteristics, including high Seebeck coefficients and power factors. These findings reveal that Sc 2 RuZ alloys represent a rare class of stable, nonmagnetic inverse Heusler semiconductors with intrinsic thermoelectric potential at room temperature, unlike many existing Heusler systems optimized for spintronics or high-temperature operation. This work expands the known design space for Heusler-based thermoelectrics and offers a theoretical basis for experimental realization of efficient, low-temperature, nonmagnetic thermoelectric materials.

alloys↗
DOE OSTI2021

Demonstration of reduced neoclassical energy transport in Wendelstein 7-X

Research on magnetic confinement of high-temperature plasmas has the ultimate goal of harnessing nuclear fusion for the production of electricity. Although the tokamak is the leading toroidal magnetic-confinement concept, it is not without shortcomings and the fusion community has therefore also pursued alternative concepts such as the stellarator. Unlike axisymmetric tokamaks, stellarators possess a three-dimensional (3D) magnetic field geometry. The availability of this additional dimension opens up an extensive configuration space for computational optimization of both the field geometry itself and the current-carrying coils that produce it. Such an optimization was undertaken in designing Wendelstein 7-X (W7-X), a large helical-axis advanced stellarator (HELIAS), which began operation in 2015 at Greifswald, Germany. A major drawback of 3D magnetic field geometry, however, is that it introduces a strong temperature dependence into the stellarator’s non-turbulent ‘neoclassical’ energy transport. Indeed, such energy losses will become prohibitive in high-temperature reactor plasmas unless a strong reduction of the geometrical factor associated with this transport can be achieved; such a reduction was therefore a principal goal of the design of W7-X. In spite of the modest heating power currently available, W7-X has already been able to achieve high-temperature plasma conditions during its 2017 and 2018 experimental campaigns, producing record values of the fusion triple product for such stellarator plasmas. The triple product of plasma density, ion temperature and energy confinement time is used in fusion research as a figure of merit, as it must attain a certain threshold value before net-energy-producing operation of a reactor becomes possible. Here we demonstrate that such record values provide evidence for reduced neoclassical energy transport in W7-X, as the plasma profiles that produced these results could not have been obtained in stellarators lacking a comparably high level of neoclassical optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

BeyondPlanck: I. Global Bayesian analysis of the Planck Low Frequency Instrument data

We describe the BEYONDPLANCK project in terms of our motivation, methodology, and main products, and provide a guide to a set of companion papers that describe each result in more detail. Building directly on experience from ESA’s Planck mission, we implemented a complete end-to-end Bayesian analysis framework for the Planck Low Frequency Instrument (LFI) observations. The primary product is a full joint posterior distribution P(ω | d), where ω represents the set of all free instrumental (gain, correlated noise, bandpass, etc.), astrophysical (synchrotron, free-free, thermal dust emission, etc.), and cosmological (cosmic microwave background – CMB – map, power spectrum, etc.) parameters. Some notable advantages of this approach compared to a traditional pipeline procedure are seamless end-to-end propagation of uncertainties; accurate modeling of both astrophysical and instrumental effects in the most natural basis for each uncertain quantity; optimized computational costs with little or no need for intermediate human interaction between various analysis steps; and a complete overview of the entire analysis process within one single framework. As a practical demonstration of this framework, we focus in particular on low-ℓ CMB polarization reconstruction with Planck LFI. In this process, we identify several important new effects that have not been accounted for in previous pipelines, including gain over-smoothing and time-variable and non-1/f correlated noise in the 30 and 44 GHz channels. Modeling and mitigating both previously known and newly discovered systematic effects, we find that all results are consistent with the ΛCDM model, and we constrained the reionization optical depth to τ = 0.066 ± 0.013, with a low-resolution CMB-based χ 2 probability to exceed of 32%. This uncertainty is about 30% larger than the official pipelines, arising from taking a more complete instrumental model into account. The marginal CMB solar dipole amplitude is 3362.7 ± 1.4 μK, where the error bar was derived directly from the posterior distribution without the need of any ad hoc instrumental corrections. We are currently not aware of any significant unmodeled systematic effects remaining in the Planck LFI data, and, for the first time, the 44 GHz channel is fully exploited in the current analysis. We argue that this framework can play a central role in the analysis of many current and future high-sensitivity CMB experiments, including LiteBIRD, and it will serve as the computational foundation of the emerging community-wide COSMOGLOBE effort, which aims to combine state-of-the-art radio, microwave, and submillimeter data sets into one global astrophysical model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Neoclassical transport of tungsten ion bundles in total-f neoclassical gyrokinetic simulations of a whole-volume JET-like plasma

Here, the application of a bundling technique to model the diverse charge states of tungsten impurity species in total-f gyrokinetic simulations is demonstrated. The gyrokinetic bundling method strategically groups tungsten ions of similar charge, optimizing computational efficiency. The initial radial configuration of these bundles and their respective charges are derived from a coronal approximation and the quasi-neutrality of the plasma. A low-density JET H-mode like plasma is simulated using the neoclassical version of XGC across the entire plasma volume, spanning from the magnetic axis to the divertor. An accumulation of tungsten is observed at the pedestal top, as a result of low-Z tungsten ions moving inward from the scrape-off-layer into the core region and high-Z tungsten ions moving outward from the core into the pedestal. This organization of the fluxes cannot be captured by a single tungsten-ion simulation. Large up-down poloidal asymmetries of tungsten form in the pedestal and strongly influence the direction of neoclassical fluxes. The temperature screening effect and its correlation with asymmetries are analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Stochastic Model Predictive Control With Gaussian Wind Direction Preview for Wake Steering

This article addresses the problem of wake steering control for wind farms that explicitly consider the tradeoff between farm-level power generation and yaw duty cycle under variable and uncertain wind conditions. A novel stochastic model predictive control (MPC) algorithm is presented, which utilizes a stochastic model of the freestream wind field components in a receding horizon framework to compute optimal yaw set points that maximize the expected value of the farm power while constraining the yaw actuation. Different configurations of the algorithm are evaluated using a steady-state wind farm simulator. The proposed stochastic MPC algorithm can plan control actions over a future prediction horizon based on probabilistic estimates of the incoming wind magnitude and direction.

17 WIND ENERGY↗
DOE OSTI2025

Ipopt Interface to Re::Solve Linear Solver

The software provides Ipopt optimization package an interface to the Re::Solve linear solver library. Re::Solve features GPU-resident direct and iterative linear solvers that could be used to accelerate optimization computations.

Alam, Maksudul [Oak Ridge National Laboratory (ORN↗
DOE OSTI2024

Efficacy of commercial recombinant HVT vaccines against a North American clade 2.3.4.4b H5N1 highly pathogenic avian influenza virus in chickens

The outbreak of clade 2.3.4.4b H5 highly pathogenic avian influenza (HPAI) in North America that started in 2021 has increased interest in applying vaccination as a strategy to help control and prevent the disease in poultry. Two commercially available vaccines based on the recombinant herpes virus of turkeys (rHVT) vector were tested against a recent North American clade 2.3.4.4b H5 HPAI virus isolate: A/turkey/Indiana/22-003707-003/2022 H5N1 in specific pathogen free white leghorn (WL) chickens and commercial broiler chickens. One rHVT-H5 vaccine encodes a hemagglutinin (HA) gene designed by the computationally optimized broadly reactive antigen method (COBRA-HVT vaccine). The other encodes an HA gene of a clade 2.2 virus (2.2-HVT vaccine). There was 100% survival of both chicken types COBRA-HVT vaccinated groups and in the 2.2-HVT vaccinated groups there was 94.8% and 90% survival of the WL and broilers respectively. Compared to the 2.2-HVT vaccinated groups, WL in the COBRA-HVT vaccinated group shed significantly lower mean viral titers by the cloacal route and broilers shed significantly lower titers by the oropharyngeal route than broilers. Virus titers detected in oral and cloacal swabs were otherwise similar among both vaccine groups and chicken types. To assess antibody-based tests to identify birds that have been infected after vaccination (DIVA-VI), sera collected after the challenge were tested with enzyme-linked lectin assay-neuraminidase inhibition (ELLA-NI) for N1 neuraminidase antibody detection and by commercial ELISA for detection of antibodies to the NP protein. As early as 7 days post challenge (DPC) 100% of the chickens were positive by ELLA-NI. ELISA was less sensitive with a maximum of 75% positive at 10DPC in broilers vaccinated with 2.2-HVT. Both vaccines provided protection from challenge to both types of chickens and ELLA-NI was sensitive at identifying antibodies to the challenge virus therefore should be evaluated further for DIVA-VI.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Robust multiresonant nonlocal metasurfaces by rational design

Abstract Dielectric metasurfaces supporting optical resonances have become a promising platform for quantum and nonlinear optics. However, resonant metasurfaces remain limited in their capacity to independently control the behavior of many distinct resonances despite efforts in computational optimization and inverse design. In this work, we overcome longstanding limitations by introducing a generalized rational design paradigm based on symmetry. Specifically, we use symmetry-broken metasurfaces with periodic “quadromer” lattices comprised of four nanostructures per unit cell to enable extensive control of multiple optical resonances. The rationally designed metasurfaces are readily fabricable, and we experimentally demonstrate metasurfaces that support up to four high Q-factor resonances with deliberately chosen free-space polarizations, spectral separations, and mode profiles. Our design paradigm may unlock new applications for multiresonant metasurfaces in quantum and nonlinear optics, optical sensing, and augmented reality displays.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Mitigating Phase Unbalance for Distribution Systems with High Penetrations of Solar PV (Final Technical Report)

Distribution system operators have traditionally limited unbalance among phases by maintaining similar loadings on each phase. High penetrations of distributed solar PV continually change the net loading on each phase, resulting in time-varying phase unbalances that can damage three-phase devices such as three-phase motors, violate grid codes, and increase technical losses. This project has developed several control strategies for the reactive power outputs of solar PV inverters in order to mitigate power quality issues related to phase unbalance. These control strategies include a decentralized approach that is solely based on local measurements, distributed and grouped approaches that consider subsets of loads and PV generators, and a centralized approach that leverages measurements from a variety of locations in order to compute optimal reactive power setpoints for each inverter. Variants of the controllers handle challenges relevant to practical implementations, including noisy measurements, delayed communications, and reactive power limits. Moreover, the project developed theory that provides convergence guarantees for systems with multiple interacting controllers as well as “balanceability” certificates that ensure satisfaction of phase unbalance requirements with variable loading. The controllers were integrated with NRECA’s Open Modeling Framework (omf.coop) and evaluated using actual distribution system models obtained from several NRECA member utilities. Application of the controllers results in significant improvements to phase unbalance in these test cases with decreases from base case levels of over 3% to under 0.5%, which is within the 2% IEC phase unbalance standard.

14 SOLAR ENERGY↗
DOE OSTI2021

Ultra-compact Imaging Technology (CRADA Final Report)

As part of the Cyclotron Road program, Synopic (formerly FlatCam LLC) sought to investigate the benefits of its depth-sensitive imaging techniques. The basic principle involved using a thin mask specially designed to encode light reaching a sensor, in conjunction with optimized computational algorithms, in such a manner that allowed for both miniaturization of imaging systems and improved computation toward high resolution, three-dimensional imaging. By thoroughly exploring the resolution, mask/optical design, and algorithmic capabilities of our depth sensitive technology, we expand the potential applications for (but not limited to) medical, consumer and industrial purposes. The project aimed to develop new imaging systems by building on previous work and using 1) materials capable of modulating and 2) sensors capable of measuring visible and longer wavelengths. Preliminary research was conducted to design, fabricate and characterize imaging systems with the goal of improving resolution, enhancing single capture, three-dimensional imaging, and extending depth of field of captured images. The primary goal was to determine whether adapting the depth-sensitive imaging system is feasible, and early results are promising.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗