PE Exam Exam  >  PE Exam Notes  >  Chemical Engineering for PE  >  Cheatsheet: Phase Equilibria

Cheatsheet: Phase Equilibria

1. Fundamental Concepts

1.1 Phase Rule

Parameter Definition
Gibbs Phase Rule F = C - P + 2, where F = degrees of freedom, C = number of components, P = number of phases
Degrees of Freedom Number of intensive variables that can be independently varied without changing the number of phases
Phase Homogeneous, physically distinct portion of a system separated by definite boundaries
Component Chemically independent constituents of a system

1.2 Equilibrium Criteria

Criterion Mathematical Expression
Temperature Equality Tα = Tβ = Tγ = ... for all phases
Pressure Equality Pα = Pβ = Pγ = ... for all phases (no curved interfaces)
Chemical Potential Equality μiα = μiβ = μiγ = ... for each component i in all phases
Fugacity Equality fiα = fiβ = fiγ = ... for each component i in all phases

2. Vapor-Liquid Equilibrium (VLE)

2.1 Ideal Systems - Raoult's Law

Equation Expression
Raoult's Law yiP = xiPisat, where yi = vapor mole fraction, xi = liquid mole fraction
K-value (Ideal) Ki = yi/xi = Pisat/P
Relative Volatility αij = Ki/Kj = (yi/xi)/(yj/xj) = Pisat/Pjsat
Bubble Point (T) ΣxiPisat(T) = P
Dew Point (T) Σyi/Ki = ΣyiP/Pisat(T) = 1

2.2 Non-Ideal Systems - Modified Raoult's Law

Equation Expression
Modified Raoult's Law yiφiP = xiγiPisat, where γi = activity coefficient, φi = fugacity coefficient
Low Pressure Form yiP = xiγiPisat (assumes φi ≈ 1)
Activity ai = γixi
K-value (Non-Ideal) Ki = (γiPisat)/(φiP)

2.3 Vapor Pressure Correlations

Correlation Equation
Antoine Equation log10Psat = A - B/(C + T), where T in °C or K, Psat in mmHg or kPa
Clausius-Clapeyron d(ln Psat)/dT = ΔHvap/(RT²)
Integrated C-C ln(P2sat/P1sat) = -(ΔHvap/R)(1/T2 - 1/T1)

2.4 Activity Coefficient Models

2.4.1 Margules Equations

Model Equation
Two-Suffix Margules ln γ1 = Ax2²; ln γ2 = Ax1²
Three-Suffix Margules ln γ1 = x2²[A12 + 2(A21 - A12)x1]; ln γ2 = x1²[A21 + 2(A12 - A21)x2]

2.4.2 Van Laar Equation

Component Equation
Component 1 ln γ1 = A12/[1 + (A12x1)/(A21x2)]²
Component 2 ln γ2 = A21/[1 + (A21x2)/(A12x1)]²

2.4.3 Wilson Equation

Expression Equation
Binary System ln γ1 = -ln(x1 + Λ12x2) + x212/(x1 + Λ12x2) - Λ21/(x2 + Λ21x1)]
Wilson Parameter Λij = (VjL/ViL)exp[-(λij - λii)/(RT)]

2.4.4 NRTL Equation

Parameter Definition
Binary Activity Coefficient ln γ1 = x2²[τ21(G21/(x1 + x2G21))² + (τ12G12)/(x2 + x1G12)²]
G Parameters G12 = exp(-α12τ12); G21 = exp(-α21τ21)
τ Parameters τ12 = (g12 - g22)/(RT); τ21 = (g21 - g11)/(RT)
α Parameter Non-randomness parameter, 0.2 ≤ α ≤ 0.47 for most systems

2.4.5 UNIQUAC Equation

  • Uses molecular size (r) and surface area (q) parameters
  • Two contributions: combinatorial (entropic) and residual (enthalpic)
  • ln γi = ln γiC + ln γiR
  • Combinatorial: ln γiC = ln(Φi/xi) + (z/2)qiln(θii) + li - (Φi/xi)Σxjlj
  • z = 10 (coordination number)

2.5 Azeotropes

Type Characteristics
Minimum Boiling xi = yi for all i; positive deviation from Raoult's Law; γi > 1
Maximum Boiling xi = yi for all i; negative deviation from Raoult's Law; γi <>
Azeotrope Condition γ1P1sat = γ2P2sat at azeotropic composition

3. Liquid-Liquid Equilibrium (LLE)

3.1 Fundamentals

Concept Description
Equilibrium Condition xiαγiα = xiβγiβ for each component in phases α and β
Mutual Solubility Concentrations of components in each phase at equilibrium
Distribution Coefficient KD,i = xiβ/xiα
Plait Point Point on phase diagram where two liquid phases become identical in composition

3.2 Ternary LLE Diagrams

  • Triangular diagram: three components at vertices
  • Binodal curve: boundary between one-phase and two-phase regions
  • Tie lines: connect equilibrium compositions of two phases
  • All tie lines meet at plait point
  • Lever rule applies for calculating phase amounts

3.3 Extraction Calculations

Parameter Definition
Distribution Ratio KD = concentration in extract/concentration in raffinate
Selectivity β = (KD,solute)/(KD,carrier)
Extraction Factor E = KD(S/F), where S = solvent flow, F = feed flow

4. Vapor-Liquid-Liquid Equilibrium (VLLE)

4.1 Heterogeneous Azeotropes

Property Description
Three Phases Vapor phase in equilibrium with two immiscible liquid phases
Total Pressure P = P1sat + P2sat for partially miscible binary at low mutual solubility
Vapor Composition y1 = P1sat/(P1sat + P2sat)
Degrees of Freedom F = C - P + 2 = 2 - 3 + 2 = 1 for binary VLLE

5. Solid-Liquid Equilibrium (SLE)

5.1 Solubility

Equation Expression
Ideal Solubility ln xi = -(ΔHfus/R)(1/T - 1/Tm), where Tm = melting point
Non-Ideal Solubility ln(γixi) = -(ΔHfus/R)(1/T - 1/Tm)
Eutectic Point Lowest temperature at which liquid phase exists for binary system

5.2 Crystallization

  • Supersaturation drives crystallization: S = C/Csat
  • Yield = (m0 - msat)/(m0), where m0 = initial mass, msat = mass in saturated solution
  • Material balance includes solid and liquid phases
  • Enthalpy balance includes heat of crystallization

6. High-Pressure VLE

6.1 Equations of State Approach

Method Expression
Fugacity from EOS ln(fi/P) = ln φi = ∫V[(∂P/∂ni)T,V,nj - RT/V]dV - ln Z
Equilibrium Condition φiVyiP = φiLxiP
K-value from EOS Ki = yi/xi = φiLiV

6.2 Cubic Equations of State

6.2.1 van der Waals EOS

Parameter Expression
Equation (P + a/V²)(V - b) = RT
Parameter a a = 27R²Tc²/(64Pc)
Parameter b b = RTc/(8Pc)

6.2.2 Redlich-Kwong EOS

Parameter Expression
Equation P = RT/(V - b) - a/(T0.5V(V + b))
Parameter a a = 0.42748R²Tc2.5/Pc
Parameter b b = 0.08664RTc/Pc

6.2.3 Soave-Redlich-Kwong (SRK) EOS

Parameter Expression
Equation P = RT/(V - b) - aα/(V(V + b))
Parameter a a = 0.42748R²Tc²/Pc
Parameter b b = 0.08664RTc/Pc
α Function α = [1 + m(1 - Tr0.5)]², where Tr = T/Tc
Parameter m m = 0.480 + 1.574ω - 0.176ω², where ω = acentric factor

6.2.4 Peng-Robinson (PR) EOS

Parameter Expression
Equation P = RT/(V - b) - aα/(V(V + b) + b(V - b))
Parameter a a = 0.45724R²Tc²/Pc
Parameter b b = 0.07780RTc/Pc
α Function α = [1 + κ(1 - Tr0.5)]²
Parameter κ κ = 0.37464 + 1.54226ω - 0.26992ω²

6.3 Mixing Rules

Rule Expression
van der Waals One-Fluid am = ΣΣyiyjaij; bm = Σyibi
Combining Rule aij = (aiaj)0.5(1 - kij), where kij = binary interaction parameter
Geometric Mean bij = (bi + bj)/2

6.4 Compressibility Factor

Definition Expression
Z Factor Z = PV/(RT)
Cubic Form (SRK/PR) Z³ - Z² + (A - B - B²)Z - AB = 0
Reduced Parameters A = aαP/(R²T²); B = bP/(RT)

7. Fugacity and Fugacity Coefficient

7.1 Definitions

Property Definition
Fugacity Corrected pressure for non-ideal behavior: dGi = RTd(ln fi)
Fugacity Coefficient φi = fi/(yiP) for vapor; φi = fi/(xiP) for liquid
Activity Coefficient γi = fi/(xifi0), where fi0 = standard state fugacity

7.2 Pure Component Fugacity

Phase Expression
Vapor ln(f/P) = ∫0P[(Z - 1)/P]dP
Liquid at Saturation fisat = φisatPisat
Poynting Correction fi = φisatPisatexp[ViL(P - Pisat)/(RT)]

7.3 Mixture Fugacity

  • ln φi = (1/RT)∫V[(∂P/∂ni)T,V,nj - RT/V]dV - ln Z
  • For SRK: ln φi = (bi/bm)(Z - 1) - ln(Z - B) - (A/B)[2Σyjaij/am - bi/bm]ln(1 + B/Z)
  • For PR: Similar form with modified ln term

8. Flash Calculations

8.1 Isothermal Flash

Equation Expression
Material Balance zi = xiL + yiV, where zi = feed composition, L + V = 1
Rachford-Rice Σ[zi(Ki - 1)/(1 + V(Ki - 1))] = 0
Liquid Composition xi = zi/[1 + V(Ki - 1)]
Vapor Composition yi = Kizi/[1 + V(Ki - 1)]

8.2 Solution Procedure

  1. Calculate K-values at given T and P
  2. Check: if all Ki > 1, all vapor; if all Ki < 1,="" all="">
  3. Solve Rachford-Rice equation for V (vapor fraction)
  4. Calculate xi and yi from V and Ki
  5. Update K-values if using non-ideal models and iterate

9. Henry's Law

9.1 Gas Solubility

Equation Expression
Henry's Law fi = xiHi, where Hi = Henry's constant
Low Pressure Form yiP = xiHi
Temperature Dependence d(ln H)/dT = -ΔHsol/(RT²)
Combined Approach For dilute solute: use Henry's Law; for solvent: use Raoult's Law

9.2 Applications

  • Gas absorption: CO₂, O₂, N₂ in water
  • Sparingly soluble gases in liquids
  • Valid for dilute solutions (xi → 0)
  • Hi units must be consistent with P and x

10. Excess Properties

10.1 Definitions

Property Expression
Excess Gibbs Energy GE = G - Gid = RT Σxiln γi
Excess Enthalpy HE = H - Hid
Excess Entropy SE = S - Sid
Excess Volume VE = V - Vid

10.2 Gibbs-Duhem Equation

Form Expression
General Σxid(ln γi) = 0 at constant T and P
Binary x1d(ln γ1) + x2d(ln γ2) = 0
Consistency Test 01ln(γ12)dx1 = 0 for thermodynamically consistent data

10.3 Activity Coefficient Relationships

Relationship Expression
From GE ln γi = [∂(nGE/RT)/∂ni]T,P,nj
Temperature Dependence [∂ln γi/∂T]P,x = -HiE/(RT²)
Infinite Dilution γi = limit as xi → 0

11. Critical Properties and Corresponding States

11.1 Critical Point

Property Definition
Critical Temperature Highest temperature at which liquid and vapor can coexist
Critical Pressure Vapor pressure at critical temperature
Critical Volume Molar volume at critical point
Critical Criteria (∂P/∂V)T = 0 and (∂²P/∂V²)T = 0

11.2 Reduced Properties

Property Definition
Reduced Temperature Tr = T/Tc
Reduced Pressure Pr = P/Pc
Reduced Volume Vr = V/Vc
Acentric Factor ω = -log10(Prsat)Tr=0.7 - 1.000

11.3 Generalized Correlations

Correlation Expression
Compressibility Z = Z(0) + ωZ(1), where Z(0) and Z(1) are functions of Tr and Pr
Fugacity Coefficient ln φ = ln φ(0) + ω ln φ(1)
Lee-Kesler Three-parameter corresponding states using ω

12. Multi-Component VLE

12.1 Bubble Point Calculation

Type Equation
Bubble Point T ΣKi(T)xi = 1 at specified P
Bubble Point P P = ΣxiKiP at specified T (rearranged form)
Vapor Composition yi = Kixi

12.2 Dew Point Calculation

Type Equation
Dew Point T Σyi/Ki(T) = 1 at specified P
Dew Point P Σyi/Ki = 1 at specified T
Liquid Composition xi = yi/Ki

12.3 K-value Estimation

Method Expression
Wilson Approximation Ki = (Pc,i/P)exp[5.373(1 + ωi)(1 - Tc,i/T)]
DePriester Charts Graphical correlation of K vs T and P for light hydrocarbons

13. Phase Diagrams

13.1 Binary VLE Diagrams

Type Characteristics
T-x-y (Isobaric) Temperature vs composition at constant pressure; upper curve = bubble point, lower curve = dew point
P-x-y (Isothermal) Pressure vs composition at constant temperature; upper curve = dew point, lower curve = bubble point
x-y Diagram y vs x at constant T or P; 45° line = x = y; equilibrium curve above = more volatile in vapor

13.2 Ternary Phase Diagrams

  • Triangular coordinates: each vertex = pure component (100%)
  • Distance from side = concentration of opposite component
  • Constant composition lines: parallel to opposite side
  • Sum of perpendicular distances = constant (height of triangle)

14. Thermodynamic Consistency Tests

14.1 Area Test

Test Criterion
Redlich-Kister |∫01ln(γ12)dx1| < 0.01="" for="" consistent="">
Herington Test D - J < 10,="" where="" d="">+ - A-|/(A+ + A-), J = 150|Tmax - Tmin|/Tmin

14.2 Point Test

  • Check Gibbs-Duhem at individual data points
  • Slope test: d(ln γ1)/dx1 vs d(ln γ2)/dx1 must satisfy Gibbs-Duhem
  • Infinite dilution: γ1 and γ2 must be consistent with GE model

15. Special Topics

15.1 Supercritical Extraction

Parameter Description
Supercritical Fluid T > Tc and P > Pc; properties between liquid and gas
Solvent Power Increases with density (pressure); tuneable by adjusting P and T
Common Solvents CO₂ (Tc = 31.1°C, Pc = 73.8 bar); ethane; propane

15.2 Polymer Solutions

Model Application
Flory-Huggins Activity coefficients for polymer-solvent systems
UNIFAC-FV Free-volume modification for polymers
SAFT Statistical Associating Fluid Theory for complex molecules

15.3 Electrolyte Solutions

Model Description
Debye-Hückel Activity coefficients for dilute electrolytes based on ionic strength
Pitzer Equations Virial expansion for concentrated electrolytes
NRTL-Electrolyte Extension of NRTL for systems with ions
The document Cheatsheet: Phase Equilibria is a part of the PE Exam Course Chemical Engineering for PE.
All you need of PE Exam at this link: PE Exam
Explore Courses for PE Exam exam
Get EduRev Notes directly in your Google search
Related Searches
Cheatsheet: Phase Equilibria, shortcuts and tricks, Sample Paper, Summary, pdf , Extra Questions, mock tests for examination, study material, practice quizzes, Exam, ppt, Semester Notes, Viva Questions, past year papers, Cheatsheet: Phase Equilibria, Objective type Questions, Important questions, video lectures, MCQs, Free, Previous Year Questions with Solutions, Cheatsheet: Phase Equilibria;