Particle Size Effect on Extraction Rate

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

Introduction & Context

Particle size optimization is a critical design parameter in solid-liquid extraction processes, such as tea brewing, pharmaceutical leaching, and hydrometallurgical recovery. In these unit operations, the extraction rate is governed by the competing phenomena of internal solute diffusion and external convective mass transfer. Selecting the appropriate particle size (\(d_{p}\)) requires balancing the kinetics of extraction—where smaller particles reduce diffusion time—against the hydraulic constraints of the packed bed, where smaller particles significantly increase pressure drop and the risk of bed clogging.

Methodology & Formulas

The optimization process relies on defining a feasible window for the particle diameter based on kinetic requirements and hydraulic limits. The following mathematical framework is utilized to determine the optimal particle size.

1. Kinetic Extraction Limit

The characteristic time required for extraction is dominated by internal diffusion. For diffusion from spherical particles, the time scale \(\tau\) to achieve a target extraction fraction \(F\) is approximated by:

\[ \tau = \frac{d_{p}^{2}}{C \cdot D_{eff}} \]

Where \(C\) is a dimensionless constant dependent on the target extraction fraction (e.g., \(C = \pi^{2}\) for 99% extraction from spheres), and \(D_{eff}\) is the effective diffusivity of the solute within the solid matrix.

2. Hydraulic Pressure Drop Limit

To ensure the system remains within pump capacity, the pressure drop across the bed is calculated using the Kozeny-Carman equation for laminar flow in packed beds:

\[ \frac{\Delta P}{L} = \frac{150 \cdot \mu \cdot u \cdot (1 - \epsilon)^{2}}{\epsilon^{3} \cdot d_{p}^{2}} \]

This equation dictates the minimum allowable particle size to prevent excessive resistance to flow. Valid for particle Reynolds number \(Re_{p} \leq 10\).

3. Mass Transfer and Flow Regimes

To validate the extraction regime, the system evaluates the dimensionless Reynolds (\(Re_{p}\)), Schmidt (\(Sc\)), and Sherwood (\(Sh\)) numbers:

\[ Re_{p} = \frac{\rho \cdot u \cdot d_{p}}{\mu} \]

\[ Sc = \frac{\nu}{D_{m}} = \frac{\mu}{\rho \cdot D_{m}} \]

For flow around isolated spheres (external mass transfer), the Sherwood number correlation is:

\[ Sh = \frac{k_{c} \cdot d_{p}}{D_{m}} = 2.0 + 0.6 \cdot Re_{p}^{1/2} \cdot Sc^{1/3} \]

The external mass transfer coefficient \(k_{c}\) is derived from the Sherwood number. For diffusion-controlled extraction, the internal diffusion resistance dominates when the characteristic diffusion time \(\tau_{diff} = d_{p}^{2}/(C \cdot D_{eff})\) is much greater than the external mass transfer time \(\tau_{ext} = d_{p}/(6 \cdot k_{c})\).

Regime/Constraint	Condition/Threshold
Practical Particle Size Range	Typically 0.1 mm ≤ \(d_{p}\) ≤ 10 mm for packed beds
Flow Regime (Laminar, Kozeny-Carman valid)	\(Re_{p} \leq 10\)
Bed Voidage (Typical packed bed)	0.35 ≤ ε ≤ 0.5
Diffusion Control Dominance	\(\tau_{diff} \gg \tau_{ext}\) or \(\frac{d_{p}^{2}}{C \cdot D_{eff}} \gg \frac{d_{p}}{6 \cdot k_{c}}\)

Reducing particle size increases the specific surface area available for solvent contact, which directly influences mass transfer through several mechanisms:

Increased surface area to volume ratio: The surface area per unit mass scales as \(1/d_{p}\), reducing the diffusion path length for the solute within the solid matrix.
Reduced internal mass transfer resistance: The characteristic diffusion time scales as \(d_{p}^{2}\), so halving the particle diameter reduces diffusion time by a factor of four.
External mass transfer coefficient (\(k_{c}\)): For a given fluid velocity, smaller particles typically experience lower Reynolds numbers, which may decrease \(k_{c}\) according to correlations like \(Sh = 2.0 + 0.6 Re_{p}^{1/2} Sc^{1/3}\). However, the increased surface area often outweighs any reduction in \(k_{c}\).
Faster approach to equilibrium: The combined effect of reduced diffusion time and increased interfacial area leads to more rapid extraction kinetics.

While finer particles accelerate extraction, they introduce significant mechanical and process challenges that must be managed:

Increased pressure drop: According to the Kozeny-Carman equation, \(\Delta P \propto 1/d_{p}^{2}\), so halving particle diameter increases pressure drop by a factor of four, potentially leading to channeling, bed compaction, or exceeding pump capacity.
Solid-liquid separation difficulties: Fine particles can blind filter media, reduce filtration rates, increase centrifuge loading, and require more extensive downstream processing.
Energy costs: The specific energy for size reduction increases significantly as target particle size decreases, following Kick's Law (\(E \propto 1/d_{p}\)).
Heat generation: Fine milling can generate substantial heat, potentially degrading heat-sensitive compounds in pharmaceutical or food applications.
Handling issues: Very fine powders may fluidize, creating dust hazards, or exhibit poor flow properties in hoppers and feeders.

Determining the ideal particle size requires balancing extraction kinetics against downstream processing constraints through systematic engineering analysis:

Define process objectives: Specify target extraction efficiency, cycle time, allowable pressure drop, and downstream separation requirements.
Characterize material properties: Measure effective diffusivity \(D_{eff}\), particle density, and mechanical properties affecting size reduction.
Establish kinetic constraints: Use the characteristic time equation \(\tau = d_{p}^{2}/(C \cdot D_{eff})\) to determine the maximum particle size that achieves the target extraction within the available time.
Establish hydraulic constraints: Apply the Kozeny-Carman equation \(\Delta P/L = [150 \mu u (1-\epsilon)^{2}]/(\epsilon^{3} d_{p}^{2})\) to determine the minimum particle size that keeps pressure drop within pump capacity.
Evaluate additional factors:
- Check that \(Re_{p} \leq 10\) for laminar flow assumption validity
- Compare internal vs. external mass transfer resistances
- Consider practical limits: 0.1-10 mm for typical packed beds
Perform economic optimization: Balance milling energy costs against benefits of faster extraction and downstream processing costs.
Validate experimentally: Conduct laboratory-scale tests across the calculated feasible range to confirm performance.

Worked Example: Determining Feasible Particle Size for Tea Extraction

A process engineer is designing a packed bed leaching system to extract soluble compounds from tea leaves. The goal is to identify a particle size that ensures complete extraction within a target time while keeping the pressure drop across the bed low enough for standard pumping equipment.

Knowns:

Target extraction time scale, \(\tau_{target} = 1500.0\ \text{s}\)
Effective solute diffusivity in leaf matrix, \(D_{eff} = 1.0 \times 10^{-10}\ \text{m}^{2}/\text{s}\)
Constant for 99% extraction from spheres, \(C = \pi^{2} \approx 9.87\)
Maximum allowable pressure drop, \(\Delta P_{max} = 50000.0\ \text{Pa}\)
Bed height, \(L = 0.5\ \text{m}\)
Solvent (water) dynamic viscosity, \(\mu = 0.00032\ \text{Pa} \cdot \text{s}\)
Superficial velocity, \(u = 0.005\ \text{m}/\text{s}\)
Bed void fraction, \(\epsilon = 0.45\)

Step-by-Step Calculation:

Find maximum particle size from kinetics constraint. For diffusion-controlled extraction, the characteristic time is \(\tau = d_{p}^{2} / (C \cdot D_{eff})\). Setting \(\tau = \tau_{target}\) and solving for the particle diameter gives the upper limit:

\[ d_{p,max} = \sqrt{\tau_{target} \cdot C \cdot D_{eff}} = \sqrt{1500.0 \cdot 9.87 \cdot 1.0 \times 10^{-10}} = \sqrt{1.4805 \times 10^{-6}} = 1.217 \times 10^{-3}\ \text{m} = 1.22\ \text{mm} \]

This is the largest particle size that will meet the 25-minute extraction target for 99% extraction.

Find minimum particle size from pressure drop constraint. Using the Kozeny-Carman equation for laminar flow, rearranged to solve for minimum diameter:

\[ d_{p,min} = \sqrt{ \frac{150 \cdot \mu \cdot u \cdot (1-\epsilon)^{2}}{\epsilon^{3} \cdot (\Delta P_{max} / L)} } \]

Calculate step-by-step:

Numerator: \(150 \cdot \mu \cdot u \cdot (1-\epsilon)^{2} = 150 \cdot 0.00032 \cdot 0.005 \cdot (0.55)^{2} = 7.260 \times 10^{-5}\ \text{Pa} \cdot \text{s} \cdot \text{m}\)

Denominator: \(\epsilon^{3} \cdot (\Delta P_{max} / L) = (0.45)^{3} \cdot (50000.0 / 0.5) = 0.091125 \cdot 100000 = 9112.5\ \text{Pa}/\text{m}\)

\[ d_{p,min} = \sqrt{ \frac{7.260 \times 10^{-5}}{9112.5} } = \sqrt{7.965 \times 10^{-9}} = 8.925 \times 10^{-5}\ \text{m} = 0.0893\ \text{mm} \]

Any particle smaller than this would cause excessive pressure drop exceeding 50 kPa across the 0.5 m bed.

Check flow regime constraint:

Using the geometric mean of the limits \(d_{p} = \sqrt{d_{p,min} \cdot d_{p,max}} = \sqrt{8.925 \times 10^{-5} \cdot 1.217 \times 10^{-3}} = 3.29 \times 10^{-4}\ \text{m}\)

Assuming water density \(\rho = 1000\ \text{kg}/\text{m}^{3}\):

\[ Re_{p} = \frac{\rho \cdot u \cdot d_{p}}{\mu} = \frac{1000 \cdot 0.005 \cdot 3.29 \times 10^{-4}}{0.00032} = 5.14 \leq 10 \]

The flow is laminar, validating use of the Kozeny-Carman equation.

Final Answer:

The feasible particle diameter \(d_{p}\) must satisfy \(8.93 \times 10^{-5}\ \text{m} < d_{p} < 1.22 \times 10^{-3}\ \text{m}\) (0.0893 mm < \(d_{p}\) < 1.22 mm). A reasonable design value within this range is the geometric mean \(d_{p,opt} = 3.3 \times 10^{-4}\ \text{m}\) or approximately 0.33 mm.

"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