How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE
|
|
- Cecilia Powell
- 5 years ago
- Views:
Transcription
1 How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE Tamal Das PhD Candidate IKP, NTNU 24th May 2017
2 Outline 1 Differential algebraic equations (DAE) 2 Extended Kalman filter (EKF) for DAE 3 Moving horizon estimation (MHE) Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 1 / 24
3 Differential algebraic equations (DAE) Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 2 / 24
4 Types of DAEs: Fully implicit and Semi-explicit Fully implicit DAE Notation x : Differential variables z : Algebraic variables y = [ x T z T ] F ( y, dy ) = 0 dt Semi-explicit DAE y (0) = y 0 dx = f (x, z) dt 0 = g (x, z) Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 3 / 24
5 DAEs in Chemical Engineering In ChemEng, mathematical models consist of the following: Dynamic conservation laws [Differential equations] Mass balance Energy balance Static conservation laws [Algebraic equations] Constitutive equations [Algebraic equations] Equation of state Pressure drop equation Heat transfer equation Physical and operating constraints [Algebraic equations] Desired operation Design constraints Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 4 / 24
6 Index of a DAE Index 1-DAE g z In ChemEng, commonly encountered DAEs are semi-explicit dx dt = f (x, z) ; 0 = g (x, z) Index of DAEs determines how hard it is to solve them is non-singular ) det 0 g z ( g z is full rank matrix Higher index (2,3,... -) DAEs g z is singular ) det = 0 g z ( g z is not full rank matrix Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 5 / 24
7 Examples of DAE indices F, c 0 dc dt = (c 0 c) τ Index 1-DAE c 0 = γ (t) is specified; g = c 0 γ (t) g z = g c 0 = 1 (Full rank) CSTR with volume V Residence time τ = V F Differential variable (x) : c Algebraic variable (z) : c 0 F, c Higher index (2,3,... -) DAEs c = γ (t) is specified; g = c γ (t) g z = g c 0 = 0 (not full rank) Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 6 / 24
8 How to determine DAE index? Index = Number of differentiations of g (x, z) w.r.t time necessary to get a differential equation for each algebraic variable Index 1-DAE c 0 = γ (t) is specified Derivative 1: dc 0 dt Index = 1 = γ (t) Higher index (2,3,... -) DAEs c = γ (t) is specified Derivative 1: dc Rearrange: c 0 c τ dt = γ (t) = γ (t) Derivative 2: dc 0 dt dc dt = τγ Rearrange: dc 0 dt = γ + τγ Index = 2 Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 7 / 24
9 Why is Index of a DAE important? Index 1-DAEs are easy whereas higher index-daes are hard Higher index DAEs are hard due to initialization of variables (among others) In ChemEng, DAEs of both kinds are encountered dx = f (x, z) ; 0 = g (x, z) dt ( ) g det 0 (Index 1-DAE), or z ( ) g det = 0 (Higher index DAE) z Where does one encounter higher index DAEs in ChemEng Making simplifications in process models Converting simulation problems into control problems Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 8 / 24
10 Initialization issues of DAEs Index 1-DAE (Easy) Initialize differential variables x (0) Solve g (x (0), z (0)) = 0 for algebraic variables z (0) Higher index DAEs (Hard) x (0) and z (0) are interdependent x (0) and z (0) cannot be specified independently If specified independently, leads to inconsistent initial conditions DAE solver will crash Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 9 / 24
11 DAE example: Two component dynamic flash calculation F, z 1, z 2 P, T M V M L V, y 1, y 2 P out Volume V L, x 1, x 2 M T = F V L M 1 = Fz 1 Vy 1 Lx 1 M 2 = Fz 2 Vy 2 Lx 2 M T = M L + M V M 1 = M L x 1 + M V y 1 M 2 = M L x 2 + M V y 2 x 1 + x 2 = y 1 + y 2 7 Conservation equations Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 10 / 24
12 DAE example: Two component dynamic flash calculation F, z 1, z 2 P, T M V M L V, y 1, y 2 P out Volume V L, x 1, x 2 y 1 = K 1 x 1 y 2 = K 2 x 2 K 1 = f K 1 (T, P) K 2 = f K 2 (T, P) ρ L = f L (x 1 ) ρ V = f V (T, P) L = h L (M L, ρ L ) V = h V (P P out ) V = M L ρ L + M V ρ V } 8 Constitutive equations Design Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 11 / 24
13 DAE example: Two component dynamic flash calculation P, T V, y 1, y 2 P out 3 differential equations 3 differential variables M V 13 algebraic equations F, z 1, z 2 M L Volume V 13 algebraic variables We can find an output set for all the algebraic and differential variables L, x 1, x 2 Index 1 Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 12 / 24
14 DAE example: Modified two component dynamic flash F, z 1, z 2 P, T M L V, y 1, y 2 M L = F V L M 1 = Fz 1 Vy 1 Lx 1 M 2 = Fz 2 Vy 2 Lx 2 M 1 = M L x 1 M 2 = M L x 2 x 1 + x 2 = y 1 + y 2 6 Conservation equations L, x 1, x 2 Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 13 / 24
15 DAE example: Modified two component dynamic flash F, z 1, z 2 P, T M L V, y 1, y 2 y 1 = K 1 x 1 y 2 = K 2 x 2 K 1 = f K 1 (T, P) K 2 = f K 2 (T, P) L = h L (M L, ρ L ) ρ L = f L (x 1 ) 6 Constitutive equations L, x 1, x 2 Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 14 / 24
16 DAE example: Modified two component dynamic flash F, z 1, z 2 V, y 1, y 2 P, T 3 differential equations M L L, x 1, x 2 3 differential variables 9 algebraic equations 9 algebraic variables We cannot find an output set for algebraic variable V. We need to differentiate twice to get there Index 2 Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 15 / 24
17 Solving higher index DAE problems Reduce index to 1 or 0 Find consistent initial conditions for all the variables of the new system Solve the new system as a Index 1 DAE system or 0 index system (ODE system) The procedure 1 is explained on the MATLAB website with added support functions: Index / order reduction ( reducedaeindex ) Finding consistent initial guesses ( decic ) DAE solvers: ode15i, ode15s, or ode23t 1 Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 16 / 24
18 Solving index 1-DAE: Robertson DAE in Simulink A B + C B C dc A = 0.04C A C B C C dt dc B = 0.04C A 10 4 C B C C CB 2 dt 0 = 1 (C A + C B + C C ) Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 17 / 24
19 Solving index 1-DAE problem in Simulink: Implementation Interpreted MATLAB Fcn derivatives 1 s Integrator z Solve f(z) = 0 f(z) Algebraic Constraint Interpreted MATLAB Fcn Algebraics Differential states Algebraic states Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 18 / 24
20 Solving index 1-DAE problem in Simulink: Result Robertson DAE problem with a Conservation Law, solved by ODE15S C A C C Concentration [-] C B t (sec) Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 19 / 24
21 Extended Kalman filter (EKF) for DAE Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 20 / 24
22 Continuous-discrete EKF algorithm for DAE dx = f (x, z) + w dt x aug = 0 = g (x, z) ; y = h (x, z) + v [ ] d ˆx dt = f (ˆx) ; L = I ( ) 1 ( ) g g z x ( ) ( ) dp f aug f aug T dt = x aug P + P x aug + LQL T ˆx aug,+ k K k = P k ( h x aug ) ( T ( h R k + = ˆx aug, k + K k ( y k h k ( P + k = ( I K k ( h x aug )) P k x aug ˆx aug, k )) ) P k ( I K k ( [ ] [ ] x ; f aug f = z g ( h ) ) T 1 x aug h x aug )) T + K k R k K T k Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 21 / 24
23 Moving horizon estimation (MHE) Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 22 / 24
24 MHE is a dynamic optimization problem At each sample time, a set of past few measurements are used to provide a state estimate This problem is discretized to a nonlinear program (NLP) [for nonlinear estimation] NLP is solved using IPOPT solver in CasADi CasADi is interfaced with Simulink using system block to provide online state estimation Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 23 / 24
25 Conclusions DAE systems are hard: particularly higher index DAEs Two approaches: simultaneous (Mass matrix solver) sequential (Simulink with Algebraic constraint block) Estimation for DAE systems: CD-DAE-EKF in Simulink MPC in Simulink using system block and CasADi MHE in Simulink using system block and CasADi For simulation files and this presentation, follow 2 Thank you for your attention 2 folk.ntnu.no/tamald Tamal Das How to simulate in Simulink: DAE, DAE-EKF, MPC & MHE 24 / 24
10.34: Numerical Methods Applied to Chemical Engineering. Lecture 19: Differential Algebraic Equations
10.34: Numerical Methods Applied to Chemical Engineering Lecture 19: Differential Algebraic Equations 1 Recap Differential algebraic equations Semi-explicit Fully implicit Simulation via backward difference
More informationDifferential Equation Types. Moritz Diehl
Differential Equation Types Moritz Diehl Overview Ordinary Differential Equations (ODE) Differential Algebraic Equations (DAE) Partial Differential Equations (PDE) Delay Differential Equations (DDE) Ordinary
More informationChemical Engineering Applications in Scilab
Chemical Engineering Applications in Scilab Prashant Dave Indian Institute of Technology Bombay (IIT Bombay) Introduction In Chemical Engineering the type of problems that occur are Modeling and Simulation
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 informationEfficient Moving Horizon Estimation of DAE Systems. John D. Hedengren Thomas F. Edgar The University of Texas at Austin TWMCC Feb 7, 2005
Efficient Moving Horizon Estimation of DAE Systems John D. Hedengren Thomas F. Edgar The University of Teas at Austin TWMCC Feb 7, 25 Outline Introduction: Moving Horizon Estimation (MHE) Eplicit Solution
More information1.2 Derivation. d p f = d p f(x(p)) = x fd p x (= f x x p ). (1) Second, g x x p + g p = 0. d p f = f x g 1. The expression f x gx
PDE-constrained optimization and the adjoint method Andrew M. Bradley November 16, 21 PDE-constrained optimization and the adjoint method for solving these and related problems appear in a wide range of
More informationMath 225 Differential Equations Notes Chapter 1
Math 225 Differential Equations Notes Chapter 1 Michael Muscedere September 9, 2004 1 Introduction 1.1 Background In science and engineering models are used to describe physical phenomena. Often these
More information10.34 Numerical Methods Applied to Chemical Engineering. Quiz 2
10.34 Numerical Methods Applied to Chemical Engineering Quiz 2 This quiz consists of three problems worth 35, 35, and 30 points respectively. There are 4 pages in this quiz (including this cover page).
More informationLecture 2. Introduction to Differential Equations. Roman Kitsela. October 1, Roman Kitsela Lecture 2 October 1, / 25
Lecture 2 Introduction to Differential Equations Roman Kitsela October 1, 2018 Roman Kitsela Lecture 2 October 1, 2018 1 / 25 Quick announcements URL for the class website: http://www.math.ucsd.edu/~rkitsela/20d/
More informationNonlinear Modeling, Estimation and Predictive Control in APMonitor
Nonlinear Modeling, Estimation and Predictive Control in APMonitor John D. Hedengren, Reza Asgharzadeh Shishavan Department of Chemical Engineering, Brigham Young University, Provo, UT 84602 Kody M. Powell,
More informationSensor Fusion, 2014 Lecture 1: 1 Lectures
Sensor Fusion, 2014 Lecture 1: 1 Lectures Lecture Content 1 Course overview. Estimation theory for linear models. 2 Estimation theory for nonlinear models 3 Sensor networks and detection theory 4 Nonlinear
More informationThe Kalman Filter ImPr Talk
The Kalman Filter ImPr Talk Ged Ridgway Centre for Medical Image Computing November, 2006 Outline What is the Kalman Filter? State Space Models Kalman Filter Overview Bayesian Updating of Estimates Kalman
More informationMath 322. Spring 2015 Review Problems for Midterm 2
Linear Algebra: Topic: Linear Independence of vectors. Question. Math 3. Spring Review Problems for Midterm Explain why if A is not square, then either the row vectors or the column vectors of A are linearly
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 informationSolutions of separable equations often are in an implicit form; that is, not solved for the dependent variable.
Section 2.2 Separable DEs Key Terms: Separable DEs; the DE can be rearranged so that all terms containing the independent variable are on one side of the equal sign and all terms containing the dependent
More informationDeterministic Global Optimization Algorithm and Nonlinear Dynamics
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
More informationNumerical approximation for optimal control problems via MPC and HJB. Giulia Fabrini
Numerical approximation for optimal control problems via MPC and HJB Giulia Fabrini Konstanz Women In Mathematics 15 May, 2018 G. Fabrini (University of Konstanz) Numerical approximation for OCP 1 / 33
More informationClick to 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
More informationPrelim Examination 2010 / 2011 (Assessing Units 1 & 2) MATHEMATICS. Advanced Higher Grade. Time allowed - 2 hours
Prelim Examination 00 / 0 (Assessing Units & ) MATHEMATICS Advanced Higher Grade Time allowed - hours Read Carefully. Calculators may be used in this paper.. Candidates should answer all questions. Full
More informationProjective geometry and spacetime structure. David Delphenich Bethany College Lindsborg, KS USA
Projective geometry and spacetime structure David Delphenich Bethany College Lindsborg, KS USA delphenichd@bethanylb.edu Affine geometry In affine geometry the basic objects are points in a space A n on
More informationNumerical Algorithms for ODEs/DAEs (Transient Analysis)
Numerical Algorithms for ODEs/DAEs (Transient Analysis) Slide 1 Solving Differential Equation Systems d q ( x(t)) + f (x(t)) + b(t) = 0 dt DAEs: many types of solutions useful DC steady state: state no
More informationNonlinear Stochastic Modeling and State Estimation of Weakly Observable Systems: Application to Industrial Polymerization Processes
Nonlinear Stochastic Modeling and State Estimation of Weakly Observable Systems: Application to Industrial Polymerization Processes Fernando V. Lima, James B. Rawlings and Tyler A. Soderstrom Department
More informationParameter Estimation in a Moving Horizon Perspective
Parameter Estimation in a Moving Horizon Perspective State and Parameter Estimation in Dynamical Systems Reglerteknik, ISY, Linköpings Universitet State and Parameter Estimation in Dynamical Systems OUTLINE
More informationNonlinear Model Predictive Control of a combined power and district heating production portfolio
Nonlinear Model Predictive Control of a combined power and district heating production portfolio Joakim Børlum Petersen Control and Automation, MSc. thesis Department of Electronic Systems June 2017 i
More informationCOMPARISON OF TWO METHODS TO SOLVE PRESSURES IN SMALL VOLUMES IN REAL-TIME SIMULATION OF A MOBILE DIRECTIONAL CONTROL VALVE
COMPARISON OF TWO METHODS TO SOLVE PRESSURES IN SMALL VOLUMES IN REAL-TIME SIMULATION OF A MOBILE DIRECTIONAL CONTROL VALVE Rafael ÅMAN*, Heikki HANDROOS*, Pasi KORKEALAAKSO** and Asko ROUVINEN** * Laboratory
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 informationMATLAB Ordinary Differential Equation (ODE) solver for a simple example 1. Introduction
MATLAB Ordinary Differential Equation (ODE) solver for a simple example 1. Introduction Differential equations are a convenient way to express mathematically a change of a dependent variable (e.g. concentration
More informationAn initialization subroutine for DAEs solvers: DAEIS
Computers and Chemical Engineering 25 (2001) 301 311 www.elsevier.com/locate/compchemeng An initialization subroutine for DAEs solvers: DAEIS B. Wu, R.E. White * Department of Chemical Engineering, Uniersity
More informationMath 5198 Mathematics for Bioscientists
Math 5198 Mathematics for Bioscientists Lecture 1: Course Conduct/Overview Stephen Billups University of Colorado at Denver Math 5198Mathematics for Bioscientists p.1/22 Housekeeping Syllabus CCB MERC
More informationA First Course on Kinetics and Reaction Engineering Example S5.1
Example S5.1 Problem Purpose This example illustrates the use of the MATLAB template file SolvBVDif.m to solve a second order boundary value ordinary differential equation. Problem Statement Solve the
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
More information1.1 Motivation: Why study differential equations?
Chapter 1 Introduction Contents 1.1 Motivation: Why stu differential equations?....... 1 1.2 Basics............................... 2 1.3 Growth and decay........................ 3 1.4 Introduction to Ordinary
More informationLinear and nonlinear optimization with CasADi
Numerical optimization course computer exercise inear and nonlinear optimization with CasADi Prof. Dr. Moritz Diehl 1, Prof. Dr. Angelika Altmann-Dieses 2, Adrian Bürger 1,2 1 Systems Control and Optimization
More informationROBUST CONSTRAINED ESTIMATION VIA UNSCENTED TRANSFORMATION. Pramod Vachhani a, Shankar Narasimhan b and Raghunathan Rengaswamy a 1
ROUST CONSTRINED ESTIMTION VI UNSCENTED TRNSFORMTION Pramod Vachhani a, Shankar Narasimhan b and Raghunathan Rengaswamy a a Department of Chemical Engineering, Clarkson University, Potsdam, NY -3699, US.
More informationMTAEA Implicit Functions
School of Economics, Australian National University February 12, 2010 Implicit Functions and Their Derivatives Up till now we have only worked with functions in which the endogenous variables are explicit
More informationF, {z} Isothermal Flash Tank
Exam III dministered: Thursday, November 7, 24 22 points For each problem part: points i not attempted or no work shown, point or partial credit, i work is shown, 2 points or correct numerical value o
More informationIndex Reduction and Discontinuity Handling using Substitute Equations
Mathematical and Computer Modelling of Dynamical Systems, vol. 7, nr. 2, 2001, 173-187. Index Reduction and Discontinuity Handling using Substitute Equations G. Fábián, D.A. van Beek, J.E. Rooda Abstract
More informationStreamline calculations. Lecture note 2
Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational
More informationDifferential-Algebraic Equations (DAEs)
Differential-Algebraic Equations (DAEs) L. T. Biegler Chemical Engineering Department Carnegie Mellon University Pittsburgh, PA 15213 biegler@cmu.edu http://dynopt.cheme.cmu.edu Introduction Simple Examples
More informationALADIN An Algorithm for Distributed Non-Convex Optimization and Control
ALADIN An Algorithm for Distributed Non-Convex Optimization and Control Boris Houska, Yuning Jiang, Janick Frasch, Rien Quirynen, Dimitris Kouzoupis, Moritz Diehl ShanghaiTech University, University of
More informationModel Predictive Control
Model Predictive Control Linear time-varying and nonlinear MPC Alberto Bemporad http://cse.lab.imtlucca.it/~bemporad 2018 A. Bemporad - "Model Predictive Control" 1/30 Course structure Linear model predictive
More informationPhysical Modelling with Simscape Rick Hyde
Physical Modelling with Simscape Rick Hyde 1 2013 The MathWorks, Inc. Outline Part 1: Introduction to Simscape Review approaches to modelling Overview of Simscape-based libraries Introduction to physical
More informationSolving Differential Equations Using MATLAB
Solving Differential Equations Using MATLAB Abraham Asfaw aasfaw.student@manhattan.edu November 28, 2011 1 Introduction In this lecture, we will follow up on lecture 2 with a discussion of solutions to
More informationDefect-based a-posteriori error estimation for implicit ODEs and DAEs
1 / 24 Defect-based a-posteriori error estimation for implicit ODEs and DAEs W. Auzinger Institute for Analysis and Scientific Computing Vienna University of Technology Workshop on Innovative Integrators
More informationx+ y = 50 Dividing both sides by 2 : ( ) dx By (7.2), x = 25m gives maximum area. Substituting this value into (*):
Solutions 7(b 1 Complete solutions to Exercise 7(b 1. Since the perimeter 100 we have x+ y 100 [ ] ( x+ y 50 ividing both sides by : y 50 x * The area A xy, substituting y 50 x gives: A x( 50 x A 50x x
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationMATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem
MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem Pulemotov, September 12, 2012 Unit Outline Goal 1: Outline linear
More informationProblem 1 Cost of an Infinite Horizon LQR
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS H O M E W O R K # 5 Ahmad F. Taha October 12, 215 Homework Instructions: 1. Type your solutions in the LATEX homework
More informationIntroduction to Differential Equations
Math0 Lecture # Introduction to Differential Equations Basic definitions Definition : (What is a DE?) A differential equation (DE) is an equation that involves some of the derivatives (or differentials)
More informationABSTRACT The use of dynamic modelling and simulation of thermal power systems has grown recent years and as the performance of computers have increased optimisation of such systems has become possible.
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
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 informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationEigenvalues and eigenvectors System Theory: electricmotor
Eigenvalues and eigenvectors System Theory: electricmotor Alireza Abouhossein Ph.D. Alireza.abouhossein@unibo.it 1 Example-1 Can you calculate A^2, A -1 -I? without multiplying A or finding the inverse
More informationState and Parameter Estimation for Natural Gas Pipeline Networks using Transient State Data
JOURNAL OF LATEX CLASS FILES, VOL. 14, NO. 8, NOVEMBER 216 1 State and Parameter Estimation for Natural Gas Pipeline Networks using Transient State Data Kaarthik Sundar Anatoly Zlotnik arxiv:183.7156v2
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationEfficient Numerical Methods for Nonlinear MPC and Moving Horizon Estimation
Efficient Numerical Methods for Nonlinear MPC and Moving Horizon Estimation Moritz Diehl, Hans Joachim Ferreau, and Niels Haverbeke Optimization in Engineering Center (OPTEC) and ESAT-SCD, K.U. Leuven,
More informationThe Euler Equations and Non-Local Conservative Riccati Equations
The Euler Equations and Non-Local Conservative Riccati Equations Peter Constantin Department of Mathematics The University of Chicago November 8, 999 The purpose of this brief note is to present an infinite
More informationUse of Differential Equations In Modeling and Simulation of CSTR
Use of Differential Equations In Modeling and Simulation of CSTR JIRI VOJTESEK, PETR DOSTAL Department of Process Control, Faculty of Applied Informatics Tomas Bata University in Zlin nám. T. G. Masaryka
More informationEconomics 101A (Lecture 3) Stefano DellaVigna
Economics 101A (Lecture 3) Stefano DellaVigna January 24, 2017 Outline 1. Implicit Function Theorem 2. Envelope Theorem 3. Convexity and concavity 4. Constrained Maximization 1 Implicit function theorem
More informationNonlinear Model Predictive Control for Combined Cycle Power Plants. Anna Johnsson
Nonlinear Model Predictive Control for Combined Cycle Power Plants Anna Johnsson MASTERS THESIS in Automatic Control 2013 Supervisor: Stéphane Velut, Modelon Killian Link, Siemens Department of Automatic
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Pseudo-Time Integration 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 10 Outline 1
More informationIntegration, Separation of Variables
Week #1 : Integration, Separation of Variables Goals: Introduce differential equations. Review integration techniques. Solve first-order DEs using separation of variables. 1 Sources of Differential Equations
More informationNonlinear State Estimation Methods Overview and Application to PET Polymerization
epartment of Biochemical and Chemical Engineering Process ynamics Group () onlinear State Estimation Methods Overview and Polymerization Paul Appelhaus Ralf Gesthuisen Stefan Krämer Sebastian Engell The
More informationMath 3313: Differential Equations Second-order ordinary differential equations
Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties
More informationSolutions Definition 2: a solution
Solutions As was stated before, one of the goals in this course is to solve, or find solutions of differential equations. In the next definition we consider the concept of a solution of an ordinary differential
More informationMathematical Modelling Using Simulink
Experiment Two Mathematical Modelling Using Simulink Control Systems Laboratory Dr. Zaer Abo Hammour Dr. Zaer Abo Hammour Control Systems Laboratory 1. Mathematical Model Definition A mathematical model
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More information10.34 Numerical Methods Applied to Chemical Engineering
10.34 Numerical Methods Applied to Chemical Engineering Quiz 1 This quiz consists of three problems worth 0, 40 and 40 points respectively. The problem statements are found on pages, 3 and 5 in this exam
More informationModule 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 3 : Analytical Solutions of Linear ODE-IVPs
Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 3 : Analytical Solutions of Linear ODE-IVPs 3 Analytical Solutions of Linear ODE-IVPs Before developing numerical
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 informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems. CDS 110b
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 110b R. M. Murray Kalman Filters 14 January 2007 Reading: This set of lectures provides a brief introduction to Kalman filtering, following
More information. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,
1.2. Direction Fields: Graphical Representation of the ODE and its Solution Section Objective(s): Constructing Direction Fields. Interpreting Direction Fields. Definition 1.2.1. A first order ODE of the
More informationMEM 255 Introduction to Control Systems: Modeling & analyzing systems
MEM 55 Introduction to Control Systems: Modeling & analyzing systems Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline The Pendulum Micro-machined capacitive accelerometer
More informationEconomic Model Predictive Control Historical Perspective and Recent Developments and Industrial Examples
1 Economic Model Predictive Control Historical Perspective and Recent Developments and Industrial Examples Public Trial Lecture Candidate: Vinicius de Oliveira Department of Chemical Engineering, Faculty
More informationMathematical Model of Catalytic Chemical Reactor. Sintayehu Agegnehu Matintu* School of Mathematical and Statistical Sciences, Hawassa University,
Mathematical Model of Catalytic Chemical Reactor Sintayehu Agegnehu Matintu* School of Mathematical and Statistical Sciences, Hawassa University, P.O.Box 05, Hawassa, Ethiopia Abstract This study concerns
More informationAsymptotic and Perturbation Methods: Lecture Notes.
Asymptotic and Perturbation Methods: Lecture Notes. Blerta Shtylla University of Utah Mathematics Department Fall 2008 1 Lecture 1 and Lecture 2 Week of August 26, 2008 Introduction: For many problems
More informationCombined Estimation and Optimal Control of Batch Membrane Processes
processes Article Combined Estimation and Optimal Control of Batch Membrane Processes Martin Jelemenský 1, *, Daniela Pakšiová 1, Radoslav Paulen 2, Abderrazak Latifi 3 and Miroslav Fikar 1 1 Faculty of
More informationSOLUTIONS TO ADDITIONAL EXERCISES FOR II.1 AND II.2
SOLUTIONS TO ADDITIONAL EXERCISES FOR II.1 AND II.2 Here are the solutions to the additional exercises in betsepexercises.pdf. B1. Let y and z be distinct points of L; we claim that x, y and z are not
More informationCh 4 Differentiation
Ch 1 Partial fractions Ch 6 Integration Ch 2 Coordinate geometry C4 Ch 5 Vectors Ch 3 The binomial expansion Ch 4 Differentiation Chapter 1 Partial fractions We can add (or take away) two fractions only
More informationControl Systems Lab - SC4070 Control techniques
Control Systems Lab - SC4070 Control techniques Dr. Manuel Mazo Jr. Delft Center for Systems and Control (TU Delft) m.mazo@tudelft.nl Tel.:015-2788131 TU Delft, February 16, 2015 (slides modified from
More informationAn Ensemble Kalman Filter for Systems Governed by Differential Algebraic Equations (DAEs)
Preprints of the 8th IFAC Symposium on Advanced Control of Chemical Processes The International Federation of Automatic Control An Ensemble Kalman Filter for Systems Governed by Differential Algebraic
More informationChapter III. Unconstrained Univariate Optimization
1 Chapter III Unconstrained Univariate Optimization Introduction Interval Elimination Methods Polynomial Approximation Methods Newton s Method Quasi-Newton Methods 1 INTRODUCTION 2 1 Introduction Univariate
More informationThe Method of Undetermined Coefficients.
The Method of Undetermined Coefficients. James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University May 24, 2017 Outline 1 Annihilators 2 Finding The
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 4, pp. 373-38, 23. Copyright 23,. ISSN 68-963. ETNA A NOTE ON THE RELATION BETWEEN THE NEWTON HOMOTOPY METHOD AND THE DAMPED NEWTON METHOD XUPING ZHANG
More informationTransduction Based on Changes in the Energy Stored in an Electrical Field
Lecture 6-1 Transduction Based on Changes in the Energy Stored in an Electrical Field Electric Field and Forces Suppose a charged fixed q 1 in a space, an exploring charge q is moving toward the fixed
More informationReal-Time Implementation of Nonlinear Predictive Control
Real-Time Implementation of Nonlinear Predictive Control Michael A. Henson Department of Chemical Engineering University of Massachusetts Amherst, MA WebCAST October 2, 2008 1 Outline Limitations of linear
More informationNonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México
Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October
More informationA Crash Course on Kalman Filtering
A Crash Course on Kalman Filtering Dan Simon Cleveland State University Fall 2014 1 / 64 Outline Linear Systems Probability State Means and Covariances Least Squares Estimation The Kalman Filter Unknown
More informationODE solvers in Julia. Gabriel Ingesson. October 2, 2015
ODE solvers in Julia Gabriel Ingesson October 2, 2015 Motivation General: Numerical methods for solving ODE s is important for system simulations. Simulation is important for controller design. Personal:
More informationChapter III. Stability of Linear Systems
1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,
More informationDefect-based a posteriori error estimation for index-1 DAEs
Defect-based a posteriori error estimation for index-1 DAEs W. Auzinger, H. Lehner, E. Weinmüller Institute for Analysis and Scientific Computing Vienna University of Technology Wiedner Hauptstrasse 8
More informationChapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey
Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix
More informationON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS
8th International Congress on Computational Mechanics, Volos, 1-15 July 015 ON THE INTEGRATION OF EQUATIONS OF MOTION: FEM AND MOLECULAR DYNAMICS PROBLEMS E.G. Kakouris, V.K. Koumousis Institute of Structural
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More informationA First Course on Kinetics and Reaction Engineering. Class 20 on Unit 19
A First Course on Kinetics and Reaction Engineering Class 20 on Unit 19 Part I - Chemical Reactions Part II - Chemical Reaction Kinetics Where We re Going Part III - Chemical Reaction Engineering A. Ideal
More informationCHAPTER 3 : MATHEMATICAL MODELLING PRINCIPLES
CHAPTER 3 : MATHEMATICAL MODELLING PRINCIPLES When I complete this chapter, I want to be able to do the following. Formulate dynamic models based on fundamental balances Solve simple first-order linear
More informationThe Plan. Initial value problems (ivps) for ordinary differential equations (odes) Review of basic methods You are here: Hamiltonian systems
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 2 Dianne P. O Leary c 2008 The Plan
More information