Reference ID: MET-A031 | Process Engineering Reference Sheets Calculation Guide
Introduction & Context
The Fixed Bed Extractor design is a critical unit operation in the food and beverage industry, specifically for the production of instant coffee. This process utilizes a percolation battery—a series of interconnected columns—to achieve high‑efficiency solid‑liquid extraction and aligns with the overall coffee solubles extraction process design. By operating in a quasi‑continuous manner, the system ensures that fresh solvent (hot water) contacts the most exhausted coffee grounds, while the most concentrated extract is produced by contacting fresh coffee grounds. This counter‑current approach maximizes the extraction yield and concentration, which is essential for downstream processing such as evaporation and spray drying.
Methodology & Formulas
The design of a percolation battery relies on balancing hydraulic constraints with mass transfer kinetics. The following formulas define the operational parameters based on the physical properties of the bed and the desired throughput.
1. Cycle and Temporal Parameters
The cycle time \( t_{c} \) is determined by the mass of the coffee charge per column \( M_{\text{bed}} \) and the total coffee feed rate \( \dot{m}_{\text{coffee}} \):
\[ t_{c} = \frac{M_{\text{bed}}}{\dot{m}_{\text{coffee}}} \]
The switching frequency \( f \) and total extraction time \( t_{\text{extract}} \) for a system with \( N_{\text{extract}} \) active columns are defined as:
\[ f = \frac{1}{t_{c}} \]
\[ t_{\text{extract}} = N_{\text{extract}} \cdot t_{c} \]
2. Hydraulic and Velocity Calculations
The solvent-to-feed ratio \( S/F \) ensures the mass balance of the system. The superficial velocity \( u_{s} \) and interstitial velocity \( u_{i} \) are calculated based on the column cross-sectional area \( A_{\text{bed}} \), water density \( \rho_{\text{water}} \), and bed porosity \( \epsilon \):
\[ S/F = \frac{\dot{m}_{\text{water}}}{\dot{m}_{\text{coffee}}} \]
\[ u_{s} = \frac{\dot{m}_{\text{water,col}}}{\rho_{\text{water}} \cdot A_{\text{bed}}} \]
\[ u_{i} = \frac{u_{s}}{\epsilon} \]
3. Pressure Drop and Mass Transfer
The pressure drop \( \Delta P \) across the bed is governed by Darcy's Law for laminar flow through porous media, where \( \mu \) is the dynamic viscosity, \( L \) is the bed length, and \( k \) is the permeability:
\[ \Delta P = \left( \frac{\mu \cdot L}{k} \right) \cdot u_{s} \]
The extraction yield \( Y(t) \) is modeled using the Weibull distribution, where \( k_{\text{ext}} \) is the rate constant, \( \beta \) is the shape factor, and \( Y_{\infty} \) is the maximum theoretical yield:
\[ Y(t) = Y_{\infty} \cdot \left[ 1 - \exp\!\left( -(k_{\text{ext}} \cdot t)^{\beta} \right) \right] \]
Operational Validity Criteria
Parameter
Symbol
Valid Range
Solvent-to-Feed Ratio
\( S/F \)
4.0 – 8.0
Interstitial Velocity
\( u_{i} \)
0.5 – 2.5 mm/s
Pressure Drop
\( \Delta P \)
0.5 – 2.5 bar
Bed Aspect Ratio
\( L/D \)
3.0 – 6.0
Operating Temperature
\( T \)
140 – 180 °C
To ensure uniform extraction efficiency and prevent channeling, process engineers must optimize the following parameters:
Bed porosity and particle size distribution of the solid matrix.
Design of the liquid distributor plate to ensure even wetting across the cross-sectional area.
Superficial velocity of the solvent relative to the bed voidage.
Viscosity and surface tension characteristics of the solvent at operating temperatures.
The aspect ratio is critical for maintaining plug flow and minimizing axial dispersion. Engineers should consider:
Maintaining a minimum ratio to prevent wall effects from dominating the flow profile.
Accounting for the pressure drop limitations across the bed depth.
Ensuring the residence time distribution meets the required extraction kinetics.
Evaluating the mechanical integrity of the bed to prevent compaction under high solvent flow rates.
Accurate pressure drop estimation is essential for pump sizing and preventing bed fluidization. Standard practice involves:
Applying Darcy's Law for laminar flow through the packed bed, which relates pressure drop to superficial velocity, fluid viscosity, bed length, and permeability.
For higher Reynolds number regimes where inertial losses become significant, the Ergun equation provides a more comprehensive correlation that accounts for both viscous and inertial contributions.
Adjusting the calculation based on the sphericity factor of the solid particles when using the Ergun equation.
Incorporating a safety factor to account for potential fines migration or bed settling over time.
Verifying the calculated pressure drop against the structural load limits of the support grid.
Worked Example: Fixed Bed Extractor Design (Instant Coffee Percolation Battery)
Scenario: A 6-column battery operates with 4 columns on active extraction, 1 draining, and 1 filling/discharging. The column geometry and feed rates are specified below.
Column geometry: Bed length \(L = 5.0\) m, bed diameter \(D = 1.0\) m → cross-sectional area \(A_{\text{bed}} = 0.785\) m²
Coffee charge per column: \(M_{\text{bed}} = 1500\) kg
Pressure drop \(\Delta P\) (Darcy's Law):
\(\Delta P = \left( \dfrac{\mu \cdot L}{k} \right) \cdot u_{s}\)
→ \(\Delta P = 222283\) Pa = 2.223 bar (within valid range 0.5–2.5 bar)
L/D ratio check:
\(L/D = 5.0\) (within valid range 3–6)
Final Answer: The extraction yield is 0.449 (44.9%) and the extract concentration is 0.101 (10.1% w/w). All operating parameters (interstitial velocity, pressure drop, S/F ratio, L/D ratio) fall within the recommended empirical ranges.
"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