# Industrial Solvent Recovery & Azeotrope Breaking: Packed Column Design and Membrane Separation in Multipurpose API Plants

# Thermodynamic Modeling, Vapor-Liquid Equilibrium (VLE), Structured Packing Hydrodynamics, and Reboiler Safety in Commercial Pharmaceutical Distillation

# Executive Summary

In Active Pharmaceutical Ingredient (API) and fine chemical synthesis, raw materials do not dictate the volumetric bulk of plant operations—solvents do. Between extraction, crystallization, chromatography, and vessel cleanings, organic solvents account for $80\text{ to }90\%$ of the total cumulative mass deployed in an API facility.

Disposing of spent mother liquors via high-temperature thermal oxidation (hazardous incineration) is economically crippling and environmentally unsustainable:

• Cost Burden: Incineration costs range from \300 to \800 per metric ton, while virgin pharma-grade solvents continue to escalate in price.
• Scope 1 & 3 Carbon Emissions: Every ton of incinerated solvent releases approximately $2.2\text{–}2.8\text{ tons of } CO_2$.

Recovering and purifying spent solvents to cGMP pharmaceutical specifications ($>99.5\text{ wt\%}$ purity, $<0.05\text{ wt\%}$ water) cuts API manufacturing costs by $20\text{ to }35\%$. However, pharmaceutical solvent mixtures rarely behave ideally. They frequently form complex minimum-boiling azeotropes, liquid-liquid immiscibility gaps, and thermally sensitive reboiler residues.

This engineering guide details the thermodynamic equations of state, column hydraulic sizing protocols (HETP, F-factor), modern azeotrope-breaking technologies (Pressure-Swing Distillation and Pervaporation), and a fully worked commercial scale-up case study.

# 1. Vapor-Liquid Equilibrium (VLE) & Non-Ideal Activity Coefficients

For an ideal vapor-liquid mixture, Raoult's law applies. In pharmaceutical solvents (such as Isopropanol/Water, Ethyl Acetate/Water, Toluene/Methanol), strong intermolecular forces (hydrogen bonding, dipole interactions) create severe non-ideality:

Where:

• $P$: Operating column pressure ($\text{kPa}$)
• $y_i, x_i$: Vapor and liquid mole fractions of component $i$
• $P_i^{sat}(T)$: Pure component vapor pressure at temperature $T$ (via Antoine equation)
• $\phi_i$: Vapor phase fugacity coefficient ($\approx 1.0$ at low to moderate pressures)
• $\gamma_i$: Liquid phase activity coefficient

  y_IPA (Vapor Mole Fraction)
  1.0 ┌──────────────────────────────────────────┐
      │                                   • Azeotrope (x = y = 0.68)
      │                           • • • •        │
      │                    • • •                 ▼
      │              • •                         No Separation Possible!
      │          • •                             (Relative Volatility α = 1.0)
      │      • •
      │   • •         Ideal Raoult Curve
      │ •       -----------------
  0.0 └──────────────────────────────────────────┘
      0.0                                      1.0
                 x_IPA (Liquid Mole Fraction)

# 1.1 The Azeotropic Condition

At the azeotropic composition ($x_i = y_i$), the relative volatility ($\alpha_{ij}$) collapses to unity:

At this point, conventional fractional distillation reaches an absolute thermodynamic pinch: no further enrichment is possible regardless of column height or reflux ratio.

# 2. Common Pharmaceutical Azeotropes & Separation Strategies

# 3. Advanced Azeotrope Breaking Technologies

                  Option A: Pressure-Swing Distillation (PSD)
┌────────────────────────┐                    ┌────────────────────────┐
│ Low-Pressure Column    │ Overhead Azeotrope │ High-Pressure Column   │ Pure Bottoms
│ (e.g., 0.2 bar Vacuum) │───────────────────►│ (e.g., 6.0 bar Gauge)  │ Component B
└────────────────────────┘                    └────────────────────────┘
                                                         │ Recycled Azeotrope
                                                         ▼ (Back to Column 1)

                  Option B: Hybrid Distillation-Pervaporation
┌────────────────────────┐ Overhead Vapor     ┌────────────────────────┐
│ Fractional Distillation│ (88 wt% IPA / H2O) │ Zeolitic Pervaporation │ Pure Anhydrous IPA
│ Column                 │───────────────────►│ Membrane Skid (NaA)    │ (99.8 wt%)
└────────────────────────┘                    └────────────────────────┘
                                                         │
                                                         ▼ Permeate: Pure Water

# 3.1 Pressure-Swing Distillation (PSD)

Pressure-swing distillation exploits the fact that azeotropic composition is pressure-dependent. By operating two columns in series at different pressures (e.g., Column 1 at $0.2\text{ bar}$ vacuum and Column 2 at $6\text{ bar}$ pressure), the feed to the second column lies outside its new azeotropic envelope, enabling pure product recovery from the bottoms of both towers.

# 3.2 Hydrophilic Pervaporation Membranes

Polymeric or inorganic Zeolite NaA (Linde Type A) ceramic membranes feature molecular pore diameters of $3.8\text{ Å}$.

• Water molecules ($d_k \approx 2.8\text{ Å}$) pass rapidly through the crystalline lattice under vacuum suction.
• Larger organic solvent molecules (IPA: $d_k \approx 4.7\text{ Å}$, Ethanol: $d_k \approx 4.3\text{ Å}$) are completely rejected.
• Energy Advantage: Pervaporation requires energy only to vaporize the small water permeate fraction ($10\text{–}15\%$), rather than boiling the entire solvent mass multiple times.

# 3.3 Solution-Diffusion Kinetics in Zeolite NaA Membranes

Water transport across a hydrophilic Zeolite NaA molecular sieve membrane is modeled by the solution-diffusion mechanism:

$J_w = \frac{P_{w,mem}}{L_{mem}} \cdot \left( \gamma_{w,feed} \cdot x_{w,feed} \cdot P_w^{sat}(T) - p_{w,permeate} \right)$

Where:

• $J_w$: Water permeate flux ($\text{kg/m}^2\cdot\text{h}$)
• $P_{w,mem} / L_{mem}$: Membrane hydraulic permeance ($\text{kg/m}^2\cdot\text{h}\cdot\text{kPa}$)
• $p_{w,permeate}$: Vacuum permeate pressure maintained downstream (typically $\le 15\text{ mbar}$)

Because the crystalline cages of Zeolite NaA ($3.8\text{ \AA}$) are strictly smaller than the kinetic diameter of Isopropanol ($4.7\text{ \AA}$), the membrane separation factor ($\alpha_{sep} = (y_w/y_{IPA}) / (x_w/x_{IPA})$) exceeds $10,000$, yielding ultra-pure permeate ($>99.5\text{ wt\% Water}$) and leaving pure anhydrous IPA behind.

# 4. Packed Column Hydraulic Design: Sizing Structured Packing

In multipurpose pharmaceutical plants, corrugated sheet metal structured packing (e.g., Sulzer Mellapak 250Y or 500X) has almost entirely replaced bubble-cap and sieve trays due to its low pressure drop per theoretical stage ($\Delta P / \text{HETP} \le 0.5\text{ mbar}$), which prevents thermal decomposition of sensitive solvent residues.

# 4.1 F-Factor and Column Diameter Sizing

The column diameter ($D_{col}$) is determined by the gas capacity F-factor:

Where:

• $u_v$: Superficial vapor velocity ($u_v = V / (\pi D_{col}^2 / 4)$, in $\text{m/s}$)
• $\rho_v$: Vapor density ($\text{kg/m}^3$)

At the flooding point ($F_{flood}$), liquid downflow is held up by ascending vapor drag. Distillation columns are safely operated at $70\text{–}80\%$ of flood velocity:

# 4.2 Height Equivalent to a Theoretical Plate (HETP)

The total packed height ($Z_{pack}$) required to achieve separation is:

• Structured Packing (250 $\text{m}^2/\text{m}^3$): $\text{HETP} \approx 0.35\text{ to }0.45\text{ m}$.
• Structured Packing (500 $\text{m}^2/\text{m}^3$): $\text{HETP} \approx 0.20\text{ to }0.30\text{ m}$.
• Random Packing (1-inch SS Pall Rings): $\text{HETP} \approx 0.55\text{ to }0.75\text{ m}$.

# 4.3 Fenske-Underwood-Gilliland (FUG) Shortcut Sizing Equations

Before executing rigorous multi-stage ASPEN simulations, process engineers utilize the classical FUG shortcut equations to establish baseline column parameters:

1. Minimum Theoretical Stages ($N_{min}$) via Fenske Equation:
   $N_{min} = \frac{\ln\left[ \left( \frac{x_{D,LK}}{x_{D,HK}} \right) \cdot \left( \frac{x_{B,HK}}{x_{B,LK}} \right) \right]}{\ln(\alpha_{avg})}$

2. Minimum Reflux Ratio ($R_{min}$) via Underwood Equations:
   $\sum_{i=1}^n \frac{\alpha_i \cdot x_{F,i}}{\alpha_i - \theta} = 1 - q$
   $R_{min} = \sum_{i=1}^n \frac{\alpha_i \cdot x_{D,i}}{\alpha_i - \theta} - 1$
   Where $\theta$ is the Underwood root lying between $\alpha_{LK}$ and $\alpha_{HK}$, and $q$ is the thermal feed condition ($q=1.0$ for saturated liquid).

# 5. Reboiler Safety: Peroxides, Foaming, and Thermal Runaway

Solvent distillation reboilers are the site of numerous catastrophic industrial accidents:

[ Reboiler Sump: Long Thermal Residence Time ]
                       │
      ┌────────────────┴────────────────┐
      ▼                                 ▼
Peroxide Concentration (Ethers)    Residue Thermal Decomposition
• THF, 2-MeTHF, Dioxane            • Reactive intermediates / catalysts
• Boiling to dryness concentrates  • Auto-catalytic runaway exotherm
  explosive hydroperoxides!        • Sudden overpressurization!

# 6. Worked Industrial Case Study: Designing a 1,200 kg/h IPA Recovery Skid

# 6.1 Process Specifications

• Feed Stream: Spent crystallization mother liquor: $1,200\text{ kg/h}$ containing $80\text{ wt\%}$ Isopropanol (IPA) and $20\text{ wt\%}$ Water.
• Target Products:
  1. Anhydrous IPA: $\ge 99.8\text{ wt\%}$ purity (maximum $0.1\text{ wt\%}$ water).
  2. Stripped Aqueous Effluent: $\le 0.5\text{ wt\%}$ IPA (suitable for biological ETP).
• System Design Selected: Fractional Packed Column (Atmospheric) producing overhead azeotropic vapor ($87.7\text{ wt\%}$ IPA), coupled directly to an inline Zeolite NaA Pervaporation Membrane Skid.

       Feed (1,200 kg/h: 80% IPA / 20% H2O)
                   │
                   ▼
     ┌───────────────────────────┐
     │ Packed Distillation Tower │ ──► Overhead Azeotrope (1,095 kg/h: 87.7% IPA)
     │ (SS316L, Mellapak 250Y)   │                │
     └─────────────┬─────────────┘                ▼
                   │ Bottoms              ┌──────────────────────────┐
                   ▼ (105 kg/h Water)     │ NaA Pervaporation Skid   │ ──► Permeate: Water (135 kg/h)
                                          └────────────┬─────────────┘
                                                       ▼ Retentate Product
                                            Anhydrous IPA: 960 kg/h (>99.8 wt%)

# 6.2 Distillation Column Sizing Calculations

# Step 1: Mass Balance

• Total IPA in Feed: $1,200\text{ kg/h} \times 0.80 = 960\text{ kg/h}$.
• Overhead Azeotrope Rate ($D$):

• Water Evaporated with Overhead: $1,094.6 - 960 = 134.6\text{ kg/h}$.
• Bottoms Effluent Rate ($B$):

# Step 2: Minimum Reflux & Stage Requirements

Using ASPEN Plus VLE simulation with NRTL-RK fluid package:

• Minimum Reflux Ratio: $R_{min} = 1.45$.
• Operating Reflux Ratio: $R = 1.3 \times R_{min} \approx \mathbf{1.90}$.
• Required Theoretical Stages: $N_{theor} = 18$ stages (including reboiler).

# Step 3: Column Internal Diameter ($D_{col}$)

At the column top ($P = 1.013\text{ bar}$, $T_{overhead} = 80.3^\circ\text{C}$):

• Total vapor traffic: $V = D \times (1 + R) = 1,094.6\text{ kg/h} \times (1 + 1.90) = 3,174.3\text{ kg/h} = 0.8818\text{ kg/s}$.
• Vapor density: $\rho_v \approx 2.05\text{ kg/m}^3$.
• Volumetric vapor flow rate:

• Structured Packing (Mellapak 250Y): $F_{flood} \approx 2.4\text{ Pa}^{0.5}$.
• Operating at $75\%$ of flood:

• Allowable vapor velocity:

• Column Cross-Sectional Area ($A_{col}$):

• Internal Column Diameter:

# Step 4: Packed Height Sizing

• Selecting Mellapak 250Y ($\text{HETP} = 0.40\text{ m}$):

• Divided into two $3.6\text{ m}$ packed beds with an intermediate liquid redistributor and feed tray.

# Step 5: Pervaporation Membrane Sizing

The overhead vapor from the column ($1,095\text{ kg/h}$ at $87.7\text{ wt\%}$ IPA) passes into the Zeolite NaA membrane module:

• Water to be removed: $134.6\text{ kg/h} - 1.0\text{ kg/h (residual)} = 133.6\text{ kg/h}$.
• High-performance NaA ceramic membrane water flux: $J_{water} \approx 3.2\text{ kg/m}^2\cdot\text{h}$ at $105^\circ\text{C}$ and $15\text{ mbar}$ vacuum permeate pressure.
• Required Membrane Surface Area ($A_{mem}$):

# 7. Economic Payback & Utility Savings