Introduction & Context
Steam distillation is a critical separation technique in process engineering used to isolate heat-sensitive organic compounds from non-volatile or thermally unstable matrices. By introducing steam into a mixture of immiscible liquids, the system total vapor pressure is increased, allowing the mixture to boil at a temperature significantly lower than the boiling point of the pure organic component or water alone.
This calculation is essential for designing distillation columns, sizing heat exchangers, and ensuring safety in batch processing. It is primarily used in the essential oil, fragrance, and pharmaceutical industries where thermal degradation must be avoided during extraction.
Methodology & Formulas
The boiling point of an immiscible mixture is reached when the sum of the partial vapor pressures of the individual components equals the external system pressure. The vapor pressure of each component is determined using the Antoine equation.
The governing relationship for the system is defined as:
\[ P_{atm} = P_{water}^{sat}(T) + P_{organic}^{sat}(T) \]
The boiling point of an immiscible mixture is reached when the sum of the partial vapor pressures of the individual components equals the external system pressure. Understanding the vapor composition in steam distillation is essential for accurate calculations, as it determines how each component contributes to the total pressure. The vapor pressure of each component is determined using the Antoine equation.
\[ \log_{10}(P^{sat}) = A - \frac{B}{T + C} \]
Which is rearranged to solve for pressure as:
\[ P^{sat} = 10^{\left( A - \frac{B}{T + C} \right)} \]
Because the temperature \(T\) is implicit within the Antoine equations, a numerical approach such as the Bisection Method is employed to find the root where the calculated total pressure matches the atmospheric pressure. The error function to be minimized is:
\[ \epsilon = P_{atm} - \left( P_{water}^{sat}(T) + P_{organic}^{sat}(T) \right) \]
| Condition |
Requirement |
| Phase State |
Liquids must be mutually immiscible (forming two distinct phases). |
| Temperature Range |
\(T_{min} \leq T \leq T_{max}\) (Must remain within Antoine constant validity limits). |
| Convergence Criteria |
\(|P_{atm} - (P_{water}^{sat} + P_{organic}^{sat})| \leq \text{Tolerance}\) |
| Thermal Stability |
Component must not undergo chemical decomposition at \(T_{boil}\). |
The boiling point of a steam distillation system is reached when the sum of the vapor pressures of the immiscible components equals the total atmospheric pressure. Because the components exert pressure independently, the mixture boils at a temperature lower than the boiling point of either pure component. The calculation follows these steps:
- Identify the vapor pressure of water at various temperatures.
- Identify the vapor pressure of the organic compound at those same temperatures.
- Sum the two pressures for each temperature point.
- Locate the temperature where the sum of the pressures equals the ambient system pressure.
Worked Example: Steam Distillation Boiling Point of Water and d‑Limonene
Scenario: A process engineer must determine the boiling point temperature of an immiscible mixture of water and d‑limonene at atmospheric pressure. The vapor phase is assumed ideal, and the total pressure is the sum of the pure‑component saturation pressures (Dalton’s law).
Known Parameters
- Atmospheric pressure, \(P_{\text{atm}} = 1.01325\ \text{bar}\)
- Antoine constants for water (valid 0 – 100 °C, pressure in bar): \(A_{\text{water}}=5.196\), \(B_{\text{water}}=1730.63\), \(C_{\text{water}}=233.43\)
- Antoine constants for d‑limonene (valid 50 – 120 °C, pressure in bar): \(A_{\text{limonene}}=4.626\), \(B_{\text{limonene}}=1535.0\), \(C_{\text{limonene}}=237.3\)
- Lower temperature bound for bisection: \(T_{\text{low}} = 50.0\ \text{°C}\)
- Upper temperature bound for bisection: \(T_{\text{high}} = 100.0\ \text{°C}\)
- Sum of saturation pressures at \(50.0\ \text{°C}\): \(P_{\text{sum,low}} = 0.315\ \text{bar}\)
- Sum of saturation pressures at \(100.0\ \text{°C}\): \(P_{\text{sum,high}} = 2.202\ \text{bar}\)
- Number of bisection iterations: \(20\)
Step‑by‑Step Calculation
- Define vapor‑pressure equations.
Water: \(P_{\text{water}}^{\text{sat}}(T) = 10^{\,A_{\text{water}} - \frac{B_{\text{water}}}{T + C_{\text{water}}}}\)
d‑Limonene: \(P_{\text{limonene}}^{\text{sat}}(T) = 10^{\,A_{\text{limonene}} - \frac{B_{\text{limonene}}}{T + C_{\text{limonene}}}}\)
- Verify that a root exists between \(T_{\text{low}}\) and \(T_{\text{high}}\).
At \(T_{\text{low}} = 50.0\ \text{°C}\): \(P_{\text{sum,low}} = 0.315\ \text{bar} < P_{\text{atm}}\).
At \(T_{\text{high}} = 100.0\ \text{°C}\): \(P_{\text{sum,high}} = 2.202\ \text{bar} > P_{\text{atm}}\).
Therefore a boiling point exists in the interval \([50.0,100.0]\ \text{°C}\).
- Perform bisection.
Starting with \(T_a = T_{\text{low}}\) and \(T_b = T_{\text{high}}\), for 20 iterations:
\(T_{\text{mid}} = \frac{T_a + T_b}{2}\)
Evaluate \(P_{\text{water}}^{\text{sat}}(T_{\text{mid}})\) and \(P_{\text{limonene}}^{\text{sat}}(T_{\text{mid}})\).
If the sum is less than \(P_{\text{atm}}\), set \(T_a = T_{\text{mid}}\); otherwise set \(T_b = T_{\text{mid}}\).
- Final converged value.
After iteration \(i = 19\), the mid‑point is \(T_{\text{mid}} = 78.143\ \text{°C}\).
At this temperature:
\(P_{\text{water}}^{\text{sat}} = 0.438\ \text{bar}\)
\(P_{\text{limonene}}^{\text{sat}} = 0.575\ \text{bar}\)
\(P_{\text{total}} = 0.438 + 0.575 = 1.013\ \text{bar}\)
Error = \(P_{\text{atm}} - P_{\text{total}} = 1.01325 - 1.013 = 0.000\ \text{bar}\)
Final Answer
The steam distillation boiling point of the water + d‑limonene mixture at 1.01325 bar is 78.143 °C.
