Structuring of Reactive Distillation Columns for Non-Ideal Mixtures using MINLP-Techniques
|
|
- Maryann O’Brien’
- 6 years ago
- Views:
Transcription
1 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, Universität Dortmund, Emil-Figge-Straße 70, Dortmund, Germany Gerhard Schembecer Process Design enter B.V., J.-v.-Fraunhofer-Straße 20, Dortmund, Germany Abstract This contribution deals with the optimisation of the design of reactive distillation columns. For a highly non-ideal reacting mixture, an MINLP-model is presented and different solution approaches are discussed. Due to the efficient modelling and implementation of the superstructure and an appropriate setting of non-standard solver options, optimisation problems with up to 10,000 variables and constraints could be solved robustly without providing problem specific initial values or scaling factors. Keywords: reactive distillation, non-ideal mixtures, mathematical programming 1. Introduction A current trend in process design is towards integrated processes, i.e. the integration of reactive and separating functionalities into a single apparatus as, e.g., a reactive distillation column. ompared with the classical serial arrangement of unit operations, this advanced concept has the potential to decrease the dimensions of the equipment and to increase the degree of heat integration. Furthermore, it provides the opportunity to overcome chemical and thermodynamical boundaries, such as chemical equilibria or distillation boundaries due to azeotropes. Separations of non-ideal mixtures with simultaneous chemical reactions belong to the most difficult design problems and should be solved in an integrated fashion. The design of a reactive distillation column constitutes a constrained combinatorial optimisation problem which is amenable to MINLP-techniques. In practice, such problems often are hard to solve due to non-linear and integrality constraints. The optimisation of reactive distillation columns has previously been addressed by Grossmann and co-worers (Jacson and Grossmann, 2001) and Stichlmair and coworers (Frey and Stichlmair, 2000, Poth et al., 2003). This wor differs from the aforementioned in the modelling of the superstructure and in the applied solution methodology. A highly non-ideal chemical system is considered, and the robust solution by means of an efficient implementation of the model and a specific treatment of the integers but without a priori (heuristic) nowledge on the solution is highlighted. In * Author to whom correspondence should be adressed: guido.sand@udo.edu
2 Sect. 2, an MINLP-model for the production of methyl-tertiary-butyl-ether (MTBE) is presented. Then three algorithmic approaches to the solution of the problem without an elaborated initialisation and scaling strategy are compared. In Sect. 4, the (locally) optimal solution of the example problem is presented and discussed: It is demonstrated graphically that the chemical and thermodynamical boundaries are overcome. 2. Mixed-Integer Model of a Non-Ideal Reactive Separation Process The inetically controlled production of MTBE from isobutene and methanol (IB + MeOH MTBE) in the presence of n-butane at a pressure of 8 bar is considered here. The column is represented by a superstructure of trays each of which may be active or inactive. The model for the costing comprises revenues, running costs, and annualised investment costs. 2.1 Tray column The model of a tray is based upon the MESH-equations which are extended by source terms for the reaction (material/enthalpy balance). The equilibrium condition is extended by the Murphree efficiency for the vapour phase: E ideal ( y y ) ( y y ) i = i, i, 1 i, i, 1,. (1) The continuous variable y i, denotes the mole-fraction of the component i of the vapour stream from tray. The assumed efficiency E = 0.7 for each tray determines the ratio between the real change of the concentrations and the maximum possible change given by the thermodynamic equilibrium. The non-ideal liquid phase is modelled by activities which appear in the reaction inetics and in the thermodynamic equilibrium condition. The activity-coefficients are calculated by Wilson s approach, describing three binary azeotropes. The temperature dependencies of the vapour pressures are modelled by the Antoine equation. The formation and the decomposition of MTBE is heterogeneously catalysed and characterised by the following inetics: form eff form 1 Pr = H ai= IB, ai= MeOH, dec eff dec 2 Pr = H ai= MTBE, ai= MeOH,, (2). (3) The production rates Pr [mole/s] are calculated from variable efficient hold-ups H eff [ltr], the inetic constants and polynomials of activities a i,. The efficient hold-up is smaller than the total hold-up reflecting the amount and the efficiency of the catalyst on a tray. The temperature dependencies of the inetic constants are modelled by an extended Arrhenius approach. The thermodynamic properties and the inetics were taen from Beßling (1998) and Rehfinger and Hoffmann (1990). The geometry of the column is constrained by an upper bound on the efficient hold-up and upper and lower bounds on the vapour velocity. Both values do not appear explicitly but are bounded implicitly by the following constraints (Dia represents the column diameter and V represents the vapour flow rate from tray ):
3 eff H 2 Dia 103, 7, (4) ,7 10 V Dia 19,2 10 V. (5) The right hand side in the height-constraint Eq. (4) results from the maximum weir height of a tray and the volume specific catalyst efficiency, and the constant factors in the velocity-constraints Eq. (5) are calculated from the minimum and maximum F-factor and the mean molar volume of the vapour phase. 2.2 Superstructure The superstructure consists of N serially arranged trays, each tray has a (free) efficient hold-up H eff and (free) feeds of both raw materials. A tray can be inactivated by means of a binary variable ϕ, except for the bottom and the top tray (reboiler/condenser). The eff reactive functionality of a tray is controlled by its effective hold-up H IR + (0 means no reaction). The separating functionality is controlled by switching its efficiency E IR + by means of the activation variable ϕ : E 0.7 =ϕ. (6) For deactivated trays, the efficient hold-up and the feed are forced to zero by big-m inequalities such that the vapour and the liquid streams pass a deactivated tray unaffected. Additionally, the sum of the feeds to the trays must be equal to the two given total feeds of methanol and isobutene/n-butane. Ambiguous activation patterns are avoided by the logic constraint: ϕ + 1 ϕ 0 < N 1. (7) Note, that the binary variables appear linearly throughout. 2.3 ost model The column is designed to produce MTBE with a purity of 99% (liquid mole fraction x =1 MTBE 0.99) with minimum annual cost ann [ /a]. These costs are calculated from the annualised investment costs for the shell shell, the internals inter, the reboiler reb, the condenser cond and the catalyst cat, in addition to running costs for cooling and heating (Q =n, Q =1 [W]) and the feed streams (F MeOH, F IB/but [mole/s]), reduced by revenues for the bottom and top distillate streams (D =1, D =N MeOH, D =N IB/but [mole/s]). ann = 13.53H Dia ϕ Q 825.1D shell reb = 1 = = = Q 138.4D = N MeOH = N 2 3 ( Dia Dia 0.021Dia ) cond T Q T Q ,6 ln, = N F 410.8D 0.6 = N MeOH IB/but = N inter cat F 3 eff H IB/but (8)
4 H [m] denotes the column height, and T =1 [K] and T ln, =N [K] denote the (logarithmic) temperature differences of the reboiler and the condenser, respectively. The objective function is taen from ooboo (2003). 2.4 Efficient implementation Preliminary tests of different NLP-solvers (cf. Sect. 4) indicated that ONOPT3 is the most promising choice. The difficulties in the algorithmic solution process can be reduced by efficiently implementing the model equations (see rules of thumb in ONOPT, 2003): Degeneracy as well as very small or large function values and gradients should be avoided, and highly non-linear expressions should be re-written. To this end, unnecessary bounds are removed (e.g. bounds for activities are not stated explicitly but implicitly by means of Wilson s approach), algorithmic bounds are added (e.g. Q =1 1 to avoid Q =1 = 0) and auxiliary variables are introduced to simplify nonlinear functions of expressions (e.g. the reformulation of the Arrhenius approach: e -E/(RT) e aux. aux T = -E/R). A specific type of algorithmic bounds is introduced into Eqs. (2) and (3) to prevent the activities a i=meoh, from assuming a value of zero due to insufficient numerical precision: Instead of imposing small lower bounds on a i=meoh, (which would be valid globally), a small value ε = 10-2 is introduced into Eqs. (2) and (3) only, they are implemented as Pr Pr form eff form ( ai= MeOH, + ε) = H ai= IB, 2 eff dec ( a + ε) = H a dec i= MeOH, i= MTBE, 3. Solution Algorithms, (9). (10) The solution procedure is based on the decomposition of the MINLP problem into an IP-master-problem (optimisation of the number of trays) and NLP-sub-problems (optimisation of continuous variables for a fixed number of trays). For the IP-part three algorithmic approaches were compared; the NLP-sub-problems are always solved by the generalized reduced gradient method implemented in ONOPT NLP-sub-problems ONOPT was adjusted to the problem at hand by iteratively setting non-standard options based on evaluations of the results of tests without applying any heuristic nowledge on good solutions. This procedure ensures a conceptual partition of the model from the solution procedure. All variables were initialised by 1 and no scaling factors were supplied. To chec the reproducibility of the results, 91 problems with different numbers of trays N (10 N 100) were solved. The models comprise 1,021 10,111 variables (8 98 binaries) and 1,051 10,411 constraints. With the standard solver options, the robustness turned out to be insufficient: Only 35% of the problems were solved to (local) optimality. This value was increased to 95% by decreasing the feasibility tolerance (rtnwma = 10-6 ), increasing the maximum number of stalled iterations (lfstal = 10 5 ), by the introduction of slac variables for the constraints in the first iteration step (lslac = true), by optimising the step-length in Newton iterations
5 (lmmxsf = 1) and by automatically scaling the variables and constraints in each iteration step (lfscal = 1). The solution times were 6 1,779 PU-s (Windows-P, 1,5 GHz). 3.2 IP-master-problems The first approach for the IP-part was the branch and bound-algorithm implemented in SBB (2003). The integrality requirements are totally relaxed in the root node and reinforced for one variable after the other from layer to layer. At each node of the search tree, a lower bound on the (locally) optimal solution is generated by solving an NLP (integer-relaxed MINLP). For a superstructure with N = 100 trays, a (locally) optimal solution was found after the exploitation of 58 nodes (PU-time 44.5 min). The annual cost is ann = 1,050 /a and 49 trays are active. Note that in the worst case the number of NLP-sub-problems that has to be solved is in the order of O(2 N ). Secondly, a problem specific approach was applied which completely enumerates the subspace of binary variables ϕ which is feasible wrt. Eq. (7). In fact, this is exactly the same as a solution of problems for different numbers of trays as treated in Sect The solution (identical with the first one) was obtained after solving 91 NLPs. However, it should be noted that in general the number of NLPs to be solved is only in the order of O(N). Under the assumption that for a fixed number of trays the minimum annual cost ann is unimodal wrt. the number of trays, the optimal number of trays can be found efficiently by an interval reduction algorithm: The search for the optimal number of trays is considered as a one-dimensional unconstrained optimisation problem, which must be solved by direct methods (without gradient information). Initially, two problems with N 0 = {10;100} were solved under the assumption that the optimal solution is within the interval [10;100]. The range of the interval is reduced iteratively by a factor of 1/3 as follows: Two new NLPs with N 1 = {40;70} are solved and the corresponding values ann are compared. From ann (40) < ann (70) it follows that the optimum must be in the interval [10;70]. This procedure is repeated iteratively until the interval range is 1, i.e. a local optimum is obtained. Here, the number of NLPs to be solved was only 17 (identical solution as above), in general it is only in the order of O(lg 3/2 N). In comparison, in the disjunctive programming approach presented by Jacson and Grossmann (2001), the number of NLPs to be solved is up to two orders of magnitude larger: For a model with 38 binary variables they solved 1,270 NLPs. Butane Feed MeOH MeOH MTBE Feed IB/Butane IB Azeotropes Separating Trays Integrated Trays Figure 1: Liquid mole fractions, chemical equilibrium lines and distillation boundary
6 4. Optimal Structure The (locally) optimal column has 49 trays with separating functionality with reactive functionality on trays Feed trays (with non-marginal feed) are for MeOH and for IB/butane. The diameter is Dia = 0.39 m and the reflux ratio is 4,5. From the graphical representation of the liquid mole fractions on the trays in Figure 1, it can be seen that the boundaries given by the chemical equilibrium and the distillation boundary (given by the azeotropes) are overcome. 5. Summary and Perspectives An MINLP-model for the optimal structuring of a reactive distillation column for the synthesis of MTBE was presented. The superstructure is modelled by linear constraints and the highly non-linear equations are implemented efficiently. Three approaches to solve the IP-master-problems which require a relatively small number of NLP-subproblems to be solved were presented. urrent research deals with the systematic integration of heuristic nowledge into the model and the algorithm in order to increase the solution efficiency while maintaining its robustness. Acnowledgement The research reported here was supported by the Deutsche Forschungsgemeinschaft as part of the research unit Integrated Reaction and Separation Processes at the Universität Dortmund. References Beßling, B., 1998, Zur Reativdestillation in der Prozeßsynthese, Dr.-Ing.-Dissertation, Fachbereich hemietechni, Universität Dortmund (in German). ONOPT, 2003, ooboo, 2003, Process Design enter B.V., Breda, The Netherlands. Frey, T. and J. Stichlmair, 2000, MINLP Optimization of Reactive Distillation olumns, in: ESAPE 10, omputer Aided hemical Engineering 8, , Ed. S. Pierucci, Elsevier, Amsterdam. Jacson, J.R. and I.E. Grossmann, 2001, A Disjunctive Programming Approach for the Optimal Design of Reactive Distillation olumns, omp. hem. Engg. 25, Poth, N., D. Brusis and J. Stichlmair, 2003, Rigorous Optimization of Reactive Distillation in GAMS with the Use of External Functions, in: ESAPE 13, omputer Aided hemical Engineering 14, , Ed. A. Kraslawsi and I. Turunen, Elsevier, Amsterdam. Rehfinger, A. and U. Hoffmann, 1990, Kinetics of Methyl Tertiary Butyl Ether Liquid Phase Synthesis atalyzed by Ion Exchange Resin I. Intrinsic Rate Expression in Liquid Phase Activities, hem. Eng. Sci., 45, SBB, 2003, Users Manual, ARKI onsulting & Development,
Optimal Design of a Reactive Distillation Column
Optimal Design of a Reactive Distillation Column Edwin Zondervan, Mayank Shah * and André B. de Haan Eindhoven University of Technology, Department of Chemistry and Chemical Engineering, P.O. Box 513,
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 informationModel-based Conceptual Design and Tool Support for the Development of Continuous Chemical Processes
Ian David Lockhart Bogle and Michael Fairweather (Editors), Proceedings of the 22nd European Symposium on Computer Aided Process Engineering, 17-20 June 2012, London. 2012 Elsevier B.V. All rights reserved
More informationOptimal structure of reactive and non-reactive stages in reactive distillation processes
Optimal structure of reactive and non-reactive stages in reactive distillation processes Alejandro O. Domancich 1,2, Nélida B. Brignole 1,2, Patricia M. Hoch 2,3 1 LIDeCC - Departamento de Ciencias e Ingeniería
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 informationExperimental and Simulation Study on the Reactive Distillation Process for the Production of Ethyl Acetate
Experimental and Simulation Study on the Reactive Distillation Process for the Production of Ethyl Acetate Jongkee Park, Na-Hyun Lee, So-Jin Park, and Jungho Cho, Separation Process Research Center, Korea
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 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 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 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 informationStructural Multiplicity and Redundancy in Chemical Process Synthesis with MINLP
Structural Multiplicity and Redundancy in Chemical Process Synthesis with MINLP Tivadar Farkas, Endre Rev and Zoltan Lelkes * Chemical Engineering Department, Budapest University of Technology and Economics
More informationBOUNDARY VALUE DESIGN METHOD FOR COMPLEX DEMETHANIZER COLUMNS
Distillation Absorption 2010 A.B. de Haan, H. Kooijman and A. Górak (Editors) All rights reserved by authors as per DA2010 copyright notice BOUNDARY AUE DESIGN METHOD FOR COMPEX DEMETHANIZER COUMNS Muneeb
More 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 informationDesign and Control Properties of Arrangements for Distillation of Four Component Mixtures Using Less Than N-1 Columns
D. M. MÉNDEZ-VALENCIA et al., Design and Control Properties of Arrangements, Chem. Biochem. Eng. Q. 22 (3) 273 283 (2008) 273 Design and Control Properties of Arrangements for Distillation of Four Component
More informationShort-cut distillation columns based on orthogonal polynomials of a discrete variable
Short-cut distillation columns based on orthogonal polynomials of a discrete variable Peter Lory Institut für Wirtschaftsinformatik, Universität Regensburg, D 93040 Regensburg, Germany Dedicated to Professor
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 informationProcess Integration Methods
Process Integration Methods Expert Systems automatic Knowledge Based Systems Optimization Methods qualitative Hierarchical Analysis Heuristic Methods Thermodynamic Methods Rules of Thumb interactive Stochastic
More informationDEVELOPMENT OF A ROBUST ALGORITHM TO COMPUTE REACTIVE AZEOTROPES
Brazilian Journal of Chemical Engineering ISSN 0104-6632 Printed in Brazil www.abeq.org.br/bjche Vol. 23, No. 03, pp. 395-403, July - September, 2006 DEVELOPMENT OF A ROBUST ALGORITHM TO COMPUTE REACTIVE
More informationThermally Coupled Distillation Systems: Study of an Energy-efficient Reactive Case
F. O. BARROSO-MUÑOZ et al., Thermally Coupled Distillation Systems: Study of, Chem. Biochem. Eng. Q. 21 (2) 115 120 (2007) 115 Thermally Coupled Distillation Systems: Study of an Energy-efficient Reactive
More informationStudy of arrangements for distillation of quaternary mixtures using less than n-1 columns
Instituto Tecnologico de Aguascalientes From the SelectedWorks of Adrian Bonilla-Petriciolet 2008 Study of arrangements for distillation of quaternary mixtures using less than n-1 columns J.G. Segovia-Hernández,
More informationOptimal Synthesis of Complex Distillation ColumnsUsing Rigorous Models
Optimal Synthesis of Complex Distillation ColumnsUsing Rigorous Models Ignacio E. Grossmann 1, Pío A. Aguirre 2 and Mariana Barttfeld 2 1 Department of Chemical Engineering, Carnegie Mellon University,
More informationStructural considerations and modeling in the synthesis of heat integrated thermally coupled distillation sequences
Structural considerations and modeling in the synthesis of heat integrated thermally coupled distillation sequences José. aballero* ; Ignacio E. Grossmann ** *epartment of hemical Engineering. University
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 informationAuthor's personal copy
Computers and Chemical Engineering 33 (2009) 788 793 Contents lists available at ScienceDirect Computers and Chemical Engineering journal homepage: www.elsevier.com/locate/compchemeng Impact of mathematical
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 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 informationDETERMINATION OF OPTIMAL ENERGY EFFICIENT SEPARATION SCHEMES BASED ON DRIVING FORCES
DETERMINATION OF OPTIMAL ENERGY EFFICIENT SEPARATION SCHEMES BASED ON DRIVING FORCES Abstract Erik Bek-Pedersen, Rafiqul Gani CAPEC, Department of Chemical Engineering, Technical University of Denmark,
More 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 informationCASE-BASED REASONING FOR SEPARATION PROCESS SYNTHESIS
CASE-BASED REASONING FOR SEPARATION PROCESS SYNTHESIS Timo Seuranen*, Elina Pajula**, Markku Hurme* * Helsinki University of Technology, Laboratory of Plant Design, P.O. Box 6100, FIN-02015 HUT, Finland
More informationSynthesis and Design of Reactive Distillation Columns
Synthesis and Design of Reactive Distillation Columns thesis submitted to the University of Manchester Institute of Science and Technology for the degree of Doctor of Philosophy by Ramona Manuela DRGOMIR
More informationImproved Big-M Reformulation for Generalized Disjunctive Programs
Improved Big-M Reformulation for Generalized Disjunctive Programs Francisco Trespalacios and Ignacio E. Grossmann Department of Chemical Engineering Carnegie Mellon University, Pittsburgh, PA 15213 Author
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 informationControl Study of Ethyl tert-butyl Ether Reactive Distillation
3784 Ind. Eng. Chem. Res. 2002, 41, 3784-3796 Control Study of Ethyl tert-butyl Ether Reactive Distillation Muhammad A. Al-Arfaj Department of Chemical Engineering, King Fahd University of Petroleum and
More informationDesign and Optimisation of Batch Reactors
Design and Optimisation of Batch Reactors Jinzhong Zhang Supervisor: Professor Robin Smith 2-1 XVII PIRC Annual Research Meeting 2000 In recent years, a methodology has been developed for the systematic
More informationMixed-Integer Nonlinear Decomposition Toolbox for Pyomo (MindtPy)
Mario R. Eden, Marianthi Ierapetritou and Gavin P. Towler (Editors) Proceedings of the 13 th International Symposium on Process Systems Engineering PSE 2018 July 1-5, 2018, San Diego, California, USA 2018
More 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 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 informationSynthesis of separation processes by using case-based reasoning
Computers and Chemical Engineering 25 (2001) 775 782 www.elsevier.com/locate/compchemeng Synthesis of separation processes by using case-based reasoning Elina Pajula *, Timo Seuranen, Tuomas Koiranen,
More informationSimplified Methodology for the Design and Optimization of Thermally Coupled Reactive Distillation Systems
pubs.acs.org/iecr Simplified Methodology for the Design and Optimization of Thermally Coupled Reactive Distillation Systems Fernando Israel Goḿez-Castro,, Vicente Rico-Ramírez,*, Juan Gabriel Segovia-Hernańdez,
More informationPRELIMINARY DESIGN OF REACTIVE DISTILLATION COLUMNS
PRELIMINARY DESIGN OF REACTIVE DISTILLATION COLUMNS R. THERY, X. M. MEYER, X. JOULIA and M. MEYER Process System Engineering Department, Laboratoire de Génie Chimique INPT-ENSIACET, Toulouse, France A
More informationFinding Rigorous Thermodynamic Model parameters for the Simulation of a Reactive Distillation Column: Case of ETBE Production
Finding Rigorous Thermodynamic Model parameters for the Simulation of a Reactive Distillation Column: Case of ETBE Production Gisela N. Durruty, M. Soledad Diaz, Patricia M. Hoch * Universidad Nacional
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 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 informationA novel design of reactive distillation configuration for 2-methoxy-2-methylheptane process
A novel design of reactive distillation configuration for 2-methoxy-2-methylheptane process Arif Hussain 1, Muhammad Abdul Qyyum 1, Le Quang Minh 1, Hong Jimin 1, and Moonyong Lee 1,* 1 Process System
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 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 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 informationImpact of mass transfer coefficient correlations on prediction of reactive distillation column behaviour
Chemical Engineering Journal 140 (2008) 381 390 Impact of mass transfer coefficient correlations on prediction of reactive distillation column behaviour Z. Švandová, J. Markoš,L. Jelemenský Institute of
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 informationOn the Design of Optimal Solvent Mixtures using Generalised Disjunctive Programming
Ian David Lockhart Bogle and Michael Fairweather (Editors), Proceedings of the 22nd European Symposium on Computer Aided Process Engineering, 17-20 June 2012, London. 2012 Elsevier B.V. All rights reserved.
More informationMultiperiod Blend Scheduling Problem
ExxonMobil Multiperiod Blend Scheduling Problem Juan Pablo Ruiz Ignacio E. Grossmann Department of Chemical Engineering Center for Advanced Process Decision-making University Pittsburgh, PA 15213 1 Motivation
More informationNovel Control Structures for Heterogeneous Reactive Distillation
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 informationINTRODUCTION TO CHEMICAL PROCESS SIMULATORS
INTRODUCTION TO CHEMICAL PROCESS SIMULATORS DWSIM Chemical Process Simulator A. Carrero, N. Quirante, J. Javaloyes October 2016 Introduction to Chemical Process Simulators Contents Monday, October 3 rd
More informationSoftware for Integer and Nonlinear Optimization
Software for Integer and Nonlinear Optimization Sven Leyffer, leyffer@mcs.anl.gov Mathematics & Computer Science Division Argonne National Laboratory Roger Fletcher & Jeff Linderoth Advanced Methods and
More informationOptimal Operation of Batch Reactive Distillation Process Involving Esterification Reaction System
1387 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 43, 2015 Chief Editors: Sauro Pierucci, Jiří J. Klemeš Copyright 2015, AIDIC Servizi S.r.l., ISBN 978-88-95608-34-1; ISSN 2283-9216 The Italian
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 informationbased on crystallization technology
Optimal synthesis of p-xylene separation processes based on crystallization technology Ricardo M. Lima, Ignacio E. Grossmann Department of Chemical Engineering, Carnegie Mellon University Pittsburgh, PA
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 informationOptimization of Heat Exchanger Network with Fixed Topology by Genetic Algorithms
Optimization of Heat Exchanger Networ with Fixed Topology by Genetic Algorithms Roman Bochene, Jace Jeżowsi, Sylwia Bartman Rzeszow University of Technology, Department of Chemical and Process Engineering,
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 informationReal-Time Feasibility of Nonlinear Predictive Control for Semi-batch Reactors
European Symposium on Computer Arded Aided Process Engineering 15 L. Puigjaner and A. Espuña (Editors) 2005 Elsevier Science B.V. All rights reserved. Real-Time Feasibility of Nonlinear Predictive Control
More informationA Short Method To Calculate Residue Curve Maps in Multireactive and Multicomponent Systems
pubs.acs.org/iecr A Short Method To Calculate Residue Curve Maps in Multireactive and Multicomponent Systems Marcelino Carrera-Rodríguez, Juan Gabriel Segovia-Hernandez,*, and Adrian Bonilla-Petriciolet
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 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 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 informationREACTIVE distillation is a separation process that combines
Simulation of Reactive Distillation: Comparison of Equilibrium and Nonequilibrium Stage Models Asfaw Gezae Daful Abstract In the present study, two distinctly different approaches are followed for modeling
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 informationOPTIMAL SOLVENT CYCLE DESIGN IN SUPERCRITICAL FLUID PROCESSES
Latin American Applied Research 33:11-15 (3) OPTIMAL SOLVENT CYCLE DESIGN IN SUPERCRITICAL FLUID PROCESSES M. S. DIAZ, S. ESPINOSA and E. A. BRIGNOLE Planta Piloto de Ingenieria Quimica - PLAPIQUI (UNS-CONICET)
More informationRecovery of Aromatics from Pyrolysis Gasoline by Conventional and Energy-Integrated Extractive Distillation
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 Recovery of Aromatics from Pyrolysis Gasoline by Conventional
More informationTitle: Logic hybrid simulation-optimization algorithm for distillation design. Computers and Chemical Engineering
Title: Logic hybrid simulation-optimization algorithm for distillation design. Author: J.A. Caballero PII: S0098-1354(14)00097-0 DOI: http://dx.doi.org/doi:10.1016/j.compchemeng.2014.03.016 Reference:
More informationInnovative Design of Diphenyl Carbonate Process via One Reactive Distillation with a Feed-Splitting Arrangement
445 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 61, 2017 Guest Editors: Petar S Varbanov, Rongxin Su, Hon Loong Lam, Xia Liu, Jiří J Klemeš Copyright 2017, AIDIC Servizi S.r.l. ISBN 978-88-95608-51-8;
More informationRobust Process Control by Dynamic Stochastic Programming
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany MARC C. STEINBACH Robust Process Control by Dynamic Stochastic Programming ZIB-Report 04-20 (May 2004) ROBUST
More informationREACTIVE DIVIDING-WALL COLUMNS: TOWARDS ENHANCED PROCESS INTEGRATION
Distillation bsorption 1.. de aan,. Kooijman and. Górak (Editors) ll rights reserved by authors as per D1 copyright notice RETIVE DIVIDING-WLL OLUMNS: TOWRDS ENNED PROESS INTEGRTION nton. Kiss, J. J. Pragt,.
More informationVapor-liquid Separation Process MULTICOMPONENT DISTILLATION
Vapor-liquid Separation Process MULTICOMPONENT DISTILLATION Outline: Introduction to multicomponent distillation Phase Equilibria in Multicomponent Distillation (Pg. 737) Bubble-point and dew-point calculation
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 informationSynthesis of Azeotropic Separation Systems by Case-Based Reasoning
Synthesis of Azeotropic Separation Systems by Case-Based Reasoning Timo Seuranen 1, Elina Pajula 2, Markku Hurme 1 1 Helsinki University of Technology, Laboratory of Plant Design, P.O. Box 6100, FIN-02015
More informationIndicator Constraints in Mixed-Integer Programming
Indicator Constraints in Mixed-Integer Programming Andrea Lodi University of Bologna, Italy - andrea.lodi@unibo.it Amaya Nogales-Gómez, Universidad de Sevilla, Spain Pietro Belotti, FICO, UK Matteo Fischetti,
More informationReduction of Energy Consumption and Greenhouse Gas Emissions in a Plant for the Separation of Amines
1462 Chem. Eng. Technol. 28, 31, No. 1, 1462 1469 Dionicio Jantes-Jaramillo 1 Juan Gabriel Segovia- Hernández 1 Salvador Hernández 1 1 Universidad de Guanajuato, Facultad de Química, Guanajuato, Mexico.
More informationDesign and Optimization of Thermally Coupled Distillation Schemes for the Separation of Multicomponent Mixtures
724 Ind. Eng. Chem. Res. 2006, 45, 724-732 Design and Optimization of Thermally Coupled Distillation Schemes for the Separation of Multicomponent Mixtures Christopher Jorge Calzon-McConville, Ma. Bibiana
More informationDevelopment of a Systematic Synthesis and Design Methodology to Achieve Process Intensification in (Bio)Chemical Processes
Development of a Systematic Synthesis and Design Methodology to Achieve Process Intensification in (Bio)Chemical Processes Philip Lutze a, Rafiqul Gani b, John M. Woodley a a PROCESS, Center for Process
More informationMixed Integer Programming Solvers: from Where to Where. Andrea Lodi University of Bologna, Italy
Mixed Integer Programming Solvers: from Where to Where Andrea Lodi University of Bologna, Italy andrea.lodi@unibo.it November 30, 2011 @ Explanatory Workshop on Locational Analysis, Sevilla A. Lodi, MIP
More informationPressure Swing Distillation with Aspen Plus V8.0
Pressure Swing Distillation with Aspen Plus V8.0 1. Lesson Objectives Aspen Plus property analysis RadFrac distillation modeling Design Specs NQ Curves Tear streams Understand and overcome azeotrope Select
More informationSeparation Trains Azeotropes. S,S&L Chapter 9.5 Terry A. Ring Chemical Engineering University of Utah
Separation Trains Azeotropes S,S&L Chapter 9.5 Terry A. Ring Chemical Engineering University of Utah Distillation Recycle Loops Closing Recycle Loops Matrix Mathematics Without Recycle Loop [P] x = y Straight
More informationA Systematic Approach towards Synthesis and Design of Pervaporation-Assisted Separation Processes
1534 Research Article Chemie A Systematic Approach towards Synthesis and Design of Pervaporation-Assisted Separation Processes Bettina Scharzec, Thomas Waltermann, and Mirko Skiborowski* DOI: 10.1002/cite.201700079
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 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 informationAzeotropic Distillation Methods. Dr. Stathis Skouras, Gas Processing and LNG RDI Centre Trondheim, Statoil, Norway
Azeotropic Distillation Methods Dr. Stathis Skouras, Gas Processing and LNG RDI Centre Trondheim, Statoil, Norway Schedule Tuesday 1/12/2015: 09.45 12.30: Lecture - Natural Gas Processing 14.00 17.00:
More informationBatchExtractiveDistillationwithLightEntrainer
V. VARGA et al., Batch Extractive Distillation with Light Entrainer, Chem. Biochem. Eng. Q. 20 (1) XXX XXX (2006) 1 BatchExtractiveDistillationwithLightEntrainer V. Varga,*, ** E. Rev,*,+ V. Gerbaud,**
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 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 informationOptimal (Solvent) Mixture Design through a Decomposition Based CAMD methodology
European Symposium on Computer-Aided Process Engineering- 14 A. Barbosa-P6voa and H. Matos (Editors) 9 2004 Elsevier B.V. All rights reserved. 217 Optimal (Solvent) Mixture Design through a Decomposition
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 informationLagrangean Decomposition for Mean-Variance Combinatorial Optimization
Lagrangean Decomposition for Mean-Variance Combinatorial Optimization Frank Baumann, Christoph Buchheim, and Anna Ilyina Fakultät für Mathematik, Technische Universität Dortmund, Germany {frank.baumann,christoph.buchheim,anna.ilyina}@tu-dortmund.de
More informationDevelopment of an algorithm for solving mixed integer and nonconvex problems arising in electrical supply networks
Development of an algorithm for solving mixed integer and nonconvex problems arising in electrical supply networks E. Wanufelle 1 S. Leyffer 2 A. Sartenaer 1 Ph. Toint 1 1 FUNDP, University of Namur 2
More informationHeterogeneous Azeotropic Distillation Operational Policies and Control
Heterogeneous Azeotropic Distillation Operational Policies and Control Claudia J. G. Vasconcelos * and Maria Regina Wolf-Maciel State University of Campinas, School of Chemical Engineering, Campinas/SP,
More informationOptimization of Batch Processes
Integrated Scheduling and Dynamic Optimization of Batch Processes Yisu Nie and Lorenz T. Biegler Center for Advanced Process Decision-making Department of Chemical Engineering Carnegie Mellon University
More informationCork Institute of Technology. Autumn 2005 CE 2.7 Separation Process & Particle Technology (Time: 3 Hours) Section A
Cork Institute of Technology Bachelor of Engineering (Honours) in Chemical and Process Engineering Stage 2 (Bachelor of Engineering in Chemical and Process Engineering Stage 2) (NFQ Level 8) Autumn 2005
More informationThe Role of Process Integration in Process Synthesis
The Role of Process Integration in Process Synthesis Jeffrey J. Siirola Purdue University, West Lafayette, Indiana Carnegie Mellon University, Pittsburgh, Pennsylvania "Process integration" and "process
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 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 information