DOCSLIB.ORG

Comparison of Simulation Tools for Dynamic Models

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Comparison of Simulation Tools for Dynamic Models Comparison of Simulation Tools for Dynamic Models Sveinung M. Sund1, Marianne Plouvier2, and Bernt Lie1 1University of South-Eastern Norway, [email protected], 2IMT Mines Albi, France Abstract els to experimental data, to reuse such models for pro- cess design, stability analysis, control design, estimator Macroscopic models are used extensively in process en- design, etc. The various simulation tools have different gineering, and can often be posed as DAE (Differential support for and capabilities wrt. analysis and synthesis. Algebraic Equation) models. Three generic tools for solv- An important question is then: is it wise to aim at a one- ing such DAEs are compared: OpenModelica, Julia, and language/tool solution, or is it better to use different tools MATLAB. To make the comparison concrete, a simple for simulation and design/synthesis as long as the lan- non-linear process model from the literature was extended guages can interact? by removing simplifying assumptions; the more complex model was posed as DAEs. Some implementation details 1.2 Previous Work of DAE models in OpenModelica, Julia, and MATLAB Specialized process engineering tools are well developed, are given. Selected simulation results are given, with re- e.g., from AspenTech1, Process Systems Enterprise2, etc. sulting execution time. The three tools gave identical sim- More generic tools are also popular, such as Modelica ulation results. The tools are then compared wrt. cost, based tools (OpenModelica3, JModelica.org4, Dymola5, ease of use, documentation, numeric quality, Eco-system etc.), other high level tools (DAE Tools6, APMonitor7), , and possibility for reuse of models/library. Overall, Julia script based tools (MATLAB8, Python9, Julia10, etc.), and appears may appear as the best choice. However, Model- computer algebra systems (Mathematica11, Maple12, etc.). ica is found to be easier to use, so an ideal solution would The possibility to integrate OpenModelica with Python to probably be some tight integration of Modelica with Julia. enable more extended analysis was studied in (Lie et al., Keywords : process modeling, dynamics, DAEs, simula- 2016). tion tool. A simple model of a nonlinear, open loop unstable re- actor is given in (Seborg et al., 2011), and is used in var- 1 Introduction ious control studies (Henson and Seborg, 1997). In a stu- 1.1 Background dent project/exam paper at University College of South- east Norway13, this model has been extended into a DAE Macroscopic models are used extensively in process engi- model by removing some of the simplifying assumptions. neering, and can in general be posed as a set of differen- Both of these models are suitable for testing basic capabil- tial equations stemming from conservation/balance laws ities of simulation tools for handling small-scale models. with the addition of algebraic relations describing thermo- dynamics, transport laws, reaction engineering, etc. The 1.3 Organization of Paper result is a set of Differential Algebraic Equations (DAE) In this paper, we compare the suitability of selected tools with a relatively simple structure. Such DAE models can (OpenModelica, Julia, MATLAB) for solving macro- be transformed/reduced to Ordinary Differential Equa- scopic DAE models in process engineering. The paper tions (ODEs) with some effort. The resulting ODEs are often simpler than the DAEs, thus more efficient solvers 1www.aspentech.com may be available. At the same time, reducing DAEs to 2www.psenterprise.com/products/gproms ODEs also eliminates many variable which may be of in- 3https://openmodelica.org/ terest to study. 4www.jmodelica.org 5www.3ds.com/products-services/catia/ In education, as well as in engineering practice, it is im- products/dymola portant to choose simulation tools for doing experiments 6http://daetools.com on such models. Important aspects are, e.g.: cost, ease 7http://apmonitor.com of use, how well the tool is documented, numeric qual- 8https://mathworks.com ity, Eco-system (packages/functions for plotting, analysis, 9www.anaconda.com 10 random numbers, control packages, access to time series https://julialang.org 11www.wolfram.com database, etc.), possibility of reuse of the model (model 12www.maplesoft.com library, extract information for control design, etc.). 13Course FM1015 Modelling of Dynamic Systems, www.usn.no Dynamic models are used for design and operation of /academics/study-and-subjectplans/#/ systems, and it is therefore of interest to fit such mod- subjects/FM1015_2018H_1 https://doi.org/10.3384/ecp18153177 177 Proceedings of The 59th Conference on Simulation and Modelling (SIMS 59), 26-28 September 2018, Oslo Metropolitan University, Norway is organized as follows. In Section 2, an overview of a suitable DAE model structure is given, with details of the case study model. In Section 3, some general character- istics and specific implementation details of the models in the chosen languages are given. In Section 4, selected simulation results are given, together with some resulting execution times. In Section 5, the findings are discussed, Figure 1. Cooled liquid reactor with reaction aA ! B. before some conclusions are drawn in Section 6. An ap- pendix gives some model details. 2.2 Process Overview 2 Model Overview We consider a liquid reactor of constant volume V, with ˙ 2.1 Mechanistic Model Structure influent volumetric flow rate Vi, influent concentration cA;i of reactant A carried via an inert solvent S,14 and influent Mechanistic models at the macroscopic level can typically temperature Ti, Fig. 1. be described using balance laws of form It is of interest to convert species A into species B through an exothermic reaction dx = f (x;z;u;q) (1) dt aA ! B; (4) where differential variable x is the “balanced” property the products are carried out of the reactor via solvent S. To (amount, momentum, energy), while z is some auxil- keep control of reactor temperature T, heat rate Q˙ is added iary algebraic variable, u is some input variable and q by flowing a liquid at temperature Tc through the tube side is some model constant/parameter. Such models need to of a coil/heat exchanger. With a high flow rate of the cool- be complemented with transport/thermodynamic/reaction ing liquid, Tc is constant through the heat exchanger, and laws (constitutive laws) of form the heat rate can be described as 0 = g(x;z;u;q): (2) Q˙ = U A(Tc − T): (5) Together, Eqs. 1 and 2 form a set of differential-algebraic where U A is a parameter. If Tc < T, Q˙ < 0 and the reactor equations, DAE. The set of differential variables and alge- is cooled. braic variables is referred to as the model descriptor, (x;z). The rate of reaction r is given as Such DAE models are conceptually simple to formulate, but may contain relatively many descriptor variables and a r = kcA (6) many parameters q. DAE models can be manipulated into sets of ordinary where a is the reaction order and k is given by the Arrhe- differential equations, ODE, of form nius expression ( ) dx˜ ( ) E=R = f˜ x˜;u;q˜ (3) k = k exp − . (7) dt 0 T wherex ˜ is a state of the system. For DAE models of index The operation of the reactor is influenced by inputs V˙ , 1, dimx ≥ dimx ˜, while for DAE models of index larger i c ; , T and Tc, and it is of interest to study how these in- than 1, it is possible that dimx ≤ dimx ˜. Normally dimq ≥ A i i fluence the outputs cA, cB and T. ( ) dimq˜. Although the case study has inputs u = V˙i;cA;i;Ti;Tc Tools for analyzing ODE models (solvers, stability con- and outputs y = (cA;cB;T), in a control problem one may cepts, etc.) are better developed than tools for DAE mod- choose to control the temperature y = T by manipulat- els. ODE models are smaller and faster to solve, but man- (ing the input) u = Tc. In that case, the additional inputs ual model simplification may introduce (model) errors. V˙i;cA;i;Ti may be considered disturbances. DAE models, on the other hand, hold more information For practical purposes, the cooling liquid temperature (the algebraic variables), and are simpler to formulate, so ◦ 15 Tc 2 [4;90] C or so. In reality, Tc is not directly con- there are important advantages in keeping models in DAE trollable. Instead, an influent “cooling” temperature T form. c;i could be varied while temperature Tc in the heat exchanger It is of interest to compare how easy it is to formulate is varied via manipulating the cooling liquid flow rate. DAE models and solve them in popular modeling/simu- lation tools. For such a comparison, it is convenient to 14Inert implies non-reacting. introduce a case study. 15If vaporized, higher temperatures may be achieved. https://doi.org/10.3384/ecp18153177 178 Proceedings of The 59th Conference on Simulation and Modelling (SIMS 59), 26-28 September 2018, Oslo Metropolitan University, Norway ˙ • ≡ 2.3 Balance Equations VB;i 0 (no feed of species B). Equations 12, 13, and 17 For the given process, the following balance equations are are based on the assumption of ideal solution, and are not relevant, generally valid. Equations 14–16, on the other hand, are valid by definition. dmS In addition to the ideal solution assumption, we need = m˙ − m˙ (8) dt S;i S;e a model of liquid specific enthalpy. For an ideal (in- dn compressible) liquid with temperature independent spe- A = n˙ − n˙ + n˙ (9) dt A;i A;e A;g cific heat capacity: dnB • ◦ • ◦ = n˙ − n˙ ; + n˙ ; (10) ˆ ˆ − dt B;i B e B g HS = HS + cˆp;S (T T ) (18) ( ) dU ◦ ◦ = H˙i − H˙e + Q˙ (11) 16 Hˆ ;T dt and similar for species A and B. Here, S are specific enthalpy at standard state and a chosen reference 6 dx ◦ wherex ˙ indicates flow of variable x (in general,x ˙ = dt ), m temperature T .
Load more

Recommended publications
  • Click to Edit Master Title Style
    Click to edit Master title style MINLP with Combined Interior Point and Active Set Methods Jose L. Mojica Adam D. Lewis John D. Hedengren Brigham Young University INFORM 2013, Minneapolis, MN Presentation Overview NLP Benchmarking Hock-Schittkowski Dynamic optimization Biological models Combining Interior Point and Active Set MINLP Benchmarking MacMINLP MINLP Model Predictive Control Chiller Thermal Energy Storage Unmanned Aerial Systems Future Developments Oct 9, 2013 APMonitor.com APOPT.com Brigham Young University Overview of Benchmark Testing NLP Benchmark Testing 1 1 2 3 3 min J (x, y,u) APOPT , BPOPT , IPOPT , SNOPT , MINOS x Problem characteristics: s.t. 0 f , x, y,u t Hock Schittkowski, Dynamic Opt, SBML 0 g(x, y,u) Nonlinear Programming (NLP) Differential Algebraic Equations (DAEs) 0 h(x, y,u) n m APMonitor Modeling Language x, y u MINLP Benchmark Testing min J (x, y,u, z) 1 1 2 APOPT , BPOPT , BONMIN x s.t. 0 f , x, y,u, z Problem characteristics: t MacMINLP, Industrial Test Set 0 g(x, y,u, z) Mixed Integer Nonlinear Programming (MINLP) 0 h(x, y,u, z) Mixed Integer Differential Algebraic Equations (MIDAEs) x, y n u m z m APMonitor & AMPL Modeling Language 1–APS, LLC 2–EPL, 3–SBS, Inc. Oct 9, 2013 APMonitor.com APOPT.com Brigham Young University NLP Benchmark – Summary (494) 100 90 80 APOPT+BPOPT APOPT 70 1.0 BPOPT 1.0 60 IPOPT 3.10 IPOPT 50 2.3 SNOPT Percentage (%) 6.1 40 Benchmark Results MINOS 494 Problems 5.5 30 20 10 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 Not worse than 2 times slower than
    [Show full text]
  • Click to Edit Master Title Style
    Click to edit Master title style APMonitor Modeling Language John Hedengren Brigham Young University Advanced Process Solutions, LLC http://apmonitor.com Overview of APM Software as a service accessible through: MATLAB, Python, Web-browser interface Linux / Windows / Mac OS / Android platforms Solvers 1 1 2 3 3 APOPT , BPOPT , IPOPT , SNOPT , MINOS Problem characteristics: min J (x, y,u, z) Large-scale x s.t. 0 f , x, y,u, z Nonlinear Programming (NLP) t Mixed Integer NLP (MINLP) 0 g(x, y,u, z) Multi-objective 0 h(x, y,u, z) n m m Real-time systems x, y u z Differential Algebraic Equations (DAEs) 1 – APS, LLC 2 – EPL 3 – SBS, Inc. Oct 14, 2012 APMonitor.com Advanced Process Solutions, LLC Overview of APM Vector / matrix algebra with set notation Automatic Differentiation st nd Exact 1 and 2 Derivatives Large-scale, sparse systems of equations Object-oriented access Thermo-physical properties Database of preprogrammed models Parallel processing Optimization with uncertain parameters Custom solver or model connections Oct 14, 2012 APMonitor.com Advanced Process Solutions, LLC Unique Features of APM Initialization with nonlinear presolve minJ(x, y,u) x s.t. 0 f ,x, y,u min J (x, y,u) t 0 g(x, y,u) 0h(x, y,u) x minJ(x, y,u) x s.t. 0 f ,x, y,u s.t. 0 f , x, y,u t 0 g(x, y,u) t 0 h(x, y,u) minJ(x, y,u) x s.t. 0 f ,x, y,u t 0g(x, y,u) 0h(x, y,u) 0 g(x, y,u) minJ(x, y,u) x s.t.
    [Show full text]
  • Treball (1.484Mb)
    Treball Final de Màster MÀSTER EN ENGINYERIA INFORMÀTICA Escola Politècnica Superior Universitat de Lleida Mòdul d’Optimització per a Recursos del Transport Adrià Vall-llaura Salas Tutors: Antonio Llubes, Josep Lluís Lérida Data: Juny 2017 Pròleg Aquest projecte s’ha desenvolupat per donar solució a un problema de l’ordre del dia d’una empresa de transports. Es basa en el disseny i implementació d’un model matemàtic que ha de permetre optimitzar i automatitzar el sistema de planificació de viatges de l’empresa. Per tal de poder implementar l’algoritme s’han hagut de crear diversos mòduls que extreuen les dades del sistema ERP, les tracten, les envien a un servei web (REST) i aquest retorna un emparellament òptim entre els vehicles de l’empresa i les ordres dels clients. La primera fase del projecte, la teòrica, ha estat llarga en comparació amb les altres. En aquesta fase s’ha estudiat l’estat de l’art en la matèria i s’han repassat molts dels models més importants relacionats amb el transport per comprendre’n les seves particularitats. Amb els conceptes ben estudiats, s’ha procedit a desenvolupar un nou model matemàtic adaptat a les necessitats de la lògica de negoci de l’empresa de transports objecte d’aquest treball. Posteriorment s’ha passat a la fase d’implementació dels mòduls. En aquesta fase m’he trobat amb diferents limitacions tecnològiques degudes a l’antiguitat de l’ERP i a l’ús del sistema operatiu Windows. També han sorgit diferents problemes de rendiment que m’han fet redissenyar l’extracció de dades de l’ERP, el càlcul de distàncies i el mòdul d’optimització.
    [Show full text]
  • GEKKO Documentation Release 1.0.1
    GEKKO Documentation Release 1.0.1 Logan Beal, John Hedengren Aug 31, 2021 Contents 1 Overview 1 2 Installation 3 3 Project Support 5 4 Citing GEKKO 7 5 Contents 9 6 Overview of GEKKO 89 Index 91 i ii CHAPTER 1 Overview GEKKO is a Python package for machine learning and optimization of mixed-integer and differential algebraic equa- tions. It is coupled with large-scale solvers for linear, quadratic, nonlinear, and mixed integer programming (LP, QP, NLP, MILP, MINLP). Modes of operation include parameter regression, data reconciliation, real-time optimization, dynamic simulation, and nonlinear predictive control. GEKKO is an object-oriented Python library to facilitate local execution of APMonitor. More of the backend details are available at What does GEKKO do? and in the GEKKO Journal Article. Example applications are available to get started with GEKKO. 1 GEKKO Documentation, Release 1.0.1 2 Chapter 1. Overview CHAPTER 2 Installation A pip package is available: pip install gekko Use the —-user option to install if there is a permission error because Python is installed for all users and the account lacks administrative priviledge. The most recent version is 0.2. You can upgrade from the command line with the upgrade flag: pip install--upgrade gekko Another method is to install in a Jupyter notebook with !pip install gekko or with Python code, although this is not the preferred method: try: from pip import main as pipmain except: from pip._internal import main as pipmain pipmain(['install','gekko']) 3 GEKKO Documentation, Release 1.0.1 4 Chapter 2. Installation CHAPTER 3 Project Support There are GEKKO tutorials and documentation in: • GitHub Repository (examples folder) • Dynamic Optimization Course • APMonitor Documentation • GEKKO Documentation • 18 Example Applications with Videos For project specific help, search in the GEKKO topic tags on StackOverflow.
    [Show full text]
  • Component-Oriented Acausal Modeling of the Dynamical Systems in Python Language on the Example of the Model of the Sucker Rod String
    Component-oriented acausal modeling of the dynamical systems in Python language on the example of the model of the sucker rod string Volodymyr B. Kopei, Oleh R. Onysko and Vitalii G. Panchuk Department of Computerized Mechanical Engineering, Ivano-Frankivsk National Technical University of Oil and Gas, Ivano-Frankivsk, Ukraine ABSTRACT Typically, component-oriented acausal hybrid modeling of complex dynamic systems is implemented by specialized modeling languages. A well-known example is the Modelica language. The specialized nature, complexity of implementation and learning of such languages somewhat limits their development and wide use by developers who know only general-purpose languages. The paper suggests the principle of developing simple to understand and modify Modelica-like system based on the general-purpose programming language Python. The principle consists in: (1) Python classes are used to describe components and their systems, (2) declarative symbolic tools SymPy are used to describe components behavior by difference or differential equations, (3) the solution procedure uses a function initially created using the SymPy lambdify function and computes unknown values in the current step using known values from the previous step, (4) Python imperative constructs are used for simple events handling, (5) external solvers of differential-algebraic equations can optionally be applied via the Assimulo interface, (6) SymPy package allows to arbitrarily manipulate model equations, generate code and solve some equations symbolically. The basic set of mechanical components (1D translational “mass”, “spring-damper” and “force”) is developed. The models of a sucker rods string are developed and simulated using these components. The comparison of results of the sucker rod string simulations with practical dynamometer cards and Submitted 22 March 2019 Accepted 24 September 2019 Modelica results verify the adequacy of the models.
    [Show full text]
  • Nonlinear Modeling, Estimation and Predictive Control in Apmonitor
    Brigham Young University BYU ScholarsArchive Faculty Publications 2014-11-05 Nonlinear Modeling, Estimation and Predictive Control in APMonitor John Hedengren Brigham Young University, [email protected] Reza Asgharzadeh Shishavan Brigham Young University Kody M. Powell University of Utah Thomas F. Edgar University of Texas at Austin Follow this and additional works at: https://scholarsarchive.byu.edu/facpub Part of the Chemical Engineering Commons BYU ScholarsArchive Citation Hedengren, John; Asgharzadeh Shishavan, Reza; Powell, Kody M.; and Edgar, Thomas F., "Nonlinear Modeling, Estimation and Predictive Control in APMonitor" (2014). Faculty Publications. 1667. https://scholarsarchive.byu.edu/facpub/1667 This Peer-Reviewed Article is brought to you for free and open access by BYU ScholarsArchive. It has been accepted for inclusion in Faculty Publications by an authorized administrator of BYU ScholarsArchive. For more information, please contact [email protected], [email protected]. Computers and Chemical Engineering 70 (2014) 133–148 Contents lists available at ScienceDirect Computers and Chemical Engineering j ournal homepage: www.elsevier.com/locate/compchemeng Nonlinear modeling, estimation and predictive control in APMonitor a,∗ a b b John D. Hedengren , Reza Asgharzadeh Shishavan , Kody M. Powell , Thomas F. Edgar a Department of Chemical Engineering, Brigham Young University, Provo, UT 84602, United States b The University of Texas at Austin, Austin, TX 78712, United States a r t i c l e i n f o a b s t r a c t Article history: This paper describes nonlinear methods in model building, dynamic data reconciliation, and dynamic Received 10 August 2013 optimization that are inspired by researchers and motivated by industrial applications.
    [Show full text]
  • Datascientist Manual
    DATASCIENTIST MANUAL . 2 « Approcherait le comportement de la réalité, celui qui aimerait s’epanouir dans l’holistique, l’intégratif et le multiniveaux, l’énactif, l’incarné et le situé, le probabiliste et le non linéaire, pris à la fois dans l’empirique, le théorique, le logique et le philosophique. » 4 (* = not yet mastered) THEORY OF La théorie des probabilités en mathématiques est l'étude des phénomènes caractérisés PROBABILITY par le hasard et l'incertitude. Elle consistue le socle des statistiques appliqué. Rubriques ↓ Back to top ↑_ Notations Formalisme de Kolmogorov Opération sur les ensembles Probabilités conditionnelles Espérences conditionnelles Densités & Fonctions de répartition Variables aleatoires Vecteurs aleatoires Lois de probabilités Convergences et théorèmes limites Divergences et dissimilarités entre les distributions Théorie générale de la mesure & Intégration ------------------------------------------------------------------------------------------------------------------------------------------ 6 Notations [pdf*] Formalisme de Kolmogorov Phé nomé né alé atoiré Expé riéncé alé atoiré L’univérs Ω Ré alisation éléméntairé (ω ∈ Ω) Evé némént Variablé aléatoiré → fonction dé l’univér Ω Opération sur les ensembles : Union / Intersection / Complémentaire …. Loi dé Augustus dé Morgan ? Independance & Probabilités conditionnelles (opération sur des ensembles) Espérences conditionnelles 8 Variables aleatoires (discret vs reel) … Vecteurs aleatoires : Multiplét dé variablés alé atoiré (discret vs reel) … Loi marginalés Loi
    [Show full text]
  • Optimization of Battery Storage Profitability with Wind Energy
    Optimization of Battery Storage Profitability with Wind Energy Joseph Bloxham, Baichuan Liu, John Hedengren, Logan Beal Date 4/20/2018 Table of Contents Table of Contents ................................................................................................................ ii List of Figures ...................................................................................................................... iii List of Tables ....................................................................................................................... iv Optimization of Battery Storage Profitability with Wind Energy ....................................... 1 Abstract .............................................................................................................................. 1 Introduction and Background............................................................................................. 1 Modeling Framework ......................................................................................................... 2 Software ......................................................................................................................... 2 Model Overview ............................................................................................................. 2 Model Assumptions ....................................................................................................... 3 Input Data .....................................................................................................................
    [Show full text]
  • Numerical Optimization of Industrial Processes
    Czech Technical University in Prague Faculty of Electrical Engineering Department of Control Engineering Bachelor’s thesis Numerical optimization of industrial processes Lenka Caletková Supervisor: Ing. Jan Šulc Study Programme: Cybernetics and Robotics Field of Study: Systems and control Date: May 2013 iv Acknowledgements I would like to show my greatest appreciation to my supervisor Ing. Jan Šulc, who was always willing to help me with this project - even in his extra time - during the last two semesters. Without his encouragement and guidance this project would not have materialized. Many thanks belong to all my friends, who supported me in my tough times. Last but not least, I want to thank my supportive sister Jana and my great parents, who enabled me to study at the CTU, always stood by my side and gave me all their love. Abstract This thesis handles the conception of numerical optimizations and presents Python libraries, which are able to solve optimization problems. It compares their usability, stability and computing possibilities. All these aspects were tested on two specific industrial processes - control of the tunnel ventilation and heating systems control in buildings. A functional code in Python programming language for performance optimization of tunnel ventilation represents the outcome of this thesis. Abstrakt Tato bakalářská práce se věnuje pojmu numerické optimalizace a prezentuje knihovny pro- gramovacího jazyka Python, které je možné použít pro řešení optimalizačních úloh. Zároveň porovnává jejich použitelnost, stabilitu a výpočetní schopnosti. Všechny tyto aspekty byly testovány na dvou konkrétních průmyslových procesech - řízení ventilace v tunelech a řízení tepelných systémů v budovách. Výstup této práce představuje funkční kód v jazyce Python, který má za úkol optimalizovat výkon při ventilaci tunelu.
    [Show full text]
  • Software Implementation for Optimization of Production Planning Within Refineries
    Vol.VII(LXX) 1 - 10 Economic Insights – Trends and Challenges No. 1/2018 Software Implementation for Optimization of Production Planning within Refineries Cătălin Popescu Faculty of Economic Sciences, Petroleum-Gas University of Ploieşti, Bd. Bucureşti 39, 100680, Ploieşti, Romania e-mail: [email protected] Abstract The paper will analyse the advantages of using the RPMS program in order to optimize the production plan and to decide on the structure of raw materials and final products within refineries. Also, the paper will include some real-time refining activities to be analysed with the aim of optimizing production. Key words: optimization; RPM; production planning; decision point JEL Classification: C44; C61; C63; L23; L86 Introduction The most popular manufacturing optimization software used today by most refineries are RPMS (Honeywell License), PIMS (Aspen License), and GRTMPS (Harperly Systems License), all based on the LP-Linear Programming principles. These allow modelling of the refinery by using different tables and aim to find the optimal solution taking into account a number of limitations introduced in LP, either linear or non-linear. The constant change in demand-to-offer ratio on the global energy market, as well as the diversity of raw materials and final products on the market require significant economic and financial assessment capacities in the shortest possible time. The development of information systems within companies has allowed shorten the time to assess opportunities and risks in making decisions in order to conclude contracts for the purchase and sale of petroleum products, by taking into account all the expenditures related to the production and transport of goods.
    [Show full text]
  • Edit Master Title Style
    Click to edit Master title style APOPT and BPOPT: New NLP and MINLP Solvers for Differential Algebraic Systems with Benchmark Testing John D. Hedengren Jose L. Mojica Brigham Young University Kody M. Powell Wesley Cole Thomas F. Edgar The University of Texas at Austin Presentation Overview Overview of Benchmarking Testing. NLP Benchmarking Hock-Schittkowski Dynamic optimization Biological models MINLP Benchmarking MacMINLP Case studies Gravity Drained tank (MIDAE control problem) Energy Storage for smart grid Future Developments Oct 17, 2012 APMonitor.com | APOPT.com BYU – UT Austin Overview of Benchmark Testing NLP Benchmark Testing 1 1 2 3 3 min J (x, y,u) APOPT , BPOPT , IPOPT , SNOPT , MINOS x Problem characteristics: s.t. 0 f , x, y,u t Hock Schittkowski, Dynamic Opt, SBML 0 g(x, y,u) Nonlinear Programming (NLP) Differential Algebraic Equations (DAEs) 0 h(x, y,u) n m APMonitor Modeling Language x, y u MINLP Benchmark Testing min J (x, y,u, z) 1 1 2 4 APOPT , BPOPT , BONMIN , DICOPT x s.t. 0 f , x, y,u, z Problem characteristics: t MacMINLP, Industrial Test Set 0 g(x, y,u, z) Mixed Integer Nonlinear Programming (MINLP) 0 h(x, y,u, z) Mixed Integer Differential Algebraic Equations (MIDAEs) x, y n u m z m APMonitor & AMPL Modeling Language 1–APS, LLC 2–EPL, 3–SBS, Inc., 4-CMU Oct 17, 2012 APMonitor.com | APOPT.com BYU – UT Austin NLP Benchmark – Hock-Schittkowski (116) 100 90 80 70 60 APOPT+BPOPT APOPT 50 1.0 BPOPT 40 1.0 Percentage (%) IPOPT 3.10 30 IPOPT 2.3 SNOPT 20 6.1 MINOS 5.5 10 0 0 0.5
    [Show full text]
  • Linear Programming 1 Linear Programming
    Linear programming 1 Linear programming Linear programming (LP, or linear optimization) is a method to achieve the best outcome (such as maximum profit or lowest cost) in a mathematical model whose requirements are represented by linear relationships. Linear programming is a special case of mathematical programming (mathematical optimization). More formally, linear programming is a technique for the optimization of a linear objective function, subject to linear equality and linear inequality constraints. Its feasible region is a convex polyhedron, which is a set defined as the intersection of finitely many half spaces, each of which is defined by a linear inequality. Its objective function is a real-valued affine function defined on this polyhedron. A linear programming algorithm finds a point in the polyhedron where this A pictorial representation of a simple linear program with two variables and six inequalities. function has the smallest (or largest) value if such a point exists. The set of feasible solutions is depicted in light Linear programs are problems that can be expressed in canonical form: red and forms a polygon, a 2-dimensional polytope. The linear cost function is represented by the red line and the arrow: The red line is a level set of the cost function, and the arrow indicates the direction in which we are optimizing. A closed feasible region of a problem with three variables is a convex polyhedron. The surfaces giving a fixed value of the objective function are planes (not shown). The linear programming problem is to find a point on the polyhedron that is on the plane with the highest possible value.
    [Show full text]