Skip to main content
Log in

Optimization Problem Coupled with Differential Equations: A Numerical Algorithm Mixing an Interior-Point Method and Event Detection

  • Published:
Journal of Optimization Theory and Applications Aims and scope Submit manuscript

    We’re sorry, something doesn't seem to be working properly.

    Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.

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.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Subscribe and save

Springer+ Basic
$34.99 /Month
  • Get 10 units per month
  • Download Article/Chapter or eBook
  • 1 Unit = 1 Article or 1 Chapter
  • Cancel anytime
Subscribe now

Buy Now

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Similar content being viewed by others

References

  1. 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)

    Article  MathSciNet  MATH  Google Scholar 

  2. 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)

    MathSciNet  MATH  Google Scholar 

  3. 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

  4. Rabier, P.J., Griewank, A.: Generic aspects of convexification with applications to thermodynamic equilibrium. Arch. Ration. Mech. Anal. 118(4), 349–397 (1992)

    Article  MathSciNet  MATH  Google Scholar 

  5. 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

  6. Park, T., Barton, P.I.: State event location in differential-algebraic models. ACM Trans. Model. Comput. Simul. (TOMACS) 6, 137–165 (1996)

    Article  MATH  Google Scholar 

  7. 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)

    Article  MathSciNet  MATH  Google Scholar 

  8. Singer, A.B., Barton, P.I.: Global optimization with nonlinear ordinary differential equations. J. Glob. Optim. 34, 159–190 (2006)

    Article  MathSciNet  MATH  Google Scholar 

  9. 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)

    Article  Google Scholar 

  10. Rockafellar, R.T.: Convex Analysis. Princeton University Press, Princeton (1970)

    MATH  Google Scholar 

  11. 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)

    Article  Google Scholar 

  12. 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)

    Article  MathSciNet  MATH  Google Scholar 

  13. Benson, H.Y., Shanno, D.F.: Interior-point methods for nonconvex nonlinear programming: Regularization and warmstarts. Comput. Optim. Appl. 40(2), 143–189 (2008)

    Article  MathSciNet  MATH  Google Scholar 

  14. Fiacco, A.V., McCormick, G.P.: Nonlinear Programming: Sequential Unconstrained Minimization Techniques. Wiley, New York (1968)

    MATH  Google Scholar 

  15. Yildirim, E.A., Wright, S.J.: Warm-start strategies in interior-point methods for linear programming. SIAM J. Optim. 12(3), 782–810 (2002)

    Article  MathSciNet  MATH  Google Scholar 

  16. Gear, C.W., Østerby, O.: Solving ordinary differential equations with discontinuities. ACM Trans. Math. Softw. 10(1), 23–44 (1984)

    Article  MATH  Google Scholar 

  17. Guglielmi, N., Hairer, E.: Computing breaking points in implicit delay differential equations. Adv. Comput. Math. 29(3), 229–247 (2008)

    Article  MathSciNet  MATH  Google Scholar 

  18. 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)

    MATH  Google Scholar 

  19. 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)

    Article  MathSciNet  Google Scholar 

  20. Deuflhard, P.: Newton Methods for Nonlinear Problems: Affine Invariance and Adaptive Algorithms. Springer-Verlag, Berlin (2004)

    MATH  Google Scholar 

  21. Mao, G., Petzold, L.R.: Efficient integration over discontinuities for differential-algebraic systems. Comput. Math. Appl. 43(1–2), 65–79 (2002)

    Article  MathSciNet  MATH  Google Scholar 

  22. 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)

    Google Scholar 

  23. 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)

    Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to A. Caboussat.

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

Reprints 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

Download citation

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10957-010-9714-1

Keywords