Lecture 1: Basics Concepts
|
|
- Harry James
- 6 years ago
- Views:
Transcription
1 National Technical University of Athens School of Chemical Engineering Department II: Process Analysis and Plant Design Lecture 1: Basics Concepts Instructor: Α. Kokossis Laboratory teaching staff: Α. Nikolakopoulos
2 Review Basic Concepts Process Simulation i. Unit simulation vs flowsheeting ii. iii. Mathematical representation Simulation types - Sequential modelling - Simultaneous approach - Equation-oriented approach
3 1. Introduction
4 Simulation of individual processes i. Simulation of process models (unit operations) Operating parameters (pressure, temperature) Feed/input Flows Composition Pressure Temperature Process Products/output Flows Composition Pressure Temperature Design parameters (Sizing, geometry)
5 Flowsheeting ii. Industrial flowsheeting Process flowsheeting as integration of process models Computational work: develop calculations around - Mass and energy balances - Scale-up and sizing - Costing with respect to the integrated processes of an industrial plant at steady state. The calculations assume individual models for the process units and the focus shifts to combine the individual calculations as a whole.
6 The role of flowsheeting in Process Design Needs Flowsheet Synthesis Analysis Mass and energy balances Costing/sizing Decision variables (Initial values) Parameter optimization Process selection Economic Evaluation Final Flowsheet
7 Potential in process design Advantages i. Significant reductions in the development time Use of unit operation libraries Property models and thermodynamics Correlations and supporting functions ii. Significant improvements of the existing process Experimentation with alternative processes (not systematic though!) what if studies and analysis Parametric optimization studies Archiving and re-use of design technology
8 Generic process representation a. Process unit representation Input/output structure Distinction between input/output variables and parameters input, x i Units Output, y i parameters, u i i.e.: f i (x i, y i, u i ) = 0, then, y i = g i (x i, u i ) fixed calculated
9 Open and closed systems Open systems: I/O dependencies are described by analytical expressions. They could be - algebraic equations - differential equations (simple, with partial derivatives) General formulation: f j (x 1, x 2, x m ) = 0, y k = g (x 1, x 2, x m ) ; x i (0) = 0 j = 1,..., m k =1,..., n i = 1,..., m The remaining of the models assumes algebraic sets of equations
10 Systems representation Equations variables - parameters residuals NE = n f 1 (x 1, x 2,, x m ) = 0 = r 1 f 2 (x 1, x 2,, x m ) = 0 = r 2... f n (x 1, x 2,, x m ) = 0 = r n NV = m
11 Design variables and selection of specs a. In simulation NE = NV. If NE>NV, the remaining NV NE = NP variables (design variables) have to be selected and fixed (set them as parameters) Distillation: product purity Reaction: conversion Flowsheet: recycle/purge stream b. Note: The parameters of one model may be variables in another. The role of inputs and parameters may be swapped across models. Examples of typical design variables set as parameters: c. In closed models, the acceptable combinations of parameters is not clear and is usually provided as user options (e.g. distillation: RR over Reflux Boilup or Reboiler Heat Duty)
12 Degree of freedom analysis When we define additional equations, NΕ > NV the systems becomes overdefined, while for NV > NE there are degrees of freedom (i.e. an optimization problem) In the past a degree of freedom analysis has been an important task in design. Nowadays, such analysis can be performed by the flowsheeting software.
13 Closed systems Open systems: minority of the process models currently used. Closed: the relationship between outputs and inputs is indirect. Examples algorithms and procedural code widely used as Fortran code for instance. Compare with Open systems: direct relationships through analytical expressions (e.g. equations) Can map open and closed systems with each other? Open Closed Open? Closed (apparent mapping) Sets of equations: defined as the residuals from closed-system calculations For every r r y i i ( x, y, z ) r, i i i i xi u Δx i i Δyi i (closed models) (Jacobian element)
14 Open and closed systems Relationships and mappings between different technologies inputs parameters outputs x 1, x 2,..., x n u 1, u 2,..., u n y 1, y 2,..., y n f 1 = r 1 f 2 = r 2... αij βij γij f n = r n f n equations that correspond to the open model x x + x y + y (closed model) i i i It calculates r i as a results of the changes of ri ri ri It calculates,, α ij, β ij, γ ij x y u j j j ( ) x, y, u i i i
15 2. Flowsheeting technology
16 Types and technologies i. Sequential modular approach ii. Simultaneous modular approach iii. Equation-oriented approach Let us use an illustration as a drive to explain the technologies
17 Illustration of a simple flowsheet Typical reactor-separation recycle structure A simple flowsheet Mixer Reactor Separator Flowsheet includes: Reactor, separator 2 mixers 1 recycle Mixer
18 Input-output representation Systems representation Four (4) subsystems y 31 x 12 y 21 x 11 y 11 Mixer Reactor x21 x31 Separator Two indices used for each flow The first denotes the process x 42 y 32 The second denotes its role as input or output x 41 Mixer y 41
19 Mathematical representation and models Sets of equations (constraints) for the flowsheet 1. Constraints for each unit (balances: mass, energy etc; economics) mixer1: f1 ( x 11,x 12,y11 ) = 0 reactor: f2 ( x 21,x 21, u2 ) = 0 f ( x,y,y, u ) = 0 seperator: mixer2: f4 ( x 41,x 42,y41 ) = 0 x 12 y 31 Generalization: x ij y ij u i input j in unit i output j from unit i parameter of unit i y 21 x 11 y 11 Mixer Reactor x21 x31 x 41 y 32 x 42 Mixer Separator y 41 Any other constraints;
20 Mathematical representation (2) Connectivity constraints 1. Connectivity x12 y31 = 0 x21 y11 = 0 x31 y21 = 0 x42 y32 = 0 x 12 y 21 x 11 y 11 Mixer Reactor x21 x31 y 31 Separator y 32 x 42 Notes: unit constraints relate to mass and energy balances, constraints for physical properties, correlations, economics, upper/lower limits Unit parameters (design parameters) relate to specifications (e.g. temperature and pressure of a reactor) required to impose user preferences and/or make the balances a square set of equations and variables x 41 Mixer y 41
21 Sequential modular approach Basic Idea 1. Simulate each process in sequence following the physical flow of streams (though not necessarily) 2. Whenever needed, guess streams from forward or backward units. These are termed tear streams and their guesses are the pivots of convergence 3. Tear stream variables are calculated over cycles of calculations 4. Overall convergence is attained when - each and every process model converges - each and every tear stream converges Commercial software Aspen Plus (Evans et al, 1979) HYSIS (Mahoney & Santollani, 1994) SUPERPO (Futtellugen)
22 Sequential modular approach (illustration) Example S2 y 31 x x 11 y 11 S3 y 21 S4 Mixer Reactor Separator S1 x21 x31 1 y 32 x 42 S5 x 41 S6 Mixer 4 y 41 S7
23 Sequential modular approach (illustration) Use tearing based on natural flows Using such tearing, the sequence of calculations is ( ) = = ( ) f x,x,y 0 y g x,x ( u ) = = ( u ) f x,x, 0 y g x, ( u ) = = ( u ) y = g ( x, u ) f x,y,y, 0 y g x, x 12 (tear at ) convergence (no tearing) convergence (no tearing) convergence (no tearing) convergence ( ) = = ( ) f x,x,y 0 y g x,x Convergence in x 12? If yes, stop. Otherwise, iterate
24 The sequential modular approach as an algorithm Sequential modular steps y 31 x 12 y 21 x 11 y 11 Mixer Reactor x21 x31 Separator x 41 y 32 x 42 Mixer y 41
25 Tearing and its importance in convergence Remarks In our case, the single tearing stream has been the recycle calculated as: ( ) x = f x,y k ' k 1 k There is no unique option for tearing. Indeed, there exist degrees of freedom around the tearing options and they can be exploited for better convergence Following natural flows in the flowsheet is not the best advice for the selection of tear streams Internally, the tearing is handled without any request by the user to define the tearing variables. However, the selection of tearing variables can be critical. Tearing may destabilize convergence. The case is pronounced when the number of recycles increase.
26 Tearing options are possible even for simple flowsheets Alternative tearing scenarios y 31 x 12 y 21 x 11 y 11 Mixer Reactor x21 x31 Separator x 41 y 32 x 42 Mixer y 41 Tearing around reactors, separators, etc
27 Simultaneous modular approach Simulation trade-offs and challenges: The sequential modular approach faces difficulties to converge in large problems with several recycle streams. The simulation approach may have to iterate around distant (process unit) blocks The approach holds back on flowsheet calculations to ensure (robust) convergence for the individual blocks Basic idea Simultaneously converge all individual subsystems by linearizing (simplifying) around nominal operating points the I/Os of individual blocks Commercial software Flowpack (predecessor of today s packages) Rosen (1962) Umeda et al. (1972), Reklaitis (1979), Mahalec and Motard (1979), Biegler & Hughes (1982) All units are written as in sequential modular approach. i.e. : = ( u ) α y g x, y x ij i ik i ij ijk ik k ij k= 1
28 Equation-oriented simulation Motivation: Linearization is already made within numerical routines (e.g. Taylor approximations typically replace nonlinear sets of equations in block balances) Advances in numerical computing to handle larger problems Basic Idea Tackle the full problem as a single set of equations, building custom-made functions to handle block-diagonal structures Commercial packages gproms Aspen Plus / HYSIS (though not truly through analytical models)
29 Equation-oriented approach Sets of equations 1. Process unit models 2. Connectivity constraints 3. Design specifications Typical challenges Break down the large system into smaller tasks Degree of freedom analysis Handle multiplicity in solutions Non-linear aspects of the problem constraints Initialization Management of side calculations related to 1 st and 2 nd order derivatives (Jacobian, Hessian)
30 Pros and cons: sequential modular vs equation-oriented Sequential modular Equation-oriented Separate simulation for each unit Modelling is not as flexible Robust convergence Initialization is useful Small memory storage Applications are easy to set up (little left to the user) Simultaneous solution of equations Modelling is very flexible Lack of robustness Initialization is critical Larger memory storage Applications take longer times and require more sophistication (much left to the user)
31 Hybrid applications Combine Sequential modular (closed systems) used initially to establish good reference (nominal) points with Equation-oriented simulations where previous (nominal) solutions can be used as Initial points Handle stiff sets of non-linear equations
Pressure Swing Distillation with Aspen Plus V8.0
Pressure Swing Distillation with Aspen Plus V8.0 1. Lesson Objectives Aspen Plus property analysis RadFrac distillation modeling Design Specs NQ Curves Tear streams Understand and overcome azeotrope Select
More informationTEACH YOURSELF THE BASICS OF ASPEN PLUS
TEACH YOURSELF THE BASICS OF ASPEN PLUS TEACH YOURSELF THE BASICS OF ASPEN PLUS RALPH SCHEFFLAN Chemical Engineering and Materials Science Department Stevens Institute of Technology A JOHN WILEY & SONS,
More informationProcess Unit Control System Design
Process Unit Control System Design 1. Introduction 2. Influence of process design 3. Control degrees of freedom 4. Selection of control system variables 5. Process safety Introduction Control system requirements»
More informationy b where U. matrix inverse A 1 ( L. 1 U 1. L 1 U 13 U 23 U 33 U 13 2 U 12 1
LU decomposition -- manual demonstration Instructor: Nam Sun Wang lu-manualmcd LU decomposition, where L is a lower-triangular matrix with as the diagonal elements and U is an upper-triangular matrix Just
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 informationIntroduction to Chemical Engineering Computing, 2 nd edition, Bruce A. Finlayson, Wiley (2012)
Introduction to Chemical Engineering Computing, 2 nd edition, Bruce A. Finlayson, Wiley (2012) The following web sites can be used to obtain material from Wiley. It includes the keys to the problems (useful
More informationApplication of Decomposition Methodology to Solve Integrated Process Design and Controller Design Problems for Reactor-Separator-Recycle Systems
Proceedings of the 9th International Symposium on Dynamics and Control of Process Systems (DYCOPS 2010), Leuven, Belgium, July 5-7, 2010 Mayuresh Kothare, Moses Tade, Alain Vande Wouwer, Ilse Smets (Eds.)
More informationChapter 4. Problem SM.7 Ethylbenzene/Styrene Column
Background Chapter 4. Problem SM.7 Ethylbenzene/Styrene Column In Problem SM.6 of the HYSYS manual, a modified form of successive substitution, called the Wegstein method, was used to close the material
More informationSlide 1. A look into the CAPE-OPEN kitchen of COCO. Jasper van Baten, AmsterCHEM
Slide 1 A look into the CAPE-OPEN kitchen of COCO Jasper van Baten, AmsterCHEM Slide 2 CAPE-OPEN to CAPE-OPEN (COCO): Simulation environment (COFE) Thermodynamic property package (TEA) Collection of unit
More informationSystems Engineering Spring Group Project #1: Process Flowsheeting Calculations for Acetic Anhydride Plant. Date: 2/25/00 Due: 3/3/00
10.551 Systems Engineering Spring 2000 Group Project #1: Process Flowsheeting Calculations for Acetic Anhydride Plant Date: 2/25/00 Due: 3/3/00 c Paul I. Barton, 14th February 2000 At our Nowhere City
More informationTHE PROCESS SIMULATION COURSE THE CULMINATION OF CORE UNDERGRADUATE COURSEWORK IN CHEMICAL ENGINEERING
THE PROCESS SIMULATION COURSE THE CULMINATION OF CORE UNDERGRADUATE COURSEWORK IN CHEMICAL ENGINEERING Mordechai Shacham, Ben-Gurion University, Beer-Sheva, Israel Michael B. Cutlip, University of Connecticut,
More informationMatrix Assembly in FEA
Matrix Assembly in FEA 1 In Chapter 2, we spoke about how the global matrix equations are assembled in the finite element method. We now want to revisit that discussion and add some details. For example,
More informationGLoBES. Patrick Huber. Physics Department VT. P. Huber p. 1
GLoBES Patrick Huber Physics Department VT P. Huber p. 1 P. Huber p. 2 General Long Baseline Experiment Simulator GLoBES is a software package designed for Simulation Analysis Comparison of neutrino oscillation
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 informationINTRODUCTION TO CHEMICAL ENGINEERING COMPUTING
INTRODUCTION TO CHEMICAL ENGINEERING COMPUTING BRUCE A. FINLÄYSON, PH.D. University of Washington Seattle, Washington iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC., PUBLICATION Microsoft product screen
More informationOptimization of Batch Processes
Integrated Scheduling and Dynamic Optimization of Batch Processes Yisu Nie and Lorenz T. Biegler Center for Advanced Process Decision-making Department of Chemical Engineering Carnegie Mellon University
More informationA Definition for Plantwide Controllability. Process Flexibility
A Definition for Plantwide Controllability Surya Kiran Chodavarapu and Alex Zheng Department of Chemical Engineering University of Massachusetts Amherst, MA 01003 Abstract Chemical process synthesis typically
More informationIntegrated Knowledge Based System for Process Synthesis
17 th European Symposium on Computer Aided Process Engineering ESCAPE17 V. Plesu and P.S. Agachi (Editors) 2007 Elsevier B.V. All rights reserved. 1 Integrated Knowledge Based System for Process Synthesis
More informationSystems of Nonlinear Equations
Systems of Nonlinear Equations L. T. Biegler Chemical Engineering Department Carnegie Mellon University Pittsburgh, PA 15213 Biegler@cmu.edu http://dynopt.cheme.cmu.edu Obtaining Derivatives Broyden vs.
More informationDividing wall columns for heterogeneous azeotropic distillation
Dividing wall columns for heterogeneous azeotropic distillation Quang-Khoa Le 1, Ivar J. Halvorsen 2, Oleg Pajalic 3, Sigurd Skogestad 1* 1 Norwegian University of Science and Technology (NTNU), Trondheim,
More informationA Flexible Framework and Model Library for Process Simulation, Optimization and Control
Mario R. Eden, Marianthi Ierapetritou Gavin P. Towler (Editors) Proceedings of the 13 th International Symposium on Process Systems Engineering PSE 2018 July 1-5, 2018, San Diego, California, USA 2018
More informationOptimization of the Sulfolane Extraction Plant Based on Modeling and Simulation
Korean J. Chem. Eng., 17(6), 712-718 (2000) Optimization of the Sulfolane Extraction Plant Based on Modeling and Simulation Yu-Jung Choi, Tae-In Kwon and Yeong-Koo Yeo Department of Chemical Engineering,
More informationDesign and Control Properties of Arrangements for Distillation of Four Component Mixtures Using Less Than N-1 Columns
D. M. MÉNDEZ-VALENCIA et al., Design and Control Properties of Arrangements, Chem. Biochem. Eng. Q. 22 (3) 273 283 (2008) 273 Design and Control Properties of Arrangements for Distillation of Four Component
More informationShortcut Distillation. Agung Ari Wibowo, S.T., M.Sc Politeknik Negeri Malang Malang - Indonesia
Shortcut Distillation Agung Ari Wibowo, S.T., M.Sc Politeknik Negeri Malang Malang - Indonesia The use of separation column in HYSYS The column utilities in HYSYS can be used to model a wide variety of
More informationThe Chemical Kinetics Time Step a detailed lecture. Andrew Conley ACOM Division
The Chemical Kinetics Time Step a detailed lecture Andrew Conley ACOM Division Simulation Time Step Deep convection Shallow convection Stratiform tend (sedimentation, detrain, cloud fraction, microphysics)
More informationProcess Design Decisions and Project Economics Prof. Dr. V. S. Moholkar Department of Chemical Engineering Indian Institute of Technology, Guwahati
Process Design Decisions and Project Economics Prof. Dr. V. S. Moholkar Department of Chemical Engineering Indian Institute of Technology, Guwahati Module - 2 Flowsheet Synthesis (Conceptual Design of
More informationThermally Coupled Distillation Systems: Study of an Energy-efficient Reactive Case
F. O. BARROSO-MUÑOZ et al., Thermally Coupled Distillation Systems: Study of, Chem. Biochem. Eng. Q. 21 (2) 115 120 (2007) 115 Thermally Coupled Distillation Systems: Study of an Energy-efficient Reactive
More informationSemi-empirical Process Modelling with Integration of Commercial Modelling Tools
European Symposium on Computer Arded Aided Process Engineering 15 L. Puigjaner and A. Espuña (Editors) 2005 Elsevier Science B.V. All rights reserved. Semi-empirical Process Modelling with Integration
More informationAC : A MODULE FOR TEACHING BATCH OPERATIONS
AC 2010-726: A MODULE FOR TEACHING BATCH OPERATIONS Richard Turton, West Virginia University Richard Turton received his B.S. degree from the University of Nottingham and his M.S. and Ph.D. degrees from
More informationCopyrighted Material. 1.1 Large-Scale Interconnected Dynamical Systems
Chapter One Introduction 1.1 Large-Scale Interconnected Dynamical Systems Modern complex dynamical systems 1 are highly interconnected and mutually interdependent, both physically and through a multitude
More informationSIMULATION ANALYSIS OF FULLY THERMALLY COUPLED DISTILLATION COLUMN
Int. J. Chem. Sci.: 14(3), 2016, 1621-1632 ISSN 0972-768X www.sadgurupublications.com SIMULATION ANALYSIS OF FULLY THERMALLY COUPLED DISTILLATION COLUMN ASMITA PRAVIN PATIL * and S. M. JADHAV Chemical
More informationSIMULTANEOUS OPTIMAL DESIGN AND OPERATION OF MULTIPURPOSE BATCH DISTILLATION COLUMNS
SIMULTANEOUS OPTIMAL DESIGN AND OPERATION OF MULTIPURPOSE BATCH DISTILLATION COLUMNS Kian Huat Low and Eva Sorensen 1 Centre for Process Systems Engineering Department of Chemical Engineering University
More informationEffect of Two-Liquid Phases on the Dynamic and the Control of Trayed Distillation Columns
Effect of Two-Liquid Phases on the Dynamic and the Control of Trayed Distillation Columns Gardênia Cordeiro, Karoline Brito, Brenda Guedes, Luis Vasconcelos, Romildo Brito Federal University of Campina
More information2x (x 2 + y 2 + 1) 2 2y. (x 2 + y 2 + 1) 4. 4xy. (1, 1)(x 1) + (1, 1)(y + 1) (1, 1)(x 1)(y + 1) 81 x y y + 7.
Homework 8 Solutions, November 007. (1 We calculate some derivatives: f x = f y = x (x + y + 1 y (x + y + 1 x = (x + y + 1 4x (x + y + 1 4 y = (x + y + 1 4y (x + y + 1 4 x y = 4xy (x + y + 1 4 Substituting
More informationNONLINEAR EQUATIONS AND TAYLOR S THEOREM
APPENDIX C NONLINEAR EQUATIONS AND TAYLOR S THEOREM C.1 INTRODUCTION In adjustment computations it is frequently necessary to deal with nonlinear equations. For example, some observation equations relate
More informationA First Course on Kinetics and Reaction Engineering Supplemental Unit S4. Numerically Fitting Models to Data
Supplemental Unit S4. Numerically Fitting Models to Data Defining the Problem Many software routines for fitting a model to experimental data can only be used when the model is of a pre-defined mathematical
More informationHeterogeneous Azeotropic Distillation Operational Policies and Control
Heterogeneous Azeotropic Distillation Operational Policies and Control Claudia J. G. Vasconcelos * and Maria Regina Wolf-Maciel State University of Campinas, School of Chemical Engineering, Campinas/SP,
More informationDeterministic Models
Deterministic Models Perfect foreight, nonlinearities and occasionally binding constraints Sébastien Villemot CEPREMAP June 10, 2014 Sébastien Villemot (CEPREMAP) Deterministic Models June 10, 2014 1 /
More informationAspen Dr. Ziad Abuelrub
Aspen Plus Lab Pharmaceutical Plant Design Aspen Dr. Ziad Abuelrub OUTLINE 1. Introduction 2. Getting Started 3. Thermodynamic Models & Physical Properties 4. Pressure Changers 5. Heat Exchangers 6. Flowsheet
More informationDISTILLATION. Keywords: Phase Equilibrium, Isothermal Flash, Adiabatic Flash, Batch Distillation
25 DISTILLATION Keywords: Phase Equilibrium, Isothermal Flash, Adiabatic Flash, Batch Distillation Distillation refers to the physical separation of a mixture into two or more fractions that have different
More informationIntroduction to Linear Programming (LP) Mathematical Programming (MP) Concept (1)
Introduction to Linear Programming (LP) Mathematical Programming Concept LP Concept Standard Form Assumptions Consequences of Assumptions Solution Approach Solution Methods Typical Formulations Massachusetts
More informationThe Abinit project. Coding is based on modern software engineering principles
The Abinit project Abinit is a robust, full-featured electronic-structure code based on density functional theory, plane waves, and pseudopotentials. Abinit is copyrighted and distributed under the GNU
More informationDirect Methods. Moritz Diehl. Optimization in Engineering Center (OPTEC) and Electrical Engineering Department (ESAT) K.U.
Direct Methods Moritz Diehl Optimization in Engineering Center (OPTEC) and Electrical Engineering Department (ESAT) K.U. Leuven Belgium Overview Direct Single Shooting Direct Collocation Direct Multiple
More informationSimulation of Butyl Acetate and Methanol Production by Transesterification Reaction via Conventional Distillation Process
Simulation of Butyl Acetate and Methanol Production by Transesterification Reaction via Conventional Distillation Process Nikhil V. Sancheti Department of Chemical Engineering L.I.T., Nagpur, Maharashtra,
More informationReactor Design within Excel Enabled by Rigorous Physical Properties and an Advanced Numerical Computation Package
Reactor Design within Excel Enabled by Rigorous Physical Properties and an Advanced Numerical Computation Package Mordechai Shacham Department of Chemical Engineering Ben Gurion University of the Negev
More informationDETERMINATION OF OPTIMAL ENERGY EFFICIENT SEPARATION SCHEMES BASED ON DRIVING FORCES
DETERMINATION OF OPTIMAL ENERGY EFFICIENT SEPARATION SCHEMES BASED ON DRIVING FORCES Abstract Erik Bek-Pedersen, Rafiqul Gani CAPEC, Department of Chemical Engineering, Technical University of Denmark,
More informationSolving linear systems (6 lectures)
Chapter 2 Solving linear systems (6 lectures) 2.1 Solving linear systems: LU factorization (1 lectures) Reference: [Trefethen, Bau III] Lecture 20, 21 How do you solve Ax = b? (2.1.1) In numerical linear
More informationMath and Numerical Methods Review
Math and Numerical Methods Review Michael Caracotsios, Ph.D. Clinical Associate Professor Chemical Engineering Department University of Illinois at Chicago Introduction In the study of chemical engineering
More informationLecture 4: Linear Algebra 1
Lecture 4: Linear Algebra 1 Sourendu Gupta TIFR Graduate School Computational Physics 1 February 12, 2010 c : Sourendu Gupta (TIFR) Lecture 4: Linear Algebra 1 CP 1 1 / 26 Outline 1 Linear problems Motivation
More informationSteady State Multiplicity of Novel Reactive Distillation Combined Distillation Column with Side Reactors for the Production of Cyclohexyl Acetate
7 A publication of EMIL ENGINEERING TRANSACTIONS VOL. 6, 07 Guest Editors: Petar S Varbanov, Rongxin Su, Hon Loong Lam, Xia Liu, Jiří J Klemeš Copyright 07, AIDIC Servizi S.r.l. ISBN 978-88-95608-5-8;
More informationAggregate Models based on Improved Group Methods for Simulation and Optimization of Distillation Systems
Aggregate Models based on Improved Group Methods for Simulation and Optimization of Distillation Systems Ravindra S. Kamath 1, Ignacio E. Grossmann 1,*, and Lorenz T. Biegler 1 1 Chemical Engineering Department,
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationDEVELOPMENT OF A CAPE-OPEN 1.0 SOCKET. Introduction
DEVELOPMENT OF A CAPE-OPEN 1. SOCKET Eric Radermecker, Belsim S.A., Belgium Dr. Ulrika Wising, Belsim S.A., Belgium Dr. Marie-Noëlle Dumont, LASSC, Belgium Introduction VALI is an advanced data validation
More informationAN INTRODUCTION TO ASPEN Plus AT CITY COLLEGE
AN INTRODUCTION TO ASPEN Plus AT CITY COLLEGE Prepared by I. H. RINARD Department of Chemical Engineering November 1997 (Revised October 1999) ASPEN Plus provided by Aspen Technology, Inc. Cambridge, Massachusetts
More informationComputer Aided Design of Thermal Systems (ME648)
Computer Aided Design of Thermal Systems (ME648) PG/Open Elective Credits: 3-0-0-9 Updated Syallabus: Introduction. Basic Considerations in Design. Modelling of Thermal Systems. Numerical Modelling and
More informationProcess Design Decisions and Project Economics Prof. Dr. V. S. Moholkar Department of Chemical Engineering Indian Institute of Technology, Guwahati
Process Design Decisions and Project Economics Prof. Dr. V. S. Moholkar Department of Chemical Engineering Indian Institute of Technology, Guwahati Module - 2 Flowsheet Synthesis (Conceptual Design of
More informationMULTIOBJECTIVE OPTIMIZATION CONSIDERING ECONOMICS AND ENVIRONMENTAL IMPACT
MULTIOBJECTIVE OPTIMIZATION CONSIDERING ECONOMICS AND ENVIRONMENTAL IMPACT Young-il Lim, Pascal Floquet, Xavier Joulia* Laboratoire de Génie Chimique (LGC, UMR-CNRS 5503) INPT-ENSIGC, 8 chemin de la loge,
More informationFacing New Products Demand through Simultaneous Structural and Operational Decisions in the Design of the Control Recipe
Facing New Products Demand through Simultaneous Structural and Operational Decisions in the Design of the Control Recipe Marta Moreno-Benito, Antonio Espuña* Chemical Engineering Department, Univesitat
More informationSolving Deterministic Models
Solving Deterministic Models Shanghai Dynare Workshop Sébastien Villemot CEPREMAP October 27, 2013 Sébastien Villemot (CEPREMAP) Solving Deterministic Models October 27, 2013 1 / 42 Introduction Deterministic
More informationSYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES
SYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES by Mark H. Holmes Yuklun Au J. W. Stayman Department of Mathematical Sciences Rensselaer Polytechnic Institute, Troy, NY, 12180 Abstract
More information1111: Linear Algebra I
1111: Linear Algebra I Dr. Vladimir Dotsenko (Vlad) Lecture 6 Dr. Vladimir Dotsenko (Vlad) 1111: Linear Algebra I Lecture 6 1 / 14 Gauss Jordan elimination Last time we discussed bringing matrices to reduced
More informationLecture Notes: Geometric Considerations in Unconstrained Optimization
Lecture Notes: Geometric Considerations in Unconstrained Optimization James T. Allison February 15, 2006 The primary objectives of this lecture on unconstrained optimization are to: Establish connections
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More informationNEURAL NETWORK BASED APPROACH APPLIED TO FOR MODELING AND OPTIMIZATION AN INDUSTRIAL ISOPRENE UNIT PRODUCTION
NEURAL NETWORK BASED APPROACH APPLIED TO FOR MODELING AND OPTIMIZATION AN INDUSTRIAL ISOPRENE UNIT PRODUCTION Rita M. B. Alves and Claudio A. O. Nascimento CESQ- Center of Engineering of Chemical Systems
More informationCryptanalysis of Achterbahn
Cryptanalysis of Achterbahn Thomas Johansson 1, Willi Meier 2, and Frédéric Muller 3 1 Department of Information Technology, Lund University P.O. Box 118, 221 00 Lund, Sweden thomas@it.lth.se 2 FH Aargau,
More informationSurrogate Model Generation Using the Concepts of Self-Optimizing Control
Surrogate Model Generation Using the Concepts of Self-Optimizing Control Julian Straus, Sigurd Skogestad Department of Chemical Engineering Norwegian University of Science and Technology Trondheim, Norway
More informationFitting. PHY 688: Numerical Methods for (Astro)Physics
Fitting Fitting Data We get experimental/observational data as a sequence of times (or positions) and associate values N points: (x i, y i ) Often we have errors in our measurements at each of these values:
More informationEfficient Handling of the Implicit Constraints Problem for the ASPEN MINLP Synthesizer
2006 Ind. Eng. Chem. Res. 1993,32, 2006-2011 Efficient Handling of the Implicit Constraints Problem for the ASPEN MINLP Synthesizer U. M. Diwekar' and E. S. Rubin Environmental Institute h Department of
More informationApplied Computational Fluid Dynamics. in Marine Engineering
Applied Computational Fluid Dynamics in Marine Engineering Objectives Understand basic CFD theory Learn how to set up and run simulations in Star CCM+ and interpret results Learn about limitations and
More informationNew CAPE-OPEN capabilities
ProSimPlus New CAPE-OPEN capabilities #535b CAPE-OPEN Unit Operations: Development and Usage (TD001) Alain VACHER, Stéphane DECHELOTTE, Olivier BAUDOUIN ProSim SA Stratège Bâtiment A BP 2738 31312 LABEGE
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationEffects of Feed Stage on the Product Composition of a Reactive Distillation Process: A Case Study of trans-2-pentene Metathesis
Journal homepage: http://www.journalijar.com INTERNATIONAL JOURNAL OF ADVANCED RESEARCH RESEARCH ARTICLE Effects of Feed Stage on the Product Composition of a Reactive Distillation Process: A Case Study
More informationA Sequential and Hierarchical Approach for the Feasibility Analysis and the Preliminary Synthesis and Design of Reactive Distillation Processes
A Sequential and Hierarchical Approach for the Feasibility Analysis and the Preliminary Synthesis and Design of Reactive Distillation Processes Raphaële Thery, Xuân-Mi Meyer 1, Xavier Joulia Laboratoire
More informationBOUNDARY VALUE DESIGN METHOD FOR COMPLEX DEMETHANIZER COLUMNS
Distillation Absorption 2010 A.B. de Haan, H. Kooijman and A. Górak (Editors) All rights reserved by authors as per DA2010 copyright notice BOUNDARY AUE DESIGN METHOD FOR COMPEX DEMETHANIZER COUMNS Muneeb
More informationAspen Polymers. Conceptual design and optimization of polymerization processes
Aspen Polymers Conceptual design and optimization of polymerization processes Aspen Polymers accelerates new product innovation and enables increased operational productivity for bulk and specialty polymer
More informationTheoretical Models of Chemical Processes
Theoretical Models of Chemical Processes Dr. M. A. A. Shoukat Choudhury 1 Rationale for Dynamic Models 1. Improve understanding of the process 2. Train Plant operating personnel 3. Develop control strategy
More informationDME(10 TPD) Process Simulation Using Aspen Plus Release Dr. Jungho Cho, Professor Department of Chemical Engineering Dong Yang University
DME(10 TPD) Process Simulation Using Aspen Plus Release 12.1 Dr. Jungho Cho, Professor Department of Chemical Engineering Dong Yang University Overall Flowsheet for DME Production Unit 18 TO FLARE 17 DA-103
More informationCONTROL PROPERTIES ANALYSIS OF ALTERNATE SCHEMES TO THERMALLY COUPLED DISTILLATION SCHEMES
8th International IFAC Symposium on Dynamics and Control of Process Systems Preprints Vol.1, June 6-8, 2007, Cancún, Mexico CONTROL PROPERTIES ANALYSIS OF ALTERNATE SCHEMES TO THERMALLY COUPLED DISTILLATION
More informationLecture 8: Determinants I
8-1 MATH 1B03/1ZC3 Winter 2019 Lecture 8: Determinants I Instructor: Dr Rushworth January 29th Determinants via cofactor expansion (from Chapter 2.1 of Anton-Rorres) Matrices encode information. Often
More informationAutomation in Complex Systems MIE090
Automation in Complex Systems MIE090 Exam Monday May 29, 2017 You may bring the course book and the reprints (defined in the course requirements), but not the solution to problems or your own solutions
More informationAzeotropic distillation Example 1
Azeotropic distillation Example 1 Separate water from iso-butanol. The phase behavior for this mixture is interesting. There is a minimum boiling azeotrope formed as well as a liquid-liquid phase separation
More informationNonrobust and Robust Objective Functions
Nonrobust and Robust Objective Functions The objective function of the estimators in the input space is built from the sum of squared Mahalanobis distances (residuals) d 2 i = 1 σ 2(y i y io ) C + y i
More informationMaximum Likelihood Estimation
Maximum Likelihood Estimation Prof. C. F. Jeff Wu ISyE 8813 Section 1 Motivation What is parameter estimation? A modeler proposes a model M(θ) for explaining some observed phenomenon θ are the parameters
More informationRedacted for Privacy
AN ABSTRACT OF THE THESIS OF Nopasit Chaiwanakupt for the degree of Master of Science in Chemical Engineering presented on September 23, 1994 Title: Simulation and Optimization of an Offshore Natural Gas
More informationAmmonia Synthesis with Aspen Plus V8.0
Ammonia Synthesis with Aspen Plus V8.0 Part 1 Open Loop Simulation of Ammonia Synthesis 1. Lesson Objectives Become comfortable and familiar with the Aspen Plus graphical user interface Explore Aspen Plus
More information4 Limit and Continuity of Functions
Module 2 : Limits and Continuity of Functions Lecture 4 : Limit at a point Objectives In this section you will learn the following The sequential concept of limit of a function The definition of the limit
More informationMultiperiod Blend Scheduling Problem
ExxonMobil Multiperiod Blend Scheduling Problem Juan Pablo Ruiz Ignacio E. Grossmann Department of Chemical Engineering Center for Advanced Process Decision-making University Pittsburgh, PA 15213 1 Motivation
More informationCEE 618 Scientific Parallel Computing (Lecture 3)
1 / 36 CEE 618 Scientific Parallel Computing (Lecture 3) Linear Algebra Basics using LAPACK Albert S. Kim Department of Civil and Environmental Engineering University of Hawai i at Manoa 2540 Dole Street,
More informationregressing the vapor-liquid equilibrium data in Mathuni et al. and Rodriguez et al., respectively. The phase equilibrium data of the other missing pai
Plant-Wide Control of Thermally Coupled Reactive Distillation to Co-Produce Diethyl Carbonate and Propylene Glycol San-Jang Wang, Shueh-Hen Cheng, and Pin-Hao Chiu Abstract Diethyl carbonate is recognized
More informationEsterification in a PFR with Aspen Plus V8.0
Esterification in a PFR with Aspen Plus V8.0 1. Lesson Objectives Use Aspen Plus to determine whether a given reaction is technically feasible using a plug flow reactor. 2. Prerequisites Aspen Plus V8.0
More informationFull Furnace Simulations and Optimization with COILSIM1D
Full Furnace Simulations and Optimization with COILSIM1D Kevin M. Van Geem 1 Alexander J. Vervust 1 Ismaël Amghizar 1 Andrés E. Muñoz G. 2 Guy B. Marin 1 1 Laboratory for Chemical Technology, Ghent University
More informationCBE495 LECTURE IV MODEL PREDICTIVE CONTROL
What is Model Predictive Control (MPC)? CBE495 LECTURE IV MODEL PREDICTIVE CONTROL Professor Dae Ryook Yang Fall 2013 Dept. of Chemical and Biological Engineering Korea University * Some parts are from
More informationNOTICE WARNING CONCERNING COPYRIGHT RESTRICTIONS: The copyright law of the United States (title 17, U.S. Code) governs the making of photocopies or
NOTICE WARNING CONCERNING COPYRIGHT RESTRICTIONS: The copyright law of the United States (title 17, U.S. Code) governs the making of photocopies or other reproductions of copyrighted material. Any copying
More informationLinear Algebra I Lecture 10
Linear Algebra I Lecture 10 Xi Chen 1 1 University of Alberta January 30, 2019 Outline 1 Gauss-Jordan Algorithm ] Let A = [a ij m n be an m n matrix. To reduce A to a reduced row echelon form using elementary
More informationDESIGN AND CONTROL OF BUTYL ACRYLATE REACTIVE DISTILLATION COLUMN SYSTEM. I-Lung Chien and Kai-Luen Zeng
DESIGN AND CONTROL OF BUTYL ACRYLATE REACTIVE DISTILLATION COLUMN SYSTEM I-Lung Chien and Kai-Luen Zeng Department of Chemical Engineering, National Taiwan University of Science and Technology, Taipei
More informationANALYSIS OF LOAD PATTERNS IN RUBBER COMPONENTS FOR VEHICLES
ANALYSIS OF LOAD PATTERNS IN RUBBER COMPONENTS FOR VEHICLES Jerome Merel formerly at Hutchinson Corporate Research Center Israël Wander Apex Technologies Pierangelo Masarati, Marco Morandini Dipartimento
More informationA Robust Controller for Scalar Autonomous Optimal Control Problems
A Robust Controller for Scalar Autonomous Optimal Control Problems S. H. Lam 1 Department of Mechanical and Aerospace Engineering Princeton University, Princeton, NJ 08544 lam@princeton.edu Abstract Is
More informationAn IDL Based Image Deconvolution Software Package
An IDL Based Image Deconvolution Software Package F. Városi and W. B. Landsman Hughes STX Co., Code 685, NASA/GSFC, Greenbelt, MD 20771 Abstract. Using the Interactive Data Language (IDL), we have implemented
More informationReactors. Reaction Classifications
Reactors Reactions are usually the heart of the chemical processes in which relatively cheap raw materials are converted to more economically favorable products. In other cases, reactions play essential
More information