Short-cut distillation columns based on orthogonal polynomials of a discrete variable
|
|
- Meredith Lamb
- 6 years ago
- Views:
Transcription
1 Short-cut distillation columns based on orthogonal polynomials of a discrete variable Peter Lory Institut für Wirtschaftsinformatik, Universität Regensburg, D Regensburg, Germany Dedicated to Professor Roland Bulirsch on the occasion of his 65th birthday Summary. Increasing competition in the process industries enforces the application of mathematical simulation techniques both in the design phase and in the operating phase of a plant. The basic apparatus for separation processes is the distillation column. Its rigorous (tray by tray) mathematical modelling results in a system of simultaneous nonlinear equations (algebraic in the steady-state case, differential-algebraic in the dynamic case). For high (however realistic) numbers of trays and components, these systems may become rather large (thousands of equations). In addition, realistic plant models often include several distillation columns. As a consequence, the numerical solution of these models may become difficult and time consuming. This has led to attempts to model the distillation columns less rigorously with the aim to achieve a considerable reduction in the number of equations. The name short-cut distillation columns is common for models of this type. The present paper uses a discrete weighted residual method for the development of short-cut models. It suggests a Galerkin method based on orthogonal polynomials of a discrete variable, the tray number. It is a remarkable advantage of this technique that even very coarse models satisfy all global balances exactly. Keywords: Mathematical models, distillation columns, discrete weighted residual method, discrete Galerkin method, Hahn polynomials. 1
2 1. Introduction Separation processes are ubiquitous in chemical engineering, see e.g. Holland (1981) or King (1980). The basic apparatus for these purposes is the distillation column. A tray distillation column consists of numerous stacked equilibrium stages that operate at different pressures and temperatures. The following system of equations is a steady-state model of a distillation column with N trays. The index i refers to the number of a tray (stage) starting at the bottom (i = 1) and ending at the top (i = N). The mixture consists if n c components. The index j denotes its number. L i+1 + V i 1 L i V i + F i = 0, (1.1a) L i+1 x i+1,j + V i 1 y i 1,j L i x i,j V i y i,j + F i z i,j = 0 (j = 1,..., n c ), (1.1b) n c y i,j = 1, (1.1c) j=1 y i,j K i,j (T i, P i, x i, y i ) x i,j = 0 (j = 1,..., n c ), (1.1d) L i = Φ(M i ), (1.1e) L i+1 h L i+1 + V i 1 h V i 1 L i h L i V i h V i + F i h F i + Q i = 0, (1.1f) h L i = h L i (T i, P i, x i ), h V i = h V i (T i, P i, y i ), (1.1g) P i P i+1 = Ψ(M i+1, V i ). (1.1h) These equations determine the following state variables: M i : molar liquid holdup on stage i, L i : molar liquid flowrate descending from stage i to stage i 1, V i : molar vapour flowrate ascending from stage i to stage i + 1, x i,j : mole fraction of component j in the liquid at stage i, y i,j : mole fraction of component j in the vapour at stage i, h L i, h V i : molar enthalpies of liquid and vapour at stage i, T i : temperature at stage i, P i : pressure at stage i. The following quantities are given: the external heat input Q i to stage i, the feed flowrate F i to stage i as well as its mole fractions z i,j and its molar enthalpy h F i, the liquid flowrate L N+1 to the top tray together with its mole fractions x N+1,j and its molar enthalpy h L N+1, the vapour flowrate V 0 to the bottom tray as well as its composition y 0,j and the molar 2
3 enthalpy h V 0, the external pressure at the column top P N+1. The row vectors x i collect the mole fractions of the liquid at stage i: x i = (x i,1,..., x i,nc ). The vectors y i are defined analogously. Equation (1.1a) is the total material balance on stage i, while (1.1b) are individual material balances for the n c components. Equation (1.1c) expresses the fact, that the mole fractions in vapour sum up to 1 (see below). Equations (1.1d) describe the thermodynamic equilibrium relation between the mole fractions in the vapour phase and in the liquid phase. The so-called K values K i,j are determined by a physical property system based on given T i, P i, x i,j, y i,j. The flowrate L i of the liquid leaving stage i is a known function Φ of the holdup M i on this stage (1.1e). Equation (1.1f) describes the total enthalpy balance over stage i, with the liquid and vapour molar enthalpies defined by (1.1g). Both values are again given by a physical property system. Equation (1.1h) expresses the pressure drop between two consecutive stages. For more details, see Eich-Soellner et al. (1997). It is assumed that both the mole fractions x N+1,j of the liquid flow to the top tray and the mole fractions y 0,j of the vapour flow to the bottom tray and the mole fractions z i,j of each feed sum up to 1. Then it is a consequence of (1.1a) and (1.1b) that n c x i,j = 1 (i = 1,..., N). (1.2) j=1 As a consequence, these equations can substitute the group of equations with j = n c in (1.1b). However, the present paper stresses the exact satisfaction of the global balances (see Section 3). Thus, it is reasonable to enforce all the local mass balances in (1.1b) explicitly. The system (1.1) consists of N(2n c + 7) equations for the same number of unknown variables. In practical applications, the number of trays is within 10 N 100, and the number of components is less than 50. As a consequence, it is not uncommon that the the system (1.1) comprises thousands of equations. Additionally, a realistic plant may include several distillation columns. Consequently, the numerical solution of these models may become difficult and time consuming. This has led to attempts to reduce the number of equations in distillation column models without losing the required level of detail. The name short cut distillation column is common for a model of this type. Most of these approaches use engineering experience and have a low level of mathematical sophistication. The paper of Eckert & Vanek (1985) is a representative for short-cut models of this type. Other authors, in particular Stewart et al. (1985), use discrete weighted residual methods for the development of short-cut models. The article of Canu & Ray (1991) gives a survey of this and subsequent papers. These authors use collocation and cannot guarantee the exact satisfaction of the global balances (see Section 3). The present paper uses a discrete weighted residual method, too. However, the weighting functions are chosen according to the discrete Galerkin method. As shown in 3
4 Section 3, this has the advantage that the global balances are satisfied exactly (up to roundoff errors and up to the errors due to the nonlinear equation solver, of course). In other words, the reduction in the dimensionality of the model does not diminish the exactness of the global balances. The following section (Section 2) describes the discrete Galerkin method, in particular the use of the Hahn polynomials as basis and test functions within the frame of the discrete weighted residual method. Section 3 shows that the proposed method produces results that satisfy the global balances exactly. Section 4 summarises first numerical experiments. 2. The discrete Galerkin method At those trays, where a feed of material or energy enters the column, the solutions obviously change abruptly. As a consequence, it is convenient for the short-cut approach to partition the column in sections, where no feed is present (e.g. rectifying or stripping sections). In the following, it is assumed that such a section is given. Its modelling equations are given by (1.1) with the exception that feeds can occur only at the bottom or top tray. Let the trays in the rigorous model for such a section be numbered from i = 1 (bottom) to i = N (top). Then the feeds of this section can be added to the corresponding flowrates (L N+1, L N+1 x N+1,j, V 0, V 0 y 0,j ) and enthalpies (L N+1 h L N+1, V 0 h V 0 ) that enter the section. As a result, the transformed (equivalent) rigorous model of this column section does not contain any feeds F i, F i z i,j, F i h F i or Q i explicitly. This is assumed in the following. It has been suggested to built a short-cut model by transforming the discrete variable i into a continuous one. The paper of Canu & Ray (1991) gives a review of these attempts. This approach yields a differential equation in the continuous variable, that has to be discretized for its numerical solution. However, it seems much more natural to retain the discrete nature of the original problem. This is pointed out in the papers of Deuflhard & Wulkow (1989, 1995) in the context of a different field of application, namely the computational treatment of polyreaction kinetics. The present paper also uses a discrete version of the classical Galerkin method. For that purpose, the unknown functions are approximated by linear combinations of certain basis functions. E.g. for the flowrates: n n L(s) = Λ k H k (s), V (s) = Φ k H k (s), s {1, 2,..., N}. (2.1) k=0 k=0 The discrete Galerkin method is a special case of the discrete weighted residual method and needs a scalar product for its formulation: The discrete Galerkin equations for (1.1a) are N f, g := f(s)g(s). (2.2) s=1 L(s + 1) + V (s 1) L(s) V (s), H l (s) = 0 (l = 0, 1,..., n). (2.3) 4
5 The Galerkin equations for (1.1b) (1.1h) are analogous. Altogether, a nonlinear algebraic system of (n + 1)(2n c + 7) equations for the same number of coefficients in the representations of the state variables arises. The number n has to be selected as small as possible without sacrificing too much accuracy. It is computationally convenient to use orthogonal polynomials as the basis functions H k (s). The scalar product (2.2) suggests the use of Hahn polynomials, which are defined by Q k (s; α, β, N) = k ν=0 ( k) ν ( s) ν (k + α + β + 1) ν ( N) ν (α + 1) ν ν! see e.g. Askey (1975). Here, the definitions (a) 0 = 1, (a) ν = a (a + 1)... (a + ν 1) for ν > 0, α, β > 1, k = 0,..., N ; (2.4) have been used. For α = β = 0, these polynomials satisfy the following orthogonality relation: 0 for k l N Q k (s; 0, 0, N) Q l (s; 0, 0, N) = (2.5) N + 1 s=0 for k = l, π k where π k = ( 1)k ( N) k (2k + 1). (N + 2) k As a consequence, the polynomials H k (s) = Q k (s 1; 0, 0, N 1) (2.6) are orthogonal with respect to the scalar product (2.2). The present paper uses these polynomials both as basis and as test (weighting) functions. The paper of Karlin & McGregor (1961) lists some interesting relations for Hahn polynomials. Please note that they use a slightly different notation: their polynomials Q k (s; α, β, N) are orthogonal on s = 0, 1,..., N Global balances According to Section 2, it is assumed that the rigorous model of a column section has been transformed to an equivalent model with no explicit feeds. Straightforward calculations show that such a rigorous model (1.1) satisfies the following global balances for material and enthalpy: L N+1 + V 0 L 1 V N = 0, (3.1a) L N+1 x N+1,j + V 0 y 0,j L 1 x 1,j V N y N,j = 0 (j = 1,..., n c ), (3.1b) L N+1 h L N+1 + V 0 h V 0 L 1 h L 1 V N h V N = 0, (3.1c) 5
6 It is a very attractive feature of the Galerkin method as proposed in Section 2 that even for arbitrary n the results of the short-cut model satisfy the same balances. Indeed, the Galerkin equations (2.3) contain the equation L(s + 1) + V (s 1) L(s) V (s), H 0 (s) = 0 (3.2) with the Hahn polynomial H 0 of degree 0 as test function. This polynomial is a constant, and from Equation (3.2) follows immediately L(N + 1) + V (0) L(1) V (N) = 0. (3.3) This is the mass balance (3.1a). The other global balances follow analogously. Remark: Stewart et al. (1985) have developed short-cut models for distillation columns on the basis of discrete weighted residual methods, too. These authors use collocation with the zeros of the Hahn polynomials as collocation points. See Canu & Ray (1991) for a survey of this and subsequent papers. However, these methods cannot guarantee the exact satisfaction of the global balances. Intuitively, the Galerkin approach seems to be more adequate than collocation because the equations of the former method merge the different stages to a higher degree. 4. Numerical experiments The test problems of Fletcher (1990) have been used as the basis for first numerical experiments. These problems are simplified versions of the system (1.1). They do not include holdups and pressures as variables. Consequently, Equations (1.1e) and (1.1h) are dropped. The K-values are given by the relationship ( ) K i,j = 1 P i exp b j a j + c j + T i, (4.1) where the Antoine constants a j, b j, c j (j = 1,..., n c ) and the stage pressures P i are data values. The molar enthalpies are given by n c h L i = j=1 x i,j (α j + α jt i + α j T 2 i ), h V i = n c j=1 y i,j (β j + β jt i + β j T 2 i ), (4.2) where the liquid enthalpy constants α j, α j, α j and the vapour enthalpy constants β j, β j, β j (j = 1,..., n c ) are again data values. In the following, these and other necessary data are given in the schematic form of Table 1. Here, f 1,..., f nc are the individual flowrates of the feed components. In the case of a rectifying section, this means a vapour feed to the bottom tray; in the case of a stripping section it is a liquid feed to the top tray. In any case, T f is the temperature of this feed. All computations have been performed with the function fsolve of the Optimization Toolbox of MATLAB with a precision of roughly 16 significant decimal digits. The termination criteria both on the independent variables and on the function (options(2) and options(3)) were 1.E 8. 6
7 N, n c ; a 1, b 1, c 1,..., a nc, b nc, c nc ; α 1, α 1, α 1,..., α nc, α n c, α n c ; β 1, β 1, β 1,..., β nc, β n c, β n c ; f 1, f 2..., f nc, T f ; P 1, P 2,..., P N ; L N+1, T N+1, x N+1,1, x N+1,2,..., x N+1,nc ; V 0, T 0, y 0,1, y 0,2,..., y 0,nc. Table 1: Scheme of numerical data Example 1 Example 1 is the stripping section (without the reboiler) of the hydrocarbon-20 problem of Fletcher (1990). Its data are given in Table 2. Figures 1 and 2 give the results for the mole fractions and for the temperatures at the different stages E Table 2: Data of Example 1 Table 3 summarises the results of the short-cut models. For the most important components of the solution vector, it gives the maxima of the relative errors between the short-cut solutions for n = 2, 3, 4 and the solution of the rigorous model. In Example 1, the components x i,1 and y i,1 are very small in comparison with the others, and the absolute errors are given in these cases. The number n is the maximum degree of the Hahn polynomials that are involved in the corresponding short-cut model. 7
8 Mole fractions Number of tray Figure 1: Example 1 ( refers to liquid, refers to vapour) Temperatures Number of tray Figure 2: Example 1 8
9 Example 1 Example 2 n = 2 n = 3 n = 4 n = 2 n = 3 n = 4 x i, % % % x i, % % % % % % x i, % % % y i, % % % y i, % % % % % % y i, % % % L i % % % % % % V i % % % % % % T i % % % % % % Table 3: Errors between short-cut solutions and the rigorous solution Example 2 Example 2 is an enlarged version of the rectifying section (without the condenser) of the methanol-8 problem of Fletcher (1990). Its data are given in Table 4. Figures 3 and 4 show the results for the mole fractions and for the temperatures at the different stages. Table 3 presents the results of the short-cut models. Here, for all components the relative errors are given Table 4: Data of Example 2 Final remarks The numerical experiments show that the proposed approach for short-cut distillation columns is very promising. In accordance with the theoretical results of Section 3, 9
10 Mole fractions Number of tray Figure 3: Example 2 ( refers to liquid, refers to vapour) Temperatures Number of tray Figure 4: Example 2 the balances (3.1) were statisfied with a relative precision of at least 2.2 E 12 in all rigorous and short-cut computations. This feature is particularly attractive when shortcut models are used as units in large flowsheets, see e.g. Eich-Soellner et al. (1997). 10
11 References Askey, R., Orthogonal Polynomials and Special Functions, Philadelphia, Society for Industrial and Applied Mathematics. Canu P. & Ray W. H., Discrete weighted residual methods applied to polymerization reactions. Computers Chem. Engng. 15, Deuflhard P. & Wulkow M., Computational treatment of polyreaction kinetics by orthogonal polynomials of a discrete variable. IMPACT Comput. Sci. Eng. 1, Deuflhard P. & Wulkow M., Simulationsverfahren für die Polymerchemie. In: Mathematik in der Praxis (A. Bachem et al., Eds.). Berlin, Springer. Pp Eckert E. & Vanek T., Simulation of stage separation columns based on the block short-cut approach. Chemical Engineering Science 40, Eich-Soellner E., Lory P., Burr P. & Kröner A., Stationary and dynamic flowsheeting in the chemical engineering industry. Surveys on Mathematics for Industry 7, Fletcher R., Distillation column test problem. In: Computational Solution of Nonlinear Systems of Equations (E. L. Allgower et al., Eds.). Lectures in Applied Mathematics 26, Providence, American Mathematical Society. Pp Holland, C. D., McGraw-Hill. Fundamentals of Multicomponent Distillation, New York, Karlin S. & McGregor J. L., The Hahn polynomials, formulas and an application. Scripta Mathematica 26, King, C. J., Separation Processes, New York, McGraw-Hill. Stewart W. E., Levien K. L. & Morari M., Simulation of fractionation by orthogonal collocation. Chemical Engineering Science 40,
BOUNDARY 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 informationAll Rights Reserved. Armando B. Corripio, PhD, P.E., Multicomponent Distillation Column Specifications... 2
Multicomponent Distillation All Rights Reserved. Armando B. Corripio, PhD, P.E., 2013 Contents Multicomponent Distillation... 1 1 Column Specifications... 2 1.1 Key Components and Sequencing Columns...
More informationMass Transfer Operations I Prof. Bishnupada Mandal Department of Chemical Engineering Indian Institute of Technology, Guwahati
Mass Transfer Operations I Prof. Bishnupada Mandal Department of Chemical Engineering Indian Institute of Technology, Guwahati Module - 5 Distillation Lecture - 6 Fractional Distillation: McCabe Thiele
More information[Thirumalesh*, 4.(8): August, 2015] ISSN: (I2OR), Publication Impact Factor: 3.785
IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY CASE STUDY ON MULTICOMPONENT DISTILLATION AND DISTILLATION COLUMN SEQUENCING Thirumalesh. B.S*, Ramesh.V * Assistant Professor
More informationChapter 10. Vapor/Liquid Equilibrium: Introduction
Chapter 10 Vapor/Liquid Equilibrium: Introduction Preceding chapters have dealt largely with pure substances or with constant-composition mixtures. e.g., air. However, composition changes are the desired
More informationLatin American Applied Research an International Journal of Chemical Engineering. 1997, vol. 27, 1-2, (Accepted: March 20, 1996).
Latin American Applied Research an International Journal of Chemical Engineering. 1997, vol. 27, 1-2, 51-60 (Accepted: March 20, 1996). NEW METHOD FOR DESIGNING DISTILLATION COLUMNS OF MULTICOMPONENT MIXTURES
More informationPRESSURE SWING BATCH DISTILLATION FOR HOMOGENOUS AZEOTROPIC SEPARATION
PRESSURE SWING BATCH DISTILLATION FOR HOMOGENOUS AZEOTROPIC SEPARATION Jens-Uwe Repke 1, Andreas Klein 1, David Bogle 2,Günter Wozny 1 1 Technical University Berlin, Department of Process- and Plant Dynamics,
More informationDynamic Effects of Diabatization in Distillation Columns
Downloaded from orbit.dtu.dk on: Dec 27, 2018 Dynamic Effects of Diabatization in Distillation Columns Bisgaard, Thomas Publication date: 2012 Document Version Publisher's PDF, also known as Version of
More informationComparison of Conventional and Middle Vessel Batch Reactive Distillation Column: Application to Hydrolysis of Methyl Lactate to Lactic Acid
A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 35, 2013 Guest Editors: Petar Varbanov, Jiří Klemeš, Panos Seferlis, Athanasios I. Papadopoulos, Spyros Voutetakis Copyright 2013, AIDIC Servizi
More informationMODULE 5: DISTILLATION
MOULE 5: ISTILLATION LECTURE NO. 3 5.2.2. Continuous distillation columns In contrast, continuous columns process a continuous feed stream. No interruptions occur unless there is a problem with the column
More informationRigorous column simulation - SCDS
Rigorous column simulation - SCDS Problem & solution principle: A rectification can be realistically displayed by performing a rigorous column simulation. It is possible to simulate ideal and real mixtures.
More informationDynamics and Control of Energy Integrated Distillation Column Networks
200 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 30-July 02, 200 ThA4.4 Dynamics and Control of Energy Integrated Distillation Column Networks Sujit S. Jogwar and Prodromos
More informationEXTENDED SMOKER S EQUATION FOR CALCULATING NUMBER OF STAGES IN DISTILLATION
EXTENDED SMOKER S EQUATION FOR CALCULATING NUMBER OF STAGES IN DISTILLATION Santanu Bandyopadhyay Energy Systems Engineering and Department of Mechanical Engineering, Indian Institute of Technology, Bombay,
More informationAn Efficient Design of Multi Component Distillation Column by Approximate & Rigorous Method
An Efficient Design of Multi Component Distillation Column by Approximate & Rigorous Method Syed Mujahed Ali Rizwan Senior Lecturer in Chemistry Challenger College, Moinabad, Hyderabad. Abstract: In this
More informationWhat is "In" and What is "Out" in Engineering Problem Solving
18 th European Symposium on Computer Aided Process Engineering ESCAPE 18 Bertrand Braunschweig and avier Joulia (Editors) 008 Elsevier B.V./Ltd. All rights reserved. M. Shacham, M. B. Cutlip and N. Brauner,"What
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 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 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 informationComputer Aided Design Module for a Binary Distillation Column
Computer Aided Design Module for a Binary Distillation Column K. R. Onifade Department of Chemical Engineering, Federal University of Technology Minna, Nigeria Abstract A Computer Aided Design (CAD) module
More informationINTEGRATION OF DESIGN AND CONTROL FOR ENERGY INTEGRATED DISTILLATION
INTEGRATION OF DESIGN AND CONTROL FOR ENERGY INTEGRATED DISTILLATION Hongwen Li, Rafiqul Gani, Sten Bay Jørgensen CAPEC, Department of Chemical Engineering Technical University of Denmark, Lyngby, Denmark
More informationOptimization of Batch Distillation Involving Hydrolysis System
273 Optimization of Batch Distillation Involving Hydrolysis System Elmahboub A. Edreder 1, Iqbal M. Mujtaba 2, Mansour Emtir 3* 1 Libyan Petroleum Institute, P.O. Box 6431, Tripoli, Libya 2 School of Engineering
More informationMULTI-LOOP CONTROL STRUCTURE FOR DIVIDING-WALL DISTILLATION COLUMNS. Abstract. Introduction
MULTI-LOOP CONTROL STRUCTURE FOR DIVIDING-WALL DISTILLATION COLUMNS Salvador Tututi-Avila, Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute, Troy, NY, Arturo Jiménez-Gutiérrez,
More informationA comparative study on the recovery of 1,2-dichloroethane and the removal of benzene contained in the byproducts of VCM process
Korean J. Chem. Eng., 27(3), 970-976 (2010) DOI: 10.1007/s11814-010-0206-7 RAPID COMMUNICATION A comparative study on the recovery of 1,2-dichloroethane and the removal of benzene contained in the byproducts
More informationMass Transfer Operations I Prof. Bishnupada Mandal Department of Chemical Engineering Indian Institute of Technology, Guwahati
Mass Transfer Operations I Prof. Bishnupada Mandal Department of Chemical Engineering Indian Institute of Technology, Guwahati Module - 5 Distillation Lecture - 5 Fractional Distillation Welcome to the
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 3, March ISSN
International Journal of Scientific & Engineering Research Volume 6 Issue 3 March-25 69 ISSN 2229-558 Modeling and Simulation of a Distillation olumn using ASPEN PLUS for separating methanol/water mixture
More informationShortcut Design Method for Columns Separating Azeotropic Mixtures
3908 Ind. Eng. Chem. Res. 2004, 43, 3908-3923 Shortcut Design Method for Columns Separating Azeotropic Mixtures Guilian Liu, Megan Jobson,*, Robin Smith, and Oliver M. Wahnschafft Department of Process
More informationReduced Order Modeling of High Purity Distillation Columns for Nonlinear Model Predictive Control
Reduced Order Modeling of High Purity Distillation Columns for Nonlinear Model Predictive Control Suabtragool Khowinij, Shoujun Bian and Michael A Henson Department of Chemical Engineering University of
More informationRATE-BASED MODELING OF TWO COMMERCIAL SCALE H 2 S STRIPPING COLUMNS
Distillation Absorption A.B. de Haan, H. Kooijman and A. Górak (Editors) All rights reserved by authors as per DA copyright notice RATE-BASED MODELING OF TWO COMMERCIAL SCALE H 2 S STRIPPING COLUMNS Ross
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 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 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 informationReflections on the use of the McCabe and Thiele method
From the Selectedorks of João F Gomes January 2007 Reflections on the use of the McCabe and Thiele method Contact Author Start Your Own Selectedorks Notify Me of New ork Available at: http://works.bepress.com/joao_gomes/42
More informationMODELING AND SIMULATION OF DISTILLATION COLUMN
MODELING AND SIMULATION OF DISTILLATION COLUMN Ronak Prajapati, Parth Patel 2, Vishal Parmar 3, Nirali Tharwala 4,2,3 student, Chemical Engineering, Department of Chemical Engineering, Shroff S R Rotary
More informationChemSep Case Book: Handling Missing Components
ChemSep Case Book: Handling Missing Components Ross Taylor and Harry Kooijman In their book Conceptual Design of Distillation Systems M.F. Doherty and M.F. Malone (McGraw-Hill, 2001) describe a process
More informationAvailable online Journal of Scientific and Engineering Research, 2018, 5(2): Research Article
Available online www.jsaer.com, 2018, 5(2):32-39 Research Article ISSN: 2394-2630 CODEN(USA): JSERBR A Shortcut Method for Modeling Multi-Component Distillation in an Existing Column Philomena K. Igbokwe
More informationSteady State Multiplicity and Stability in a Reactive Flash
Steady State Multiplicity and Stability in a Reactive Flash Iván E. Rodríguez, Alex Zheng and Michael F. Malone Department of Chemical Engineering University of Massachusetts Amherst, MA 01003 Abstract
More informationTHERMAL INTEGRATION OF A DISTILLATION COLUMN THROUGH SIDE-EXCHANGERS
THERMAL INTEGRATION OF A DISTILLATION COLUMN THROUGH SIDE-EXCHANGERS Santanu Bandyopadhyay Energy Systems Engineering and Department of Mechanical Engineering, Indian Institute of Technology, Bombay, Powai,
More informationTHERMODYNAMIC ANALYSIS OF MULTICOMPONENT DISTILLATION-REACTION PROCESSES FOR CONCEPTUAL PROCESS DESIGN
THERMODYNAMIC ANALYSIS OF MULTICOMPONENT DISTILLATION-REACTION PROCESSES FOR CONCEPTUAL PROCESS DESIGN Oliver Ryll, Sergej Blagov 1, Hans Hasse Institute of Thermodynamics and Thermal Process Engineering,
More informationCHEMICAL 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,
More informationMinimum Energy Demand and Split Feasibility for a Class of Reactive Distillation Columns
Minimum Energy emand and Split easibility for a Class of Reactive istillation Columns Rosa Y. Urdaneta, Jürgen ausa + and olfgang Marquardt * Lehrstuhl für Prozesstechnik RTH Aachen University, Templergraben
More informationDesign and Analysis of Divided Wall Column
Proceedings of the 6th International Conference on Process Systems Engineering (PSE ASIA) 25-27 June 2013, Kuala Lumpur. Design and Analysis of Divided Wall Column M. Shamsuzzoha, a* Hiroya Seki b, Moonyong
More informationPerformance of esterification system in reaction-distillation column
Performance of esterification system in reaction-distillation column Proceedings of European Congress of Chemical Engineering (ECCE-6) Copenhagen, 16-20 September 2007 Performance of esterification system
More informationA First Course on Kinetics and Reaction Engineering Unit 4. Reaction Rates and Temperature Effects
Unit 4. Reaction Rates and Temperature Effects Overview This course is divided into four parts, I through IV. Part II is focused upon modeling the rates of chemical reactions. Unit 4 is the first unit
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 informationMake distillation boundaries work for you!
Make distillation boundaries work for you! Michaela Tapp*, Simon T. Holland, Diane Hildebrandt and David Glasser Centre for Optimization, Modeling and Process Synthesis School of Process and Materials
More informationPlacement and Integration of Distillation column Module 06 Lecture 39
Module 06: Integration and placement of equipment Lecture 39: Placement and Integration of Distillation Column Key word: Pinch Technology, Pinch technology has established that good process integration
More informationMultivariable model predictive control design of reactive distillation column for Dimethyl Ether production
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Multivariable model predictive control design of reactive distillation column for Dimethyl Ether production To cite this article:
More informationIncorporating Reality Into Process Simulation. Nathan Massey Chemstations, Inc. January 10, 2002
Incorporating Reality Into Process Simulation Nathan Massey Chemstations, Inc. January 10, 2002 Levels of Reality in Process Simulation 1. Data Accuracy and Comprehensiveness Physical Properties Phase
More informationPRINCIPLES AND MODERN APPLICATIONS OF MASS TRANSFER OPERATIONS
PRINCIPLES AND MODERN APPLICATIONS OF MASS TRANSFER OPERATIONS Jaime Benitez iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC., PUBLICATION Preface Nomenclature xiii xv 1. FUNDAMENTALS OF MASS TRANSFER 1
More informationA NON-LINEAR WAVE MODEL WITH VARIABLE MOLAR FLOWS FOR DYNAMIC BEHAVIOUR AND DISTURBANCE PROPAGATION IN DISTILLATION COLUMNS
A NON-LINEAR WAVE MODEL WITH VARIABLE MOLAR FLOWS FOR DYNAMIC BEHAVIOUR AND DISTURBANCE PROPAGATION IN DISTILLATION COLUMNS Nick Hankins Department of Engineering Science, The University of Oxford, Parks
More informationComparison of distillation arrangement for the recovery process of dimethyl sulfoxide
Korean J. Chem. Eng., 24(3), 438-444 (2007) SHORT COMMUNICATION Comparison of distillation arrangement for the recovery process of dimethyl sulfoxide Jungho Cho and Dong Min Kim* Department of Chemical
More informationReview of Chemical Equilibrium Introduction
Review of Chemical Equilibrium Introduction Copyright c 2015 by Nob Hill Publishing, LLC This chapter is a review of the equilibrium state of a system that can undergo chemical reaction Operating reactors
More informationRate-based design of integrated distillation sequences
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 Rate-based design of integrated distillation sequences
More informationGRAPHICAL REPRESENTATION OF THE CHANGES OF SECTOR FOR PARTICULAR CASES IN THE PONCHON SAVARIT METHOD
GRAPICAL REPRESENTATION OF TE CANGES OF SECTOR FOR PARTICULAR CASES IN TE PONCON SAVARIT METOD J.A. REYES-LABARTA *, M.D. SERRANO and A. MARCILLA *Chemical Engng Department, University of Alicante, Apdo.
More informationA Rate-Based Equation-Oriented Parallel Column Model: Application to Dividing Wall Columns
A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 69, 2018 Guest Editors: Elisabetta Brunazzi, Eva Sorensen Copyright 2018, AIDIC Servizi S.r.l. ISBN 978-88-95608-66-2; ISSN 2283-9216 The Italian
More informationBatch extractive distillation of mixture methanol-acetonitrile using aniline as a asolvent
48 Pol. J. Chem. Polish Tech., Journal Vol. 14, of No. Chemical 3, 2012 Technology, 14, 3, 48 53, 10.2478/v10026-012-0083-4 Batch extractive distillation of mixture methanol-acetonitrile using aniline
More information,, Seong-Bo Kim,Hai-SongBae, and Jeong-Sik Han
Jungho Cho, So-Jin Park,, Myung-Jae Choi,, Seong-Bo Kim,Hai-SongBae, and Jeong-Sik Han Department of Chemical Engineering, Dong-Yang University, Kyoungbuk, 750-711, Korea *Department of Chemical Engineering,
More informationA 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
More informationRigorous calculation of LNG flow reliefs using the GERG-2004 equation of state
Rigorous calculation of LNG reliefs using the GERG-2004 equation of state Luigi Raimondi Process Simulation Services www.xpsimworld.com Via Galvani 105, 20025 Legnano (MI) - Italy The design of process
More informationThe most important nomenclature used in this report can be summarized in:
Notation The most important nomenclature used in this report can be summarized in: V Vapor flow rate V T Vapor flow rate in the top L Liquid flow rate D Distillation product B Bottom product q Liquid fraction
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 informationNonlinear Dynamic Modelling of Methanol-Water Continuous Distillation Column النمذجة الذيناميكية الالخطية لعمود التقطير الميثانول الماء المستمر
Nonlinear Dynamic Modelling of Methanol-Water Continuous Distillation Column النمذجة الذيناميكية الالخطية لعمود التقطير الميثانول الماء المستمر Lecturer / Asseel M. Rasheed Electrical & Electronics Eng.
More informationNeural Networks for Two-Group Classification Problems with Monotonicity Hints
Neural Networks for Two-Group Classification Problems with Monotonicity Hints P. Lory 1, D. Gietl Institut für Wirtschaftsinformatik, Universität Regensburg, D-93040 Regensburg, Germany Abstract: Neural
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 informationReview of Chemical Equilibrium Introduction
Review of Chemical Equilibrium Introduction Copyright c 2015 by Nob Hill Publishing, LLC This chapter is a review of the equilibrium state of a system that can undergo chemical reaction Operating reactors
More informationCHAPTER SIX THERMODYNAMICS Vapor-Liquid Equilibrium in a Binary System 6.2. Investigation of the Thermodynamic Properties of Pure Water
CHAPTER SIX THERMODYNAMICS 6.1. Vapor-Liquid Equilibrium in a Binary System 6.2. Investigation of the Thermodynamic Properties of Pure Water 2 6.1. VAPOR-LIQUID EQUILIBRIUM IN A BINARY SYSTEM Keywords:
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 informationProblem Solving Using CAPE-OPEN Software: Residue Curves Case Study
Contemporary Engineering Sciences, Vol. 11, 2018, no. 78, 3857-3864 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ces.2018.88409 Problem Solving Using CAPE-OPEN Software: Residue Curves Case Study
More informationData Reconciliation of Streams with Low Concentrations of Sulphur Compounds in Distillation Operation
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 Data Reconciliation of Streams with Low Concentrations
More informationSpecific Problems of Using Unisim Design in the Dynamic Simulation of the Propylene-Propane Distillation Column
Specific Problems of Using Unisim Design in the Dynamic Simulation of the Propylene-Propane Distillation Column CRISTIAN PATRASCIOIU*, MARIAN POPESCU, NICOLAE PARASCHIV Petroleum - Gas University of Ploieºti,
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 informationModeling, Simulation and Optimization of a Cross Flow Molten Carbonate Fuel Cell
Modeling, Simulation and Optimization of a Cross Flow Molten Carbonate Fuel Cell P. Heidebrecht 1, M. Mangold 2, M. Gundermann 1, A. Kienle 2,3, K. Sundmacher 1,2 1 Otto-von-Guericke-University Magdeburg,
More informationCOPYRIGHTED MATERIAL INTRODUCTION CHAPTER 1
CHAPTER 1 INTRODUCTION The development of the chemical industry over the last two centuries has provided modern civilization with a whole host of products that improve the well-being of the human race.
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 informationDistillation. Sep-tek. Ch.11 (continued) Distillation (Multistage with reflux) Sigurd Skogestad. Separation of liquid mixtures by repeated evaporation
Sep-tek. Ch.11 (continued) Distillation (Multistage with reflux) Sigurd Skogestad Distillation Separation of liquid mixtures by repeated evaporation multi-stage with reflux Old name: Rectification Basis:
More informationComprehend and execute the 10 elements of effective problem
Lecture 8, 3/9/2012 Chapter 7: A GENERAL Strategy for Solving Material Balance Problems Objectives: Comprehend and execute the 10 elements of effective problem Drive a flow chart and Place labels on the
More informationSolving Equations and Optimizing Functions
Solving Equations and Optimizing Functions At the end of this lecture, you will be able to: solve for the roots of a polynomial using polyroots. obtain approximate solutions to single equations from tracing
More informationChemical 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
More informationApproximate Methods Fenske-Underwood-Gilliland (FUG) Method Selection of Two Key Components
Lecture 3. Approximate Multicomponent Methods () [Ch. 9] Approximate Methods Fenske-Underwood-Gilliland (FUG) Method Selection of Two Key Components Column Operating Pressure Fenske Equation for Minimum
More informationConceptual Design of Reactive Distillation Columns with Non-Reactive Sections
Conceptual esign of Reactive istillation Columns with Non-Reactive Sections R. M. ragomir, M. Jobson epartment of Process Integration, UMIST, PO ox 88, M60 Q, Manchester, UK Abstract Reactive distillation
More informationSelecting the Appropriate Numerical Software for a Chemical Engineering Course
Selecting the Appropriate Numerical Software for a Chemical Engineering Course Mordechai Shacham, shacham@bgumail.bgu.ac.il Chemical Engineering Department, Ben-Gurion University of the Negev, Beer-Sheva
More informationEXPERIMENTAL SIMULATION OF DISTILLATION COLUMN PROFILE MAPS. Tshepo Sehole David Modise
EXPERIMENTAL SIMULATION OF DISTILLATION COLUMN PROFILE MAPS Tshepo Sehole David Modise A thesis submitted to the Faculty of Engineering and the Built Environment, University of the Witwatersrand, Johannesburg,
More informationTHEORETICAL AND EXPERIMENTAL STUDIES ON STARTUP STRATEGIES FOR A HEAT-INTEGRATED DISTILLATION COLUMN SYSTEM
HEOREICAL AND EXPERIMENAL SUDIES ON SARUP SRAEGIES FOR A HEA-INEGRAED DISILLAION COLUMN SYSEM M. Wendt, R. Königseder, P. Li, G. Wozny Institut für Prozess und Anlagen echnik, echnische Universität erlin
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 informationSeparation of ternary heteroazeotropic mixtures in a closed multivessel batch distillation decanter hybrid
Chemical Engineering and Processing 43 (2004) 291 304 Separation of ternary heteroazeotropic mixtures in a closed multivessel batch distillation decanter hybrid S. Skouras, S. Skogestad Norwegian University
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 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 informationChapter 2 Equilibria, Bubble Points, Dewpoints, Flash Calculations, and Activity Coefficients
Chapter 2 Equilibria, Bubble Points, Dewpoints, Flash Calculations, and Activity Coefficients 2.1 Vapour Pressure Calculations The basis for all phase equilibrium calculations are the vapour pressures
More informationSTARTUP ANALYSIS OF MASS- AND HEAT-INTEGRATED TWO-COLUMN-SYSTEMS
STARTUP ANALYSIS OF MASS- AND HEAT-INTEGRATED TWO-COLUMN-SYSTEMS Varbanov, P., Klein, A., Repke, J.-U., Wozny, G. Technical University Berlin, Department of Process and Plant-Dynamics, Straße des 17, Juni
More informationA 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
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... Introduce the topic of numerical methods Consider the Error analysis and sources of errors Introduction A numerical method which
More informationStructuring of Reactive Distillation Columns for Non-Ideal Mixtures using MINLP-Techniques
Structuring of Reactive Distillation olumns for Non-Ideal Mixtures using MINLP-Techniques Guido Sand *, Sabine Barmann and Sebastian Engell Process ontrol Laboratory, Dept. of Biochemical and hemical Engineering,
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 informationChapter 3 PROPERTIES OF PURE SUBSTANCES
Thermodynamics: An Engineering Approach Seventh Edition Yunus A. Cengel, Michael A. Boles McGraw-Hill, 2011 Chapter 3 PROPERTIES OF PURE SUBSTANCES Copyright The McGraw-Hill Companies, Inc. Permission
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 informationUNIVERSITY OF MINNESOTA DULUTH DEPARTMENT OF CHEMICAL ENGINEERING ChE CONTINUOUS BINARY DISTILLATION
UNIVERSITY OF MINNESOTA DULUTH DEPARTMENT OF CHEMICAL ENGINEERING ChE 3211-4211 CONTINUOUS BINARY DISTILLATION OBJECTIVE The objective of this experiment is to determine the overall column efficiency for
More informationProcess Control Exercise 2
Process Control Exercise 2 1 Distillation Case Study Distillation is a method of separating liquid mixtures by means of partial evaporation. The volatility α of a compound determines how enriched the liquid
More informationModeling and Simulation of Ethyl Acetate Reactive Distillation Column Using Aspen Plus
nternational Journal Of Scientific & Engineering Research, Volume 3, ssue 8, August-2012 1 Modeling and Simulation of Ethyl Acetate Reactive Distillation Column Using Aspen Plus Sohail Rasool Lone & Syed
More informationTeaching Assistant: Tracy Bucholz CPE office hours: M, W, F 8:30-9:30am
ChE 317: INTRODUCTION TO CHEMICAL ENGINEERING ANALYSIS T. Th 8-9:30 CPE 2.218 Personnel Contact Information Course Instructor: Prof. Lynn Loo lloo@che.utexas.edu CPE 4.422 471-6300 office hours: T, Th
More informationTHERMODYNAMIC INSIGHT ON EXTRACTIVE DISTILLATION WITH ENTRAINER FORMING NEW AZEOTROPES
Distillation Absorption 2010 A.. de Haan, H. Kooijman and A. Górak (Editors) All rights reserved by authors as per DA2010 copyright notice THERMODYNAMIC INSIGHT ON EXTRTIVE DISTILLATION WITH ENTRAINER
More information