Vapor-Liquid Equilibria for Water+Hydrochloric Acid+Magnesium Chloride and Water+Hydrochloric Acid+Calcium Chloride Systems at Atmospheric Pressure

Similar documents
Energy, Entropy, and Availability Balances Phase Equilibria. Nonideal Thermodynamic Property Models. Selecting an Appropriate Model

General Thermodynamics for Process Simulation. Dr. Jungho Cho, Professor Department of Chemical Engineering Dong Yang University

VAPOR LIQUID EQUILIBRIUM DATA GENERATION FOR ACETIC ACID AND p-xylene AT ATMOSPHERIC PRESSURE

Supplementary Notes for Chapter 9 Mixture Thermodynamics

If two volatile and miscible liquids are combined to form a solution, Raoult s law is not obeyed. Use the experimental data in Table 9.

Journal of Polytechnic, 2017; 20 (1) : Journal of Polytechnic, 2017; 20 (1) :

Chapter 18, Part 1. Fundamentals of Atmospheric Modeling

Gasometric Determination of NaHCO 3 in a Mixture

3. Be able to derive the chemical equilibrium constants from statistical mechanics.

10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Lecture 8. Chapter 7. - Thermodynamic Web - Departure Functions - Review Equations of state (chapter 4, briefly)

Solution Thermodynamics

Application of Activity Coefficient Models in VLE Determination for Azeotropic System Using Othmer Type Ebulliometer

I wish to publish my paper on The International Journal of Thermophysics. A Practical Method to Calculate Partial Properties from Equation of State

PETE 310 Lectures # 24 & 25 Chapter 12 Gas Liquid Equilibrium

( ) 1/ 2. ( P SO2 )( P O2 ) 1/ 2.

Influence Of Operating Conditions To The Effectiveness Of Extractive Distillation Columns

Estimation of the composition of the liquid and vapor streams exiting a flash unit with a supercritical component

MODELING THE HIGH-PRESSURE BEHAVIOR OF BINARY MIXTURES OF CARBON DIOXIDE+ALKANOLS USING AN EXCESS FREE ENERGY MIXING RULE

Assessing the use of NMR chemical shifts for prediction of VLE in nonideal binary liquid mixtures

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Name: SID: Discussion Session:

Department of Physical Chemistry, Faculty of Chemistry, University of Bucharest, Bd. Regina Elisabeta 4-12, , Bucharest, Romania

4.2 Chemical Driving Force

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014

Lecture. Polymer Thermodynamics 0331 L Chemical Potential

CinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure

Phase equilibria for the oxygen-water system up to elevated temperatures and pressures

CHEMICAL ENGINEERING

Non-Commercial Use Only

New Perturbation Model for Prediction of Amino Acid and Peptide Activity Coefficients

LNG CARGO TRANSFER CALCULATION METHODS AND ROUNDING-OFFS

NAME and Section No. it is found that 0.6 mol of O

A SHORT METHOD TO CALCULATE REACTIVE RESIDUE CURVE MAPS

CALCULATION OF VAPOUR - LIQUID - LIQUID EQUILIBRIA

ARTICLE IN PRESS. Fluid Phase Equilibria 275 (2008) Contents lists available at ScienceDirect. Fluid Phase Equilibria

Appendix II Summary of Important Equations

Thermodynamics II. Department of Chemical Engineering. Prof. Kim, Jong Hak

Electrochemical Equilibrium Electromotive Force

A Self-Consistent Gibbs Excess Mixing Rule for Cubic Equations of State: derivation and fugacity coefficients

y i x P vap 10 A T SOLUTION TO HOMEWORK #7 #Problem

UNIFAC. Documentation. DDBSP Dortmund Data Bank Software Package

Adsorption: A gas or gases from a mixture of gases or a liquid (or liquids) from a mixture of liquids is bound physically to the surface of a solid.

KINETICS OF GAS HYDRATE FORMATION FROM PYROLYSIS GAS IN WATER-IN-OIL EMULSION SYSTEM

Determination of Structure and Formation Conditions of Gas Hydrate by Using TPD Method and Flash Calculations

Prediction of steady state input multiplicities for the reactive flash separation using reactioninvariant composition variables

Equation of State Modeling of Phase Equilibrium in the Low-Density Polyethylene Process

(1) The saturation vapor pressure as a function of temperature, often given by the Antoine equation:

INTRODUCTION TO CHEMICAL PROCESS SIMULATORS

Available online at Energy Procedia 4 (2011) Energy Procedia 00 (2010) GHGT-10

Determination of activity coefficients of dimethyl ether in butyl acetate and 2-propanol

Prediction of Ultrasonic Velocity in Binary Mixtures of a Nuclear Extractant and Monocarboxylic Acids using Several Theoretical Models

a for save as PDF Chemistry 163B Introduction to Multicomponent Systems and Partial Molar Quantities

salt concentration. Department of Chemical Engineering, University of Burgos, Burgos. Spain

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

Assignment 4. Adsorption Isotherms

Investigation of High-Pressure Phase Equilibrium with the Observation of the Cloud Point

The ChemSep Book. Harry A. Kooijman Consultant. Ross Taylor Clarkson University, Potsdam, New York University of Twente, Enschede, The Netherlands

Three-Phase Distillation in Packed Towers: Short-Cut Modelling and Parameter Tuning

Experimental and Modeling Studies for a Reactive Batch Distillation Column

Chem 2A Exam 1. First letter of your last name

Answers Problem Set 2 Chem 314A Williamsen Spring 2000

Design Equations. ν ij r i V R. ν ij r i. Q n components. = Q f c jf Qc j + Continuous Stirred Tank Reactor (steady-state and constant phase)

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram

Simulation and Experimental System Terner Aseton-Butanol-Ethanol with Batch Distillation

PI-21 Data correlation in homogeneous azeotropic mixtures using NRTL model and stochastic optimization methods. A. Briones-Ramírez

Thermodynamic Model of Charging the Gas/Water Interface

Electrical double layer: revisit based on boundary conditions

829. An adaptive method for inertia force identification in cantilever under moving mass

DETERMINATION OF CO 2 MINIMUM MISCIBILITY PRESSURE USING SOLUBILITY PARAMETER

Similar Constructing Method for Solving the Boundary Value Problem of the Compound Kummer Equation

Solid-Liquid Equilibrium for the System Activated Carbon and Aqueous Solution of both Phenol and Dinitro-o-Cresol.

BINARY VAPOUR-LIQUID EQUILIBRIUM FOR SYSTEMS OF INDUSTRIAL IMPORTANCE

Chapter 8 Solutions Engineering and Chemical Thermodynamics 2e

The influence of non-ideal vapor-liquid-equilibrium on vaporization of multicomponent hydrocarbon fuels

Outline. Unit Eight Calculations with Entropy. The Second Law. Second Law Notes. Uses of Entropy. Entropy is a Property.

STAT 511 FINAL EXAM NAME Spring 2001

Modeling Interfacial Tension in Liquid Liquid Systems Containing Electrolytes

Non-Ideality Through Fugacity and Activity

between standard Gibbs free energies of formation for products and reactants, ΔG! R = ν i ΔG f,i, we

PREDICTION OF GAS HYDRATE PHASE BEHAVIOR IN THE PRESENCE OF ALCOLHOLS AND GLYCOLS WITH PRSV EQUATION OF STATE AND THE VDW-P MODEL

Solutions Review Worksheet

Lecture 12. Transport in Membranes (2)

Chapter 3 Thermochemistry of Fuel Air Mixtures

Chemical Equilibrium. Chapter 6 Spontaneity of Reactive Mixtures (gases) Taking into account there are many types of work that a sysem can perform

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014

Be true to your work, your word, and your friend.

Relationship between Refractive Index and Molar Concentration of Multi-Component Solutions Zhu Xingyu 1, a, Mai Tiancheng 2, b and Zhao Zilong 2, c

Module 3: The Whole-Process Perspective for Thermochemical Hydrogen

1) Silicon oxide has a typical surface potential in an aqueous medium of ϕ,0

Ph.D. Qualifying Examination in Kinetics and Reactor Design

A New Thermodynamic Function for Phase-Splitting at Constant Temperature, Moles, and Volume

Process Modeling. Improving or understanding chemical process operation is a major objective for developing a dynamic process model

Outlet temperature of a WGS reactor (Stage I) for the conversion of CO, applied for the abatement of CO to a fixed value.

( ) Phase equilibrium Some basic principles for phase calculations

Prediction of the flash point of ternary ideal mixtures

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors

Modeling of Phase and Chemical Equilibria for Systems Involved in Biodiesel Production

Transcription:

Chnese J. Chem. Eng., 4() 76 80 (006) RESEARCH OES Vapor-Lqud Equlbra for Water+Hydrochlorc Acd+Magnesum Chlorde and Water+Hydrochlorc Acd+Calcum Chlorde Systems at Atmospherc Pressure ZHAG Yng( 张颖 ) and ZHOU Rongq( 周荣琪 )* Department of Chemcal Engneerng, snghua Unversty, Bejng 00084, Chna Abstract Vapor-lqud equlbra for water+hydrochlorc acd+magnesum chlorde and water+hydrochlorc acd+calcum chlorde systems at atmospherc pressure were measured usng a Othmer-type equlbrum stll. he epermental data are correlated usng a modfed Messner s method. Satsfactory agreements are obtaned between the epermental and the calculated results. Keywords hydrochlorc acd, magnesum chlorde, calcum chlorde, vapor-lqud equlbra IRODUCIO Aqueous soluton of a volatle strong electrolyte, such as HCl-H O system, has a mamum azeotrope so that separaton by dstllaton s dffcult. he addton of salts n ths system s an effectve method for changng the relatve volatlty and modfyng the azeotropc ponts, and the recovery of waste hydrochlorc acd by etractve dstllaton usng magnesum chlorde or calcum chlorde has receved consderable ndustral attenton []. Vapor-Lqud equlbra (VLE) data for H O+HCl +MgCl and H O+HCl+CaCl systems have been reported by Sako [] wth MgCl mole fracton of 0.0 0.058 and CaCl of 0.0 0.067. However, more concentrated salt solutons are often used to get greater salt effect n the separaton process. In ths artcle, VLE for H O+HCl+MgCl system wth MgCl mole fracton of 0.05 0.090 and H O+HCl+CaCl system wth CaCl mole fracton of 0.05 0. at atmospherc pressure were measured. he epermental results are correlated wth a modfed Messner s method [3]. EXPERIMEAL. Materals MgCl and CaCl (mnmum purty 98%) and HCl (mamum mpurty 0.% ) were purchased from Bejng Chemcal Reagents Co, all the reagents were used wthout further purfcaton. Dstllated water was used throughout the epermental work.. Apparatus and procedure he Othmer-type equlbrum stll [] was used to obtan VLE data. he pressure was controlled at (0.3±0.)kPa through the buffer tank by manual operaton and the equlbrum temperature was measured wthn ± 0.K by a standard mercury thermometer. After crculaton for about 4h, the equlbrum temperature was measured and samples of the lqud and vapor phases were wthdrawn from the equlbrum stll and the vapor phase samplng cell to determne the compostons of both phases. HCl n both phases was analyzed by acd-base ttraton wth 00mol m -3 sodum hydrode soluton. MgCl and CaCl were ttrated wth 0mol m -3 EDA (Ethylene damne tetraacetc acd, mnmum purtes 98%) soluton. he accuraces of analyss were wthn 0.000 mole fracton. o check the consstency of the epermental apparatus, the VLE for the H O+HCl system was measured under atmospherc pressure and was n good agreement wth the lterature []. 3 RESULS AD DISCUSSIO VLE data for H O ()+HCl ()+MgCl (3) and H O ()+HCl ()+CaCl (3) systems at 0.3kPa are gven n ables, where, s the lqud mole fracton of component, y s the vapor mole fracton of component, and s the equlbrum temperature. he VLE for the volatle strong electrolyte s v v y p = γ H () t ± ± Receved 005-07-06, accepted 006-0-9. * o whom correspondence should be addressed.

Vapor-Lqud Equlbra for H O+HCl+MgCl and H O+HCl+CaCl Systems at Atmospherc Pressure 77 able Epermental VLE data for H O ()+HCl ()+ MgCl (3) system at 0.3kPa 3, K y 0.0430 0.06 378.5 0.008 0.073 0.053 38. 0.043 0.0880 0.059 38. 0.059 0.0 0.05 38.6 0.938 0.05 0.033 379.8 0.0095 0.0388 0.030 38. 0.053 0.053 0.036 38.0 0.0536 0.069 0.03 38.3 0.0987 0.0778 0.03 38.0 0.0 0.030 0.046 38.6 0.009 0.053 0.0455 38.5 0.030 0.038 0.045 383.3 0.0758 0.0468 0.045 383.4 0.30 0.0566 0.0446 383.0 0.953 0.0050 0.0609 384.5 0.0088 0.0099 0.0606 385.0 0.06 0.00 0.060 385.5 0.0677 0.096 0.0607 385. 0.350 0.037 0.060 384. 0.066 0.0007 0.075 387.7 0.003 0.003 0.0739 388. 0.00 0.0079 0.079 388.0 0.0439 0.03 0.0740 387.9 0.09 0.03 0.073 386.3 0.95 0.0003 0.0885 393.4 0.003 0.000 0.0903 393.9 0.037 0.0065 0.0886 39.3 0.087 0.000 0.089 390.0 0.95 able Epermental VLE data for H O ()+HCl ()+ CaCl (3)system at 0.3kPa 3, K y 0.039 0.053 377.6 0.0063 0.0606 0.05 379.6 0.046 0.0788 0.050 38.3 0.066 0.095 0.053 38.8 0.50 0.006 0.046 38.7 0.789 0.07 0.036 380.0 0.08 0.0386 0.0356 380.9 0.030 0.05 0.0353 38.9 0.0665 0.0635 0.035 38. 0.8 0.073 0.0350 38.8 0.87 0.0080 0.0563 38.6 0.0083 0.063 0.0555 38.5 0.036 0.035 0.055 383.6 0.0775 0.049 0.0550 383.5 0.58 0.0503 0.0543 38.8 0.05 0.003 0.0747 385.6 0.00 0.0075 0.074 386.0 0.076 0.0 0.0750 386. 0.0573 0.085 0.0743 385.9 0.008 0.073 0.0739 385.0 0.789 0.00 0.093 39.0 0.0098 0.0055 0.0943 390.5 0.0636 0.00 0.0938 389. 0.3 0.04 0.0935 388.0 0.936 0.0005 0.08 397.4 0.03 0.005 0. 396.6 0.0457 0.0040 0.03 394.5 0.37 0.0078 0.090 39.3 0.0 where ± s the lqud phase mean onc mole fracton, γ ± s the mean onc actvty coeffcent, and H s Henry s law constant of the electrolyte. he VLE for the solvent s S t ap y p = () where p S s the vapor pressure of pure solvent at the system temperature, and a s the lqud phase actvty of the solvent whch s normalzed for the pure solvent at the system temperature, he γ ± n Eq.() s calculated by a modfed Messner s method. he Messner s method [4] lead to large errors n calculaton at hgh temperature and onc strength, and three-suff Margules terms are added by Sako [3] to descrbe the short-range nteractons between the chemcal speces for onc strengths hgher than the upper lmt of applcablty of the Messner s method. (Append ) he actvty of water s calculated from the Gbbs-Duhem equaton. (Append ) he Henry s law constant n Eq.() and Margules constants n Eq.(A), whch are to be determned from the epermental data, are assumed to be functons only of temperature. h0 ln H = + hln + h + h3 (3) A0 A = + Aln + A + A3 (4) B0 B = + Bln + B + B3 (5) Chnese J. Ch. E. 4() 76 (006)

78 Chnese J. Ch. E. (Vol. 4, o.) he parameters for the H O ()+HCl () system and H O ()+salt (3) system, H, A, A 3, B and B 3 were determned from the vapor pressure data of the electrolyte solutons [5] by the mnmzaton of the followng objectve functon. F = [( pep pcal) / pep] (6) and for the H O ()+salt (3) system, B 3 =0, as the system s not so complcated. he values of these parameters and mean absolute devatons n vapor pressures are reported n ables 3 6. he other parameters for the ternary systems, A 3, B 3 and Q 3, were determned from the ternary VLE data n ables by the mnmzaton of the followng objectve functon F = [( p p ) / p ] + ( y y ) (7) ep cal ep ep cal the temperature dependence of the three parameters s neglected because the temperature range of the data used s narrow. he values of these parameters and the devaton δp/p and δy are reported n ables 7 8. he results ndcate a good agreement between the able 3 Parameters of Henry s law constants System h 0 h h h 3 H O+HCl -.6088 0 9.8909 0 -.06 0 - -5.0300 0 able 4 Parameters of A for bnary system System A 0 A A A 3 H O+HCl 3.37 0 5 -.5993 0 3.996 9.566 0 3 H O+MgCl.098 0 5 9.06 0-9.0969 0-5.6698 0 3 H O+CaCl -.6865 0 5.6739 0 3.605 9.6640 0 3 able 5 Parameters of B for bnary systems System B 0 B B B 3 H O+HCl.8405 0 6.0565 0 4 -.4875 0-6.94 0 4 able 6 Mean absolute devatons for bnary systems System Range of, K Range of δp/p 0 H O+HCl 333 393 0.0 0.55 0.98 H O+MgCl 333 393 0.085 0.0795 0.84 H O+CaCl 333 393 0.076 0.43 0.70 δ / = Σ ( )/. ote: pp pep pcal pep able 7 Parameters for ternary systems System A 3 B 3 Q 3 H O+HCl+MgCl.034 0.0956 0 5.380 0 H O+HCl+CaCl -5.4709 0.75 0 3.785 0 able 8 Mean absolute devatons for ternary systems System H O+HCl+MgCl H O+HCl+CaCl Present Sako [] Present Sako [] Range of onc Strength, mol kg - 5.39 7. 6.4.50 5.00.06 5.45 3.0 δp/p 0 0.67.0 0.3 0.74 δy 0 0.78.34 0.39 0.7 δ / = Σ ( )/ ; ote: pp p p p ep cal ep δ y = Σ y y.,ep,cal Aprl, 006

Vapor-Lqud Equlbra for H O+HCl+MgCl and H O+HCl+CaCl Systems at Atmospherc Pressure 79 epermental and calculated data. Although the data reported by Sako [] could not be drectly compared wth those reported n ths work because of the dfference n epermental condtons, we were able to use the model and parameters n ths artcle to reproduce the VLE data at condtons reported by these authors. he devatons n able 8 ndcate that the results are also n good agreement wth the epermental data. When MgCl or CaCl s added n the hydrochlorc acd soluton, the concentraton of the chlorc on, the chemcal potental of HCl, and the pressure of HCl ncreases n the vapor phase. And accordng to the hydraton theory, the actvty of solvent decreases sgnfcantly due to a fracton of the water beng bound n the hydraton shells of the ons, and ths bound water of dmnshed actvty s less avalable as solvent for the ons, so actvty of the ons s rased. As a result, HCl becomes more volatle and water becomes less volatle when the salt s added. OMECLAURE A, B bnary Margules constant a actvty H Henry s law constant, kpa I onc strength, mol kg - p equlbrum pressure, kpa p t total pressure, kpa Q ternary Margules constant equlbrum temperature, K lqud phase mole fracton ± mean onc mole fracton n lqud phase y vapor phase mole fracton α constant calculated by Eq.(A4) γ onc actvty coeffcent of component on mole fracton bass γ± mean onc actvty coeffcent on mole fracton bass ν number of moles of on produced by dssocaton of one mole of electrolyte Superscrpts S saturated Subscrpts cal calculated results ep epermental results component (=,,3) ± mean onc REFERECES Furter, W.F., Salt effect n dstllaton: A lterature revew Ⅱ, the Canadan J. Chem. Eng., 55, 9 39(977) Sako,., Hakuta,., Yoshtome, H., Salt effect on vapor-lqud equlbrum for volatle strong electrolyte-water systems, J. Chem.Eng. Jpn., 7, 38 388(984) 3 Sako,., Hakuta,., Yoshtome, H., Vapor-lqud equlbra of strong electrolyte aqueous solutons, Flud Phase Equlbra, 38, 63 8(987) 4 Joseph F.Z., Handbook of aqueous electrolyte thermodynamcs: heory & applcaton, AIChE, ew York (986) 5 Sako,., Hakuta,., Yoshtome, H., Vapor pressures of bnary (H O-HCl, H O-MgCl and H O-CaCl ) and ternary aqueous (H O-MgCl -CaCl ) solutons, J. Chem. Eng. Data, 30, 4 8(985) 6 eruya, K., Hosako, S., akano,., akamor, I., Estmaton of water actvtes n multcomponent electrolyte solutons, J. Chem. Eng. Jpn., 9, 5(976) APPEDIX he actvty coeffcent of the electrolyte s calculated by combnng the Messner and Margules equatons ν ν ± ± (MEI) (MAR) ln γ = ln γ + ln γ, =, 3 (A) where γ ± ( MEI) s the mean onc actvty coeffcent of electrolyte determned by Messner s equaton [4], and γ ( MAR ) s calculated by Margules s equaton. For the bnary system ( ) ( ) ν ν 3 ± γ ± (MEI) A ± B ± ± ln γ = ln + + +.5 0.5 (A) For the ternary system ν ln γ = α ( ) + α ( ) + ± (MAR) 3 3± ± 3 ± 3± ± α ± ( ± ) α ± ± ( ± ) α3 ± 3± α3 ± 3± + [0.5( α3 + + α3 + α+ α3+ α3) Q3] ( ) α (A3) ± 3± ± ± 3± where Q 3 s the ternary Margules constant. and α j n the above equaton are related to the bnary Margules constants α = A + B, α = A + 0.5B α3 = A3 + B3, α3 = A3 + 0.5B3 α = A + B, α = A + 0.5B (A4) 3 3 3 3 3 3 ± s the lqud phase mean onc mole fracton, whch s defned by Chnese J. Ch. E. 4() 76 (006)

80 Chnese J. Ch. E. (Vol. 4, o.) /( + ν ) ± = Σ (A5) ν+ ν / ν ± + = ( ), =,3 (A6) APPEDIX he actvty of water s calculated from the Gbbs-Duhem equaton, for the bnary system + ( ν ) γ ± ln a = ln v d ln γ γ ± = ± (A7) and for the ternary system, the actvty of water s calculated by a smplfed epresson [6] ln a = I ln a (A8) k ( k ) It k where I t s the total onc strength n the ternary system, I k s the onc strength of electrolyte k n the ternary mture, and a (k) s the actvty of water n the bnary soluton contanng electrolyte k at the same temperature and onc strength as n the ternary soluton. Aprl, 006