Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

QUBO formulations for training machine learning models

Abstract Training machine learning models on classical computers is usually a time and compute intensive process. With Moore’s law nearing its inevitable end and an ever-increasing demand for large-scale data analysis using machine learning, we must leverage non-conventional computing paradigms like quantum computing to train machine learning models efficiently. Adiabatic quantum computers can approximately solve NP-hard problems, such as the quadratic unconstrained binary optimization (QUBO), faster than classical computers. Since many machine learning problems are also NP-hard, we believe adiabatic quantum computers might be instrumental in training machine learning models efficiently in the post Moore’s law era. In order to solve problems on adiabatic quantum computers, they must be formulated as QUBO problems, which is very challenging. In this paper, we formulate the training problems of three machine learning models—linear regression, support vector machine (SVM) and balanced k-means clustering—as QUBO problems, making them conducive to be trained on adiabatic quantum computers. We also analyze the computational complexities of our formulations and compare them to corresponding state-of-the-art classical approaches. We show that the time and space complexities of our formulations are better (in case of SVM and balanced k-means clustering) or equivalent (in case of linear regression) to their classical counterparts.

97 MATHEMATICS AND COMPUTING↗

Quantum annealing algorithms for Boolean tensor networks

Abstract Quantum annealers manufactured by D-Wave Systems, Inc., are computational devices capable of finding high-quality heuristic solutions of NP-hard problems. In this contribution, we explore the potential and effectiveness of such quantum annealers for computing Boolean tensor networks. Tensors offer a natural way to model high-dimensional data commonplace in many scientific fields, and representing a binary tensor as a Boolean tensor network is the task of expressing a tensor containing categorical (i.e., $$\{0, 1\}$$ { 0 , 1 } ) values as a product of low dimensional binary tensors. A Boolean tensor network is computed by Boolean tensor decomposition, and it is usually not exact. The aim of such decomposition is to minimize the given distance measure between the high-dimensional input tensor and the product of lower-dimensional (usually three-dimensional) tensors and matrices representing the tensor network. In this paper, we introduce and analyze three general algorithms for Boolean tensor networks: Tucker, Tensor Train, and Hierarchical Tucker networks. The computation of a Boolean tensor network is reduced to a sequence of Boolean matrix factorizations, which we show can be expressed as a quadratic unconstrained binary optimization problem suitable for solving on a quantum annealer. By using a novel method we introduce called parallel quantum annealing, we demonstrate that Boolean tensor’s with up to millions of elements can be decomposed efficiently using a DWave 2000Q quantum annealer.

97 MATHEMATICS AND COMPUTING↗

Solving the homogeneous Bethe-Salpeter equation with a quantum annealer

The homogeneous Bethe-Salpeter equation (hBSE), describing a bound system in a genuinely relativistic quantum-field theory framework, was solved for the first time by using a D-Wave quantum annealer. After applying standard techniques of discretization, the hBSE, in ladder approximation, can be formally transformed in a generalized eigenvalue problem (GEVP), with two square matrices: one symmetric and the other nonsymmetric. The latter matrix poses the challenge of obtaining a suitable formal approach for investigating the GEVP by means of a quantum annealer, i.e., to recast it as a quadratic unconstrained binary optimization problem. A broad numerical analysis of the proposed algorithms, applied to matrices of dimension up to 64, was carried out by using both the simulated-annealing package and the D-Wave . The numerical results very nicely compare with those obtained with standard classical algorithms, and also show interesting scalability features. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A QUBO formulation for top-τ eigencentrality nodes

The efficient calculation of the centrality or “hierarchy” of nodes in a network has gained great relevance in recent years due to the generation of large amounts of data. The eigenvector centrality (aka eigencentrality) is quickly becoming a good metric for centrality due to both its simplicity and fidelity. In this work we lay the foundations for solving the eigencentrality problem of ranking the importance of the nodes of a network with scores from the eigenvector of the network, using quantum computational paradigms such as quantum annealing and gate-based quantum computing. The problem is reformulated as a quadratic unconstrained binary optimization (QUBO) that can be solved on both quantum architectures. The results focus on correctly identifying a given number of the most important nodes in numerous networks given by the sparse vector solution of our QUBO formulation of the problem of identifying the top- τ highest eigencentrality nodes in a network on both the D-Wave and IBM quantum computers.

97 MATHEMATICS AND COMPUTING↗

Using Machine Learning for Quantum Annealing Accuracy Prediction

Quantum annealers, such as the device built by D-Wave Systems, Inc., offer a way to compute solutions of NP-hard problems that can be expressed in Ising or quadratic unconstrained binary optimization (QUBO) form. Although such solutions are typically of very high quality, problem instances are usually not solved to optimality due to imperfections of the current generations quantum annealers. In this contribution, we aim to understand some of the factors contributing to the hardness of a problem instance, and to use machine learning models to predict the accuracy of the D-Wave 2000Q annealer for solving specific problems. We focus on the maximum clique problem, a classic NP-hard problem with important applications in network analysis, bioinformatics, and computational chemistry. By training a machine learning classification model on basic problem characteristics such as the number of edges in the graph, or annealing parameters, such as the D-Wave’s chain strength, we are able to rank certain features in the order of their contribution to the solution hardness, and present a simple decision tree which allows to predict whether a problem will be solvable to optimality with the D-Wave 2000Q. We extend these results by training a machine learning regression model that predicts the clique size found by D-Wave.

97 MATHEMATICS AND COMPUTING↗

Approximation in LQG control of a thermoelastic rod

Control and estimator gains are computed for linear-quadratic-Gaussian (LQG) optimal control of the axial vibrations of a thermoelastic rod. The computations are based on a modal approximation of the partial differential equations representing the rod, and convergence of the approximations to control and estimator gains is the main issue.

Gibson, J. S.↗

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.↗