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
📑 Key Components
🎯 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.