fast nonlinear model predictive control algorithms and
play

Fast Nonlinear Model Predictive Control Algorithms and Applications - PowerPoint PPT Presentation

Fast Nonlinear Model Predictive Control Algorithms and Applications in Process Engineering Moritz Diehl, Optimization in Engineering Center (OPTEC) & Electrical Engineering Department (ESAT) K.U. Leuven, Belgium INRIA-Rocquencourt, May


  1. Fast Nonlinear Model Predictive Control Algorithms and Applications in Process Engineering Moritz Diehl, Optimization in Engineering Center (OPTEC) & Electrical Engineering Department (ESAT) K.U. Leuven, Belgium INRIA-Rocquencourt, May 30, 2007

  2. Outline of the Talk � K.U.Leuven‘s Optimization in Engineering Center OPTEC � Nonlinear Model Predictive Control (NMPC) � How to solve dynamic optimization problems? � Four crucial features for fast NMPC algorithms � Application to a Distillation Column

  3. OPTEC Aim: Connect Optimization Methods & Applications Applications: Smart problem formulations allow efficient solution (e.g. convexity) Methods: New developments are inspired and driven by application needs

  4. Optimization in Engineering Center OPTEC Five year project, from 2005 to 2010, 500.000 Euro per year, about 20 professors, 10 postdocs, and 60 PhD students involved in OPTEC research Promoted by four departments: � Electrical Engineering � Mechanical Engineering � Chemical Engineering � Computer Science Many real world applications at OPTEC...

  5. Quarterly Stevin Lecture: Everyone Invited! � Quarterly „Simon Stevin Lecture on Optimization in Engineering“: • Dec 6: Larry Biegler , CMU Pittsburgh • Apr 18: Stephen Boyd , Stanford • July 9: Steve Wright , Madison, Wisconsin • Oct 24: Manfred Morari , ETH Zurich • Dec X: David Mayne , Imperial, London Simon Stevin, 1548-1620) � „K.U. Leuven Seminar on Optimization in Engineering“ : • Jan. 31: Mario Milanese (Torino): MPC of semi-active damping • Feb. 8: Philippe Toint (Namur): large scale optimization methods Lecture and following • Feb. 22: Peter Kuehl (Heidelberg): Robust optimal feedback control Reception in Arenberg Castle, Leuven • March 1: Yurii Nesterov (UCL)/ Florian Jarre (Duesseldorf): new optimization algorithms

  6. Outline of the Talk � K.U.Leuven‘s Optimization in Engineering Center OPTEC � Nonlinear Model Predictive Control (NMPC) � How to solve dynamic optimization problems? � Four crucial features for fast NMPC algorithms � Online MPC of a combustion engine

  7. First Principle Dynamic System Models E.g. some equations modelling a distillation column (in Stuttgart) � Nonlinear differential algebraic equations (DAE) � often in modeling languages like gPROMS, SIMULINK, Modelica � typical order of magnitude: some hundreds to thousands variables � difficulties: stiffness, discontinuities, high index Can we use these models directly for optimization and feedback control?

  8. Principle of Optimal Feedback Control / Nonlinear MPC: x 0 u 0 x 0 u 0

  9. Nonlinear Model Predictive Control When We Drive a Car Always look a bit into the future! Brain predicts and optimizes: e.g. slow down before curve Main challenge for NMPC: fast and reliable real-time optimization!

  10. NMPC applications, with decreasing timescales Polymeri- Distillation Combined Cycle sation column (with Power Plant (with reactor (with Univ. Stuttgart) Univ. Pavia) BASF) Chromatographic PET plant: Plant wide Separation (with control project with Univ. Dortmund) Politecnico di Milano Robot arms (with Looping kites for power Columbia Univ. & generation, with TU Delft, INRIA Grenoble) Politecnico di Torino Car Engines: EU Project with Univ. Linz, Stuttgart, Politecnico di Milano

  11. Outline of the Talk � K.U.Leuven‘s Optimization in Engineering Center OPTEC � Nonlinear Model Predictive Control (NMPC) � How to solve dynamic optimization problems? � Four crucial features for fast NMPC algorithms � Application to a Distillation Column

  12. Optimal Control Family Tree

  13. Optimal Control Family Tree

  14. Optimal Control Family Tree

  15. Optimal Control Problem in Simplest Form

  16. Simplest Approach: Single Shooting

  17. Nonlinear Program (NLP) in Single Shooting � After control discretization, obtain NLP: � Solve with NLP solver, e.g. Sequential Quadratic Programming (SQP)

  18. Sequential Quadratic Programming (SQP)

  19. Toy Problem with One ODE for Illustration Mildly nonlinear and unstable system.

  20. Single Shooting

  21. Single Shooting: First Iteration

  22. Single Shooting: Second Iteration

  23. Single Shooting: Third Iteration

  24. Single Shooting: 4th Iteration

  25. Single Shooting: 5th Iteration

  26. Single Shooting: 6th Iteration

  27. Single Shooting: 7th Iteration (Solution)

  28. Single Shooting: Pros and Cons

  29. Alternative: Direct Multiple Shooting [Bock, Plitt 1981]

  30. Nonlinear Program in Multiple Shooting

  31. SQP for Multiple Shooting Summarize NLP:

  32. Toy Example: Multiple Shooting Initialization

  33. Multiple Shooting: First Iteration

  34. Multiple Shooting: Second Iteration

  35. Multiple Shooting: 3 rd Iteration (already solution!)

  36. Multiple Shooting: 3 rd Iteration (already solution!) Single shooting converged much slower!

  37. The MUSCOD-II Developer Team [Heidelberg, Leuven, Madrid]

  38. Outline of the Talk � K.U.Leuven‘s Optimization in Engineering Center OPTEC � Nonlinear Model Predictive Control (NMPC) � How to solve dynamic optimization problems? � Four crucial features for fast NMPC algorithms � Online MPC of a combustion engine

  39. NMPC Computation from 1998 to 2006 � 1998: 5th order distillation model allows sampling times of only 5 minutes [Allgower, Findeisen, 1998] � 2001: 206th order distillation model, sampling times of 20 seconds [D. et al. ‚01] � 2006: 5th order engine model, sampling times of 10-20 milliseconds [Ferreau et al. ‘06] 5*60*1000 / 20 = 15 000 times faster, due to Moore‘s law + Algorithm Development

  40. NMPC Computation from 1998 to 2006 � 1998: 5th order distillation model allows sampling times of only 5 minutes [Allgower, Findeisen, 1998] � 2001: 206th order distillation model, sampling times of 20 seconds [D. et al. ‚01] cf. [Biegler] � 2006: 5th order engine model, sampling times of 10-20 milliseconds [Ferreau et al. ‘06], [Albersmeyer, Findeisen `06] 5*60*1000 / 20 = 15 000 times faster, due to Moore‘s law + Algorithm Development

  41. Four Crucial Features for Fast NMPC � Direct, simultaneous optimal control: Multiple Shooting � Efficient derivative generation for ODE/DAE solvers � Initialization by „ Initial Value Embedding “ � Real-Time Iterations for fast tracking of optimal solutions

  42. x 0

  43. x 0

  44. x 0

  45. x 0

  46. x 0

  47. x 0

  48. x 0

  49. Never simulate a nonlinear system open-loop! Conventional: : Conventional Initial Initial Value Value Embedding: Embedding:

  50. x 0 x 0

  51. x 0 x 0

  52. x 0 x 0

  53. x 0 u 0

  54. Real-Time Iterations minimize feedback delay

  55. x 0 u 0

  56. Outline of the Talk � K.U.Leuven‘s Optimization in Engineering Center OPTEC � Nonlinear Model Predictive Control (NMPC) � How to solve dynamic optimization problems? � Four crucial features for fast NMPC algorithms � Application to a Distillation Column

  57. Transient in 15 minutes instead of 2 hours!

  58. Conclusions � Recent progress makes Nonlinear MPC with first principles models in millisecond range possible (now 15 000 x faster than 1998) � Emerging consensus for NMPC algorithms: • employ direct, simultaneous methods • use Initial Value Embedding (first order predictor) • perform Real-Time Iterations to trace NMPC problem solution while data change • Use SQP type method to track active set changes

  59. An Invitation � 13th Czech-French-German Conference on Optimization (CFG07), Heidelberg, Germany, September 17-21, 2007. (inv. speakers: Fletcher, Scherer, Trelat, Waechter,...) Traditionally strong in optimal control.

  60. 4 PhD Positions in Numerical Optimization: • Sequential Convex Programming Algorithms for Nonlinear SDP • Large Scale & PDE Constrained Real-Time Optimization Algorithms • Fast Model Predictive Control Applications in Mechatronic Systems • Shape Optimization of Mechanical Parts under Inertia Loading (deadline: June 21, 2007)

Recommend


More recommend