Abstract
The numerical analysis of a dynamic constrained optimization problem is presented. It consists of a global minimization problem that is coupled with a system of ordinary differential equations. The activation and the deactivation of inequality constraints induce discontinuity points in the time evolution. A numerical method based on an operator splitting scheme and a fixed point algorithm is advocated. The ordinary differential equations are approximated by the Crank-Nicolson scheme, while a primal-dual interior-point method with warm-starts is used to solve the minimization problem. The computation of the discontinuity points is based on geometric arguments, extrapolation polynomials and sensitivity analysis. Second order convergence of the method is proved when an inequality constraint is activated. Numerical results for atmospheric particles confirm the theoretical investigations.
Similar content being viewed by others
References
Amundson, N.R., Caboussat, A., He, J.W., Seinfeld, J.H.: Primal-dual interior-point algorithm for chemical equilibrium problems related to modeling of atmospheric organic aerosols. J. Optim. Theory Appl. 130(3), 375–407 (2006)
Amundson, N.R., Caboussat, A., He, J.W., Landry, C., Seinfeld, J.H.: A dynamic optimization problem related to organic aerosols. C. R. Acad. Sci. 344(8), 519–522 (2007)
Landry, C.: Numerical analysis of optimization-constrained differential equations: Applications to atmospheric chemistry. PhD thesis, Ecole Polytechnique Fédérale de Lausanne (2009). Available at http://library.epfl.ch/theses/?nr=4345. Accessed 18 May 2010
Rabier, P.J., Griewank, A.: Generic aspects of convexification with applications to thermodynamic equilibrium. Arch. Ration. Mech. Anal. 118(4), 349–397 (1992)
Caboussat, A.: Primal-dual interior-point method for thermodynamic gas-particle partitioning. Comput. Optim. Appl. (2009). http://dx.doi.org/10.1007/s10589-009-9262-5
Park, T., Barton, P.I.: State event location in differential-algebraic models. ACM Trans. Model. Comput. Simul. (TOMACS) 6, 137–165 (1996)
Shampine, L.F., Gladwell, I., Brankin, R.W.: Reliable solution of special event location problems for ODEs. ACM Trans. Math. Softw. 17, 11–25 (1991)
Singer, A.B., Barton, P.I.: Global optimization with nonlinear ordinary differential equations. J. Glob. Optim. 34, 159–190 (2006)
Esposito, J.M., Kumar, V.: A state event detection algorithm for numerically simulating hybrid systems with model singularities. ACM Trans. Model. Comput. Simul. (TOMACS) 17, 1–22 (2007)
Rockafellar, R.T.: Convex Analysis. Princeton University Press, Princeton (1970)
McDonald, C.M., Floudas, C.A.: GLOPEQ: A new computational tool for the phase and chemical equilibrium problem. Comput. Chem. Eng. 21(1), 1–23 (1996)
Gondzio, J., Grothey, A.: A new unblocking technique to warmstart interior point methods based on sensitivity analysis. SIAM J. Optim. 19(3), 1184–1210 (2008)
Benson, H.Y., Shanno, D.F.: Interior-point methods for nonconvex nonlinear programming: Regularization and warmstarts. Comput. Optim. Appl. 40(2), 143–189 (2008)
Fiacco, A.V., McCormick, G.P.: Nonlinear Programming: Sequential Unconstrained Minimization Techniques. Wiley, New York (1968)
Yildirim, E.A., Wright, S.J.: Warm-start strategies in interior-point methods for linear programming. SIAM J. Optim. 12(3), 782–810 (2002)
Gear, C.W., Østerby, O.: Solving ordinary differential equations with discontinuities. ACM Trans. Math. Softw. 10(1), 23–44 (1984)
Guglielmi, N., Hairer, E.: Computing breaking points in implicit delay differential equations. Adv. Comput. Math. 29(3), 229–247 (2008)
Hairer, E., Nørsett, S.P., Wanner, G.: Solving Ordinary Differential Equations I: Nonstiff Problems. 2nd edn. Springer Series in Computational Mathematics, vol. 8. Springer-Verlag, Berlin (1993)
Landry, C., Caboussat, A., Hairer, E.: Solving optimization-constrained differential equations with discontinuity points, with application to atmospheric chemistry. SIAM J. Sci. Comput. 31(5), 3806–3826 (2009)
Deuflhard, P.: Newton Methods for Nonlinear Problems: Affine Invariance and Adaptive Algorithms. Springer-Verlag, Berlin (2004)
Mao, G., Petzold, L.R.: Efficient integration over discontinuities for differential-algebraic systems. Comput. Math. Appl. 43(1–2), 65–79 (2002)
Caloz, G., Rappaz, J.: Numerical analysis for nonlinear and bifurcation problems. In: Ciarlet, P.G., Lions, J.L. (eds.) Handbook of Numerical Analysis, vol. V, pp. 487–637. Elsevier, Amsterdam (1997)
Rappaz, J.: Numerical approximation of PDEs and Clément’s interpolation. In: Operator Theory: Advances and Applications, vol. 168, pp. 237–250. Birkhäuser, Basel (2006)
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by R. Glowinski.
This work was partially supported by US Environmental Protection Grant X-83234201. The first author gratefully acknowledges the support of the Institute of Analysis and Scientific Computing (EPFL); part of this work has been achieved during his sabbatical leave in 2009–2010 at EPFL. The authors thank Prof. Ernst Hairer (University of Geneva) and Prof. Jiwen He (University of Houston) for fruitful discussions.
Rights and permissions
About this article
Cite this article
Caboussat, A., Landry, C. & Rappaz, J. Optimization Problem Coupled with Differential Equations: A Numerical Algorithm Mixing an Interior-Point Method and Event Detection. J Optim Theory Appl 147, 141–156 (2010). https://doi.org/10.1007/s10957-010-9714-1
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10957-010-9714-1