SEARCH · Engineering Papers
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.
NREL researcher Peter Graf discusses distributed reinforcement learning with ADMM-RL as part of the Innovative Optimization and Control Methods for Highly Distributed Autonomous Systems workshop. NREL hosted the workshop to exchange ideas on state-of-the-art control and optimization strategies as well as get feedback on work being done at NREL.
The integration of aircraft maneuver characteristics into an optimal taxiway scheduling solution is challenging due to the uncertainties that are intrinsic to ramp area aircraft trajectories. To address the challenge, we build a stochastic model of ramp area aircraft trajectories that is used to generate a probabilistic measure of conflict within the Charlotte Douglas International Airport (CLT) ramp area. Parameters of the conflict distributions are estimated and passed to a Mixed Integer Linear Program that solves for an optimal taxiway schedule constrained to be conflict free in the presence of trajectory uncertainties. Here we extend our previous research by accounting for departing and arriving aircraft whereas our prior formulation only accounted for departing aircraft.
As large electricity consumers, water distribution system (WDS) pumping stations have the potential to become meaningful participants in demand response (DR) programs. The authors propose an optimization framework for assessing the DR capacity of a WDS and identifying the optimal bidding strategy for maximizing WDS revenue in the DR spot market. The proposed mixed integer linear programming (MILP) model overcomes computational constraints of previous DR optimization models by adopting a preprocessing procedure to minimize the number of binary variables and implementing a convex relaxation technique to linearize the hydraulic equations. The proposed MILP model also explicitly accounts for varying levels of risk tolerance of WDS operators by varying the recovery period over which pumping returns to business-as-usual operation. The optimization framework is implemented on a skeletonized 48-node WDS model that includes 7 pumps, 6 tanks, and 39 pipes. Using a simulated DR event and water consumption profile, the authors derive the optimal DR supply curves (i.e., compensation price versus load curtailment quantity) and revenue potential of the WDS under six scenarios for DR participation.
A method for optimally shifting the imaginary parts of the open-loop poles of a multivariable control system to the desirable closed-loop locations is presented. The optimal solution with respect to a quadratic performance index is obtained by solving a linear matrix Liapunov equation.
Applications for kinetic equations such as optimal design and inverse problems often involve finding unknown parameters through gradient-based optimization algorithms. Based on the adjoint-state method, we derive two different frameworks for approximating the gradient of an objective functional constrained by the nonlinear Boltzmann equation. While the forward problem can be solved by the DSMC method, it is difficult to efficiently solve the high-dimensional continuous adjoint equation obtained by the “optimize-then-discretize” approach. This challenge motivates us to propose an adjoint DSMC method following the “discretize-then-optimize” approach for Boltzmann-constrained optimization. We also analyze the properties of the two frameworks and their connections. Here, several numerical examples are presented to demonstrate their accuracy and efficiency.
A self-healing power grid can mitigate the impact of power outages and recover the affected area swiftly. To properly size the optimal amount and location of de-energized loads in the presence of both conventional and renewable sources, a precise optimization-based load recovery tool is proposed herein. The steady-state and transient behaviors of loads have been taken into account. The proposed tool contains three separate problems that interact with each other during the restoration period. The first-level problem determines the dispatching of generation units. The second-level problem comprises the non-linear algebraic equations, whose task is to provide the initial conditions for the third-level problem. Finally, the third-level problem intends to determine the optimal location and amount of load pickup through modeling the governor of generators as well as loads' dynamic behaviors. To examine the effectiveness of the proposed tool, a modified IEEE 39-system under a partial blackout condition has been tested. The results are verified by a commercial time-domain simulation software (PSS/E) to demonstrate the accuracy of the proposed approach. Ultimately, the impact of wind turbine frequency control on the load recovery process has been studied under different penetration levels.
Optimal design of hybrid renewable mini-grids requires both economic and power quality considerations. Existing modeling approaches address these considerations via separate or loosely coupled models. Here, we extend REopt—a techno-economic optimization model developed at the National Renewable Energy Laboratory—to consider both within a single model. REopt formulates the design problem as a mixed-integer linear program that solves for a site's optimal technology mix, sizing, and operation to minimize life cycle cost. REopt has traditionally assumed a single node system. In the work presented here, we expand the REopt platform to consider multiple connected nodes with associated voltage constraints. In order to do this, we model power flow using a fixed-point linear approximation method. Additionally, we then use the model to explore design considerations of mini-grids in Sub-Saharan Africa. Specifically, we evaluate under what combinations of transmission line distance and capacity it is technically viable and economically preferable to build multiple isolated mini-grids versus an interconnected, centralized system.
This paper studies the joint fleet sizing and charging system planning problem for a company operating a fleet of autonomous electric vehicles (AEVs) for passenger and goods transportation. Most of the relevant published papers focus on intracity scenarios and adopt heuristic approaches, e.g., agent based simulation, which do not guarantee optimality. In contrast, we propose a mixed integer linear programming model for intercity scenarios. This model incorporates comprehensive considerations of 1) limited AEV driving range; 2) optimal AEV routing and relocating operations; 3) time-varying origin-destination transport demands; and 4) differentiated operation cost structure of passenger and goods transportation. The proposed model can be computational expensive when the scale of the transportation network is large. We then exploit the structure of this program to expedite its solution. Numerical experiments are conducted to validate the proposed method. Our experimental results show that AEVs in passenger and goods transportation have remarkable planning and operation differences. We also demonstrate that intelligent routing and relocating operations, charging system and vehicle parameters, e.g., charging power, battery capacity, driving speed etc., can significantly affect the economic efficiency and the planning results of an AEV fleet.
Aspects of optimization theory and switching theory are discussed, taking into account the necessary conditions for extrema, a solution subject to constraints, the calculus of variations, the Pontryagin maximum principle, the canonical transformation, Contensou's domain of maneuverability, optimal switching, a junction with singular arc, and linearized singular control. Equations of motion are considered along with aerodynamic and propulsive forces, the general properties of optimal trajectories, flight in a horizontal plane, optimal coasting flight, supersonic cruise, the supersonic turn, supersonic maneuvers in a vertical plane, energy state approximation, a modified Chapman's formulation for optimal reentry trajectories, optimal planar reentry trajectories, and an optimal glide of reentry vehicles. Orbital aerodynamic maneuvers are examined, giving attention to aerodynamic capture, a change in the apogee, a change in the eccentricity, a change in the perigee, an orbital maneuver, an aerodynamic maneuver, and a combined maneuver.
Carbon Capture and Storage (CCS) is a critical technology for reducing anthropogenic CO2 emissions, but its large-scale deployment is complicated by uncertainties in geological storage performance. These uncertainties pose significant financial and operational risks, as underperforming storage sites can lead to costly infrastructure modifications, inefficient pipeline routing, and economic shortfalls. To address this challenge, we propose a novel optimization workflow that is based on mixed-integer linear programming and explicitly integrates probabilistic modeling of storage uncertainty into CCS infrastructure design. This workflow generates multiple infrastructure scenarios by sampling storage capacity distributions, optimally solving each scenario using a mixed-integer linear programming model, and aggregating results into a heatmap to identify core infrastructure components that have a low likelihood of underperforming. A risk index parameter is introduced to balance trade-offs between cost, CO2 processing capacity, and risk of underperformance, allowing stakeholders to quantify and mitigate uncertainty in CCS planning. Applying this workflow to a CCS dataset from the US Department of Energy’s Carbon Utilization and Storage Partnership project reveals key insights into infrastructure resilience. Reducing the risk index from 15% to 0% is observed to lead to an 83.7% reduction in CO2 processing capacity and a 77.1% decrease in project profit, quantifying the trade-off between risk tolerance and project performance. Furthermore, our results highlight critical breakpoints, where small adjustments in the risk index produce disproportionate shifts in infrastructure performance, providing actionable guidance for decision-makers. Unlike prior approaches that aimed to cheaply repair underperforming infrastructure, our workflow constructs robust CCS networks from the ground up, ensuring cost-effective infrastructure under storage uncertainty. These findings demonstrate the practical relevance of incorporating uncertainty-aware optimization into CCS planning, equipping decision-makers with a tool to make informed project planning decisions.
This paper presents a study of fuzzy optimal variance control problem for dynamical systems subject to actuator amplitude and rate constraints. Using Takagi-Sugeno fuzzy modeling and dynamic Parallel Distributed Compensation technique, the stability and the constraints can be cast as a multi-objective optimization problem in the form of Linear Matrix Inequalities. By utilizing the formulations and solutions for the input and output variance constraint problems, we develop a fuzzy full-state feedback controller. The stability and performance of the proposed controller is demonstrated through its application to the airfoil flutter suppression.
The NASA Engine Performance Program (NEPP) can configure and analyze almost any type of gas turbine engine that can be generated through the interconnection of a set of standard physical components. In addition, the code can optimize engine performance by changing adjustable variables under a set of constraints. However, for engine cycle problems at certain operating points, the NEPP code can encounter difficulties: nonconvergence in the currently implemented Powell's optimization algorithm and deficiencies in the Newton-Raphson solver during engine balancing. A project was undertaken to correct these deficiencies. Nonconvergence was avoided through a cascade optimization strategy, and deficiencies associated with engine balancing were eliminated through neural network and linear regression methods. An approximation-interspersed cascade strategy was used to optimize the engine's operation over its flight envelope. Replacement of Powell's algorithm by the cascade strategy improved the optimization segment of the NEPP code. The performance of the linear regression and neural network methods as alternative engine analyzers was found to be satisfactory. This report considers two examples-a supersonic mixed-flow turbofan engine and a subsonic waverotor-topped engine-to illustrate the results, and it discusses insights gained from the improved version of the NEPP code.
The linearity of the equations of motion of dynamical systems such as spacecraft or manipulators when controlled by acceleration commands is exploited in generating exact fast maneuvers. A simulated robot example thereof is examined, with advantages and limitations of the proposed design approach discussed in its context.
Idaho National Laboratory (INL) is developing a first-of-a-kind leadout instrumented capsule experiment design to enable in-situ measurement capabilities in the Advanced Test Reactor (ATR) core. The Ceramic Advanced Thermal Evolution Research (CRATER) experiment supports the aLEU program objective to accelerate fuel performance irradiation testing for identifying alternative high-assay, low enriched uranium (HALEU) fuel systems. CRATER is a fueled, instrumented capsule experiment to measure in-situ temperature and thermal conductivity of ceramic fuels. Two ceramic fuel types will be used, uranium mono-nitride (UN) and uranium mono-carbide (UC), with a third metallic fuel used for comparison (UMo). The three fuel specimens will use a stainless-steel cladding. Programmatic objectives include linear heat generation rates (LHGR) of 210 ± 25 Watts per cm. and an inner clad temperature of 300-450 °C. The evolution of fuel thermal conductivity during irradiation has never been successfully measured in-situ for these systems and this experiment is designed to use advances in measurement sciences to characterize how thermal transport properties evolve while in reactor. Neutronic simulations of the experiment and its surrounding reactor environment were conducted using the Monte Carlo N-Particle Transport code (MCNP) and result in optimized fuel enrichment to meet target linear heat generation rates (LHGRs) influencing fuel temperatures, and fuel burnup requirements. Fabrication research and development (R&D) efforts are underway to produce annular right cylinder UC and UN pellets using carbothermic reduction and nitridation (or hydride-dehydride-nitride) synthesis methods, followed double-action die cold isostatic pressing.
Robust optimization (RO) has attracted great attention from the optimization community over the past decade. RO is dedicated to solving optimization problems subject to uncertainty: design constraints must be satisfied for all the values of the uncertain parameters within a given uncertainty set. Uncertainty sets may be modeled as deterministic sets (boxes, polyhedra, ellipsoids), in which case the RO problem may be reformulated via worst-case analysis, or as families of distributions. The challenge of RO is to reformulate or approximate robust constraints so that the uncertain optimization problem is transformed into a tractable deterministic optimization problem. Most reformulation methods assume linearity of the robust constraints or uncertainty sets of favorable shape, which represents only a fraction of real-world applications. This survey addresses nonlinear RO and includes problem formulations and applications, solution approaches, and available software with code samples.
Recent severe outages highlight the urgency of improving grid resiliency in the U.S. Microgrid formation schemes are proposed to restore critical loads after outages occur. Most distribution networks have unbalanced configurations that are not represented in sufficient detail by single-phase models. This study provides a microgrid formation plan that adopts a three-phase network model to represent unbalanced distribution networks. The problem formulation has a quadratic objective function with mixed-integer linear constraints. The three-phase network model enables us to examine the three-phase power outputs of distributed generators (DGs), preventing unbalanced operation that might trip DGs. Because the DG unbalanced operation constraint is non-convex, an iterative process is presented that checks whether the unbalanced operation limits for DGs are satisfied after each iteration of optimization. We also develop a relatively conservative linear approximation on the unbalanced operation constraint to handle larger networks. Compared with the iterative solution process, the conservative linear approximation is able to accelerate the solution process at the cost of sacrificing optimality to a limited extent. Simulation in the IEEE 34 node and IEEE 123 test feeders indicate that the proposed method yields more practical microgrid formations results. In addition, this paper explores the coordinated operation of DGs and energy storage (ES) installations. The unbalanced three-phase outputs of ESs combined with the relatively balanced outputs of DGs could supply unbalanced loads. In conclusion, the case study also validates the DG-ES coordination.
This Technical Report presents work completed by The Applied Research Laboratory at The Pennsylvania State University, in conjunction with Sandia National Labs, on the optimization of the power conversion chain (PCC) design to maximize the Average Annual Electric Power (AAEP) output of an Oscillating Water Column (OWC) device. The design consists of two independent stages. First, the design of a floating OWC, a Backward Bent Duct Buoy (BBDB), and second the design of the PCC. The pneumatic power output of the BBDB in random waves is optimized through the use of a hydrodynamically coupled, linear, frequency-domain, performance model that links the oscillating structure to internal air-pressure fluctuations. The PCC optimization is centered on the selection and sizing of a Wells Turbine and electric power generation equipment. The optimization of the PCC involves the following variables: the type of Wells Turbine (fixed or variable pitched, with and without guide vanes), the radius of the turbine, the optimal vent pressure, the sizing of the power electronics, and number of turbines. Also included in this Technical Report are further details on how rotor thrust and torque are estimated, along with further details on the type of variable frequency drive selected.
The dawn of quantum computing brings on a revolution in the way combinatorially complex power system problems such as Unit Commitment are solved. The Unit Commitment problem complexity is expected to increase in the future because of the trend toward the increase of penetration of intermittent renewables. Even though quantum computing has proven effective for solving a host of problems, its applications for power systems’ problems have been rather limited. Here, in this paper, a quantum unit commitment is innovatively formulated and the quantum version of the decomposition and coordination alternate direction method of multipliers (ADMM) is established. The above is achieved by devising quantum algorithms and by exploiting the superposition and entanglement of quantum bits (qubits) for solving subproblems, which are then coordinated through ADMM to obtain feasible solutions. The main contributions of this paper include: 1) the innovative development of a quantum model for Unit Commitment; 2) development of decomposition and coordination-supported framework which paves the way for the utilization of limited quantum resources to potentially solve the large-scale discrete optimization problems; 3) devising the novel quantum distributed unit commitment (QDUC) to solve the problem in a larger scale than currently available quantum computers are capable of solving. The QDUC results are compared with those from its classical counterpart, which validate the efficacy of quantum computing.