Search NASA⌕ Search

SEARCH · Search NASA

Results for “quadratic 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 289 records · Page 16

Project Wish: The Emerald City, phase 3

Phase 3 of Project Wish saw the evolution of the Emerald City (E-City) from a collection of specialized independent analyses and ideas to a working structural design integrated with major support systems and analyses. Emphasis was placed on comparing and contrasting the closed and open cycle gas core nuclear rocket engines to further determine the optimum propulsive system for the C-City. Power and thermal control requirements were then defined and the question of how to meet these requirements was addressed. Software was developed to automate the mission/system/configuration analysis so changes dictated by various subsystems constraints could be managed efficiently and analyzed interactively. In addition, the liquid hydrogen propellant tank was statically designed for minimum mass and shape optimization using a finite element modeling package called SDRC I-DEAS while spoke and shaft cross-sectional areas were optimized on ASTROS (Automated Structural Optimization System). A structural dynamic analysis also conducted using ASTROS enabled a study of the displacements, accelerations, modes and frequencies of the C-City. Finally, the attitude control system design began with an initial mass moment of inertia analysis and was then designed and optimized using linear quadratic regulator control theory.

Source record↗

Active control of acoustic pressure fields using smart material technologies

An overview describing the use of piezoceramic patches in reducing noise in a structural acoustics setting is presented. The passive and active contributions due to patches which are bonded to an Euler-Bernoulli beam or thin shell are briefly discussed and the results are incorporated into a 2-D structural acoustics model. In this model, an exterior noise source causes structural vibrations which in turn lead to interior noise as a result of nonlinear fluid/structure coupling mechanism. Interior sound pressure levels are reduced via patches bonded to the flexible boundary (a beam in this case) which generate pure bending moments when an out-of-phase voltage is applied. Well-posedness results for the infinite dimensional system are discussed and a Galerkin scheme for approximating the system dynamics is outlined. Control is implemented by using linear quadratic regulator (LQR) optimal control theory to calculate gains for the linearized system and then feeding these gains back into the nonlinear system of interest. The effectiveness of this strategy for this problem is illustrated in an example.

Banks, H. T.↗

Project WISH: The Emerald City

Phase 3 of Project WISH saw the evolution of the Emerald City (E-City) from a collection of specialized independent analyses and ideas to a working structural design integrated with major support systems and analyses. Emphasis was placed on comparing and contrasting the closed and open cycle gas core nuclear rocket engines to further determine the optimum propulsive system for the E-City. Power and thermal control requirements were then defined and the question of how to meet these requirements was addressed. Software was developed to automate the mission/system/configuration analysis so changes dictated by various subsystem constraints could be managed efficiently and analyzed interactively. In addition, the liquid hydrogen propellant tank was statically designed for minimum mass and shape optimization using a finite element modeling package called SDRC I-DEAS. Spoke and shaft cross-sectional areas were optimized on ASTROS (Automated Structural Optimization System) for mass minimization. A structural dynamic analysis of the optimal structure also conducted using ASTROS enabled a study of the modes, frequencies, displacements, and accelerations of the E-City. Finally, the attitude control system design began with an initial mass moment of inertia analysis and was then designed and optimized using linear quadratic regulator control theory.

Oz, Hayrani↗

Shape Control of Plates with Piezo Actuators and Collocated Position/Rate Sensors

This paper treats the control problem of shaping the surface deformation of a circular plate using embedded piezo-electric actuator and collocated rate sensors. An explicit Linear Quadratic Gaussian (LQG) optimizer stability augmentation compensator is derived as well as the optimal feed-forward control. Corresponding performance evaluation formulas are also derived.

Balakrishnan, A. V.↗

Shape Control of Plates with Piezo Actuators and Collocated Position/Rate Sensors

This paper treats the control problem of shaping the surface deformation of a circular plate using embedded piezo-electric actuators and collocated rate sensors. An explicit Linear Quadratic Gaussian (LQG) optimizer stability augmentation compensator is derived as well as the optimal feed-forward control. Corresponding performance evaluation formulas are also derived.

Balakrishnan, A. V.↗

Preliminary Optimal Orbit Design for the Laser Interferometer Space Antenna (LISA)

In this paper we present a preliminary optimal orbit analysis for the Laser Interferometer Space Antenna (LISA). LISA is a NASA/ESA mission to study gravitational waves and test predictions of general relativity. The nominal formation consists of three spacecraft in heliocentric orbits at 1 AU and trailing the Earth by twenty degrees. This configuration was chosen as a trade off to reduce the noise sources that will affect the instrument and to reduce the fuel to achieve the final orbit. We present equations for the nominal orbit design and discuss several different measures of performance for the LISA formation. All of the measures directly relate the formation dynamics to science performance. Also, constraints on the formation dynamics due to spacecraft and instrument limitations are discussed. Using the nominal solution as an initial guess, the formation is optimized using Sequential Quadratic Programming to maximize the performance while satisfying a set of nonlinear constraints. Results are presented for each of the performance measures.

Hughes, Steven P.↗

Spacecraft Station-Keeping Trajectory and Mission Design Tools

Two tools were developed for designing station-keeping trajectories and estimating delta-v requirements for designing missions to a small body such as a comet or asteroid. This innovation uses NPOPT, a non-sparse, general-purpose sequential quadratic programming (SQP) optimizer and the Two-Level Differential Corrector (T-LDC) in LTool (Libration point mission design Tool) to design three kinds of station-keeping scripts: vertical hovering, horizontal hovering, and orbiting. The T-LDC is used to differentially correct several trajectory legs that join hovering points. In a vertical hovering, the maximum and minimum range points must be connected smoothly while maintaining the spacecrafts range from a small body, all within the law of gravity and the solar radiation pressure. The same is true for a horizontal hover. A PatchPoint is an LTool class that denotes a space-time event with some extra information for differential correction, including a set of constraints to be satisfied by T-LDC. Given a set of PatchPoints, each with its own constraint, the T-LDC differentially corrects the entire trajectory by connecting each trajectory leg joined by PatchPoints while satisfying all specified constraints at the same time. Vertical and horizontal hover both are needed to minimize delta-v spent for station keeping. A Python I/F to NPOPT has been written to be used from an LTool script. In vertical hovering, the spacecraft stays along the line joining the Sun and a small body. An instantaneous delta-v toward the anti- Sun direction is applied at the closest approach to the small body for station keeping. For example, the spacecraft hovers between the minimum range (2 km) point and the maximum range (2.5 km) point from the asteroid 1989ML. Horizontal hovering buys more time for a spacecraft to recover if, for any reason, a planned thrust fails, by returning almost to the initial position after some time later via a near elliptical orbit around the small body. The mapping or staging orbit may be similarly generated using T-LDC with a set of constraints. Some delta-v tables are generated for several different asteroid masses.

Chung, Min-Kun J.↗

QSPIN: A High Level Java API for Quantum Computing Experimentation

QSPIN is a high level Java language API for experimentation in QC models used in the calculation of Ising spin glass ground states and related quadratic unconstrained binary optimization (QUBO) problems. The Java API is intended to facilitate research in advanced QC algorithms such as hybrid quantum-classical solvers, automatic selection of constraint and optimization parameters, and techniques for the correction and mitigation of model and solution errors. QSPIN includes high level solver objects tailored to the D-Wave quantum annealing architecture that implement hybrid quantum-classical algorithms [Booth et al.] for solving large problems on small quantum devices, elimination of variables via roof duality, and classical computing optimization methods such as GPU accelerated simulated annealing and tabu search for comparison. A test suite of documented NP-complete applications ranging from graph coloring, covering, and partitioning to integer programming and scheduling are provided to demonstrate current capabilities.

Quantu↗

Pterodactyl: Development and Comparison of Control Architectures for a Mechanically Deployed Entry Vehicle

The Pterodactyl project, seeks to advance the current state-of-the-art for entry vehicles by developing novel guidance and control technologies for Deployable Entry Vehicles (DEVs) that can be applied to various entry vehicle configurations. This paper details the efforts on the NASA-funded Pterodactyl project to investigate and implement multiple control techniques for an asymmetric mechanical DEV. We design multiple control architectures for a Pterodactyl Baseline Vehicle (PBV) and evaluate their performance in achieving varying guidance commands during entry. The control architectures studied are (i) propulsive control systems such as reaction control systems and (ii) non-propulsive control systems such as aerodynamic control surfaces and internal moving masses. For each system, state-feedback integral controllers based on linear quadratic regulator (LQR) optimal control methods are designed to track guidance commands of either (i) bank angle or (ii) angle of attack and sideslip angle as determined by the desired guidance trajectory. All control systems are compared for a lunar return reference mission and by providing a comparative analysis of these systems, configurations, and performance, the efforts detailed in this paper and the Pterodactyl project as a whole will help entry vehicle system designers determine suitable control architectures for integration with DEVs and other entry vehicle types.

Okolo, Wendy A.↗

Pterodactyl: Guidance and Control of a Symmetric Deployable Entry Vehicle using an Aerodynamic Control System

The NASA-funded Pterodactyl project seeks to advance the state-of-the-art for varying entry vehicle types by developing unconventional guidance and control technologies for Deployable Entry Vehicles (DEVs) that can be applied to different entry vehicle configurations. Prior work by the authors [1–5] involved developing both traditional and novel integrated guidance and control solutions for a Pterodactyl Baseline Vehicle (PBV), a variant of an asymmetric DEV called the Lifting Nano ADEPT (LNA) [6]. In the prior studies, two different guidance schemes were designed and implemented for the PBV: (i) traditional bank angle guidance developed using the Fully Numerical Predictor-Corrector Entry Guidance (FNPEG) and (ii) novel angle of attack and sideslip (α - β) guidance developed using FNPEG with Uncoupled Range Control [4]. Using Linear Quadratic Regulator (LQR) optimal control methods with state-feedback integral control designs, these guidance trajectories were designed to be tracked using (i) a conventional propulsive entry vehicle control hardware architecture - reaction control systems (RCS) and (ii) novel non-propulsive entry vehicle control systems - aerodynamic flap control system (FCS) and moving mass control system (MMCS) [1]. The novel FCS and MMCS architectures were designed to track α - β guidance commands while the RCS was designed to track bank angle commands. It was discovered that the asymmetric DEV, the PBV, experienced a non-zero induced roll moment due to sideslip that the FCS and MMCS architectures had limited capability to trim out. These two architectures were designed to provide independent angle of attack and sideslip commands with limited consideration for roll moment generation to trim. As a result, for the PBV, the FCS and MMCS configurations as designed, were limited in providing the control authority needed to track an α - β guidance trajectory [1]. These results form the motivation for the work presented in this paper - utilizing an aerodynamic control system to track α - β guidance commands for a symmetric DEV with the expectation that a symmetric entry vehicle will have zero or significantly reduced roll moment due to sideslip that the FCS can handle when tracking an α - β guidance trajectory. To demonstrate the feasibility of a novel guidance and control architecture on a DEV, we utilize a symmetric DEV, the PBV-II, for (i) the novel α - β guidance development using FNPEG with Uncoupled Range Control and (ii) LQR control design using eight aerodynamic control surfaces. This paper demonstrates that the novel uncoupled α - β guidance tracking can be achieved using aerodynamic control surfaces on a symmetric deployable entry vehicle configuration.

Wendy A Okolo↗

QBTNs - Quantum Boolean Tensor Networks

We develop algorithms and software that uses the D-Wave 2000Q quantum annealer to solve several types of Boolean tensor factorization problems. Boolean tensor factorization refers to the problem of representing a high-dimensional tensor filled with Boolean values as a product of smaller Boolean core tensors and Boolean matrices. We consider different tensor factorization models, including Boolean Tensor Train, Boolean Tucker, and Boolean Hierarchical Tucker. As an exact decomposition of a given type may not exist in the general case, the objective is to minimize the difference between the input high-dimensional tensor and the product of the lower-dimensional tensors of the proposed factorization, using a specified tensor norm. In our approach, we reduce the Boolean tensor factorization problem to a sequence of quadratic unconstrained binary optimization problems suitable for the D-Wave 2000Q quantum annealer. Although current quantum technology is still fairly restricted in the problems it can tackle, we show that complex tensor factorization problems as the ones addressed by us can be solved efficiently and accurately.

Alexandrov, Boian↗

Formation Tetrahedron Design for Phase I of the Magnetospheric Multiscale Mission

The Magnetospheric Multiscale Mission (MMS) is a NASA mission intended to make fundamental advancements in our understanding of the Earth's Magnetosphere. There are three processes that MMS will study including magnetic reconnection, charged particle acceleration, and turbulence. There are four phases in the nominal mission and this work addresses some of the outstanding issues in phase I. The nominal phase I orbit is 1.2 x 12 R(sub e) highly elliptic orbit with four spacecraft nominally forming a regular tetrahedron. In this paper we investigate the relative dynamics of the four MMS spacecraft about an assumed reference orbit. There are several tetrahedron dimensions required in Phase I of the mission and in this work we design optimal tetrahedrons for the 10 km baseline. The performance metric used in the optimization process is directly related to the science return, and is based on an extension of previous work performed by Glassmeier. The optimizer we use is a commercially available Sequential Quadratic Programming (SQP) routine. Multiple optimal solutions are found, and we characterize how the performance of the formation varies between different regions of the reference orbit.

Hughes, Steven P.↗

Least H 2 norm updating of quadratic interpolation models for derivative-free trust-region algorithms

One particular class of derivative-free optimization algorithms is trust-region algorithms based on quadratic models given by the under-determined interpolation. Different techniques in updating the quadratic model from iteration to iteration will give different interpolation models. We propose a new way to update the quadratic model by minimizing the $H^{2}$ norm of the difference between neighboring quadratic models. The motivation for applying the $H^{2}$ norm is given. The theoretical properties of our new updating technique are also presented. We propose the projection in the sense of $H^{2}$ norm and the interpolation error analysis of our model function. We obtain the coefficients of the quadratic model function using the Karush–Kuhn–Tucker (KKT) conditions. Numerical results show the advantages of our model on the test set considered, and the derivative-free algorithms based on our least $H^{2}$ norm updating quadratic model functions can solve test problems with fewer function evaluations than the algorithm based on the least Frobenius norm updating model and the other compared methods.

derivative-free optimization↗

Results of an integrated structure-control law design sensitivity analysis

Next generation air and space vehicle designs are driven by increased performance requirements, demanding a high level of design integration between traditionally separate design disciplines. Interdisciplinary analysis capabilities have been developed, for aeroservoelastic aircraft and large flexible spacecraft control for instance, but the requisite integrated design methods are only beginning to be developed. One integrated design method which has received attention is based on hierarchal problem decompositions, optimization, and design sensitivity analyses. This paper highlights a design sensitivity analysis method for Linear Quadratic Cost, Gaussian (LQG) optimal control laws, which predicts change in the optimal control law due to changes in fixed problem parameters using analytical sensitivity equations. Numerical results of a design sensitivity analysis for a realistic aeroservoelastic aircraft example are presented. In this example, the sensitivity of the optimally controlled aircraft's response to various problem formulation and physical aircraft parameters is determined. These results are used to predict the aircraft's new optimally controlled response if the parameter was to have some other nominal value during the control law design process. The sensitivity results are validated by recomputing the optimal control law for discrete variations in parameters, computing the new actual aircraft response, and comparing with the predicted response. These results show an improvement in sensitivity accuracy for integrated design purposes over methods which do not include changess in the optimal control law. Use of the analytical LQG sensitivity expressions is also shown to be more efficient that finite difference methods for the computation of the equivalent sensitivity information.

Gilbert, Michael G.↗

Robust perturbation guidance for the advanced launch system

Real-time optimal trajectory generation for constructing guidance laws seems feasible. However, these deterministic guidance schemes may be sensitive to system and environmental uncertainties. To overcome this deficiency, it is suggested that the accessory problem about an optimal path, which generates a quadratic cost criterion in the perturbed states and control away from the optimal trajectory, be used as the cost criterion. A modified quadratic cost is minimized with respect to the perturbed control but maximized with respect to additive uncertainties in the linearized dynamics and in the linearized measurements. In this way, a robust neighboring optimum scheme is constructed. A perturbed guidance scheme is stated. To illustrate the theory, the second stage of a launch vehicle is considered. An optimal path is generated, and about this path a robust neighboring optimum guidance law is constructed and tested.

Speyer, Jason L.↗

Advanced Energy Storage Management in Distribution Network

With increasing penetration of distributed generation (DG) in the distribution networks (DN), the secure and optimal operation of DN has become an important concern. In this paper, an iterative mixed integer quadratic constrained quadratic programming model to optimize the operation of a three phase unbalanced distribution system with high penetration of Photovoltaic (PV) panels, DG and energy storage (ES) is developed. The proposed model minimizes not only the operating cost, including fuel cost and purchasing cost, but also voltage deviations and power loss. The optimization model is based on the linearized sensitivity coefficients between state variables (e.g., node voltages) and control variables (e.g., real and reactive power injections of DG and ES). To avoid slow convergence when close to the optimum, a golden search method is introduced to control the step size and accelerate the convergence. The proposed algorithm is demonstrated on modified IEEE 13 nodes test feeders with multiple PV panels, DG and ES. Numerical simulation results validate the proposed algorithm. Various scenarios of system configuration are studied and some critical findings are concluded.

linearization↗

Minimum fuel horizontal flightpaths in the terminal area

The problem of minimum fuel airplane trajectories from arbitrary initial states to be fixed final state is considered. There are four state variables (two position coordinates, heading, and constrained velocity) and two constrained controls (thrust and bank angle). The fuel optimality of circular and straight line flightpaths is examined. Representative extremals (trajectories satisfying the necessary conditions of the minimum principle) of various types are computed and used to evaluate trajectories generated by an on line algorithm. Attention is paid to the existence of Darboux points (beyond which an extremal ceases to be globally optimal). One fuel flow rate model includes a term quadratic in thrust; hence, the optimal thrust is continuous and nonsingular. The other fuel flow rate model is linear in thrust, and consequently the optimal thrust is discontinuous and singular.

Kreindler, E.↗

Trajectory optimization for real-time guidance. I - Time-varying LQR on a parallel processor

A key algorithmic element of a real-time trajectory optimization hardware/software implementation, the quadratic program (QP) solver element, is presented. The purpose of the effort is to make nonlinear trajectory optimization fast enough to provide real-time commands during guidance of a vehicle such as an aeromaneuvering orbiter. Many methods of nonlinear programming require the solution of a QP at each iteration. In the trajectory optimization case the QP has a special dynamic programming structure, a LQR-like structure. QP algorithm speed is increased by taking advantage of this special structure and by parallel implementation.

Psiaki, Mark L.↗