Reference ID: MET-831F | Process Engineering Reference Sheets Calculation Guide
Introduction & Context
The Hildebrand solubility parameter (δ) is a fundamental thermodynamic property used in process engineering to quantify the solvation power or polarity of a solvent. In the context of Supercritical Fluid (SCF) extraction, this parameter is critical for predicting the efficiency of solute extraction from a solid or liquid matrix. Because SCFs exhibit gas-like diffusivity and liquid-like density, their solvation power is highly tunable via pressure and temperature adjustments.
This calculation is typically employed in the design and optimization of supercritical extraction processes, such as the decaffeination of coffee, the extraction of essential oils, or the removal of contaminants from polymers. By calculating δ, engineers can determine if a specific SCF (such as CO2) is suitable for dissolving a target compound based on the principle of "like dissolves like."
Methodology & Formulas
The calculation relies on an empirical correlation that relates the solubility parameter to the critical properties of the fluid and its current density. The density of the SCF is the primary lever for controlling solvation power; as the fluid density increases near the critical point, the solubility parameter increases accordingly. For processes that incorporate co‑solvents to further boost solubility, refer to the co‑solvent enhancement calculation.
The governing equation for the Hildebrand solubility parameter is defined as:
\(k\): Empirical scaling constant (dimensionless). For common SCFs like carbon dioxide, \(k \approx 1.25\).
\(P_{c}\): Critical pressure of the SCF [MPa]
\(\rho_{SCF}\): Density of the supercritical fluid at operating conditions [g/cm3]
\(\rho_{ref}\): Reference liquid density of the SCF at its normal boiling point [g/cm3]
To ensure thermodynamic consistency and the validity of the empirical model, the following conditions must be satisfied:
Parameter
Condition / Constraint
Operating Temperature
\(T_{op} > T_{c}\)
Operating Pressure
\(P_{op} > P_{c}\)
Reduced Pressure
\(0.9 < \frac{P_{op}}{P_{c}} < 2.5\)
Density Ratio
\(\frac{\rho_{SCF}}{\rho_{ref}} \leq 2.0\)
Note: The accuracy of the resulting solubility parameter is entirely dependent on the precision of the supercritical fluid density (\(\rho_{SCF}\)). It is strongly recommended to obtain this value from a reliable real-fluid equation of state (e.g., Peng-Robinson) or an NIST-validated thermophysical property database rather than assuming ideal gas behavior.
To calculate the cohesive energy density for a supercritical fluid (SCF) system, you must evaluate the energy required to vaporize the fluid at a specific temperature and pressure. Follow these steps:
Calculate the molar enthalpy of vaporization at the target state.
Subtract the product of pressure and molar volume (RT) to account for the work of expansion.
Divide the resulting cohesive energy by the molar volume of the fluid.
Take the square root of the final value to obtain the Hildebrand solubility parameter.
In supercritical fluids, the solubility parameter is highly dependent on density because the intermolecular distance changes significantly with small variations in pressure or temperature. Since the cohesive energy density is defined as the square of the solubility parameter, even minor fluctuations in fluid density lead to non-linear changes in the solvent power of the SCF.
While several models exist, the following equations are generally preferred for process engineering applications:
Peng-Robinson (PR) equation of state for non-polar or slightly polar fluids.
Soave-Redlich-Kwong (SRK) for hydrocarbon-based SCF systems.
PC-SAFT for complex fluids or systems involving polymers where molecular shape and association effects are significant.
For polar SCFs, the total solubility parameter should be partitioned into individual components using the Hansen solubility parameter framework. You must sum the squares of the individual contributions:
Dispersion forces (\(\delta_d\))
Polar interactions (\(\delta_p\))
Hydrogen bonding (\(\delta_h\))
The total parameter is calculated as \(\delta_{total} = \sqrt{\delta_d^2 + \delta_p^2 + \delta_h^2}\).
Worked Example: Solubility Parameter Calculation for Supercritical CO₂
A process engineer is designing a supercritical fluid extraction system using carbon dioxide (CO₂) to recover non-polar essential oils from plant material. To evaluate the solvent's solvation power, the Hildebrand solubility parameter (δ) must be calculated at the proposed operating conditions.
Critical temperature of CO₂, \( T_{c} = 31.0 \, ^\circ\text{C} \)
Reference liquid density of CO₂ (at normal boiling point), \( \rho_{ref} = 0.85 \, \text{g/cm}^3 \)
Supercritical fluid density of CO₂ at \( T_{op} \) and \( P_{op} \), \( \rho_{SCF} = 0.75 \, \text{g/cm}^3 \) (from thermophysical property data)
Step-by-Step Calculation:
Compute the density ratio:
\[ \frac{\rho_{SCF}}{\rho_{ref}} = \frac{0.75}{0.85} = 0.882 \]
(This value is rounded from 0.88235 for display.)
Compute the square root of the critical pressure:
\[ \sqrt{P_{c}} = \sqrt{7.38 \, \text{MPa}} = 2.717 \, \text{MPa}^{0.5} \]
(This value is rounded from 2.71662 for display.)
Apply the empirical correlation to calculate the Hildebrand solubility parameter:
\[ \delta = k \cdot \sqrt{P_{c}} \cdot \left( \frac{\rho_{SCF}}{\rho_{ref}} \right) \]
where \( k = 1.25 \) is the empirical scaling constant for CO₂.
Substituting the values:
\[ \delta = 1.25 \cdot 2.717 \cdot 0.882 = 2.996 \, \text{MPa}^{0.5} \]
(Using the displayed rounded intermediates yields δ ≈ 2.996. The precise calculation gives 2.9960 MPa⁰·⁵.)
Final Answer: The Hildebrand solubility parameter for supercritical CO₂ at the given conditions is \( \delta = 2.996 \, \text{MPa}^{0.5} \).
"Un projet n'est jamais trop grand s'il est bien conçu."— André Citroën
"La difficulté attire l'homme de caractère, car c'est en l'étreignant qu'il se réalise."— Charles de Gaulle
