Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimization methods”

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

Navigator function for the conformal bootstrap

Current numerical conformal bootstrap techniques carve out islands in theory space by repeatedly checking whether points are allowed or excluded. We propose a new method for searching theory space that replaces the binary information "allowed"/"excluded" with a continuous "navigator" function that is negative in the allowed region and positive in the excluded region. Such a navigator function allows one to efficiently explore high-dimensional parameter spaces and smoothly sail towards any islands they may contain. The specific functions we introduce have several attractive features: they are well-defined in large regions of parameter space, can be computed with standard methods, and evaluation of their gradient is immediate due to an SDP gradient formula that we provide. The latter property allows for the use of efficient quasi-Newton optimization methods, which we illustrate by navigating towards the 3d Ising island.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Minimizing Cr-Evaporation from Balance Of Plant Components By Utilizing Cost-Effective Alumina-Forming Austenitic Steels

A solid oxide fuel cell (SOFC) is a clean and efficient energy conversion device. The development of intermediate-temperature SOFCs has made it preferable to use metallic interconnects (MICs) to greatly reduce the cost and significantly increase the efficiency compared to ceramic interconnect materials. However, gaseous chromium species will evaporate from the chromium-containing layer formed on the surface of commonly used MICs and balance of plant (BoP) components. Volatile chromium species have been shown to form solid deposits which poison the cathodes of SOFCs, causing drastic cell performance degradation and thereby limiting commercialization. In order to alleviate the Cr poisoning and achieve long-term high performance of SOFC stacks, various Al 2 O 3 -forming austenitic (AFA) stainless steels applied at different temperatures are evaluated in this project. Based on our Phase I results (500 h operation), it is shown that on the AFAs, an alumina-based protective layer forms under high temperature that is invulnerable to water vapor effects and suppresses the diffusion of chromium and manganese which can prevent the generation of spinels on the alloy surface. Therefore, 310S, OC4 and OC5 at 800 °C and 625, OC11 and OC11LZ at 900 °C are selected to be tested in the Phase II long-term operation. From the results, we find that besides the lower Cr evaporation rates and better oxidation resistance of AFAs than benchmark alloys after short-term (500 h) operation, AFAs also possess the sturdy and compact alumina layer after a long-term operation (5000 h). A protective oxide layer is of great importance for the long-term high-temperature operation of structural materials. The high-temperature oxidation behavior of AFAs and commercial 310S and 625 alloys in Air + 10% H 2 O at 800 °C and 900 °C is systematically investigated. Severe breakaway oxidation and minor spallation are observed for the 310S and 625 after the long-term operation, while the AFAs show high oxidation resistance. It is found that the formation of a continuous alumina layer could greatly prevent the volatilization of chromium vapors. In addition, the corresponding evolving models are discussed. Chromium evaporation from BoP components in high-temperature environment could severely deteriorate the electrochemical performance of SOFC. Several methods were applied to evaluate the Cr evaporation rates of BoP after 500 h exposure at 800 °C to 900 °C in air with 10% H 2 O. An optimal method was designed to exclude the effect of silicon (Si) deposits from quartz tube and sodium (Na) deposits from the sodium carbonate on the oxidation process and the chemical interaction between Cr gaseous species and alumina tube which could provide further quantitative correlation of the evaporated Cr species quantities and degradation rates of SOFC. Based on the great performance of AFAs after long-term operation, AFAs are assembled with Anode-supported cells (ASC) to investigate the the Cr deposition of anode-supported cell under a constant current density of 0.5 A cm -2 at 800 °C with AFA alloys compared with commercial alloys. In addition, the anodic and cathodic processes are deconvoluted by distribution of relaxation times (DRT) method which are comprehensively discussed. It is found out that the voltage exhibited a slight decrease of 5.09 % and 1.54 % in the presence of OC11 and OC11LZA alloy, respectively. However, a considerable decrease of 22.14 % and 12.06 % was determined in the presence of 310S and 625 alloy, respectively. AFAs possessing the low-cost, low Cr evaporation rates and the high oxidation resistance will make it of great potential to replace the existing BoP components.

01 COAL, LIGNITE, AND PEAT↗

Minimizing Cr-Evaporation from Balance of Plant Components by Utilizing Cost-Effective Alumina-Forming Austenitic Steels

A solid oxide fuel cell (SOFC) is a clean and efficient energy conversion device. The development of intermediate-temperature SOFCs has made it preferable to use metallic interconnects (MICs) to greatly reduce the cost and significantly increase the efficiency compared to ceramic interconnect materials. However, gaseous chromium species will evaporate from the chromium-containing layer formed on the surface of commonly used MICs and balance of plant (BoP) components. Volatile chromium species have been shown to form solid deposits which poison the cathodes of SOFCs, causing drastic cell performance degradation and thereby limiting commercialization. In order to alleviate the Cr poisoning and achieve long-term high performance of SOFC stacks, various Al2O3-forming austenitic (AFA) stainless steels applied at different temperatures are evaluated in this project. Based on our Phase I results (500 h operation), it is shown that on the AFAs, an alumina-based protective layer forms under high temperature that is invulnerable to water vapor effects and suppresses the diffusion of chromium and manganese which can prevent the generation of spinels on the alloy surface. Therefore, 310S, OC4 and OC5 at 800 °C and 625, OC11 and OC11LZ at 900 °C are selected to be tested in the Phase II long-term operation. From the results, we find that besides the lower Cr evaporation rates and better oxidation resistance of AFAs than benchmark alloys after short-term (500 h) operation, AFAs also possess the sturdy and compact alumina layer after a long-term operation (5000 h). A protective oxide layer is of great importance for the long-term high-temperature operation of structural materials. The high-temperature oxidation behavior of AFAs and commercial 310S and 625 alloys in Air + 10% H 2 O at 800 °C and 900 °C is systematically investigated. Severe breakaway oxidation and minor spallation are observed for the 310S and 625 after the long-term operation, while the AFAs show high oxidation resistance. It is found that the formation of a continuous alumina layer could greatly prevent the volatilization of chromium vapors. In addition, the corresponding evolving models are discussed. Chromium evaporation from BoP components in high-temperature environment could severely deteriorate the electrochemical performance of SOFC. Several methods were applied to evaluate the Cr evaporation rates of BoP after 500 h exposure at 800 °C to 900 °C in air with 10% H 2 O. An optimal method was designed to exclude the effect of silicon (Si) deposits from quartz tube and sodium (Na) deposits from the sodium carbonate on the oxidation process and the chemical interaction between Cr gaseous species and alumina tube which could provide further quantitative correlation of the evaporated Cr species quantities and degradation rates of SOFC. Based on the great performance of AFAs after long-term operation, AFAs are assembled with Anode-supported cells (ASC) to investigate the the Cr deposition of anode-supported cell under a constant current density of 0.5 A cm -2 at 800 °C with AFA alloys compared with commercial alloys. In addition, the anodic and cathodic processes are deconvoluted by distribution of relaxation times (DRT) method which are comprehensively discussed. It is found out that the voltage exhibited a slight decrease of 5.09 % and 1.54 % in the presence of OC11 and OC11LZA alloy, respectively. However, a considerable decrease of 22.14 % and 12.06 % was determined in the presence of 310S and 625 alloy, respectively. AFAs possessing the low-cost, low Cr evaporation rates and the high oxidation resistance will make it of great potential to replace the existing BoP components.

36 MATERIALS SCIENCE↗

Green’s functions on a renormalized lattice: An improved method for the integer quantum Hall transition

Highlights: • Speed up of numerical calculations by exact renormalization group steps. • Practical approach to handle delicate scaling analysis. • Absence of evidence to prefer marginal corrections over power-law corrections. We introduce a performance-optimized method to simulate localization problems on bipartite tight-binding lattices. It combines an exact renormalization group step to reduce the sparseness of the original problem with the recursive Green’s function method. We apply this framework to investigate the critical behavior of the integer quantum Hall transition of a tight-binding Hamiltonian defined on a simple square lattice. In addition, we employ an improved scaling analysis that includes two irrelevant exponents to characterize the shift of the critical energy as well as the corrections to the dimensionless Lyapunov exponent. We compare our findings with the results of a conventional implementation of the recursive Green’s function method, and we put them into broader perspective in view of recent development in this field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-F-T-336: A Quick Auto-Planning (QAP) Method for Patient Intensity Modulated Radiotherapy (IMRT)

Purpose: The aim of this study is to develop a quick auto-planning system that permits fast patient IMRT planning with conformal dose to the target without manual field alignment and time-consuming dose distribution optimization. Methods: The planning target volume (PTV) of the source and the target patient were projected to the iso-center plane in certain beameye- view directions to derive the 2D projected shapes. Assuming the target interior was isotropic for each beam direction boundary analysis under polar coordinate was performed to map the source shape boundary to the target shape boundary to derive the source-to-target shape mapping function. The derived shape mapping function was used to morph the source beam aperture to the target beam aperture over all segments in each beam direction. The target beam weights were re-calculated to deliver the same dose to the reference point (iso-center) as the source beam did in the source plan. The approach was tested on two rectum patients (one source patient and one target patient). Results: The IMRT planning time by QAP was 5 seconds on a laptop computer. The dose volume histograms and the dose distribution showed the target patient had the similar PTV dose coverage and OAR dose sparing with the source patient. Conclusion: The QAP system can instantly and automatically finish the IMRT planning without dose optimization.

60 APPLIED LIFE SCIENCES↗

Dynamic ADMM for Real-Time Optimal Power Flow

This paper considers distribution networks featuring distributed energy resources (DERs), and develops a dynamic optimization method to maximize given operational objectives in real time while adhering to relevant network constraints. The design of the dynamic algorithm is based on suitable linearization of the AC power flow equations, and it leverages the so-called alternating direction method of multipliers (ADMM). The steps of the ADMM, however, are suitably modified to accommodate appropriate measurements from the distribution network and the DERs. With the aid of these measurements, the resultant algorithm can enforce given operational constraints in spite of inaccuracies in the representation of the AC power flows, and it avoids ubiquitous metering to gather the state of noncontrollable resources. Optimality and convergence of the proposed algorithm are established in terms of tracking of the solution of a convex surrogate of the AC optimal power flow problem.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Hydropower Optimization Using Artificial Neural Network Surrogate Models of a High-Fidelity Hydrodynamics and Water Quality Model

Hydropower operations optimization subject to environmental constraints is limited by challenges associated with dimensionality and spatial and temporal resolution. The need for high-fidelity hydrodynamic and water quality models within optimization schemes is driven by improved computational capabilities, increased requirements to meet specific points of compliance with greater resolution, and the need to optimize operations of not just single reservoirs but systems of reservoirs. This study describes an important advancement for computing hourly power generation schemes for a hydropower reservoir using high-fidelity models, surrogate modeling techniques, and optimization methods. The predictive power of the high-fidelity hydrodynamic and water quality model CE-QUAL-W2 is successfully emulated by an artificial neural network, then integrated into a genetic algorithm optimization approach to maximize hydropower generation subject to constraints on dam operations and water quality. This methodology is applied to a multipurpose reservoir near Nashville, Tennessee, USA. The model successfully reproduced high-fidelity reservoir information while enabling 6.8% and 6.6% increases in hydropower production value relative to actual operations for dissolved oxygen (DO) limits of 5 and 6 mg/L, respectively, while witnessing an expected decrease in power generation at more restrictive DO constraints. Exploration of simultaneous temperature and DO constraints revealed capability to address multiple water quality constraints at specified locations. Here, the reduced computational requirements of the new modeling approach demonstrated an ability to provide decision support for reservoir operations scheduling while maintaining high-fidelity hydrodynamic and water quality information as part of the optimization decision support routines.

13 HYDRO ENERGY↗

Modeling design and control problems involving neural network surrogates

Here, we consider nonlinear optimization problems that involve surrogate models represented by neural networks. We demonstrate first how to directly embed neural network evaluation into optimization models, highlight a difficulty with this approach that can prevent convergence, and then characterize stationarity of such models. We then present two alternative formulations of these problems in the specific case of feedforward neural networks with ReLU activation: as a mixed-integer optimization problem and as a mathematical program with complementarity constraints. For the latter formulation we prove that stationarity at a point for this problem corresponds to stationarity of the embedded formulation. Each of these formulations may be solved with state-of-the-art optimization methods, and we show how to obtain good initial feasible solutions for these methods. We compare our formulations on three practical applications arising in the design and control of combustion engines, in the generation of adversarial attacks on classifier networks, and in the determination of optimal flows in an oil well network.

97 MATHEMATICS AND COMPUTING↗

Manipulating Two-Photon Absorption of Molecules through Efficient Optimization of Entangled Light

We report how the unique temporal and spectral features of pulsed entangled photons from a parametric downconversion source can be utilized for manipulating electronic excitations through the optimization of their spectral phase. A new comprehensive optimization protocol based on Bayesian optimization has been developed in this work to selectively excite electronic states accessible by two-photon absorption. Using our optimization method, the entangled two-photon absorption probability for a thiophene dendrimer can be enhanced by up to a factor of 20, while classical light turns out to be nonoptimizable. Moreover, the optimization involving photon entanglement enables selective excitation that would not be possible otherwise. In addition to optimization, we have explored entangled two-photon absorption in the small entanglement time limit showing that entangled light can excite molecular electronic states that are vanishingly small for classical light. As a result, we demonstrate these opportunities with an application to a thiophene dendrimer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiphysics Design Optimization and Additive Manufacturing of Nuclear Components (Final CRADA Report - Executive Summary)

Westinghouse Electric Company (WEC) actively participated in the advancement of the nuclear fuel and reactor design space and requested the help of Oak Ridge National Laboratory (ORNL) in the creation of a new design tool set. This report details the creation of a collection of software tool sets that are linked together to collectively assist WEC design engineers in developing novel ideas outside the normal scope of traditional nuclear fuel and reactor design formulas. Specifically, Siemens HEEDS, a design space exploration and parametric optimization software, monitored and changed parameters in a collection of softwares to meet the team’s objective. The HEEDS parametric optimization method, SHERPA, was developed to control the Siemens NX CAD platform to adjust the native CAD of a hexahedral spacer grid. This new geometry can be used to execute a topological design optimization by the NX Topology software add-in. The resulting geometry is additively manufacturable. This topological optimization occurred twice—once on the spacer grid’s spring, and once on the dimple geometry. These new geometries were imported by Siemens’ STAR-CCM+, a multiphysics structural and fluid dynamic computational solver in which the spring geometry is deflected to match the rod insertion configuration. Along with the dimple geometry, this new deflected spring was used to complete a hydraulic assessment of a single-unit cell comprising one rod, one spring, and two dimples. The HEEDS SHERPA algorithm ranks the design based on the final mass of the unit cell and the hydraulic pressure drop performance. The ORNL team demonstrated the ability to use this software and provided engineering judgement to apply modern aerospace aerodynamic design. The effort has been focused on thinking outside the conventional design space and redesigning a spacer grid to perform beyond the WEC set objectives. Furthermore, the ORNL team also demonstrated that the HEEDS optimization routine can independently develop a design that meets the WEC design goals. Although these designs were at a low technology readiness level, their demonstration confirmed the team’s capability to create novel advanced nuclear concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Dynamic ADMM for Real-Time Optimal Power Flow: Preprint

This paper considers distribution networks featuring distributed energy resources (DERs), and develops a dynamic optimization method to maximize given operational objectives in real time while adhering to relevant network constraints. The design of the dynamic algorithm is based on suitable linearizations of the AC power flow equations, and it leverages the so-called alternating direction method of multipliers (ADMM). The steps of the ADMM, however, are suitably modified to accommodate appropriate measurements from the distribution network and the DERs. With the aid of these measurements, the resultant algorithm can enforce given operational constraints in spite of inaccuracies in the representation of the AC power flows, and it avoids ubiquitous metering to gather the state of non-controllable resources. Optimality and convergence of the propose algorithm are established in terms of tracking of the solution of a convex surrogate of the AC optimal power flow problem.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Lax-Oleinik-Type Formulas and Efficient Algorithms for Certain High-Dimensional Optimal Control Problems

Two of the main challenges in optimal control are solving problems with state-dependent running costs and developing efficient numerical solvers that are computationally tractable in high dimension. In this paper, we provide analytical solutions to certain optimal control problems whose running cost depends on the state variable and with constraints on the control. We also provide Lax-Oleinik-type representation formulas for the corresponding Hamilton-Jacobi partial differential equations with state-dependent Hamiltonians. Additionally, we present an efficient, grid-free numerical solver based on our representation formulas, which is shown to scale linearly with the state dimension, and thus, to overcome the curse of dimensionality. Using existing optimization methods and the min-plus technique, we extend our numerical solvers to address more general classes of convex and nonconvex initial costs. We demonstrate the capabilities of our numerical solvers using implementations on a central processing unit (CPU) and a field-programmable gate array (FPGA). In several cases, our FPGA implementation obtains over a 10 times speedup compared to the CPU, which demonstrates the promising performance boosts FPGAs can achieve. Furthermore, our numerical results show that our solvers have the potential to serve as a building block for solving broader classes of high-dimensional optimal control problems in real-time.

97 MATHEMATICS AND COMPUTING↗

Development of multi-objective core optimization framework and application to sodium-cooled fast test reactors

The optimization of a Sodium-cooled Fast Reactor (SFR) core is a challenging process, due to the large number of design parameters, the nonlinearities among inputs and outputs, and the complicated correlation among output parameters. This study attempts to develop a generalized framework for the SFR core optimization by coupling the sensitivity analysis, advanced optimization algorithm, and optionally the surrogate modeling. The framework is built based on the fast reactor modeling capability of the Argonne Reactor Computation (ARC) suite and the sensitivity analysis and optimization modules embedded in the DAKOTA code, both have been integrated within the NEAMS Workbench. The genetic algorithm is selected as the optimization method for its robustness, while the option of surrogate modeling is also explored to alleviate the computational burden caused by employing the ARC direct core physics simulation and thus enhance the efficiency of the optimization. Finally, the normalized deviations of performance parameters of the near-optimal solution from those of the ideal core are calculated and used as criteria to down select the final core design. The developed framework is applied to the Advanced Burner Test Reactor (ABTR) core, and optimal solutions are determined by balancing various objectives simultaneously, including peak fast flux, core volume, power, reactivity swing, plutonium mass feed, while at the same time satisfying the predefined constraints due to safety or economics considerations. The optimal ABTR core design obtained using the direct physical simulation and surrogated model are compared and discussed. It is found that using the accurately constructed surrogate models can significantly reduce the required computational time while maintaining satisfactory accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Convex Relaxation of Grid-Connected Energy Storage System Models With Complementarity Constraints in DC OPF

Including complementarity constraints in energy storage system (ESS) models in optimization problems ensure an optimal solution will not produce a physically unrealizable control strategy where there is simultaneous charging and discharging. However, the current approaches to impose complementarity constraints require the use of non-convex optimization methods. Here, we propose a convex relaxation for a common ESS model that has terms for both charging and discharging based on a penalty reformulation for use in a model predictive control (MPC) based optimal power flow (DC OPF) problem. In this approach, the complementarity constraints are omitted and a penalty term is added to the optimization objective function. For the DC OPF problem, we provide analysis for the conditions under which the convex relaxation of the complementarity constraint ensures that a solution with simultaneous ESS charging and discharging operation is suboptimal. Simulation results demonstrating ESS behavior with and without the penalty reformulation are provided for an MPC-based DC OPF problem on multiple IEEE test systems.

25 ENERGY STORAGE↗

First-principles study of structural stability, electronic, optical and elastic properties of binary intermetallic: PtZr

Structural, electronic, optical and elastic properties of PtZr have been studied using the full-potential linearized augmented plane wave (FP-LAPW) method within density functional theory (DFT). The energy against volume and enthalpy vs. pressure variation in three different structures i.e. B{sub 1}, B{sub 2} and B{sub 3} for PtZr has been presented. The equilibrium lattice parameter, bulk modulus and its pressure derivative have been obtained using optimization method for all the three phases. Furthermore, electronic structure was discussed to reveal the metallic character of the present compound. The linear optical properties are also studied under zero pressure for the first time. Results on elastic properties are obtained using generalized gradient approximation (GGA) for exchange correlation potentials. Ductile nature of PtZr compound is predicted in accordance with Pugh’s criteria.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Grain Size and Phase Purity Characterization of U 3 Si 2 Pellet Fuel

Characterization of U 3 Si 2 fresh fuel pellets is important for quality assurance and validation of the finished product. Grain size measurement methods, phase identification methods using scanning electron microscopes equipped with energy dispersive spectroscopy and x-ray diffraction, and phase quantification methods via image analysis have been developed and implemented on U 3 Si 2 pellet samples. A wide variety of samples have been characterized including representative pellets from an initial irradiation experiment, and samples produced using optimized methods to enhance phase purity from an extended fabrication effort. The average grain size for initial pellets was between 16 and 18 µm. The typical average grain size for pellets from the extended fabrication was between 20 and 30 µm with some samples exhibiting irregular grain growth. Pellets from the latter half of extended fabrication had a bimodal grain size distribution consisting of coarsened grains (>80 µm) surrounded by the typical (20-30 µm) grain structure around the surface. Phases identified in initial uranium silicide pellets included: U 3 Si 2 as the main phase composing about 80 vol. %, Si rich phases (USi and U 5 Si 4 ) composing about 13 vol. %, and UO 2 composing about 5 vol. %. Initial batches from the extended U 3 Si 2 pellet fabrication had similar phases and phase quantities. The latter half of the extended fabrication pellet batches did not contain Si rich phases, and had between 1-5% UO 2 : achieving U 3 Si 2 phase purity between 95 vol. % and 98 vol. % U 3 Si 2 . The amount of UO 2 in sintered U 3 Si 2 pellets is correlated to the length of time between U 3 Si 2 powder fabrication and pellet formation. These measurements provide information necessary to optimize fabrication efforts and a baseline for future work on this fuel compound.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Optimization of a Nucleic Acids united-RESidue 2-Point model (NARES-2P) with a maximum-likelihood approach

Coarse-grained models are useful tools to investigate the structural and thermodynamic properties of biomolecules. They are obtained by merging several atoms into one interaction site. Such simplified models try to capture as much as possible information of the original biomolecular system in all-atom representation but the resulting parameters of these coarse-grained force fields still need further optimization. In this paper, a force field optimization method, which is based on maximum-likelihood fitting of the simulated to the experimental conformational ensembles and least-squares fitting of the simulated to the experimental heat-capacity curves, is applied to optimize the Nucleic Acid united-RESidue 2-point (NARES-2P) model for coarse-grained simulations of nucleic acids recently developed in our laboratory. The optimized NARES-2P force field reproduces the structural and thermodynamic data of small DNA molecules much better than the original force field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗