Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimality”

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 1,567 records · Page 87

Monotonically convergent optimization in quantum control using Krotov's method

The non-linear optimization method developed by A. Konnov and V. Krotov [Autom. Remote Cont. (Engl. Transl.) 60, 1427 (1999)] has been used previously to extend the capabilities of optimal control theory from the linear to the non-linear Schroedinger equation [S. E. Sklarz and D. J. Tannor, Phys. Rev. A 66, 053619 (2002)]. Here we show that based on the Konnov-Krotov method, monotonically convergent algorithms are obtained for a large class of quantum control problems. It includes, in addition to nonlinear equations of motion, control problems that are characterized by non-unitary time evolution, nonlinear dependencies of the Hamiltonian on the control, time-dependent targets, and optimization functionals that depend to higher than second order on the time-evolving states. We furthermore show that the nonlinear (second order) contribution can be estimated either analytically or numerically, yielding readily applicable optimization algorithms. We demonstrate monotonic convergence for an optimization functional that is an eighth-degree polynomial in the states. For the ''standard'' quantum control problem of a convex final-time functional, linear equations of motion and linear dependency of the Hamiltonian on the field, the second-order contribution is not required for monotonic convergence but can be used to speed up convergence. We demonstrate this by comparing the performance of first- and second-order algorithms for two examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Inverse Biot–Savart Optimization for Superconducting Accelerator Magnets

Superconducting (SC) magnets for accelerator concepts are often synthesized by numerically optimizing magnetic field waveforms, a process that requires a subsequent solution of a constrained inverse problem to identify suitable SC magnet windings. When the desired field distribution is intuitive, the inverse process is facilitated by seeding preconceived coil distributions into design optimization methods for refinement. With more complex magnetic field distributions, an initial design may be unknown, and topology optimization tools are required to synthesize current distributions without a priori guidance from a subject matter expert. In this work, we develop a constrained inverse Biot-Savart topology optimization methodology that synthesizes optimal distributions of current density in racetrack-like SC coils. The problem structure is exploited through a computationally efficient quadratic programming formulation, and the method is applied to recently published magnetic field waveforms for a recirculating proton phase shifter, a proton therapy gantry, and dipole magnets with sharp field transitions. The method and results herein identify novel winding configurations that can help magnet designers bring accelerator concepts to fruition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Bidding strategy for microgrid in day-ahead market based on hybrid stochastic/robust optimization

In this paper, we propose an optimal bidding strategy in the day-ahead market of a microgrid consisting of intermittent distributed generation (DG), storage, dispatchable DG and price responsive loads. The microgrid coordinates the energy consumption or production of its components and trades electricity in both the day-ahead and real-time markets to minimize its operating cost as a single entity. The bidding problem is challenging due to a variety of uncertainties, including power output of intermittent DG, load variation, day-ahead and real-time market prices. A hybrid stochastic/robust optimization model is proposed to minimize the expected net cost, i.e., expected total cost of operation minus total benefit of demand. This formulation can be solved by mixed integer linear programming. The uncertain output of intermittent DG and day-ahead market price are modeled via scenarios based on forecast results, while a robust optimization is proposed to limit the unbalanced power in real-time market taking account of the uncertainty of real-time market price. Numerical simulations on a microgrid consisting of a wind turbine, a PV panel, a fuel cell, a micro-turbine, a diesel generator, a battery and a responsive load show the advantage of stochastic optimization in addition to robust optimization.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

REopt: A Platform for Energy System Integration and Optimization: Preprint

REopt is NREL's energy planning platform offering concurrent, multi-technology integration and optimization capabilities to help clients meet their cost savings and energy performance goals. The REopt platform provides techno-economic decision-support analysis throughout the energy planning process, from agency-level screening and macro planning to project development to energy asset operation. REopt employs an integrated approach to optimizing a site?s energy costs by considering electricity and thermal consumption, resource availability, complex tariff structures including time-of-use, demand and sell-back rates, incentives, net-metering, and interconnection limits. Formulated as a mixed integer linear program, REopt recommends an optimally-sized mix of conventional and renewable energy, and energy storage technologies; estimates the net present value associated with implementing those technologies; and provides the cost-optimal dispatch strategy for operating them at maximum economic efficiency. The REopt platform can be customized to address a variety of energy optimization scenarios including policy, microgrid, and operational energy applications. This paper presents the REopt techno-economic model along with two examples of recently completed analysis projects.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Binary Optimal Control of Single-Flux-Quantum Pulse Sequences

We introduce a binary, relaxed gradient, trust-region method for optimizing pulse sequences for single flux quanta (SFQ) control of a quantum computer. The pulse sequences are optimized with the goal of realizing unitary gate transformations. Each pulse has a fixed amplitude and duration. Here we model this process as an binary optimal control problem, constrained by Schrödinger’s equation, where the binary variables indicate whether each pulse is on or off. We introduce a first-order trust-region method, which takes advantage of a relaxed gradient to determine an optimal pulse sequence that minimizes the gate infidelity, while also suppressing leakage to higher energy levels. The proposed algorithm has a computational complexity of O(p log(p)), where p is the number of pulses in the sequence. We present numerical results for the H and X gates, where the optimized pulse sequences give gate fidelity’s better than 99.9%, in ≈ 25 trust-region iterations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimizing an Integrated Renewable-Electrolysis System

Hydrogen is a versatile energy carrier that is used in a wide variety of chemical and industrial processes. Producing hydrogen using electrolysis can enable integration of multiple sectors including electricity, heating, and industrial sectors; however, the cost of producing hydrogen from electrolysis remains a challenge for encouraging greater adoption. With growing amounts of renewable generation on the California grid, there is downward pressure on wholesale electricity prices, particularly during the afternoon from photovoltaics (PV). These lower, or even potentially negative prices, challenge the business cases for new and existing PV plants. In addition, as the grid transitions to less flexible generation, there is greater need for system flexibility. To help improve the economics for both solar PV and hydrogen production using electrolyzers, we explore the benefit of combining PV and electrolysis systems. The optimal breakeven hydrogen production cost for six unique market participation configurations is calculated at six candidate locations where PV is already installed. The six market configurations include islanded, separated, retail, net energy metering (NEM), hybrid retail/wholesale and wholesale. Using the Revenue Operation and Device Optimization Model (RODeO) model, the optimal breakeven hydrogen price over the lifetime of the equipment is calculated. The cost includes production, storage, and compression in preparation for gaseous delivery trucks. Revenue streams include the sale of hydrogen, low carbon fuel standard (LCFS) credits, renewable electricity sold to the grid, and Renewable Energy Credits (REC). The costs included are the electricity costs, capital and fixed operation and maintenance cost (FOM) for the electrolyzer, PV, and storage and compression systems as well as taxes and financing costs. In addition, cost reductions are achieved through retail and wholesale rate optimization, by which electricity is purchased at the lowest price and sold, if possible, at the highest price. For all locations, the breakeven hydrogen production cost results show that, in the order of decreasing cost, the system configurations are islanded (highest), separated, NEM, retail, hybrid retail/wholesale, and wholesale (lowest). The resulting system design balances between the capital and maintenance cost components, the operation costs (i.e., electricity costs) and the additional market revenues. The integration of solar PV and electrolysis is shown to provide a mutually beneficial relationship. For PV, integration with electrolysis offers the potential to hedge against wholesale market price volatility, and integration with electrolysis may offer the potential to defer or avoid transmission investment to deliver power to the point-of-use and instead use it on-site. When compared with SMR without considering any renewable hydrogen premiums, this study finds that PV + Electrolysis systems with current costs are likely not competitive; however, with cost reductions for electrolysis equipment consistent with DOE projections, it was found that systems with wholesale market access would be competitive, largely on account of both low capital costs and low-cost electricity. The electrolysis units can provide greater flexibility than is required based on retail rate optimization, so there is an opportunity for a utility or CAISO to increase system flexibility with PV + Electrolysis systems in return for commensurate compensation. In this way, there are potentially several solutions that fall between the hybrid configuration and the wholesale configuration that could provide sufficient compensation for a PV + Electrolysis unit to compete with SMR while also providing greater flexibility to the grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Effectively using multifidelity optimization for wind turbine design

Abstract. Wind turbines are complex multidisciplinary systems that are challenging to design because of the tightly coupled interactions between different subsystems. Computational modeling attempts to resolve these couplings so we can efficiently explore new wind turbine systems early in the design process. Low-fidelity models are computationally efficient but make assumptions and simplifications that limit the accuracy of design studies, whereas high-fidelity models capture more of the actual physics but with increased computational cost. This paper details the use of multifidelity methods for optimizing wind turbine designs by using information from both low- and high-fidelity models to find an optimal solution at reduced cost. Specifically, a trust-region approach is used with a novel corrective function built from a nonlinear surrogate model. We find that for a diverse set of design problems – with examples given in rotor blade geometry design, wind turbine controller design, and wind power plant layout optimization – the multifidelity method finds the optimal design using 38 %–58 % of the computational cost of the high-fidelity-only optimization. The success of the multifidelity method in disparate applications suggests that it could be more broadly applied to other wind energy or otherwise generic applications.

17 WIND ENERGY↗

Automated Multi-Robot Assembly of Compliance Optimized Structures

Autonomous assembly of large structures is one of the fundamental challenges on the way towards NASA’s objectives of deep space exploration. In this work, we propose an algorithmic framework to optimize the assembly process of a prescribed target structure by constraining the assembly effort as well as maintaining structural soundness throughout the process. This framework uses structural topology optimization with assembly effort metrics to generate checkpoints for robotic traversal algorithms. Assembly effort is quantified by the Wasserstein metric between consecutive structural configurations during the assembly process. The robotic assembly task is split into two subtasks, where we first optimize for a set of key frames, then perform reconfiguration between consecutive frames. Key frames are optimized by adopting topology optimization techniques to reduce assembly effort and maintain structural integrity during the assembly process, while reconfiguration between key frames is performed using a path planning algorithm with a minimum weight maximum matching approach on a bipartite graph. We employ a Crystalline robot model in which each structural element is capable of locomotion through the structure and locking into place with neighboring elements after reaching its destination. An example assembly of a two-dimensional cantilever beam under volume constraints and structural compliance considerations is presented to demonstrate the approach. Finally, we conclude by discussing possible future extensions to this work, including adoption of better metrics, extension to three-dimensional large-scale problems, and exacting finer control of structural integrity during the path-planning phase.

robotic assembly↗

Economic Dispatch Optimization of Multi-Unit SMR Site

With the increase of integration of variable renewable into the grid, grid resilience can be compromised without flexible electricity suppliers. Flexible operation of NPP can expand the scope of integration of new nuclear power plants in regions with high share of renewables. It also provides opportunity to increase profits through participation in ancillary services. In this research, NuScale modules are used as the case study. NuScale is an integrated pressurized water reactor with thermal power of 250 MWth. Having multiple units can provide higher flexibility in changing power levels and supplying the grid with electricity all year long. Additionally, load-following might perturbate the ordered sequence of refueling (every 24 months for NuScale). To maintain the 2 year refueling cycles, the loading patterns will need to be adjusted in a year-to-year basis, resulting in additional unnecessary costs. Another solution is to find the optimal outage schedule. This will also shift the focus of the planned outage to a more optimized predictive outages and maintenance. Results from short-term optimization (i.e., daily participation in energy and ancillary services markets) shows that average commitment to regulation down ancillary service markets during late night and early morning hours can increase revenues. Participating in upward reserves during pre-work commuting hours and early evening hours can also generate additional revenues due to high ramp up of electricity demand. Having multiple units at a single site can also help operate flexibly. Results from long-term optimization show that daily fluctuation in electric power is usually handled by few units while other units operate on base-load profile. It is possible to combine base-load operation with flexible operation to extend component remaining useful life in some units and find optimal outage schedule for the units.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Equivalent Area Targets for Inverse Design Optimization With Changing Low-Boom Cruise Conditions

This paper provides theoretical and numerical justifications for inverse design optimization of reversed equivalent area (A e,r ) of low-boom supersonic aircraft when the cruise condition changes during the optimization iterations. A modified linear theory for steady flow around a supersonic projectile is used to establish the accuracy of A e,r -based body-of-revolution approximation of a low-boom supersonic aircraft for undertrack sonic boom analysis using computational fluid dynamics off-body pressure. As a result, designing a low-boom shape of the off-body pressure at three body lengths below the aircraft is equivalent to designing a low-boom A e,r shape. A Bezier curve with eight control points is used to define an A e,r target for inverse design optimization of A e,r of a supersonic aircraft. To make an A e,r target matchable by a low-boom supersonic aircraft, the A e,r target must have the minimal perceived level of decibels (PLdB) for sonic boom and satisfy two constraints defined by the cruise condition, effective length, and two nondimensional parameters. Numerical results are used to verify that, for different cruise conditions, all optimized A e,r targets with PLdB below 70 can be approximately generated by scaling one of them if the targets have the same nondimensional parameters and approximately the same effective length.

multidisciplinary optimization↗

Exploiting Power Flow Manifold to Solve AC Optimal Power Flow

AC optimal power flow has proven difficult to solve with interior point methods on GPUs. This is largely due to challenging linear algebra problems that current state of the art massively parallel linear solvers struggle with. However, the advent of Riemannian optimization techniques and the fact that the power flow equations form a smooth manifold present an alternative approach. In this talk, we present the basics of Riemannian optimization techniques in which optimization is done directly on a manifold. Then we present computational results showing that Riemannian techniques are capable of producing solutions of comparable quality as interior point methods.

AC optimal power flow↗

Calibration of the E3SM Land Model Using Surrogate‐Based Global Optimization

Abstract Calibration of the Energy Exascale Earth System Model (E3SM), land model (ELMv0) is challenging because of its model complexity, strong model nonlinearity, and significant computational requirements. Therefore, only a limited number of simulations can be allowed in any attempt to find a near‐optimal solution within an affordable time. The goal of this study is to calibrate some of the ELMv0 parameters to improve model projection of carbon fluxes. We propose a computationally efficient global optimization procedure using sparse‐grid based surrogates. We first use advanced sparse grid (SG) interpolation to construct a surrogate system of the ELMv0, and then calibrate the surrogate model in the optimization process. As the surrogate model is a polynomial whose evaluation is fast, it can be efficiently evaluated a sufficiently large number of times in the optimization, which facilitates the global search. We calibrate eight parameters against five years of net ecosystem exchange, total leaf area index, and latent heat flux data from the U.S. Missouri Ozark flux tower. The calibrated model is then used for predicting the three variables in the following 4 years. The results indicate that an accurate surrogate model can be created for the ELMv0 with a relatively small number of SG points, i.e., a few ELMv0 simulations that can be fully parallel. And, the application of the optimized parameters leads to a better model performance and a higher predictive capability than the default parameter values in the ELMv0.

54 ENVIRONMENTAL SCIENCES↗

Selective Heating Through Y-Junction Waveguide Designed by Acoustic Shape Optimization

Unlike phononic crystals or systems designed by topology optimization, waveguides designed by shape optimization do not have voids or internal defects, making the fabrication process more suitable for additive manufacturing. By designing a Y-junction waveguide through shape optimization, an ultrasonic wave can be controlled so that it propagates to a predetermined location just by adjusting its frequency. These demultiplexed ultrasonic waves can be used to transport signals or stimulate nearby materials. As an example, the ultrasonic wave is converted to heat at different locations, which can be applied to mechanisms that can take advantage of heating. First, shape optimization is performed on a cylindrical structure to selectively propagate ultrasonic waves of a particular frequency while attenuating others, which is analyzed through a finite element model. The numerical study results are compared with experimental measurements from samples fabricated through additive manufacturing methods. After verifying the concept, the Y-junction waveguide is fabricated to demultiplex the wave and selectively heat different locations. Here the results show that the method of combining shape optimization with additive manufacturing is exceptionally simple and capable of demultiplexing ultrasonic waves, which can replace complex electrical components with single-material waveguides.

36 MATERIALS SCIENCE↗

The IDAES process modeling framework and model library—Flexibility for process simulation and optimization

Abstract Energy systems and manufacturing processes of the 21st century are becoming increasingly dynamic and interconnected, which require new capabilities to effectively model and optimize their design and operations. Such next generation computational tools must leverage state‐of‐the‐art techniques in optimization and be able to rapidly incorporate new advances. To address these requirements, we have developed the Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform, which builds on the strengths of both process simulators (model libraries) and algebraic modeling languages (advanced solvers). This paper specifically presents the IDAES Core Modeling Framework (IDAES‐CMF), along with a case study demonstrating the application of the framework to solve process optimization problems. Capabilities provided by this framework include a flexible, modifiable, open‐source platform for optimization of process flowsheets utilizing state‐of‐the‐art solvers and solution techniques, fully open and extensible libraries of dynamic unit operations models and thermophysical property models, and integrated support for superstructure‐based conceptual design and optimization under uncertainty.

42 ENGINEERING↗

Bayesian optimization to design a novel x-ray shaping device

In radiation therapy, x-ray dose must be precisely sculpted to the tumor, while simultaneously avoiding surrounding organs at risk. This requires modulation of x-ray intensity in space and/or time. Typically, this is achieved using a multi leaf collimator (MLC) - a complex mechatronic device comprising over one hundred individually powered tungsten ‘leaves’ that move in or out of the radiation field as required. Here, an all-electronic x-ray collimation concept with no moving parts is presented, termed “SPHINX”: Scanning Pencil-beam High-speed Intensity-modulated X-ray source. SPHINX utilizes a spatially distributed bremsstrahlung target and collimator array in conjunction with magnetic scanning of a high energy electron beam to generate a plurality of small x-ray “beamlets.” A simulation framework was developed in Topas Monte Carlo incorporating a phase space electron source, transport through user defined magnetic fields, bremsstrahlung x-ray production, transport through a SPHINX collimator, and dose in water. This framework was completely parametric, meaning a simulation could be built and run for any supplied geometric parameters. This functionality was coupled with Bayesian optimization to find the best parameter set based on an objective function which included terms to maximize dose rate for a user defined beamlet width while constraining inter-channel cross talk and electron contamination. Designs for beamlet widths of 5, 7, and 10 mm 2 were generated. Each optimization was run for 300 iterations and took approximately 40 h on a 24-core computer. For the optimized 7-mm model, a simulation of all beamlets in water was carried out including a linear scanning magnet calibration simulation. Finally, a back-of-envelope dose rate formalism was developed and used to estimate dose rate under various conditions. The optimized 5–, 7–, and 10-mm models had beamlet widths of 5.1 , 7.2 , and 10.1 mm 2 and dose rates of 3574, 6351, and 10 015 Gy/C, respectively. The reduction in dose rate for smaller beamlet widths is a result of both increased collimation and source occlusion. For the simulation of all beamlets in water, the scanning magnet calibration reduced the offset between the collimator channels and beam centroids from 2.9 ±1.9 mm to 0.01 ±0.03 mm. A slight reduction in dose rate of approximately 2% per degree of scanning angle was observed. Based on a back-of-envelope dose rate formalism, SPHINX in conjunction with next-generation linear accelerators has the potential to achieve substantially higher dose rates than conventional MLC-based delivery, with delivery of an intensity modulated 100 x 100 mm 2 field achievable in 0.9 to 10.6 s depending on the beamlet widths used. Bayesian optimization was coupled with Monte Carlo modeling to generate SPHINX geometries for various beamlet widths. A complete Monte Carlo simulation for one of these designs was developed, including electron beam transport of all beamlets through scanning magnets, x-ray production and collimation, and dose in water. These results demonstrate that SPHINX is a promising candidate for sculpting radiation dose with no moving parts, and has the potential to vastly improve both the speed and robustness of radiotherapy delivery. A multi-beam SPHINX system may be a candidate for delivering magavoltage FLASH RT in humans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Reliability‐based layout optimization in offshore wind energy systems

Abstract Existing methods for optimizing wind array layouts typically use power or cost objectives and rarely consider reliability‐based objectives. Component and system failure rates, however, are dependent on location‐specific wind conditions, are influenced by array layout and wake interactions, and have a direct and significant impact on capital costs, operational costs, and power production. Although wind power plant models exist that calculate wind loads with sufficient resolution to capture component loading dynamics from wind conditions, they are computationally expensive and thus not suitable for research applications requiring many evaluations, particularly optimization. This study describes the development of computationally efficient, reliability‐based layout optimization methods, enabling us to explore the relationship between component reliability and layout optimization. These methods include the surrogate modeling of the planet bearing life based on varying wind conditions simulated in FAST.Farm and the formulation of reliability‐based objectives based on failure cost and power production models. Through demonstration of this method, we explore how wind conditions, objective functions, and capacity density influence reliability‐based layout optimization. Results indicate that considering reliability alongside power production can reduce failure costs associated with replacement costs and downtime while maintaining or improving power production. Our conclusions highlight the opportunity for wind power plant developers to integrate reliability and operational expenditures alongside performance and capital expenditure objectives in plant design and development to improve plant performance and costs.

17 WIND ENERGY↗

Blade planform design optimization to enhance turbine wake control

Abstract This study considers optimizing the planform of wind turbine blades to ultimately enhance wind plant controls, namely, wake steering strategies. Adjoint‐enabled unsteady actuator line simulations are carried out to obtain gradients for optimization of several different performance objectives with respect to blade chord length at 10 locations along blade span. We demonstrate different blade design optimizations that can maximize time‐averaged lateral wake deflection, entrainment of kinetic energy, or total power of multiple turbines. Our optimized designs can produce a 4+% increase in wake deflection, a 4× increase in vertical kinetic energy entrainment, or a 3.6% increase in power when compared with the baseline case. While lateral wake deflection is only modestly sensitive to chord changes, we find that increasing the outboard chord length can dramatically increase kinetic energy entrainment, resulting in faster wake recovery and gains in net power. While this work develops only a few case studies emphasizing relative performance improvements and general trends, these results show the promise of a framework that combines mid‐fidelity computation with adjoint‐based optimization for control and design problems.

WIND ENERGY↗

Aero‐servo‐elastic co‐optimization of large wind turbine blades with distributed aerodynamic control devices

Abstract This work introduces automated wind turbine optimization techniques based on full aero‐servo‐elastic models and investigates the potential of trailing edge flaps to reduce the levelized cost of energy (LCOE) of wind turbines. The Wind Energy with Integrated Servo‐control (WEIS) framework is improved to conduct the presented research. Novel methods for the generic implementation and tuning of trailing edge flap devices and their controller are also introduced. Primary flap and controller parameters are optimized to demonstrate potential maximum blade tip deflection reductions of 21 % . Concurrent design optimization (i.e., co‐design) of a novel segmented wind turbine blade with trailing edge flaps and its controller is then conducted to demonstrate blade cost savings of 5 % . Additionally, rotor diameter co‐design optimization is demonstrated to reduce the LCOE by 1.3 % without significant load increases to the tower. These results demonstrate the efficacy of control co‐design optimization using trailing edge flaps, and the entirety of this work provides a foundation for numerous control co‐design‐oriented studies for distributed aerodynamic control devices.

17 WIND ENERGY↗