Vapor-Liquid Equilibria, Azeotropy and Phase Rule Applications

Why this unit matters

Liquid–liquid, dissolved gases, solids and electrolyte effects. The preceding prerequisite is Vapor–liquid systems; the destination is Reaction equilibrium. Engineering students should be able to move from a physical scenario to a transparent mathematical model and back to an interpretation of what would be observed in equipment. Thermodynamic design becomes unreliable when an input is precise to many decimal places while the phase model, standard states or material identities remain unspecified.

Four connected lessons

1. Liquid–liquid partitioning

When two poorly miscible solvents form distinct phases, equilibrium distribution depends on chemical potential equality and phase composition.

Useful calculation: A dilute partition model K=c_org/c_aq allows a first approximation to solute extraction.

Physical interpretation: A single extraction stage removes more solute when solvent volume increases, but decreasing marginal returns may favor several stages.

Failure to watch: Distribution factors depend on temperature, solvent choice and solute speciation; a dilute K cannot describe tie-lines at high loading.

Decision: Estimate extraction recovery before laboratory validation.

2. Gas solubility and Henry law

At dilute concentrations, dissolved-gas concentration can be approximately proportional to gas partial pressure at a given temperature.

Useful calculation: With the explicitly defined convention c=k_H p, solute concentration c increases linearly with gas partial pressure p.

Physical interpretation: A gas absorption contactor at higher partial pressure may dissolve more gas if equilibrium rather than mass-transfer rate is limiting.

Failure to watch: Different handbooks use reciprocal Henry coefficients, so writing down the chosen convention and units is mandatory.

Decision: Screen absorber conditions while noting reaction and salinity effects.

3. Solid–liquid and crystallization equilibrium

Cooling or evaporation can push a dissolved species beyond its saturation concentration and create a thermodynamic driving force to crystallize.

Useful calculation: Relative supersaturation σ=(c−c_sat)/c_sat when c_sat>0.

Physical interpretation: A crystallizer with c slightly above saturation may require seeding to control nucleation and product crystal size.

Failure to watch: Thermodynamic saturation predicts final equilibrium but not crystal growth kinetics, metastable zone width or fouling rates.

Decision: Distinguish equilibrium yield from throughput and morphology.

4. Osmosis and electrolyte solution limits

Solute activities control solvent chemical potential; dilute ideal solutions exhibit osmotic pressure roughly proportional to total dissolved particle concentration.

Useful calculation: π≈i c R T in an ideal dilute model; i approximates the number of dissolved particles per solute formula unit.

Physical interpretation: Reverse osmosis requires pressure exceeding the membrane-separated osmotic pressure difference, with further losses and concentration polarization.

Failure to watch: Real electrolytes require activity corrections; ideal van’t Hoff predictions lose reliability at high ionic strength.

Decision: Evaluate membrane driving-force directions with unit-consistent estimates.

Unit assessment

• Problem 1A — Liquid–liquid partitioning. Given an engineering case of your choice, identify a defensible thermodynamic boundary and relevant equation, and justify its variables.
• Problem 1B — Data reliability. Describe a measurement or property-data uncertainty that affects liquid–liquid partitioning and a way to quantify its consequence.
• Problem 1C — Counterexample. Construct an operating case where the assumption “Distribution factors depend on temperature, solvent choice and solute speciation; a dilute K cannot describe tie-lines at high loading.” becomes important and would invalidate a simplistic answer.
• Problem 2A — Gas solubility and Henry law. Given an engineering case of your choice, identify a defensible thermodynamic boundary and relevant equation, and justify its variables.
• Problem 2B — Data reliability. Describe a measurement or property-data uncertainty that affects gas solubility and henry law and a way to quantify its consequence.
• Problem 2C — Counterexample. Construct an operating case where the assumption “Different handbooks use reciprocal Henry coefficients, so writing down the chosen convention and units is mandatory.” becomes important and would invalidate a simplistic answer.
• Problem 3A — Solid–liquid and crystallization equilibrium. Given an engineering case of your choice, identify a defensible thermodynamic boundary and relevant equation, and justify its variables.
• Problem 3B — Data reliability. Describe a measurement or property-data uncertainty that affects solid–liquid and crystallization equilibrium and a way to quantify its consequence.
• Problem 3C — Counterexample. Construct an operating case where the assumption “Thermodynamic saturation predicts final equilibrium but not crystal growth kinetics, metastable zone width or fouling rates.” becomes important and would invalidate a simplistic answer.
• Problem 4A — Osmosis and electrolyte solution limits. Given an engineering case of your choice, identify a defensible thermodynamic boundary and relevant equation, and justify its variables.
• Problem 4B — Data reliability. Describe a measurement or property-data uncertainty that affects osmosis and electrolyte solution limits and a way to quantify its consequence.
• Problem 4C — Counterexample. Construct an operating case where the assumption “Real electrolytes require activity corrections; ideal van’t Hoff predictions lose reliability at high ionic strength.” becomes important and would invalidate a simplistic answer.

Assessment evidence

Unit 6 · Lesson 1: Liquid–liquid partitioning

Learning objectives. Explain liquid–liquid partitioning in physical terms; identify the relevant state or balance variables; apply a quantitative relation only within its validity range; and critically review a proposed process-engineering decision.

The physical question

When two poorly miscible solvents form distinct phases, equilibrium distribution depends on chemical potential equality and phase composition.

Governing relation and interpretation

Key relation: A dilute partition model K=c_org/c_aq allows a first approximation to solute extraction.

Engineering example and model boundary

A single extraction stage removes more solute when solvent volume increases, but decreasing marginal returns may favor several stages.

Worked investigation

An extraction removes a dye from 8 L of aqueous solution into 2 L organic solvent with partition ratio K=5 at dilute conditions.

The simple equilibrium model predicts the transferred fraction (K V_org)/(V_aq+K V_org), here 10/18, before accounting for emulsions or kinetics.

Where the relation breaks down

Distribution factors depend on temperature, solvent choice and solute speciation; a dilute K cannot describe tie-lines at high loading.

Engineering decision and data checklist

Estimate extraction recovery before laboratory validation.

Guided problems

1. Definition and boundary. Describe the relevant system and identify which properties must be uniform before using the governing relation for liquid–liquid partitioning.

2. Derivation. State which conservation principle, potential differential or equilibrium condition supports the key relation. Explain each approximation that is needed to reach the displayed form.

3. Units and sign. Determine the natural units of each term. Describe a sign error or reference-state error that would produce a misleading yet numerical-looking result.

4. Sensitivity. Change one physically adjustable input by approximately 10% and predict the direction of the response before performing an exact calculation.

5. Independent validation. Propose a laboratory measurement or alternate property relation that would help detect a model error.

6. Engineering judgment. Explain which additional information would be mandatory before the calculation could support purchase, construction or operation of real equipment.

Mini-case: two competing models

Suppose a project team has a quick spreadsheet estimate and a detailed simulation package that disagree. Model A is simpler and transparent; Model B includes correlations for effects missing from Model A. First, compare both results under the same physical conditions and reference convention. Then determine whether they are truly independent or share the same underlying property parameters. A difference between them is not automatically proof that Model B is better. Ask which result fits experimental observations, and whether the disagreement is material to the stated decision. Apply this comparison specifically to Liquid–liquid partitioning using the assumptions discussed in this lesson.

Lesson synthesis

Subject focus: The key distinguishing result of this lesson is the link between liquid–liquid partitioning and the application “Estimate extraction recovery before laboratory validation.” The governing relation is A dilute partition model K=c_org/c_aq allows a first approximation to solute extraction.; the stated failure mode is Distribution factors depend on temperature, solvent choice and solute speciation; a dilute K cannot describe tie-lines at high loading.

Unit 6 · Lesson 2: Gas solubility and Henry law

Learning objectives. Explain gas solubility and henry law in physical terms; identify the relevant state or balance variables; apply a quantitative relation only within its validity range; and critically review a proposed process-engineering decision.

The physical question

At dilute concentrations, dissolved-gas concentration can be approximately proportional to gas partial pressure at a given temperature.

Governing relation and interpretation

Key relation: With the explicitly defined convention c=k_H p, solute concentration c increases linearly with gas partial pressure p.

Engineering example and model boundary

A gas absorption contactor at higher partial pressure may dissolve more gas if equilibrium rather than mass-transfer rate is limiting.

Worked investigation

A dissolved gas exhibits c=0.7 mol/m³ at p=70 kPa. Define kH as 0.01 mol/(m³ kPa) and predict c at 90 kPa only within dilute conditions.

If a handbook instead defines p=H c, its constant is the reciprocal of this concentration-over-pressure convention.

Where the relation breaks down

Different handbooks use reciprocal Henry coefficients, so writing down the chosen convention and units is mandatory.

Engineering decision and data checklist

Screen absorber conditions while noting reaction and salinity effects.

Guided problems

1. Definition and boundary. Describe the relevant system and identify which properties must be uniform before using the governing relation for gas solubility and henry law.

2. Derivation. State which conservation principle, potential differential or equilibrium condition supports the key relation. Explain each approximation that is needed to reach the displayed form.

3. Units and sign. Determine the natural units of each term. Describe a sign error or reference-state error that would produce a misleading yet numerical-looking result.

4. Sensitivity. Change one physically adjustable input by approximately 10% and predict the direction of the response before performing an exact calculation.

5. Independent validation. Propose a laboratory measurement or alternate property relation that would help detect a model error.

6. Engineering judgment. Explain which additional information would be mandatory before the calculation could support purchase, construction or operation of real equipment.

Mini-case: two competing models

Suppose a project team has a quick spreadsheet estimate and a detailed simulation package that disagree. Model A is simpler and transparent; Model B includes correlations for effects missing from Model A. First, compare both results under the same physical conditions and reference convention. Then determine whether they are truly independent or share the same underlying property parameters. A difference between them is not automatically proof that Model B is better. Ask which result fits experimental observations, and whether the disagreement is material to the stated decision. Apply this comparison specifically to Gas solubility and Henry law using the assumptions discussed in this lesson.

Lesson synthesis

Subject focus: The key distinguishing result of this lesson is the link between gas solubility and henry law and the application “Screen absorber conditions while noting reaction and salinity effects.” The governing relation is With the explicitly defined convention c=k_H p, solute concentration c increases linearly with gas partial pressure p.; the stated failure mode is Different handbooks use reciprocal Henry coefficients, so writing down the chosen convention and units is mandatory.

Unit 6 · Lesson 3: Solid–liquid and crystallization equilibrium

Learning objectives. Explain solid–liquid and crystallization equilibrium in physical terms; identify the relevant state or balance variables; apply a quantitative relation only within its validity range; and critically review a proposed process-engineering decision.

The physical question

Cooling or evaporation can push a dissolved species beyond its saturation concentration and create a thermodynamic driving force to crystallize.

Governing relation and interpretation

Key relation: Relative supersaturation σ=(c−c_sat)/c_sat when c_sat>0.

Engineering example and model boundary

A crystallizer with c slightly above saturation may require seeding to control nucleation and product crystal size.

Worked investigation

A solution is cooled until c_sat decreases to 0.12 mol/L while actual c is 0.15 mol/L. Relative supersaturation is (0.15−0.12)/0.12=0.25.

The thermodynamic driving force alone says little about nucleation rate; crystal sizes depend on kinetics and mixing.

Where the relation breaks down

Thermodynamic saturation predicts final equilibrium but not crystal growth kinetics, metastable zone width or fouling rates.

Engineering decision and data checklist

Distinguish equilibrium yield from throughput and morphology.

Guided problems

1. Definition and boundary. Describe the relevant system and identify which properties must be uniform before using the governing relation for solid–liquid and crystallization equilibrium.

2. Derivation. State which conservation principle, potential differential or equilibrium condition supports the key relation. Explain each approximation that is needed to reach the displayed form.

3. Units and sign. Determine the natural units of each term. Describe a sign error or reference-state error that would produce a misleading yet numerical-looking result.

4. Sensitivity. Change one physically adjustable input by approximately 10% and predict the direction of the response before performing an exact calculation.

5. Independent validation. Propose a laboratory measurement or alternate property relation that would help detect a model error.

6. Engineering judgment. Explain which additional information would be mandatory before the calculation could support purchase, construction or operation of real equipment.

Mini-case: two competing models

Suppose a project team has a quick spreadsheet estimate and a detailed simulation package that disagree. Model A is simpler and transparent; Model B includes correlations for effects missing from Model A. First, compare both results under the same physical conditions and reference convention. Then determine whether they are truly independent or share the same underlying property parameters. A difference between them is not automatically proof that Model B is better. Ask which result fits experimental observations, and whether the disagreement is material to the stated decision. Apply this comparison specifically to Solid–liquid and crystallization equilibrium using the assumptions discussed in this lesson.

Lesson synthesis

Subject focus: The key distinguishing result of this lesson is the link between solid–liquid and crystallization equilibrium and the application “Distinguish equilibrium yield from throughput and morphology.” The governing relation is Relative supersaturation σ=(c−c_sat)/c_sat when c_sat>0.; the stated failure mode is Thermodynamic saturation predicts final equilibrium but not crystal growth kinetics, metastable zone width or fouling rates.

Unit 6 · Lesson 4: Osmosis and electrolyte solution limits

Learning objectives. Explain osmosis and electrolyte solution limits in physical terms; identify the relevant state or balance variables; apply a quantitative relation only within its validity range; and critically review a proposed process-engineering decision.

The physical question

Solute activities control solvent chemical potential; dilute ideal solutions exhibit osmotic pressure roughly proportional to total dissolved particle concentration.

Governing relation and interpretation

Key relation: π≈i c R T in an ideal dilute model; i approximates the number of dissolved particles per solute formula unit.

Engineering example and model boundary

Reverse osmosis requires pressure exceeding the membrane-separated osmotic pressure difference, with further losses and concentration polarization.

Worked investigation

A dilute nonelectrolyte concentration is 100 mol/m³ at 300 K. Ideal osmotic pressure is approximately 249 kPa when i=1.

Real seawater is concentrated and ionic; this dilute estimate is a demonstration of direction and scale, not a membrane-design tool.

Where the relation breaks down

Real electrolytes require activity corrections; ideal van’t Hoff predictions lose reliability at high ionic strength.

Engineering decision and data checklist

Evaluate membrane driving-force directions with unit-consistent estimates.

Guided problems

1. Definition and boundary. Describe the relevant system and identify which properties must be uniform before using the governing relation for osmosis and electrolyte solution limits.

2. Derivation. State which conservation principle, potential differential or equilibrium condition supports the key relation. Explain each approximation that is needed to reach the displayed form.

3. Units and sign. Determine the natural units of each term. Describe a sign error or reference-state error that would produce a misleading yet numerical-looking result.

4. Sensitivity. Change one physically adjustable input by approximately 10% and predict the direction of the response before performing an exact calculation.

5. Independent validation. Propose a laboratory measurement or alternate property relation that would help detect a model error.

6. Engineering judgment. Explain which additional information would be mandatory before the calculation could support purchase, construction or operation of real equipment.

Mini-case: two competing models

Suppose a project team has a quick spreadsheet estimate and a detailed simulation package that disagree. Model A is simpler and transparent; Model B includes correlations for effects missing from Model A. First, compare both results under the same physical conditions and reference convention. Then determine whether they are truly independent or share the same underlying property parameters. A difference between them is not automatically proof that Model B is better. Ask which result fits experimental observations, and whether the disagreement is material to the stated decision. Apply this comparison specifically to Osmosis and electrolyte solution limits using the assumptions discussed in this lesson.

Lesson synthesis

Subject focus: The key distinguishing result of this lesson is the link between osmosis and electrolyte solution limits and the application “Evaluate membrane driving-force directions with unit-consistent estimates.” The governing relation is π≈i c R T in an ideal dilute model; i approximates the number of dissolved particles per solute formula unit.; the stated failure mode is Real electrolytes require activity corrections; ideal van’t Hoff predictions lose reliability at high ionic strength.

Laboratory 11: Ideal-vapor dew pressure

Model: 1/Pdew=yA/PsatA+(1−yA)/PsatB

Controllable variables: Vapor mole fraction A (fraction) and Temperature (K).

Governing relation

1/Pdew=yA/PsatA+(1−yA)/PsatB. This relation is used for an educational screening calculation, not certified physical-property or process-safety design. Every output must be checked against dimensional consistency, equilibrium assumptions and an independent reference when applied to real chemicals.

Stated limitations

Low-pressure ideal vapor/liquid, hypothetical species data.

Step-by-step investigation

1. Record the default inputs and the readout with its units. Sketch the qualitative expected trend before inspecting the chart.
2. Decrease Vapor mole fraction A to its minimum while keeping Temperature fixed. Record the change and explain which term in the governing equation is responsible.
3. Increase Vapor mole fraction A to its maximum and compare the response with the first observation. Determine whether the dependence is linear, logarithmic, exponential or nonlinear for a reason you can explain.
4. Reset, then alter Temperature and repeat. Estimate a finite-difference sensitivity between two nearby parameter values, including its units.
5. Export the response curve to CSV. Explain how a plotted sensitivity helps select a range for a laboratory measurement or a process-control test.
6. Identify the first physical mechanism omitted by the model that would become important if the system were pressurized, heated, mixed or scaled up substantially.

Questions

Reasoning: Why does the dew-pressure expression use reciprocal saturation pressures?

Extension: Compare bubble and dew pressure at numerically equal xA=yA.

Uncertainty: How would uncertainty in temperature influence your interpretation of the response curve? Which single additional measurement would best distinguish model error from input uncertainty?

Evidence to submit

Export a parameter sweep from the interactive model and attach it with your explanation of 1/Pdew=yA/PsatA+(1−yA)/PsatB. Demonstrate the predicted change in Vapor mole fraction A at two different settings of Temperature. Conclude by explaining whether the laboratory's approximation (Low-pressure ideal vapor/liquid, hypothetical species data.) would remain credible in a real chemical process.

Laboratory 12: Margules nonideal bubble model

Model: P=Σ xi γi Psati; lnγ1=A(1−x)²

Controllable variables: Liquid mole fraction A (fraction) and Margules interaction A (dimensionless).

Governing relation

P=Σ xi γi Psati; lnγ1=A(1−x)². This relation is used for an educational screening calculation, not certified physical-property or process-safety design. Every output must be checked against dimensional consistency, equilibrium assumptions and an independent reference when applied to real chemicals.

Stated limitations

One-parameter two-suffix Margules, fixed T=350 K and low-pressure vapor.

Step-by-step investigation

1. Record the default inputs and the readout with its units. Sketch the qualitative expected trend before inspecting the chart.
2. Decrease Liquid mole fraction A to its minimum while keeping Margules interaction A fixed. Record the change and explain which term in the governing equation is responsible.
3. Increase Liquid mole fraction A to its maximum and compare the response with the first observation. Determine whether the dependence is linear, logarithmic, exponential or nonlinear for a reason you can explain.
4. Reset, then alter Margules interaction A and repeat. Estimate a finite-difference sensitivity between two nearby parameter values, including its units.
5. Export the response curve to CSV. Explain how a plotted sensitivity helps select a range for a laboratory measurement or a process-control test.
6. Identify the first physical mechanism omitted by the model that would become important if the system were pressurized, heated, mixed or scaled up substantially.

Questions

Reasoning: When can activity corrections create extrema in bubble pressure?

Extension: Set A=0 and verify recovery of Raoult-law results.

Uncertainty: How would uncertainty in margules interaction a influence your interpretation of the response curve? Which single additional measurement would best distinguish model error from input uncertainty?

Evidence to submit

Export a parameter sweep from the interactive model and attach it with your explanation of P=Σ xi γi Psati; lnγ1=A(1−x)². Demonstrate the predicted change in Liquid mole fraction A at two different settings of Margules interaction A. Conclude by explaining whether the laboratory's approximation (One-parameter two-suffix Margules, fixed T=350 K and low-pressure vapor.) would remain credible in a real chemical process.

Unit 6 Assessment — Other equilibria

Topic 1: Liquid–liquid partitioning

1. State the governing physical assumption for the relation A dilute partition model K=c_org/c_aq allows a first approximation to solute extraction., and describe a failure case.
2. Derive or interpret the connection between liquid–liquid partitioning and the following equipment context: A single extraction stage removes more solute when solvent volume increases, but decreasing marginal returns may favor several stages.
3. Create two credible numerical inputs using correct units, and explain the predicted response direction before computing it.
4. Assess this warning using a counterexample: Distribution factors depend on temperature, solvent choice and solute speciation; a dilute K cannot describe tie-lines at high loading.
5. Propose an independent reference or experimental measurement that would distinguish numerical accuracy from physical accuracy.

Answer guidance: A correct explanation must state: When two poorly miscible solvents form distinct phases, equilibrium distribution depends on chemical potential equality and phase composition. The central decision is: Estimate extraction recovery before laboratory validation. Check unit conventions, the sign of the outcome and whether the selected thermodynamic state exists in the assumed phase region.

Topic 2: Gas solubility and Henry law

1. State the governing physical assumption for the relation With the explicitly defined convention c=k_H p, solute concentration c increases linearly with gas partial pressure p., and describe a failure case.
2. Derive or interpret the connection between gas solubility and henry law and the following equipment context: A gas absorption contactor at higher partial pressure may dissolve more gas if equilibrium rather than mass-transfer rate is limiting.
3. Create two credible numerical inputs using correct units, and explain the predicted response direction before computing it.
4. Assess this warning using a counterexample: Different handbooks use reciprocal Henry coefficients, so writing down the chosen convention and units is mandatory.
5. Propose an independent reference or experimental measurement that would distinguish numerical accuracy from physical accuracy.

Answer guidance: A correct explanation must state: At dilute concentrations, dissolved-gas concentration can be approximately proportional to gas partial pressure at a given temperature. The central decision is: Screen absorber conditions while noting reaction and salinity effects. Check unit conventions, the sign of the outcome and whether the selected thermodynamic state exists in the assumed phase region.

Topic 3: Solid–liquid and crystallization equilibrium

1. State the governing physical assumption for the relation Relative supersaturation σ=(c−c_sat)/c_sat when c_sat>0., and describe a failure case.
2. Derive or interpret the connection between solid–liquid and crystallization equilibrium and the following equipment context: A crystallizer with c slightly above saturation may require seeding to control nucleation and product crystal size.
3. Create two credible numerical inputs using correct units, and explain the predicted response direction before computing it.
4. Assess this warning using a counterexample: Thermodynamic saturation predicts final equilibrium but not crystal growth kinetics, metastable zone width or fouling rates.
5. Propose an independent reference or experimental measurement that would distinguish numerical accuracy from physical accuracy.

Answer guidance: A correct explanation must state: Cooling or evaporation can push a dissolved species beyond its saturation concentration and create a thermodynamic driving force to crystallize. The central decision is: Distinguish equilibrium yield from throughput and morphology. Check unit conventions, the sign of the outcome and whether the selected thermodynamic state exists in the assumed phase region.

Topic 4: Osmosis and electrolyte solution limits

1. State the governing physical assumption for the relation π≈i c R T in an ideal dilute model; i approximates the number of dissolved particles per solute formula unit., and describe a failure case.
2. Derive or interpret the connection between osmosis and electrolyte solution limits and the following equipment context: Reverse osmosis requires pressure exceeding the membrane-separated osmotic pressure difference, with further losses and concentration polarization.
3. Create two credible numerical inputs using correct units, and explain the predicted response direction before computing it.
4. Assess this warning using a counterexample: Real electrolytes require activity corrections; ideal van’t Hoff predictions lose reliability at high ionic strength.
5. Propose an independent reference or experimental measurement that would distinguish numerical accuracy from physical accuracy.

Answer guidance: A correct explanation must state: Solute activities control solvent chemical potential; dilute ideal solutions exhibit osmotic pressure roughly proportional to total dissolved particle concentration. The central decision is: Evaluate membrane driving-force directions with unit-consistent estimates. Check unit conventions, the sign of the outcome and whether the selected thermodynamic state exists in the assumed phase region.