Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear 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 919 records · Page 51

Sizing Energy Storage System for Energy Arbitrage in Extreme Fast Charging Station

This paper proposes a non-linear programming (NLP) model to optimally size the energy storage system (ESS) and obtain an optimal energy management for energy arbitrage of an extreme fast charging station (XFCS) for electric vehicles (EVs), with minimized total cost of XFCS operation and ESS investment. Different from most reported work on sizing the ESS for EV charging stations, this paper proposes a pragmatic approach to model the ESS life degradation and accurately count the ESS cycles. Moreover, this work incorporates the peak demand charges in the operational cost of the charging station which are often overlooked in the literature. The proposed model is formulated and solved using AIMMS. Finally, a thorough sensitivity analysis is performed to offer insights into how different input parameters impact the ESS sizing and savings from the energy arbitrage perspective.

25 ENERGY STORAGE↗

Optimal Regulation of Virtual Power Plants

This paper develops a real-time algorithmic framework for aggregations of distributed energy resources (DERs) in distribution networks to provide regulation services in response to transmission-level requests. Leveraging online primal-dual-type methods for time-varying optimization problems and suitable linearizations of the nonlinear AC power-flow equations, we believe this work establishes the system-theoretic foundation to realize the vision of distribution-level virtual power plants. The optimization framework controls the output powers of dispatchable DERs such that, in aggregate, they respond to automatic-generation-control and/or regulation-services commands. This is achieved while concurrently regulating voltages within the feeder and maximizing customers' and utility's performance objectives. Convergence and tracking capabilities are analytically established under suitable modeling assumptions. Simulations are provided to validate the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Direct handling of equality constraints in multilevel optimization

In recent years there have been several hierarchic multilevel optimization algorithms proposed and implemented in design studies. Equality constraints are often imposed between levels in these multilevel optimizations to maintain system and subsystem variable continuity. Equality constraints of this nature will be referred to as coupling equality constraints. In many implementation studies these coupling equality constraints have been handled indirectly. This indirect handling has been accomplished using the coupling equality constraints' explicit functional relations to eliminate design variables (generally at the subsystem level), with the resulting optimization taking place in a reduced design space. In one multilevel optimization study where the coupling equality constraints were handled directly, the researchers encountered numerical difficulties which prevented their multilevel optimization from reaching the same minimum found in conventional single level solutions. The researchers did not explain the exact nature of the numerical difficulties other than to associate them with the direct handling of the coupling equality constraints. The coupling equality constraints are handled directly, by employing the Generalized Reduced Gradient (GRG) method as the optimizer within a multilevel linear decomposition scheme based on the Sobieski hierarchic algorithm. Two engineering design examples are solved using this approach. The results show that the direct handling of coupling equality constraints in a multilevel optimization does not introduce any problems when the GRG method is employed as the internal optimizer. The optimums achieved are comparable to those achieved in single level solutions and in multilevel studies where the equality constraints have been handled indirectly.

Renaud, John E.↗

Measure this, not that: Optimizing the cost and model-based information content of measurements

Model-based design of experiments (MBDoE) is a powerful framework for selecting and calibrating science-based mathematical models from data. Here, this work extends popular MBDoE workflows by proposing a convex mixed integer (non)linear programming (MINLP) to optimize the selection of measurements. The solver MindtPy is modified to support calculating the D-optimality objective and its gradient via an external package, scipy, using the grey-box module in Pyomo. The new approach is demonstrated in two case studies: estimating highly correlated kinetics from a batch reactor and estimating transport parameters in a large-scale rotary packed bed for CO 2 capture. Both case studies show how examining the Pareto optimal trade-offs between information content measured by A- and D-optimality versus measurement budget offers practical guidance for selecting measurements for scientific experiments.

97 MATHEMATICS AND COMPUTING↗

Comparative Evaluation of Different Optimization Algorithms for Structural Design Applications

Non-linear programming algorithms play an important role in structural design optimization. Fortunately, several algorithms with computer codes are available. At NASA Lewis Research Centre, a project was initiated to assess the performance of eight different optimizers through the development of a computer code CometBoards. This paper summarizes the conclusions of that research. CometBoards was employed to solve sets of small, medium and large structural problems, using the eight different optimizers on a Cray-YMP8E/8128 computer. The reliability and efficiency of the optimizers were determined from the performance of these problems. For small problems, the performance of most of the optimizers could be considered adequate. For large problems, however, three optimizers (two sequential quadratic programming routines, DNCONG of IMSL and SQP of IDESIGN, along with Sequential Unconstrained Minimizations Technique SUMT) outperformed others. At optimum, most optimizers captured an identical number of active displacement and frequency constraints but the number of active stress constraints differed among the optimizers. This discrepancy can be attributed to singularity conditions in the optimization and the alleviation of this discrepancy can improve the efficiency of optimizers.

Patnaik, Surya N.↗

Minimum feature size control in level set topology optimization via density fields

A level set topology optimization approach that uses an auxiliary density field to nucleate holes during the optimization process and achieves minimum feature size control in optimized designs is explored. The level set field determines the solid-void interface and the density field describes the distribution of a fictitious porous material using the solid isotropic material with penalization. These fields are governed by two sets of independent optimization variables which are initially coupled using a penalty for hole nucleation. The strength of the density field penalization and projection is gradually increased during the optimization process to promote a 0-1 density distribution. In addition, a second penalty regulates the evolution of the density field in the void phase. The treatment of the density field combined with the second penalty mitigate the appearance of small design features. The minimum feature size of optimized designs is controlled by the radius of the linear filter applied to the density optimization variables. The structural response is predicted by the extended finite element method, the sensitivities by the adjoint method, and the optimization variables are updated by a gradient-based optimization algorithm. Numerical examples investigate the robustness of this approach with respect to algorithmic parameters and mesh refinement. The results show the applicability of the combined density level set topology optimization approach for both optimal hole nucleation and for minimum feature size control in 2D and 3D. This comes, however, at the cost of a more complex problem formulation and additional computational cost due to an increased number of optimization variables.

42 ENGINEERING↗

Resiliency of Distribution Systems Incorporating Asynchronous Information for System Restoration

Resiliency of distribution systems under extreme operating conditions is critical, especially when the utility is not available. With the large-scale deployment of distributed resources, it becomes possible to restore critical loads with local non-utility resources. Distribution system operators (DSOs) need to determine the critical loads to be restored, considering limited resources and distribution facilities. Several studies on resiliency have been conducted for the restoration of distribution systems. However, the inherent asynchronous characteristic on the information availability has not been incorporated. With incomplete and asynchronous information, decisions may be made that result in underutilization of generation resources. In this paper, a new distribution system restoration approach is proposed, considering uncertain devices and associated asynchronous information. It uses a two-module architecture that efficiently optimizes restoration actions using a binary linear programming model and evaluates their feasibility with unbalanced optimal power flow. Networked microgrids are included in the model. The IEEE 123-node test feeder is used for validation. Results show that asynchronous messages may affect the restoration actions significantly and the impacts can be mitigated by the proposed decision support tool for the DSOs.

Bedoya, Juan C.↗

Design of High-Accuracy Multiple Flyby Trajectories Using Constrained Optimization

The trajectory optimization technique described in this paper provides several distinct advantages over previous formulations. First, fully numerically integrated trajectory modeling is used. That is, no approximations to the trajectory are made and the inclusion of any level of complicated force models desired is allowed. Second, only trajectory propagation is used so there is no requirement for optimization. This is accomplished by the novel method of splitting the trajectory into independent legs, which are then subjected to constrained optimization. Third, each of the trajectory legs may be specified by any convenient set of parameters particularly useful for that leg. Any of these parameters may then be subject to constraints. Fourth, the nonlinear optimization problem is solved by solving a sequence of linear problems which converges to the optimal nonlinear solution. Fifth, the robustness of this formulation requires little or no user interaction with the optimization once a feasible problem has been posed.

flyby↗

Spacecraft Mission Design for the Mitigation of the 2017 PDC Hypothetical Asteroid Threat

This paper presents a detailed mission design analysis results for the 2017 Planetary Defense Conference (PDC) Hypothetical Asteroid Impact Scenario, documented at https://cneos.jpl.nasa.gov/ pd/cs/pdc17/. The mission design includes campaigns for both reconnaissance (flyby or rendezvous) of the asteroid (to characterize it and the nature of the threat it poses to Earth) and mitigation of the asteroid, via kinetic impactor deflection, nuclear explosive device (NED) deflection, or NED disruption. Relevant scenario parameters are varied to assess the sensitivity of the design outcome, such as asteroid bulk density, asteroid diameter, momentum enhancement factor, spacecraft launch vehicle, and mitigation system type. Different trajectory types are evaluated in the mission design process from purely ballistic to those involving optimal midcourse maneuvers, planetary gravity assists, and/or lowthrust solar electric propulsion. The trajectory optimization is targeted around peak deflection points that were found through a novel linear numerical technique method. The optimization process includes constrain parameters, such as Earth departure date, launch declination, spacecraft/asteroid relative velocity and solar phase angle, spacecraft dry mass, minimum/maximum spacecraft distances from Sun and Earth, and Earth/spacecraft communications line of sight. Results show that one of the best options for the 2017 PDC deflection is solar electric propelled rendezvous mission with a single spacecraft using NED for the deflection

Barbee, Brent W.↗

Spacecraft Mission Design for the Mitigation of the 2017 PDC Hypothetical Asteroid Threat

This paper presents a detailed mission design analysis results for the 2017 Planetary Defense Conference (PDC) Hypothetical Asteroid Impact Scenario, documented at https:cneos.jpl.nasa.govpdcspdc17. The mission design includes campaigns for both reconnaissance (flyby or rendezvous) of the asteroid (to characterize it and the nature of the threat it poses to Earth) and mitigation of the asteroid, via kinetic impactor deflection, nuclear explosive device (NED) deflection, or NED disruption. Relevant scenario parameters are varied to assess the sensitivity of the design outcome, such as asteroid bulk density, asteroid diameter, momentum enhancement factor, spacecraft launch vehicle, and mitigation system type. Different trajectory types are evaluated in the mission design process from purely ballistic to those involving optimal midcourse maneuvers, planetary gravity assists, and/or low-thrust solar electric propulsion. The trajectory optimization is targeted around peak deflection points that were found through a novel linear numerical technique method. The optimization process includes constrain parameters, such as Earth departure date, launch declination, spacecraft, asteroid relative velocity and solar phase angle, spacecraft dry mass, minimum/maximum spacecraft distances from Sun and Earth, and Earth-spacecraft communications line of sight. Results show that one of the best options for the 2017 PDC deflection is solar electric propelled rendezvous mission with a single spacecraft using NED for the deflection.

Barbee, Brent W.↗

The discrete complementary variational principle and optimal control systems

A discrete complementary variational principle is developed and applied to linear and nonlinear discrete-time optimal control systems. Using the variational approach, a primal-dual relationship is established. This relationship provides a precise measure of system suboptimality independent of any a priori knowledge of the optimal solution.

Chan, W. L.↗

A method to stabilize linear systems using eigenvalue gradient information

Formal optimization methods and eigenvalue gradient information are used to develop a stabilizing control law for a closed loop linear system that is initially unstable. The method was originally formulated by using direct, constrained optimization methods with the constraints being the real parts of the eigenvalues. However, because of problems in trying to achieve stabilizing control laws, the problem was reformulated to be solved differently. The method described uses the Davidon-Fletcher-Powell minimization technique to solve an indirect, constrained minimization problem in which the performance index is the Kreisselmeier-Steinhauser function of the real parts of all the eigenvalues. The method is applied successfully to solve two different problems: the determination of a fourth-order control law stabilizes a single-input single-output active flutter suppression system and the determination of a second-order control law for a multi-input multi-output lateral-directional flight control system. Various sets of design variables and initial starting points were chosen to show the robustness of the method.

Wieseman, C. D.↗

SU-E-T-592: Relationship Between Dose of Distribution and Area of Segment Fields Among Different Intensity-Modulated Radiotherapy Planning in Cervix Cancer

Purpose: In premise of uninfluenced to dose distribution of tumor target and organ at risk(OAR) in cervical cancer,area of segment fields was changed to increase efficacy and optimize treatment method by designing different plan of intensity modulated radiotherapy(IMRT). Methods: 12 cases of cervical cancer were confirmed in pathology and treated with step and shoot IMRT. Dose of PTV was 50Gy/25fractions. Every patient was designed 9 treatment plans of IMRT by Pinnacle 8.0m planning system,each plan was used with 9 beams of uniform distribution and fixing incidence direction(200°,240°,280°,320°,0°,40°,80°,120°and 160°respectively),and designed for delivery on Elekta Synergy linear accelerator. All plans were optimized with the direct machine parameter optimization(DMPO) algorithm using the same set of optimization objectives. Number of maximum segment field was defined at 80 and minimum MU in each segment was 5MU,and minimal segment area was 2*1cm{sup 2},2*2cm{sup 2},3*3cm{sup 2},4*4cm{sup 2},5*5cm{sup 2},6*6cm{sup 2},7*7cm{sup 2},8*8cm{sup 2}and 9*9cm{sup 2},respectively.Coverage,homogeneity and conformity of PTV,sparing of OAR, MU and number of segment were compared. Results: In this group, mean volume of PTV was 916.8±228.7 cm{sup 3}. Compared with the area of minimal segment field increased from 2*1cm{sup 2} to 9*9 cm{sup 2},the number of mean MU was decreased from 1405±170 to 490±47 and the number of segment field was reduced from 76±4 to 39±7 respectively(p<0.05). When the limit of minimal segment area was increased from 2*1cm{sup 2} to 7*7 cm{sup 2},dose distribution of PTV,OAR,CI,HI and V{sub 2} {sub 3} were not different (p>0.05),but when the minimal segment area was 8*8 cm{sup 2} and 9*9 cm{sup 2},they were changed compared with 7*7 cm{sup 2} and below(p<0.05). Conclusion: The minimal segment field of IMRT plan designed by Pinnacle 8.0m planning system in cervical carcinoma should be enlarge reasonably and minimal segment area of 7*7 cm{sup 2} was recommend.

60 APPLIED LIFE SCIENCES↗

Linear nanometric tunnel junction sensors with exchange pinned sensing layer

Highly sensitive nanosensors with high spatial resolution provide the necessary features for high accuracy imaging of isolated magnetic nanoparticles. In this work, we report the fabrication and characterization of MgO-barrier magnetic tunnel junction nanosensors, with two exchange-pinned electrodes. The perpendicular magnetization configuration for field sensing is set using a two-step annealing process, where the second annealing temperature was optimized to yield patterned sensors responses with improved linearity. The optimized circular nanosensors show sensitivities up to 0.1%/Oe, larger than previously reported for nanometric sensors and comparable to micrometric spin-valves. Our strategy avoids the use of external permanent biasing or demagnetizing fields (large for smaller structures) to achieve a linear response, enabling the control of the linear operation range using only the stack and thus providing a small footprint device.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Control-structure interaction

The feasibility of using conventional proportional-integral-derivative (PID) control and an alternative optimal control to perform the pointing and tracking functions of the Space Station solar dynamic power module is investigated. A very large state model of 6 rigid body modes and 272 flexible modes is used in conjunction with classical linear-quadratic-Gaussian (LQG) optimal control to produce a full-order controller that satisfies the requirements. The results are compared with a classically designed PID controller that was implemented for a much smaller (6 rigid body, 40 flexible modes) model. The conventional control design approach is shown to be very much influenced by the order reduction of the plant model, i.e., the number of retained elastic modes from the full-order model, suggesting that for a complex, large space structure, such as the Space Station Freedom solar dynamic module, application of conventional control system design methods may not be adequate. The use of LQG control is recommended, and method for solving the large matrix. Riccati equation that arises from the optimal formulation is provided.

Cheng, Joseph K.↗