The simultaneous prediction of vapor-liquid equilibrium and excess enthalpy. Kwon, Jung Hun. Thermodynamics and properties lab.

Similar documents
A modification of Wong-Sandler mixing rule for the prediction of vapor-liquid equilibria in binary asymmetric systems

CH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions

EQUATION OF STATE DEVELOPMENT

Status and results of group contribution methods

Generalized binary interaction parameters in the Wong Sandler mixing rules for mixtures containing n-alkanols and carbon dioxide

Vapour Liquid Equilibrium in Asymmetric Mixtures of n-alkanes with Ethane

Chemical and Engineering Thermodynamics

Calorimetry in the near-critical and supercritical regions. Nitrous oxide + hydrocarbon mixtures*

MODELING OF PHASE EQUILIBRIA FOR BINARY AND TERNARY MIXTURES OF CARBON DIOXIDE, HYDROGEN AND METHANOL

510 Subject Index. Hamiltonian 33, 86, 88, 89 Hamilton operator 34, 164, 166

Modified Raoult's Law and Excess Gibbs Free Energy

CH2351 Chemical Engineering Thermodynamics II Unit I, II Phase Equilibria. Dr. M. Subramanian

PETE 310 Lectures # 36 to 37

Goal. Prerequisites. Performance Objectives. chenvle2.mcd Page 1

Overall: 75 ECTS: 7.0

CHAPTER SIX THERMODYNAMICS Vapor-Liquid Equilibrium in a Binary System 6.2. Investigation of the Thermodynamic Properties of Pure Water

Equations of State. Equations of State (EoS)

CONTENTS. Notes to Students Acknowledgments ABOUT THE AUTHORS UNIT I FIRST AND SECOND LAWS 1

Prediction of phase equilibria in waterralcoholralkane systems

Modeling Vapor Liquid Equilibrium of Binary and Ternary Systems of CO 2 + Hydrocarbons at High-Pressure Conditions

Prediction of Vapor/Liquid Equilibrium Behavior from Quantum Mechanical Data

2. Derive ideal mixing and the Flory-Huggins models from the van der Waals mixture partition function.

Optimization of the Sulfolane Extraction Plant Based on Modeling and Simulation

Cooling Temperatures of Binary Mixed Refrigerants: Vapor-Liquid-Liquid Equilibrium versus Vapor-Liquid Equilibrium

Comparison of different mixing rules for prediction of density and residual internal energy of binary and ternary Lennard Jones mixtures

Agung Ari Wibowo, S.T, M.Sc THERMODYNAMICS MODEL

Numerical Aspects of the SAFT Equation of State

Modeling of the solubility of Naproxen and Trimethoprim in different solvents at different temperature

COSMO-RS Theory. The Basics

Virial Expansion. Silicon Valley FIG March 24, Chen-Hanson Ting

Vapor liquid equilibria of carbon dioxide with diethyl oxalate, ethyl laurate, and dibutyl phthalate binary mixtures at elevated pressures

A new segmental local composition model for calculation of thermodynamic properties of binary polymer solutions

New correlation for hydrogen-natural gas mixture compressibility factor

Measurement and Calculation of Physico-Chemical Properties of Binary Mixtures Containing Xylene and 1- Alkanol

Chapter 2 Equilibria, Bubble Points, Dewpoints, Flash Calculations, and Activity Coefficients

Global Phase Diagrams and Critical Phenomena of Binary Mixtures. Ji Lin Wang

An augmented hard core equation of state generalized in terms of T,, P, and o

Excess molar volume of acetonitrile + alcohol systems at K. Part II: Correlation by cubic equation of state

ME 201 Thermodynamics

Overview. Types of Solutions. Intermolecular forces in solution. Concentration terms. Colligative properties. Osmotic Pressure 2 / 46

POLYMER SOLUTIONS AND SOLVENT-SUPERCRITICAL CARBON DIOXIDE STUDIES APPLYING PC-SAFT

Chem 112 Dr. Kevin Moore

P1: IML/FFX P2: IML/FFX QC: IML/FFX T1: IML AT029-FM AT029-Manual AT029-Manual-v8.cls December 11, :59. Contents

Chemical, Biochemical, and Engineering Thermodynamics

= = 10.1 mol. Molar Enthalpies of Vaporization (at Boiling Point) Molar Enthalpy of Vaporization (kj/mol)

Solubility of Supercritical CO 2 in Polystyrene during Foam Formation via Statistical Associated Fluid Theory (SAFT) Equation of State

Solubility of solids in supercritical fluids using equations of state - excess Gibbs free energy models.

Vapor liquid equilibrium of carbon dioxide with ethyl caproate, ethyl caprylate and ethyl caprate at elevated pressures

Polymath Examples: Nonlinear Algebraic Equation and Regression Problems


RESEARCH NOTE. DOUGLAS HENDERSON? Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah USA

On the Boyle temperature

4.3 CONCLUSION: HOW TO CHOOSE A MODEL The right questions Ionic liquids What is the property of interest?

PHASE EQUILIBRIUM CALCULATIONS OF HIGHLY POLAR SYSTEMS

Vapor liquid equilibria for the binary system 2,2 dimethylbutane + 1,1 dimethylpropyl methyl ether (TAME) at , , and 338.

PRACTICAL DATA CORRELATION OF FLASHPOINTS OF BINARY MIXTURES BY A RECIPROCAL FUNCTION: THE CONCEPT AND NUMERICAL EXAMPLES

THE PROPERTIES OF GASES AND LIQUIDS

A crossover equation of state for associating fluids

Copyright Warning & Restrictions

Fluid data consistency testing incorporation into the ThermoDataEngine

The successes and limitations of this thermodynamic model in predicting the solubility of organic solutes in pure solvents

물 - 에탄올공비증류공정의최적화 공주대학교화학공학부조정호

Outline of the Course

ABSTRACT. A Systematic Study of Short and Long Range Interactions in Associating Fluids Using Molecular Theory. Wael A. Fouad

Simple Mixtures. Chapter 7 of Atkins: Section

Chapter 10. Vapor/Liquid Equilibrium: Introduction

Chemistry 111 Syllabus

A Thermodynamic Study of Methane Hydrates Formation In Glass Beads

Modeling and Computation of Phase Equilibrium. Using Interval Methods. Department of Chemical and Biomolecular Engineering, University of Notre Dame

MATSCI 204 Thermodynamics and Phase Equilibria Winter Chapter #4 Practice problems

Exam 3 Solutions. ClO g. At 200 K and a total pressure of 1.0 bar, the partial pressure ratio for the chlorine-containing compounds is p ClO2

Intermolecular Forces, Liquids, Solids. IM Forces and Physical Properties

CALCULATION OF THE COMPRESSIBILITY FACTOR AND FUGACITY IN OIL-GAS SYSTEMS USING CUBIC EQUATIONS OF STATE

Mutual diffusion in the ternary mixture of water + methanol + ethanol: Experiments and Molecular Simulation

Chemical Engineering Thermodynamics

SOFTWARE INTELIGENT PACKAGE FOR PHASE EQULIBRIA (PHEQ) IN SYSTEMS APPLIED IN CHEMISTRY AND CHEMICAL ENGINEERING

Note: items marked with * you should be able to perform on a closed book exam. Chapter 10 Learning Objective Checklist

Introduction (1) where ij denotes the interaction energy due to attractive force between i and j molecules and given by; (2)

Finitely Concentrated Partial Molar Excess Properties of Solvent/Polymer [poly(4-methylstyrene) (PMS), poly(vinylbenzyl chloride) (PVBC)] Systems

PVTpetro: A COMPUTATIONAL TOOL FOR ISOTHERM TWO- PHASE PT-FLASH CALCULATION IN OIL-GAS SYSTEMS

A perturbed hard sphere chain equation of state for applications to hydrofluorocarbons, hydrocarbons and their mixtures

CH.8 Polymers: Solutions, Blends, Membranes, and Gels

arxiv: v1 [cond-mat.soft] 20 Apr 2015

ChBE BIBLE. Robert A. Pinnick. 28 April 2006

- L. _/cit ...B N - Q. - PAD"t..t... QatJ11.81L. ><"'g. Below is a blank flow diagram of a generalized distillation column: (10 Points)

INTERFACIAL PROPERTIES AND PHASE EQUILIBRIA FOR MIXTURES RELEVANT IN THE OIL AND GAS INDUSTRY

Surface Tension Prediction for Liquid Mixtures

Week 8 Intermolecular Forces

Prediction of methanol content in natural gas with the GC-PR-CPA model Hajiw, Martha; Chapoy, Antonin; Coquelet, Christophe; Lauermann, Gerhard

Model Framework of predictive PR78 Validation Results of ppr78 EOS Discussion Conclusion for Validation Work Further Improvements

Evaluation of the SPUNG Equation of State for use in Carbon Capture and Storage Modelling

Imperfect Gases. NC State University

Vapor liquid equilibria of carbon dioxide with ethyl benzoate, diethyl succinate and isoamyl acetate binary mixtures at elevated pressures

Py x P P P. Py x P. sat. dq du PdV. abs Q S. An Innovative Approach in the G U TS PV P P G U TS PV T H U PV H U PV. abs. Py x P. sat.

Chapter 2 Experimental sources of intermolecular potentials

Mutual diffusion of binary liquid mixtures containing methanol, ethanol, acetone, benzene, cyclohexane, toluene and carbon tetrachloride

Index to Tables in SI Units

Volumetric, Viscometric and Refractive Indices Properties of Binary Mixtures of Acetyl Acetone with 1-Butanol at Different Temperatures

,, Seong-Bo Kim,Hai-SongBae, and Jeong-Sik Han

A Generalized Correlation for Pool Boiling Heat Transfer Coefficient Based on Corresponding State Rule for Pure Compounds and Binary Mixtures

Transcription:

The simultaneous prediction of vapor-liquid equilibrium and excess enthalpy Kwon, Jung Hun.

2 Contents 1 A comparison of cubic EOS mixing rules for the simultaneous description of excess enthalpies and VLE Hasan Orbey, Stanley I. Sandler (1996) 2 Calculation of excess enthalpies with the CCOR EOS H.-M Lin, M.-J Lee (1994) 3 A new method to predict thermodynamic properties of mixtures by means of a cubic EOS YOSHINORI ADACHI AND HIDEZUMI SUGIE (1998) 4 Compilation of group-contribution of VLE and excess enthalpy Kazuo Kojima (1997)

3 1 The simultaneous correlation VLE and H E Hasan Orbey, Stanley I. Sandler (1996) An attempt to represent simultaneously VLE and H E behavior Non ideal binary mixtures Representational method of VLE and H E 1 The activity coefficient model 2 Equation of state Cubic EOS and various types of mixing rules VLACM (van Laar activity coefficient) model VDWEOS (van der Waals EOS) model WSEOS (Wong Sandler ) model HVOS (Orbey Sandler ) model

4 The comparison of models Filled circles : Experimental data The dotted lines : VLACM The solid line : VDWEOS The long-dashed line : WSEOS The short-dashed line : HVOS Correlation of the excess enthalpy of 2-propanol + water system at 298.15K

5 The prediction of VLE and H E Filled circles : Experimental data The dotted lines : VLACM The solid line : VDWEOS The long-dashed line : WSEOS The short-dashed line : HVOS Prediction of the excess enthalpy of acetone + water system at 293.15K and 305.15K

6 Exception Filled circles : Experimental data The dotted lines : VLACM The solid line : VDWEOS The long-dashed line : WSEOS The short-dashed line : HVOS Prediction of the excess enthalpy of the methanol +benzene system at 298.15K with parameters obtained from the correlation of VLE data at 293.15K

7 The behavior of EOS model Filled circles :298.15K data Filled triangles : 348.15K data The dotted lines : VLACM The solid line : VDWEOS The long-dashed line : WSEOS The short-dashed line : HVOS Excess enthalpy of 2-propanol+ water system obtained from simultaneous correlation with VLE data at 298.15K, by prediction at 348.15K with the same parameter system at 293.15K and 305.15K

8 Conclusion ⑴ All of the EOS models are good or better than activity coefficient models for simultaneous correlation and prediction of VLE and H E properties. ⑵ Fitting parameters used here only to excess enthalpy data or VLE data are compatible but another property did not give satisfactory results. ⑶ The VDWEOS model led to extrapolations of reasonable accuracy for the excess enthalpy and VLE at another temperature.

9 2 Calculation of H E with the CCOR EOS The Cubic chain - of - rotators (CCOR) equation of state H.-M Lin, M.-J Lee (1994) P = RT[1 + 0.77( b / V )] + C V 0.42b V ( V bd + C)( V 0.42b) R 0.055RT ( b / V ) a V 0.42b V ( V + c) 1 The first term : simplification of the Carnahan-Starling equation for a hard-sphere fluid 2 The second term : contribution to pressure due to the rotational motion of a polyatomic molecule 3 The last two term : the attractive force contributions

10 CCOR EOS (Lin et al.,1983 ; Kim et al.,1986) 1 Applications : To calculate excess enthalpies for a variety of binary mixtures 2 Characteristics Non-polar/polar and polar/polar system The conditions in both gas and liquid phase 3 Deviations Mixtures containing long-chain molecules such as n-hexadecane Asymmetric mixture

11 Exception of CCOR EOS System : methanol + chlorobenzene [ AAD = 48 % ] ethanol + chlorobenzene [ AAD = 63.1 % ]

12 Conclusion ⑴ ⑵ ⑶ ⑷ The CCOR equation was applied to calculate excess enthalpies for a variety of mixtures in both gas and liquid states. These systems include polar/non-polar and polar/polar mixtures. However, the equation did not work as well for highly asymmetric mixtures. The substantial discrepancies between calculated and experimental data were obtained for some system.

13 3 The prediction of H E by mixing rules YOSHINORI ADACHI AND HIDEZUMI SUGIE (1998) 1 2 3 The effect of the binary interaction parameters of the modified conventional mixing rules on excess thermodynamic properties of mixture. A new method to predict thermodynamic properties of mixtures by means of a cubic EOS. The Modified Conventional Mixing Rules a m = xi x j ( 1 kij )( aia j 0.5 ) k ij = l ij + m ij ( xi x j 4 Excess molar properties ) M E + = M ' x M ' i

14 The Modified Conventional Mixing Rule circles : Experimental data : modified conventional mixing rules ( VDW with l ij =-0.253 and m ij =-0.060 ) : conventional mixing rules ( VDW with k ij =-0.269 ) Excess molar enthalpy of water- acetone system at 323.15K and 0.101325MPa

15 The effect of binary interaction parameters * Effect of l value on H E at fixed m * Effect of m value on H E at fixed l l = - 0.253 m = - 0.060 The effect of binary interaction parameter ( l, m ) on Excess enthalpy calculation

16 Conclusion ⑴ The modified conventional mixing rules give satisfactory excess enthalpy and excess volume predictions even for strongly non-ideal solution. ⑵ The calculated result for H E are more sensitive to the binary interaction parameters than those for V E.

4 Group Contribution prediction of VLE and H E Kazuo Kojima (1997) 1. VLE calculations At low temperature Activity coefficient equations At high temperature Equations of state 2. Group Contribution Methods ASOG model, UNIFAC model, The group interaction parameters can be estimated by regressing available VLE data. 3. The activity coefficient of component i COM RES ln γ i = lnγ i + lnγ i 1 Combinatorial part : the difference in size and shape of molecule 2 Residual part : intermolecular force 17

18 ASOG model lnγ COM i Si Si = 1 + ln N N S x S j j j j j x j lnγ RES i = (ln Γ k k ln Γ (i) k ) ln Γ k = ln X l m = N X a + (1 m l= 1 l l, k l= 1 N i= 1 i k, i m x x v v i= 1 i I = 1 l, i a l, k m a X l i= 1 i, k X i ) ak, l = exp( mk, l + n k, l T )

19 Group interaction parameters Available ASOG group interaction parameters [ temperature range : 298 ~ 423K ]

20 Group interaction parameters Available modified UNIFAC (Dortmund) group interaction parameters [ temperature range : 293~ 398K ]

The ASOG model for predicting the excess enthalpy 21 Available ASOG-H E group interaction parameters [ temperature range : 288~ 348K ]

22 Conclusion ⑴ Group contribution methods (ASOG, UNIFAC, modified UNIFAC) for predicting the vapor-liquid equilibrium and the excess enthalpy for mixtures of non-electrolyte solutions are presents. ⑵ Group contribution methods can be used in absence of experimental data. ⑶ The group contribution method have gained widespread popularity in chemical engineering computations.