Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fast 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 901 records · Page 50

Active operator learning with predictive uncertainty quantification for partial differential equations

With the increased prevalence of neural operators being used to provide rapid solutions to partial differential equations (PDEs), understanding the accuracy of model predictions and the associated error levels is necessary for deploying reliable surrogate models in scientific applications. Existing uncertainty quantification (UQ) frameworks employ ensembles or Bayesian methods, which can incur substantial computational costs during both training and inference. Here, we propose a lightweight predictive UQ method tailored for Deep operator networks (DeepONets) that also generalizes to other operator networks. Numerical experiments on linear and nonlinear PDEs demonstrate that the framework’s uncertainty estimates are unbiased and provide accurate out-of-distribution uncertainty predictions with a sufficiently large training dataset. Our framework provides fast inference and uncertainty estimates that can efficiently drive outer-loop analyses that would be prohibitively expensive with conventional solvers. We demonstrate how predictive uncertainties can be used in the context of Bayesian optimization and active learning problems to yield improvements in accuracy and data-efficiency for outer-loop optimization procedures. In the active learning setup, we extend the framework to Fourier Neural Operators (FNO) and describe a generalized method for other operator networks. To enable real-time deployment, we introduce an inference strategy based on precomputed trunk outputs and a sparse placement matrix, reducing evaluation time by more than a factor of five. Our method provides a practical route to uncertainty-aware operator learning in time-sensitive settings.

97 MATHEMATICS AND COMPUTING↗

A compact neutron spectrometer system

A portable neutron spectrometer, nSpec, has been developed for passive interrogation of fast neutrons ( 1 MeV) from neutron sources. The design of the system includes four EJ-301 liquid scintillation cells, a custom stand, and tungsten caps for gamma-ray shielding. Data acquisition settings were optimized to provide a dynamic range of up to > 15 MeV neutron energy. Additionally, the detector response matrix was characterized using a time-of-flight method based on gamma-ray tagged neutron emission from 252 Cf. The spectrum unfolding was characterized using measured and simulated spectra of AmBe and 252 Cf with the GRAVEL and MLEM methods. Monte Carlo methods were used to examine the uncertainty propagation and determine sufficient counts for a measurement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SCORE (Synthesis of COnsists as Rolling Energy microgrids) (Final Technical Report)

The SCORE toolset is an open source, web-based application to assess the impact of new powering technologies on railroad performance, specifically technologies capable of both putting power into the rail (motive force) and taking power from the rail (regeneration). SCORE’s primary goal is generate trade studies to analyze different powering and train make-up options to minimize energy usage and greenhouse gas generation. At the core of generating these trade spaces is calculating the optimal powering policy for the train given the makeup of the train, the route, and time constraints. This paper presents details on the algorithms used in SCORE to compute this powering policies that is fast and accurate, discusses its implementation in an Energy-Longitudinal Train Dynamics (E-LTD) model, compares it to naïve approaches, and demonstrates its use across a variety of train/route pairs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Optimization of Scrap Melting Using an Electric Arc in Steel Manufacturing

Steel industry is crucial to the national economy and security. Around 67% of crude steel in the U.S is produced in electric arc furnaces (EAF), which is energy intensive. Around 140 EAFs operate in the U.S., consuming about 8.6x10 7 MMBtu/year of electricity. One of major challenges for EAFs includes maximizing the efficiency of the electrical energy provided in the form of electric arcs to melt various scrap mixes. To address this issue, a computational fluid dynamics (CFD) methodology is chosen to analyze scrap melting using the electric arc. Due to complex furnace phenomena and the wide variety of potential scenarios, high performance computing (HPC) is essential to yield comprehensive and detailed CFD analyses and systematic parametric studies for optimized EAF operation. The objectives are to 1) simulate scrap melting using electric arc, 2) evaluate electrode/arc position for optimum scrap melting and 3) establish reduced order model for CFD data-base for fast model calculation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Factorization and the synthesis of optimal feedback kernels for differential-delay systems

A combination of ideas from the theories of operator Riccati equations and Volterra factorizations leads to the derivation of a novel, relatively simple set of hyperbolic equations which characterize the optimal feedback kernel for the finite-time regulator problem for autonomous differential-delay systems. Analysis of these equations elucidates the underlying structure of the feedback kernel and leads to the development of fast and accurate numerical methods for its computation. Unlike traditional formulations based on the operator Riccati equation, the gain is characterized by means of classical solutions of the derived set of equations. This leads to the development of approximation schemes which are analogous to what has been accomplished for systems of ordinary differential equations with given initial conditions.

Milman, Mark M.↗

Minimizing distortion and internal forces in truss structures by simulated annealing

Inaccuracies in the length of members and the diameters of joints of large truss reflector backup structures may produce unacceptable levels of surface distortion and member forces. However, if the member lengths and joint diameters can be measured accurately it is possible to configure the members and joints so that root-mean-square (rms) surface error and/or rms member forces is minimized. Following Greene and Haftka (1989) it is assumed that the force vector f is linearly proportional to the member length errors e(sub M) of dimension NMEMB (the number of members) and joint errors e(sub J) of dimension NJOINT (the number of joints), and that the best-fit displacement vector d is a linear function of f. Let NNODES denote the number of positions on the surface of the truss where error influences are measured. The solution of the problem is discussed. To classify, this problem was compared to a similar combinatorial optimization problem. In particular, when only the member length errors are considered, minimizing d(sup 2)(sub rms) is equivalent to the quadratic assignment problem. The quadratic assignment problem is a well known NP-complete problem in operations research literature. Hence minimizing d(sup 2)(sub rms) is is also an NP-complete problem. The focus of the research is the development of a simulated annealing algorithm to reduce d(sup 2)(sub rms). The plausibility of this technique is its recent success on a variety of NP-complete combinatorial optimization problems including the quadratic assignment problem. A physical analogy for simulated annealing is the way liquids freeze and crystallize. All computational experiments were done on a MicroVAX. The two interchange heuristic is very fast but produces widely varying results. The two and three interchange heuristic provides less variability in the final objective function values but runs much more slowly. Simulated annealing produced the best objective function values for every starting configuration and was faster than the two and three interchange heuristic.

Kincaid, Rex K.↗

An efficient and robust grid optimization algorithm

The development of an efficient and robust grid optimization is presented. This algorithm is developed by combining the best characteristics of algebraic, elliptic, and hyperbolic grid generation techniques. This development is based on the following observations and evaluations: (1) algebraic systems are fast and economical; (2) precise spacing control is always achieved; (3) grid generation by elliptic systems is always smooth; and (4) the hyperbolic system preserves the orthogonality at the solid boundary and the point distribution in the field. Computational examples representing practical internal flow configurations are presented to demonstrate the algorithm.

Soni, Bharat K.↗

Transient Vibration Prediction for Rotors on Ball Bearings Using Load-dependent Non-linear Bearing Stiffness

Rolling-element bearing forces vary nonlinearly with bearing deflection. Thus an accurate rotordynamic transient analysis requires bearing forces to be determined at each step of the transient solution. Analyses have been carried out to show the effect of accurate bearing transient forces (accounting for non-linear speed and load dependent bearing stiffness) as compared to conventional use of average rolling-element bearing stiffness. Bearing forces were calculated by COBRA-AHS (Computer Optimized Ball and Roller Bearing Analysis - Advanced High Speed) and supplied to the rotordynamics code ARDS (Analysis of Rotor Dynamic Systems) for accurate simulation of rotor transient behavior. COBRA-AHS is a fast-running 5 degree-of-freedom computer code able to calculate high speed rolling-element bearing load-displacement data for radial and angular contact ball bearings and also for cylindrical and tapered roller beatings. Results show that use of nonlinear bearing characteristics is essential for accurate prediction of rotordynamic behavior.

Fleming, David P.↗

Optimization with Neural Network Feasibility Surrogates: Formulations and Application to Security-Constrained Optimal Power Flow

In many areas of constrained optimization, representing all possible constraints that give rise to an accurate feasible region can be difficult and computationally prohibitive for online use. Satisfying feasibility constraints becomes more challenging in high-dimensional, non-convex regimes which are common in engineering applications. A prominent example that is explored in the manuscript is the security-constrained optimal power flow (SCOPF) problem, which minimizes power generation costs, while enforcing system feasibility under contingency failures in the transmission network. In its full form, this problem has been modeled as a nonlinear two-stage stochastic programming problem. In this work, we propose a hybrid structure that incorporates and takes advantage of both a high-fidelity physical model and fast machine learning surrogates. Neural network (NN) models have been shown to classify highly non-linear functions and can be trained offline but require large training sets. In this work, we present how model-guided sampling can efficiently create datasets that are highly informative to a NN classifier for non-convex functions. We show how the resultant NN surrogates can be integrated into a non-linear program as smooth, continuous functions to simultaneously optimize the objective function and enforce feasibility using existing non-linear solvers. Overall, this allows us to optimize instances of the SCOPF problem with an order of magnitude CPU improvement over existing methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Optimizing High Performance Markov Clustering for Pre-Exascale Architectures

HipMCL is a high-performance distributed memory implementation of the popular Markov Cluster Algorithm (MCL) and can cluster large-scale networks within hours using a few thousand CPU-equipped nodes. It relies on sparse matrix computations and heavily makes use of the sparse matrix-sparse matrix multiplication kernel (SpGEMM). The existing parallel algorithms in HipMCL are not scalable to Exascale architectures, both due to their communication costs dominating the runtime at large concurrencies and also due to their inability to take advantage of accelerators that are increasingly popular. In this work, we systematically remove scalability and performance bottlenecks of HipMCL. We enable GPUs by performing the expensive expansion phase of the MCL algorithm on GPU. Additionally, we propose a CPU-GPU joint distributed SpGEMM algorithm called pipelined Sparse SUMMA and integrate a probabilistic memory requirement estimator that is fast and accurate. Furthermore, we develop a new merging algorithm for the incremental processing of partial results produced by the GPUs, which improves the overlap efficiency and the peak memory usage. We also integrate a recent and faster algorithm for performing SpGEMM on CPUs. We validate our new algorithms and optimizations with extensive evaluations. With the enabling of the GPUs and integration of new algorithms, HipMCL is up to 12.4x faster, being able to cluster a network with 70 million proteins and 68 billion connections just under 15 minutes using 1024 nodes of ORNL's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗

Advanced electrolytes for fast-charging high-voltage lithium-ion batteries in wide-temperature range

LiNixMnyCo1-x-yO2 (NMC) cathode materials with Ni=0.8 have attracted great interest for application in high energy-density lithium (Li)-ion batteries (LIBs) because of their high specific capacities at high voltages. However, the practical application of Ni-rich NMC in LIBs under high charge voltages (e.g. 4.4 V and above) still faces big challenges due to the severe capacity fading, which is directly related to the instability of the cathode/electrolyte interface. In this work, we develop new localized high-concentration electrolytes (LHCEs) based on carbonate solvents, which are compatible with graphite (Gr) anode and have a high oxidation potential over 4.9 V vs. Li/Li+. The optimal LHCE enables the Gr||NMC811 cells with 2.8 mAh cm-2 cathode areal capacity loading to achieve an excellent cycling stability over 600 cycles with a capacity retention of 94.2 % under a charge cut-off voltage of 4.4 V at 25 ?C, good rate capabilities under fast charging or discharging up to 3C rate, and superior low-temperature discharge performance down to -30 ?C with a capacity retention of 85.6% at C/5 rate. The findings in this work shed light a very promising strategy to develop new electrolytes for practical high-energy LIBs with Ni-rich NMC cathodes.

Zhang, Xianhui↗

In-route inductive versus stationary conductive charging for shared automated electric vehicles: A university shuttle service

In-route inductive charging technology, as applied to automated electric vehicles, can help realize a fully automated system of both vehicles and chargers. This study presents a planning optimization analysis for fixed-route automated shuttles supported by in-route inductive charging technology. A techno-economic feasibility of inductive charging was assessed in comparison with stationary charging, including Level 2 AC chargers, and DC fast chargers (DCFCs). This analysis considered both present-day and future vehicle operations and overall system costs. A real project with two circulator Navya Arma shared automated electric vehicles (SAEVs) at the University of Michigan was investigated using real-world collected energy and travel data. The outcomes show that the proper design of quasi-dynamic inductive chargers at designated stops allows SAEVs to realize unlimited driving range and be cost-competitive to DCFC technology. Considering present-day costs and vehicles, low-speed SAEVs can realize charge-sustaining operation at a minimum cost either by implementing a 50-kW inductive charger at two stops with one segment per position and a 29-kWh onboard battery, or by installing a 100-kW inductive charger at one stop with one segment per position and a 28-kWh onboard battery. Overall, considering future costs and vehicles, either a 40-kW charger at one stop with a 29-kWh battery or a 50-kW charger at the north stop with a 14-kWh battery would enable charge-sustaining operation. In addition, quasi-dynamic inductive solution can reduce the onboard battery by about 15% while providing unlimited driving range, but stationary scenarios require about 112% additional battery capacity to support a 12-h driving range.

33 ADVANCED PROPULSION SYSTEMS↗

Collaborative Research: Improved Efficiency and Coupling of the Radiation Code in the ACME Earth System Model (Final Report)

The complexity of radiative transfer, its importance to the exchange of energy in the climate system, and its high computational cost establishes importance of an accurate and efficient radiative transfer parameterization for climate simulation. Previous work by the proposing team at AER led to the development of the radiation code, RRTMG, which has been widely accepted by the global modeling community as a fast and accurate advancement over the previous generation of radiation codes. It has been in use in the NCAR CESM for many years, and it has been implemented in the initial version the DOE E3SM model. However, its computational cost remains high relative to other components in part due to its complexity and to its inefficient use of modern optimization strategies, and this project helped address this limitation for the code’s application in E3SM. Under other funding, the Investigators of this project led an effort to develop a high-performance broadband radiation code, called RTE+RRTMGP, which is a completely restructured code that will take advantage of modern computational capabilities to enhance its performance while retaining the strengths and accuracy of the original code. Designed to perform over a range of computer architectures, RTE+RRTMGP makes extensive use of Fortran 2003 features to improve both its efficiency and its use of memory. RTE+RRTMGP is expected to be adopted in the next generation of RTE+RRTMGP. This project allowed for advancements to the code’s computational capabilities and adding new and enhanced features. This included revising RTE+RRTMGP to run on GPU processors, analysis of and improvements to the code’s timing, modifying the gas optics in RRTMGP to increase the code’s accuracy, extensive validation of computed fluxes, heating rates and forcings, generating a cloud optical property and vertical sampling capabilities for RTE+RRTMGP, implementing a fast and accurate longwave scattering capability, and groundwork for the inclusion of a capability to specify solar variability. Many of the accomplishment in this project necessitated significant collaboration with the E3SM development team. The result of this project was optimization of a key physical component (radiative transfer calculations) of E3SM, directly supporting E3SM’s overarching global modeling objectives. More broadly, this project provided overall advancements in the use of radiative transfer calculations in atmospheric modeling and simulation, particularly for climate.

58 GEOSCIENCES↗

Feasible region-based heuristics for optimal transmission switching

In this paper, we develop a optimal transmission switching (OTS) heuristic based on DC optimal power flow (OPF) and assess the efficacy of the approach when implemented within AC OPF. Traditional formulations of the OTS problem can result in hundreds or thousands of binary variables for large networks, making the OTS problem challenging to solve on fast timescales even for relatively small networks. Here, we identify which constraints and therefore which variables are constraining the DC OPF feasible region, and rank them based on their impact on the cost function. We develop a heuristic algorithm which iteratively removes these constraints and solves a series of standard DC OPF problems. The heuristic is tested on a variety of PGlib networks and the results show that the algorithm can provide substantial cost decreases without having to solve any mixed integer programs. Additionally, we provide insights about the OTS problem, including identifying scenarios outside congestion where OTS can prove useful. Lastly, the performance of the DC-based heuristic is shown when the line switching decisions are implemented within AC OPF.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Fast Automatic Knot Placement Method for Accurate B-spline Curve Fitting

The choice of knot vector has immense influence on the resulting accuracy of a B-spline approximation of a curve. However, despite the significance of this problem and the various solutions that were proposed in the literature, optimizing the number and placement of knots remains a difficult task. This paper presents a novel method for the approximation of a curve by a B-spline of arbitrary order, which automatically determines a knot vector that achieves high approximation quality. Additionally, at the core of our approach is a feature function that characterizes the amount and spatial distribution of geometric details in the input curve by estimating its derivatives. Knots are then selected in such a way as to evenly distribute the feature contents across their intervals. A comparison to the state of the art for a wide variety of curves shows that our method is faster and achieves more accurate reconstruction results, while typically reducing the number of necessary knots.

97 MATHEMATICS AND COMPUTING↗

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because it supports stiff fast dispersive waves. The preconditioner is formulated using an operator-split approximate block factorization (Schur complement), informed by physics insight. We use a vector-potential formulation (instead of a magnetic field one) to allow a clean segregation of the problematic $\nabla$ x $\nabla$ x operator in the electron Ohm's law subsystem. This segregation allows the formulation of an effective damped block-Jacobi smoother for multigrid. We demonstrate by analysis that our proposed block-Jacobi iteration is convergent and has the smoothing property. The resulting HMHD solver is verified linearly with wave propagation examples, and nonlinearly with the GEM challenge reconnection problem by comparison against another HMHD code. We demonstrate the excellent algorithmic and parallel performance of the algorithm up to 16384 MPI tasks in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enabling high areal capacity for Co-free high voltage spinel materials in next-generation Li-ion batteries

The rapidly growing technological demand for lithium-ion batteries has prompted the development of novel cathode materials with high energy density, low cost, and improved safety. High voltage spinel, LiNi 0.5 Mn 1.5 O 4 (LNMO), is one of the most promising candidates yet to be commercialized. The two primary obstacles for this material are the inferior electronic conductivity and fast capacity degradation in full cells due to the high operating voltage. By systematically addressing these limitations, we successfully develop a thick LNMO electrode with areal capacity loadings up to 3 mAh·cm –2 . The optimized thick electrode is paired with a commercial graphite anode at both the coin cell and pouch cell level, achieving a full cell capacity retention as high as 72% and 78%, respectively, after 300 cycles. We attribute this superior cycling stability to careful optimizations of cell components and testing conditions, with a specific focus improving electronic conductivity and high voltage compatibility. These results suggest precise control of materials quality, electrode architecture and electrolyte optimization can soon support the development of a cobalt-free battery system based on a thick LNMO cathode (>4 mAh·cm 2 ), which will eventually meet the needs of next-generation Li-ion batteries with reduced cost, improved safety, and assured sustainability.

25 ENERGY STORAGE↗

Molecular Design Considerations for Azobenzene Anolytes

Realization of high-density batteries requires the development of anolytes that display highly negative reduction potentials, solubility, and persistence in the charged state. Azobenzenes have garnered interest as potential anolytes for redox flow batteries. Here, we report the synthesis of a family of substituted azobenzene derivatives and evaluation of their solution-phase electrochemical properties. Systematic synthetic derivatization of this scaffold allows (1) access to anolytes of varying solubility, including intrinsically liquid derivatives that represent potential high-density charge carriers; (2) systematic variation of the reduction potential, and in some cases redox inventory, that provides azobenzenes with highly negative reduction potentials; and (3) control of the lifetime of the azobenzene radical anions that result from one-electron reduction. Electrokinetic experiments demonstrated that fast electron transfer occurs for all derivatives examined. Spectroscopic characterization of monoreduced azobenzene derivatives establishes that decomposition of the azobenzene radical anion proceeds via bimolecular disproportionation. Together, these results provide an experimental basis for the optimization of azobenzene anolytes for electrochemical storage applications, including redox flow batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗