Deterministic Global Optimization Algorithm and Nonlinear Dynamics
|
|
- Suzan Norris
- 5 years ago
- Views:
Transcription
1 Deterministic Global Optimization Algorithm and Nonlinear Dynamics C. S. Adjiman and I. Papamichail Centre for Process Systems Engineering Department of Chemical Engineering and Chemical Technology Imperial College London Global Optimization Theory Institute, Argonne National Laboratory September 8-10, 2003 Financial support: Engineering and Physical Sciences Research Council 1
2 Outline Motivation Dynamic Optimization: Problem and Methods Convex Relaxation of the Dynamic Information Deterministic Global Optimization Algorithm Convergence of the Algorithm Case Studies Conclusions and Perspectives 2
3 Applications of Dynamic Optimization Chemical systems Determination of Kinetic Constants from Time Series Data (Parameter Estimation) Optimal Control of Batch and semi-batch Chemical Reactors Safety Analysis of Industrial Processes Biological and ecological systems Economic and other dynamic systems 3
4 Dynamic Optimization Problem min p J(x(t i, p), p ; i = 0, 1,..., NP ) subject to: g i (x(t i, p), p) 0, i = 0, 1,..., NP p L p p U Solution approaches Simultaneous approach: full discretization Sequential approach: p 0 NLP Algorithm p p opt ẋ = f(t, x, p), t [t 0, t NP ] x(t 0, p) = x 0 (p) Integrator x,x p 4
5 Usually nonconvex NLP Problems Multiple optimum solutions exist Commercially available numerical solvers guarantee local optimum solutions Poor economic performance Algorithm classes (combined with simultaneous or sequential approach) Stochastic algorithms Deterministic algorithms 5
6 Global optimization methods Simultaneous approaches Application of global optimization algorithms for NLPs Issues: problem size; quality of discretization. Smith and Pantelides (1996): spatial BB + reformulation Esposito and Floudas (2000): αbb algorithm 6
7 Global optimization methods: Sequential approaches Stochastic algorithms Luus et al. (1990): direct search procedure. Banga and Seider (1996), Banga et al. (1997): randomly directed search. Deterministic algorithms New techniques for convex relaxation of time-dependent parts of problem Lack of analytical forms for the constraints / objective function Esposito and Floudas (2000): extension of αbb to handle nonlinear dynamics. Singer and Barton (2002): convex relaxation of integral objective function with linear dynamics Papamichail and Adjiman (2002): convex relaxations of nonlinear dynamics. 7
8 Outline Motivation Dynamic Optimization: Problem and Methods Convex Relaxation of the Dynamic Information Deterministic Global Optimization Algorithm Convergence of the Algorithm Case Studies Conclusions and Perspectives 8
9 Reformulated NLP Problem subject to: min ˆx,p J(ˆx i, p ; i = 0, 1,..., NP ) g i (ˆx i, p) 0, i = 0, 1,..., NP ˆx i = x(t i, p), i = 0, 1,..., NP p [p L, p U ] ẋ = f(t, x, p), t [t 0, t NP ] x(t 0, p) = x 0 (p) 9
10 Convex Relaxation of J and g i (1) f(z) = f CT (z) + bt i=1 b iz Bi,1z Bi,2 + ut i=1 f UT,i(z i ) + nt i=1 f NT,i(z) Underestimating Bilinear Terms (McCormick, 1976) w = z 1 z 2 over [z L 1, zu 1 ] [zl 2, zu 2 ] w z L 1 z 2 + z L 2 z 1 z L 1 zl 2 w z U 1 z 2 + z U 2 z 1 z U 1 zu 2 w z L 1 z 2 + z U 2 z 1 z L 1 zu 2 w z U 1 z 2 + z L 2 z 1 z U 1 zl 2 Underestimating Univariate Concave Terms f UT (z) = f UT (z L ) + f UT (z U ) f UT (z L ) z U z L (z z L ) over [z L, z U ] R 10
11 Convex Relaxation of J and g i (2) Underestimating General Nonconvex Terms in C 2 (Maranas and Floudas, 1994; Androulakis et al., 1995) f NT (z) = f NT (z) + m i=1 α i (z L i z i )(z U i z i ) over [z L, z U ] R m f NT (z) is always less than f NT f NT (z) is convex if α i is big enough H f NT (z) = H fnt (z) + 2 diag(α i ) Rigorous α calculations using the scaled Gerschgorin method (Adjiman et al., 1998) z [z L, z U ] H fnt (z) [H fnt ] = H fnt ([z L, z U ]) 11
12 Convex Relaxation of ˆx i = x(t i, p) ˆx i x(t i, p) 0 x(t i, p) ˆx i 0 Constant bounds: x(t i ) ˆx i x(t i ) Affine bounds: M(t i )p + N(t i ) ˆx i M(t i )p + N(t i ) α-based bounds (Esposito and Floudas, 2000): ˆx ik x k (t i, p) + x k (t i, p) + r j=1 r j=1 α kij (pl j p j)(p U j p j) 0 α + kij (pl j p j)(p U j p j) ˆx ik 0 12
13 Illustrative example ẋ(t) = x(t) 2 + p, t [0, 1] x(0) = 9 p [ 5, 5] 10 8 p L = 5 p U =5 x p=5 p=2 p=0 p= 2 How do we find bounding trajectories? t p= 5 13
14 Differential Inequalities (1) Consider the following parameter dependent ODE: ẋ = f(t, x, p), t [t 0, t NP ] x(t 0, p) = x 0 (p) p [p L, p U ] R r where x and ẋ R n, f : (t 0, t NP ] R n [p L, p U ] R n and x 0 : [p L, p U ] R n. Let x = (x 1, x 2,..., x n ) T and x k = (x 1, x 2,..., x k 1, x k+1,..., x n ) T. The notation f(t, x, p) = f(t, x k, x k, p) is used. Theory on diff. inequalities (Walter, 1970) has been extended. 14
15 Differential Inequalities (2) Bounds on the solutions of the parameter dependent ODE: ẋ k = inf f k (t, x k, [x k, x k ], [p L, p U ]) ẋ k = sup f k (t, x k, [x k, x k ], [p L, p U ]) x(t 0 ) = inf x 0 ([p L, p U ]) x(t 0 ) = sup x 0 ([p L, p U ]) t [t 0, t NP ] and k = 1,..., n x(t) is a subfunction and x(t) is a superfunction for the solution of the ODE, i.e., x(t) x(t, p) x(t), p [p L, p U ], t [t 0, t NP ] 15
16 Quasi-monotonicity Definition 1: Let g(x) be a mapping g : D R with D R n. Again the notation g(x) = g(x k, x k ) is used. The function g is called unconditionally partially isotone (antitone) on D with respect to the variable x k if g(x k, x k ) g( x k, x k ) for x k x k (x k x k ) and for all (x k, x k ), ( x k, x k ) D. Definition 2: Let f(t, x) = (f 1 (t, x),..., f 2 (t, x)) T and each f k (t, x k, x k ) be unconditionally partially isotone on I 0 R R n 1 with respect to any component of x k, but not necessarily with respect to x k. Then f is quasi-monotone increasing on I 0 R n with respect to x (Walter, 1970) 16
17 Example: Constant bounds ẋ(t) = x(t) 2 5 ẋ(t) = x(t) x(0) = 9 x(0) = p L = 5 p U =5 x 4 2 x 0 2 x t 17
18 Parameter Dependent Bounds Let f(t, x, p) f(t, x, p) x [x(t), x(t)], p [p L, p U ], t [t 0, t NP ] and x 0 (p) x 0 (p) p [p L, p U ], where f : [t 0, t NP ] R n [p L, p U ] R n and x 0 : [p L, p U ] R n. If f is quasi-monotone increasing w.r.t. x and x(t, p) is the solution of the ODE: ẋ = f(t, x, p), t [t 0, t NP ] x(t 0, p) = x 0 (p) p [p L, p U ] then x(t, p) x(t, p), p [p L, p U ], t [t 0, t NP ]. 18
19 Affine Bounds Let f(t, x, p) = A(t)x + B(t)p + C(t) and x 0 (p) = Dp + E, where A(t), B(t) and C(t) are continuous on [t 0, t NP ]. Then the analytical solution is (Zadeh and Desoer, 1963): x(t, p) = { Φ(t, t 0 )D + +Φ(t, t 0 )E + where Φ(t, t 0 ) is the transition matrix: t t 0 Φ(t, τ)b(τ)dτ t } t 0 Φ(t, τ)c(τ)dτ, Φ(t, t 0 ) = A(t)Φ(t, t 0 ) t [t 0, t NP ] Φ(t 0, t 0 ) = I and I is the identity matrix. x(t, p) = M(t)p + N(t). p 19
20 M(t i ), N(t i ) calculation 1. Apply x(t i, p) = M(t i )p + N(t i ) for r + 1 values of p 2. Calculate x(t i, p) for the r + 1 values of p from the integration of the linear ODE 3. Solve n linear systems to find the r+1 unknowns for each one of the n dimensions of x 20
21 Underestimating IVP Overestimating IVPs Example: Affine bounds ẋ = (x + x)x + xx + v t [0, 1] x(0, v) = 9 ẋ 1 = 2xx 1 + x 2 + v t [0, 1] x 1 (0, v) = 9 ẋ 2 = 2xx 2 + x 2 + v t [0, 1] x 2 (0, v) = 9 21
22 Example: Affine bounds for p = p L = 5 p U =5 x(t,p=0) 4 2 x 0 x t 22
23 α calculation for x k (t i, p) [H xk (t i ) ] H x k (t i ) (p) = 2 x k (t i, p), p [p L, p U ] R r i = 0, 1,..., NS, k = 1, 2,..., n 1. 1st and 2nd order sensitivity equations 2. Create bounds using Differential Inequalities 3. Construct the interval Hessian matrix 4. Calculate α using the scaled Gerschgorin method 23
24 Example: All bounds x(t=1,p) 3 6 x(t=1,p) p p 24
25 Outline Motivation Dynamic Optimization: Problem and Methods Convex Relaxation of the Dynamic Information Deterministic Global Optimization Algorithm Convergence of the Algorithm Case Studies Conclusions and Perspectives 25
26 Spatial BB Algorithm (Horst and Tuy, 1996) Initialisation Upper and Lower Bounds Upper Bound Branch NO Lower Bound Subregion Selection u J - J l J R l R < ε r YES Terminate 26
27 Global optimization algorithm Step 1 Initialization: empty list of subregions, bounds on solution Step 2 First upper bound calculation (local optimization) Step 3 First lower bound calculation, including relaxation. Add subregions to list. Step 4 Subregion selection Terminate if list is empty Choose region with lowest lower bound otherwise Step 5 Check for convergence (relative tolerance, max iter) Step 6 Branch with standard rule (least reduced axis) Step 7 Upper bound for new regions (not needed at every iteration) Step 8 Lower bound calculation for each region. Go to Step 4. 27
28 Lower bound calculation for region R (Step 8) Let J L be the lower bound on the parent region of R. Let J U be the best known upper bound. Obtain bounds on the differential variables If affine bounds are used, obtain necessary matrices Form convex relaxation of problem for region R If a feasible solution with objective function JR L is obtained, then If affine bounds are used and JR L < JL, set JR L = JL If JR L JU, add R to the list of subregions 28
29 Outline of proof of convergence Three main properties are needed: Bound improvement after branching Constant bounds improve after branching α-based bounds improve after branching Affine bounds do not improve after branching. Ensure improvement through test in Step 8: if J L R < JL, set J L R = JL Bound improving selection operation Consequence of region selection criterion (Step 6) Consistent bounding operation Maximum distance between objective function and its relaxation converges to zero. Maximum distance between any constraint and its relaxation converges to zero. 29
30 Key elements of proof Bounds on the solutions of the ODE are such that: ẋ k = inf f k (t, x k, [x k, x k ], [p L, p U ]) inf f k (t, [x, x], [p L, p U ]) ẋ k = sup f k (t, x k, [x k, x k ], [p L, p U ]) sup f k (t, [x, x], [p L, p U ]) t [t 0, t NP ] and k = 1,..., n Inclusion monotonicity of interval operations ensures consistency of bounding operation with constant bounds. Similar approach can be taken to show α-based underestimators yield a consistent bounding operation: Interval Hessian matrix obtained through differential inequalities has desired properties. Hence, α and α-based bounds have desired properties. 30
31 Implementation of algorithm MATLAB 5.3 implementation NLPs: Use fmincon function (Optimization Toolbox) IVP solution: ode45 (Runge-Kutta based on Dormand-Prince pair) Interval calculations: INTLAB with directed outward rounding. Runs performed on an Ultra 60 workstation. 31
32 Case Study I: Parameter estimation for 1 st order reaction (Tjoa and Biegler, 1991) A k 1 B k 2 C subject to: min k 1,k j=1 i=1 (x i (t j ) x exp i (t j )) 2 ẋ 1 = k 1 x 1 t [0, 1] ẋ 2 = k 1 x 1 k 2 x 2 x 1 (0) = 1 x 2 (0) = 0 0 k k 2 10 Up to 8 affine underestimators and 8 affine overestimators can be constructed. α-values 0.5. Convexity is identified. 32
33 Results for case study I Obj. fun. = e-06; k 1 = ; k 2 = Underestimation Branching CPU time ɛ r Iter. scheme strategy (sec) C e-02 3,501 2,828 C e-03 34,508 22,959 C & A e C & A e C & α e C & α e C & α e C & α e C & A & α e C & A & α e
34 Case Study II: Parameter estimation for catalytic cracking of gas oil (Tjoa and Biegler, 1991) A k 1 Q k 3 k 2 S subject to: min k 1,k 2,k j=1 i=1 (x i (t j ) x exp i (t j )) 2 ẋ 1 = (k 1 + k 3 )x 21 t [0, 0.95] ẋ 2 = k 1 x 2 1 k 2x 2 x 1 (0) = 1 x 2 (0) = 0 0 k k k
35 Results for case study II Obj. fun. = e 03; k = ( , , ) T Underestimation Branching CPU time ɛ r Iter. scheme strategy (sec) C e-02 10,000 16,729 C e , ,816 C & A e ,597 C & A e ,478 C & α e ,415 C & α e ,524 C & α e ,116 C & α e , affine underestimators + 64 affine overestimators 35
36 Case Study III: Parameter estimation for reversible gas phase reaction (Bellman, 1967): 2NO + O 2 2NO 2 min k 1,k 2 14 j=1 (x(t = t j, k 1, k 2 ) x exp (t j )) 2 s.t. ẋ = k 1 (126.2 x)(91.9 x) 2 k 2 x 2 t [0, 39] x(t = 0, k 1, k 2 ) = 0 0 k k x is the difference of the pressure of the system from the initial pressure. k 1 and k 2 are the rate constants of the forward and reverse reactions. 36
37 Results for case study III Obj. fun.= ; k 1 = e-06; k 2 = e-04 Underestimation Branching CPU time ɛ r Iter. scheme strategy (sec) constant e-02 75,441 58,513 only constant bounds used 32 affine underestimators affine overestimators stiff systems, expensive to integrate quality of bounds not better than constant bounds quality of α-based bounds poor due to wrapping effect 37
38 Case Study IV: Optimal control with end-point constraint (Goh and Teo, 1988) Problem formulation using control vector parameterization: min u 1,u 2 x 2 (t = 1, u 1, u 2 ) s.t. ẋ 1 = u 1 (1 t) + u 2 t t ẋ 2 = x (u 1(1 t) + u 2 t) 2 x 1 (t = 0, u 1, u 2 ) = 1 x 2 (t = 0, u 1, u 2 ) = 0 x 1 (t = 1, u 1, u 2 ) 1 x 1 (t = 1, u 1, u 2 ) 1 1 u u 2 1 [0, 1] 16 affine underestimators + 2 affine overestimators 38
39 Results for case study IV Obj. fun.= e-01; u 1 = ; u 2 = Underestimation Branching CPU time ɛ r Iter. scheme strategy (sec) C e C e-03 1,062 1,106 C & A e C & A e C & α 1 or e α-based bounds recognize convexity of problem at root node. 39
40 Conclusions Three types of rigorous convex relaxations have been developed Convergence of the algorithm has been proved A BB global optimization algorithm has been applied successfully to case studies in parameter estimation and optimal control References: Papamichail and Adjiman, J. Glob Opt, Papamichail and Adjiman, Comp Chem Eng,
41 Perspectives Basic theoretical developments of recent years make global optimization of problems with nonlinear IVPs in the constraints possible. Practical applicability limited by cost of constructing underestimators and overestimators, quality of estimators for highly nonlinear systems (e.g. oscillatory) and long time horizons. Further research needed to identify classes of IVPs for which current estimators are effective, develop new estimators for other problem classes, establish basic theory for DAE systems. 41
Deterministic Global Optimization for Dynamic Systems Using Interval Analysis
Deterministic Global Optimization for Dynamic Systems Using Interval Analysis Youdong Lin and Mark A. Stadtherr Department of Chemical and Biomolecular Engineering University of Notre Dame, Notre Dame,
More informationImplementation of an αbb-type underestimator in the SGO-algorithm
Implementation of an αbb-type underestimator in the SGO-algorithm Process Design & Systems Engineering November 3, 2010 Refining without branching Could the αbb underestimator be used without an explicit
More informationDeterministic Global Optimization of Nonlinear Dynamic Systems
Deterministic Global Optimization of Nonlinear Dynamic Systems Youdong Lin and Mark A. Stadtherr Department of Chemical and Biomolecular Engineering University of Notre Dame, Notre Dame, IN 46556, USA
More informationNonlinear and robust MPC with applications in robotics
Nonlinear and robust MPC with applications in robotics Boris Houska, Mario Villanueva, Benoît Chachuat ShanghaiTech, Texas A&M, Imperial College London 1 Overview Introduction to Robust MPC Min-Max Differential
More information1 Introduction A large proportion of all optimization problems arising in an industrial or scientic context are characterized by the presence of nonco
A Global Optimization Method, BB, for General Twice-Dierentiable Constrained NLPs I. Theoretical Advances C.S. Adjiman, S. Dallwig 1, C.A. Floudas 2 and A. Neumaier 1 Department of Chemical Engineering,
More informationMixed Integer Non Linear Programming
Mixed Integer Non Linear Programming Claudia D Ambrosio CNRS Research Scientist CNRS & LIX, École Polytechnique MPRO PMA 2016-2017 Outline What is a MINLP? Dealing with nonconvexities Global Optimization
More informationMcCormick Relaxations: Convergence Rate and Extension to Multivariate Outer Function
Introduction Relaxations Rate Multivariate Relaxations Conclusions McCormick Relaxations: Convergence Rate and Extension to Multivariate Outer Function Alexander Mitsos Systemverfahrenstecnik Aachener
More informationA Decomposition-based MINLP Solution Method Using Piecewise Linear Relaxations
A Decomposition-based MINLP Solution Method Using Piecewise Linear Relaxations Pradeep K. Polisetty and Edward P. Gatzke 1 Department of Chemical Engineering University of South Carolina Columbia, SC 29208
More informationA global optimization method, αbb, for general twice-differentiable constrained NLPs I. Theoretical advances
Computers Chem. Engng Vol. 22, No. 9, pp. 1137 1158, 1998 Copyright 1998 Elsevier Science Ltd All rights reserved. Printed in Great Britain PII: S0098-1354(98)00027-1 0098 1354/98 $19.00#0.00 A global
More informationRobust Parameter Estimation in Nonlinear Dynamic Process Models
European Symposium on Computer Arded Aided Process Engineering 15 L. Puigjaner and A. Espuña (Editors) 2005 Elsevier Science B.V. All rights reserved. Robust Parameter Estimation in Nonlinear Dynamic Process
More informationGlobal Optimization by Interval Analysis
Global Optimization by Interval Analysis Milan Hladík Department of Applied Mathematics, Faculty of Mathematics and Physics, Charles University in Prague, Czech Republic, http://kam.mff.cuni.cz/~hladik/
More informationMixed-Integer Nonlinear Programming
Mixed-Integer Nonlinear Programming Claudia D Ambrosio CNRS researcher LIX, École Polytechnique, France pictures taken from slides by Leo Liberti MPRO PMA 2016-2017 Motivating Applications Nonlinear Knapsack
More informationSolving Mixed-Integer Nonlinear Programs
Solving Mixed-Integer Nonlinear Programs (with SCIP) Ambros M. Gleixner Zuse Institute Berlin MATHEON Berlin Mathematical School 5th Porto Meeting on Mathematics for Industry, April 10 11, 2014, Porto
More informationNonlinear convex and concave relaxations for the solutions of parametric ODEs
OPTIMAL CONTROL APPLICATIONS AND METHODS Optim. Control Appl. Meth. 0000; 00:1 22 Published online in Wiley InterScience (www.interscience.wiley.com). Nonlinear convex and concave relaxations for the solutions
More informationSolving Global Optimization Problems by Interval Computation
Solving Global Optimization Problems by Interval Computation Milan Hladík Department of Applied Mathematics, Faculty of Mathematics and Physics, Charles University in Prague, Czech Republic, http://kam.mff.cuni.cz/~hladik/
More informationA Branch-and-Refine Method for Nonconvex Mixed-Integer Optimization
A Branch-and-Refine Method for Nonconvex Mixed-Integer Optimization Sven Leyffer 2 Annick Sartenaer 1 Emilie Wanufelle 1 1 University of Namur, Belgium 2 Argonne National Laboratory, USA IMA Workshop,
More informationHot-Starting NLP Solvers
Hot-Starting NLP Solvers Andreas Wächter Department of Industrial Engineering and Management Sciences Northwestern University waechter@iems.northwestern.edu 204 Mixed Integer Programming Workshop Ohio
More informationChapter 2 Optimal Control Problem
Chapter 2 Optimal Control Problem Optimal control of any process can be achieved either in open or closed loop. In the following two chapters we concentrate mainly on the first class. The first chapter
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationCHAPTER 2: QUADRATIC PROGRAMMING
CHAPTER 2: QUADRATIC PROGRAMMING Overview Quadratic programming (QP) problems are characterized by objective functions that are quadratic in the design variables, and linear constraints. In this sense,
More informationExact Computation of Global Minima of a Nonconvex Portfolio Optimization Problem
Frontiers In Global Optimization, pp. 1-2 C. A. Floudas and P. M. Pardalos, Editors c 2003 Kluwer Academic Publishers Exact Computation of Global Minima of a Nonconvex Portfolio Optimization Problem Josef
More informationDevelopment of an algorithm for solving mixed integer and nonconvex problems arising in electrical supply networks
Development of an algorithm for solving mixed integer and nonconvex problems arising in electrical supply networks E. Wanufelle 1 S. Leyffer 2 A. Sartenaer 1 Ph. Toint 1 1 FUNDP, University of Namur 2
More informationA Test Set of Semi-Infinite Programs
A Test Set of Semi-Infinite Programs Alexander Mitsos, revised by Hatim Djelassi Process Systems Engineering (AVT.SVT), RWTH Aachen University, Aachen, Germany February 26, 2016 (first published August
More informationA software toolbox for the dynamic optimization of nonlinear processes
European Symposium on Computer Arded Aided Process Engineering 15 L. Puigjaner and A. Espuña (Editors) 2005 Elsevier Science B.V. All rights reserved. A software toolbox for the dynamic optimization of
More informationN. L. P. NONLINEAR PROGRAMMING (NLP) deals with optimization models with at least one nonlinear function. NLP. Optimization. Models of following form:
0.1 N. L. P. Katta G. Murty, IOE 611 Lecture slides Introductory Lecture NONLINEAR PROGRAMMING (NLP) deals with optimization models with at least one nonlinear function. NLP does not include everything
More informationTechnische Universität Ilmenau Institut für Mathematik
Technische Universität Ilmenau Institut für Mathematik Preprint No. M 18/03 A branch-and-bound based algorithm for nonconvex multiobjective optimization Julia Niebling, Gabriele Eichfelder Februar 018
More informationOptimisation: From theory to better prediction
Optimisation: From theory to better prediction Claire Adjiman Claire Adjiman 23 May 2013 Christodoulos A. Floudas Memorial Symposium 6 May 2017 Why look for the global solution? Deterministic global optimisation
More informationFrom structures to heuristics to global solvers
From structures to heuristics to global solvers Timo Berthold Zuse Institute Berlin DFG Research Center MATHEON Mathematics for key technologies OR2013, 04/Sep/13, Rotterdam Outline From structures to
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Mathematics 234 (2) 538 544 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationA Two-Stage Algorithm for Multi-Scenario Dynamic Optimization Problem
A Two-Stage Algorithm for Multi-Scenario Dynamic Optimization Problem Weijie Lin, Lorenz T Biegler, Annette M. Jacobson March 8, 2011 EWO Annual Meeting Outline Project review and problem introduction
More informationmin f(x). (2.1) Objectives consisting of a smooth convex term plus a nonconvex regularization term;
Chapter 2 Gradient Methods The gradient method forms the foundation of all of the schemes studied in this book. We will provide several complementary perspectives on this algorithm that highlight the many
More informationOptimization-based Domain Reduction in Guaranteed Parameter Estimation of Nonlinear Dynamic Systems
9th IFAC Symposium on Nonlinear Control Systems Toulouse, France, September 4-6, 2013 ThC2.2 Optimization-based Domain Reduction in Guaranteed Parameter Estimation of Nonlinear Dynamic Systems Radoslav
More informationReliable Modeling Using Interval Analysis: Chemical Engineering Applications
Reliable Modeling Using Interval Analysis: Chemical Engineering Applications Mark A. Stadtherr Department of Chemical Engineering University of Notre Dame Notre Dame, IN 46556 USA SIAM Workshop on Validated
More informationHeuristics for nonconvex MINLP
Heuristics for nonconvex MINLP Pietro Belotti, Timo Berthold FICO, Xpress Optimization Team, Birmingham, UK pietrobelotti@fico.com 18th Combinatorial Optimization Workshop, Aussois, 9 Jan 2014 ======This
More informationPositive Semidefiniteness and Positive Definiteness of a Linear Parametric Interval Matrix
Positive Semidefiniteness and Positive Definiteness of a Linear Parametric Interval Matrix Milan Hladík Charles University, Faculty of Mathematics and Physics, Department of Applied Mathematics, Malostranské
More informationNonlinear Optimization for Optimal Control
Nonlinear Optimization for Optimal Control Pieter Abbeel UC Berkeley EECS Many slides and figures adapted from Stephen Boyd [optional] Boyd and Vandenberghe, Convex Optimization, Chapters 9 11 [optional]
More informationFast Probability Generating Function Method for Stochastic Chemical Reaction Networks
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. 71 (2014) 57-69 ISSN 0340-6253 Fast Probability Generating Function Method for Stochastic Chemical Reaction
More informationICRS-Filter: A Randomized Direct Search Algorithm for. Constrained Nonconvex Optimization Problems
ICRS-Filter: A Randomized Direct Search Algorithm for Constrained Nonconvex Optimization Problems Biyu Li 1, Viet H. Nguyen 1, Chieh. L. Ng 1, E. A. del Rio-Chanona 1 Vassilios S. Vassiliadis 1, Harvey
More informationAn exact reformulation algorithm for large nonconvex NLPs involving bilinear terms
An exact reformulation algorithm for large nonconvex NLPs involving bilinear terms Leo Liberti CNRS LIX, École Polytechnique, F-91128 Palaiseau, France. E-mail: liberti@lix.polytechnique.fr. Constantinos
More informationDisconnecting Networks via Node Deletions
1 / 27 Disconnecting Networks via Node Deletions Exact Interdiction Models and Algorithms Siqian Shen 1 J. Cole Smith 2 R. Goli 2 1 IOE, University of Michigan 2 ISE, University of Florida 2012 INFORMS
More informationTools for dynamic model development. Spencer Daniel Schaber
Tools for dynamic model development by Spencer Daniel Schaber B.Ch.E., University of Minnesota (2008) M.S.C.E.P., Massachusetts Institute of Technology (2011) Submitted to the Department of Chemical Engineering
More informationAlpha-helical Topology and Tertiary Structure Prediction of Globular Proteins Scott R. McAllister Christodoulos A. Floudas Princeton University
Alpha-helical Topology and Tertiary Structure Prediction of Globular Proteins Scott R. McAllister Christodoulos A. Floudas Princeton University Department of Chemical Engineering Program of Applied and
More informationDeterministic Global Optimization for Error-in-Variables Parameter Estimation
Deterministic Global Optimization for Error-in-Variables Parameter Estimation Chao-Yang Gau and Mark A. Stadtherr Department of Chemical Engineering University of Notre Dame 182 Fitzpatrick Hall Notre
More informationOptimizing Economic Performance using Model Predictive Control
Optimizing Economic Performance using Model Predictive Control James B. Rawlings Department of Chemical and Biological Engineering Second Workshop on Computational Issues in Nonlinear Control Monterey,
More informationSolving Models With Off-The-Shelf Software. Example Of Potential Pitfalls Associated With The Use And Abuse Of Default Parameter Settings
An Example Of Potential Pitfalls Associated With The Use And Abuse Of Default Parameter Settings Ric D. Herbert 1 Peter J. Stemp 2 1 Faculty of Science and Information Technology, The University of Newcastle,
More informationNumerical Optimal Control Overview. Moritz Diehl
Numerical Optimal Control Overview Moritz Diehl Simplified Optimal Control Problem in ODE path constraints h(x, u) 0 initial value x0 states x(t) terminal constraint r(x(t )) 0 controls u(t) 0 t T minimize
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationLecture Note 13:Continuous Time Switched Optimal Control: Embedding Principle and Numerical Algorithms
ECE785: Hybrid Systems:Theory and Applications Lecture Note 13:Continuous Time Switched Optimal Control: Embedding Principle and Numerical Algorithms Wei Zhang Assistant Professor Department of Electrical
More informationGlobal optimization of minimum weight truss topology problems with stress, displacement, and local buckling constraints using branch-and-bound
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2004; 61:1270 1309 (DOI: 10.1002/nme.1112) Global optimization of minimum weight truss topology problems with stress,
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationTime-Optimal Automobile Test Drives with Gear Shifts
Time-Optimal Control of Automobile Test Drives with Gear Shifts Christian Kirches Interdisciplinary Center for Scientific Computing (IWR) Ruprecht-Karls-University of Heidelberg, Germany joint work with
More information1. Introduction. We consider the general global optimization problem
CONSTRUCTION OF VALIDATED UNIQUENESS REGIONS FOR NONLINEAR PROGRAMS IN WHICH CONVEX SUBSPACES HAVE BEEN IDENTIFIED R. BAKER KEARFOTT Abstract. In deterministic global optimization algorithms for constrained
More information1 Convexity, Convex Relaxations, and Global Optimization
1 Conveity, Conve Relaations, and Global Optimization Algorithms 1 1.1 Minima Consider the nonlinear optimization problem in a Euclidean space: where the feasible region. min f(), R n Definition (global
More informationHow to simulate in Simulink: DAE, DAE-EKF, MPC & MHE
How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE Tamal Das PhD Candidate IKP, NTNU 24th May 2017 Outline 1 Differential algebraic equations (DAE) 2 Extended Kalman filter (EKF) for DAE 3 Moving horizon
More informationTheory, Solution Techniques and Applications of Singular Boundary Value Problems
Theory, Solution Techniques and Applications of Singular Boundary Value Problems W. Auzinger O. Koch E. Weinmüller Vienna University of Technology, Austria Problem Class z (t) = M(t) z(t) + f(t, z(t)),
More informationConvex Optimization. Newton s method. ENSAE: Optimisation 1/44
Convex Optimization Newton s method ENSAE: Optimisation 1/44 Unconstrained minimization minimize f(x) f convex, twice continuously differentiable (hence dom f open) we assume optimal value p = inf x f(x)
More informationOn rigorous integration of piece-wise linear systems
On rigorous integration of piece-wise linear systems Zbigniew Galias Department of Electrical Engineering AGH University of Science and Technology 06.2011, Trieste, Italy On rigorous integration of PWL
More informationImproved quadratic cuts for convex mixed-integer nonlinear programs
Improved quadratic cuts for convex mixed-integer nonlinear programs Lijie Su a,b, Lixin Tang a*, David E. Bernal c, Ignacio E. Grossmann c a Institute of Industrial and Systems Engineering, Northeastern
More informationSTRONG VALID INEQUALITIES FOR MIXED-INTEGER NONLINEAR PROGRAMS VIA DISJUNCTIVE PROGRAMMING AND LIFTING
STRONG VALID INEQUALITIES FOR MIXED-INTEGER NONLINEAR PROGRAMS VIA DISJUNCTIVE PROGRAMMING AND LIFTING By KWANGHUN CHUNG A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN
More informationFeasibility Pump for Mixed Integer Nonlinear Programs 1
Feasibility Pump for Mixed Integer Nonlinear Programs 1 Presenter: 1 by Pierre Bonami, Gerard Cornuejols, Andrea Lodi and Francois Margot Mixed Integer Linear or Nonlinear Programs (MILP/MINLP) Optimize
More informationReview for Exam 2 Ben Wang and Mark Styczynski
Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note
More informationApplications and algorithms for mixed integer nonlinear programming
Applications and algorithms for mixed integer nonlinear programming Sven Leyffer 1, Jeff Linderoth 2, James Luedtke 2, Andrew Miller 3, Todd Munson 1 1 Mathematics and Computer Science Division, Argonne
More informationComparing Convex Relaxations for Quadratically Constrained Quadratic Programming
Comparing Convex Relaxations for Quadratically Constrained Quadratic Programming Kurt M. Anstreicher Dept. of Management Sciences University of Iowa European Workshop on MINLP, Marseille, April 2010 The
More informationQUALIFYING EXAM IN SYSTEMS ENGINEERING
QUALIFYING EXAM IN SYSTEMS ENGINEERING Written Exam: MAY 23, 2017, 9:00AM to 1:00PM, EMB 105 Oral Exam: May 25 or 26, 2017 Time/Location TBA (~1 hour per student) CLOSED BOOK, NO CHEAT SHEETS BASIC SCIENTIFIC
More informationLarge Scale Semi-supervised Linear SVMs. University of Chicago
Large Scale Semi-supervised Linear SVMs Vikas Sindhwani and Sathiya Keerthi University of Chicago SIGIR 2006 Semi-supervised Learning (SSL) Motivation Setting Categorize x-billion documents into commercial/non-commercial.
More informationResearch Article A New Global Optimization Algorithm for Solving Generalized Geometric Programming
Mathematical Problems in Engineering Volume 2010, Article ID 346965, 12 pages doi:10.1155/2010/346965 Research Article A New Global Optimization Algorithm for Solving Generalized Geometric Programming
More informationQuasiconvex relaxations based on interval arithmetic
Linear Algebra and its Applications 324 (2001) 27 53 www.elsevier.com/locate/laa Quasiconvex relaxations based on interval arithmetic Christian Jansson Inst. für Informatik III, Technical University Hamburg-Harburg,
More informationSlovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS
Slovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS of the 18 th International Conference on Process Control Hotel Titris, Tatranská
More informationDistributionally Robust Convex Optimization
Distributionally Robust Convex Optimization Wolfram Wiesemann 1, Daniel Kuhn 1, and Melvyn Sim 2 1 Department of Computing, Imperial College London, United Kingdom 2 Department of Decision Sciences, National
More informationLecture 13: Constrained optimization
2010-12-03 Basic ideas A nonlinearly constrained problem must somehow be converted relaxed into a problem which we can solve (a linear/quadratic or unconstrained problem) We solve a sequence of such problems
More informationPrimal-Dual Interior-Point Methods for Linear Programming based on Newton s Method
Primal-Dual Interior-Point Methods for Linear Programming based on Newton s Method Robert M. Freund March, 2004 2004 Massachusetts Institute of Technology. The Problem The logarithmic barrier approach
More informationNonlinear Model Predictive Control Tools (NMPC Tools)
Nonlinear Model Predictive Control Tools (NMPC Tools) Rishi Amrit, James B. Rawlings April 5, 2008 1 Formulation We consider a control system composed of three parts([2]). Estimator Target calculator Regulator
More informationThe use of second-order information in structural topology optimization. Susana Rojas Labanda, PhD student Mathias Stolpe, Senior researcher
The use of second-order information in structural topology optimization Susana Rojas Labanda, PhD student Mathias Stolpe, Senior researcher What is Topology Optimization? Optimize the design of a structure
More informationSemidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization
Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization Instructor: Farid Alizadeh Author: Ai Kagawa 12/12/2012
More informationAn Inexact Newton Method for Nonlinear Constrained Optimization
An Inexact Newton Method for Nonlinear Constrained Optimization Frank E. Curtis Numerical Analysis Seminar, January 23, 2009 Outline Motivation and background Algorithm development and theoretical results
More informationConvex Optimization M2
Convex Optimization M2 Lecture 8 A. d Aspremont. Convex Optimization M2. 1/57 Applications A. d Aspremont. Convex Optimization M2. 2/57 Outline Geometrical problems Approximation problems Combinatorial
More informationOptimization of Batch Processes
Integrated Scheduling and Dynamic Optimization of Batch Processes Yisu Nie and Lorenz T. Biegler Center for Advanced Process Decision-making Department of Chemical Engineering Carnegie Mellon University
More informationLCPI method to find optimal solutions of nonlinear programming problems
ISSN 1 746-7233, Engl, UK World Journal of Modelling Simulation Vol. 14 2018) No. 1, pp. 50-57 LCPI method to find optimal solutions of nonlinear programming problems M.H. Noori Skari, M. Ghaznavi * Faculty
More informationEuler s Method applied to the control of switched systems
Euler s Method applied to the control of switched systems FORMATS 2017 - Berlin Laurent Fribourg 1 September 6, 2017 1 LSV - CNRS & ENS Cachan L. Fribourg Euler s method and switched systems September
More informationMixed-Integer Nonlinear Decomposition Toolbox for Pyomo (MindtPy)
Mario R. Eden, Marianthi Ierapetritou and Gavin P. Towler (Editors) Proceedings of the 13 th International Symposium on Process Systems Engineering PSE 2018 July 1-5, 2018, San Diego, California, USA 2018
More informationStochastic geometric optimization with joint probabilistic constraints
Stochastic geometric optimization with joint probabilistic constraints Jia Liu a,b, Abdel Lisser 1a, Zhiping Chen b a Laboratoire de Recherche en Informatique (LRI), Université Paris Sud - XI, Bât. 650,
More informationScenario Grouping and Decomposition Algorithms for Chance-constrained Programs
Scenario Grouping and Decomposition Algorithms for Chance-constrained Programs Siqian Shen Dept. of Industrial and Operations Engineering University of Michigan Joint work with Yan Deng (UMich, Google)
More informationThe Pooling Problem - Applications, Complexity and Solution Methods
The Pooling Problem - Applications, Complexity and Solution Methods Dag Haugland Dept. of Informatics University of Bergen Norway EECS Seminar on Optimization and Computing NTNU, September 30, 2014 Outline
More informationModels for Control and Verification
Outline Models for Control and Verification Ian Mitchell Department of Computer Science The University of British Columbia Classes of models Well-posed models Difference Equations Nonlinear Ordinary Differential
More informationComputer Sciences Department
Computer Sciences Department Valid Inequalities for the Pooling Problem with Binary Variables Claudia D Ambrosio Jeff Linderoth James Luedtke Technical Report #682 November 200 Valid Inequalities for the
More informationResearch Article Complete Solutions to General Box-Constrained Global Optimization Problems
Journal of Applied Mathematics Volume, Article ID 47868, 7 pages doi:.55//47868 Research Article Complete Solutions to General Box-Constrained Global Optimization Problems Dan Wu and Youlin Shang Department
More informationOptimization of Heat Exchanger Network with Fixed Topology by Genetic Algorithms
Optimization of Heat Exchanger Networ with Fixed Topology by Genetic Algorithms Roman Bochene, Jace Jeżowsi, Sylwia Bartman Rzeszow University of Technology, Department of Chemical and Process Engineering,
More informationECE7850 Lecture 9. Model Predictive Control: Computational Aspects
ECE785 ECE785 Lecture 9 Model Predictive Control: Computational Aspects Model Predictive Control for Constrained Linear Systems Online Solution to Linear MPC Multiparametric Programming Explicit Linear
More informationAn Inexact Newton Method for Optimization
New York University Brown Applied Mathematics Seminar, February 10, 2009 Brief biography New York State College of William and Mary (B.S.) Northwestern University (M.S. & Ph.D.) Courant Institute (Postdoc)
More informationThe HJB-POD approach for infinite dimensional control problems
The HJB-POD approach for infinite dimensional control problems M. Falcone works in collaboration with A. Alla, D. Kalise and S. Volkwein Università di Roma La Sapienza OCERTO Workshop Cortona, June 22,
More informationCourse Outline. FRTN10 Multivariable Control, Lecture 13. General idea for Lectures Lecture 13 Outline. Example 1 (Doyle Stein, 1979)
Course Outline FRTN Multivariable Control, Lecture Automatic Control LTH, 6 L-L Specifications, models and loop-shaping by hand L6-L8 Limitations on achievable performance L9-L Controller optimization:
More informationCyclic short-term scheduling of multiproduct batch plants using continuous-time representation
Computers and Chemical Engineering (00) Cyclic short-term scheduling of multiproduct batch plants using continuous-time representation D. Wu, M. Ierapetritou Department of Chemical and Biochemical Engineering,
More informationLecture: Convex Optimization Problems
1/36 Lecture: Convex Optimization Problems http://bicmr.pku.edu.cn/~wenzw/opt-2015-fall.html Acknowledgement: this slides is based on Prof. Lieven Vandenberghe s lecture notes Introduction 2/36 optimization
More informationEnhancing Energy Efficiency among Communication, Computation and Caching with QoI-Guarantee
Enhancing Energy Efficiency among Communication, Computation and Caching with QoI-Guarantee Faheem Zafari 1, Jian Li 2, Kin K. Leung 1, Don Towsley 2, Ananthram Swami 3 1 Imperial College London 2 University
More informationOn Rigorous Upper Bounds to a Global Optimum
Noname manuscript No. (will be inserted by the editor) On Rigorous Upper Bounds to a Global Optimum Ralph Baker Kearfott Received: date / Accepted: date Abstract In branch and bound algorithms in constrained
More informationGLOBAL OPTIMIZATION WITH GAMS/BARON
GLOBAL OPTIMIZATION WITH GAMS/BARON Nick Sahinidis Chemical and Biomolecular Engineering University of Illinois at Urbana Mohit Tawarmalani Krannert School of Management Purdue University MIXED-INTEGER
More informationWorst case analysis of mechanical structures by interval methods
Worst case analysis of mechanical structures by interval methods Arnold Neumaier (University of Vienna, Austria) (joint work with Andrzej Pownuk) Safety Safety studies in structural engineering are supposed
More informationAn RLT Approach for Solving Binary-Constrained Mixed Linear Complementarity Problems
An RLT Approach for Solving Binary-Constrained Mixed Linear Complementarity Problems Miguel F. Anjos Professor and Canada Research Chair Director, Trottier Institute for Energy TAI 2015 Washington, DC,
More informationConvex optimization problems. Optimization problem in standard form
Convex optimization problems optimization problem in standard form convex optimization problems linear optimization quadratic optimization geometric programming quasiconvex optimization generalized inequality
More informationCE 191: Civil & Environmental Engineering Systems Analysis. LEC 17 : Final Review
CE 191: Civil & Environmental Engineering Systems Analysis LEC 17 : Final Review Professor Scott Moura Civil & Environmental Engineering University of California, Berkeley Fall 2014 Prof. Moura UC Berkeley
More information