Optimization-Constrained Differential Equations with Active Set Changes

Foundational theory is established for nonlinear differential equations with embedded nonlinear optimization problems exhibiting active set changes. Existence, uniqueness, and continuation of solutions are shown, followed by lexicographically smooth (implying Lipschitzian) parametric dependence. The sensitivity theory found here accurately characterizes sensitivity jumps resulting from active set changes via an auxiliary nonsmooth sensitivity system obtained by lexicographic directional differentiation. The results in this article hold under easily verifiable regularity conditions (linear independence of constraints and strong second-order sufficiency), which are shown to imply generalized differentiation index one of a nonsmooth differential-algebraic equation system obtained by replacing the optimization problem with its optimality conditions and recasting the complementarity conditions as nonsmooth algebraic equations. The theory in this article is computationally relevant, allowing for implementation of dynamic optimization strategies (i.e., open-loop optimal control), and recovers (and rigorously formalizes) classical results in the absence of active set changes. Along the way, contributions are made to the theory of piecewise differentiable functions.

[1]  P. I. Barton,et al.  Efficient sensitivity analysis of large-scale differential-algebraic systems , 1997 .

[2]  C. W. Gear,et al.  The index of general nonlinear DAEs , 1995 .

[3]  Yurii Nesterov,et al.  Lexicographic differentiation of nonsmooth functions , 2005, Math. Program..

[4]  Paul I. Barton,et al.  Generalized Derivatives for Hybrid Systems , 2017, IEEE Transactions on Automatic Control.

[5]  Andreas Griewank,et al.  Evaluating derivatives - principles and techniques of algorithmic differentiation, Second Edition , 2000, Frontiers in applied mathematics.

[6]  Paul I. Barton,et al.  A vector forward mode of automatic differentiation for generalized derivative evaluation , 2015, Optim. Methods Softw..

[7]  S. Scholtes Introduction to Piecewise Differentiable Equations , 2012 .

[8]  Paul I. Barton,et al.  Computationally relevant generalized derivatives: theory, evaluation and applications , 2018, Optim. Methods Softw..

[9]  P. I. Barton,et al.  Parametric sensitivity functions for hybrid discrete/continuous systems , 1999 .

[10]  Uwe Naumann,et al.  The Art of Differentiating Computer Programs - An Introduction to Algorithmic Differentiation , 2012, Software, environments, tools.

[11]  Paul I. Barton,et al.  Generalized Sensitivity Analysis of Nonlinear Programs , 2018, SIAM J. Optim..

[12]  Paul I. Barton,et al.  Generalized derivatives of optimal control problems with nonsmooth differential-algebraic equations embedded , 2016, 2016 IEEE 55th Conference on Decision and Control (CDC).

[13]  Jan Vlcek,et al.  A bundle-Newton method for nonsmooth unconstrained minimization , 1998, Math. Program..

[14]  Ernst Hairer,et al.  Solving Optimization-Constrained Differential Equations with Discontinuity Points, with Application to Atmospheric Chemistry , 2009, SIAM J. Sci. Comput..

[15]  Shengtai Li,et al.  Adjoint Sensitivity Analysis for Differential-Algebraic Equations: The Adjoint DAE System and Its Numerical Solution , 2002, SIAM J. Sci. Comput..

[16]  Uri M. Ascher,et al.  Computer methods for ordinary differential equations and differential-algebraic equations , 1998 .

[17]  W. Rheinboldt,et al.  Theoretical and numerical analysis of differential-algebraic equations , 2002 .

[18]  Paul I. Barton,et al.  Generalized Derivatives for Solutions of Parametric Ordinary Differential Equations with Non-differentiable Right-Hand Sides , 2014, J. Optim. Theory Appl..

[19]  Andreas Griewank,et al.  On stable piecewise linearization and generalized algorithmic differentiation , 2013, Optim. Methods Softw..

[20]  F. Facchinei,et al.  Finite-Dimensional Variational Inequalities and Complementarity Problems , 2003 .

[21]  Paul I. Barton,et al.  Nonsmooth DAEs with Applications in Modeling Phase Changes , 2018 .

[22]  Vladimir M. Veliov,et al.  On the Time-Discretization of Control Systems , 1997 .

[23]  Anthony V. Fiacco,et al.  Introduction to Sensitivity and Stability Analysis in Nonlinear Programming , 2012 .

[24]  Liqun Qi,et al.  A nonsmooth version of Newton's method , 1993, Math. Program..

[25]  F. Doyle,et al.  Dynamic flux balance analysis of diauxic growth in Escherichia coli. , 2002, Biophysical journal.

[26]  Paul I. Barton,et al.  Dependence of solutions of nonsmooth differential-algebraic equations on parameters , 2017 .

[27]  A dynamic optimization problem related to organic aerosols , 2007 .

[28]  Daniel Ralph,et al.  Sensitivity analysis of composite piecewise smooth equations , 1997, Math. Program..

[29]  Jong-Shi Pang,et al.  Differential variational inequalities , 2008, Math. Program..

[30]  J. Seinfeld,et al.  Primal-Dual Interior-Point Method for an Optimization Problem Related to the Modeling of Atmospheric Organic Aerosols , 2006 .

[31]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[32]  Kamil A. Khan Branch-locking AD techniques for nonsmooth composite functions and nonsmooth implicit functions , 2018, Optim. Methods Softw..

[33]  J. M. Schumacher,et al.  Complementarity systems in optimization , 2004, Math. Program..

[34]  F. Clarke Optimization And Nonsmooth Analysis , 1983 .

[35]  Paul I. Barton,et al.  Generalized Derivatives of Differential–Algebraic Equations , 2016, J. Optim. Theory Appl..

[36]  B. Palsson,et al.  Constraining the metabolic genotype–phenotype relationship using a phylogeny of in silico methods , 2012, Nature Reviews Microbiology.

[37]  Volker Mehrmann,et al.  Differential-Algebraic Equations: Analysis and Numerical Solution , 2006 .

[38]  Chantal Landry,et al.  Optimization Problem Coupled with Differential Equations: A Numerical Algorithm Mixing an Interior-Point Method and Event Detection , 2010, J. Optim. Theory Appl..