Search NASA⌕ Search

SEARCH · Search NASA

Results for “Partial 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 OSTI2023

Learning heterogeneous reaction kinetics from X-ray movies pixel-by-pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries 1 and electrocatalysts 2. Experimental characterizations of such materials by operando microscopy produce rich image datasets 3–6, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation 7. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

Chueh, William↗
DOE OSTI2026

Design of a GammaT-Jump System for Fermilab Booster

A $\gamma_{T}$ scheme may be required for the PIP-II era performance or ACE-MIRT era performance of the Booster. PIP-II era operations of the Fermilab proton complex will require the Fermilab Booster to increase beam intensity from 4.5e12 to 6.7e12 protons, while also increasing its ramp from 15 Hz to 20 Hz. These changes pose particular challenges for transition-crossing in the Booster, where longitudinal beam quality must be controlled in order to facilitate slip-stacking in the Recycler Ring later in the Main Injector cycle. Two novel $\gamma_{T}$-jump schemes are proposed, termed Double $\gamma_{T}$ jump and Partial $\gamma_{T}$-jump , which optimize the magnitude of the $\gamma_{T}$-jump within optics and power supply constraints.

Eldred, J. [Fermilab] (ORCID:000000034432072X)↗
DOE OSTI2023

Improving the Accuracy of Variational Quantum Eigensolvers with Fewer Qubits Using Orbital Optimization

Near-term quantum computers will be limited in the number of qubits on which they can process information as well as the depth of the circuits that they can coherently carry out. To date, experimental demonstrations of algorithms such as the Variational Quantum Eigensolver (VQE) have been limited to small molecules using minimal basis sets for this reason. In this work we propose incorporating an orbital optimization scheme into quantum eigensolvers wherein a parametrized partial unitary transformation is applied to the basis functions set in order to reduce the number of qubits required for a given problem. The optimal transformation is found by minimizing the ground state energy with respect to this partial unitary matrix. Through numerical simulations of small molecules up to 16 spin orbitals, we demonstrate that this method has the ability to greatly extend the capabilities of near-term quantum computers with regard to the electronic structure problem. Finally, we find that VQE paired with orbital optimization consistently achieves lower ground state energies than traditional VQE when using the same number of qubits and even frequently achieves lower ground state energies than VQE methods using more qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Machine Learning with Gradient-Based Optimization of Nuclear Waste Vitrification with Uncertainties and Constraints

Gekko is an optimization suite in Python that solves optimization problems involving mixed-integer, nonlinear, and differential equations. The purpose of this study is to integrate common Machine Learning (ML) algorithms such as Gaussian Process Regression (GPR), support vector regression (SVR), and artificial neural network (ANN) models into Gekko to solve data based optimization problems. Uncertainty quantification (UQ) is used alongside ML for better decision making. These methods include ensemble methods, model-specific methods, conformal predictions, and the delta method. An optimization problem involving nuclear waste vitrification is presented to demonstrate the benefit of ML in this field. ML models are compared against the current partial quadratic mixture (PQM) model in an optimization problem in Gekko. GPR with conformal uncertainty was chosen as the best substitute model as it had a lower mean squared error of 0.0025 compared to 0.018 and more confidently predicted a higher waste loading of 37.5 wt% compared to 34 wt%. The example problem shows that these tools can be used in similar industry settings where easier use and better performance is needed over classical approaches. Future works with these tools include expanding them with other regression models and UQ methods, and exploration into other optimization problems or dynamic control.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Hardware Design and Demonstration of a 100kW, 99% Efficiency Dual Active Half Bridge Converter Based on 1700V SiC Power MOSFET

High efficiency, high density and galvanically isolated power converters are attractive for numerous medium voltage high power applications. This paper presents a 1.5kVdc, 100kW bidirectional Dual Active Half Bridge (DAHB) using newly developed 1.7kV SiC MOSFET modules. The DAHB achieved an efficiency of 98.6% at full load of 100kW and a maximum of 99% at 45kW, operating as an isolated DC/DC converter. The converter also achieved an efficiency of 97.8% operating as an isolated DC/AC inverter. An optimized PCB-based busbar design has significantly reduced the voltage overshoot across the device, making the design suitable for 1500Vdc input application such as 1500V PV inverters. Typical partial discharge inception voltage (PDIV) of the optimized PCB busbar is 1.7kVpeak with total charge <10pC. The power density of the DAHB converter is 1.8MVA/m 3 which is much higher than traditional two-stage industry products. The DAHB converter can be used as a building block for a 4.16kV/1MVA utility scale PV systems with input parallel and output series configuration.

Xu, Wei↗
DOE OSTI2025

Monitoring Sulfuric Acid and Temperature Using Raman Spectroscopy and Multivariate Chemometrics

Multivariate regression models were optimized for the quantification of sulfuric acid (H 2 SO 4 ) [0–8 M] and temperature (20 °C–80 °C) in the presence of ammonium sulfate ((NH 4 ) 2 SO 4 [0–0.6 M]) using Raman spectroscopy. Optical vibrational spectroscopy is a useful nondestructive technique for the in situ analysis of complex chemical systems notoriously difficult to monitor in situ and in real-time. Multivariate analysis, a chemometrics method, can be paired with these nondestructive optical methods for determining analyte concentration and speciation in complex solutions, such as dissociated species in polyprotic acids, e.g., H 2 SO 4 . The effect of temperature is often overlooked although it can have a major influence on speciation and the corresponding Raman spectra. Here, in this study, partial least squares regression models were optimized for the quantification of H 2 SO 4 and its two deprotonated forms as a function of temperature. Measuring bisulfate as a function of temperature is particularly challenging owing to changes in the second dissociation constant. A designed training set effectively minimized the sample set size and trained a robust predictive model with percent root mean square error of <3% for H 2 SO 4 . The practical strategy employed here was demonstrated to be effective for building chemometric models that directly account for dynamic temperatures with static samples and is shown to be amenable to flow cell analysis applications with a simple calibration transfer for process monitoring applications.

D-optimal design↗
DOE OSTI2020

Analytical solution for partial heating on the exterior of the pipe with application to measuring fluid flow rate

An analytical model using Green’s functions for partial external heating of a pipe is developed, which results in an exact mathematical solution for the radial and axial temperature distribution in the pipe wall. Partial heating consists of a constant heat flux function imposed over a small section of the exterior of a pipe and for a limited time duration. The solution comprises steady-state and transient parts, and an algebraic identity is used to decrease the number of summation terms in the slowly-converging steady- state part. Intrinsic verification principles are used to verify the solution. As an example application, this transient solution is applied toward the development of a simple, noninvasive method for in-field measurement of the flow rate in pipes. Furthermore, to simulate this application, a pulse of energy is imposed to the wall of the pipe, and the developed mathematical solution is used to find the flow rate of the fluid inside the pipe. An optimal experiment is designed to find the best measurement location and time.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Plastic work constrained elastoplastic topology optimization

An elastoplastic topology optimization framework for limiting plastic work generation while maximizing stiffness is presented. The kinematics and constitutive model are based on finite strain linear isotropic hardening plasticity, and the balance laws are solved using a total Lagrangian finite element formulation. Aggregation of the specific plastic work combined with an adaptive normalization scheme efficiently constrains the maximum specific plastic work. The optimization problem is regularized using an augmented partial differential equation filter, and is solved by the method of moving asymptotes where path-dependent sensitivities are derived using the adjoint method. The numerical examples show a clear dependence on the optimized maximum stiffness structures for different levels of constrained specific plastic work. It is also shown that due to the history dependency of the plasticity, the load path significantly influences the structural performance and optimized topology.

42 ENGINEERING↗
DOE OSTI2022

Stochastic Learning Approach for Binary Optimization: Application to Bayesian Optimal Design of Experiments

Here, we present a novel stochastic approach to binary optimization suited for optimal experimental design (OED) for Bayesian inverse problems governed by mathematical models such as partial differential equations. The OED utility function, namely, the regularized optimality criterion, is cast into a stochastic objective function in the form of an expectation over a multivariate Bernoulli distribution. The probabilistic objective is then solved by using a stochastic optimization routine to find an optimal observational policy. This formulation (a) is generally applicable to binary optimization problems with soft constraints and is ideal for OED and sensor placement problems; (b) does not require differentiability of the original objective function (e.g., a utility function in OED applications) with respect to the design variable, and thus it enables direct employment of sparsity-enforcing penalty functions such as $\ell_0$, without needing to utilize a continuation procedure or apply a rounding technique; (c) exhibits much lower computational cost than traditional gradient-based relaxation approaches; and (d) can be applied to both linear and nonlinear OED problems with proper choice of the utility function. The proposed approach is analyzed from an optimization perspective with detailed convergence analysis of the optimization approach and is also analyzed from a machine learning perspective with correspondence to policy gradient reinforcement learning. The approach is demonstrated numerically by using an idealized two-dimensional Bayesian linear inverse problem and validated by extensive numerical experiments carried out for sensor placement in a parameter identification setup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Multi-objective, robust constraints enforced global topology optimizer for optical devices

A method for optimization of photonic devices is disclosed. The method includes receiving a set of unconstrained latent variables; mapping the set of unconstrained latent variables to a constrained space to generate a constrained device; calculating the permittivity across each element of the constrained device; determining a permittivity-constrained width gradient based at least partially on the permittivity across each element; and optimizing the set of unconstrained latent variables by at least partially using the permittivity-constrained width gradient.

Fan, Jonathan↗
DOE OSTI2021

Electrochemical reduction of acetonitrile to ethylamine

Electrifying chemical manufacturing using renewable energy is an attractive approach to reduce the dependence on fossil energy sources in chemical industries. Primary amines are important organic building blocks; however, the synthesis is often hindered by the poor selectivity because of the formation of secondary and tertiary amine byproducts. Herein, we report an electrocatalytic route to produce ethylamine selectively through an electroreduction of acetonitrile at ambient temperature and pressure. Among all the electrocatalysts, Cu nanoparticles exhibit the highest ethylamine Faradaic efficiency (~96%) at -0.29 V versus reversible hydrogen electrode. Under optimal conditions, we achieve an ethylamine partial current density of 846 mA cm -2 . A 20-hour stable performance is demonstrated on Cu at 100 mA cm -2 with an 86% ethylamine Faradaic efficiency. Moreover, the reaction mechanism is investigated by computational study, which suggests the high ethylamine selectivity on Cu is due to the moderate binding affinity for the reaction intermediates.

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

Application of helium line intensity ratio spectroscopy to xenon plasma in E × B Penning discharge

Here, we propose the application of helium line intensity ratio spectroscopy in a low-pressure (0.3 mTorr) xenon E × B discharge with an electron temperature of ~2 eV and a density of 10 10 –10 11 cm –3 . We successfully identified the helium atom line emissions at 388.9, 447.1, 501.6, 504.8, and 706.5 nm with helium pressures of up to ~20 mTorr. The measured electron temperature, density, and I–V characteristics of the discharge remained almost constant in all helium pressures in the present experiment, indicating the suitability of the helium gas as a diagnostic gas. The results of helium line intensity ratio spectroscopy using the line emissions at 388.9, 447.1, and 504.8 nm showed fair agreement with the Langmuir probe measurement. Considering the trade-off relationship between the disturbance introduced by the helium gas and the signal-to-noise ratio, we conclude that a helium pressure of approximately 4 mTorr (approximately 13 times the partial pressure of xenon) represents the optimal pressure range for the application of the helium line emission intensity ratio method to this xenon plasma. It is found that the use of the line emissions at 501.6 and 706.5 nm result in a significant disturbance in the helium line intensity ratio method due to the radiation trapping effect.

helium line intensity ratio spectroscopy↗
DOE OSTI2025

Dynamic Modeling, Trajectory Optimization, and Linear Control of Cable-Driven Parallel Robots for Automated Panelized Building Retrofits

The construction industry faces a growing need for automation to reduce costs, improve accuracy and productivity, and address labor shortages. One area that stands to benefit significantly from automation is panelized prefabricated building envelope retrofits, which can improve a building’s energy efficiency in heating and cooling interior spaces. In this paper, we propose using cable-driven parallel robots (CDPRs), which can effectively lift and handle large objects, to install these panels. However, implementing CDPRs presents significant challenges because of their nonlinear dynamics, complex trajectory planning, and precise control requirements. To tackle these challenges, this work focuses on a new application of established control and trajectory optimization theories in a CDPR simulation of a building envelope retrofit under real-world conditions. We first model the dynamics of CDPRs, highlighting the critical role of damping in system behavior. Building on this dynamic model, we formulate a trajectory optimization problem to generate feasible and efficient motion plans for the robot under operational and environmental constraints. Given the high precision required in the construction industry, accurately tracking the optimized trajectory is essential. However, challenges such as partial observability and external vibrations complicate this task. To address these issues, a Linear Quadratic Gaussian control framework is applied, enabling the robot to track the optimized trajectories with precision. Simulation results show that the proposed controller enables precise end effector positioning with errors under 4 mm, even in the presence of external wind disturbances. Through comprehensive simulations, our approach allows for an in-depth exploration of the system’s nonlinear dynamics, trajectory optimization, and control strategies under controlled yet highly realistic conditions. The results demonstrate the feasibility of CDPRs for automating panel installation and provide insights into their practical deployment.

CDPR↗
DOE OSTI2024

PETSc/TAO Users Manual V.3.21

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations (PDEs) and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for implementing large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and optimizers that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started ). The library is organized hierarchically, enabling users to employ the abstraction level most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Modeling and Characterization of Surface Discharges in Insulating Material for Spacers: Electrode Shape, Discharge Mode, and Revision of the Creepage Concept

In the design of MV AC and DC spacers, the predominant factors are surface and interface conditions. Design is generally carried out on specifications and standards which are based on long-term experience and lab testing. However, the diffusion of power electronics with a trend to increase electric field, switching frequency, and rise time to achieve higher power density calls for an innovative, global approach to optimized insulation system design. A new methodology, based on field simulation, discharge modeling, and partial discharge inception measurements, called the three-leg approach, can form the basis to optimize insulation design for any type of supply voltage waveform. This paper focuses on the influence of the type of electrode on the inception and phenomenology of surface discharges and, as a consequence, on the interpretation of the results used for application of the three-leg approach. It is demonstrated that a typical electrode system used for insulating material testing can generate both gas and surface discharges at the triple point, when the electrodes have a smooth profile that is used to avoid corona or flashover. Hence, testing partial discharge may not provide a straightforward indication of the surface discharge inception and, thus, be partially misleading for insulation design. Another takeover is that such analysis must benefit from PD testing tools endowed with analytics able to provide automatic identification of the type of defect generating PD, i.e., internal, surface, and corona, since design and remedy actions can be taken, and adequate insulating materials developed, only knowing the type of source generating PD. Hence, testing partial discharge may not provide a straightforward indication of surface discharge inception and, thus, be partially misleading for insulation design. In addition to the importance of the three-leg approach to favor reliable and optimized design of insulation systems, there is a clear need to have a PD testing tool endowed with analytics. It should preferably be able to provide automatic identification of the type of defect generating PD, i.e., internal, surface, and corona.

36 MATERIALS SCIENCE↗
DOE OSTI2022

On the energy landscape of symmetric quantum signal processing

Symmetric quantum signal processing provides a parameterized representation of a real polynomial, which can be translated into an efficient quantum circuit for performing a wide range of computational tasks on quantum computers. For a given polynomial f , the parameters (called phase factors) can be obtained by solving an optimization problem. However, the cost function is non-convex, and has a very complex energy landscape with numerous global and local minima. It is therefore surprising that the solution can be robustly obtained in practice, starting from a fixed initial guess &#x03A6; 0 that contains no information of the input polynomial. To investigate this phenomenon, we first explicitly characterize all the global minima of the cost function. We then prove that one particular global minimum (called the maximal solution) belongs to a neighborhood of &#x03A6; 0 , on which the cost function is strongly convex under the condition &#x2016; f &#x2016; &#x221E; = O ( d &#x2212; 1 ) with d = d e g ( f ) . Our result provides a partial explanation of the aforementioned success of optimization algorithms.

Wang, Jiasu↗
DOE OSTI2023

Simulation results for a low energy nuclear recoil yields measurement in liquid xenon using the MiX detector

Measuring the scintillation and ionization yields of liquid xenon in response to ultra-low energy nuclear recoil events is necessary to increase the sensitivity of liquid xenon experiments to light dark matter. Neutron capture on xenon can be used to produce nuclear recoil events with energies below 0.3 keV NR via the asymmetric emission of γ rays during nuclear de-excitation. The feasibility of an ultra-low energy nuclear recoil measurement using neutron capture was investigated for the Michigan Xenon (MiX) detector, a small dual-phase xenon time projection chamber that is optimized for a high scintillation gain. Simulations of the MiX detector, a partial neutron moderator, and a pulsed neutron generator indicate that a population of neutron capture events can be isolated from neutron scattering events. Additionally, the rate of neutron captures in the MiX detector was optimized by varying the thickness of the partial neutron moderator, neutron pulse width, and neutron pulse frequency.

detector calibration↗