15.3 Estimation Methods and Corresponding States Principles

15.3 Estimation Methods and Corresponding States Principles

In chemical process design and simulation, engineers frequently encounter novel reaction intermediates, specialty solvents, or complex multicomponent hydrocarbon blends where experimental physical property data are incomplete or unavailable. In such cases, the NCEES PE Chemical Exam tests the ability to predict critical constants ($T_c, P_c, V_c$), compressibility factors ($Z$), vapor pressures, and phase equilibrium activity coefficients using Corresponding States Principles (CSP) and Group Contribution Methods.

---

1. The Principle of Corresponding States (CSP)

Two-Parameter Corresponding States (van der Waals)

Originally deduced by van der Waals in 1873, the two-parameter Principle of Corresponding States asserts that all fluids, when compared at the same reduced temperature ($T_r$) and reduced pressure ($P_r$), possess the same compressibility factor ($Z$) and identical dimensionless thermodynamic properties:

$$T_r = \frac{T}{T_c}, \quad P_r = \frac{P}{P_c}, \quad Z = f(T_r, P_r)$$

While two-parameter CSP provides excellent accuracy for simple, spherical, non-polar fluids (argon, krypton, xenon, and to a good approximation methane), it breaks down significantly for non-spherical, elongated, polar, or hydrogen-bonding molecules.

Three-Parameter Corresponding States: Pitzer's Acentric Factor ($\omega$)

In 1955, Kenneth Pitzer introduced a third parameter—the acentric factor ($\omega$)—to quantify the degree of molecular non-sphericity and acentricity of intermolecular force fields. Pitzer defined $\omega$ based on the reduced saturated vapor pressure at a reduced temperature of $T_r = 0.700$:

$$\omega \equiv -\log_{10}\left( P_r^{\text{sat}} \right)_{T_r = 0.700} - 1.000$$

Physical Basis of $\omega$:

• For simple spherical molecules (argon, krypton, xenon), experimental measurements confirm that $(P_r^{\text{sat}})_{T_r=0.700} \approx 0.100$. Taking $-\log_{10}(0.100) - 1.000 = -(-1.000) - 1.000 = 0.000$. Thus, $\omega \equiv 0.000$ for simple spherical fluids.
• For non-spherical or elongated hydrocarbons, intermolecular attractive forces are stronger and vapor pressures at $T_r = 0.700$ are lower than $0.100$ ($(P_r^{\text{sat}})_{T_r=0.7} < 0.100$). Consequently, $-\log_{10}(P_r^{\text{sat}}) > 1.000$, yielding positive acentric factors ($\omega > 0$).

Representative values:

• Argon ($\text{Ar}$): $\omega = 0.000$
• Methane ($\text{CH}_4$): $\omega = 0.011$
• Ethane ($\text{C}_2\text{H}_6$): $\omega = 0.099$
• Propane ($\text{C}_3\text{H}_8$): $\omega = 0.152$
• n-Butane ($\text{C}_4\text{H}_{10}$): $\omega = 0.200$
• n-Octane ($\text{C}_8\text{H}_{18}$): $\omega = 0.398$
• Benzene ($\text{C}_6\text{H}_6$): $\omega = 0.210$
• Water ($\text{H}_2\text{O}$): $\omega = 0.344$

Pitzer Generalized Compressibility Formulation

Under three-parameter CSP, the compressibility factor $Z$ is expanded as a linear function of $\omega$:

$$Z = Z^{(0)}(T_r, P_r) + \omega Z^{(1)}(T_r, P_r)$$

Where:

• $Z^{(0)}$ = compressibility factor of a simple spherical fluid ($\omega = 0$).
• $Z^{(1)}$ = reference correction function (tabulated alongside $Z^{(0)}$ in Lee-Kesler tables in the NCEES handbook).

Generalized Virial Equation (Pitzer-Curl-Tsonopoulos)

At low to moderate pressures ($P_r < 0.5$), the virial equation truncated at the second coefficient provides rapid, closed-form compressibility calculations:

$$Z = 1 + \frac{B P}{R T} = 1 + \left( \frac{B P_c}{R T_c} \right) \frac{P_r}{T_r}$$ $$\frac{B P_c}{R T_c} = B^{(0)}(T_r) + \omega B^{(1)}(T_r)$$

Where the generalized functions are:

$$B^{(0)}(T_r) = 0.1445 - \frac{0.330}{T_r} - \frac{0.1385}{T_r^2} - \frac{0.0121}{T_r^3} - \frac{0.000607}{T_r^8}$$ $$B^{(1)}(T_r) = 0.0637 + \frac{0.331}{T_r^2} - \frac{0.423}{T_r^3} - \frac{0.008}{T_r^8}$$

---

2. Multicomponent Gas Mixtures: Pseudo-Critical Properties

To apply generalized corresponding states correlations to gas mixtures, the blend is treated as a single pseudo-pure fluid with effective pseudo-critical properties.

Kay's Rule

The simplest and most widely tested mixing rule on the PE exam is Kay's Rule, which calculates mole-fraction-weighted linear averages of pure-component critical constants:

$$T_{pc} = \sum_{i=1}^{C} y_i T_{c,i}$$ $$P_{pc} = \sum_{i=1}^{C} y_i P_{c,i}$$ $$\omega_m = \sum_{i=1}^{C} y_i \omega_i$$

The pseudo-reduced temperature ($T_{pr}$) and pseudo-reduced pressure ($P_{pr}$) are defined as:

$$T_{pr} = \frac{T}{T_{pc}}, \quad P_{pr} = \frac{P}{P_{pc}}$$

Kay's rule provides excellent accuracy when components have comparable critical properties ($T_{c,i}$ and $P_{c,i}$ within a factor of $2$).

Prausnitz-Gunn Modification

For mixtures containing molecules with widely disparate sizes (e.g., hydrogen/hydrocarbon blends), the Prausnitz-Gunn rule improves pressure predictions by weighting molar volume:

$$T_{pc} = \sum y_i T_{c,i}, \quad V_{pc} = \sum y_i V_{c,i}, \quad Z_{pc} = \sum y_i Z_{c,i}$$ $$P_{pc} = \frac{Z_{pc} R T_{pc}}{V_{pc}}$$

---

3. Group Contribution Methods: The Joback Method

Group contribution methods decompose a molecule into distinct structural building blocks (e.g., $-\text{CH}_3, -\text{CH}_2-, -\text{OH}, -\text{COOH}$). Each functional group is assigned a fixed empirical contribution determined from extensive experimental regression. Summing these group values allows prediction of properties for unmeasured compounds.

The Joback Method (Joback and Reid, 1987)

The Joback method is the standard group contribution framework featured in the NCEES PE Chemical Reference Handbook for estimating critical properties and thermochemical parameters.

   2-Butanone (Methyl Ethyl Ketone): CH3 - C(=O) - CH2 - CH3
   
   +-------------+    +-------------+    +-------------+    +-------------+
   | -CH3 group  |    | -C(=O)- grp |    | -CH2- group |    | -CH3 group  |
   | (Terminal)  |----+  (Ketone)   |----+ (Methylene) |----+ (Terminal)  |
   +-------------+    +-------------+    +-------------+    +-------------+
   Total Contribution = Sum( N_k * Delta_Property_k )

1. Normal Boiling Point ($T_b$)

$$T_b = 198.2 + \sum_k N_k \Delta T_{b,k} \quad (\text{K})$$

2. Critical Temperature ($T_c$)

$$T_c = T_b \left[ 0.584 + 0.965 \sum_k N_k \Delta T_{c,k} - \left( \sum_k N_k \Delta T_{c,k} \right)^2 \right]^{-1} \quad (\text{K})$$

3. Critical Pressure ($P_c$)

$$P_c = \left[ 0.113 + 0.0032 N_A - \sum_k N_k \Delta P_{c,k} \right]^{-2} \quad (\text{bar})$$

Where $N_A$ is the total number of atoms in the molecule.

4. Critical Volume ($V_c$)

$$V_c = 17.5 + \sum_k N_k \Delta V_{c,k} \quad (\text{cm}^3/\text{mol})$$

5. Standard Enthalpy of Formation ($\Delta H_f^\circ$ at $298.15\text{ K}$, Ideal Gas)

$$\Delta H_f^\circ(298\text{ K}) = 68.29 + \sum_k N_k \Delta H_{f,k} \quad (\text{kJ/mol})$$

Representative Joback Group Contributions

---

4. UNIFAC Group Contribution for Liquid Phase Equilibria

While Joback predicts pure-component critical constants, UNIFAC (Universal Quasi-Chemical Functional Group Activity Coefficients) predicts liquid-phase activity coefficients ($\gamma_i$) in non-ideal liquid mixtures directly from molecular fragments without requiring binary experimental VLE/LLE data.

In the UNIFAC framework, the natural logarithm of the activity coefficient is decomposed into two distinct contributions:

$$\ln \gamma_i = \ln \gamma_i^C + \ln \gamma_i^R$$

• Combinatorial Part ($\ln \gamma_i^C$): Governed by entropic and steric differences in molecular size and surface area. Computed directly from pure-component group volume parameters ($R_k$) and surface area parameters ($Q_k$): $$\ln \gamma_i^C = \ln\left( \frac{\Phi_i}{x_i} \right) + 1 - \frac{\Phi_i}{x_i} - \frac{z}{2} q_i \left[ \ln\left( \frac{\Phi_i}{\theta_i} \right) + 1 - \frac{\Phi_i}{\theta_i} \right]$$ Where $\Phi_i$ is segment fraction, $\theta_i$ is area fraction, and $z = 10$ is the lattice coordination number.
• Residual Part ($\ln \gamma_i^R$): Governed by energetic intermolecular pairwise interactions between functional groups: $$\ln \gamma_i^R = \sum_k \nu_k^{(i)} \left[ \ln \Gamma_k - \ln \Gamma_k^{(i)} \right]$$ Where $\Gamma_k$ is the group activity coefficient in the mixture and $\Gamma_k^{(i)}$ is the group activity coefficient in a reference solution of pure component $i$, evaluated using temperature-dependent binary group interaction parameters ($a_{mn}$ and $a_{nm}$).

---

5. Summary Comparison Table: Estimation Methodologies

---

6. Comprehensive Worked Numerical Example

Problem Statement

An engineer is characterizing 2-butanone (methyl ethyl ketone, $\text{CH}_3-\text{CO}-\text{CH}_2-\text{CH}_3$, chemical formula $\text{C}_4\text{H}_8\text{O}$) and designing a natural gas separator. Complete the following three engineering evaluations:

1. Joback Critical Constant Estimation: Using the Joback method, estimate the normal boiling point $T_b$ in $\text{K}$, critical temperature $T_c$ in $\text{K}$, critical pressure $P_c$ in $\text{bar}$, critical volume $V_c$ in $\text{cm}^3/\text{mol}$, and standard heat of formation $\Delta H_f^\circ(298\text{ K})$ in $\text{kJ/mol}$ for 2-butanone.
2. Acentric Factor Calculation: Experimental measurements on a newly synthesized fluorinated hydrocarbon indicate $T_c = 360.0\text{ K}$ and $P_c = 40.0\text{ bar}$. At $T = 252.0\text{ K}$, the saturated vapor pressure is measured as $P^{\text{sat}} = 2.45\text{ bar}$. Determine the Pitzer acentric factor $\omega$ of this fluid.
3. Kay's Rule Gas Blend Processing: A natural gas blend contains $70.0\text{ mol}\%$ methane ($T_{c1} = 190.6\text{ K}, P_{c1} = 46.0\text{ bar}$) and $30.0\text{ mol}\%$ propane ($T_{c2} = 369.8\text{ K}, P_{c2} = 42.5\text{ bar}$). The separator operates at $T = 320.0\text{ K}$ and $P = 44.0\text{ bar}$. Calculate the mixture pseudo-critical temperature $T_{pc}$, pseudo-critical pressure $P_{pc}$, and pseudo-reduced properties ($T_{pr}, P_{pr}$).

---

Step 1: Joback Property Estimation for 2-Butanone

Decompose 2-butanone ($\text{CH}_3-\text{CO}-\text{CH}_2-\text{CH}_3$):

• $2 \times (-\text{CH}_3)$ group
• $1 \times (-\text{CH}_2-)$ group
• $1 \times (-\text{C}(=\text{O})-)$ (ketone) group
• Total atoms: $N_A = 4\text{ (C)} + 8\text{ (H)} + 1\text{ (O)} = 13\text{ atoms}$.

Sum the group contributions:

• Boiling Point ($\sum \Delta T_b$): $$\sum \Delta T_b = 2(23.58) + 1(22.88) + 1(76.75) = 47.16 + 22.88 + 76.75 = 146.79\text{ K}$$ $$T_b = 198.2 + 146.79 = \mathbf{344.99\text{ K}} \quad (71.8^\circ\text{C}, \text{experimental: } 352.8\text{ K})$$
• Critical Temperature ($\sum \Delta T_c$): $$\sum \Delta T_c = 2(0.0141) + 1(0.0189) + 1(0.0380) = 0.0282 + 0.0189 + 0.0380 = 0.0851$$ $$T_c = 344.99 \times \left[ 0.584 + 0.965(0.0851) - (0.0851)^2 \right]^{-1}$$ $$\text{Bracket} = 0.584 + 0.08212 - 0.00724 = 0.65888$$ $$T_c = \frac{344.99}{0.65888} = \mathbf{523.6\text{ K}} \quad (250.4^\circ\text{C}, \text{experimental: } 535.5\text{ K})$$
• Critical Pressure ($\sum \Delta P_c$): $$\sum \Delta P_c = 2(-0.0012) + 1(0.0000) + 1(0.0031) = -0.0024 + 0.0031 = 0.0007$$ $$P_c = \left[ 0.113 + 0.0032(13) - 0.0007 \right]^{-2} = \left[ 0.113 + 0.0416 - 0.0007 \right]^{-2} = [0.1539]^{-2} = \mathbf{42.22\text{ bar}} \quad (\text{experimental: } 42.1\text{ bar})$$
• Critical Volume ($\sum \Delta V_c$): $$\sum \Delta V_c = 2(65) + 1(56) + 1(61) = 130 + 56 + 61 = 247\text{ cm}^3/\text{mol}$$ $$V_c = 17.5 + 247 = \mathbf{264.5\text{ cm}^3/\text{mol}} \quad (\text{experimental: } 267\text{ cm}^3/\text{mol})$$
• Standard Heat of Formation ($\sum \Delta H_f$): $$\sum \Delta H_f = 2(-76.45) + 1(-20.64) + 1(-133.22) = -152.90 - 20.64 - 133.22 = -306.76\text{ kJ/mol}$$ $$\Delta H_f^\circ(298\text{ K}) = 68.29 + (-306.76) = \mathbf{-238.47\text{ kJ/mol}} \quad (\text{experimental: } -238.5\text{ kJ/mol})$$

(Notice the remarkable predictive accuracy: calculated $P_c = 42.22\text{ bar}$ vs experimental $42.1\text{ bar}$, and calculated $\Delta H_f^\circ = -238.5\text{ kJ/mol}$ matches experimental within $0.1\%$).

---

Step 2: Pitzer Acentric Factor Calculation

Check reduced temperature:

$$T_r = \frac{T}{T_c} = \frac{252.0\text{ K}}{360.0\text{ K}} = 0.700$$

Because $T_r$ is exactly $0.700$, Pitzer's definition applies directly:

$$P_r^{\text{sat}} = \frac{P^{\text{sat}}}{P_c} = \frac{2.45\text{ bar}}{40.0\text{ bar}} = 0.06125$$ $$\log_{10}(P_r^{\text{sat}}) = \log_{10}(0.06125) = -1.21289$$ $$\omega = -\log_{10}(P_r^{\text{sat}}) - 1.000 = -(-1.21289) - 1.000 = 1.21289 - 1.000 = \mathbf{0.213}$$

---

Step 3: Kay's Rule for Natural Gas Blend

Calculate pseudo-critical properties using Kay's rule:

$$T_{pc} = y_1 T_{c1} + y_2 T_{c2} = (0.700 \times 190.6\text{ K}) + (0.300 \times 369.8\text{ K}) = 133.42 + 110.94 = \mathbf{244.36\text{ K}} \approx 244.4\text{ K}$$ $$P_{pc} = y_1 P_{c1} + y_2 P_{c2} = (0.700 \times 46.0\text{ bar}) + (0.300 \times 42.5\text{ bar}) = 32.20 + 12.75 = \mathbf{44.95\text{ bar}} \approx 45.0\text{ bar}$$

Compute pseudo-reduced properties at separator operating conditions ($T = 320.0\text{ K}, P = 44.0\text{ bar}$):

$$T_{pr} = \frac{T}{T_{pc}} = \frac{320.0\text{ K}}{244.36\text{ K}} = \mathbf{1.31}$$ $$P_{pr} = \frac{P}{P_{pc}} = \frac{44.0\text{ bar}}{44.95\text{ bar}} = \mathbf{0.98}$$

(These pseudo-reduced values can now be used directly with generalized compressibility charts to find $Z \approx 0.84$ for vessel sizing).

---

7. Critical PE Exam Traps & Pitfalls