MODELING LANGUAGES in MATHEMATICAL OPTIMIZATION Applied Optimization Volume 88

Total Page:16

File Type:pdf, Size:1020Kb

MODELING LANGUAGES in MATHEMATICAL OPTIMIZATION Applied Optimization Volume 88 MODELING LANGUAGES IN MATHEMATICAL OPTIMIZATION Applied Optimization Volume 88 Series Editors: Panos M. Pardalos University o/Florida, U.S.A. Donald W. Hearn University o/Florida, U.S.A. MODELING LANGUAGES IN MATHEMATICAL OPTIMIZATION Edited by JOSEF KALLRATH BASF AG, GVC/S (Scientific Computing), 0-67056 Ludwigshafen, Germany Dept. of Astronomy, Univ. of Florida, Gainesville, FL 32611 Kluwer Academic Publishers Boston/DordrechtiLondon Distributors for North, Central and South America: Kluwer Academic Publishers 101 Philip Drive Assinippi Park Norwell, Massachusetts 02061 USA Telephone (781) 871-6600 Fax (781) 871-6528 E-Mail <[email protected]> Distributors for all other countries: Kluwer Academic Publishers Group Post Office Box 322 3300 AlI Dordrecht, THE NETHERLANDS Telephone 31 78 6576 000 Fax 31 786576474 E-Mail <[email protected]> .t Electronic Services <http://www.wkap.nl> Library of Congress Cataloging-in-Publication Kallrath, Josef Modeling Languages in Mathematical Optimization ISBN-13: 978-1-4613-7945-4 e-ISBN-13:978-1-4613 -0215 - 5 DOl: 10.1007/978-1-4613-0215-5 Copyright © 2004 by Kluwer Academic Publishers Softcover reprint of the hardcover 1st edition 2004 All rights reserved. No part ofthis pUblication may be reproduced, stored in a retrieval system or transmitted in any form or by any means, electronic, mechanical, photo-copying, microfilming, recording, or otherwise, without the prior written permission ofthe publisher, with the exception of any material supplied specifically for the purpose of being entered and executed on a computer system, for exclusive use by the purchaser of the work. Permissions for books published in the USA: P"§.;r.m.i.§J?i..QD.§.@:w.k~p" ..,.g.Qm Permissions for books published in Europe: [email protected] Printed on acid-free paper. Contents List of Figures xiii List of Tables xv Preface xvii Contributing Authors xxiii Introduction xxix Josef Kallrath Part I Theoretical and Practical Concepts of Modeling Languages 1 Mathematical Optimization and the Role of Modeling Languages 3 Josef Kallrath 1.1 Mathematical Optimization 3 1.2 Classes of Problems in Mathematical Optimization 6 1.2.1 A Deterministic Standard MINLP Problem 7 1.2.2 Constraint Satisfaction Problems 9 1.2.3 Multi-Objective Optimization 9 1.2.4 Multi-Level Optimization 10 1.2.5 Semi-Infinite Programming 11 1.2.6 Optimization Involving Differential Equations 12 1.2.7 Safety Programming 12 1.2.8 Optimization Under Uncertainty 13 1.2.8.1 Approaches to Optimization Under Uncertainty 13 1.2.8.2 Stochastic Optimization 15 1.2.8.3 Beyond Stochastic Programming 16 1.3 The History of Modeling Languages in Optimization 17 1.4 Conventions and Abbreviations 24 2 Models and the History of Modeling 25 Hermann Schichl 2.1 The History of Modeling 25 2.2 Models 28 2.3 Mathematical Models 29 2.4 The Modeling Process 31 2.4.1 The Importance of Good Modeling Practice 32 2.4.2 Making Mathematical Models Accessible for Computers 35 VI MODELING LANGUAGES IN MATHEMATICAL OPTIMIZATION 3 Mathematical Model Building 37 Arnold Neumaier 3.1 Why Mathematical Modeling? 37 3.2 A List of Applications 38 3.3 Basic numerical tasks 39 3.4 The Modeling Diagram 40 3.5 General Rules 41 3.6 Conflicts 42 3.7 Attitudes 43 4 Theoretical Concepts and Design of Modeling Languages 45 Hermann Schichl 4.1 Modeling Languages 45 4.1.1 Algebraic Modeling Languages 46 4.1.2 Non-algebraic Modeling Languages 51 4.1.3 Integrated Modeling Environments 54 4.1.4 Model-Programming Languages 55 4.1.5 Other Modeling Tools 55 4.2 Global Optimization 56 4.2.1 Problem Description 57 4.2.2 Algebraic Modeling Languages and Global Optimization 58 4.3 A Vision - What the Future Needs to Bring 59 4.3.1 Data Handling 59 4.3.2 Solver Views 60 4.3.3 GUI 60 4.3.4 Object Oriented Modeling - Derived Models 61 4.3.5 Hierarchical Modeling 61 4.3.6 Building Blocks 62 4.3.7 Open Model Exchange Format 62 5 The Importance of Modeling Languages for Solving Real-World Problems 63 Josef Kallrath and Josef Liesenfeld 5.1 Modeling Languages and Real World Problems 64 5.2 Requirements from Practitioners towards Modeling Languages and Modeling Systems 65 Part II The Modeling Languages in Detail 6 The Modeling Language AIMMS 71 Johannes Bisschop and Marcel Roelofs 6.1 AIMMS Design Philosophy, Features and Benefits 72 6.2 AIMMS Outer Approximation (AOA) Algorithm 78 6.2.1 Problem Statement 79 6.2.2 Basic Algorithm 79 6.2.3 Open Solver Approach 81 Contents vii 6.204 Alternative Uses of the Open Approach 83 6.3 Units of Measurement 84 6.3.1 Unit Analysis 86 6.3.2 Unit-Based Scaling 87 6.3.3 Unit Conventions 87 604 Time-Based Modeling 88 604.1 Calendars 90 604.2 Horizons 90 604.3 Data Conversion of Time-Dependent Identifiers 92 6.5 The AIMMS Excel Interface 93 6.5.1 Excel as the Main Application 94 6.5.2 AIMMS as the Main Application 94 6.6 Multi-Agent Support 96 6.6.1 Basic Agent Concepts 96 6.6.2 Examples of Motivation 97 6.6.3 Agent-Related Concepts in AIMMS 99 6.604 Agent Construction Support 101 6.7 Future Developments 101 Appendix 103 A AIMMS Features Overview 103 A.l Language Features 103 A.2 Mathematical Programming Features 103 A.3 End-User Interface Features 103 AA Connectivity and Deployment Features 104 B Application Examples 104 7 Design Principles and New Developments 105 in the AMPL Modeling Language Robert Fourer, David M. Gay and Brian W. Kernighan 7.1 Background and Early History 106 7.2 The McDonald's Diet Problem 107 7.3 The Airline Fleet Assignment Problem 111 704 Iterative Schemes 114 704.1 Flow of Control 114 704.2 Named Subproblems 115 704.3 Debugging 118 7.5 Other Types of Models 118 7.5.1 Piecewise-Linear Terms 119 7.5.2 Complementarity Problems 120 7.5.3 Combinatorial Optimization 123 7.504 Stochastic Programming 126 7.6 Communicating with Other Systems 128 7.6.1 Relational Database Access 128 7.6.2 Internet Optimization Services 129 7.6.3 Communication with Solvers via Suffixes 133 7.7 Updated AMPL Book 135 7.8 Concluding Remarks 135 V111 MODELING LANGUAGES IN MATHEMATICAL OPTIMIZATION 8 General Algebraic Modeling System (GAMS) 137 Michael R. Bussieck & Alex Meeraus 8.1 Background and Motivation 138 8.2 Design Goals and Changing Focus 139 8.3 A User's View of Modeling Languages 141 8.3.1 Academic Research Models 141 8.3.2 Domain Expert Models 145 8.3.3 Black Box Models 147 8.4 Summary and Conclusion 149 Appendix 150 A Selected Language Features 150 B GAMS External Functions 151 C Secure Work Files 151 D GAMS versus FORTRAN Matrix Generators 152 E Sample GAMS Problem 153 9 The LINGO Algebraic Modeling Language 159 Kevin Cunningham and Linus Schrage 9.1 History 159 9.2 Design Philosophy 160 9.2.1 Simplified Syntax for Small Models 160 9.2.2 Close Coupled Solvers 162 9.2.3 Interface to Excel 163 9.2.4 Model Class Identification 166 9.2.5 Automatic Linearization and Global Optimization 167 9.2.6 Debugging Models 168 9.2.7 Programming Interface 170 9.3 Future Directions 171 10 The LPL Modeling Language 173 Tony Hiirlimann 10.1 History 173 10.2 Some Basic Ideas 175 10.3 Highlights 176 10.4 The Cutting Stock Problem 178 10.5 Liquid Container 181 10.6 Model Documentation 182 10.7 Conclusion 183 11 The MINOPT Modeling Language 185 Carl A. Schweiger and Christodoulos A. Floudas 11.1 Introduction 186 11.1.1 Motivation 186 11.1.2 MINOPT Overview 188 Contents ix 11.2 Model Types and Solution Algorithms 189 11.2.1 Mixed-Integer Nonlinear Program (MINLP) 190 11.2.1.1 Generalized Benders Decomposition (GBD) 190 11.2.1.2 Outer ApproximationlEquality Relaxation! Augmented Penalty (OA/ER/AP) 192 11.2.2 Nonlinear Program with Differential and Algebraic Constraints (NLPIDAE) 195 11.2.3 Mixed-Integer Nonlinear Program with Differential and Al- gebraic Constraints (MINLPIDAE) 197 11.2.4 Optimal Control Problem (OCP) and Mixed Integer Optimal Control 200 11.2.5 External Solvers 204 11.3 Example Problems 204 11.3.1 Language Overview 205 11.3.2 MINLP Problem-Nonconvex Portfolio Optimization Prob- lem 205 11.3.3 Optimal Control Problem-Dow Batch Reactor 207 11.4 Summary 209 12 Mosel: A Modular Environment for Modeling and Solving Optimization 211 Problems Yves Colombani, Bob Daniel & Susanne Heipcke 12.1 Introduction 211 12.1.1 Solver Modules 212 12.1.2 Other Modules 212 12.1.3 User Modules 213 12.1.4 Contents of this Chapter 213 12.2 The Mosel Language 213 12.2.1 Example Problem 214 12.2.2 Types and Data Structures 215 12.2.3 Initialization of DatalData File Access 217 12.2.4 Language Constructs 218 12.2.4.1 Selections 218 12.2.4.2 Loops 219 12.2.5 Set Operations 219 12.2.6 Subroutines 220 12.3 Mosel Libraries 220 12.4 Mosel Modules 222 12.4.1 Available Modules 223 12.4.2 QP Example with Graphical Output 224 12.4.3 Example of a Solution Algorithm 225 12.5 Writing User Modules 230 12.5.1 Defining a New Subroutine 231 12.5.2 Creating a New Type 233 12.5.2.1 Module Context 234 12.5.2.2 Type Creation and Deletion 235 12.5.2.3 Type Transformation to and from String 235 12.5.2.4 Overloading of Arithmetic Operators 236 12.6 Summary 238 x MODELING LANGUAGES IN MATHEMATICAL OPTIMIZATION 13 The MPL Modeling System 239 Bjarni Kristjansson & Denise Lee 13.1 Maximal Software and Its History 239 13.2 Algebraic Modeling Languages 241 13.2.1 Comparison of Modeling Languages 242 13.2.1.1 Modeling Language 243 13.2.1.2 Multiple Platforms 243 13.2.1.3 Open Design 243 13.2.1.4 Indexing 244 13 .2.1.5 Scalability 244 ] 3.2.].6 Memory Management 244 13.2.1.7 Speed 245 13.2.1.8 Robustness 245 ] 3.2.].9 Deployment 245 13.2.1.1 0 Pricing 245 13.3 MPL Modeling System 245 13.3.] MPL Integrated Model Development Environment 246
Recommended publications
  • How to Use Your Favorite MIP Solver: Modeling, Solving, Cannibalizing
    How to use your favorite MIP Solver: modeling, solving, cannibalizing Andrea Lodi University of Bologna, Italy [email protected] January-February, 2012 @ Universit¨atWien A. Lodi, How to use your favorite MIP Solver Setting • We consider a general Mixed Integer Program in the form: T maxfc x : Ax ≤ b; x ≥ 0; xj 2 Z; 8j 2 Ig (1) where matrix A does not have a special structure. A. Lodi, How to use your favorite MIP Solver 1 Setting • We consider a general Mixed Integer Program in the form: T maxfc x : Ax ≤ b; x ≥ 0; xj 2 Z; 8j 2 Ig (1) where matrix A does not have a special structure. • Thus, the problem is solved through branch-and-bound and the bounds are computed by iteratively solving the LP relaxations through a general-purpose LP solver. A. Lodi, How to use your favorite MIP Solver 1 Setting • We consider a general Mixed Integer Program in the form: T maxfc x : Ax ≤ b; x ≥ 0; xj 2 Z; 8j 2 Ig (1) where matrix A does not have a special structure. • Thus, the problem is solved through branch-and-bound and the bounds are computed by iteratively solving the LP relaxations through a general-purpose LP solver. • The course basically covers the MIP but we will try to discuss when possible how crucial is the LP component (the engine), and how much the whole framework is built on top the capability of effectively solving LPs. • Roughly speaking, using the LP computation as a tool, MIP solvers integrate the branch-and-bound and the cutting plane algorithms through variations of the general branch-and-cut scheme [Padberg & Rinaldi 1987] developed in the context of the Traveling Salesman Problem (TSP).
    [Show full text]
  • Julia, My New Friend for Computing and Optimization? Pierre Haessig, Lilian Besson
    Julia, my new friend for computing and optimization? Pierre Haessig, Lilian Besson To cite this version: Pierre Haessig, Lilian Besson. Julia, my new friend for computing and optimization?. Master. France. 2018. cel-01830248 HAL Id: cel-01830248 https://hal.archives-ouvertes.fr/cel-01830248 Submitted on 4 Jul 2018 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. « Julia, my new computing friend? » | 14 June 2018, IETR@Vannes | By: L. Besson & P. Haessig 1 « Julia, my New frieNd for computiNg aNd optimizatioN? » Intro to the Julia programming language, for MATLAB users Date: 14th of June 2018 Who: Lilian Besson & Pierre Haessig (SCEE & AUT team @ IETR / CentraleSupélec campus Rennes) « Julia, my new computing friend? » | 14 June 2018, IETR@Vannes | By: L. Besson & P. Haessig 2 AgeNda for today [30 miN] 1. What is Julia? [5 miN] 2. ComparisoN with MATLAB [5 miN] 3. Two examples of problems solved Julia [5 miN] 4. LoNger ex. oN optimizatioN with JuMP [13miN] 5. LiNks for more iNformatioN ? [2 miN] « Julia, my new computing friend? » | 14 June 2018, IETR@Vannes | By: L. Besson & P. Haessig 3 1. What is Julia ? Open-source and free programming language (MIT license) Developed since 2012 (creators: MIT researchers) Growing popularity worldwide, in research, data science, finance etc… Multi-platform: Windows, Mac OS X, GNU/Linux..
    [Show full text]
  • Solving Mixed Integer Linear and Nonlinear Problems Using the SCIP Optimization Suite
    Takustraße 7 Konrad-Zuse-Zentrum D-14195 Berlin-Dahlem fur¨ Informationstechnik Berlin Germany TIMO BERTHOLD GERALD GAMRATH AMBROS M. GLEIXNER STEFAN HEINZ THORSTEN KOCH YUJI SHINANO Solving mixed integer linear and nonlinear problems using the SCIP Optimization Suite Supported by the DFG Research Center MATHEON Mathematics for key technologies in Berlin. ZIB-Report 12-27 (July 2012) Herausgegeben vom Konrad-Zuse-Zentrum f¨urInformationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Telefon: 030-84185-0 Telefax: 030-84185-125 e-mail: [email protected] URL: http://www.zib.de ZIB-Report (Print) ISSN 1438-0064 ZIB-Report (Internet) ISSN 2192-7782 Solving mixed integer linear and nonlinear problems using the SCIP Optimization Suite∗ Timo Berthold Gerald Gamrath Ambros M. Gleixner Stefan Heinz Thorsten Koch Yuji Shinano Zuse Institute Berlin, Takustr. 7, 14195 Berlin, Germany, fberthold,gamrath,gleixner,heinz,koch,[email protected] July 31, 2012 Abstract This paper introduces the SCIP Optimization Suite and discusses the ca- pabilities of its three components: the modeling language Zimpl, the linear programming solver SoPlex, and the constraint integer programming frame- work SCIP. We explain how these can be used in concert to model and solve challenging mixed integer linear and nonlinear optimization problems. SCIP is currently one of the fastest non-commercial MIP and MINLP solvers. We demonstrate the usage of Zimpl, SCIP, and SoPlex by selected examples, give an overview of available interfaces, and outline plans for future development. ∗A Japanese translation of this paper will be published in the Proceedings of the 24th RAMP Symposium held at Tohoku University, Miyagi, Japan, 27{28 September 2012, see http://orsj.or.
    [Show full text]
  • [20Pt]Algorithms for Constrained Optimization: [ 5Pt]
    SOL Optimization 1970s 1980s 1990s 2000s 2010s Summary 2020s Algorithms for Constrained Optimization: The Benefits of General-purpose Software Michael Saunders MS&E and ICME, Stanford University California, USA 3rd AI+IoT Business Conference Shenzhen, China, April 25, 2019 Optimization Software 3rd AI+IoT Business Conference, Shenzhen, April 25, 2019 1/39 SOL Optimization 1970s 1980s 1990s 2000s 2010s Summary 2020s SOL Systems Optimization Laboratory George Dantzig, Stanford University, 1974 Inventor of the Simplex Method Father of linear programming Large-scale optimization: Algorithms, software, applications Optimization Software 3rd AI+IoT Business Conference, Shenzhen, April 25, 2019 2/39 SOL Optimization 1970s 1980s 1990s 2000s 2010s Summary 2020s SOL history 1974 Dantzig and Cottle start SOL 1974{78 John Tomlin, LP/MIP expert 1974{2005 Alan Manne, nonlinear economic models 1975{76 MS, MINOS first version 1979{87 Philip Gill, Walter Murray, MS, Margaret Wright (Gang of 4!) 1989{ Gerd Infanger, stochastic optimization 1979{ Walter Murray, MS, many students 2002{ Yinyu Ye, optimization algorithms, especially interior methods This week! UC Berkeley opened George B. Dantzig Auditorium Optimization Software 3rd AI+IoT Business Conference, Shenzhen, April 25, 2019 3/39 SOL Optimization 1970s 1980s 1990s 2000s 2010s Summary 2020s Optimization problems Minimize an objective function subject to constraints: 0 x 1 min '(x) st ` ≤ @ Ax A ≤ u c(x) x variables 0 1 A matrix c1(x) B . C c(x) nonlinear functions @ . A c (x) `; u bounds m Optimization
    [Show full text]
  • Numericaloptimization
    Numerical Optimization Alberto Bemporad http://cse.lab.imtlucca.it/~bemporad/teaching/numopt Academic year 2020-2021 Course objectives Solve complex decision problems by using numerical optimization Application domains: • Finance, management science, economics (portfolio optimization, business analytics, investment plans, resource allocation, logistics, ...) • Engineering (engineering design, process optimization, embedded control, ...) • Artificial intelligence (machine learning, data science, autonomous driving, ...) • Myriads of other applications (transportation, smart grids, water networks, sports scheduling, health-care, oil & gas, space, ...) ©2021 A. Bemporad - Numerical Optimization 2/102 Course objectives What this course is about: • How to formulate a decision problem as a numerical optimization problem? (modeling) • Which numerical algorithm is most appropriate to solve the problem? (algorithms) • What’s the theory behind the algorithm? (theory) ©2021 A. Bemporad - Numerical Optimization 3/102 Course contents • Optimization modeling – Linear models – Convex models • Optimization theory – Optimality conditions, sensitivity analysis – Duality • Optimization algorithms – Basics of numerical linear algebra – Convex programming – Nonlinear programming ©2021 A. Bemporad - Numerical Optimization 4/102 References i ©2021 A. Bemporad - Numerical Optimization 5/102 Other references • Stephen Boyd’s “Convex Optimization” courses at Stanford: http://ee364a.stanford.edu http://ee364b.stanford.edu • Lieven Vandenberghe’s courses at UCLA: http://www.seas.ucla.edu/~vandenbe/ • For more tutorials/books see http://plato.asu.edu/sub/tutorials.html ©2021 A. Bemporad - Numerical Optimization 6/102 Optimization modeling What is optimization? • Optimization = assign values to a set of decision variables so to optimize a certain objective function • Example: Which is the best velocity to minimize fuel consumption ? fuel [ℓ/km] velocity [km/h] 0 30 60 90 120 160 ©2021 A.
    [Show full text]
  • Derivative-Free Optimization: a Review of Algorithms and Comparison of Software Implementations
    J Glob Optim (2013) 56:1247–1293 DOI 10.1007/s10898-012-9951-y Derivative-free optimization: a review of algorithms and comparison of software implementations Luis Miguel Rios · Nikolaos V. Sahinidis Received: 20 December 2011 / Accepted: 23 June 2012 / Published online: 12 July 2012 © Springer Science+Business Media, LLC. 2012 Abstract This paper addresses the solution of bound-constrained optimization problems using algorithms that require only the availability of objective function values but no deriv- ative information. We refer to these algorithms as derivative-free algorithms. Fueled by a growing number of applications in science and engineering, the development of derivative- free optimization algorithms has long been studied, and it has found renewed interest in recent time. Along with many derivative-free algorithms, many software implementations have also appeared. The paper presents a review of derivative-free algorithms, followed by a systematic comparison of 22 related implementations using a test set of 502 problems. The test bed includes convex and nonconvex problems, smooth as well as nonsmooth prob- lems. The algorithms were tested under the same conditions and ranked under several crite- ria, including their ability to find near-global solutions for nonconvex problems, improve a given starting point, and refine a near-optimal solution. A total of 112,448 problem instances were solved. We find that the ability of all these solvers to obtain good solutions dimin- ishes with increasing problem size. For the problems used in this study, TOMLAB/MULTI- MIN, TOMLAB/GLCCLUSTER, MCS and TOMLAB/LGO are better, on average, than other derivative-free solvers in terms of solution quality within 2,500 function evaluations.
    [Show full text]
  • TOMLAB –Unique Features for Optimization in MATLAB
    TOMLAB –Unique Features for Optimization in MATLAB Bad Honnef, Germany October 15, 2004 Kenneth Holmström Tomlab Optimization AB Västerås, Sweden [email protected] ´ Professor in Optimization Department of Mathematics and Physics Mälardalen University, Sweden http://tomlab.biz Outline of the talk • The TOMLAB Optimization Environment – Background and history – Technology available • Optimization in TOMLAB • Tests on customer supplied large-scale optimization examples • Customer cases, embedded solutions, consultant work • Business perspective • Box-bounded global non-convex optimization • Summary http://tomlab.biz Background • MATLAB – a high-level language for mathematical calculations, distributed by MathWorks Inc. • Matlab can be extended by Toolboxes that adds to the features in the software, e.g.: finance, statistics, control, and optimization • Why develop the TOMLAB Optimization Environment? – A uniform approach to optimization didn’t exist in MATLAB – Good optimization solvers were missing – Large-scale optimization was non-existent in MATLAB – Other toolboxes needed robust and fast optimization – Technical advantages from the MATLAB languages • Fast algorithm development and modeling of applied optimization problems • Many built in functions (ODE, linear algebra, …) • GUIdevelopmentfast • Interfaceable with C, Fortran, Java code http://tomlab.biz History of Tomlab • President and founder: Professor Kenneth Holmström • The company founded 1986, Development started 1989 • Two toolboxes NLPLIB och OPERA by 1995 • Integrated format for optimization 1996 • TOMLAB introduced, ISMP97 in Lausanne 1997 • TOMLAB v1.0 distributed for free until summer 1999 • TOMLAB v2.0 first commercial version fall 1999 • TOMLAB starting sales from web site March 2000 • TOMLAB v3.0 expanded with external solvers /SOL spring 2001 • Dash Optimization Ltd’s XpressMP added in TOMLAB fall 2001 • Tomlab Optimization Inc.
    [Show full text]
  • Process Optimization
    Process Optimization Mathematical Programming and Optimization of Multi-Plant Operations and Process Design Ralph W. Pike Director, Minerals Processing Research Institute Horton Professor of Chemical Engineering Louisiana State University Department of Chemical Engineering, Lamar University, April, 10, 2007 Process Optimization • Typical Industrial Problems • Mathematical Programming Software • Mathematical Basis for Optimization • Lagrange Multipliers and the Simplex Algorithm • Generalized Reduced Gradient Algorithm • On-Line Optimization • Mixed Integer Programming and the Branch and Bound Algorithm • Chemical Production Complex Optimization New Results • Using one computer language to write and run a program in another language • Cumulative probability distribution instead of an optimal point using Monte Carlo simulation for a multi-criteria, mixed integer nonlinear programming problem • Global optimization Design vs. Operations • Optimal Design −Uses flowsheet simulators and SQP – Heuristics for a design, a superstructure, an optimal design • Optimal Operations – On-line optimization – Plant optimal scheduling – Corporate supply chain optimization Plant Problem Size Contact Alkylation Ethylene 3,200 TPD 15,000 BPD 200 million lb/yr Units 14 76 ~200 Streams 35 110 ~4,000 Constraints Equality 761 1,579 ~400,000 Inequality 28 50 ~10,000 Variables Measured 43 125 ~300 Unmeasured 732 1,509 ~10,000 Parameters 11 64 ~100 Optimization Programming Languages • GAMS - General Algebraic Modeling System • LINDO - Widely used in business applications
    [Show full text]
  • Specifying “Logical” Conditions in AMPL Optimization Models
    Specifying “Logical” Conditions in AMPL Optimization Models Robert Fourer AMPL Optimization www.ampl.com — 773-336-AMPL INFORMS Annual Meeting Phoenix, Arizona — 14-17 October 2012 Session SA15, Software Demonstrations Robert Fourer, Logical Conditions in AMPL INFORMS Annual Meeting — 14-17 Oct 2012 — Session SA15, Software Demonstrations 1 New and Forthcoming Developments in the AMPL Modeling Language and System Optimization modelers are often stymied by the complications of converting problem logic into algebraic constraints suitable for solvers. The AMPL modeling language thus allows various logical conditions to be described directly. Additionally a new interface to the ILOG CP solver handles logic in a natural way not requiring conventional transformations. Robert Fourer, Logical Conditions in AMPL INFORMS Annual Meeting — 14-17 Oct 2012 — Session SA15, Software Demonstrations 2 AMPL News Free AMPL book chapters AMPL for Courses Extended function library Extended support for “logical” conditions AMPL driver for CPLEX Opt Studio “Concert” C++ interface Support for ILOG CP constraint programming solver Support for “logical” constraints in CPLEX INFORMS Impact Prize to . Originators of AIMMS, AMPL, GAMS, LINDO, MPL Awards presented Sunday 8:30-9:45, Conv Ctr West 101 Doors close 8:45! Robert Fourer, Logical Conditions in AMPL INFORMS Annual Meeting — 14-17 Oct 2012 — Session SA15, Software Demonstrations 3 AMPL Book Chapters now free for download www.ampl.com/BOOK/download.html Bound copies remain available purchase from usual
    [Show full text]
  • Theory and Experimentation
    Acta Astronautica 122 (2016) 114–136 Contents lists available at ScienceDirect Acta Astronautica journal homepage: www.elsevier.com/locate/actaastro Suboptimal LQR-based spacecraft full motion control: Theory and experimentation Leone Guarnaccia b,1, Riccardo Bevilacqua a,n, Stefano P. Pastorelli b,2 a Department of Mechanical and Aerospace Engineering, University of Florida 308 MAE-A building, P.O. box 116250, Gainesville, FL 32611-6250, United States b Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, Torino 10129, Italy article info abstract Article history: This work introduces a real time suboptimal control algorithm for six-degree-of-freedom Received 19 January 2015 spacecraft maneuvering based on a State-Dependent-Algebraic-Riccati-Equation (SDARE) Received in revised form approach and real-time linearization of the equations of motion. The control strategy is 18 November 2015 sub-optimal since the gains of the linear quadratic regulator (LQR) are re-computed at Accepted 18 January 2016 each sample time. The cost function of the proposed controller has been compared with Available online 2 February 2016 the one obtained via a general purpose optimal control software, showing, on average, an increase in control effort of approximately 15%, compensated by real-time implement- ability. Lastly, the paper presents experimental tests on a hardware-in-the-loop six- Keywords: Spacecraft degree-of-freedom spacecraft simulator, designed for testing new guidance, navigation, Optimal control and control algorithms for nano-satellites in a one-g laboratory environment. The tests Linear quadratic regulator show the real-time feasibility of the proposed approach. Six-degree-of-freedom & 2016 The Authors.
    [Show full text]
  • AIMMS 4 Release Notes
    AIMMS 4 Release Notes Visit our web site www.aimms.com for regular updates Contents Contents 2 1 System Requirements 3 1.1 Hardware and operating system requirements ......... 3 1.2 ODBC and OLE DB database connectivity issues ........ 4 1.3 Viewing help files and documentation .............. 5 2 Installation Instructions 7 2.1 Installation instructions ....................... 7 2.2 Solver availability per platform ................... 8 2.3 Aimms licensing ............................ 9 2.3.1 Personal and machine nodelocks ............. 9 2.3.2 Installing an Aimms license ................ 12 2.3.3 Managing Aimms licenses ................. 14 2.3.4 Location of license files ................... 15 2.4 OpenSSL license ............................ 17 3 Project Conversion Instructions 20 3.1 Conversion of projects developed in Aimms 3.13 and before 20 3.2 Conversionof Aimms 3 projects to Aimms 4 ........... 22 3.2.1 Data management changes in Aimms 4.0 ........ 24 4 Getting Support 27 4.1 Reporting a problem ......................... 27 4.2 Known and reported issues ..................... 28 5 Release Notes 29 What’s new in Aimms 4 ........................ 29 Chapter 1 System Requirements This chapter discusses the system requirements necessary to run the various System components of your Win32 Aimms 4 system successfully. When a particular requirements requirement involves the installation of additional system software compo- nents, or an update thereof, the (optional) installation of such components will be part of the Aimms installation procedure. 1.1 Hardware and operating system requirements The following list provides the minimum hardware requirements to run your Hardware Aimms 4 system. requirements 1.6 Ghz or higher x86 or x64 processor XGA display adapter and monitor 1 Gb RAM 1 Gb free disk space Note, however, that performance depends on model size and type and can vary.
    [Show full text]
  • Spatially Explicit Techno-Economic Optimisation Modelling of UK Heating Futures
    Spatially explicit techno-economic optimisation modelling of UK heating futures A thesis submitted in fulfilment of the requirements for the Degree of Doctor of Philosophy UCL Energy Institute University College London Date: 3rd April 2013 Candidate: Francis Li Supervisors: Robert Lowe, Mark Barrett I, Francis Li, confirm that the work presented in this thesis is my own. Where information has been derived from other sources, I confirm that this has been indicated in the thesis. Francis Li, 3rd April 2013 2 For Eugenia 3 Abstract This thesis describes the use of a spatially explicit model to investigate the economies of scale associated with district heating technologies and consequently, their future technical potential when compared against individual building heating. Existing energy system models used for informing UK technology policy do not employ high enough spatial resolutions to map district heating potential at the individual settlement level. At the same time, the major precedent studies on UK district heating potential have not explored future scenarios out to 2050 and have a number of relevant low-carbon heat supply technologies absent from their analyses. This has resulted in cognitive dissonance in UK energy policy whereby district heating is often simultaneously acknowledged as both highly desirable in the near term but ultimately lacking any long term future. The Settlement Energy Demand System Optimiser (SEDSO) builds on key techno-economic studies from the last decade to further investigate this policy challenge. SEDSO can be distinguished from other models used for investigating UK heat decarbonisation by employing a unique combination of extensive spatial detail, technical modelling which captures key cost-related nonlinearities, and a least-cost constrained optimisation approach to technology selection.
    [Show full text]