📦 Resource excel

Thermodynamic Property Estimation Workbook (NRTL, UNIFAC, Peng-Robinson)

The Thermodynamic Property Estimation Workbook is an Excel-based computational resource that implements thermodynamic models—NRTL (Non-Random Two-Liquid), UNIFAC (Universal Quasi-Chemical Functional-Group Activity Coefficients), and Peng-Robinson Equation of State—to estimate phase equilibrium properties (e.g., activity coefficients, vapor pressures, bubble/dew points, and fugacities) for process simulation and digital twin development. It integrates empirical, semi-empirical, and cubic EOS approaches to enable rapid, spreadsheet-based property prediction without requiring commercial simulation software. Designed for engineers and researchers, it bridges fundamental thermodynamics with practical process design and operational digital twin calibration.

📖 Overview

This workbook provides a modular, transparent, and auditable platform for estimating key thermodynamic properties essential in chemical process design and operation. The NRTL model handles highly non-ideal liquid-phase behavior using binary interaction parameters and local composition theory, making it especially effective for polar and hydrogen-bonding mixtures. UNIFAC extends predictability to unseen mixtures by decomposing molecules into functional groups and applying group contribution rules—enabling estimation of activity coefficients without experimental binary data. The Peng-Robinson (PR) equation of state, implemented in its cubic form, computes vapor-phase fugacities, compressibility factors, and PVT relationships for hydrocarbon and supercritical systems, with built-in alpha-function temperature dependence and mixing rules (e.g., van der Waals one-fluid mixing) for multi-component mixtures. The workbook typically includes input sheets for component selection, mixture composition, temperature/pressure conditions, and parameter databases (e.g., UNIFAC group interaction parameters, PR critical properties, and acentric factors); calculation sheets perform iterative phase-equilibrium solving (e.g., bubble-point T/P calculations via successive substitution or Newton-Raphson); and output sheets deliver tabulated results, convergence diagnostics, and visualization-ready data. Its Excel foundation supports integration with VBA macros, sensitivity analysis tools, and live links to digital twin platforms (e.g., via Power Query or API-connected add-ins), enabling real-time property updates based on sensor-derived compositions and operating conditions.

📑 Key Components

1 NRTL Activity Coefficient Calculator
2 UNIFAC Group Contribution Engine
3 Peng-Robinson EOS Solver with Mixing Rules

🎯 Applications

  • Process design and feasibility screening for separation units (e.g., distillation, extraction)
  • Digital twin calibration and real-time thermodynamic property reconciliation
  • Teaching and conceptual validation of phase equilibrium thermodynamics

📐 Key Formulas

NRTL Activity Coefficient (γ_i)

ln γ_i = ∑_j (x_j τ_ji G_ji) / ∑_k x_k G_ki + ∑_j x_j (G_ij / ∑_k x_k G_kj) (τ_ij − ∑_m x_m τ_mj G_mj / ∑_n x_n G_nj)

Calculates the activity coefficient of component i in a liquid mixture using non-randomness parameter α_ij, interaction energy parameter τ_ij, and local composition factor G_ij.

UNIFAC Residual and Combinatorial Contributions

ln γ_i = ln γ_i^C + ln γ_i^R

Decomposes activity coefficient into combinatorial (size/shape) and residual (energy) terms based on functional group volumes, surface areas, and group interaction parameters.

Peng-Robinson EOS

P = RT/(v−b) − a(T)/(v(v+b)+b(v−b))

Cubic equation of state relating pressure (P), molar volume (v), temperature (T), and composition-dependent parameters a(T) and b; used to compute phase densities, fugacities, and saturation properties.

🔗 Related Concepts

Phase Equilibrium Activity Coefficient Models Cubic Equations of State Group Contribution Methods Digital Twin Calibration

📚 References

#thermodynamics #process-simulation #digital-twin #excel-tool #phase-equilibrium