A validated integration algorithm for nonlinear ODEs using Taylor models and ellipsoidal calculus

This paper presents a novel algorithm for bounding the reachable set of parametric nonlinear differential equations. This algorithm is based on a first-discretize-then-bound approach to enclose the reachable set via propagation of a Taylor model with ellipsoidal remainder, and it accounts for truncation errors that are inherent to the discretization. In contrast to existing algorithms that proceed in two phases-an a priori enclosure phase, followed by a tightening phase-the proposed algorithm first predicts a continuous-time enclosure and then seeks a maximal step-size for which validity of the predicted enclosure can be established. It is shown that this reversed approach leads to a natural step-size control mechanism, which no longer relies on the availability of an a priori enclosure. Also described in the paper is an open-source implementation of the algorithm in ACADO Toolkit. A simple numerical case study is presented to illustrate the performance and stability of the algorithm.

[1]  Benoît Chachuat,et al.  Discretize-then-relax approach for convex/concave relaxations of the solutions of parametric ODEs , 2011 .

[2]  Moritz Diehl,et al.  ACADO toolkit—An open‐source framework for automatic control and dynamic optimization , 2011 .

[3]  Nedialko S. Nedialkov,et al.  Validated solutions of initial value problems for ordinary differential equations , 1999, Appl. Math. Comput..

[4]  Benoît Chachuat,et al.  Convergence analysis of Taylor models and McCormick-Taylor models , 2013, J. Glob. Optim..

[5]  W. Walter.  Differential and Integral Inequalities , 1970 .

[6]  Aaas News,et al.  Book Reviews , 1893, Buffalo Medical and Surgical Journal.

[7]  J. Lasserre.  Moments, Positive Polynomials And Their Applications , 2009 .

[8]  M. Stadtherr,et al.  Validated solutions of initial value problems for parametric ODEs , 2007 .

[9]  Claire S. Adjiman,et al.  A Rigorous Global Optimization Algorithm for Problems with Ordinary Differential Equations , 2002, J. Glob. Optim..

[10]  Boris Houska,et al.  Robust Optimization of Dynamic Systems (Robuuste optimalisatie van dynamische systemen) , 2011 .

[11]  Franco Blanchini,et al.  Set-theoretic methods in control , 2007 .

[12]  Garth P. McCormick,et al.  Computability of global solutions to factorable nonconvex programs: Part I — Convex underestimating problems , 1976, Math. Program..

[13]  Ramon E. Moore.  Methods and applications of interval analysis , 1979, SIAM studies in applied mathematics.

[14]  Paul I. Barton,et al.  Improved relaxations for the parametric solutions of ODEs using differential inequalities , 2012, Journal of Global Optimization.

[15]  Nedialko S. Nedialkov,et al.  On Taylor Model Based Integration of ODEs , 2007, SIAM J. Numer. Anal..

[16]  Martin Berz,et al.  Performance of Taylor Model Methods for Validated Integration of ODEs , 2004, PARA.

[17]  Benoît Chachuat,et al.  Bounding the Solutions of Parametric ODEs: When Taylor Models Meet Differential Inequalities , 2012 .

[18]  Benoît Chachuat,et al.  Branch-and-Lift Algorithm for Deterministic Global Optimization in Nonlinear Optimal Control , 2014, J. Optim. Theory Appl..

[19]  Jan Van Impe,et al.  Robust Optimization of Nonlinear Dynamic Systems with Application to a Jacketed Tubular Reactor , 2012 .

[20]  PAUL I. BARTON,et al.  Bounding the Solutions of Parameter Dependent Nonlinear Ordinary Differential Equations , 2005, SIAM J. Sci. Comput..

[21]  A. M. Sahlodin,et al.  Convex/concave relaxations of parametric ODEs using Taylor models , 2011, Comput. Chem. Eng..

[22]  Paul I. Barton,et al.  Nonlinear convex and concave relaxations for the solutions of parametric ODEs , 2013 .

[23]  Paul I. Barton,et al.  McCormick-Based Relaxations of Algorithms , 2009, SIAM J. Optim..

[24]  P. I. Barton,et al.  Global methods for dynamic optimization and mixed-integer dynamic optimization , 2006 .