Very Low Thrust Trajectory Optimization

We derive local error estimates for the discretization of optimal control problems governed by linear advection-diffusion partial differential equations (PDEs) using the streamline upwind/Petrov Galerkin (SUPG) stabilized finite element method. We show that if the SUPG method is used to solve optimization problems governed by an advection-dominated PDE, the convergence properties of the SUPG method is substantially different from the convergence properties of the SUPG method applied for the solution of an advection-dominated PDE. The reason is that the solution of the optimal control problem involves another advection-dominated PDE, the so-called adjoint equation, whose advection field is just the negative of the advection of the governing PDEs. For the solution of the optimal control problem, a coupled system involving both the original governing PDE as well as the adjoint PDE must be solved. We show that in the presence of a boundary layer, the local error between the solution of the SUPG discretized optimal control problem and the solution of the infinite dimensional problem is only of first order even if the error is computed locally in a region away from the boundary layer. In the presence of interior layers, we prove optimal convergence rates for the local error in a region away from the layer between the solution of the SUPG discretized optimal control problems and the solution of the infinite dimensional problem. Numerical examples are presented to illustrate some of the theoretical results.

This article presents a unified approach to simultaneous solution of optimization and optimal control problems in batch distillation, operating under different modes of operation: variable, constant, or optimal reflux. The simplified, computationally efficient short‐cut method and a novel algorithm to solve the optimal control problems in batch distillation is the basis of this unified approach. The short‐cut method identifies the feasible region of operation essential for optimization and optimal control problems, and provides analytical partial derivatives of the model parameters crucial to the solution.The new algorithm for the solution of optimal combination of the maximum principle and NLP optimization techniques. It circumvents the problems associated with the maximum principle approach (iterative solution of a two‐point boundary value problem, undounded control variables, and inability to handle the simultaneous optimization and optimal control problem), and the coupled ODE discretization‐NLP optimization scheme for nonlinear models (higher system nonlinearities, multiplicity of solutions, sensitivity of convergence to initial guesses). This algorithm reduces the dimensionality of the problem, and the nature of the algorithm allows a common platform to optimal solutions of different operating conditions. This article also shows that different categories of the optimal control problems in batch distillation essentially involve the solution of the maximum distillate problem.

Mathematical formulation and numerical solutions of an optimal Dirichlet boundary control problem for the Boussinesq equations are considered. The solution of the optimal control problem is obtained by adjusting of the temperature on the boundary. We analyze finite element approximations. A gradient method for the solution of the discrete optimal control problem is presented and analyzed. Finally, the results of some computational experiments are presented.

The most effective numerical techniques for the solution of trajectory optimization and optimal control problems combine a nonlinear iteration procedure with some type of parametric approximation to the trajectory dynamics. Early methods attempted to parameterize the dynamics using a small number of variables because the iterative search procedures could not successfully solve larger problems. With the development of more robust nonlinear programming algorithms, it is now feasible and desirable to consider formulations of the trajectory optimization problem incorporating a large number of variables and constraints. The purpose of this paper is to address the manner in which a trajectory is parameterized and the design of the nonlinear programming algorithm to effectively deal with this formulation.

1 A Survey on Computational Optimal Control.- Issues in the Direct Transcription of Optimal Control Problems to Sparse Nonlinear Programs.- Optimization in Control of Robots.- Large-scale SQP Methods and their Application in Trajectory Optimization.- Solving Optimal Control and Pursuit-Evasion Game Problems of High Complexity.- 2 Theoretical Aspects of Optimal Control and Nonlinear Programming.- Continuation Methods In Boundary Value Problems.- Second Order Optimality Conditions for Singular Extremals.- Synthesis of Adaptive Optimal Controls for Linear Dynamic Systems.- Control Applications of Reduced SQP Methods.- Time Optimal Control of Mechanical Systems.- 3 Algorithms for Optimal Control Calculations.- Second Order Algorithm for Time Optimal Control of a Linear System.- An SQP-type Solution Method for Constrained Discrete-Time Optimal Control Problems.- Numerical Methods for Solving Differential Games, Prospective Applications to Technical Problems.- Construction of the Optimal Feedback Controller for Constrained Optimal Control Problems with Unknown Disturbances.- Repetitive Optimization for Predictive Control of Dynamic Systems under Uncertainty.- Optimal Control of Multistage Systems Described by High-Index Differential-Algebraic Equations.- A New Class of a High Order Interior Point Method for the Solution of Convex Semiinfinite Optimization Problems.- A Structured Interior Point SQP Method for Nonlinear Optimal Control Problems.- 4 Software for Optimal Control Calculations.- Automated Approach for Optimizing Dynamic Systems.- ANDECS: A Computation Environment for Control Applications of Optimization.- Application of Automatic Differentiation to Optimal Control Problems.- OCCAL: A mixed symbolic-numeric Optimal Control CALculator.- 5 Applications of Optimal Control.- A Robotic Satellite with Simplified Design.- Nonlinear Control under Constraints of a Biological System.- An Object-Oriented Approach to Optimally Describe and Specify a SCADA System Applied to a Power Network.- Near-Optimal Flight Trajectories Generated by Neural Networks.- Performance of a Feedback Method with Respect to Changes in the Air-Density during the Ascent of a Two-Stage-To-Orbit Vehicle.- Linear Optimal Control for Reentry Flight.- Steady-State Modelling of Turbine Engine with Controllers.- Shortest Paths for Satellite Mounted Robot Manipulators.- Optimal Control of the Industrial Robot Manutec r3.

An enhanced scheme of transcribing the system dynamics for the numerical solution of optimal control problems is proposed. This new scheme is based on the standard method of direct collocation that converts an optimal control problem into a nonlinear programming problem via simultaneous state and control discretization. When compared with the standard method, the enhanced scheme has the advantage of higher solution accuracy with minimal additional computational effort. It is particularly suited for systems with states that are related to each other in a special form. For such systems, the ensuing nonlinear programming problem has the same number of constraints as those using the standard method. Numerical results on several optimal control problems using the enhanced scheme are presented, together with comparisons with the results obtained from the standard scheme.

Optimal control problems are common in aerospace engineering. A Python software program called PySCP is described for solving multiple-phase optimal control problems using sequential convex programming methods. By constructing a series of approximated second-order cone programming subproblems, PySCP approaches to the solution of the original optimal control problem in an iterative way. The key components of the software are described in detail, including convexification, discretization, and the adaptive trust region method. The convexification of the first-order differential dynamic equation is implemented using successive linearization. Six discretization methods, including zero-order hold, first-order hold, Runge-Kutta, and three hp pseudospectral collocation methods, are implemented so that different types of optimal control problems can be tackled efficiently. Adaptive trust region method is employed, and robust convergence is achieved. Both free-final-time problem and fixed-final-time problem can be solved by the software. The application of the software is demonstrated on three optimal control problems with varying complexity. PySCP provides researchers a useful toolkit to solve a wide variety of optimal control problems using sequential convex programming.

Existing methods for nonlinear robust control often use scenario‐based approaches to formulate the control problem as large nonlinear optimization problems. The optimization problems are challenging to solve due to their size, especially if the control problems include time‐varying uncertainty. This paper draws from local reduction methods used in semi‐infinite optimization to solve robust optimal control problems with parametric and time‐varying uncertainty. By iteratively adding interim worst‐case scenarios to the problem, methods based on local reduction provide a way to manage the total number of scenarios. We show that the local reduction method for optimal control problems consists of solving a series of simplified optimal control problems to find worst‐case constraint violations. In particular, we present examples where local reduction methods find worst‐case scenarios that are not on the boundary of the uncertainty set. We also provide bounds on the error if local solvers are used. The proposed approach is illustrated with two case studies with parametric and additive time‐varying uncertainty. In the first case study, the number of scenarios obtained from local reduction is 101, smaller than in the case when all extreme scenarios are considered. In the second case study, the number of scenarios obtained from the local reduction is two compared to 512 extreme scenarios. Our approach was able to satisfy the constraints both for parametric uncertainty and time‐varying disturbances, whereas approaches from literature either violated the constraints or became computationally expensive.

Three algorithms for efficient solution of optimal control problems for high-dimensional systems are presented. Each bypasses the intermediate (and, unnecessary) step of open-loop model reduction. Each also bypasses the solution of the full Riccati equation corresponding to the LQR problem, which is numerically intractable for large n. Motivation for this effort comes from the field of model-based flow control, where open-loop model reduction often fails to capture the dynamics of interest (governed by the Navier–Stokes equation). Our minimum control energy method is a simplified expression for the well-known minimum-energy stabilizing control feedback that depends only on the left eigenvectors corresponding to the unstable eigenvalues of the system matrix A. Our Adjoint of the Direct-Adjoint method is based on the repeated iterative computation of the adjoint of a forward problem, itself defined to be the direct-adjoint vector pair associated with the LQR problem. Our oppositely-shifted subspace iteration (OSSI, the main new result of the present paper) method is based on our new subspace iteration method for computing the Schur vectors corresponding, notably, to the \(m\ll n\) central eigenvalues (near the imaginary axis) of the Hamiltonian matrix related to the Riccati equation of interest. Prototype OSSI implementations are tested on a low-order control problem to illustrate its behavior.

We have already briefly touched nonlinear programming in the last two chapters, in which we introduced the multiple shooting discretized optimal control problem, a Nonlinear Program (NLP), and presented the outer convexification and relaxation approach that allows to compute approximations to local Mixed–Integer Optimal Control Problem (MIOCP) solutions by solving a reformulated and discretized but possibly much larger Optimal Control Problem (OCP). This chapter is concerned with theory and numerical methods for the solution of NLPs and equips the reader with definitions and algorithms required for the following chapters of this thesis. We present optimality conditions characterizing locally optimal solutions and introduce Sequential Quadratic Programming (SQP) methods for the iterative solution of NLPs. The evaluation of the matching condition constraints of the discretized optimal control problem requires the solution of Initial Value Problems (IVPs). To this end, we present one step methods for non–stiff Ordinary Differential Equations (ODEs) and discuss the efficient and consistent computation of sensitivities of IVPs solutions. The familiar reader may wish to continue with chapter 4 rightaway, and refer back to this chapter should the need arise.KeywordsSequential Quadratic ProgrammingActive ConstraintAutomatic DifferentiationDiscretization GridStrict ComplementarityThese keywords were added by machine and not by the authors. This process is experimental and the keywords may be updated as the learning algorithm improves.

A neural network based optimal control synthesis is presented for solving optimal control problems with control and state constraints. The optimal control problem is transcribed into nonlinear programming problem which is implemented with adaptive critic neural network. The proposed simulation methods is illustrated by the optimal control problem of feeding adaptation of filter feeders of Daphnia. Results show that adaptive critic based systematic approach holds promise for obtaining the optimal control with control and state constraints.

This paper describes a novel offline approach to estimate the parameters of a nonlinear system using the optimal control theory. In the conventional methods, only the output is considered in the identification process and the problem is stated as a nonlinear optimization problem. Nevertheless, the proposed method here is based on the least square error of both output and input. To do this, the problem is first stated as a tracking problem with some unknown parameters. The necessary and sufficiency conditions of optimality conditions are derived using the fundamental theorem of variational calculus. The problem is deployed as an effort optimal control problem which finally leads to a two point boundary value problem. To solve the resultant equations, a secondary optimal control problem is deployed to be solved using the homotopy algorithm. So the parameters are gradually converged to the real values. The proposed algorithm is analyzed for a one-link manipulator with the full measurements. Thereafter, the proposed method is applied on a real one-link manipulator and nonlinear model of permanent magnet dc motor. Application of the method for the identification of a few test beds at the AUT proves the validity and efficiency of the proposed approach.

This paper deals with the reconstruction of the desired demand in an optimal control problem, stated over a tree-shaped transportation network which is governed by a linear hyperbolic conservation law. As desired demands typically undergo fluctuations due to seasonality or unexpected events making short-term adjustments necessary, such an approach can exemplary be used for forecasting from past data. We suggest to model this problem as a so-called inverse optimal control problem, i.e., a hierarchical optimization problem whose inner problem is the optimal control problem and whose outer problem is the reconstruction problem. In order to guarantee the existence of solutions in the function space framework, the hyperbolic conservation law is interpreted in weak sense allowing for control functions in Lebesgue spaces. For the computational treatment of the model, we transfer the hierarchical problem into a nonsmooth single-level one by plugging the uniquely determined solution of the inner optimal control problem into the outer reconstruction problem before applying techniques from nonsmooth optimization. Some numerical experiments are presented to visualize various features of the model including different types of noise in the demand and strategies of how to observe the network in order to obtain good reconstructions of the desired demand.

This paper presents an efficient and a novel implementation of a combined multiple shooting and collocation (CMSC) algorithm for the solution of nonlinear optimal control problems. The implemented algorithm is a modification of the approach proposed in [17, 18]. The new implementation is done under the JModelica.org framework along with CasADi and Ipopt. The framework uses a symbolic pre-calculation of functions and derivatives. Besides the integration of various components of JModelica.org, Ipopt, and CasADi, the implementation facilitates simpler modeling of optimal control problems along with a choice of options for various linear algebra algorithms. The paper gives a description of the algorithm and elaborates the components of the framework. Numerical experimentations show that the new implementation is efficient in com

Numerical solutions of multi-phase optimal control problems are considered, where a phase is defined as a subinterval in which the right-hand sides of the differential equations are continuous. At the phase boundaries, the state, control, and design vectors may have jumps; in addition, the dimension of these vectors as well as the dimension of the right-hand sides may change. Two direct methods for solving such problems are presented : a multiple-shooting approach and a collocation method. Both methods avoid the adjoint differential equations. By these transcriptions, the optimal control problem is converted into a nonlinear programming problem (NLP), which is solved using standard sequential quadratic programming (SQP) algorithms. The methods are applied to determine optimal ascent and reentry trajectories of a two-stage-to-orbit vehicle. Both methods are embedded into an advanced user interface which allows a user to edit most of the necessary input in a simple way. This utility is described briefly.