Search NASA⌕ Search

SEARCH · Search NASA

Results for “derivative free 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Data-driven Optimal Control Strategy for Virtual Synchronous Generator via Deep Reinforcement Learning Approach

This paper aims at developing a data-driven optimal control strategy for virtual synchronous generator (VSG) in the scenario where no expert knowledge or requirement for system model is available. Firstly, the optimal and adaptive control problem for VSG is transformed into a reinforcement learning task. Specifically, the control variables, i.e., virtual inertia and damping factor, are defined as the actions. Meanwhile, the active power output, angular frequency and its derivative are considered as the observations. Moreover, the reward mechanism is designed based on three preset characteristic functions to quantify the control targets: (1) maintaining the deviation of angular frequency within special limits; (2) preserving well-damped oscillations for both the angular frequency and active power output; (3) obtaining slow frequency drop in the transient process. Next, to maximize the cumulative rewards, a decentralized deep policy gradient algorithm, which features model-free and faster convergence, is developed and employed to find the optimal control policy. With this effort, a data-driven adaptive VSG controller can be obtained. By using the proposed controller, the inverter-based distributed generator can adaptively adjust its control variables based on current observations to fulfill the expected targets in model-free fashion. Finally, simulation results validate the feasibility and effectiveness of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Polybenzimidazole-derived carbon molecular sieve hollow fiber membranes with tailored oxygen selective transport

Polymer membrane carbonization is a promising strategy to obtain carbon molecular sieve (CMS) membrane materials with tailorable cavity size distributions for high gas selectivities beyond that attained from traditional polymer membrane materials. Despite their already demonstrated exceptional separation performance characteristics for several gas separation applications (e.g. CO 2 /CH 4 and H 2 /CO 2 ), further development of CMS membranes having high O 2 /N 2 perm-selectivity for energy-efficient high-purity O 2 production has proven challenging. Herein, we explore CMS hollow fiber membranes (HFMs) derived from highly rigid and tightly-packed polybenzimidazole (PBI) materials. Nearly defect-free PBI-derived CMS HFMs were fabricated and evaluated for O 2 /N 2 separation performance for the first time. The micro-structural and O 2 selective gas separation characteristics of PBI-CMS HFMs pyrolyzed at different pyrolysis conditions were intensively studied. With rigorous carbonization protocol optimization, we obtained highly O 2 permselective PBI-CMS HFMs having O 2 /N 2 selectivities approaching 14, enabling high purity O 2 production. The PBI-CMS HFM O 2 /N 2 separation performance dependence on operating conditions (e.g. temperature and feed pressure) are reported.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A General Purpose Real-Time Optimization Strategy Applied to Minimizing Simultaneous Heating and Cooling in Buildings

Supervisory control strategies in heating, ventilating, and air-conditioning (HVAC) systems hold the potential for significant energy savings but they are difficult to set up and often depend on unreliable sensor measurements. In this paper, we propose a sensor-free approach to minimizing the energy penalties due to simultaneous heating and cooling in variable-air-volume (VAV) systems with reheat. A general-purpose real-time optimizer is described that is easy to set-up and integrate with existing control logic. A cost function is derived that quantifies the simultaneous heating and cooling penalty that does not require any sensor measurements. The optimizer then minimizes this cost function by automatically adjusting the supply temperature setpoint in the air-handling unit. The paper describes the approach and shows results from application from field trials on a 27-zone building.

Salsbury, Timothy↗
DOE OSTI2020

Gel‐Derived Amorphous Bismuth–Nickel Alloy Promotes Electrocatalytic Nitrogen Fixation via Optimizing Nitrogen Adsorption and Activation

Abstract To achieve the electrochemical nitrogen reduction reaction (NRR) for efficient and sustainable NH 3 production, catalysts should exhibit high selectivity and activity with optimal adsorption energy. Herein we developed a three‐dimensional (3D) amorphous BiNi alloy toward a significantly enhanced NRR compared with its crystalline and metal counterparts. Ni alloying enables the chemisorption of nitrogen and the lower free‐energy change for the *NNH formation, and the 3D alloy electrocatalyst exhibits high catalytic activity for NH 3 production with a yield rate of 17.5 μg h −1 mg cat −1 and Faradaic efficiency of 13.8 %. The enhanced electron transfer and increased electrochemical surface area were revealed in the interconnected porous scaffold, affording it sufficiently efficient and stable activity for potential practical applications. This work offers new insights into optimizing the adsorption energy of reactants and intermediates combined with tuning the crystallinity of NRR electrocatalysts.

Fang, Zhiwei↗
DOE OSTI2022

Solvent‐free enzymatic esterification of free fatty acids with glycerol for biodiesel application: Optimized using the Taguchi experimental method

Abstract Presence of free fatty acids along with glycerides poses a technical difficulty for biodiesel production. This work used a Taguchi L9 design to optimize the solvent‐free enzymatic process to result in the esterification of oleic acid with glycerol. Under optimal conditions the esterification reaction temperature of 60°C, enzyme dose of 5 wt%, glycerol: oleic acid molar ratio of 5:1, and reaction time of 3 h, a 75.235 ± 2.19% conversion of oleic acid to esters was achieved. With the addition of molecular sieves, the conversion increased to 86.73% ± 1.09%. However, using the parameters predicted by Taguchi design (60°C, 5 wt%, 5:1, and 4.5 h), 88.5% ± 1.11% of oleic acid could be converted to esters derivative. Diglycerides were the major product, and the reaction equilibrium was attained after 4 h. The immobilized enzyme could be used up to seven times with only a 10% reduction in the conversion. Thus, the process can efficiently reduce the free fatty acid content of oil to make it suitable for biodiesel production.

09 BIOMASS FUELS↗
DOE OSTI2022

Data for Solvent-Free Enzymatic Esterification of Free Fatty Acids with Glycerol for Biodiesel Application: Optimized Using the Taguchi Experimental Method

Presence of free fatty acids along with glycerides poses a technical difficulty for biodiesel production. This work used a Taguchi L9 design to optimize the solvent-free enzymatic process to result in the esterification of oleic acid with glycerol. Under optimal conditions the esterification reaction temperature of 60°C, enzyme dose of 5 wt%, glycerol: oleic acid molar ratio of 5:1, and reaction time of 3 h, a 75.235 ± 2.19% conversion of oleic acid to esters was achieved. With the addition of molecular sieves, the conversion increased to 86.73% ± 1.09%. However, using the parameters predicted by Taguchi design (60°C, 5 wt%, 5:1, and 4.5 h), 88.5% ± 1.11% of oleic acid could be converted to esters derivative. Diglycerides were the major product, and the reaction equilibrium was attained after 4 h. The immobilized enzyme could be used up to seven times with only a 10% reduction in the conversion. Thus, the process can efficiently reduce the free fatty acid content of oil to make it suitable for biodiesel production.

Conversion↗
DOE OSTI2026

SENTRA: A Modular Computational Graph Framework for Critical Mineral and Materials Supply Chains: Part I: Network Construction Latent-Quantity Estimation, and Temporal Graph Forecasting

Global supply chains for critical minerals and materials are complex, evolving networks of countries, products, production stages, and trade relationships. Existing analytical approaches are limited by fragmented data and static network representations that do not capture the dynamic production dependencies linking raw materials, intermediate products, and final goods across multiple countries. Trade and production statistics provide only a partial view of domestic production, inventories, and material flows, making it difficult to identify indirect sourcing pathways, hidden dependencies, and embedded foreign exposures. This paper introduces the Supply Chain Exposure Network Tracking and Risk Assessment (SENTRA) framework, a modular graph-based computational framework for constructing, analyzing, and forecasting dynamic supply chain networks. As the first paper in a three-part methodological series, it establishes the computational foundation of SENTRA by constructing a temporal attributed multi-relational graph whose nodes represent product–country pairs and whose edges encode observed trade and within-country value-chain relationships. Statistical estimation and constrained optimization recover latent production, final demand, and product input dependency coefficients while enforcing economic accounting constraints. Graph-derived exposure measures quantify direct, transshipment, value-chain, and multi-hop supply chain dependencies independently of the forecasting model. A temporal graph forecasting architecture based on a relational graph neural network then forecasts the evolution of the graph under mass-balance constraints with distribution-free conformal uncertainty quantification. Validation on the global aluminum supply chain shows that the learned graph representations recover economically meaningful supply chain structure, accurately forecast out-of-sample trade relationships, and produce well-calibrated prediction intervals. Subsequent papers apply this computational foundation to exposure assessment, disruption analysis, and scenario-based policy analysis, and extend the framework to multimaterial supply chain modeling and decision support.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Universal Compiling and (No-)Free-Lunch Theorems for Continuous-Variable Quantum Learning

Quantum compiling, where a parameterized quantum circuit is trained to learn a target unitary, is an important primitive for quantum computing that can be used as a subroutine to obtain optimal circuits or as a tomographic tool to study the dynamics of an experimental system. While much attention has been paid to quantum compiling on discrete-variable hardware, less has been paid to compiling in the continuous-variable paradigm. Here we motivate several, closely related, short-depth continuous-variable algorithms for quantum compilation. We analyze the trainability of our proposed cost functions and numerically demonstrate our algorithms by learning arbitrary Gaussian operations and Kerr nonlinearities. We further make connections between this framework and quantum learning theory in the continuous-variable setting by deriving no-free-lunch theorems. These generalization bounds demonstrate a linear resource reduction for learning Gaussian unitaries using entangled coherent-Fock states and an exponential resource reduction for learning arbitrary unitaries using two-mode-squeezed states.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Metal–organic frameworks as O 2 -selective adsorbents for air separations

Oxygen is a critical gas in numerous industries and is produced globally on a gigatonne scale, primarily through energy-intensive cryogenic distillation of air. The realization of large-scale adsorption-based air separations could enable a significant reduction in associated worldwide energy consumption and would constitute an important component of broader efforts to combat climate change. Certain small-scale air separations are carried out using N 2 -selective adsorbents, although the low capacities, poor selectivities, and high regeneration energies associated with these materials limit the extent of their usage. In contrast, the realization of O 2 -selective adsorbents may facilitate more widespread adoption of adsorptive air separations, which could enable the decentralization of O 2 production and utilization and advance new uses for O 2 . Here, we present a detailed evaluation of the potential of metal–organic frameworks (MOFs) to serve as O 2 -selective adsorbents for air separations. Drawing insights from biological and molecular systems that selectively bind O 2 , we survey the field of O 2 -selective MOFs, highlighting progress and identifying promising areas for future exploration. As a guide for further research, the importance of moving beyond the traditional evaluation of O 2 adsorption enthalpy, ΔH, is emphasized, and the free energy of O 2 adsorption, ΔG, is discussed as the key metric for understanding and predicting MOF performance under practical conditions. Based on a proof-of-concept assessment of O 2 binding carried out for eight different MOFs using experimentally derived capacities and thermodynamic parameters, we identify two existing materials and one proposed framework with nearly optimal ΔG values for operation under user-defined conditions. While enhancements are still needed in other material properties, the insights from the assessments herein serve as a guide for future materials design and evaluation. Computational approaches based on density functional theory with periodic boundary conditions are also discussed as complementary to experimental efforts, and new predictions enable identification of additional promising MOF systems for investigation.

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

Physics-preserving enriched Galerkin method for a fully-coupled thermo-poroelasticity model

This paper proposes a new numerical method for a fully-coupled, quasi-static thermo-poroelasticity model in a unified enriched Galerkin (EG) method framework. In our method, the mechanics sub-problem is solved using a locking-free EG method, and the flow and heat sub-problems are solved using a locally-conservative EG method. The proposed method offers mass and energy conservation properties with much lower costs than other methods with the same properties, including discontinuous Galerkin methods and mixed finite element methods. The well-posedness and optimal a priori error estimates are carefully derived. Here, several numerical tests confirm the theoretical optimal convergence rates and the mass and energy conservation properties of the new method.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Shear band velocity and activation volume during shear deformation by acoustic emission in a Zr-based bulk metallic glass

Here, an analysis is presented for the shear band velocity of plastic deformation and activation volume of the shear transformation zone (STZ) in a Zr 60 Cu 30 Al 10 bulk metallic glass (BMG). Compared to the values calculated by the stick-slip model, the values of shear band velocity calculated after optimization by the acoustic emission technique are much larger than those reported previously. Plastic flow in BMGs is found to be described by a power-law correlation between the maximum velocity of the shear band and the strain rates. Combining the free volume theory and the STZ theory, the relationship between the volume change in the STZ and the shear band velocity is derived as $Ω=A$(27–ln$v/l$). The value of A, which is a material dependent constant, is 5.43 × 10 –3 for the current BMG at room temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A scalable matrix-free spectral element approach for unsteady PDE constrained optimization using PETSc/TAO

In this work, we provide a new approach for the efficient matrix-free application of the transpose of the Jacobian for the spectral element method for the adjoint-based solution of partial differential equation (PDE) constrained optimization. This results in optimizations of nonlinear PDEs using explicit integrators where the integration of the adjoint problem is not more expensive than the forward simulation. Solving PDE constrained optimization problems entails combining expertise from multiple areas, including simulation, computation of derivatives, and optimization. The Portable, Extensible Toolkit for Scientific computation (PETSc) together with its companion package, the Toolkit for Advanced Optimization (TAO), is an integrated numerical software library that contains an algorithmic/software stack for solving linear systems, nonlinear systems, ordinary differential equations, differential algebraic equations, and large-scale optimization problems and, as such, is an ideal tool for performing PDE-constrained optimization. This paper describes an efficient approach in which the software stack provided by PETSc/TAO can be used for large-scale nonlinear time-dependent problems. Time integration can involve a range of high-order methods, both implicit and explicit. The PDE-constrained optimization algorithm used is gradient-based and seamlessly integrated with the simulation of the physical problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Differentiable lagrangian shock hydrodynamics with application to stable shock acceleration of density interfaces

We develop a gradient based optimization approach for the equations of compressible, Lagrangian hydrodynamics and demonstrate how it can be employed to automatically uncover strategies to control hydrodynamic instabilities arising from shock acceleration of density interfaces. Strategies for controlling the Richtmyer-Meshkov instability (RMI) are of great benefit for inertial confinement fusion (ICF) where shock interactions with many small imperfections in the density interface lead to instabilities which rapidly grow over time. These instabilities lead to mixing which, in the case of laser driven ICF, quenches the runaway fusion process ruining the potential for positive energy return. Here, we demonstrate that control of these instabilities can be achieved by optimization of initial conditions with ( > 100) parameters. Optimizing over a large parameter space like this is not possible with gradient-free optimization strategies. This requires computation of the gradient of the outputs of a numerical solution to the equations of Lagrangian hydrodynamics with respect to the inputs. We show that the efficient computation of these gradients is made possible via a judicious application of (i) adjoint methods, the exact formal representation of sensitivities involving partial differential equations, and (ii) automatic differentiation (AD), the algorithmic calculation of derivatives of functions. Careful regularization of multiple operators including artificial viscosity and timestep control is required. We perform design optimization of > 100 parameter energy field driving the Richtmyer Meshkov instability showing significant suppression while simultaneously enhancing the acceleration of the interface relative to a nominal baseline case.

Hydrophysics↗
DOE OSTI2020

Pore-scale view of microbial turnover: combining 14 C imaging, μCT, and zymography after adding soluble carbon to soil pores of specific sizes

The location of microorganisms and substrates within soil pore networks plays a crucial role for organic carbon (C) processing, its microbial utilization and turnover, and has direct consequences for C and nutrient cycling. An optimal approach to quantify responses to new C inputs from microorganisms residing in specific pores is the addition of new C to pores of target sizes in undisturbed soil cores. We used the matric potential approach to add 14 C labeled glucose to small (< 40 μm, root free) or large (60-180 μm, potentially inhabited by roots) pores of undisturbed soil cores. Localization of glucose-derived C via 14 C imaging was related to pore size distributions and connectivity, characterized via X-ray computed micro-tomography (μCT), and to β-glucosidase activity, characterized via zymography. After two-week incubations, 1.3 times more glucose was mineralized ( 14 CO 2 ) when it was added to the large pores, however more 14 C remained in microbial biomass when glucose was added to the small pores. Consequently, while utilizing the same amounts of easily available C, the microorganisms localized in the large pores had faster turnover compared to microorganisms in small pores. Stronger associations between β-glucosidase activity and glucose-derived C were observed when glucose was added to the large pores. We conclude that 1) the matric potential approach allows placing, albeit not exactly, soluble substrates into pores of target diameter range, and 2) microorganisms localized in large pores respond to new C inputs with faster turnover, greater growth, and more intensive enzyme production compared to those inhabiting the small pores.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

FLOWERS AEP: An Analytical Model for Wind Farm Layout Optimization

Annual energy production (AEP) is commonly used in objective functions for wind farm layout optimization. AEP is proportional to wind farm power production integrated over an annual distribution of free-stream wind conditions. Physics-based estimates of wind farm power production typically rely on low-fidelity engineering wake models that approximate the steady-state wind farm flow field. AEP estimates are then obtained by performing independent simulations for discrete wind conditions and using rectangular quadrature to account for each condition's expected frequency of occurrence. Depending on the number of simulated discrete wind conditions, this numerical integral could be hampered by poor accuracy or high computational costs. The FLOWERS AEP model instead poses an analytical integral of the engineering wake model over the variable wind conditions, yielding a closed-form, analytical function for wind farm AEP. This paper derives the analytical functions for FLOWERS AEP and its derivatives with respect to turbine position, which are useful for gradient-based wind farm layout optimization, in nondimensional form. We then analyze the benefits of the FLOWERS AEP model over conventional reference models, focusing on its low cost, adequate wake loss predictions, and smooth design space. Although the FLOWERS approach is found to predict the exact value of AEP with some error relative to the reference model (within 14% on average), it dramatically reduces computation time by an order of magnitude, produces a qualitatively similar design space at relatively low resolution, and yields comparable optimal layouts. This significant speed improvement is critical in layout optimization applications, where determining an optimal layout in an efficient manner is more important than precise AEP prediction.

17 WIND ENERGY↗
DOE OSTI2024

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Chemical Vapor Deposition for Atomically Dispersed and Nitrogen Coordinated Single Metal Site Catalysts

Abstract Atomically dispersed and nitrogen coordinated single metal sites (M‐N‐C, M=Fe, Co, Ni, Mn) are the popular platinum group‐metal (PGM)‐free catalysts for many electrochemical reactions. Traditional wet‐chemistry catalyst synthesis often requires complex procedures with unsatisfied reproducibility and scalability. Here, we report a facile chemical vapor deposition (CVD) strategy to synthesize the promising M‐N‐C catalysts. The deposition of gaseous 2‐methylimidazole onto M‐doped ZnO substrates, followed by an in situ thermal activation, effectively generated single metal sites well dispersed into porous carbon. In particular, an optimal CVD‐derived Fe‐N‐C catalyst exclusively contains atomically dispersed FeN 4 sites with increased Fe loading relative to other catalysts from wet‐chemistry synthesis. The catalyst exhibited outstanding oxygen‐reduction activity in acidic electrolytes, which was further studied in proton‐exchange membrane fuel cells with encouraging performance.

Liu, Shengwen↗
DOE OSTI2020

Automated Fitting of a Semi-empirical Multiphase Equation of State for Carbon

The equation of state (EOS) of carbon is important in high explosive, geophysical, and inertial confinement fusion applications. Within the semi-empirical Sesame framework, the EOS of each phase is represented by a sum of cold, vibrational, and thermal electronic Helmholtz free energy contributions. Each phase has ~5-10 independent parameters that are adjusted to reproduce single-phase data (e.g., thermal expansion, isothermal compression) as well as experimentally- and computationally- derived phase boundaries. Manual calibration of the full multiphase EOS is arduous. We present our progress in development of automated EOS parameterization based on minimization of an objective function. Here, this function encodes deviation of model EOS results from experimental/computational benchmarks. Optimization is implemented as a combination of global (particle swarm) and local optimization techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗