Solution Thermodynamics, Excess Properties and Activity Models

Why this unit matters

Reaction Gibbs energy, equilibrium constants, temperature and extent. The preceding prerequisite is Other equilibria; the destination is Engineering applications. 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. Standard reaction Gibbs energy

For a defined stoichiometric reaction, equilibrium is governed by the Gibbs energy change and dimensionless reaction quotient.

Useful calculation: Δ_rG=Δ_rG°+RT ln Q; K=exp(−Δ_rG°/(RT)).

Physical interpretation: When Q

Failure to watch: Q must be dimensionless using standard-state activities; substituting raw dimensional partial pressures is a convention error.

Decision: Predict reaction direction from measured compositions.

2. Temperature effects and van’t Hoff relation

Equilibrium constants shift with temperature according to the standard reaction enthalpy, under the usual assumptions for its integration.

Useful calculation: d(ln K)/dT=Δ_rH°/(RT²); for constant Δ_rH°, ln(K2/K1)=−Δ_rH°/R(1/T2−1/T1).

Physical interpretation: An exothermic reaction generally has smaller K at higher temperature; kinetics may accelerate simultaneously, creating a tradeoff.

Failure to watch: Do not infer economic reactor performance from K alone: rate, pressure drop, catalyst deactivation and heat removal matter.

Decision: Compare operating-temperature targets with conversion limits.

3. Reaction extent and equilibrium conversion

Stoichiometric extent tracks every species consistently while a nonlinear equilibrium constraint determines feasible composition.

Useful calculation: n_i=n_i0+ν_i ξ, together with Q(ξ,T,P)=K(T).

Physical interpretation: A feed containing excess reactant may shift an equilibrium reaction but changes separation and recycle requirements downstream.

Failure to watch: The allowable extent interval comes from nonnegative mole numbers; solving outside this interval is not a valid chemical state.

Decision: Build a safeguarded scalar equilibrium solver.

4. Coupled phase and reaction equilibria

Reactive absorption and synthesis processes can involve gas–liquid distribution, reactions and ionic speciation simultaneously.

Useful calculation: At equilibrium satisfy component balances, elemental balances, μ equality and reaction Δ_rG=0 for each independent reaction.

Physical interpretation: CO₂ absorption in alkaline water is not captured by a single nonreactive gas-solubility coefficient because dissolved speciation affects the equilibrium.

Failure to watch: A unique set of independently specified reactions is required; redundant equilibrium constraints may create singular numerical systems.

Decision: Plan a staged model hierarchy from ideal screening to validated software.

Unit assessment

• Problem 1A — Standard reaction Gibbs energy. 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 standard reaction gibbs energy and a way to quantify its consequence.
• Problem 1C — Counterexample. Construct an operating case where the assumption “Q must be dimensionless using standard-state activities; substituting raw dimensional partial pressures is a convention error.” becomes important and would invalidate a simplistic answer.
• Problem 2A — Temperature effects and van’t Hoff relation. 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 temperature effects and van’t hoff relation and a way to quantify its consequence.
• Problem 2C — Counterexample. Construct an operating case where the assumption “Do not infer economic reactor performance from K alone: rate, pressure drop, catalyst deactivation and heat removal matter.” becomes important and would invalidate a simplistic answer.
• Problem 3A — Reaction extent and equilibrium conversion. 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 reaction extent and equilibrium conversion and a way to quantify its consequence.
• Problem 3C — Counterexample. Construct an operating case where the assumption “The allowable extent interval comes from nonnegative mole numbers; solving outside this interval is not a valid chemical state.” becomes important and would invalidate a simplistic answer.
• Problem 4A — Coupled phase and reaction equilibria. 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 coupled phase and reaction equilibria and a way to quantify its consequence.
• Problem 4C — Counterexample. Construct an operating case where the assumption “A unique set of independently specified reactions is required; redundant equilibrium constraints may create singular numerical systems.” becomes important and would invalidate a simplistic answer.

Assessment evidence

Unit 7 · Lesson 1: Standard reaction Gibbs energy

Learning objectives. Explain standard reaction gibbs energy 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

For a defined stoichiometric reaction, equilibrium is governed by the Gibbs energy change and dimensionless reaction quotient.

Governing relation and interpretation

Key relation: Δ_rG=Δ_rG°+RT ln Q; K=exp(−Δ_rG°/(RT)).

Engineering example and model boundary

When Q

Worked investigation

At 330 K, an illustrative standard reaction Gibbs energy is +4.0 kJ/mol. Estimate the equilibrium constant and compare it with a chosen Q.

The reaction direction test requires a dimensionless quotient and a clearly specified stoichiometric orientation.

Where the relation breaks down

Q must be dimensionless using standard-state activities; substituting raw dimensional partial pressures is a convention error.

Engineering decision and data checklist

Predict reaction direction from measured compositions.

Guided problems

1. Definition and boundary. Describe the relevant system and identify which properties must be uniform before using the governing relation for standard reaction gibbs energy.

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 Standard reaction Gibbs energy using the assumptions discussed in this lesson.

Lesson synthesis

Subject focus: The key distinguishing result of this lesson is the link between standard reaction gibbs energy and the application “Predict reaction direction from measured compositions.” The governing relation is Δ_rG=Δ_rG°+RT ln Q; K=exp(−Δ_rG°/(RT)).; the stated failure mode is Q must be dimensionless using standard-state activities; substituting raw dimensional partial pressures is a convention error.

Unit 7 · Lesson 2: Temperature effects and van’t Hoff relation

Learning objectives. Explain temperature effects and van’t hoff relation 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

Equilibrium constants shift with temperature according to the standard reaction enthalpy, under the usual assumptions for its integration.

Governing relation and interpretation

Key relation: d(ln K)/dT=Δ_rH°/(RT²); for constant Δ_rH°, ln(K2/K1)=−Δ_rH°/R(1/T2−1/T1).

Engineering example and model boundary

An exothermic reaction generally has smaller K at higher temperature; kinetics may accelerate simultaneously, creating a tradeoff.

Worked investigation

A reaction has K=4 at 300 K and ΔrH°=−40 kJ/mol. Estimate K at 320 K under constant reaction enthalpy.

For an exothermic reaction K decreases with temperature; this does not automatically imply slower reaction kinetics.

Where the relation breaks down

Do not infer economic reactor performance from K alone: rate, pressure drop, catalyst deactivation and heat removal matter.

Engineering decision and data checklist

Compare operating-temperature targets with conversion limits.

Guided problems

1. Definition and boundary. Describe the relevant system and identify which properties must be uniform before using the governing relation for temperature effects and van’t hoff relation.

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 Temperature effects and van’t Hoff relation using the assumptions discussed in this lesson.

Lesson synthesis

Subject focus: The key distinguishing result of this lesson is the link between temperature effects and van’t hoff relation and the application “Compare operating-temperature targets with conversion limits.” The governing relation is d(ln K)/dT=Δ_rH°/(RT²); for constant Δ_rH°, ln(K2/K1)=−Δ_rH°/R(1/T2−1/T1).; the stated failure mode is Do not infer economic reactor performance from K alone: rate, pressure drop, catalyst deactivation and heat removal matter.

Unit 7 · Lesson 3: Reaction extent and equilibrium conversion

Learning objectives. Explain reaction extent and equilibrium conversion 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

Stoichiometric extent tracks every species consistently while a nonlinear equilibrium constraint determines feasible composition.

Governing relation and interpretation

Key relation: n_i=n_i0+ν_i ξ, together with Q(ξ,T,P)=K(T).

Engineering example and model boundary

A feed containing excess reactant may shift an equilibrium reaction but changes separation and recycle requirements downstream.

Worked investigation

For A↔B, start with 2 mol A and 0.2 mol B. The allowable forward extent is bounded by availability of A.

Solve only within physical nonnegative species inventories; roots outside that interval describe nonexistent compositions.

Where the relation breaks down

The allowable extent interval comes from nonnegative mole numbers; solving outside this interval is not a valid chemical state.

Engineering decision and data checklist

Build a safeguarded scalar equilibrium solver.

Guided problems

1. Definition and boundary. Describe the relevant system and identify which properties must be uniform before using the governing relation for reaction extent and equilibrium conversion.

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 Reaction extent and equilibrium conversion using the assumptions discussed in this lesson.

Lesson synthesis

Subject focus: The key distinguishing result of this lesson is the link between reaction extent and equilibrium conversion and the application “Build a safeguarded scalar equilibrium solver.” The governing relation is n_i=n_i0+ν_i ξ, together with Q(ξ,T,P)=K(T).; the stated failure mode is The allowable extent interval comes from nonnegative mole numbers; solving outside this interval is not a valid chemical state.

Unit 7 · Lesson 4: Coupled phase and reaction equilibria

Learning objectives. Explain coupled phase and reaction equilibria 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

Reactive absorption and synthesis processes can involve gas–liquid distribution, reactions and ionic speciation simultaneously.

Governing relation and interpretation

Key relation: At equilibrium satisfy component balances, elemental balances, μ equality and reaction Δ_rG=0 for each independent reaction.

Engineering example and model boundary

CO₂ absorption in alkaline water is not captured by a single nonreactive gas-solubility coefficient because dissolved speciation affects the equilibrium.

Worked investigation

A gas containing CO₂ contacts an alkaline solution. Compare a simple Henry-law description with a coupled dissolution-and-acid–base equilibrium model.

Consumption of dissolved molecular CO₂ by reaction can raise total uptake; phase and reaction constraints cannot be solved separately without checking consistency.

Where the relation breaks down

A unique set of independently specified reactions is required; redundant equilibrium constraints may create singular numerical systems.

Engineering decision and data checklist

Plan a staged model hierarchy from ideal screening to validated software.

Guided problems

1. Definition and boundary. Describe the relevant system and identify which properties must be uniform before using the governing relation for coupled phase and reaction equilibria.

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 Coupled phase and reaction equilibria using the assumptions discussed in this lesson.

Lesson synthesis

Subject focus: The key distinguishing result of this lesson is the link between coupled phase and reaction equilibria and the application “Plan a staged model hierarchy from ideal screening to validated software.” The governing relation is At equilibrium satisfy component balances, elemental balances, μ equality and reaction Δ_rG=0 for each independent reaction.; the stated failure mode is A unique set of independently specified reactions is required; redundant equilibrium constraints may create singular numerical systems.

Laboratory 13: Rachford–Rice binary flash

Model: Σ zi(Ki−1)/(1+β(Ki−1))=0

Controllable variables: Feed mole fraction A (fraction) and Equilibrium ratio K_A (dimensionless).

Governing relation

Σ zi(Ki−1)/(1+β(Ki−1))=0. 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

Fixed K_B=0.4 and fixed T,P; endpoint flags identify single-phase states.

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 Feed mole fraction A to its minimum while keeping Equilibrium ratio K_A fixed. Record the change and explain which term in the governing equation is responsible.
3. Increase Feed 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 Equilibrium ratio K_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: What phase-feasibility checks should precede a bisection solve?

Extension: Change feed composition until the two-phase region disappears.

Uncertainty: How would uncertainty in equilibrium ratio k_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 Σ zi(Ki−1)/(1+β(Ki−1))=0. Demonstrate the predicted change in Feed mole fraction A at two different settings of Equilibrium ratio K_A. Conclude by explaining whether the laboratory's approximation (Fixed K_B=0.4 and fixed T,P; endpoint flags identify single-phase states.) would remain credible in a real chemical process.

Laboratory 14: Gibbs energy of ideal mixing

Model: Δg_mix=RT[x ln x+(1−x)ln(1−x)]

Controllable variables: Mole fraction A (fraction) and Temperature (K).

Governing relation

Δg_mix=RT[x ln x+(1−x)ln(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

Ideal binary solution and constant T,P, with zero reference at pure endpoints.

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 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 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 is the mixing Gibbs energy negative even without heat release?

Extension: Find the composition of largest magnitude ideal mixing free energy.

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 Δg_mix=RT[x ln x+(1−x)ln(1−x)]. Demonstrate the predicted change in Mole fraction A at two different settings of Temperature. Conclude by explaining whether the laboratory's approximation (Ideal binary solution and constant T,P, with zero reference at pure endpoints.) would remain credible in a real chemical process.

Unit 7 Assessment — Reaction equilibrium

Topic 1: Standard reaction Gibbs energy

1. State the governing physical assumption for the relation Δ_rG=Δ_rG°+RT ln Q; K=exp(−Δ_rG°/(RT))., and describe a failure case.
2. Derive or interpret the connection between standard reaction gibbs energy and the following equipment context: When Q
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: Q must be dimensionless using standard-state activities; substituting raw dimensional partial pressures is a convention error.
5. Propose an independent reference or experimental measurement that would distinguish numerical accuracy from physical accuracy.

Answer guidance: A correct explanation must state: For a defined stoichiometric reaction, equilibrium is governed by the Gibbs energy change and dimensionless reaction quotient. The central decision is: Predict reaction direction from measured compositions. Check unit conventions, the sign of the outcome and whether the selected thermodynamic state exists in the assumed phase region.

Topic 2: Temperature effects and van’t Hoff relation

1. State the governing physical assumption for the relation d(ln K)/dT=Δ_rH°/(RT²); for constant Δ_rH°, ln(K2/K1)=−Δ_rH°/R(1/T2−1/T1)., and describe a failure case.
2. Derive or interpret the connection between temperature effects and van’t hoff relation and the following equipment context: An exothermic reaction generally has smaller K at higher temperature; kinetics may accelerate simultaneously, creating a tradeoff.
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: Do not infer economic reactor performance from K alone: rate, pressure drop, catalyst deactivation and heat removal matter.
5. Propose an independent reference or experimental measurement that would distinguish numerical accuracy from physical accuracy.

Answer guidance: A correct explanation must state: Equilibrium constants shift with temperature according to the standard reaction enthalpy, under the usual assumptions for its integration. The central decision is: Compare operating-temperature targets with conversion limits. Check unit conventions, the sign of the outcome and whether the selected thermodynamic state exists in the assumed phase region.

Topic 3: Reaction extent and equilibrium conversion

1. State the governing physical assumption for the relation n_i=n_i0+ν_i ξ, together with Q(ξ,T,P)=K(T)., and describe a failure case.
2. Derive or interpret the connection between reaction extent and equilibrium conversion and the following equipment context: A feed containing excess reactant may shift an equilibrium reaction but changes separation and recycle requirements downstream.
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: The allowable extent interval comes from nonnegative mole numbers; solving outside this interval is not a valid chemical state.
5. Propose an independent reference or experimental measurement that would distinguish numerical accuracy from physical accuracy.

Answer guidance: A correct explanation must state: Stoichiometric extent tracks every species consistently while a nonlinear equilibrium constraint determines feasible composition. The central decision is: Build a safeguarded scalar equilibrium solver. Check unit conventions, the sign of the outcome and whether the selected thermodynamic state exists in the assumed phase region.

Topic 4: Coupled phase and reaction equilibria

1. State the governing physical assumption for the relation At equilibrium satisfy component balances, elemental balances, μ equality and reaction Δ_rG=0 for each independent reaction., and describe a failure case.
2. Derive or interpret the connection between coupled phase and reaction equilibria and the following equipment context: CO₂ absorption in alkaline water is not captured by a single nonreactive gas-solubility coefficient because dissolved speciation affects the equilibrium.
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: A unique set of independently specified reactions is required; redundant equilibrium constraints may create singular numerical systems.
5. Propose an independent reference or experimental measurement that would distinguish numerical accuracy from physical accuracy.

Answer guidance: A correct explanation must state: Reactive absorption and synthesis processes can involve gas–liquid distribution, reactions and ionic speciation simultaneously. The central decision is: Plan a staged model hierarchy from ideal screening to validated software. Check unit conventions, the sign of the outcome and whether the selected thermodynamic state exists in the assumed phase region.