Domain s Sweep, a Computational Method to Optimise Chemical and Physical Processes
|
|
- Melvin Hopkins
- 5 years ago
- Views:
Transcription
1 Article Domain s Sweep, a Computational Method to Optimise Chemical and Physical Processes Alexandre César Balbino Barbosa Filho Department of Chemical Engineering, Federal University of Campina Grande, Campina Grande, PB 882, BRAZIL; Acbbarbosaf2@gmail.com or alexandrecb.barbosa@eq.ufcg.edu.br; Tel.: +55(82) Abstract: Many of engineering s problems are about to optimise equations, functions and process models equations which appear on common or even complex science s cases. Most applied optimisation methods nowadays can only be used on particular cases, so the main objective of this article is to define a computational method that can optimise or even find target values for specified objective functions and variables, by just using computer effort. Chemical plants involve many differential equations and they can be optimised with more facility now. On this paper, there are shown as example, an optimisation for a chemical reactor, which is found the optimal temperature and volumetric flow rate of feed, for the given objective function being the composition of the desired product. Named Domain s Sweep, it is an algorithm that evaluate given(s) mathematical function(s) and/or equation(s) by varying your(s) independent variable(s) through loops with a given step size and solving it after closing the degrees of freedom, and finally, with some condition statements, store all the optimum values of given or created objective functions with its respective independent variables. In another words, the user create an objective function and this method find the function s maximum, minimum or a certain chosen target value, even if it does not have an inflection point in the given search interval of the independent variables. Keywords: computational methods; optimisation; plant optimisation; chemical processes; physical processes; computational intelligence; maximum value; minimum value; functions; 1. Introduction Numerical methods for optimisation of equations usually uses mathematical fundamentals and theories over each specific situation and on most of them, it is necessary a given stop criteria and these methods are just used to find the maximum and the minimum value of a function. Sometimes at some cases, specially the complex ones, these methods do not converge because of the given initial estimate or even the given search interval, or because it does not have a maximum or minimum in the function that is being optimised at a certain range of independent variables values. It is commonly used numerical methods to optimise because most of the problems do not have an analytical solution, like partial differential equations of 3 independent variables and some other non-linear equations that a linearization does not satisfy the specific case because the approximation would not give a true value [1]. So the computational intelligence on some cases is a good tool to be used to solve complex problems and to facilitate its solution. The computational intelligence is evolving over the years, and it is functionally worth to deal with it on complex tasks that human being could not do by their own hands because of the hard work and the resolution time. With the developing of the CI and the search engineering, many discovered solutions are beyond the human engineer imagination [2]. CI is used to do a lot of tasks, but in some cases it is not preferred because of the difficulty of its architecture, but with the pass of the years it is 2017 by the author(s). Distributed under a Creative Commons CC BY license.
2 2 of 12 changing. One of its uses is a linear search algorithm that seeks a collection and a target item and finds whether the target is in the collection, the run time is linear and the target is found by an organized order of searching through the items of the collection, stopping when it finds the target [3]. A several number of CI algorithms have been developed and published during the last years but as they were not unified, complicated the excellent use of this computational science [4]. With the CI, human being tasks (complex or even simple tasks) had been solved by a machine with the implementation of a computing algorithm on the machine`s software, so, the CI theories show interests in knowing how the problems (specified by human being) can be represented and solved by a machine [5]. On many years, chemical and physical processes were optimised based on objective functions, which were linked to the financial business. And instead of finding the operational conditions of a chemical plant s equipment based on its best performance (for example, the maximum conversion of a chemical reactor or the complete separation of a mixture on a column distillation), plants were constructed with clearance. A degree of freedom (DOF) is the main analysis that must be done before solving any equation, as it is already known, an equation can only be solved if the DOF is equal to zero, which is the number of equations subtracted by the number of independent variables of each equation. The main idea for the algorithm presented in this paper is to close the degree of freedom of the user equation(s) by specifying values for its independent variables at order, through closed loops, and the objective of this paper is to show a computational method that find the optimum or desired value for chosen independents variables in equations, as users will, that flees the divergence problem and lead to the optimum stage of chemical and physical processes. 2. Methodology Suppose it is required to find the amount of N optimum independent variables values that leads to the maximum, minimum or even a certain specified target value of a function (must be defined by the user) on a specified interval of the independent function s variables. It s only needed a computational algorithm that does the following: I) Specify values to store functions (its value is updated according the hike of the sweep) that will be compared into those ones found by III or IV, on posterior items V, VI and VII. A minimum store function value must have an initial estimated value of a high positive number, because thus the interval of possible values will be large and no problem will occur at storing the first minimum value. On about storing a maximum store function value it is needed a high negative number as initial estimated value. The target store function value is determined by the user. II) Sweep the determined function(s) or equation(s) domain with N loops, where N is the number of function s independent variable. A sweep is an expression to refer that a loop is being done to vary a function s independent variable to evaluate its value. The sweep is always in the direction of minus infinite to positive infinite. Each loop varies a different independent variable starting at an initial chosen value to a final chosen value with a given step size. The loops must be done one inside another. III) Solve equation(s) or function(s) for independent variables in the last loop; IV) Still in the last loop, define new functions if wanted, that have the same independent variables which values were found in item III;
3 3 of 12 V) Then in the same loop, it is put a condition statement that can express that if the current calculated function s value is higher than the current maximum store function value, then, store all the independent variables value and update the maximum store function value to the current value of the calculated function. At the end of all loops, it will have stored all the optimum independent variables values that lead to the maximum value of the defined function by the user. VI) After closing the first condition statement in last loop, it is put another one that can express that if the current calculated function s value is lower than the current minimum store function value, then, store all the independent variables value and update the minimum store function value to the current value of the calculated function. At the end of all loops, it will have stored all the optimum independent variables values that lead to the minimum value of the defined function by the user. VII) Before closing the last loop, the last condition statement may do that if the currently calculated function s value is between the limits of and of the specified target store function value, then, store all the function s and independent variable s value. At the end of all loops, it will have stored the independent variables values that lead to the specified target function s value with the given tolerance limits, and it is common on a function of more than one variable to have several points that gives the desired target function value, so, it must store all these values in a vector. About the tolerance limits, it can be changed with the will of the user. On some cases it is optionally to the user to put a condition statement in any of other loops (not necessary in the last loop), but remembering that it is not occasionally done in optimisation, because the need of sweep all the function s domain. User may choose which conditions statements will be used (V, VI and VII) Algorithm s schema Since Vi is any independent variable, a step scheme algorithm of the method for a function of N independent variables should be: Step 1: Function target value = input value; Step 2: Maximum store function = ; Step 3: Minimum store function = ; Step 4: Loop 1: Vary V1 with a given step to a given maximum range Step 5: Loop 2: Vary V2 with a given step to a given maximum range And so on till loop N Step 6: Loop N: Vary VN with a given step to a given maximum range Step 7: Solve equation(s) or function(s) declared by the user and define new functions (if wanted), based on the independent variables; Step 8: If 1: (Function s value) > Maximum store function Then: Maximum store function = (Function s value); V1 max = V1; V2 max = V2; VN max = VN. End If 1 Step 9: If 2: (Function s value) < Minimum store function
4 4 of 12 Then: Minimum store function = (Function s value); V1 min = V1; V2 min = V2; VN min = VN. End If 2 Step10: If 3: (Function s value) >= 0.999x (Function target value) and <=1.001x (Function target value) Then: Store the function and independent variables values in a vector; V1 target = V1; V2 target = V2; VN target = VN. End If 3; End loop 1; End loop 2; End loop N. When the algorithm s run finishes, all the function(s) maximum, minimum and even a target value are stored with its respective independent variables. Remember that solutions found on step 8 must converge, so then computer will not enter on an infinite loop. Because of that, it can be optionally put another condition statement to warn if the problem is not converging or the function(s) value diverge, or even to ignore that and keep running. 3. Real-life Impact and applications 3.1. Chemical reactor optimisation This method is very useful in day-to-day on processes that can be described by algebraic equations. One application for example, is to find the optimum operational conditions in a chemical reactor by solving the process model equations for independent variables and creating an objective function, like the conversion of a specific reactant or even the exit s composition of a desired product and so on. Now, it is shown the optimisation of a continuous stirred tank reactor (CSTR), which the objective was to maximize the composition of the desired product (B) at the reactor s exit [6-9]. The example was performed by MATLAB (R2013a, Mathworks, Natick, MA, USA): Reactions occurring: (1) 1 (2) 2 Energy balance for reactor:
5 5 of 12 = + ( + ).. (3) Energy balance for coolant jacket: = +.. (4) Component molar balances: = (5) = + (6) = + (7) = + 2 (8) Where, =... (9) =.. ( ) (10) =. ( ). (11) =. ( ). (12) =. ( ). (13) The table below shows the meaning of each variable used on this application: Table 1. Variables meaning. Variable MW Meaning Fluid density inside the reactor Molecular weight Feed concentration of A Feed concentration of B Table 1. Cont. Variable Meaning Feed concentration of C
6 6 of 12 Feed concentration of D Enthalpy of the reaction 1 Enthalpy of the reaction 2 Enthalpy of the reaction 3 Inside reactor s fluid heat capacity Feed Coolant s jacket fluid temperature Coolant s jacket fluid density Coolant s jacket fluid heat capacity Pre-exponential factor of reaction 1 Pre-exponential factor of reaction 2 Pre-exponential factor of reaction 3 Activation energy of reaction 1 Activation energy of reaction 2 Activation energy of reaction 3 Coolant s jacket fluid volumetric flow rate D Reactor s diameter R Reactor s radius h Reactor s fluid level U Global heat transfer coefficient esp Coolant s jacket thickness Reactor s lateral area Reactor s volume Coolant s jacket volume Coolant s jacket residence time Feed concentration of A Reactor s residence time T Reactor s temperature Coolant s jacket temperature Molar concentration of A Molar concentration of B Molar concentration of C Molar concentration of D Composition of B at Reactor s exit Feed temperature Reactor s fluid volumetric flow rate Then it was specified fixed parameters of the process model equations: Table 2. Parameters for the simulation. Parameter Value Unit MW Table 2. Cont. Parameter Value Unit
7 7 of K x10 12 min x10 12 min x K K 8560 K 30. D 2 m R 1 m h 4 m U esp 0.25 m Some other terms appears in the balance equations, so it must be calculated before the simulation. Assuming an approximately cylinder volume for reactor and that the coolant jacket includes the top, base and the side of the reactor and has a specified thickness of refrigerant fluid, it follows that the geometric calculus arrive to: =2..h. = (14) =..h = (15) Then, the volume of the coolant jacket is represented by a cylindrical shell with a specified thickness: =.h..(2. + ) = (16) The refrigerant fluid used was the DOWTHERM-Q and your residence time is: = = (17) As shown in Table 1, it was adopted a value of 40. to the molecular weight of the reactant A, and a medium value of for the reactor s fluid density. Then the feed concentration value of A is given by: = =25 (18) After declared all fixed parameters and store functions, it was used the domain s sweep algorithm to architecture the optimisation of the process. It was varied through the loops, the feed temperature and the feed volumetric flow of the reactant A on an interval of 307 to 407 and 0.1 to 100, respectively, with a step size of 0.1 for both independent variables. But to solve the balance
8 8 of 12 equations it needs the value of the reactant A residence time inside the reactor, so it was written in the last loop of the code its formula: = (19) Note that the value of changes with the change of the feed volumetric flow of reactant A, so its value is updated at each iteration. Also in the last loop of the algorithm s code it was solved the balance equations using Runge-Kutta s method for ordinary differential equations (ODE), declared conditions statements to store the optimum values over the store functions and created the objective function of interest (The composition of product B). As the example was performed by MATLAB, the command to do the loop, the conditional statements and to solve the system of ODES were for and if respectively, according to its programming language. The simulation (Solving the balance equations) was done in the last for and it was done by setting initial conditions, based on feed properties: Table 3. Initial conditions for the simulation (Same values of Table 2). Dependent variable T Initial condition ( it is being varied) As follows the Domain s Sweep algorithm and with the results in hands, the objective function created was the composition of product B that is function of the steady-state value (SSV) of the concentration of reactant A, product B, C and D found in the simulation of the balance equations. So the simulation was run with a big considerably time so that the values could arrive the steady-state. After the MATLAB code was run the maximum value of the composition of product B (objective function) and its feed temperature and feed volumetric flow of the reactant A (operational conditions) attached were: = = (20) = (21) = 100 (22) Remembering that since the operational conditions are the values, which have the maximum value for the objective function, on this case, the store function created before the loops was a maximum store function that has its values updated by doing step 8 of algorithm s schema. Also in the same step, the independent variables (, ) were being stored in a store function when the conditional statement was true.
9 9 of 12 The results of the simulation of the process model are shown. First it is shown an analysis of sensibility of the composition of product B with a fixed optimum value of K for the feed temperature: Figure 1. Analysis of sensitivity: versus to show the consistence of the results. As the found optimum values that give the maximum value for the composition of B were a feed temperature of K and a feed volumetric flow of reactant A of 100 m 3, then, if it is plotted a graph of XB versus Fe, the maximum value must be at the point (100, ) and as it can be seen in Figure 1, the composition of B has a maximum value at the found optimum operational conditions ( = and =100). But it is not enough at all, on a chemical process just to do the optimisation, so it is needed a proof that the process is real and stable, so since the simulation was done when it was also doing the optimisation, it is also plotted below the dynamic plot of the process variables of the system. Figure 2 and 3 show that the temperature and the concentration of the reactants and products remains with real values in the whole simulation and get a constant and real value as they reach the steady-state, proofing, that the system is real, stable and reachable. Figure 2. Molar concentrations of reactants and products.
10 10 of 12 Figure 3. Reactor s and coolant jacket s temperature. The system is stable as shown on the graphs, and with the found optimum operational conditions, the steady-state values for the reactor s temperature, coolant jacket s temperature and the exit molar concentration of A, B, C and D with their given units are , , 7.636, 3.912, and 4.672, respectively. Also, a 3D-plot and a contour graph were done to enhance the reader s understanding, as follow in Figure 4 and 5: Figure 4. 3D-plot of the composition of product B, as a function of feed temperature and feed volumetric flow rate with a blue circle on the maximum value.
11 11 of 12 Figure 5. Contour graph showing the composition values of product B in the colorbar. 4. Conclusions Of easy understanding and application, it is a method to find the optimum value as the user choices, remembering that according to the user s code, it can find the maximum, minimum or a determined function value and the accuracy of the found value compared to the real value is linked to the step s size criteria on about varying the independent variables of the function(s). Smaller the step size, higher the accuracy of the optimum value and higher the machine processing time. In respect of time, the Domain s Sweep algorithm s delay time to finish is proportional to step s size criteria, specified sweep s range and the number of functions independent variables. Conflicts of Interest: The author declares no conflict of interest. References 1. Chapra, Steven C.; Canale, Raymond P. Numerical Methods for engineers, 6th ed.; McGraw-Hill: 1221 Avenue of the Americas, New York, NY 10020, USA, 2010; pp Downey, Allen B. Think Complexity, version 1.2.3; Green Tea Press: 9 Washburn Ave, Needham 02492, USA, 2012; pp Downey, Allen B. Think Complexity, version 1.2.3; Green Tea Press: 9 Washburn Ave, Needham 02492, USA, 2012; pp Xing, B.; Gao, W.J. Innovative Computational Intelligence: A Rough Guide to 134 Clever Algorithms, volume 62; Springer International Publishing: Switzerland, 2014; pp Konar, A. Computational Intelligence: Principles, Techniques and Applications, 1st ed.; Springer International Publishing: Switzerland, 2005; p Incropera, Frank P.; DeWitt, David P.; Bergman, Theodore L.; Lavine, Adrienne S. Fundamentals of Heat and Mass Transfer, 6th ed.; John Wiley & Sons, Inc.: New Jersey, USA, 2007; pp Seborg, Dale E.; Edgar, Thomas F.; Mellichamp, Duncan A.; Doyle III, Francis J. Process Dynamics and Control, 3rd ed.; John Wiley & Sons, Inc.: New Jersey, USA, 2011; pp
12 12 of Magalhães, Otto I. B. Desenvolvimento de um sistema de otimização dinâmica em tempo real. Master s Program, UFRJ, Rio De Janeiro, Brazil, 2010; pp Smith, J. M.; Van Ness, H. C.; Abbott, M. M. Introduction to Chemical Engineering Thermodynamics, 6th ed.; McGraw-Hill Education: New York, NY, 2001; pp
Use 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
Principles of Instability and Stability in Digital PID Control Strategies and Analysis for a Continuous Alcoholic Fermentation Tank Process Start-up
Principles of Instability and Stability in Digital PID Control Strategies and Analysis for a Continuous Alcoholic Fermentation Tank Process Start-up Alexandre César Balbino Barbosa Filho a a Faculty of
CHEMICAL REACTORS - PROBLEMS OF REACTOR ASSOCIATION 47-60
2011-2012 Course CHEMICL RECTORS - PROBLEMS OF RECTOR SSOCITION 47-60 47.- (exam jan 09) The elementary chemical reaction in liquid phase + B C is carried out in two equal sized CSTR connected in series.
Chemical Reaction Engineering Prof. JayantModak Department of Chemical Engineering Indian Institute of Science, Bangalore
Chemical Reaction Engineering Prof. JayantModak Department of Chemical Engineering Indian Institute of Science, Bangalore Module No. #05 Lecture No. #29 Non Isothermal Reactor Operation Let us continue
USE OF MATLAB TO UNDERSTAND BASIC MATHEMATICS
USE OF MATLAB TO UNDERSTAND BASIC MATHEMATICS Sanjay Gupta P. G. Department of Mathematics, Dev Samaj College For Women, Punjab ( India ) ABSTRACT In this paper, we talk about the ways in which computer
BASIC RESEARCH OF THERMAL TRANSFER SIMULATIONS
27TH DAAAM INTERNATIONAL SYMPOSIUM ON INTELLIGENT MANUFACTURING AND AUTOMATION DOI:.257/27th.daaam.proceedings.85 BASIC RESEARCH OF THERMAL TRANSFER SIMULATIONS Václav Marek This Publication has to be
A First Course on Kinetics and Reaction Engineering Unit 2. Reaction Thermochemistry
Unit 2. Reaction Thermochemistry Overview This course is divided into four parts, I through IV. Part I reviews some topics related to chemical reactions that most students will have encountered in previous
Model of Mass and Energy Transfer in a Clinker Rotary Kiln
Model of Mass and Energy Transfer in a Clinker Rotary Kiln J.A. Guirao S. Iglesias Numerical Analysis TEChnologies S.L. J. Pistono Universidad de Oviedo R. Vázquez IMASA Div. Internacional Abstract The
Development of Dynamic Models. Chapter 2. Illustrative Example: A Blending Process
Development of Dynamic Models Illustrative Example: A Blending Process An unsteady-state mass balance for the blending system: rate of accumulation rate of rate of = of mass in the tank mass in mass out
REFERENCE TEXTBOOKS MATERIAL AND ENERGY BALANCE: THERMODYNAMICS: TRANSPORT: KINETICS AND REACTOR DESIGN:
REFERENCE TEXTBOOKS MATERIAL AND ENERGY BALANCE: Elementary Principles of Chemical Processes, 4 th Edition - Richard M. Felder, Ronald W. Rousseau, Lisa G. Bullard THERMODYNAMICS: Introduction to Chemical
Review Paper Modeling and Simulation Study of the CSTR for Complex Reaction by Using Polymath
Research Journal of Chemical Sciences ISSN 2231-606X. Review Paper Modeling and Simulation Study of the CSTR for Complex Reaction by Using Polymath Abstract Kanse Nitin G. 1, Dhanke P.B. 1 and Thombare
(Refer Slide Time: 00:10)
Chemical Reaction Engineering 1 (Homogeneous Reactors) Professor R. Krishnaiah Department of Chemical Engineering Indian Institute of Technology Madras Lecture No 10 Design of Batch Reactors Part 1 (Refer
INTRODUCTION 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
Reactor 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
Chapter 5. Formulation of FEM for Unsteady Problems
Chapter 5 Formulation of FEM for Unsteady Problems Two alternatives for formulating time dependent problems are called coupled space-time formulation and semi-discrete formulation. The first one treats
Design of de-coupler for an interacting tanks system
IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 78-1676,p-ISSN: 3-3331, Volume 7, Issue 4 (Sep. - Oct. 13), PP 48-53 Design of de-coupler for an interacting tanks system Parag
Numerical Solution of Differential Equations
1 Numerical Solution of Differential Equations A differential equation (or "DE") contains derivatives or differentials. In a differential equation the unknown is a function, and the differential equation
A First Course on Kinetics and Reaction Engineering Example 11.5
Example 11.5 Problem Purpose This problem illustrates the use of the age function measured using a step change stimulus to test whether a reactor conforms to the assumptions of the ideal PFR model. Then
Example 8: CSTR with Multiple Solutions
Example 8: CSTR with Multiple Solutions This model studies the multiple steady-states of exothermic reactions. The example is from Parulekar (27) which in turn was modified from one by Fogler (1999). The
CONSIM - MS EXCEL BASED STUDENT FRIENDLY SIMULATOR FOR TEACHING PROCESS CONTROL THEORY
CONSIM - MS EXCEL BASED STUDENT FRIENDLY SIMULATOR FOR TEACHING PROCESS CONTROL THEORY S. Lakshminarayanan, 1 Rao Raghuraj K 1 and S. Balaji 1 1 Department of Chemical and Biomolecular Engineering, 4 Engineering
Chemical Engineering 140. Chemical Process Analysis C.J. Radke Tentative Schedule Fall 2013
Chemical Process Analysis C.J. Radke Tentative Schedule Fall 2013 Week 0 *8/30 1. Definition of Chemical Engineering: flow sheet, reactor trains and separation processes, raw materials, power production
Subject Guide. Facultad de Química
Bachelor s Degree in Chemistry Subject Guide 1. Information about the subject SUBJECT Chemical Engineering CODE GQUIMI01-2-008 EDUCATIONAL OFFER Bachelor s Degree in Chemistry CENTER Facultad de Química
CHAPTER 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
Module 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
International Journal of Research in Advent Technology Available Online at:
A COMPUTER PROGRAMMED DESIGN OPTIMISATION AND ANALYSIS OF COMPRESSOR IMPELLER G. Naga Malleshwar Rao 1, Dr. S.L.V. Prasad 2, Dr. S. Sudhakarbabu 3 1, 2 Professor of Mechanical Engineering, Shri Shirdi
Theoretical 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
Numerical Analysis of Plate Heat Exchanger Performance in Co-Current Fluid Flow Configuration
Numerical Analysis of Plate Heat Exchanger Performance in Co-Current Fluid Flow Configuration H. Dardour, S. Mazouz, and A. Bellagi Abstract For many industrial applications plate heat exchangers are demonstrating
Chabot College Fall Course Outline for Mathematics 25 COMPUTATIONAL METHODS FOR ENGINEERS AND SCIENTISTS
Chabot College Fall 2010 Course Outline for Mathematics 25 COMPUTATIONAL METHODS FOR ENGINEERS AND SCIENTISTS Catalog Description: MTH 25 - Computational Methods for Engineers and Scientists 3.00 units
CHE 404 Chemical Reaction Engineering. Chapter 8 Steady-State Nonisothermal Reactor Design
Textbook: Elements of Chemical Reaction Engineering, 4 th Edition 1 CHE 404 Chemical Reaction Engineering Chapter 8 Steady-State Nonisothermal Reactor Design Contents 2 PART 1. Steady-State Energy Balance
Name. Honor Code: I have neither given nor received unauthorized aid on this examination, nor have I concealed any violations of the Honor Code.
ChE 344 Fall 014 Mid Term Exam II Wednesday, November 19, 014 Open Book Closed Notes (but one 3x5 note card), Closed Computer, Web, Home Problems and In-class Problems Name Honor Code: I have neither given
CHEMICAL ENGINEEERING AND CHEMICAL PROCESS TECHNOLOGY Vol. III - Ideal Models Of Reactors - A. Burghardt
IDEAL MODELS OF REACTORS A. Institute of Chemical Engineering, Polish Academy of Sciences, Poland Keywords: Thermodynamic state, conversion degree, extent of reaction, classification of chemical reactors,
Chemical Reaction Engineering - Part 14 - intro to CSTRs Richard K. Herz,
Chemical Reaction Engineering - Part 4 - intro to CSTRs Richard K. Herz, rherz@ucsd.edu, www.reactorlab.net Continuous Stirred Tank Reactors - CSTRs Here are a couple screenshots from the ReactorLab, Division
Fundamentals of Heat Transfer Muhammad Rashid Usman
Fundamentals of Heat Transfer Muhammad Rashid Usman Institute of Chemical Engineering and Technology University of the Punjab, Lahore. Figure taken from: http://heatexchanger-design.com/2011/10/06/heat-exchangers-6/
Effect Analysis of Volume Fraction of Nanofluid Al2O3-Water on Natural Convection Heat Transfer Coefficient in Small Modular Reactor
World Journal of Nuclear Science and Technology, 2016, 6, 79-88 Published Online January 2016 in SciRes. http://www.scirp.org/journal/wjnst http://dx.doi.org/10.4236/wjnst.2016.61008 Effect Analysis of
WILEY. Differential Equations with MATLAB (Third Edition) Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg
Differential Equations with MATLAB (Third Edition) Updated for MATLAB 2011b (7.13), Simulink 7.8, and Symbolic Math Toolbox 5.7 Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg All
Development of Dynamic Models. Chapter 2. Illustrative Example: A Blending Process
Development of Dynamic Models Illustrative Example: A Blending Process An unsteady-state mass balance for the blending system: rate of accumulation rate of rate of = of mass in the tank mass in mass out
Modeling and Control of Chemical Reactor Using Model Reference Adaptive Control
Modeling and Control of Chemical Reactor Using Model Reference Adaptive Control Padmayoga.R, Shanthi.M 2, Yuvapriya.T 3 PG student, Dept. of Electronics and Instrumentation, Valliammai Engineering College,
INDUSTRIAL INSULATION for Systems Operating Above Ambient Temperature
INDUSTRIAL INSULATION for Systems Operating Above Ambient Temperature U.S. Department of Energy Energy Efficiency and Renewable Energy Office of Industrial Technologies Washington, D.C. 20585 From September,
Advanced Chemical Reaction Engineering Prof. H. S. Shankar Department of Chemical Engineering IIT Bombay. Lecture - 03 Design Equations-1
(Refer Slide Time: 00:19) Advanced Chemical Reaction Engineering Prof. H. S. Shankar Department of Chemical Engineering IIT Bombay Lecture - 03 Design Equations-1 We are looking at advanced reaction engineering;
H-Infinity Controller Design for a Continuous Stirred Tank Reactor
International Journal of Electronic and Electrical Engineering. ISSN 974-2174 Volume 7, Number 8 (214), pp. 767-772 International Research Publication House http://www.irphouse.com H-Infinity Controller
CHAPTER 2 CONTINUOUS STIRRED TANK REACTOR PROCESS DESCRIPTION
11 CHAPTER 2 CONTINUOUS STIRRED TANK REACTOR PROCESS DESCRIPTION 2.1 INTRODUCTION This chapter deals with the process description and analysis of CSTR. The process inputs, states and outputs are identified
Recovering from a Gradual Degradation in MPC Performance
Preprints of the The International Federation of Automatic Control Furama Riverfront, Singapore, July 1-13, 212 Recovering from a Gradual Degradation in MPC Performance Mohammed Jimoh*, John Howell** *School
A First Course on Kinetics and Reaction Engineering Example 13.2
Example 13. Problem Purpose This example illustrates the analysis of kinetics data from a CSTR where the rate expression must be linearized. Problem Statement A new enzyme has been found for the dehydration
Solutions for Tutorial 4 Modelling of Non-Linear Systems
Solutions for Tutorial 4 Modelling of Non-Linear Systems 4.1 Isothermal CSTR: The chemical reactor shown in textbook igure 3.1 and repeated in the following is considered in this question. The reaction
What are Numerical Methods? (1/3)
What are Numerical Methods? (1/3) Numerical methods are techniques by which mathematical problems are formulated so that they can be solved by arithmetic and logic operations Because computers excel at
THE METHOD OF THE WORKING FLUID SELECTION FOR ORGANIC RANKINE CYCLE (ORC) SYSTEM WITH VOLUMETRIC EXPANDER. * Corresponding Author ABSTRACT
Paper ID: 79, Page 1 THE METHOD OF THE WORKING FLUID SELECTION FOR ORGANIC RANKINE CYCLE (ORC) SYSTEM WITH VOLUMETRIC EXPANDER Piotr Kolasiński* 1 1 Wrocław University of Technology, Department of Thermodynamics,
Chapter 2 Solutions of Equations of One Variable
Chapter 2 Solutions of Equations of One Variable 2.1 Bisection Method In this chapter we consider one of the most basic problems of numerical approximation, the root-finding problem. This process involves
Solutions for Tutorial 3 Modelling of Dynamic Systems
Solutions for Tutorial 3 Modelling of Dynamic Systems 3.1 Mixer: Dynamic model of a CSTR is derived in textbook Example 3.1. From the model, we know that the outlet concentration of, C, can be affected
Acoustic radiation by means of an acoustic dynamic stiffness matrix in spherical coordinates
Acoustic radiation by means of an acoustic dynamic stiffness matrix in spherical coordinates Kauê Werner and Júlio A. Cordioli. Department of Mechanical Engineering Federal University of Santa Catarina
Delay Differential Equations with Constant Lags
Delay Differential Equations with Constant Lags L.F. Shampine Mathematics Department Southern Methodist University Dallas, TX 75275 shampine@smu.edu S. Thompson Department of Mathematics & Statistics Radford
Model Predictive Control Design for Nonlinear Process Control Reactor Case Study: CSTR (Continuous Stirred Tank Reactor)
IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331, Volume 7, Issue 1 (Jul. - Aug. 2013), PP 88-94 Model Predictive Control Design for Nonlinear Process
Basic Analysis of Data
Basic Analysis of Data Department of Chemical Engineering Prof. Geoff Silcox Fall 008 1.0 Reporting the Uncertainty in a Measured Quantity At the request of your supervisor, you have ventured out into
A First Course on Kinetics and Reaction Engineering Example 26.3
Example 26.3 unit. Problem Purpose This problem will help you determine whether you have mastered the learning objectives for this Problem Statement A perfectly insulated tubular reactor with a diameter
TUBE BANKS TEST PROBLEMS
TUBE BANKS TEST PROBLEMS The non-proprietary tests used to validate INSTED analysis of flow and heat transfer over tube banks are presented in this section. You may need to consult the original sources
TEACH 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,
Appendix A Course Syllabi Appendix A: Syllabi. Engineering Physics. Bachelor of Science in Engineering Physics. Self-Study Report
Appendix A Course Syllabi Appendix A: Syllabi Engineering Physics Bachelor of Science in Engineering Physics Self-Study Report New Mexico State University 159 Chemical Engineering Courses Chemical Engineering
Automation 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
Introduction to the course ``Theory and Development of Reactive Systems'' (Chemical Reaction Engineering - I)
Introduction to the course ``Theory and Development of Reactive Systems'' (Chemical Reaction Engineering - I) Prof. Gabriele Pannocchia Department of Civil and Industrial Engineering (DICI) University
BERYLLIUM IMPREGNATION OF URANIUM FUEL: THERMAL MODELING OF CYLINDRICAL OBJECTS FOR EFFICIENCY EVALUATION
BERYLLIUM IMPREGNATION OF URANIUM FUEL: THERMAL MODELING OF CYLINDRICAL OBJECTS FOR EFFICIENCY EVALUATION A Senior Scholars Thesis by NICHOLAS MORGAN LYNN Submitted to the Office of Undergraduate Research
Electrochemical Process Engineering. A Guide to the Design of Electrolytic Plant
Electrochemical Process Engineering A Guide to the Design of Electrolytic Plant Electrochemical Process Engineering A Guide to the Design of Electrolytic Plant F. Goodridge and K. Scott University of Newcastle
Fundamentals of Heat Transfer Muhammad Rashid Usman
Fundamentals of Heat Transfer Muhammad Rashid Usman Institute of Chemical Engineering and Technology University of the Punjab, Lahore. Figure taken from: http://heatexchanger-design.com/2011/10/06/heat-exchangers-6/
CHAPTER 5 Batch and Semi-Batch Operations
CHAPTER 5 Batch and Semi-Batch Operations PROCESS CYCLE TIME We must know the process cycle time for a batch or semi-batch process in order to estimate the number of reactors required for a given plant
reality is complex process
ISS0080 Automation and Process Control Lecture 5 1 Process models the desire to describe reality Model of the process, model simplication, identication. model reality is complex process Replaces the original;
6. Multiple Reactions
6. Multiple Reactions o Selectivity and Yield o Reactions in Series - To give maximum selectivity o Algorithm for Multiple Reactions o Applications of Algorithm o Multiple Reactions-Gas Phase 0. Types
Appendix F. + 1 Ma 1. 2 Ma Ma Ma ln + K = 0 (4-173)
5:39p.m. Page:949 Trimsize:8.5in 11in Appendix F F.1 MICROSOFT EXCEL SOLVER FOR NON-LINEAR EQUATIONS The Solver is an optimization package that finds a maximum, minimum, or specified value of a target
MODELING OF CONTINUOUS OSCILLATORY BAFFLED REACTOR FOR BIODIESEL PRODUCTION FROM JATROPHA OIL ABSTRACT
MODELING OF CONTINUOUS OSCILLATORY BAFFLED REACTOR FOR BIODIESEL PRODUCTION FROM JATROPHA OIL B. K. Highina, I. M. Bugaje & B. Gutti Department of Chemical Engineering University of Maiduguri, Borno State,
Extending the Patched-Conic Approximation to the Restricted Four-Body Problem
Monografías de la Real Academia de Ciencias de Zaragoza 3, 133 146, (6). Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Thomas R. Reppert Department of Aerospace and Ocean
BITS-Pilani Dubai, International Academic City, Dubai Second Semester. Academic Year
BITS-Pilani Dubai, International Academic City, Dubai Second Semester. Academic Year 2007-2008 Evaluation Com anent: Com rehensive Examination Closed Book CHE UC441/11NSTR UC 45'1 PROCESS CONTROL Date:
Axial profiles of heat transfer coefficients in a liquid film evaporator
Axial profiles of heat transfer coefficients in a liquid film evaporator Pavel Timár, Ján Stopka, Vladimír Báleš Department of Chemical and Biochemical Engineering, Faculty of Chemical and Food Technology,
Getting started with BatchReactor Example : Simulation of the Chlorotoluene chlorination
Getting started with BatchReactor Example : Simulation of the Chlorotoluene chlorination 2011 ProSim S.A. All rights reserved. Introduction This document presents the different steps to follow in order
Introduction 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
A First Course on Kinetics and Reaction Engineering Unit 33. Axial Dispersion Model
Unit 33. Axial Dispersion Model Overview In the plug flow reactor model, concentration only varies in the axial direction, and the sole causes of that variation are convection and reaction. Unit 33 describes
Mathematics with Maple
Mathematics with Maple A Comprehensive E-Book Harald Pleym Preface The main objective of these Maple worksheets, organized for use with all Maple versions from Maple 14, is to show how the computer algebra
Fault Location in Distribution Feeders with Distributed Generation using Positive Sequence Apparent Impedance
Fault Location in Distribution Feeders with Distributed Generation using Positive Sequence Apparent Impedance ARTURO SUMAN BRETAS Federal University of Rio Grande do Sul Department of Electrical Engineering
ChE 344 Winter 2013 Mid Term Exam II Tuesday, April 9, 2013
ChE 344 Winter 2013 Mid Term Exam II Tuesday, April 9, 2013 Open Course Textbook Only Closed everything else (i.e., Notes, In-Class Problems and Home Problems Name Honor Code (Please sign in the space
A First Course on Kinetics and Reaction Engineering Example 29.3
Example 29.3 Problem Purpose This problem will help you determine whether you have mastered the learning objectives for this unit. It illustrates the analysis of a parallel reactor network and the effect
Introduction to System Dynamics
SE1 Prof. Davide Manca Politecnico di Milano Dynamics and Control of Chemical Processes Solution to Lab #1 Introduction to System Dynamics Davide Manca Dynamics and Control of Chemical Processes Master
Automation in Complex Systems MIE090
Automation in Complex Systems MIE090 Exam Tuesday May 31, 2016 You may bring the course book and the reprints (defined in the course requirements), but not the solution to problems or your own solutions
Sensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transducers 2015 by IFSA Publishing, S. L. http://www.sensorsportal.com Multi-Model Adaptive Fuzzy Controller for a CSTR Process * Shubham Gogoria, Tanvir Parhar, Jaganatha Pandian B. Electronics
AC : BALANCING THEORY, SIMULATION AND PHYSICAL EXPERIMENTS IN
AC 2011-682: BALANCING THEORY, SIMULATION AND PHYSICAL EXPERIMENTS IN Anthony William Duva, Wentworth Institute of Technology Anthony W. Duva has been a faculty member in the Mechanical Engineering and
The simplified model now consists only of Eq. 5. Degrees of freedom for the simplified model: 2-1
. a) Overall mass balance: d( ρv ) Energy balance: = w + w w () d V T Tref C = wc ( T Tref ) + wc( T Tref ) w C T Because ρ = constant and ( Tref ) V = V = constant, Eq. becomes: () w = + () w w b) From
A First Course on Kinetics and Reaction Engineering Unit 22. Analysis of Steady State CSTRs
Unit 22. Analysis of Steady State CSRs Overview Reaction engineering involves constructing an accurate mathematical model of a real world reactor and then using that model to perform an engineering task
A First Course on Kinetics and Reaction Engineering Example 23.1
Example 23.1 parameter. Problem Purpose This problem illustrates the transient analysis of a CSTR following a change in an operating Problem Statement Recall the isothermal 4430 cm 3 steady-state chemostat
Lecture (9) Reactor Sizing. Figure (1). Information needed to predict what a reactor can do.
Lecture (9) Reactor Sizing 1.Introduction Chemical kinetics is the study of chemical reaction rates and reaction mechanisms. The study of chemical reaction engineering (CRE) combines the study of chemical
MATH 1231 MATHEMATICS 1B Calculus Section 4.3: - Series.
MATH 1231 MATHEMATICS 1B 2010. Calculus Section 4.3: - Series. 1. Sigma notation 2. What is a series? 3. The big question 4. What you should already know 5. Telescoping series 5. Convergence 6. n th term
Compressed Air Discharge from a Blowing Unit
The Open Mechanical Engineering Journal, 9, 3, 1-8 1 Compressed Air Discharge from a Blowing Unit Open Access G.L. Antinori and M. Spiga * Department of Industrial Engineering, University of Parma, Italy
CE 329, Fall 2015 Second Mid-Term Exam
CE 39, Fall 15 Second Mid-erm Exam You may only use pencils, pens and erasers while taking this exam. You may NO use a calculator. You may not leave the room for any reason if you do, you must first turn
Outline Calculus for the Life Sciences II. Pollution in a Lake 1. Introduction. Lecture Notes Numerical Methods for Differential Equations
Improved Improved Outline Calculus for the Life Sciences II tial Equations Joseph M. Mahaffy, mahaffy@math.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences
IJSRD - International Journal for Scientific Research & Development Vol. 1, Issue 8, 2013 ISSN (online):
IJSRD - International Journal for Scientific Research & Development Vol. 1, Issue 8, 2013 ISSN (online): 2321-0613 Development and theoretical analysis of mathematical expressions for change of entropy
Chemical, Biochemical, and Engineering Thermodynamics
Chemical, Biochemical, and Engineering Thermodynamics Fourth Edition Stanley I. Sandler University of Delaware John Wiley & Sons, Inc. Contents CHAPTER 1 INTRODUCTION 1 1.1 The Central Problems of Thermodynamics
MATHEMATICAL MODEL OF A VAPOUR ABSORPTION REFRIGERATION UNIT
ISSN 1726-4529 Int j simul model 9 (2010) 2, 86-97 Original scientific paper MATHEMATICAL MODEL OF A VAPOUR ABSORPTION REFRIGERATION UNIT Micallef, D. & Micallef, C. University of Malta, Department of
A First Course on Kinetics and Reaction Engineering Example 30.1
Example 30.1 Problem Purpose This problem will help you determine whether you have mastered the learning objectives for this unit. It illustrates the general approach to solving the design equations for
Chemical Reaction Engineering
Lecture 19 Chemical Reaction Engineering (CRE) is the field that studies the rates and mechanisms of chemical reactions and the design of the reactors in which they take place. oday s lecture Gas Phase
SYMBOLIC 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
Chemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017
Chemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017 Objective: Text: To introduce the basic concepts of fluid mechanics and heat transfer necessary for solution of engineering
Chemical Kinetics and Reaction Engineering
Chemical Kinetics and Reaction Engineering MIDTERM EXAMINATION II Friday, April 9, 2010 The exam is 100 points total and 20% of the course grade. Please read through the questions carefully before giving
Chemical and Engineering Thermodynamics
Chemical and Engineering Thermodynamics Third Edition Stanley I. Sandler University of Delaware John Wiley & Sons, Inc. New York Chichester Weinheim Brisbane Singapore Toronto Contents NOTATION xv CHAPTER1
Modeling of a Fluid Catalytic Cracking (FCC) Riser Reactor
Modeling of a Fluid Catalytic Cracking (FCC) Riser Reactor Dr. Riyad Ageli Mahfud Department of Chemical Engineering- Sabrattah Zawia University Abstract: Fluid catalytic cracking (FCC) is used in petroleum
De-Coupler Design for an Interacting Tanks System
IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331, Volume 7, Issue 3 (Sep. - Oct. 2013), PP 77-81 De-Coupler Design for an Interacting Tanks System
Process Classification
Process Classification Before writing a material balance (MB) you must first identify the type of process in question. Batch no material (mass) is transferred into or out of the system over the time period