Aroma Recovery in Distillation

Reference ID: MET-DC57 | Process Engineering Reference Sheets Calculation Guide

Introduction & Context

Aroma recovery in distillation is a critical unit operation in the food, beverage, and fragrance industries. The objective is to isolate volatile organic compounds from aqueous matrices while preserving their organoleptic properties. Because many aroma compounds are thermally sensitive, these systems are typically operated under vacuum to maintain low operating temperatures (typically below 60 °C); applying proven thermal damage and aroma retention strategies further protects product quality. This calculation blueprint provides a systematic approach to designing a distillation column for dilute aroma recovery, ensuring high product purity and recovery efficiency while adhering to strict thermal constraints.

Methodology & Formulas

The design process follows a sequential engineering approach, transitioning from mass balances to stage-wise equilibrium analysis.

1. Mass Balance and Recovery

The system is defined by the total feed flow rate \(\dot{m}_{F}\) and the mass fractions of the key aroma component in the feed (\(z_{i}\)), distillate (\(x_{D,i}\)), and bottoms (\(x_{B,i}\)). The mass balance is governed by:

\[ \dot{m}_{F} = \dot{m}_{D} + \dot{m}_{B} \]

\[ \dot{m}_{F} \cdot z_{i} = \dot{m}_{D} \cdot x_{D,i} + \dot{m}_{B} \cdot x_{B,i} \]

The recovery efficiency (\(R_{i}\)) is calculated as:

\[ R_{i} = \frac{\dot{m}_{D} \cdot x_{D,i}}{\dot{m}_{F} \cdot z_{i}} \]

2. Equilibrium and Minimum Stages

Mole fractions (\(x_{\text{mole}}\)) are derived from mass fractions using the molecular weights (\(MW_{i}\)) of the aroma compound and water. The minimum number of theoretical stages (\(N_{\text{min}}\)) is determined using the Fenske equation, assuming constant relative volatility (\(\alpha\)):

\[ N_{\text{min}} = \frac{\ln\left( \dfrac{x_{D,\text{LK}}}{x_{D,\text{HK}}} \cdot \dfrac{x_{B,\text{HK}}}{x_{B,\text{LK}}} \right)}{\ln(\alpha)} \]

3. Minimum Reflux and Actual Stages

For a saturated liquid feed (\(q=1\)), the Underwood equation is used to find the root (\(\theta\)) and the minimum reflux ratio (\(R_{\text{min}}\)):

\[ \frac{\alpha \cdot z_{\text{LK}}}{\alpha - \theta} + \frac{1 \cdot z_{\text{HK}}}{1 - \theta} = 0 \]

\[ R_{\text{min}} = \left( \sum \frac{\alpha_{i} \cdot x_{D,i}}{\alpha_{i} - \theta} \right) - 1 \]

The actual number of stages (\(N_{\text{actual}}\)) is estimated using a simplified Gilliland-type correlation, incorporating a safety factor (\(SF\)) to account for non-ideal column performance. Note: The form \(\sqrt{X_{\text{Gill}}}\) is a conservative approximation of the full Gilliland curve; for rigorous design, use the complete correlation.

\[ N_{\text{actual}} = (N_{\text{min}} + \sqrt{X_{\text{Gill}}}) \cdot SF \]

Where \(X_{\text{Gill}}\) is the Gilliland abscissa defined by the operating reflux ratio (\(R_{\text{op}}\)):

\[ X_{\text{Gill}} = \frac{R_{\text{op}} - R_{\text{min}}}{R_{\text{op}} + 1} \]

Parameter	Regime / Constraint	Threshold / Limit
Relative Volatility	Empirical Validity	\(1.2 \leq \alpha \leq 4.0\)
Bottom Purity	Product Specification	\(x_{B,i} \leq 0.001\) (mass fraction)
Reflux Ratio	Operating / Economic	\(1.05 \leq \dfrac{R_{\text{op}}}{R_{\text{min}}} \leq 2.0\)
Thermal Limit	Aroma Degradation	\(T_{\text{bottom}} < 60^{\circ}\text{C}\)

To balance aroma retention with purity, process engineers should focus on the following adjustments:

Monitor the concentration of the key volatile aroma compounds in the distillate using real-time gas chromatography.
Reduce the reflux ratio incrementally while maintaining the target aroma concentration to prevent the loss of light-end volatile fractions.
Adjust the feed tray location to ensure that high-boiling impurities are effectively separated from the desired aromatic compounds.

Thermal degradation often manifests through specific chemical and physical changes in the process stream:

A noticeable shift in the organoleptic profile, often described as cooked, burnt, or oxidized notes.
An increase in the concentration of degradation byproducts such as oxidation products (e.g., limonene oxide from terpenes) or furfural from sugar decomposition.
A decrease in the recovery yield of heat-sensitive esters and terpenes as measured by the component mass balance closure.

Operating under vacuum is critical for preserving delicate aromatics because it lowers the boiling point of the mixture. By reducing the operating pressure, you can:

Minimize the residence time at high temperatures, which prevents the thermal breakdown of volatile organic compounds.
Improve the relative volatility of key aromatic components, allowing for more efficient separation at lower energy inputs.
Prevent the formation of off-flavors that occur at atmospheric boiling temperatures.

Worked Example: Minimum Stages and Packing Height for Citrus Essence Recovery

A fruit juice processing plant produces 1000 kg/h of an aqueous distillate stream at 50°C containing 200 ppm (0.02 mass%) of limonene, the key citrus aroma compound. The objective is to recover at least 95% of the limonene into a distillate product containing 1 mass% limonene, while keeping the bottom concentration at or below 10 ppm. The column operates under vacuum at 0.2 bar to maintain temperatures below 60°C, with an average relative volatility of limonene to water of 2.5. Structured packing with an HETP of 0.4 m will be used.

Knowns (Input Parameters and Units):

Feed mass flow rate, \(F = 1000.0 \; \text{kg/h}\)
Feed limonene mass fraction, \(z_{\text{lim,mass}} = 0.0002\) (200 ppm)
Distillate limonene mass fraction (target), \(x_{D,\text{lim,mass}} = 0.01\)
Maximum bottom limonene mass fraction, \(x_{B,\text{lim,mass,spec}} = 1 \times 10^{-5}\) (10 ppm)
Molecular weight of limonene, \(MW_{\text{limonene}} = 136.0 \; \text{g/mol}\)
Molecular weight of water, \(MW_{\text{water}} = 18.0 \; \text{g/mol}\)
Average relative volatility (limonene/water), \(\alpha_{\text{avg}} = 2.5\)
Column operating pressure, \(P_{\text{col}} = 0.2 \; \text{bar}\)
Packing HETP, \(\text{HETP} = 0.4 \; \text{m}\)
Safety factor for actual stages, \(SF = 1.2\)

Step-by-Step Calculation:

Material Balance. The limonene mass flow in the feed is: \(M_{\text{lim},F} = F \cdot z_{\text{lim,mass}} = 0.200 \; \text{kg/h}\). Using the specified distillate and bottom purities, the distillate flow rate is: \(D = (M_{\text{lim},F} - F \cdot x_{B,\text{lim,mass,spec}}) / (x_{D,\text{lim,mass}} - x_{B,\text{lim,mass,spec}}) = (0.200 - 0.010) / (0.01 - 1\times10^{-5}) = 19.019 \; \text{kg/h}\). The bottoms flow rate is: \(B = F - D = 980.981 \; \text{kg/h}\). The limonene mass in the distillate is: \(M_{\text{lim},D} = D \cdot x_{D,\text{lim,mass}} = 0.1902 \; \text{kg/h}\), and in the bottoms (from specification): \(M_{\text{lim},B} = B \cdot x_{B,\text{lim,mass,spec}} = 0.00981 \; \text{kg/h}\). This gives an actual recovery of: \(R_{\text{recovery}} = M_{\text{lim},D} / M_{\text{lim},F} = 0.951\) (95.1%). The computed bottom limonene mass fraction is \(x_{B,\text{lim,mass}} = 1.0 \times 10^{-5}\), which matches the specification.
Mole Fraction Conversions. All mass-based compositions are converted to mole fractions using the molecular weights.
- Feed: Moles of limonene: \(N_{\text{lim},F} = 0.200 / 136.0 = 0.001471 \; \text{kmol/h}\). Moles of water: \(N_{\text{water},F} = (1000.0 - 0.200) / 18.0 = 55.544 \; \text{kmol/h}\). Total moles: \(N_{\text{total},F} = 55.546 \; \text{kmol/h}\). Feed mole fractions: \(z_{\text{lim,mole}} = 2.648 \times 10^{-5}\), \(z_{\text{water,mole}} = 0.99997\).
- Distillate: Moles of limonene: \(N_{\text{lim},D} = 0.1902 / 136.0 = 0.001399 \; \text{kmol/h}\). Moles of water: \(N_{\text{water},D} = (19.019 - 0.1902) / 18.0 = 1.046 \; \text{kmol/h}\). Distillate mole fractions: \(x_{D,\text{lim,mole}} = 0.001335\), \(x_{D,\text{water,mole}} = 0.99866\).
- Bottoms: Moles of limonene: \(N_{\text{lim},B} = 0.00981 / 136.0 = 7.213 \times 10^{-5} \; \text{kmol/h}\). Moles of water: \(N_{\text{water},B} = (980.981 - 0.00981) / 18.0 = 54.498 \; \text{kmol/h}\). Bottom mole fractions: \(x_{B,\text{lim,mole}} = 1.324 \times 10^{-6}\), \(x_{B,\text{water,mole}} = 0.999999\).
Minimum Theoretical Stages (Fenske Equation). The Fenske equation is applied using the distillate and bottom mole fractions of the light key (limonene) and heavy key (water):
\[ N_{\text{min}} = \frac{\ln\left[\left(\dfrac{x_{D,\text{lim,mole}}}{x_{D,\text{water,mole}}}\right) \left(\dfrac{x_{B,\text{water,mole}}}{x_{B,\text{lim,mole}}}\right)\right]}{\ln(\alpha_{\text{avg}})} \]
Substituting the mole fraction values and \(\alpha_{\text{avg}} = 2.5\):
\[ N_{\text{min}} = \frac{\ln\left[\left(\dfrac{0.001335}{0.99866}\right) \left(\dfrac{0.999999}{1.324 \times 10^{-6}}\right)\right]}{\ln(2.5)} = 7.550 \; \text{theoretical stages (including reboiler)}. \]
Minimum Reflux Ratio (Underwood Method, Saturated Liquid Feed, q=1). For a binary system with saturated liquid feed, the Underwood root \(\theta\) is found analytically. Note: For extremely dilute feeds, \(\theta\) approaches \(\alpha\) very closely; high numerical precision is required.
\[ \theta = \frac{\alpha_{\text{avg}}}{\alpha_{\text{avg}} \cdot z_{\text{lim,mole}} + z_{\text{water,mole}}} = \frac{2.5}{2.5 \times 2.648 \times 10^{-5} + 0.99997} \approx 2.4999 \]
The minimum reflux ratio is then:
\[ R_{\text{min}} = \frac{\alpha_{\text{avg}} \cdot x_{D,\text{lim,mole}}}{\alpha_{\text{avg}} - \theta} + \frac{1.0 \cdot x_{D,\text{water,mole}}}{1.0 - \theta} - 1 = 31.96 \]
(Note: The high \(R_{\text{min}}\) reflects the extreme dilution and pinch condition near \(\theta \approx \alpha\); in practice, a much lower operating reflux is used, as dictated by design heuristics for dilute systems.)
Actual Stages and Packing Height (Simplified Gilliland Correlation). An operating reflux ratio is set at \(R_{\text{op}} = 2.0 \times R_{\text{min}} = 63.91\). The Gilliland parameter is:
\[ X_{\text{Gill}} = \frac{R_{\text{op}} - R_{\text{min}}}{R_{\text{op}} + 1} = 0.492 \]
The actual number of theoretical stages, including the safety factor, is:
\[ N_{\text{actual}} = (N_{\text{min}} + \sqrt{X_{\text{Gill}}}) \times SF = (7.550 + \sqrt{0.492}) \times 1.2 = 9.902 \]
The required packing height is:
\[ H_{\text{packing}} = N_{\text{actual}} \times \text{HETP} = 9.902 \times 0.4 = 3.961 \; \text{m} \]

Final Answer: The column requires a minimum of 7.550 theoretical stages (including reboiler) and, with a safety factor of 1.2, an actual of 9.902 stages. Using structured packing with an HETP of 0.4 m, the total packing height needed is 3.961 m. The design achieves a limonene recovery of 95.1% into the distillate, meeting the 10 ppm bottom specification.

"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