Supplementary Notes for Chapter 9 Mixture Thermodynamics

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

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

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.

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

Lecture. Polymer Thermodynamics 0331 L Chemical Potential

Solution Thermodynamics

Solution Thermodynamics

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

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

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

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

NUMERICAL DIFFERENTIATION

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

Name: SID: Discussion Session:

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

CHEMICAL REACTIONS AND DIFFUSION

Non-Ideality Through Fugacity and Activity

Appendix II Summary of Important Equations

Thermodynamics General

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

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

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

4.2 Chemical Driving Force

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

McCabe-Thiele Diagrams for Binary Distillation

...Thermodynamics. If Clausius Clapeyron fails. l T (v 2 v 1 ) = 0/0 Second order phase transition ( S, v = 0)

MA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials

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

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

Lecture 12: Discrete Laplacian

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

NAME and Section No.

Module 3: Element Properties Lecture 1: Natural Coordinates

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

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

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

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Osmotic pressure and protein binding

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

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

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

Inductance Calculation for Conductors of Arbitrary Shape

Chapter 18, Part 1. Fundamentals of Atmospheric Modeling

Electrical double layer: revisit based on boundary conditions

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

Numerical Solution of Ordinary Differential Equations

CHEMICAL ENGINEERING

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

Numerical Heat and Mass Transfer

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

V T for n & P = constant

Thermodynamics. Section 4

The Geometry of Logit and Probit

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

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

This column is a continuation of our previous column

THE SUMMATION NOTATION Ʃ

Lecture 7: Boltzmann distribution & Thermodynamics of mixing

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

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

Gasometric Determination of NaHCO 3 in a Mixture

Phase equilibria Introduction General equilibrium conditions

Chapter 13: Multiple Regression

Marginal Effects in Probit Models: Interpretation and Testing. 1. Interpreting Probit Coefficients

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Part II Thermodynamic properties Properties of Pure Materials Chapter 8

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

Assignment 4. Adsorption Isotherms

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

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

Computation of Phase Equilibrium and Phase Envelopes

Difference Equations

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

Lecture 3. Ax x i a i. i i

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

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

INTRODUCTION TO CHEMICAL PROCESS SIMULATORS

( ) Phase equilibrium Some basic principles for phase calculations

More metrics on cartesian products

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

Linear Regression Analysis: Terminology and Notation

x = , so that calculated

#64. ΔS for Isothermal Mixing of Ideal Gases

Laboratory 3: Method of Least Squares

Formal solvers of the RT equation

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

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

Lecture 10 Support Vector Machines II

Review of Classical Thermodynamics

Implicit Integration Henyey Method

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

Module 3: The Whole-Process Perspective for Thermochemical Hydrogen

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.

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

Week 9 Chapter 10 Section 1-5

Singular Value Decomposition: Theory and Applications

Transcription:

Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects of mxng on heat and work nteractons and state property changes 4. Gbbs-Duhem relatonshp and thermodynamc consstency 5. Mxng functons 6. Ideal gas mxtures and deal solutons 7. Fugacty and actvty concepts 8. Fugacty coeffcents from PVTN EOS property relatonshps 9. Actvty coeffcents from G property relatonshps 10. Ideal reversble work effects n mxng or separatng components 1. Notaton and operatonal equatons for mxtures Partal molar propertes are extensve propertes that determne how derved propertes change as a functon of mole number or composton. For a general property, B, whch could be U, H, S, V, or A, the partal molar B s defned as: B B = (9-8) N TPN,, j [] There are a few mportant relatonshps n Chapter 9 that allow you to calculate partal molar propertes gven the extensve form B of the property or the ntensve form B. An mportant one of these s: B B = B x j (9-53) j 1 x j TPxj,,, Note that for a bnary, two component system that Eq.(9-53) quckly yelds the tangentntercept method of evaluatng partal molar propertes from a graph of B versus x j at constant T and P. For example, see Fgure 9.1 and the dscusson n the text. 2. PVTN EOSs for mxtures Volume and pressure explct equatons of state are commonly used to represent the volumetrc propertes of flud mxtures of both gases and lquds. The key feature that dstngushes a mxture EOS from ts pure component counterpart s the presence of compostonal dependence. Ths dependence expresses tself n the form of so-called mxng rules that ncorporate pure component EOS parameters and weght them [ ] Modfed: 10/16/2003 1

proportonally to the concentraton of each component followng a specfc mathematcal recpe. Wth the excepton of the vral equaton of state no rgorous theoretcal approach exsts to specfy a mxng rule recpe. The most common approach s to use some varaton of the Lorentz-Bertlelot rules used by van der Waals and others over a century ago (see Eqs (9-17 and 9-18). In mxtures that exhbt consderable non-dealty, a bnary or hgher order nteracton parameter δ j s ntroduced to capture specfc nteractons between molecules of type and j. Each mxture PVTN EOS wll have some prescrbed recpe for ts mxng rules. For example see the conventons followed for the RK and PR EOS on pages 324-326. In stuatons where hgh level non-deal effects are present often more complex forms for mxng rules are ntroduced. The Wong-Sandler and Chueh-Prausntz rules dscussed on pages 327-328 are examples of ths type of mxng rule. 3. General effects of mxng on heat and work nteractons and state property changes There are a few typcal classes of problems that you should be able to solve. These nclude: Calculatng the change n enthalpy as a result of mxng to determne a heat nteracton for mantanng a constant temperature or to follow some prescrbed recpe for temperature, for example, see Problem 9.2. Calculatng a change n volume as a result of mxng to determne a change n densty or the magntude of a work nteracton at constant pressure. 4. Gbbs-Duhem relatonshp and thermodynamc consstency For mxture data and correlatons for any property B, the partal molar quanttes are nterrelated through the use of the general Gbbs-Duhem relatonshp: n B B xdb = dt dp = 1 T + Px, P Tx, Note that B s a functon of T, P, and x for = 1,, n. The Gbbs-Duhem relatonshp can be appled to any partal molar property, such as: H,V, µ, ln ˆφ, ln ˆf,G, ln a, ln γ,etc. For a bnary system, the Gbbs-Duhem relatonshp s frequently used to check thermodynamc consstency of thermodynamc data; for example, actvty coeffcents. Addtonally, for a bnary mxture, f you have a measurement of the partal molar property of one component as a functon of composton you can determne the property for the other component. Modfed: 10/16/2003 2

5. Mxng functons n + + + + ( 1 ) ( 1 ) B B T,P,x, =,...,n x B T,P,x, =,...,n (9-68) mx = 1 whch s n ntensve form. A smlar expresson results for B mx n extensve form wth x replaced by N. You should be famlar wth the concept of a reference state (+). Strctly speakng, reference states are arbtrary, but n practce several common forms appear for x +, T +, and P +, for example: (1) pure at T and P of the mxture (2) x + 0 at T and P of the mxture that s n an nfnte dluton state where Henry s law behavor s followed (3) a fxed composton of 1 molal or 1 molar concentraton at T and P of the mxture that behaves n some deal manner commonly used for electrolytes (see Chapter 12) 6. Ideal gas mxtures and deal solutons In ths secton of the text, a set of defntons were used to characterze non-deal solutons n terms of a devaton from deal behavor. g deal gas mxture: G =µ = RT ln yp +λ( T ) ID deal soluton: G =µ = RT ln x +Λ ( T,P ) where λ and Λ are constants specfc to component. The key ponts to remember are that partal pressure y P s the deal gas mxture compostonally dependent varable whle for deal condensed phase soluton t s the mole fracton x. 7. Fugacty and actvty concepts For a real mxture or soluton we ntroduce the followng models: real flud mxture: G = µ = RT ln ˆf +λ + + + + real soluton: G =µ = RT ln a +µ ( T,P,x ) whch lead to the followng defntons of the fugacty coeffcent and actvty coeffcent ˆφ and γ, respectvely: and ˆ ˆf φ ---- represents a devaton from deal gas/flud mxture behavor yp ˆf f ˆ γ = ---- represents a devaton from deal soluton behavor a ˆ f + x Modfed: 10/16/2003 3

8. Fugacty coeffcent relatonshps from PVTN EOS property models There are two basc approaches one nvolvng pressure explct EOSs lke the PR or RKS, and the other for volume explct EOSs, such as Vral or Correspondng States formulatons, ncludng compressblty expansons n densty or molar volume or smlarly structured equatons. Equatons (9-142) and (9-143) provde convenent forms for pressure explct EOS models for mxture and pure components, respectvely, whle Eqs(9-129) and (9-127) work for volume explct EOS models. For example, for component n a mxture: or RT V ln ˆ P RT φ = dv RT ln Z N V T,V,Nj [] RT ln ˆφ =+ P 0 RT V P dp At ths pont you should know how to calculate the fugacty or fugacty coeffcent for a pure component usng a pressure- or volume-explct EOS or for a component n a bnary mxture usng a sutable PVTN EOS that has been properly formulated wth mxng rules for ts constants (for example, the a mx and b mx constants n the RKS EOS) n terms of pure component propertes and a bnary nteracton parameter, eg δ j. Beng able to do ths provdes a powerful tool for calculatons requred later n the course, for example to estmate the vapor pressure of a pure component you would equate ˆφ for the lqud and vapor phases at a gven temperature by estmatng the P and usng the EOS to calculate lqud and vapor volumes (denstes) untl the fugacty coeffcents for each phase were equal. 9. Actvty coeffcents from G property relatonshps The actvty coeffcent s defned n terms of the partal molar excess Gbbs free energy of mxng and can be drectly related to the fugacty usng the defnton of actvty: G RT ln γ = G = N [] T,P,Nj + where γ a ˆ ˆ x = f f x To evaluate the partal dervatve we need an expresson for G : ID G Gmx G = f T, P, N ( = 1,..., n) whch represents the dfference between the actual enthalpy of mxng and the enthalpy of mxng for an deal soluton at the same T, P and composton as the actual mxture. Modfed: 10/16/2003 4

Typcally, one has access to a model that gves G = G / N as a functon of T, P, and x, and wth Eq. (9-53) one can easly obtan γ for each component, e.g. G RT ln γ = G x j x j T,P,x j, In addton, the Gbbs-Duhem equaton can be appled to calculate the other actvty coeffcent gven a set of data for one component (e.g. f γ solvent s known as a functon of composton then γ solute can be estmated by ntegratng the Gbbs-Duhem relatonshp). Conversely, f both actvty coeffcents of a bnary mxture are known then the Gbbs- Duhem equaton can be used to check the thermodynamc consstency of a gven set of data. For example, the slope and area tests have been developed specfcally for ths purpose (see pages 357-358 and Fg. 9.3). A key ssue here s how to deal wth the standard or reference state condton (+), as that wll have a drect effect on the magntude and behavor of γ. If a symmetrc reference state s used then ˆ f + = f( pure ) and the Lews and Randall rule s followed as x goes to 1.0 wth γ approachng 1.0. Alternatvely, a unsymmetrcal reference state can be used where the nfnte dluton behavor as defned by Henry s Law determnes that γ ** approaches 1.0 as x goes to 0. Ths s frequently called the McMllan-Mayer reference condton. Another popular alternatve commonly employed for electrolytes s to use a 1 molal soluton at T and P of the mxture where the mxture follows Henry s Law n dlute solutons as x or m 0 (see Fg. 9.9 for example). The dscusson n the text from pages 360 to 365 should be carefully revewed to see how the actvty coeffcent s related to fugacty behavor for dfferent reference state condtons. Partcular attenton should be pad to Fgures 9.6 9.9. 10. Ideal reversble work assocated wth mxng or separatng components Under sothermal condtons the reversble work s equal to the net change n Gbbs free energy assocated wth the process, more generally the reversble work can be related drectly to the change n avalablty usng the methods ntroduced n Chapter 14. Examples 9.8 and 9.9 should be revewed to see how the concepts are employed. The approach s qute straghtforward f you are dealng wth a steady state process. Thngs are a bt more complcated f the system or process condtons are changng wth tme, but the general concept remans unchanged. [ ] Modfed: 10/16/2003 5