Introduction & Context
The calculation compares steam distillation and vacuum distillation for the separation of a low‑boiling organic component (e.g., limonene) from water. In process engineering, these two options are evaluated to determine the most economical operating pressure, energy consumption, and equipment sizing for a given product throughput. The methodology is applicable to batch or continuous reactors where the product rate, feed composition, and thermal properties are known, and where the choice of operating pressure directly influences the required reboiler duty and auxiliary utilities (steam or electricity).
Methodology & Formulas
1. Physical Constants and Input Variables
The following symbols are used throughout the derivations: R – universal gas constant (kJ·mol⁻¹·K⁻¹), MWwater – molecular weight of water (kg·kmol⁻¹), MWoil – molecular weight of the organic component (kg·kmol⁻¹), ΔHvap – molar heat of vaporisation of the oil (kJ·mol⁻¹), Tmax – maximum allowable temperature (°C), Psteam – supplied steam pressure (bar), hsteam – specific enthalpy of saturated steam (kJ·kg⁻¹), Csteam – steam cost (€/GJ), Celec – electricity cost (€/kWh), and Rref – external reflux ratio (dimensionless).
2. Antoine Equation for Saturation Pressure
The saturation pressure of a component at temperature T (°C) is expressed by the Antoine correlation:
\[ p^{*}_{i}(T) = 10^{\,A_{i} - \dfrac{B_{i}}{T + C_{i}}}\;\text{mmHg} \]Conversion to bar is performed with the factor 0.00133322 bar · mmHg⁻¹:
\[ p^{*}_{i,\; \text{bar}}(T) = p^{*}_{i}(T)\times 0.00133322 \]3. System Pressure for Steam Distillation
The target boiling temperature for steam distillation is selected as Tb,steam (°C) such that Tb,steam < Tmax. The corresponding saturation pressures of water and oil are:
\[ p^{*}_{\text{water}} = p^{*}_{\text{water, bar}}(T_{b,\;steam}),\qquad p^{*}_{\text{oil}} = p^{*}_{\text{oil, bar}}(T_{b,\;steam}) \]The total system pressure is the sum of the partial pressures:
\[ P_{\text{sys,steam}} = p^{*}_{\text{water}} + p^{*}_{\text{oil}} \]4. Vapor‑Phase Composition and Required Steam Flow
The mole fraction of oil in the vapor phase is:
\[ y_{\text{oil}} = \frac{p^{*}_{\text{oil}}}{P_{\text{sys,steam}}} \]The corresponding mass ratio of water to oil in the vapor is obtained from the ideal‑gas relation and molecular‑weight conversion:
\[ \text{mass\_ratio} = \frac{y_{\text{oil}}}{1 - y_{\text{oil}}}\;\times\;\frac{MW_{\text{water}}}{MW_{\text{oil}}} \]The required steam mass flow rate to achieve the desired product rate Fp (kg·h⁻¹) is:
\[ \dot{m}_{\text{steam}} = F_{p}\;\times\;\text{mass\_ratio} \]5. Steam‑Side Reboiler Duty and Cost
The thermal duty supplied by steam is:
\[ Q_{\text{steam}} = \dot{m}_{\text{steam}}\; \times\; h_{\text{steam}}\quad (\text{kJ·h}^{-1}) \]Conversion to gigajoules per hour and cost calculation:
\[ Q_{\text{steam,GJ}} = Q_{\text{steam}}\times g2GJ,\qquad C_{\text{steam}} = Q_{\text{steam,GJ}}\times C_{\text{steam}} \]6. Vacuum Distillation – System Pressure from Clausius‑Clapeyron
For a target boiling temperature Tb,vac (°C) the required system pressure is obtained from the integrated Clausius‑Clapeyron relation:
\[ \ln P_{\text{sys,vac}} = \frac{\Delta H_{\text{vap}}}{R}\!\left(\frac{1}{T_{\text{ref}}}-\frac{1}{T_{b,\;vac,K}}\right) + \ln P_{\text{ref}} \]where Tref and Pref are a reference temperature (K) and pressure (bar) for which the vapor pressure is known, and Tb,vac,K = Tb,vac + 273.15 K.
7. Vacuum Reboiler Duty
The specific latent heat per kilogram of oil is:
\[ \lambda = \frac{\Delta H_{\text{vap}}\times 1000}{MW_{\text{oil}}}\quad (\text{kJ·kg}^{-1}) \]The total reboiler duty, accounting for the external reflux ratio, is:
\[ Q_{\text{vac}} = F_{p}\;(1+R_{\text{ref}})\;\lambda\quad (\text{kJ·h}^{-1}) \]Expressed as electrical power (kW) and cost:
\[ P_{\text{vac}} = \frac{Q_{\text{vac}}}{kWh2kJ},\qquad C_{\text{vac}} = P_{\text{vac}}\; \times\; C_{\text{elec}} \]8. Vacuum Pump Power Requirement
An air leakage rate \dot{m}_{\text{air}} (kg·h⁻¹) is assumed. Its molar flow is:
\[ \dot{n}_{\text{air}} = \frac{\dot{m}_{\text{air}}}{MW_{\text{air}}} \]The volumetric flow at suction conditions (temperature Tsuction, pressure Psuction) is approximated by the ideal‑gas law using the standard molar volume (22.4 m³·kmol⁻¹ at 273.15 K and 1.013 bar):
\[ V_{\text{suction}} = \dot{n}_{\text{air}}\; 22.4\; \frac{T_{\text{suction}}}{273.15}\; \frac{1.013}{P_{\text{suction}}} \]The pump shaft power is estimated with a proportionality constant (3.5 kW·m³⁻¹·h) and the logarithmic pressure ratio:
\[ P_{\text{pump}} = 3.5\; \frac{V_{\text{suction}}\;\ln\!\left(\dfrac{P_{\text{discharge}}}{P_{\text{suction}}}\right)}{1000} \]where Pdischarge = 1 bar. The associated electricity cost is:
\[ C_{\text{pump}} = P_{\text{pump}}\; \times\; C_{\text{elec}} \]9. Total Vacuum Distillation Cost
\[ C_{\text{total,vac}} = C_{\text{vac}} + C_{\text{pump}} \]10. Validity Checks
| Check | Condition | Action |
|---|---|---|
| Steam‑Distillation Temperature Range | \(10 \le T_{b,\;steam} \le 120\) °C | Issue warning if violated |
| Steam‑Distillation System Pressure Range | \(0.02 \le P_{\text{sys,steam}} \le 2.0\) bar | Issue warning if violated |
| Vacuum‑Distillation Pressure Range | \(0.001 \le P_{\text{sys,vac}} \le 1.0\) bar | Issue warning if violated |
11. Summary of Output Quantities
- Steam‑distillation: boiling temperature Tb,steam, system pressure Psys,steam, required steam flow \(\dot{m}_{\text{steam}}\), and hourly steam cost \(C_{\text{steam}}\).
- Vacuum‑distillation: boiling temperature Tb,vac, system pressure Psys,vac, reboiler power \(P_{\text{vac}}\), pump power \(P_{\text{pump}}\), and total hourly electricity cost \(C_{\text{total,vac}}\).