Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point 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 OSTI2026

Measurement-informed Dynamic Aggregation of Distribution Systems

This paper proposes a measurement-informed dynamic aggregation methodology in order to create equivalent representations of distribution systems that are compatible with large-scale transmission analysis. By optimizing an equivalent feeder parameters using time-series measurements of active power, reactive power, and voltage at the Point of Interconnection (POI), the approach yields simplified yet dynamically accurate equivalents. Implemented in PSCAD with models of photovoltaic–battery systems, three-phase motors, and static loads, the method employs hybrid differential evolution and bounded least-squares optimization laying the foundation for for real-time state estimation and optimized sensor placement in distribution networks.

Ahmed, Kazi Ishrak [University of Tennessee, Knoxv↗
DOE OSTI2025

Machine Learning-Driven Conservative-to-Primitive Conversion in Hybrid Piecewise Polytropic and Tabulated Equations of State

We present a novel machine learning (ML)-based method to accelerate conservative-to-primitive inversion, focusing on hybrid piecewise polytropic and tabulated equations of state. Traditional root-finding techniques are computationally expensive, particularly for large-scale relativistic hydrodynamics simulations. To address this, we employ feedforward neural networks (NNC2PS and NNC2PL), trained in PyTorch (2.0+) and optimized for GPU inference using NVIDIA TensorRT (8.4.1), achieving significant speedups with minimal accuracy loss. The NNC2PS model achieves 𝐿 1 and 𝐿 ∞ errors of 4.54 × 10 −7 and 3.44 × 10−6, respectively, while the NNC2PL model exhibits even lower error values. TensorRT optimization with mixed-precision deployment substantially accelerates performance compared to traditional root-finding methods. Specifically, the mixed-precision TensorRT engine for NNC2PS achieves inference speeds approximately 400 times faster than a traditional single-threaded CPU implementation for a dataset size of 1,000,000 points. Ideal parallelization across an entire compute node in the Delta supercomputer (dual AMD 64-core 2.45 GHz Milan processors and 8 NVIDIA A100 GPUs with 40 GB HBM2 RAM and NVLink) predicts a 25-fold speedup for TensorRT over an optimally parallelized numerical method when processing 8 million data points. Moreover, the ML method exhibits sub-linear scaling with increasing dataset sizes. We release the scientific software developed, enabling further validation and extension of our findings. By exploiting the underlying symmetries within the equation of state, these findings highlight the potential of ML, combined with GPU optimization and model quantization, to accelerate conservative-to-primitive inversion in relativistic hydrodynamics simulations.

conservative-to-primitive conversion↗
DOE OSTI2023

Optimal Strategies for Hybrid Battery‐Storage Systems Design

As stationary hybrid energy‐storage systems (HESS) for power systems applications have recently drawn interest due to their enhanced performance and decreasing cost, developing systematic approaches for HESS design while considering controls is gaining traction. Herein, a method is presented to optimally design hybrid battery storage by proposing a mathematical modeling framework, formulated as a mixed integer linear programming model. The optimization is capable of handling multiple subsystems of batteries, considering their economic and technological performance. Decisions involve sizing of the batteries, optimal temporal and strategic dispatch to end uses, and energy sources for charging each battery. The applicability of the model is tested on four case studies for three battery chemistries representing distinct objectives: high‐power, high‐energy, and second life. Compared to traditionally designed battery storage with a homogeneous battery, optimally designed hybrid systems can save 12%–26% of system costs, depending on the nature of the dispatch profile. Findings point to design preference toward the second life battery supplemented with some high‐power or high‐energy battery capacity, or both. With the utilized electricity price structure, customers can experience approximately 10%–35% reduction in their bills.

Koleva, Mariya↗
DOE OSTI2025

A Fast Dynamic Internal Predictive Power Scheduling Approach for Power Management in Microgrids: Preprint

This paper presents a Dynamic Internal Predictive Power Scheduling (DIPPS) approach for optimizing power management in microgrids, particularly focusing on external power exchanges among diverse prosumers. DIPPS utilizes a dynamic objective function with a time-varying binary parameter to control the timing of power transfers to the external grid, facilitated by efficient usage of energy storage for surplus renewable power. The microgrid power scheduling problem is modeled as a mixed-integer nonlinear programming (MINLP-PS) and subsequently transformed into a mixed-integer linear programming (MILPPS) optimization through McCormick's relaxation to reduce computational complexity. A predictive window window with 6 data points is solved at an average of 0.92s, a 97.6% improvement over the 38.27s required for the MINLP-PS formulation, implying the numerical feasibility of the DIPPS approach for real-time implementation. Finally, the approach is validated against a static objective using real-world load data across three case studies with different time-varying parameters, demonstrating the ability of DIPPS to optimize power exchanges and efficiently utilize distributed resources while shifting the external power transfers to specified time durations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Rapid generation of optimal generalized Monkhorst-Pack grids

In this report computational modeling of the properties of crystalline materials has become an increasingly important aspect of materials research, consuming hundreds of millions of CPU-hours at scientific computing centers around the world each year, if not more. A routine operation in such calculations is the evaluation of integrals over the Brillouin zone. We have previously demonstrated that performing such integrals using generalized Monkhorst-Pack k-point grids can roughly double the speed of these calculations relative to the widely-used traditional Monkhorst-Pack grids. However the generation of optimal generalized Monkhorst-Pack grids is not implemented in most software packages due to the computational cost and difficulty of identifying the best grids. To address this problem, we present new algorithms that allow rapid generation of optimal generalized Monkhorst-Pack grids on the fly. We demonstrate that the grids generated by these algorithms are on average significantly more efficient than those generated using existing algorithms across a range of grid densities. For grids that correspond to a real-space supercell with at least 50 Å between lattice points, which is sufficient to converge density functional theory calculations within 1 meV/atom for nearly all materials, our algorithm finds optimized grids in an average of 0.19 s on a single processing core. To facilitate the widespread adoption of this approach, we present new open-source tools including a library designed for integration with third-party software packages.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tuning Structural and Electronic Configuration of FeN 4 via External S for Enhanced Oxygen Reduction Reaction

The Fe–N–C material represents an attractive oxygen reduction reaction electrocatalyst, and the FeN 4 moiety has been identified as a very competitive catalytic active site. Fine tuning of the coordination structure of FeN 4 has an essential impact on the catalytic performance. Herein, we construct a sulfur‐modified Fe–N–C catalyst with controllable local coordination environment, where the Fe is coordinated with four in‐plane N and an axial external S. The external S atom affects not only the electron distribution but also the spin state of Fe in the FeN 4 active site. The appearance of higher valence states and spin states for Fe demonstrates the increase in unpaired electrons. With the above characteristics, the adsorption and desorption of the reactants at FeN 4 active sites are optimized, thus promoting the oxygen reduction reaction activity. This work explores the key point in electronic configuration and coordination environment tuning of FeN 4 through S doping and provides new insight into the construction of M–N–C‐based oxygen reduction reaction catalysts.

25 ENERGY STORAGE↗
DOE OSTI2022

On the closedness and geometry of tensor network state sets

Tensor network states (TNS) are a powerful approach for the study of strongly correlated quantum matter. The curse of dimensionality is addressed by parametrizing the many-body state in terms of a network of partially contracted tensors. These tensors form a substantially reduced set of effective degrees of freedom. In practical algorithms, functionals like energy expectation values or overlaps are optimized over certain sets of TNS. Concerning algorithmic stability, it is important whether the considered sets are closed because, otherwise, the algorithms may approach a boundary point that is outside the TNS set and tensor elements diverge. Here we discuss the closedness and geometries of TNS sets, and we propose regularizations for optimization problems on non-closed TNS sets. We show that sets of matrix product states (MPS) with open boundary conditions, tree tensor network states, and the multiscale entanglement renormalization ansatz are always closed, whereas sets of translation-invariant MPS with periodic boundary conditions (PBC), heterogeneous MPS with PBC, and projected entangled pair states are generally not closed. The latter is done using explicit examples like the W state, states that we call two-domain states, and fine-grained versions thereof.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Vapor liquid equilibrium of pure and aqueous Methyl diethanolamine (MDEA), 2-Dimethylmonoethanolamine (DMMEA), N-Methylmonoethanolamine (MMEA) and 1-(2-Aminoethyl)piperazine (AEP): Experimental results and modeling

This work presents new vapor-liquid data for the pure components 2-dimethylmonoethanolamine (DMMEA) and N-methylethanolamine (MMEA) and for the binary systems of aqueous methyldiethanolamine (MDEA), DMMEA, MMEA, and 1-(2-aminoethyl)piperazine (AEP). New Antoine models are developed for the pure component systems, and new and improved NRTL models for the binary systems. In the model development, together with new data from this work, all available data on VLE, excess heat of mixing, H E , freezing point depression, and specific heat were gathered, evaluated and used selectively for parameter optimization. Emphasis has been put on a best possible representation of amine volatility to enable reasonable predictions of amine emissions from absorption plants and to form a basis for subsequent models of ternary systems with CO 2 . In conclusion, model predictions are compared with experimental results and other models, source by source, and the overall results are satisfactory.

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

Feasibility of DEIM for retrieving the initial field via dimensionality reduction

When parameter estimation is solved in a high-dimensional space, the dimensionality reduction strategy becomes the primary consideration for alleviating the tremendous computational cost. Here, the discrete empirical interpolation method (DEIM) is explored to retrieve the initial condition (IC) by combining the polynomial chaos (PC) based ensemble Kalman filter (i.e. PC-EnKF), where a non-intrusive PC expansion is considered as a surrogate model in place of the forward model in the prediction step of the ensemble Kalman filter, resulting in fewer forward model integrations but with a comparable accuracy as Monte Carlo-based approaches. The DEIM acts as a hyper-reduction tool to provide the low-dimensional input for the high-dimensional initial field, which can be reconstructed using the information on the sparse interpolation grid points that is adaptively obtained through PC-EnKF data assimilation method. Thus an innovative framework to reconstruct the IC is developed. The detailed procedure at each assimilation iteration includes: the determination of the spatial interpolation points, the estimation of the initial values on the interpolation locations using the optimal observations, and the reconstruction of IC in the full space. The current study uses the reconstruction field of initial conditions of the Navier-Stokes equations as an example to illustrate the efficacy of our method. The experimental results demonstrate the proposed algorithm achieves a satisfactory reconstruction for the initial field. The proposed method helps to extend the applicable area of DEIM in solving inverse problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Nanoscale Control of Intrinsic Magnetic Topological Insulator MnBi 2 Te 4 Using Molecular Beam Epitaxy: Implications for Defect Control

Intrinsic magnetic topological insulators have emerged as a promising platform to study the interplay between the topological surface states and ferromagnetism. This unique interplay can give rise to a variety of exotic quantum phenomena, including the quantum anomalous Hall effect and axion insulating states. Here, in this study, utilizing molecular beam epitaxy (MBE), we present a comprehensive study of the growth of MnBi 2 Te 4 thin films on Si (111), epitaxial graphene, and highly ordered pyrolytic graphite substrates. By combining a suite of in situ characterization techniques, we obtain critical insights into the nanoscale control of MnBi 2 Te 4 epitaxial growth. First, we extract the free energy landscape for the epitaxial relationship as a function of the in-plane angular distribution. Then, by employing an optimized layer-by-layer growth, we determine the chemical potential and Dirac point of the thin film at different thicknesses and how this quantity is manifested by the dopant compensation from different antisite defects. Overall, these results establish a foundation for understanding the growth kinetics of MnBi 2 Te 4 and pave the way for future applications of MBE-grown thin films in emerging topological quantum materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Photometric metallicity for 694 233 Galactic giant stars from Gaia DR3 synthetic Strömgren photometry: Metallicity distribution functions of halo substructures

We used two previous calibrations and the standardised synthetic photometry in the v, b, and y Strömgren passbands from Gaia DR3 BP/RP spectra to obtain photometric metallicities for a selected sample of 694 233 old Galactic giant stars having |b|> 20.0° and parallax uncertainties lower than 10%. The zero point of both sets of photometric metallicities has been shifted to ensure an optimal match with the spectroscopic [Fe/H] values for 44 785 stars in common with APOGEE DR17, focusing on the metallicity range where they provide the highest accuracy. The metallicities derived in this way from one calibration display a typical accuracy of ≲0.1 dex and 1σ precision ≲0.2 dex in the range -2.2 ≲ [Fe/H] ≲ - 0.4, while they show a systematic trend with [Fe/H] at higher metallicity, beyond the applicability range of the relation. Those derived from the other calibration display, in general, reduced precision, and lower accuracy in the metal-poor regime, but have a median accuracy < 0.05 dex for [Fe/H] ≥ - 0.8. These results are confirmed and, consequently, the metallicities are validated, by comparison with large sets of spectroscopic metallicities from various surveys. The newly obtained metallicities are used to derive metallicity distributions for several previously identified substructures in the Galactic halo with an unprecedented number of stars. The catalogue including the two sets of metallicities and the associated uncertainties is made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Optimal celestial bodies for dark matter detection

A wide variety of celestial bodies have been considered as dark matter detectors. Which stands the best chance of delivering the discovery of dark matter? Which is the most powerful dark matter detector? We investigate a range of objects, including the Sun, Earth, Jupiter, Brown Dwarfs, White Dwarfs, Neutron Stars, Stellar populations, and Exoplanets. We quantify how different objects are optimal dark matter detectors in different regimes by deconstructing some of the in-built assumptions in these search sensitivities, including observation potential and particle model assumptions. We find new constraints and future sensitivities across a range of dark matter annihilation final states. We quantify mediator properties leading to detectable celestial-body energy injection or Standard Model fluxes, and show how different objects can be expected to deliver corroborating signals. We discuss different search strategies, their opportunities and limitations, and the interplay of regimes where different celestial objects are optimal dark matter detectors. Deconstructing the assumptions of these searches leads us to point out a new search using the Galactic center stellar population that can provide greater sensitivity to the dark matter-nucleon scattering cross section than the Sun, despite being significantly further away in our Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

BEAM COLLIMATION IN THE PIP-II LINAC TO BOOSTER TRANSFER LINE

The new PIP-II superconducting linac will deliver a 2 mA average H⁻ beam to the existing Booster synchrotron. The injected beam is accumulated by charge exchange over approximately 300 turns; phase space painting is used to mitigate space charge effects. To limit the power load on the internal waste beam absorber from the transverse tails of the H⁻ distribution missing the foil, the beam will be collimated in both planes in the linac to Booster transfer line using compact collimators of a novel design. Both the number of parasitic hits and the fraction of the beam missing the foil are sensitive functions of the H⁻ beam centroid position with respect to the edge of the foil. The positioning of the collimation is constrained by the availability of suitable space in the transfer line lattice, by specifics of the collimator design, by the phase space orientation at the collimator, and by the betatron phase advance to the foil needed to achieve proper orientation of the spatial distribution at the injection point. In this contribution, we describe the procedure by which collimator positions were optimized. We then discuss the expected performance of the overall system.

Johnson, David E.↗
DOE OSTI2020

Evaluation of Thermolytic Hydrogen Generation Rate Models at High-Temperature/High-Hydroxide Regimes

This report describes the results of testing performed to extend the applicable ranges of temperature and hydroxide concentration for use within the Glycolate and Global Total Organic Carbon (TOC) Hydrogen Generation Rate (HGR) expressions. Seven experimental conditions (six simulants of the 242-25H Evaporator system chosen as a D-optimal set of experiments and a single test conducted at an elevated boiling point of 170 °C) were investigated in the presence of sodium glycolate and Xiameter TM AFE-1010. Glycolate was employed to study the extension of the Glycolate Thermolytic HGR expression while Xiameter TM AFE-1010 was employed to study the extension of the Global TOC Thermolytic HGR expression. The following conclusions were derived from this testing: The Glycolate Thermolytic HGR expression may be confidently used to predict thermolytic HGRs from glycolate at temperatures as high as 170 °C and hydroxide concentrations as high as 23 M.; The hydroxide and temperature-dependence predicted by the Global TOC Thermolytic HGR expression has been confirmed at temperatures as high as 170 °C and hydroxide concentrations as high as 23 M, suggesting that the Global TOC Thermolytic HGR expression may be used at these ranges.; Methane was observed from tests with Xiameter TM AFE-1010 at production rates higher than those observed for hydrogen. These rates were observed at temperatures higher than 100 °C.; Preliminary models suggest that increasing hydroxide/temperature causes an increase in Methane Generation Rate (MGR) from Xiameter TM AFE-1010. The following recommendations are based on this testing: The existing equations for thermolytic HGR from glycolate and non-glycolate organics should be used at Concentration, Storage, and Transfer Facilities (CSTF) storage and evaporation conditions, including temperatures and hydroxide concentrations exhibited in the 242-25H Evaporator.; Further investigation should be made into the influence of methylsilanes on CSTF flammability. This investigation should include: determination of the types of methylsilanes historically added to the CSTF, determination of methane formation rates from each type of methylsilane, and determination of the extent of degradation of methylsilanes in CSTF waste.; Characterization techniques should be developed by Savannah River National Laboratory (SRNL) to assist in the speciation of methylsilane-containing waste in the CSTF.; Additional testing with radioactive waste should be performed to determine the MGRs possible in radioactive waste and better inform model predictions made from testing with simulants.

08 HYDROGEN↗
DOE OSTI2024

Automated AI-driven Molecular Design for Therapeutic Discovery

In recent years, artificial intelligence and machine learning (AI/ML) approaches have revolutionized the process of designing new therapeutics, enabling scientists to rapidly respond to emerging threats from various pathogens. A prime example is the SARS-CoV-2 main protease, a key target for the development of antiviral inhibitors. In this study, we employed a novel, integrated approach that combines AI-driven iterative design of inhibitor candidates, screening based on physio-chemical properties and toxicity, physics-based computational modeling of protein-inhibitor interactions, and AI-assisted analysis of Native MS biophysical assay and characterization of designed candidates. Our deep learning 3D-scaffold model, which uses an input scaffold as a starting point, generated tens of thousands of compounds while preserving the key scaffold. To optimize these candidates, we calculated a comprehensive set of 136 descriptors, including both 2D and 3D molecular features, for compounds targeting the SARS-CoV-2 Main protease (Mpro) and a neurodegenerative disease-associated protein, cyclophilin (Cyp). The generated compounds were initially filtered based on their properties and then ranked according to their predicted binding affinity using our automated modeling and ML methods. Experimental validation of the Mpro candidates showing inhibitory activity demonstrates that our workflow can expedite the therapeutic discovery.

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

Exact signed distance fields using parallel Fast Sweeping Method

Signed distance fields are often used in multiphysics simulations to track material interfaces. We present a simple methodology based on the fast sweeping method to generate the exact signed distance from triangular meshes and linear paths on Cartesian grids. The methodology propagates the closest primitive to the boundary to the rest of the domain following the characteristics. A local upwind criterion is used to decide between the new and existing closest primitive at each grid point while capturing the correct sign of the global function. The methodology has optimal computational complexity and runs efficiently in distributed-memory architectures. We include 2D and 3D test cases along with a resolution study up to 0.512 trillion zones and 1,000 computer cores. The solution strategy can also be applied to other types of meshes or collections of primitives.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Lossy checkpoint compression in full waveform inversion: a case study with ZFPv0.5.5 and the overthrust model

This paper proposes a new method that combines checkpointing methods with error-controlled lossy compression for large-scale high-performance full-waveform inversion (FWI), an inverse problem commonly used in geophysical exploration. This combination can significantly reduce data movement, allowing a reduction in run time as well as peak memory. In the exascale computing era, frequent data transfer (e.g., memory bandwidth, PCIe bandwidth for GPUs, or network) is the performance bottleneck rather than the peak FLOPS of the processing unit. Like many other adjoint-based optimization problems, FWI is costly in terms of the number of floating-point operations, large memory footprint during backpropagation, and data transfer overheads. Past work for adjoint methods has developed checkpointing methods that reduce the peak memory requirements during backpropagation at the cost of additional floating-point computations. Combining this traditional checkpointing with error-controlled lossy compression, we explore the three-way tradeoff between memory, precision, and time to solution. We investigate how approximation errors introduced by lossy compression of the forward solution impact the objective function gradient and final inverted solution. Empirical results from these numerical experiments indicate that high lossy-compression rates (compression factors ranging up to 100) have a relatively minor impact on convergence rates and the quality of the final solution.

58 GEOSCIENCES↗
DOE OSTI2021

Enabling Long-range Exploration in Minimization of Multimodal Functions

We consider the problem of minimizing multi-modal loss functions with a large number of local optima. Since the local gradient points to the direction of the steepest slope in an infinitesimal neighborhood, an optimizer guided by the local gradient is often trapped in a local minimum. To address this issue, we develop a novel nonlocal gradient to skip small local minima by capturing major structures of the loss’s landscape in black-box optimization. The nonlocal gradient is defined by a directional Gaussian smoothing (DGS) approach. The key idea of DGS is to conducts 1D long-range exploration with a large smoothing radius along d orthogonal directions in Rd, each of which defines a nonlocal directional derivative as a 1D integral. Such long-range exploration enables the nonlocal gradient to skip small local minima. The d directional derivatives are then assembled to form the nonlocal gradient. We use the Gauss-Hermite quadrature rule to approximate the d 1D integrals to obtain an accurate estimator. The superior performance of our method is demonstrated in three sets of examples, including benchmark functions for global optimization, and two real-world scientific problems.

Zhang, Jiaxin↗