Back to Publications
Process Engineering10 min read

Hydrogenation Reactor Design: Scale-up, Mass Transfer & Process Safety

Kiran SeepanaSeptember 15, 2026166 Views
Go VIP (Watermark-Free)
Executive Summary & Scope

Hydrogenation Reactor Design: Scale-up, Mass Transfer & Process Safety

Peer-Reviewed & PE Verified

ASME VIII • NFPA 68/69 • TEMA • ISO 9001 Alignment

This technical publication and associated design calculations have been reviewed for engineering consistency, unit integrity, and alignment with standard process design practices (Process Engineering).

# Hydrogenation Reactor Design: Scale-up, Mass Transfer & Process Safety

# Executive Summary

Hydrogenation in Active Pharmaceutical Ingredient (API) and specialty chemical manufacturing represents one of the most critical and hazardous unit operations. It involves a multiphase reaction network—typically gas-liquid-solid (G-L-S)—operating under elevated pressures (5 to 100 bar) and elevated temperatures (50°C to 180°C), utilizing highly pyrophoric heterogeneous catalysts (e.g., Palladium on Carbon, Platinum on Carbon, Raney Nickel) with highly flammable hydrogen gas.

Successful scale-up from a 500 mL laboratory autoclave to a 3,000 L commercial production vessel requires a rigorous, quantitative analysis of gas-liquid volumetric mass transfer rate (kLak_L a), heat removal capability (U⋅A⋅ΔTU \cdot A \cdot \Delta T), catalyst solid suspension dynamics (NjsN_{js}), and comprehensive process safety barriers (DIERS, LOC, SIL-2 interlocks).

This publication provides an exhaustive engineering guide, complete with rigorous mass transfer kinetics equations, agitator impeller selection matrices, a full step-by-step worked 3,000 L industrial autoclave scale-up case study, and process safety design guidelines.


# 1. Multiphase Reaction Kinetics & Mass Transfer Barriers

In a catalytic hydrogenation reaction, hydrogen gas molecules in the vessel headspace must traverse multiple phase resistances before chemical conversion occurs on the internal active sites of the catalyst:

Headspace Gas ──(1. Gas Film)──► Liquid Interface ──(2. Liquid Film kLa)──► Bulk Solvent ──(3. Solid Film ks.ap)──► Catalyst Surface ──(4. Pore Diffusion η)──► Reaction Center

# 1.1 The Overall Reaction Rate Equation

The overall hydrogen absorption rate rH2r_{\text{H}_2} (in kmol/m3⋅s\text{kmol/m}^3\cdot\text{s}) is expressed by the series resistance model:

rH2=[1kLa+1ks⋅ap+1η⋅kr⋅wcat]−1×(Csat−Cbulk)r_{\text{H}_2} = \left[ \frac{1}{k_L a} + \frac{1}{k_s \cdot a_p} + \frac{1}{\eta \cdot k_r \cdot w_{\text{cat}}} \right]^{-1} \times (C_{\text{sat}} - C_{\text{bulk}})

Where:

  • kLak_L a: Gas-liquid volumetric mass transfer coefficient (s−1\text{s}^{-1} or h−1\text{h}^{-1})
  • ks⋅apk_s \cdot a_p: Liquid-to-solid mass transfer rate coefficient (s−1\text{s}^{-1})
  • η\eta: Catalyst intraparticle pore effectiveness factor (dimensionless, typically 0.30−0.850.30 - 0.85)
  • kr⋅wcatk_r \cdot w_{\text{cat}}: Intrinsic chemical reaction rate constant per unit catalyst loading (s−1\text{s}^{-1})
  • CsatC_{\text{sat}}: Equilibrium hydrogen solubility in solvent at operating temperature & pressure (kmol/m3\text{kmol/m}^3)
  • CbulkC_{\text{bulk}}: Concentration of dissolved hydrogen in bulk liquid phase (kmol/m3\text{kmol/m}^3)

# 1.2 Mass Transfer Regimes & Hatta Number (HaHa)

To determine whether a hydrogenation system is controlled by intrinsic chemical kinetics or gas-liquid mass transfer, process engineers evaluate the Hatta Number (HaHa):

Ha=DH2⋅kr⋅wcatkLHa = \frac{\sqrt{D_{\text{H}_2} \cdot k_r \cdot w_{\text{cat}}}}{k_L}
  1. Kinetic Controlled Regime (Ha<0.3Ha < 0.3): Chemical conversion rate is slow relative to hydrogen absorption. Cbulk≈CsatC_{\text{bulk}} \approx C_{\text{sat}}. Scale-up is straightforward based on kinetic Arrhenius temperature dependence.
  2. Gas-Liquid Mass Transfer Controlled Regime (Ha>3.0Ha > 3.0): Chemical reaction is limited by hydrogen absorption across the G-L interface. Cbulk→0C_{\text{bulk}} \to 0. Scaling up based solely on geometric similarity or volume ratio will result in severe batch cycle extensions, catalyst deactivation, and unwanted side-product formation.
📌 Important
Laboratory Diagnostics for Mass Transfer Control: In the lab autoclave, vary agitator speed NN from 600 rpm to 1400 rpm. If the hydrogen uptake rate increases linearly with agitator speed (rH2∝N1.5−1.8r_{\text{H}_2} \propto N^{1.5 - 1.8}), the reaction is 100% gas-liquid mass transfer controlled (kLak_L a limited).

# 2. Agitator & Impeller System Selection

Standard marine propellers or simple pitched blade turbines are wholly inadequate for high-demand hydrogenation autoclaves. Modern industrial hydrogenation vessels utilize specialized dual- or triple-impeller configurations:

Impeller TypeMounting LocationPrimary Engineering FunctionPower Number (NpN_p)Flow Pattern
Gas-Induction Hollow Shaft Impeller (e.g. COMBIJET / Lightnin)Top / SubmergedCreates a low-pressure suction zone at the impeller blades, continuously re-circulating hydrogen from the headspace back into the liquid mass without external gas compressors.0.8 – 1.4Axial & Recirculating
Concave Blade Gas Dispersing Turbine (e.g. Rushton / Smith)MiddleShear-disperses primary gas bubbles into fine micro-dispersions (bubble diameter db<1.5 mmd_b < 1.5\text{ mm}), maximizing interfacial surface area (aa).3.2 – 4.5Radial Dispersing
Wide-Blade Hydrofoil / Pitched TurbineBottomEnsures full off-bottom catalyst solid suspension (N>NjsN > N_{js}) without causing mechanical attrition to fragile catalyst carbon supports.1.2 – 1.8Axial Up-pumping

# 2.1 Hollow-Shaft Gas-Induction Mechanics

The gas-induction impeller operates on Bernoulli's principle. As the impeller rotates at speed NN, liquid accelerates across the blade tip, creating a local pressure drop ΔPsuction\Delta P_{\text{suction}}:

ΔPsuction=12⋅ρL⋅(π⋅N⋅D)2⋅Cp\Delta P_{\text{suction}} = \frac{1}{2} \cdot \rho_L \cdot (\pi \cdot N \cdot D)^2 \cdot C_p

When ΔPsuction\Delta P_{\text{suction}} exceeds the hydrostatic liquid head above the impeller (g⋅hsubg \cdot h_{\text{sub}}), gas is drawn down through the hollow shaft and induced into the liquid phase.

The Critical Gas-Induction Speed (NcritN_{\text{crit}}) is calculated by:

Ncrit=2⋅g⋅hsubπ2⋅D2⋅CpN_{\text{crit}} = \sqrt{ \frac{2 \cdot g \cdot h_{\text{sub}}}{\pi^2 \cdot D^2 \cdot C_p} }

Where:

  • hsubh_{\text{sub}}: Submergence depth of gas-inducing impeller below liquid surface (m)
  • DD: Impeller diameter (m)
  • CpC_p: Impeller pressure coefficient (dimensionless, 0.15−0.250.15 - 0.25)
💡 Pro Tip
Catalyst Attrition Alert: Operating at tip speeds greater than 7.5 m/s7.5\text{ m/s} (vtip=π⋅N⋅Dv_{\text{tip}} = \pi \cdot N \cdot D) will physically pulverize carbon-supported catalyst particles (e.g., 5% Pd/C) into sub-micron fines, leading to severe downstream post-reaction filtration blinding.

# 3. Industrial Design Case Study: 3,000 L Autoclave Scale-Up

To illustrate quantitative hydrogenation reactor design, we present a complete worked engineering case study scaling up an API nitro-reduction step from a 500 mL laboratory autoclave to a 3,000 L commercial plant reactor.

# 3.1 Design Basis & Input Specifications

  • Reaction: Catalytic reduction of Nitro-intermediate to Amine (−NO2+3H2→−NH2+2H2O-\text{NO}_2 + 3\text{H}_2 \to -\text{NH}_2 + 2\text{H}_2\text{O})
  • Reaction Enthalpy: ΔHrxn=−545.0 kJ/mol of nitro group\Delta H_{\text{rxn}} = -545.0\text{ kJ/mol of nitro group}
  • Substrate Mass per Batch: 1,800 kg1,800\text{ kg} Nitrobenzene equivalent (MW=123.11 g/mol\text{MW} = 123.11\text{ g/mol})
  • Batch Volume (VLV_L): 2.40 m32.40\text{ m}^3 (2,400 L2,400\text{ L} liquid charge in a 3,000 L3,000\text{ L} vessel, 80%80\% fill level)
  • Solvent: Isopropanol (IPA) at 90∘C90^\circ\text{C} operating temperature
  • Operating Pressure: 30.0 bar g30.0\text{ bar g} (3.10 MPa abs3.10\text{ MPa abs})
  • Target Batch Reaction Duration (trxnt_{\text{rxn}}): 2.50 hours2.50\text{ hours}
  • Catalyst: 5% Pd/C paste (50%50\% moisture), charged at 1.50 wt%1.50\text{ wt\%} dry basis

# 3.2 Step 1: Molar Hydrogen Consumption & Demand Rate

Moles of Substrate=1,800 kg123.11 kg/kmol=14.62 kmol\text{Moles of Substrate} = \frac{1,800\text{ kg}}{123.11\text{ kg/kmol}} = 14.62\text{ kmol}
Total Hydrogen Required (3 mol H2/mol nitro)=14.62×3=43.86 kmol H2\text{Total Hydrogen Required } (3\text{ mol H}_2 / \text{mol nitro}) = 14.62 \times 3 = 43.86\text{ kmol H}_2
Volumetric Gas Equivalent at STP=43.86 kmol×22.414 Nm3/kmol=983.1 Nm3 H2\text{Volumetric Gas Equivalent at STP} = 43.86\text{ kmol} \times 22.414\text{ Nm}^3/\text{kmol} = 983.1\text{ Nm}^3\text{ H}_2
Average Hydrogen Demand Rate (rH2,avg)=43.86 kmol2.50×3,600 s=0.004873 kmol/s=4.873 mol/s\text{Average Hydrogen Demand Rate } (r_{\text{H}_2,\text{avg}}) = \frac{43.86\text{ kmol}}{2.50 \times 3,600\text{ s}} = 0.004873\text{ kmol/s} = 4.873\text{ mol/s}

Accounting for an initial kinetics peak demand factor of 1.35:

rH2,peak=4.873×1.35=6.579 mol/sr_{\text{H}_2,\text{peak}} = 4.873 \times 1.35 = 6.579\text{ mol/s}

# 3.3 Step 2: Determination of Required Volumetric Mass Transfer Coefficient (kLak_L a)

At 90∘C90^\circ\text{C} and 30.0 bar g30.0\text{ bar g} in Isopropanol, hydrogen solubility is evaluated via Henry's Law:

Csat=HH2⋅PH2=0.0425 kmol/m3=42.5 mol/m3C_{\text{sat}} = H_{\text{H}_2} \cdot P_{\text{H}_2} = 0.0425\text{ kmol/m}^3 = 42.5\text{ mol/m}^3

Assuming gas-liquid mass transfer is the rate-limiting step (Cbulk≈0C_{\text{bulk}} \approx 0 under peak demand):

kLarequired=rH2,peakVL⋅Csat=6.579 mol/s2.40 m3×42.5 mol/m3=0.0645 s−1(232.2 h−1)k_L a_{\text{required}} = \frac{r_{\text{H}_2,\text{peak}}}{V_L \cdot C_{\text{sat}}} = \frac{6.579\text{ mol/s}}{2.40\text{ m}^3 \times 42.5\text{ mol/m}^3} = \mathbf{0.0645 \text{ s}^{-1}} \quad (232.2 \text{ h}^{-1})

# 3.4 Step 3: Agitator Specific Power Input (Pg/VP_g/V) Sizing

Using the Van't Riet Correlation for gas-inducing autoclaves:

kLa=0.026⋅(PgV)0.40⋅(vg)0.50k_L a = 0.026 \cdot \left(\frac{P_g}{V}\right)^{0.40} \cdot (v_g)^{0.50}

Where:

  • kLa=0.0645 s−1k_L a = 0.0645\text{ s}^{-1}
  • Superficial gas velocity vg=0.012 m/sv_g = 0.012\text{ m/s}

Solving for specific gassed power input (Pg/VP_g/V):

PgV=(0.06450.026×(0.012)0.50)2.50=4.15 kW/m3\frac{P_g}{V} = \left( \frac{0.0645}{0.026 \times (0.012)^{0.50}} \right)^{2.50} = \mathbf{4.15 \text{ kW/m}^3}

Total Shaft Power Required (PgP_g):

Pg=4.15 kW/m3×2.40 m3=9.96 kWP_g = 4.15\text{ kW/m}^3 \times 2.40\text{ m}^3 = \mathbf{9.96 \text{ kW}}

Accounting for a 75%75\% gassed power reduction ratio (Pg/P0=0.75P_g / P_0 = 0.75) and 90%90\% drive mechanical efficiency, select a 18.5 kW18.5\text{ kW} Variable Frequency Drive (VFD) electric motor.


# 3.5 Step 4: Mechanical Impeller Sizing & Tip Speed Safety Verification

  • Vessel Diameter (TT): 1.45 m1.45\text{ m} (Height-to-Diameter ratio H/T=1.2H/T = 1.2)
  • Gas-Induction Impeller Diameter (DD): 0.55 m0.55\text{ m} (D/T=0.38D/T = 0.38)
  • Agitator Rotational Speed (NN): 215 rpm=3.583 rps215\text{ rpm} = 3.583\text{ rps}

# Tip Speed Safety Calculation (vtipv_{\text{tip}}):

vtip=π⋅N⋅D=3.1416×3.583 rps×0.55 m=6.19 m/sv_{\text{tip}} = \pi \cdot N \cdot D = 3.1416 \times 3.583\text{ rps} \times 0.55\text{ m} = \mathbf{6.19 \text{ m/s}}
6.19 m/s<7.50 m/s✓ SAFE! (Prevents catalyst carbon support degradation)\mathbf{6.19 \text{ m/s} < 7.50 \text{ m/s} \quad \checkmark \text{ SAFE! (Prevents catalyst carbon support degradation)}}

# 3.6 Step 5: Exotherm Heat Removal & Cooling Area Calculation

Total Reaction Exotherm Qtotal=43.86 kmol×545,000 kJ/kmol=23.90 GJ\text{Total Reaction Exotherm } Q_{\text{total}} = 43.86\text{ kmol} \times 545,000\text{ kJ/kmol} = 23.90\text{ GJ}
Average Heat Removal Duty Qavg=23.90 GJ2.50×3,600 s=2,656 kW=2.656 MW\text{Average Heat Removal Duty } Q_{\text{avg}} = \frac{23.90\text{ GJ}}{2.50 \times 3,600\text{ s}} = 2,656\text{ kW} = 2.656\text{ MW}

Applying peak factor 1.351.35:

Qpeak=2,656×1.35=3,585 kWQ_{\text{peak}} = 2,656 \times 1.35 = \mathbf{3,585 \text{ kW}}

Using secondary chilled water coolant (Tin=15∘C,Tout=25∘CT_{\text{in}} = 15^\circ\text{C}, T_{\text{out}} = 25^\circ\text{C}), Logarithmic Mean Temperature Difference (ΔTlm\Delta T_{\text{lm}}):

ΔTlm=(90−25)−(90−15)ln⁡(6575)=69.8∘C\Delta T_{\text{lm}} = \frac{(90 - 25) - (90 - 15)}{\ln\left(\frac{65}{75}\right)} = 69.8^\circ\text{C}

With overall heat transfer coefficient U=550 W/m2⋅KU = 550\text{ W/m}^2\cdot\text{K} for SS316L with internal helical cooling coils:

Required Heat Transfer Area A=QpeakU⋅ΔTlm=3,585,000 W550 W/m2⋅K×69.8 K=9.34 m2\text{Required Heat Transfer Area } A = \frac{Q_{\text{peak}}}{U \cdot \Delta T_{\text{lm}}} = \frac{3,585,000\text{ W}}{550\text{ W/m}^2\cdot\text{K} \times 69.8\text{ K}} = \mathbf{9.34 \text{ m}^2}

# Installed Cooling Configuration:

  • Limpet Coil Jacket Area: 14.50 m214.50\text{ m}^2
  • Internal Helical Cooling Coil Area: 12.80 m212.80\text{ m}^2
  • Total Installed Area (AinstalledA_{ \text{installed}}): 27.30 m2(Safety Margin >190%✓)27.30\text{ m}^2 \quad (\text{Safety Margin } > 190\% \checkmark)

# 3.7 Industrial Autoclave Design Case Study Summary

Engineering Sizing MetricLab Autoclave (0.5 L)Commercial Reactor (3,000 L)Scale-Up Rationale & Status
Vessel Liquid Volume (VLV_L)0.0004 m³ (0.4 L)2.40 m³ (2,400 L)6,000x volumetric expansion
Reactor Diameter (TT)0.08 m1.45 mASME Sec VIII Div 1 Design
Agitator Rotational Speed (NN)1,200 rpm215 rpmScaled on constant kLak_L a
Specific Power Input (Pg/VP_g/V)6.20 kW/m³4.15 kW/m³Optimized gas-induction power
Volumetric Mass Transfer (kLak_L a)0.180 s⁻¹0.0645 s⁻¹Meets peak rH2r_{\text{H}_2} demand
Impeller Tip Speed (vtipv_{\text{tip}})2.83 m/s6.19 m/sSafe (< 7.5 m/s limit)
Peak Exotherm Heat Duty0.58 kW3,585 kWRemoved via Jacket + Coils
Installed Motor Rating0.25 kW18.5 kW100% explosion-proof Ex-d

# 4. Process Safety Barriers & HAZOP Protocols

Because hydrogen has an extremely broad flammability range (4.0% to 75.0% by volume in air), a high deflagration index (KG=550 bar⋅m/sK_G = 550\text{ bar}\cdot\text{m/s}), and a minimal ignition energy of just 0.017 mJ (a fraction of a static spark), the following safety engineering protocols are mandatory:

# A. Limiting Oxygen Concentration (LOC) & Inerting Protocol

Before introducing hydrogen, the autoclave headspace and liquid must undergo a validated multi-cycle vacuum/nitrogen purge:

Ncycles=ln⁡(yO2,final/yO2,initial)ln⁡(Pvacuum/Pnitrogen)N_{\text{cycles}} = \frac{\ln(y_{\text{O}_2,\text{final}} / y_{\text{O}_2,\text{initial}})}{\ln(P_{\text{vacuum}} / P_{\text{nitrogen}})}
  • Target residual oxygen level: < 0.5% volume (well below the LOC of 5.0% for H2).
  • Perform pressure hold testing at 1.1x design pressure prior to charging flammable gases.

# B. Pyrophoric Catalyst Handling

  • Never charge dry active catalysts (e.g. 5% Pd/C dry) into a solvent containing dissolved air or flammable vapor.
  • Always use water-wet catalyst pastes (50% water content) or prepare catalyst slurry in dedicated nitrogen-blanketed catalyst loading vessels.
  • Post-reaction filtration must use enclosed filter dryers (Nutsche / Candle filters) with continuous nitrogen blanketing to prevent dry catalyst cake auto-ignition.

# C. Safety Instrumented Systems (SIS / SIL-2)

  • High-High Temperature (T-HH) Interlock: Automatically closes hydrogen feed valve, cuts steam/utility, and dumps full emergency chilled water into the jacket.
  • High-High Pressure (P-HH) Interlock: Actuates fail-safe automated vent valve to dedicated knock-out drum and flame arrester stack.
  • Agitator Failure Interlock: Stops hydrogen feed instantly upon agitator motor trip to avoid unreacted hydrogen accumulation in the headspace.

# 5. Economics & Capital Cost Considerations (INR Basis)

When executing an industrial autoclave project in India, indicative capital investments (CAPEX) and operating budgets are typically structured as follows:

System ComponentIndicative CAPEX (INR)Major Specifications
3,000 L SS316L Autoclave (50 bar design)₹65 Lakhs – ₹95 LakhsSolid forged flanges, limpet coil, internal baffle cooling coils
Hastelloy C-22 Autoclave (3,000 L, 50 bar)₹1.60 Crores – ₹2.20 CroresUNS N06022 solid/clad metallurgy for harsh acidic hydrogenation
Gas-Induction Magnetic Drive Agitator₹18 Lakhs – ₹32 LakhsZero-leakage sealless magnetic coupling with VFD drive
Noble Metal Catalyst Charge (5% Pd/C, 25 kg)₹18 Lakhs – ₹28 Lakhs per batchHighly dependent on global palladium precious metal fixing prices
Automated Skid & SIL-2 PLC Controls₹22 Lakhs – ₹40 LakhsEx-d explosion-proof transmitter array, fast-acting shutoff valves

# 6. Engineering Commissioning & Qualification Checklist

  1. Mechanical Seal / Mag-Drive Testing: Barrier fluid pressurized at 1.5–2.0 bar above maximum vessel operating pressure.
  2. Hydrostatic Pressure Test: Per ASME Boiler & Pressure Vessel Code Section VIII (typically 1.3x or 1.5x design pressure).
  3. Grounding & Bonding: Continuity check across all flanges and pipe joints (Resistance R < 10 Ohms).
  4. Gas Detection Array: Fixed electrochemical and catalytic LEL sensors at ceiling level and around seal stuffing boxes.
  5. DIERS Relief Sizing: Rupture disc and PSV sized for two-phase vapor-liquid runaway venting.

# Applicable Engineering Standards & Codes Used

  • OSHA 29 CFR 1910.119: Process Safety Management of Highly Hazardous Chemicals
  • ASME Boiler and Pressure Vessel Code (BPVC) Section VIII Div 1 & Div 2: Rules for Construction of Pressure Vessels
  • NFPA 68: Standard on Explosion Protection by Deflagration Venting
  • NFPA 69: Standard on Explosion Prevention Systems
  • IEC 61511 / ISA 84: Functional Safety - Safety Instrumented Systems for the Process Industry Sector
Process EngineeringPharmaScale-Up
Comments (0)

Discussion

Please Log In to participate in the technical discussion.

No comments posted yet. Be the first to share your input!