Skip to main content

Lund University Publications

LUND UNIVERSITY LIBRARIES

Numerical and Symbolic Methods for Dynamic Optimization

Magnusson, Fredrik LU (2016) In PhD Thesis TFRT-1115
Abstract
Mathematical optimization is becoming increasingly important for engineering in general and control in particular. This thesis deals with numerical methods, primarily direct collocation, and symbolic methods, primarily block-triangular ordering and tearing, for numerical solution of general dynamic optimization problems involving dynamical systems modeled by large-scale differential-algebraic equations (DAE). These methods have been implemented in a software framework in the open-source JModelica.org platform, which is a software tool for simulation- and optimization-based analysis of DAEs described in the modeling language Modelica. The framework relies heavily upon the open-source, third-party software packages CasADi for symbolic... (More)
Mathematical optimization is becoming increasingly important for engineering in general and control in particular. This thesis deals with numerical methods, primarily direct collocation, and symbolic methods, primarily block-triangular ordering and tearing, for numerical solution of general dynamic optimization problems involving dynamical systems modeled by large-scale differential-algebraic equations (DAE). These methods have been implemented in a software framework in the open-source JModelica.org platform, which is a software tool for simulation- and optimization-based analysis of DAEs described in the modeling language Modelica. The framework relies heavily upon the open-source, third-party software packages CasADi for symbolic operations and algorithmic differentiation and IPOPT for solving the resulting nonconvex optimization problems.

Modelica is a standardized modeling language, which permeates the thesis. One of the many benefits of Modelica is that it is supported by several different tools, allowing implemented models to be used for different purposes. However, Modelica models are often developed for dynamic simulation and sometimes with little regard for numerics, which is enabled by the power of the available simulation software. Consequently, the models may be difficult to reuse for dynamic optimization, which is one of the challenges addressed by this thesis.

The application of direct collocation to DAE-constrained optimization problems is conventionally done by discretizing the full DAE. This often turns out to be inefficient, especially for DAEs originating from Modelica code. The thesis proposes various schemes to symbolically eliminate many of the algebraic variables in a preprocessing step before discretization to improve the efficiency of numerical methods for dynamic optimization, in particular direct collocation. These techniques are inspired by the causalization and tearing techniques often used when solving DAE initial-value problems in the Modelica community. Since sparsity is crucial for some dynamic optimization methods, we also propose a novel approach to preserving sparsity during this procedure.

A collection of five computationally challenging and industrially relevant optimal control problems is presented. The collection is used to evaluate the performance of the methods. We consider both computational time and probability of solving problems in a timely manner. We demonstrate that the proposed methods often are an order of magnitude faster than the standard way of discretizing the full DAE, and that they also increase probability of successful convergence significantly. It is also demonstrated that the methods are beneficial not only for DAEs originating from Modelica code, but also for more conventional textbook DAEs that have been developed specifically for optimization purposes. (Less)
Please use this url to cite or link to this publication:
author
supervisor
opponent
  • Assistant Professor Hedengren, John, Brigham Young University, Utah, USA
organization
publishing date
type
Thesis
publication status
published
subject
keywords
dynamic optimization, differential-algebraic equations, Modelica, direct collocation, block-triangular ordering
in
PhD Thesis TFRT-1115
pages
133 pages
publisher
Department of Automatic Control, Lund Institute of Technology, Lund University
defense location
M-huset, M:B
defense date
2016-11-18 10:15:00
ISSN
0280-5316
ISBN
978-91-7753-004-6
978-91-7753-005-3
project
LCCC
Numerical and Symbolic Algorithms for Dynamic Optimization
language
English
LU publication?
yes
id
4b50a0a7-844c-499b-9506-cfe9eee14e4b
alternative location
https://www.control.lth.se/media/2016/dynopt_thesis_web.pdf
date added to LUP
2016-10-12 10:35:07
date last changed
2018-11-21 21:26:33
@phdthesis{4b50a0a7-844c-499b-9506-cfe9eee14e4b,
  abstract     = {{Mathematical optimization is becoming increasingly important for engineering in general and control in particular. This thesis deals with numerical methods, primarily direct collocation, and symbolic methods, primarily block-triangular ordering and tearing, for numerical solution of general dynamic optimization problems involving dynamical systems modeled by large-scale differential-algebraic equations (DAE). These methods have been implemented in a software framework in the open-source JModelica.org platform, which is a software tool for simulation- and optimization-based analysis of DAEs described in the modeling language Modelica. The framework relies heavily upon the open-source, third-party software packages CasADi for symbolic operations and algorithmic differentiation and IPOPT for solving the resulting nonconvex optimization problems.<br/><br/>Modelica is a standardized modeling language, which permeates the thesis. One of the many benefits of Modelica is that it is supported by several different tools, allowing implemented models to be used for different purposes. However, Modelica models are often developed for dynamic simulation and sometimes with little regard for numerics, which is enabled by the power of the available simulation software. Consequently, the models may be difficult to reuse for dynamic optimization, which is one of the challenges addressed by this thesis.<br/><br/>The application of direct collocation to DAE-constrained optimization problems is conventionally done by discretizing the full DAE. This often turns out to be inefficient, especially for DAEs originating from Modelica code. The thesis proposes various schemes to symbolically eliminate many of the algebraic variables in a preprocessing step before discretization to improve the efficiency of numerical methods for dynamic optimization, in particular direct collocation. These techniques are inspired by the causalization and tearing techniques often used when solving DAE initial-value problems in the Modelica community. Since sparsity is crucial for some dynamic optimization methods, we also propose a novel approach to preserving sparsity during this procedure.<br/><br/>A collection of five computationally challenging and industrially relevant optimal control problems is presented. The collection is used to evaluate the performance of the methods. We consider both computational time and probability of solving problems in a timely manner. We demonstrate that the proposed methods often are an order of magnitude faster than the standard way of discretizing the full DAE, and that they also increase probability of successful convergence significantly. It is also demonstrated that the methods are beneficial not only for DAEs originating from Modelica code, but also for more conventional textbook DAEs that have been developed specifically for optimization purposes.}},
  author       = {{Magnusson, Fredrik}},
  isbn         = {{978-91-7753-004-6}},
  issn         = {{0280-5316}},
  keywords     = {{dynamic optimization; differential-algebraic equations; Modelica; direct collocation; block-triangular ordering}},
  language     = {{eng}},
  month        = {{11}},
  publisher    = {{Department of Automatic Control, Lund Institute of Technology, Lund University}},
  school       = {{Lund University}},
  series       = {{PhD Thesis TFRT-1115}},
  title        = {{Numerical and Symbolic Methods for Dynamic Optimization}},
  url          = {{https://lup.lub.lu.se/search/files/15978652/dynopt_thesis_print.pdf}},
  year         = {{2016}},
}