Back to Publications
Process Engineering31 min read

Batch Reactor Agitator Design: Impeller Selection, Power, Mixing, Scale-Up, and Engineering Correlations

Kiran SeepanaAugust 15, 202623 Views
Executive Summary & Scope

A definitive, calculation-oriented chemical engineering guide on batch reactor agitator design, covering hydrodynamics, impeller selection, power number, pumping capacity, Zwietering Njs solid suspension, heat transfer, shaft dynamics, and pilot-to-commercial scale-up.

# Batch Reactor Agitator Design: Impeller Selection, Power, Mixing, Scale-Up, and Engineering Correlations

In chemical synthesis and Active Pharmaceutical Ingredient (API) manufacturing, a batch reactor agitator is far more than a simple electric motor driving a rotating shaft.

It is the primary transport phenomena engine that dictates the rate of species blending, turbulent energy dissipation, boundary-layer heat removal, solid-liquid dissolution, and gas-liquid interfacial mass transfer across the entire vessel.

A poorly designed or improperly scaled agitation system does not merely cause extended batch cycle times—it leads to uncontrolled temperature hot-spots, localized selectivity degradation, secondary byproduct formation, severe crystal attrition, un-suspended catalyst dead zones, and catastrophic mechanical shaft failure.

The central engineering question this technical guide addresses is:

How do process, project, and equipment design engineers select, calculate, size, and scale an agitator for a batch reactor so that required mixing performance is achieved seamlessly from laboratory (1 L) to pilot (100 L) and commercial scale (5 KL to 10 KL)?

Batch Reactor Agitator System Architecture & Scale-Up Blueprint
Batch Reactor Agitator System Architecture & Scale-Up Blueprint


# 1. Introduction: Why Agitator Design Matters

In an industrial batch reactor, an agitator must accomplish multiple concurrent physical and chemical objectives across different reaction steps:

  • Homogeneous Blending: Rapidly dispersing miscible liquid feeds to avoid localized stoichiometry excesses.
  • Thermal Uniformity & Heat Transfer: Sweeping liquid past jacketed vessel walls and internal baffles to maximize process-side heat transfer coefficients (hih_i).
  • Solids Suspension & Dissolution: Maintaining heterogeneous catalysts (e.g., 5% Pd/C) or raw material crystals in a fully off-bottom suspended state.
  • Gas-Liquid Dispersion: Shearing sparged gases into micro-bubbles and recirculating the headspace gas to maximize interfacial area (aa) and overall volumetric mass transfer (kLak_L a).
  • Liquid-Liquid Emulsification: Shearing immiscible organic-aqueous phases into stable droplets during phase-transfer catalysis or extractive workups.
  • Controlled Crystallization: Providing gentle bulk circulation while keeping maximum shear rates low enough to prevent crystal attrition and uncontrolled secondary nucleation.

An impeller assembly perfectly optimized for low-viscosity solvent blending (e.g., high-speed hydrofoil) will fail completely when tasked with circulating a pseudoplastic, shear-thinning polymer slurry or dispersing high-flow hydrogen gas. Agitator design must always begin with the fundamental mixing process objective.


# 2. What Exactly Does an Agitator Need to Achieve?

Before specifying motor kilowatts or selecting impeller geometries, the process engineer must classify the dominant hydrodynamic duty:

Mixing ObjectivePrimary Fluid Dynamic RequirementKey Governing MetricCritical Engineering Hazard of Undersizing
Miscible Liquid BlendingHigh bulk volumetric pumping capacity (QQ).Blend Time (tmix<30 st_{mix} < 30\text{ s})Localized pH/concentration spikes, unreacted reagent pockets.
Heat-Transfer EnhancementHigh fluid shear and boundary layer disruption at vessel wall.Nusselt Number (NuNu), Inside Film Coefficient (hih_i)Thermal runaway, wall fouling, severe batch cooling delays.
Solid-Liquid SuspensionAxial off-bottom upward velocity exceeding particle settling velocity (utu_t).Just-Suspended Speed (NjsN_{js} via Zwietering)Catalyst settling, incomplete conversion, localized hot spots.
Gas-Liquid Mass TransferHigh turbulent shear near sparger; gas bubble breakup.Volumetric Mass Transfer Coefficient (kLak_L a), Gassed Power (PgP_g)Gas flooding, severe reaction rate starvation (Ha > 3).
Liquid-Liquid DispersionMicro-scale turbulent eddy dissipation (ϵavg\epsilon_{avg}).Sauter Mean Droplet Diameter (d32d_{32}), Weber Number (WeWe)Phase separation, low interfacial reaction rate.
Controlled CrystallizationUniform bulk suspension with minimal maximum tip shear (γ˙max\dot{\gamma}_{max}).Tip Speed (vtip<2.2 m/sv_{tip} < 2.2\text{ m/s}), Power/Volume (P/VP/V)Crystal breakage, fines generation, clogged filter dryers (ANFD).

# 3. Basic Components of a Reactor Agitation System

A standard top-entry batch reactor agitation system comprises seven integrated mechanical components:

  1. Electric Drive Motor: Standard 4-pole or 6-pole TEFC / Flameproof (Ex-d) induction motor, operated through a Variable Frequency Drive (VFD) for precision speed modulation.
  2. Speed Reduction Gearbox: High-efficiency helical or bevel-helical gearbox providing fixed torque multiplication (typical ratios 1:81:8 to 1:251:25).
  3. Flexible / Rigid Coupling: Connects the gearbox output shaft to the agitator drive shaft, absorbing minor angular and parallel misalignments.
  4. Mechanical Seal Assembly: Double-acting, pressurized liquid-lubricated or dry-gas mechanical seal with a barrier fluid reservoir (operating at 1.52.0 bar1.5\text{--}2.0\text{ bar} above maximum reactor operating pressure) to ensure hermetic containment.
  5. Agitator Shaft: Solid or heavy-wall tubular stainless steel (SS316L, Hastelloy C-22, Alloy 59) cantilevered shaft designed for combined torsional shear and cyclic bending stresses.
  6. Impeller Stages: Single or multi-stage hubs with bolted or welded blades positioned at calibrated liquid heights.
  7. Vessel Baffles: Fully welded or removable vertical wall baffles that convert rotational liquid swirl into vertical axial circulation loops.

# Mounting Configurations:

  • Top-Entry Agitator (Dominant in Pharma/API): Cantilevered shaft entering through the top dome nozzle. Allows easy maintenance, clean-in-place (CIP), and eliminates leak paths below the liquid level.
  • Bottom-Entry Agitator: Common in bioprocessing fermenters and sterile magnetic-drive mixing tanks. Requires low headroom but presents seal reliability risks in abrasive chemical slurries.
  • Side-Entry Agitator: Used almost exclusively in large petroleum and bulk storage tanks; unsuited for batch API synthesis reactors.

# 4. Important Reactor Geometry Parameters

Geometric similarity forms the mathematical foundation of chemical agitation and scale-up. The primary dimensions of an agitated vertical cylindrical vessel with a standard dished bottom are:

  • TT = Vessel Inside Diameter (m)
  • HH = Liquid Charge Height from bottom tangent line (m)
  • DD = Impeller Swept Outer Diameter (m)
  • CC = Impeller Bottom Clearance from lowest point of dish (m)
  • BB = Radial Baffle Width (m)
  • SS = Inter-Impeller Axial Spacing in multi-tier assemblies (m)
  • dsd_s = Solid Agitator Shaft Diameter (m)

# Standard Geometric Starting Ratios:

Dimension RatioStandard RangeRecommended BaselinePhysical Significance
D/TD / T (Impeller-to-Tank Ratio)0.250.500.25 - 0.500.330.33 (Axial) / 0.350.35 (Radial)Balances shear rate vs. bulk volumetric circulation flow.
H/TH / T (Aspect Ratio)0.801.500.80 - 1.501.01.21.0 - 1.2Standard dished batch reactor envelope; H/T>1.2H/T > 1.2 requires dual impellers.
C/TC / T (Bottom Clearance Ratio)0.150.350.15 - 0.350.250.330.25 - 0.33Optimum off-bottom solid suspension and dish sweeping.
B/TB / T (Baffle Width Ratio)0.080.120.08 - 0.120.100.10 (T/10T/10)Standard 4-baffle fully baffled condition; suppresses vortexing.
Bgap/TB_{gap} / T (Baffle-to-Wall Gap)0.0150.0250.015 - 0.0250.020.02 (T/50T/50)Prevents solid particle stagnation and dead-zone accumulation behind baffles.
S/DS / D (Dual Impeller Spacing)1.01.51.0 - 1.51.21.31.2 - 1.3Prevents hydrodynamic flow interference between adjacent stages.

# 5. Impeller Types: Fluid Mechanics & Selection Matrix

Impellers are classified by their primary discharge vector into Axial-Flow, Radial-Flow, and Close-Clearance (High-Viscosity) designs:

Impeller FamilySpecific GeometryPrimary Flow VectorPower Number (NpN_p)Pumping Number (NQN_Q)Operating Viscosity Range (mPa·s)Primary Industrial Application
High-Efficiency HydrofoilHE-3, A310, MaxfloPure Axial (Downward)0.280.400.28 - 0.400.500.600.50 - 0.60<2,000< 2,000Low-shear blending, rapid circulation, low-power solids suspension.
Pitched-Blade Turbine (PBT)4-Blade 45° PitchMixed (Axial + Radial)1.201.701.20 - 1.700.750.850.75 - 0.85<10,000< 10,000General synthesis, solid-liquid dissolution, crystallization, heat transfer.
Marine Propeller3-Blade True PitchPure Axial0.350.500.35 - 0.500.450.550.45 - 0.55<1,000< 1,000Small pilot vessels, high-speed portable mixers, bench autoclaves.
Rushton Disk Turbine (RDT)6 Flat Blades on DiskPure Radial4.805.504.80 - 5.500.700.800.70 - 0.80<5,000< 5,000Gas-liquid dispersion, hydrogenation, liquid-liquid high-shear emulsification.
Curved-Blade RadialSmith (CD-6), ConcavePure Radial2.803.202.80 - 3.200.800.900.80 - 0.90<5,000< 5,000High-rate gas dispersion; resists gas cavity flooding better than Rushton.
Anchor AgitatorContour Wall ScraperTangential / Laminar0.501.200.50 - 1.20<0.15< 0.155,00050,0005,000 - 50,000High-viscosity resin synthesis, cooling crystallization wall heat transfer.
Helical RibbonSingle/Double FlightTop-to-Bottom Axial1.503.501.50 - 3.500.300.500.30 - 0.5020,000500,000+20,000 - 500,000+Polymers, creams, high-density pseudoplastic pastes in laminar regime.

# 6. Axial vs. Radial Flow Hydrodynamics

Understanding the physical flow field generated by an impeller is essential for selecting the correct geometry:

# A. Axial-Flow Impellers (PBT, Hydrofoil, Propeller)

  • Fluid Motion: Fluid is propelled parallel to the agitator shaft toward the bottom dish, sweeps across the lower vessel head, travels upward along the outer baffled wall, and returns to the impeller eye.
  • Hydrodynamic Characteristic: Generates maximum volumetric liquid turnover (QQ) per unit of motor power consumed (PP).
  • Primary Utility: Blending of miscible fluids, solids suspension, and general heat-transfer circulation.

# B. Radial-Flow Impellers (Rushton Turbine, Flat Blade)

  • Fluid Motion: Fluid is discharged radially outward perpendicular to the shaft toward the vessel wall, where it splits into two distinct circulation loops: one moving upward toward the surface, and one moving downward toward the dish.
  • Hydrodynamic Characteristic: Generates high localized shear rates and turbulent energy dissipation (ϵ\epsilon) in the impeller discharge stream, but exhibits lower bulk circulation efficiency.
  • Primary Utility: Breaking up sparged gas bubbles, dispersing immiscible liquid droplets, and high-intensity gas-liquid mass transfer.

# 7. Reynolds Number for Agitated Vessels

Fluid flow inside an agitated reactor is governed by the Impeller Reynolds Number (ReRe):

Re=ρND2μRe = \frac{\rho \cdot N \cdot D^2}{\mu}

Where:

  • ReRe = Impeller Reynolds Number (Dimensionless)
  • ρ\rho = Process fluid density (kg/m3\text{kg/m}^3)
  • NN = Impeller rotational speed in revolutions per second (rev/s\text{rev/s} or s1\text{s}^{-1}) [Note: N=RPM/60N = \text{RPM} / 60]
  • DD = Impeller outer diameter (m)
  • μ\mu = Dynamic fluid viscosity (Pas\text{Pa}\cdot\text{s} or kg/(ms)\text{kg}/(\text{m}\cdot\text{s})) [Note: 1 cP=103 Pas1\text{ cP} = 10^{-3}\text{ Pa}\cdot\text{s}]

# Hydrodynamic Flow Regimes in Agitated Vessels:

  1. Laminar Regime (Re<10Re < 10): Viscous forces dominate. Fluid motion occurs solely along streamlines; zero turbulence. Power consumption is directly proportional to viscosity and independent of density (PμN2D3P \propto \mu N^2 D^3).
  2. Transitional Regime (10Re10,00010 \le Re \le 10,000): Both inertial and viscous forces interact. Complex vortex structures form near blade tips while bulk fluid exhibits laminar decay.
  3. Fully Turbulent Regime (Re>10,000Re > 10,000): Inertial forces completely dominate. The boundary layer is fully turbulent; Power Number (NpN_p) reaches an asymptotic constant value independent of fluid viscosity (PρN3D5P \propto \rho N^3 D^5).

# 8. The Impeller Power Number (NpN_p)

The Power Number (NpN_p) (also known as the Newton Number) is the fundamental dimensionless drag coefficient representing an impeller's resistance to rotational fluid motion:

Np=PρN3D5N_p = \frac{P}{\rho \cdot N^3 \cdot D^5}

Rearranging to solve for required shaft power demand (PP):

P=NpρN3D5P = N_p \cdot \rho \cdot N^3 \cdot D^5

Where:

  • PP = Net shaft power delivered to fluid (Watts, W\text{W} or J/s\text{J/s})
  • NpN_p = Turbulent Impeller Power Number (Dimensionless)
  • ρ\rho = Liquid density (kg/m3\text{kg/m}^3)
  • NN = Agitator rotational speed (rev/s\text{rev/s})
  • DD = Impeller outer diameter (m)

Critical Scaling Principle: In the turbulent regime, power escalates with the 3rd power of speed (N3N^3) and the 5th power of diameter (D5D^5)! A mere 20%20\% increase in impeller diameter at constant RPM increases required power by (1.20)5=2.488(1.20)^5 = 2.488 (+149%+149\% increase)!


# 9. The Impeller Power Curve (NpN_p vs. ReRe)

Agitator equipment manufacturers provide empirical NpN_p vs. ReRe Power Curves calibrated for specific impeller blade profiles, hub angles, and vessel baffle configurations:

  • Laminar Zone (Re<10Re < 10): Np=KpReN_p = \frac{K_p}{Re}, where KpK_p is a laminar shape factor (Kp4065K_p \approx 40 - 65 for PBT, 300400\approx 300 - 400 for anchors).
  • Turbulent Zone (Re>10,000Re > 10,000): NpN_p becomes constant (Np0.30N_p \approx 0.30 for HE-3 Hydrofoil; Np1.27N_p \approx 1.27 for standard 4-blade 45° PBT; Np5.0N_p \approx 5.0 for Rushton Turbine).

# 10. The Pumping Number (NQN_Q) and Bulk Flow Rate

The Pumping Number (NQN_Q) quantifies the volumetric liquid discharge rate pumped by the impeller per revolution:

NQ=QND3N_Q = \frac{Q}{N \cdot D^3}

Rearranging to solve for total impeller discharge capacity (QQ):

Q=NQND3Q = N_Q \cdot N \cdot D^3

Where:

  • QQ = Primary volumetric pumping rate discharged from impeller blades (m3/s\text{m}^3/\text{s})
  • NQN_Q = Impeller Pumping Number (Dimensionless, typically 0.500.850.50\text{--}0.85)
  • NN = Rotational speed (rev/s\text{rev/s})
  • DD = Impeller diameter (m)

# Flow-to-Power Efficiency Metric:

To evaluate how efficiently an impeller generates bulk fluid turnover per kilowatt consumed, engineers evaluate the Pumping Efficiency Ratio (Q/PQ / P):

QP=NQND3NpρN3D5=(NQNp)1ρN2D2\frac{Q}{P} = \frac{N_Q \cdot N \cdot D^3}{N_p \cdot \rho \cdot N^3 \cdot D^5} = \left( \frac{N_Q}{N_p} \right) \cdot \frac{1}{\rho \cdot N^2 \cdot D^2}

High-efficiency hydrofoils maximize the (NQ/Np)(N_Q / N_p) ratio (0.56/0.30=1.87\approx 0.56 / 0.30 = 1.87), delivering over 5 times more bulk pumping turnover per kilowatt than a pitched blade turbine (0.75/1.30=0.580.75 / 1.30 = 0.58) or Rushton turbine (0.75/5.0=0.150.75 / 5.0 = 0.15).


# 11. Power per Unit Volume (P/VP/V)

Power per Unit Volume (P/VP/V) is the most widely cited process engineering metric for reactor agitation intensity:

PV=NpρN3D5V\frac{P}{V} = \frac{N_p \cdot \rho \cdot N^3 \cdot D^5}{V}

Where:

  • P/VP/V = Volumetric power input (W/m3\text{W/m}^3 or kW/m3\text{kW/m}^3)
  • VV = Total agitated liquid volume (m3\text{m}^3)

# Standard Industrial P/VP/V Design Benchmarks:

  • Mild Blending & Storage (0.10.3 kW/m30.1 - 0.3\text{ kW/m}^3): Equalizing bulk temperature in solvent storage tanks.
  • Moderate Chemical Synthesis (0.51.2 kW/m30.5 - 1.2\text{ kW/m}^3): Standard single-phase homogeneous reactions, dissolving soluble salts.
  • Intense Reaction & Dissolution (1.53.0 kW/m31.5 - 3.0\text{ kW/m}^3): Viscous slurries, competitive fast reactions, high-rate heat transfer.
  • Extreme Gas-Liquid & Hydrogenation (3.58.0+ kW/m33.5 - 8.0+\text{ kW/m}^3): High-pressure 3-phase hydrogenations, aerobic fermentations.

# 12. Impeller Tip Speed (vtipv_{tip}) and Shear Rates

The Impeller Tip Speed (vtipv_{tip}) represents the maximum linear peripheral velocity of the outer blade edge:

vtip=πDNv_{tip} = \pi \cdot D \cdot N

Where:

  • vtipv_{tip} = Impeller peripheral tip speed (m/s\text{m/s})
  • DD = Impeller outer diameter (m)
  • NN = Rotational speed (rev/s\text{rev/s})

# Maximum vs. Average Shear Rates:

The maximum fluid shear rate occurs directly in the narrow trailing vortex zone at the blade tip:

γ˙maxkshearN\dot{\gamma}_{max} \approx k_{shear} \cdot N

Where kshear50100k_{shear} \approx 50 - 100 for flat blades and 204020 - 40 for hydrofoils.

The Crystallization Shear Limit: In API crystallization, exceeding a critical tip speed (vtip>2.22.5 m/sv_{tip} > 2.2\text{--}2.5\text{ m/s}) shatters fragile crystal facets, generating massive populations of sub-10 micron fines that blind centrifuge filter bags and increase ANFD drying times by 300%300\%.


# 13. Agitator Shaft Torque (TqT_q)

Torque is the mechanical rotational force transmitted from the drive motor and gearbox down the shaft to overcome liquid drag:

P=2πNTq    Tq=P2πNP = 2\pi \cdot N \cdot T_q \implies T_q = \frac{P}{2\pi \cdot N}

Where:

  • TqT_q = Shaft torque (Newton-meters, Nm\text{N}\cdot\text{m})
  • PP = Shaft power (Watts, W\text{W})
  • NN = Rotational speed (rev/s\text{rev/s})

# Operating Torque vs. Starting Torque:

During cold startup or re-suspension of settled settled solids after a power failure, the Starting Torque can reach 180% to 250%180\%\text{ to }250\% of nominal steady-state operating torque. Gearboxes and solid shafts must be mechanically rated against peak starting torque.


# 14. Motor Power Selection & Mechanical Drive Sizing

The nameplate electrical motor rating (PmotorP_{motor}) must always exceed the net shaft process power (PshaftP_{shaft}) to account for mechanical transmission losses and process viscosity variations:

PmotorPshaftFsηgearηsealηmotorP_{motor} \ge \frac{P_{shaft} \cdot F_s}{\eta_{gear} \cdot \eta_{seal} \cdot \eta_{motor}}

Where:

  • ηgear\eta_{gear} = Gearbox mechanical transmission efficiency (typically 0.940.960.94 - 0.96 for premium helical drives).
  • ηseal\eta_{seal} = Mechanical seal drag efficiency (typically 0.970.990.97 - 0.99; seal friction absorbs 0.20.8 kW0.2 - 0.8\text{ kW} depending on shaft diameter and barrier fluid pressure).
  • ηmotor\eta_{motor} = Electric motor electrical efficiency (typically 0.900.950.90 - 0.95 for IE3/IE4 premium efficiency motors).
  • FsF_s = Project Engineering Service Factor (typically 1.201.351.20 - 1.35 for standard liquids; 1.501.50 for slurry/polymer reactions).

# 15. Blend Time and Mixing Kinetics (tmixt_{mix})

Mixing Time (t95t_{95}) is the time required after a tracer injection to achieve 95%95\% chemical composition homogeneity across all points in the vessel:

In the fully turbulent regime (Re>10,000Re > 10,000), the dimensionless product of rotational speed and mixing time (NtmixN \cdot t_{mix}) is an invariant geometric constant:

Ntmix=Cmix    tmix=CmixNN \cdot t_{mix} = C_{mix} \implies t_{mix} = \frac{C_{mix}}{N}

For a standard 4-blade pitched-blade turbine (D/T=0.33D/T = 0.33, fully baffled):

Cmix5.4T2D25.4(3.0)248.6C_{mix} \approx 5.4 \cdot \frac{T^2}{D^2} \approx 5.4 \cdot (3.0)^2 \approx 48.6

Therefore, at N=2.0 rev/sN = 2.0\text{ rev/s} (120 RPM):

tmix48.62.0=24.3 secondst_{mix} \approx \frac{48.6}{2.0} = \mathbf{24.3\text{ seconds}}

# 16. General Empirical Blending Correlations

For vessels operating across transitional and turbulent regimes, Grenville and Nienow established the generalized blending correlation:

Nt95=5.4Np1/3(TD)2N \cdot t_{95} = \frac{5.4}{N_p^{1/3}} \cdot \left( \frac{T}{D} \right)^2

Substituting N=(PNpρD5)1/3N = \left( \frac{P}{N_p \cdot \rho \cdot D^5} \right)^{1/3} yields the fundamental power-mixing relationship:

t95=5.4(PρV)1/3(TD)2/3t_{95} = 5.4 \cdot \left( \frac{P}{\rho \cdot V} \right)^{-1/3} \cdot \left( \frac{T}{D} \right)^{2/3}

Key Process Lesson: Blending time scales inversely with the cube root of power per unit volume ((P/V)1/3(P/V)^{-1/3})! Doubling the motor power reduces mixing time by only (2)1/3=0.794(2)^{-1/3} = 0.794 (a modest 20.6%20.6\% reduction), while increasing impeller diameter (D/TD/T) reduces blend time far more cost-effectively.


# 17. Vessel Baffles and Vortex Suppression

Operating an unbaffled vertical reactor with a centered top-entry agitator in low-viscosity liquid creates a massive central rotational vortex.

  • Unbaffled Systems: Fluid rotates as a solid-body centrifuge. Vertical top-to-bottom circulation collapses; power draw drops by up to 75%75\%; central vortex draws air/nitrogen into the liquid, causing severe frothing and shaft whipping.
  • Fully Baffled Systems: Vertical wall baffles intercept circular flow vectors, redirecting energy into powerful axial turnover loops and maximizing motor power transfer into turbulent micro-mixing.

# Non-Negotiable Baffle Rules:

  1. Standard Baffle Count: 4 vertical wall baffles spaced at 90° intervals around the internal circumference.
  2. Standard Baffle Width: B=T/10B = T / 10 (for general synthesis) or B=T/12B = T / 12 (for high-efficiency hydrofoils).
  3. Wall Offset Clearance: Baffles must be stood off from the vessel shell by Bgap=T/50B_{gap} = T / 50 (or 2538 mm25\text{--}38\text{ mm}) to prevent solids accumulation and promote self-cleaning during CIP.

# 18. Multi-Stage Impeller Assemblies

When the liquid depth-to-tank diameter ratio exceeds H/T>1.2H/T > 1.2, a single impeller cannot establish complete top-to-bottom circulation, creating isolated stagnant zones in the upper liquid layer.

# Design Criteria for Dual-Impeller Systems:

  • Lower Impeller (C=0.250.33TC = 0.25\text{--}0.33 T): High-efficiency downward axial hydrofoil or PBT to sweep the dished bottom head.
  • Upper Impeller (Submerged 0.50.7H0.5\text{--}0.7 H below surface): Downward axial PBT or hydrofoil to draw fresh liquid feeds downward into the core.
  • Axial Spacing (SS): Must be maintained between 1.0DS1.5D1.0 D \le S \le 1.5 D. Spacing <1.0D< 1.0 D causes flow coupling where both stages act as a single inefficient thick impeller; spacing >2.0D> 2.0 D splits the fluid into two disconnected circulation loops with zero inter-stage exchange.

# 19. Solid-Liquid Suspension & The Zwietering NjsN_{js} Correlation

For heterogeneous catalytic reactions (e.g., hydrogenation using Pd/C, Raney Nickel) or dissolving solid API crystal feeds, the agitator must operate at or above the Just-Suspended Speed (NjsN_{js}), defined as the minimum rotational speed at which no solid particle remains stationary on the vessel bottom for longer than 1 to 2 seconds (ASTM / Zwietering 1-second criterion).

# The Zwietering Fundamental Equation:

Njs=Sν0.1[g(ρsρl)ρl]0.45dp0.2X0.13D0.85N_{js} = S \cdot \nu^{0.1} \cdot \left[ \frac{g \cdot (\rho_s - \rho_l)}{\rho_l} \right]^{0.45} \cdot d_p^{0.2} \cdot X^{0.13} \cdot D^{-0.85}

Where:

  • NjsN_{js} = Just-suspended impeller speed (rev/s\text{rev/s})
  • SS = Dimensionless Zwietering geometry constant (function of impeller type, D/TD/T, and C/TC/T)
  • ν\nu = Kinematic viscosity of liquid (m2/s=μ/ρl\text{m}^2/\text{s} = \mu / \rho_l)
  • gg = Gravitational acceleration (9.81 m/s29.81\text{ m/s}^2)
  • ρs\rho_s = Solid particle density (kg/m3\text{kg/m}^3)
  • ρl\rho_l = Liquid carrier density (kg/m3\text{kg/m}^3)
  • dpd_p = Mean solid particle diameter (m)
  • XX = Solids mass concentration ratio (100×mass solids/mass liquid100 \times \text{mass solids} / \text{mass liquid}, %\%)
  • DD = Impeller diameter (m)

# Zwietering Constant (SS) Reference Values (T/3T/3 Clearance, Fully Baffled):

  • High-Efficiency Hydrofoil (A310 / HE-3): S5.26.5S \approx 5.2 - 6.5
  • 4-Blade Pitched Blade Turbine (45°): S6.88.2S \approx 6.8 - 8.2
  • 6-Blade Rushton Radial Turbine: S8.511.5S \approx 8.5 - 11.5 (Very inefficient for solids suspension)

# 20. Gas-Liquid Agitation, Mass Transfer (kLak_L a), and Gassed Power

In gas-liquid reactions (e.g., catalytic hydrogenation, oxidation, halogenation), the chemical reaction rate is governed by the volumetric gas-liquid mass transfer rate (kLak_L a):

Rgas=kLa(CCL)R_{gas} = k_L a \cdot (C^* - C_L)

# The Van't Riet Correlation for Coalescing Organic Systems:

kLa=0.026(PgV)0.4vs0.5k_L a = 0.026 \cdot \left( \frac{P_g}{V} \right)^{0.4} \cdot v_s^{0.5}

Where:

  • kLak_L a = Volumetric mass transfer coefficient (s1\text{s}^{-1})
  • Pg/VP_g / V = Gassed power per unit volume (W/m3\text{W/m}^3)
  • vsv_s = Superficial gas velocity through the vessel cross section (m/s=Qgas/[π4T2]\text{m/s} = Q_{gas} / [\frac{\pi}{4} T^2])

# The Gassed Power Drop (Pg/PuP_g / P_u):

When gas is sparged beneath an impeller, low-density gas cavities form behind the blades, reducing the effective density and drag. The Gassed Power (PgP_g) drops to 40% to 70%40\%\text{ to }70\% of the Ungassed Power (PuP_u).

Severe Design Trap: If gas supply trips during a live commercial reaction, the agitator immediately transitions from gassed to ungassed conditions. If the drive motor was sized strictly for gassed power (PgP_g), the sudden jump to ungassed power (PuP_u) will trip the motor on electrical overload! Always size drive motors for 100% ungassed liquid power.


# 21. Agitation Effects on Reactor Heat Transfer

In jacketed synthesis reactors, agitation directly governs the Inside Process-Side Heat Transfer Film Coefficient (hih_i):

Q=UAwettedΔTlmQ = U \cdot A_{wetted} \cdot \Delta T_{lm}
1U=1hi+twallkwall+1hj+Rfouling\frac{1}{U} = \frac{1}{h_i} + \frac{t_{wall}}{k_{wall}} + \frac{1}{h_j} + R_{fouling}

# Generalized Nusselt Correlation for Agitated Jacketed Vessels:

Nu=hiTkf=ChtRe2/3Pr1/3(μμw)0.14Nu = \frac{h_i \cdot T}{k_f} = C_{ht} \cdot Re^{2/3} \cdot Pr^{1/3} \cdot \left( \frac{\mu}{\mu_w} \right)^{0.14}

Where:

  • NuNu = Nusselt Number (Dimensionless)
  • PrPr = Prandtl Number of fluid (Pr=Cpμ/kfPr = C_p \cdot \mu / k_f)
  • kfk_f = Thermal conductivity of fluid (W/(mK)\text{W}/(\text{m}\cdot\text{K}))
  • μw\mu_w = Fluid viscosity at vessel wall temperature
  • ChtC_{ht} = Empirical heat transfer constant (Cht0.54C_{ht} \approx 0.54 for PBT; Cht0.74C_{ht} \approx 0.74 for Rushton turbines; Cht0.36C_{ht} \approx 0.36 for anchors in laminar regime)

# 22. Agitator Design for High-Viscosity Systems

When fluid viscosity rises above 5,000 mPas5,000\text{ mPa}\cdot\text{s} (e.g., polymerizations, concentrated API slurries, fermentation broths), conventional open turbines fail completely because turbulent kinetic energy decays within millimeters of the blade tip, creating an isolated rotating liquid cylinder (the "cavern effect") surrounded by stagnant fluid.

# Close-Clearance Impeller Dynamics:

  • Anchor Agitators: Clearance from vessel wall wclearance=515 mmw_{clearance} = 5\text{--}15\text{ mm} (D/T0.950.98D/T \approx 0.95 - 0.98). Sweeps viscous boundary layers off the cooled jacket wall, preventing thermal scorching.
  • Helical Ribbon Agitators: Double-flight helical ribbons physically lift fluid up the vessel wall and pump it down the center core, enforcing positive turnover even at Re<1.0Re < 1.0.

# 23. Scale-Up of Batch Reactor Agitators: The 6 Governing Criteria

When scaling an agitation system from lab/pilot (11) to commercial scale (22), process engineers must select a single primary governing scale-up criterion:

Scale-Up CriterionMathematical Governing RuleScale Exponent (N2=N1[D1/D2]nN_2 = N_1 \cdot [D_1/D_2]^n)Typical Industrial Application
1. Constant Power / Volume (P/VP/V)(P/V)1=(P/V)2(P/V)_1 = (P/V)_2n=2/3=0.67n = 2/3 = 0.67Mass-transfer limited synthesis, fast competitive reactions.
2. Constant Tip Speed (vtipv_{tip})(vtip)1=(vtip)2(v_{tip})_1 = (v_{tip})_2n=1.00n = 1.00Shear-sensitive crystallization, mammalian cell culture.
3. Constant Blend Time (tmixt_{mix})(tmix)1=(tmix)2(t_{mix})_1 = (t_{mix})_2n=0.00n = 0.00 (N2=N1N_2 = N_1)Unrealistic in large plants; causes explosive power demand!
4. Constant Reynolds Number (ReRe)Re1=Re2Re_1 = Re_2n=2.00n = 2.00Strictly hydrodynamic study; impractical for mass transfer.
5. Constant Solids Suspension (NjsN_{js})(Njs)1=(Njs)2(N_{js})_1 = (N_{js})_2n=0.85n = 0.85 (Zwietering)Slurry dissolving, heterogeneous catalyst suspension.
6. Constant Mass Transfer (kLak_L a)(kLa)1=(kLa)2(k_L a)_1 = (k_L a)_2n=0.670.75n = 0.67 - 0.75Hydrogenation, gas-liquid sparged reactions.

# 24. Why RPM Should Never Be Scaled Directly

Consider scaling a process from a laboratory reactor (D1=0.05 mD_1 = 0.05\text{ m}) running at N1=600 RPMN_1 = 600\text{ RPM} to a commercial 5 KL5\text{ KL} reactor (D2=0.55 mD_2 = 0.55\text{ m}, scale factor λ=D2/D1=11.0\lambda = D_2 / D_1 = 11.0):

  • If scaling directly on Constant RPM (N2=600 RPMN_2 = 600\text{ RPM}):
P2=P1(N2N1)3(D2D1)5=P1(1)3(11)5=161,051×P1P_2 = P_1 \cdot \left(\frac{N_2}{N_1}\right)^3 \cdot \left(\frac{D_2}{D_1}\right)^5 = P_1 \cdot (1)^3 \cdot (11)^5 = \mathbf{161,051 \times P_1}

The power requirement explodes by over 160,000-fold, destroying the shaft and motor!

  • If scaling on Constant Power per Unit Volume (P/VP/V):
N2=N1(D1D2)2/3=600(111)0.667=121.3 RPMN_2 = N_1 \cdot \left(\frac{D_1}{D_2}\right)^{2/3} = 600 \cdot \left(\frac{1}{11}\right)^{0.667} = \mathbf{121.3\text{ RPM}}

The commercial agitator operates at 121 RPM121\text{ RPM}, achieving identical volumetric power input.

  • If scaling on Constant Tip Speed (vtipv_{tip}):
N2=N1(D1D2)1.0=600(111)=54.5 RPMN_2 = N_1 \cdot \left(\frac{D_1}{D_2}\right)^{1.0} = 600 \cdot \left(\frac{1}{11}\right) = \mathbf{54.5\text{ RPM}}

# 25. Complete Multi-Scale Scale-Up Comparison Table

The following table demonstrates the hydrodynamic shift across five geometrically similar vessel scales (D/T=0.35D/T = 0.35, H/T=1.1H/T = 1.1, ρ=950 kg/m3\rho = 950\text{ kg/m}^3, μ=1.2 cP\mu = 1.2\text{ cP}, Np=1.30N_p = 1.30) scaled on Constant Power per Volume (P/V=1.0 kW/m3P/V = 1.0\text{ kW/m}^3):

ParameterLab Scale (1 L)Bench Scale (100 L)Pilot Scale (1 KL)Plant Scale (5 KL)Production Scale (10 KL)
Vessel Diameter (TT)0.105 m0.105\text{ m}0.485 m0.485\text{ m}1.05 m1.05\text{ m}1.60 m1.60\text{ m}2.02 m2.02\text{ m}
Impeller Diameter (DD)0.037 m0.037\text{ m}0.170 m0.170\text{ m}0.368 m0.368\text{ m}0.560 m0.560\text{ m}0.707 m0.707\text{ m}
Agitator Speed (NN)812 RPM812\text{ RPM}299 RPM299\text{ RPM}175 RPM175\text{ RPM}132 RPM132\text{ RPM}113 RPM113\text{ RPM}
Net Shaft Power (PP)0.001 kW0.001\text{ kW}0.100 kW0.100\text{ kW}1.00 kW1.00\text{ kW}5.00 kW5.00\text{ kW}10.00 kW10.00\text{ kW}
Power/Volume (P/VP/V)1.00 kW/m31.00\text{ kW/m}^31.00 kW/m31.00\text{ kW/m}^31.00 kW/m31.00\text{ kW/m}^31.00 kW/m31.00\text{ kW/m}^31.00 kW/m31.00\text{ kW/m}^3
Tip Speed (vtipv_{tip})1.57 m/s1.57\text{ m/s}2.66 m/s2.66\text{ m/s}3.37 m/s3.37\text{ m/s}3.87 m/s3.87\text{ m/s}4.18 m/s4.18\text{ m/s}
Reynolds Number (ReRe)1.48×1041.48 \times 10^41.08×1051.08 \times 10^53.00×1053.00 \times 10^55.21×1055.21 \times 10^57.12×1057.12 \times 10^5
Blend Time (t95t_{95})3.6 s3.6\text{ s}9.8 s9.8\text{ s}16.7 s16.7\text{ s}22.1 s22.1\text{ s}25.8 s25.8\text{ s}
Shaft Torque (TqT_q)0.012 Nm0.012\text{ N}\cdot\text{m}3.19 Nm3.19\text{ N}\cdot\text{m}54.6 Nm54.6\text{ N}\cdot\text{m}361.7 Nm361.7\text{ N}\cdot\text{m}845.2 Nm845.2\text{ N}\cdot\text{m}

# 26. Practical Engineering Worked Example: Sizing a 5 KL Batch Reactor

# Process Design Input Data:

  • Vessel Nominal Volume: 5.0 m35.0\text{ m}^3 (5,000 L5,000\text{ L})
  • Operating Working Batch Volume: V=4.0 m3V = 4.0\text{ m}^3 (4,000 L4,000\text{ L})
  • Process Fluid Density: ρ=920 kg/m3\rho = 920\text{ kg/m}^3
  • Process Fluid Dynamic Viscosity: μ=2.4 cP=0.0024 Pas\mu = 2.4\text{ cP} = 0.0024\text{ Pa}\cdot\text{s}
  • Target Mixing Intensity: Moderate Chemical Synthesis (P/V=1.0 kW/m3P/V = 1.0\text{ kW/m}^3)

# Step 1: Establish Vessel Geometry

  • Choosing standard aspect ratio H/T=1.15H/T = 1.15:
V=π4T2H=π4T2(1.15T)=0.9032T3=4.0 m3V = \frac{\pi}{4} T^2 \cdot H = \frac{\pi}{4} T^2 \cdot (1.15 T) = 0.9032 \cdot T^3 = 4.0\text{ m}^3
T=(4.00.9032)1/3=1.64 m(Specify Standard Shell OD: 1,600 mm ID)T = \left( \frac{4.0}{0.9032} \right)^{1/3} = \mathbf{1.64\text{ m}} \quad \text{(Specify Standard Shell OD: 1,600 mm ID)}
  • Liquid height: H=1.15×1.64=1.89 mH = 1.15 \times 1.64 = \mathbf{1.89\text{ m}}
  • Impeller Diameter (D/T=0.35D/T = 0.35): D=0.35×1.60=0.56 mD = 0.35 \times 1.60 = \mathbf{0.56\text{ m}} (560 mm560\text{ mm})

# Step 2: Select Impeller Geometry & Power Number

  • Dual Impeller Selection: Lower 4-Blade PBT (Np1=1.30N_{p1} = 1.30), Upper 3-Blade Hydrofoil (Np2=0.35N_{p2} = 0.35).
  • Combined Assembly Power Number: Np,total=1.30+0.35=1.65N_{p,total} = 1.30 + 0.35 = \mathbf{1.65}

# Step 3: Calculate Required Rotational Speed (NN)

  • Target total shaft power: Pshaft=(P/V)×V=1,000 W/m3×4.0 m3=4,000 W=4.0 kWP_{shaft} = (P/V) \times V = 1,000\text{ W/m}^3 \times 4.0\text{ m}^3 = \mathbf{4,000\text{ W}} = 4.0\text{ kW}
  • Solving for NN:
P=NpρN3D5    N=(PNpρD5)1/3P = N_p \cdot \rho \cdot N^3 \cdot D^5 \implies N = \left( \frac{P}{N_p \cdot \rho \cdot D^5} \right)^{1/3}
N=(4,0001.65920(0.56)5)1/3=(4,0001,5180.05507)1/3=(47.85)1/3=3.63 rev/s218 RPMN = \left( \frac{4,000}{1.65 \cdot 920 \cdot (0.56)^5} \right)^{1/3} = \left( \frac{4,000}{1,518 \cdot 0.05507} \right)^{1/3} = (47.85)^{1/3} = \mathbf{3.63\text{ rev/s}} \approx \mathbf{218\text{ RPM}}

# Step 4: Verify Hydrodynamic Regimes & Safety Checks

  • Reynolds Number:
Re=9203.63(0.56)20.0024=4.36×105(Fully Turbulent, Re>10,000    Np is valid)Re = \frac{920 \cdot 3.63 \cdot (0.56)^2}{0.0024} = \mathbf{4.36 \times 10^5} \quad (\text{Fully Turbulent, } Re > 10,000 \implies N_p \text{ is valid})
  • Impeller Tip Speed:
vtip=πDN=3.14160.563.63=6.39 m/sv_{tip} = \pi \cdot D \cdot N = 3.1416 \cdot 0.56 \cdot 3.63 = \mathbf{6.39\text{ m/s}}
  • Shaft Operating Torque:
Tq=4,0002π3.63=175.4 NmT_q = \frac{4,000}{2\pi \cdot 3.63} = \mathbf{175.4\text{ N}\cdot\text{m}}

# Step 5: Motor & Gearbox Rating

  • Applying mechanical efficiencies (ηgear=0.95,ηseal=0.98\eta_{gear} = 0.95, \eta_{seal} = 0.98) and Service Factor Fs=1.30F_s = 1.30:
Pmotor=4.0 kW×1.300.95×0.98=5.59 kWP_{motor} = \frac{4.0\text{ kW} \times 1.30}{0.95 \times 0.98} = \mathbf{5.59\text{ kW}}
  • Selected Commercial Standard Motor: 7.5 kW7.5\text{ kW} (10 HP), 4-Pole (1440 RPM nominal), 415V / 50 Hz, Ex-d Flameproof
  • Gearbox Reduction Ratio: Ratio =1440/2181:6.6= 1440 / 218 \approx \mathbf{1:6.6} (Helical Bevel Geared Unit)

# 27. Shaft Mechanical Design, ASME Code Sizing & Critical Speed Dynamics

A top-entry batch reactor agitator shaft operates as an extended, vertical overhung cantilever beam. It is subjected to simultaneous, complex mechanical stresses:

  1. Continuous Torsional Shear Stress (\tau): Transmitted from the electric motor and gearbox down the shaft to overcome fluid drag.
  2. Cyclic Reverse Bending Stress (\sigma_b): Induced by dynamic, asymmetric fluid hydraulic side loads (FhF_h) acting on the rotating impeller blades as they pass near baffles or liquid level interfaces.
  3. Axial Tensile / Compressive Stress: From impeller weight and hydraulic down-thrust/up-thrust forces.

# A. ASME Code Governing Equations for Solid Shaft Diameter (dsd_s)

Under the ASME Code for Design of Transmission Shafting (incorporating Guest's Maximum Shear Stress Theory and Rankine's Maximum Normal Stress Theory), the minimum solid shaft diameter is sized from the Equivalent Combined Torsional Moment (TeT_e):

Te=(kmMb)2+(ktTq)2T_e = \sqrt{(k_m \cdot M_b)^2 + (k_t \cdot T_q)^2}
Me=12[kmMb+(kmMb)2+(ktTq)2]M_e = \frac{1}{2} \left[ k_m \cdot M_b + \sqrt{(k_m \cdot M_b)^2 + (k_t \cdot T_q)^2} \right]

The minimum required solid shaft diameter (dsd_s) is calculated as:

ds=[16Teπτallow]1/3=[16πτallow(kmMb)2+(ktTq)2]1/3d_s = \left[ \frac{16 \cdot T_e}{\pi \cdot \tau_{allow}} \right]^{1/3} = \left[ \frac{16}{\pi \cdot \tau_{allow}} \cdot \sqrt{(k_m \cdot M_b)^2 + (k_t \cdot T_q)^2} \right]^{1/3}

Where:

  • TqT_q = Transmitted operating torque (Nm\text{N}\cdot\text{m}):
Tq=Pshaftω=Pshaft×602πN=9550PkWNrpmT_q = \frac{P_{shaft}}{\omega} = \frac{P_{shaft} \times 60}{2\pi \cdot N} = \frac{9550 \cdot P_{kW}}{N_{rpm}}
  • FhF_h = Hydraulic lateral side force acting at the lowest impeller (N\text{N}):
Fh=fhTqD/2F_h = \frac{f_h \cdot T_q}{D / 2}

(where fhf_h is the dimensionless hydraulic side-load factor: 0.300.400.30 - 0.40 for single-phase liquid blending; 0.500.700.50 - 0.70 for gas dispersion, liquid level draw-down, or solid suspension)

  • MbM_b = Maximum bending moment at the lower support bearing (Nm\text{N}\cdot\text{m}):
Mb=FhLM_b = F_h \cdot L

(where LL is the total overhung length from the lower support bearing to the lowest impeller hub)

  • kmk_m = Combined shock and fatigue factor for bending (1.52.01.5 - 2.0 for rotating shafts with minor to heavy shock)
  • ktk_t = Combined shock and fatigue factor for torsion (1.01.51.0 - 1.5)
  • τallow\tau_{allow} = Maximum allowable shear stress for the shaft material (Pa\text{Pa} or MPa\text{MPa}):
τallow=min(0.30Syield,0.18Sultimate)×Keyway Factor (0.75)\tau_{allow} = \min\left( 0.30 \cdot S_{yield}, \, 0.18 \cdot S_{ultimate} \right) \times \text{Keyway Factor (0.75)}

(For standard annealed SS316L: τallow4550 MPa\tau_{allow} \approx 45 - 50\text{ MPa}; for Hastelloy C-22: τallow6575 MPa\tau_{allow} \approx 65 - 75\text{ MPa})


# B. Deflection Limits at the Mechanical Seal & Impeller Tip

Excessive shaft deflection is the #1 root cause of premature mechanical seal failure and batch contamination in pharmaceutical reactors:

  1. Deflection at the Mechanical Seal Face (δseal\delta_{seal}):
δseal=FhLseal2(3LLseal)6EI0.0500.080 mm(0.0020.003 in)\delta_{seal} = \frac{F_h \cdot L_{seal}^2 \cdot (3L - L_{seal})}{6 \cdot E \cdot I} \le \mathbf{0.050 - 0.080\text{ mm}} \quad (0.002 - 0.003\text{ in})

(where LsealL_{seal} is the distance from the lower support bearing to the seal faces, and I=πds464I = \frac{\pi \cdot d_s^4}{64})
2. Deflection at the Lowest Impeller Tip (δtip\delta_{tip}):

δtip=FhL33EI5.010.0 mm(Must prevent blade-to-wall/baffle collision)\delta_{tip} = \frac{F_h \cdot L^3}{3 \cdot E \cdot I} \le \mathbf{5.0 - 10.0\text{ mm}} \quad (\text{Must prevent blade-to-wall/baffle collision})

# C. Critical Whirling Speed (NcritN_{crit}) & Resonance Avoidance

A cantilevered agitator shaft possesses natural lateral vibration frequencies. If the operating rotational speed (NN) coincides with a natural frequency, violent shaft whirling / resonance occurs:

Using the Rayleigh-Ritz Method for Overhung Shafts with Concentrated Impeller Masses:

Ncrit=30πgδstaticN_{crit} = \frac{30}{\pi} \cdot \sqrt{\frac{g}{\delta_{static}}}

Where δstatic\delta_{static} is the total static lateral deflection under gravitational and concentrated blade loads:

δstatic=(mimpellers+0.25mshaft)gL33EI\delta_{static} = \frac{(m_{impellers} + 0.25 \cdot m_{shaft}) \cdot g \cdot L^3}{3 \cdot E \cdot I}
  • Rigid Shaft Design (Sub-Critical — Mandatory for Top-Entry Pharma Reactors):
NNcrit0.70(Operating speed must be at least 30% below first critical speed)\frac{N}{N_{crit}} \le \mathbf{0.70} \quad (\text{Operating speed must be at least } 30\% \text{ below first critical speed})
  • Flexible Shaft Design (Super-Critical):
NNcrit1.30(Requires rapid VFD acceleration through the resonant band 0.801.20Ncrit)\frac{N}{N_{crit}} \ge \mathbf{1.30} \quad (\text{Requires rapid VFD acceleration through the resonant band } 0.80 - 1.20 N_{crit})

# D. Comprehensive Worked Case Study: Minimum Shaft Diameter Sizing for a 5 KL Commercial Batch Reactor

# Design Problem Statement:

Size the minimum solid shaft diameter, evaluate operating torque, bending moments, seal face deflections, and verify critical speed resonance for a 5,000 L (5 KL5\text{ KL}) Pharmaceutical Synthesis Reactor (R-101).

# Input Design Specifications:

  • Working Batch Volume (VV): 4.0 m34.0\text{ m}^3 (4,000 L4,000\text{ L})
  • Liquid Density (ρ\rho): 1,050 kg/m31,050\text{ kg/m}^3, Viscosity μ=0.015 Pas\mu = 0.015\text{ Pa}\cdot\text{s}
  • Dual-Impeller Assembly: Lower 4-Blade PBT (D1=0.60 mD_1 = 0.60\text{ m}) + Upper PBT (D2=0.60 mD_2 = 0.60\text{ m}), Total Blade Mass mimp=35.0 kgm_{imp} = 35.0\text{ kg}
  • Operating Speed (NN): 120 rpm\mathbf{120\text{ rpm}} (Ns=2.0 rev/sN_s = 2.0\text{ rev/s})
  • Total Absorbed Process Power (PP): 7.50 kW\mathbf{7.50\text{ kW}} (7,500 W7,500\text{ W})
  • Overhung Shaft Length (LL): 2.40 m\mathbf{2.40\text{ m}} (distance from lower bearing in drive lantern to bottom impeller)
  • Seal Face Location (LsealL_{seal}): 0.30 m0.30\text{ m} below lower support bearing
  • Shaft Material of Construction: Forged SS316L (E=193 GPa=193×109 PaE = 193\text{ GPa} = 193 \times 10^9\text{ Pa}, ρsteel=7,950 kg/m3\rho_{steel} = 7,950\text{ kg/m}^3, τallow=48.0 MPa\tau_{allow} = 48.0\text{ MPa})
  • ASME Safety Factors: Shock factor km=1.75k_m = 1.75, Fatigue factor kt=1.25k_t = 1.25, Hydraulic side factor fh=0.40f_h = 0.40, Design Safety Factor SF=1.50SF = 1.50

# Step-by-Step Calculation:

# Step 1: Calculate Transmitted Operating Torque (TqT_q)

ω=2πN60=2π×12060=12.566 rad/s\omega = \frac{2\pi \cdot N}{60} = \frac{2\pi \times 120}{60} = 12.566\text{ rad/s}
Tq=Pω=7,500 W12.566 rad/s=596.8 NmT_q = \frac{P}{\omega} = \frac{7,500\text{ W}}{12.566\text{ rad/s}} = \mathbf{596.8\text{ N}\cdot\text{m}}

Applying the design safety factor (SF=1.50SF = 1.50 for starting load):

Tq,design=1.50×596.8=895.2 NmT_{q,design} = 1.50 \times 596.8 = \mathbf{895.2\text{ N}\cdot\text{m}}

# Step 2: Calculate Fluid Hydraulic Lateral Side Load (FhF_h)

Acting at the lower impeller radius (D/2=0.30 mD/2 = 0.30\text{ m}):

Fh=fhTqD/2=0.40×596.8 Nm0.30 m=795.7 NF_h = \frac{f_h \cdot T_q}{D / 2} = \frac{0.40 \times 596.8\text{ N}\cdot\text{m}}{0.30\text{ m}} = \mathbf{795.7\text{ N}}

# Step 3: Calculate Maximum Bending Moment at Lower Support Bearing (MbM_b)

Mb=FhL=795.7 N×2.40 m=1,909.7 NmM_b = F_h \cdot L = 795.7\text{ N} \times 2.40\text{ m} = \mathbf{1,909.7\text{ N}\cdot\text{m}}

# Step 4: Calculate ASME Equivalent Combined Torsional Moment (TeT_e)

Te=(kmMb)2+(ktTq,design)2T_e = \sqrt{(k_m \cdot M_b)^2 + (k_t \cdot T_{q,design})^2}
Te=(1.75×1909.7)2+(1.25×895.2)2=(3342.0)2+(1119.0)2=3,524.3 NmT_e = \sqrt{(1.75 \times 1909.7)^2 + (1.25 \times 895.2)^2} = \sqrt{(3342.0)^2 + (1119.0)^2} = \mathbf{3,524.3\text{ N}\cdot\text{m}}

# Step 5: Determine Minimum Required Solid Shaft Diameter (dmind_{min})

Using allowable shear stress for SS316L (τallow=48.0 MPa=48.0×106 Pa\tau_{allow} = 48.0\text{ MPa} = 48.0 \times 10^6\text{ Pa}):

dmin=[16Teπτallow]1/3=[16×3,524.3π×48.0×106]1/3=[3.739×104]1/3=0.0720 m=72.0 mmd_{min} = \left[ \frac{16 \cdot T_e}{\pi \cdot \tau_{allow}} \right]^{1/3} = \left[ \frac{16 \times 3,524.3}{\pi \times 48.0 \times 10^6} \right]^{1/3} = \left[ 3.739 \times 10^{-4} \right]^{1/3} = \mathbf{0.0720\text{ m}} = \mathbf{72.0\text{ mm}}
  • Commercial Standard Metric Shaft Diameters: 65 mm, 70 mm, 75 mm, 80 mm, 90 mm, 100 mm.
  • Selected Standard Shaft Diameter: Select ds=80.0 mmd_s = \mathbf{80.0\text{ mm}} (0.080 m0.080\text{ m}) solid precision-ground bar stock.

# Step 6: Verify Deflection at Mechanical Seal Faces (δseal\delta_{seal})

For ds=0.080 md_s = 0.080\text{ m}:

I=πds464=π×(0.080)464=2.0106×106 m4I = \frac{\pi \cdot d_s^4}{64} = \frac{\pi \times (0.080)^4}{64} = 2.0106 \times 10^{-6}\text{ m}^4

Deflection at seal position (Lseal=0.30 mL_{seal} = 0.30\text{ m}):

δseal=FhLseal2(3LLseal)6EI\delta_{seal} = \frac{F_h \cdot L_{seal}^2 \cdot (3L - L_{seal})}{6 \cdot E \cdot I}
δseal=795.7×(0.30)2×(3×2.400.30)6×(193×109)×(2.0106×106)=71.613×6.902.328×106=2.12×104 m=0.042 mm\delta_{seal} = \frac{795.7 \times (0.30)^2 \times (3 \times 2.40 - 0.30)}{6 \times (193 \times 10^9) \times (2.0106 \times 10^{-6})} = \frac{71.613 \times 6.90}{2.328 \times 10^6} = 2.12 \times 10^{-4}\text{ m} = \mathbf{0.042\text{ mm}}

Mechanical Seal Check: δseal=0.042 mm0.080 mm\delta_{seal} = \mathbf{0.042\text{ mm}} \le \mathbf{0.080\text{ mm}} maximum allowable limit. PASSED! (Seal faces will maintain parallel alignment without premature face separation or leakage).


# Step 7: Verify First Critical Whirling Speed (NcritN_{crit}) & Resonance

  1. Bare shaft volume & mass:
Vshaft=πds24L=π×(0.080)24×2.40=0.01206 m3V_{shaft} = \frac{\pi \cdot d_s^2}{4} \cdot L = \frac{\pi \times (0.080)^2}{4} \times 2.40 = 0.01206\text{ m}^3
mshaft=Vshaftρsteel=0.01206×7,950=95.9 kgm_{shaft} = V_{shaft} \cdot \rho_{steel} = 0.01206 \times 7,950 = \mathbf{95.9\text{ kg}}
  1. Effective vibrating mass:
meff=mimpellers+0.25mshaft=35.0+(0.25×95.9)=58.98 kgm_{eff} = m_{impellers} + 0.25 \cdot m_{shaft} = 35.0 + (0.25 \times 95.9) = \mathbf{58.98\text{ kg}}
  1. Static gravitational deflection (δstatic\delta_{static}):
δstatic=meffgL33EI=58.98×9.81×(2.40)33×(193×109)×(2.0106×106)=8,007.81.164×106=6.879×103 m=6.88 mm\delta_{static} = \frac{m_{eff} \cdot g \cdot L^3}{3 \cdot E \cdot I} = \frac{58.98 \times 9.81 \times (2.40)^3}{3 \times (193 \times 10^9) \times (2.0106 \times 10^{-6})} = \frac{8,007.8}{1.164 \times 10^6} = 6.879 \times 10^{-3}\text{ m} = \mathbf{6.88\text{ mm}}
  1. First critical whirling speed (NcritN_{crit}):
Ncrit=30π9.816.879×103=9.549×1,426.1=9.549×37.76=360.6 rpmN_{crit} = \frac{30}{\pi} \cdot \sqrt{\frac{9.81}{6.879 \times 10^{-3}}} = 9.549 \times \sqrt{1,426.1} = 9.549 \times 37.76 = \mathbf{360.6\text{ rpm}}
  1. Operating speed ratio:
NNcrit=120 rpm360.6 rpm=0.3330.70\frac{N}{N_{crit}} = \frac{120\text{ rpm}}{360.6\text{ rpm}} = \mathbf{0.333} \le \mathbf{0.70}

Dynamic Whirling Check: Operating speed is 33.3%33.3\% of critical speed, well below the 0.700.70 rigid-shaft ceiling. The shaft will run with zero resonance, zero harmonic whipping, and maximum bearing life.


# 28. Agitator Mechanical Seal Design

A reactor mechanical seal operates under far harsher conditions than a pump seal due to long shaft overhang, dynamic shaft runout (1.0 mm\le 1.0\text{ mm}), vapor-phase operation, and thermal cycling.

# Mechanical Seal Architecture:

  • Double Pressurized Liquid-Lubricated Seal (Standard): Employs dual back-to-back or face-to-face silicon carbide / carbon seal faces with a dedicated Thermosiphon Barrier Fluid Pot (operating with white oil, propylene glycol, or compatible solvent at 1.52.0 bar1.5\text{--}2.0\text{ bar} overpressure).
  • Dry-Running Non-Contacting Gas Seal: Uses nitrogen gas barrier film; ideal for ultra-pure API synthesis where zero barrier fluid contamination is tolerated.

# 29. API Process Impeller Selection Decision Matrix

Process OperationRecommended Impeller ConfigurationTarget P/VP/V RangeCritical Scale-Up Governing Rule
Fast Homogeneous Liquid ReactionDual Pitched Blade Turbines (4545^\circ)1.22.5 kW/m31.2 - 2.5\text{ kW/m}^3Constant Power per Unit Volume (P/VP/V)
High-Pressure Gas HydrogenationGas-Induction / Cavity-Dispersion Radial3.56.0 kW/m33.5 - 6.0\text{ kW/m}^3Constant kLak_L a and Superficial Gas Velocity (vsv_s)
Shear-Sensitive CrystallizationLow-Shear High-Efficiency Axial Hydrofoil0.20.5 kW/m30.2 - 0.5\text{ kW/m}^3Constant Tip Speed (vtip<2.0 m/sv_{tip} < 2.0\text{ m/s})
Dense Slurry Catalyst SuspensionWide-Blade Hydrofoil or Downward PBT0.81.8 kW/m30.8 - 1.8\text{ kW/m}^3Zwietering Just-Suspended Speed (NjsN_{js})
High-Viscosity Polymerization (>20 Pas> 20\text{ Pa}\cdot\text{s})Double-Flight Helical Ribbon / Anchor2.05.0 kW/m32.0 - 5.0\text{ kW/m}^3Bulk Wall Turnover & Shear-Thinning Dynamics

# 30. 15 Common Agitator Design Mistakes in Plant Scale-Up

  1. Selecting the Motor Before Sizing the Impeller: Buying a 10 HP motor based on catalogue guessing rather than calculating net shaft power from NpρN3D5N_p \rho N^3 D^5.
  2. Treating RPM as an Absolute Mixing Metric: Believing 100 RPM provides the same mixing across 100 L and 10 KL reactors.
  3. Ignoring Viscosity & Reynolds Regime: Applying turbulent power numbers (NpN_p) to viscous transitional or laminar slurries.
  4. Direct RPM Linear Scale-Up: Scaling RPM directly (N2=N1N_2 = N_1), which causes power to explode by (D2/D1)5(D_2/D_1)^5.
  5. Assuming Constant P/VP/V Is Universally Applicable: Using constant P/VP/V for crystallization scale-up, causing severe crystal attrition and filter blinding.
  6. Operating Unbaffled Vessels: Leaving a reactor unbaffled in turbulent liquid, resulting in solid-body swirl, zero axial turnover, and a vortex drawing air into the pump suction.
  7. Ignoring Gassed vs. Ungassed Power Shift: Sizing motor strictly for gassed conditions; motor trips when gas is stopped during liquid priming.
  8. Neglecting Mechanical Seal Barrier Pressure: Running barrier fluid pressure lower than reactor internal pressure, forcing toxic reaction vapor into the seal faces.
  9. Operating in the Critical Speed Resonance Band (0.81.2Ncrit0.8 - 1.2 N_{crit}): Causing catastrophic shaft whipping, bearing failure, and shattered seal faces.
  10. Inadequate Inter-Impeller Spacing (S<1.0DS < 1.0 D): Squeezing dual impellers too close together, destroying hydrodynamic efficiency.
  11. Using Single-Zone Agitators in High Aspect Ratio Vessels (H/T>1.4H/T > 1.4): Leaving stagnant unmixed fluid in the upper third of the vessel.
  12. Ignoring Starting Torque: Selecting an undersized gearbox that shears gears when restarting an agitator submerged in settled solids.
  13. Ignoring Non-Newtonian Rheology: Assuming pseudoplastic or thixotropic slurries behave like water.
  14. Over-Agitating Crystallizations: Believing higher power always improves product; destroys crystal size distribution (CSD).
  15. Accepting Generic Vendor Guarantees: Accepting "this 5 KL agitator is standard" without verifying fluid properties, NjsN_{js}, and tmixt_{mix} requirements.

# 31. Standard Agitator Design Calculation Sheet

When generating process data sheets for equipment fabrication, process engineers must complete the following standard specification sheet:

================================================================================
          PHARMACEUTICAL BATCH REACTOR AGITATOR PROCESS DATASHEET
================================================================================
PROJECT / VESSEL TAG: R-101 (5 KL API Batch Synthesis Reactor)
--------------------------------------------------------------------------------
1. PROCESS OPERATING DATA:
   - Working Batch Volume (V):           4.0 m³ (4,000 L)
   - Operating Liquid Density (r\rho):     920 kg/m³
   - Operating Dynamic Viscosity (mu):   2.4 cP (0.0024 Pa·s)
   - Design Operating Temperature:       -20 °C to +140 °C
   - Design Operating Pressure:          Full Vacuum to +6.0 barg
   - Solids Content / Particle Size:     5.0 wt% Pd/C Catalyst / 35 µm

2. VESSEL INTERNAL GEOMETRY:
   - Tank Inside Diameter (T):           1,600 mm
   - Total Liquid Height (H):            1,890 mm (Aspect Ratio H/T = 1.18)
   - Bottom Head Geometry:               ASME Torispherical Dish Head
   - Baffles:                            4 Vertical Baffles, Width B = 160 mm (T/10)
   - Baffle Wall Clearance (B_gap):      32 mm (T/50)

3. AGITATOR MECHANICAL SPECIFICATION:
   - Mounting Type:                      Top-Entry Centered Cantilever
   - Lower Impeller:                     4-Blade Pitched Blade Turbine (45°), D1 = 560 mm
   - Upper Impeller:                     3-Blade High-Efficiency Hydrofoil, D2 = 560 mm
   - Bottom Clearance (C):               480 mm (C/T = 0.30)
   - Inter-Impeller Spacing (S):         700 mm (S/D = 1.25)
   - Impeller Assembly Power No (Np):    1.65 (Turbulent)
   - Impeller Pumping Number (Nq):       0.82

4. HYDRODYNAMIC CALCULATION SUMMARY:
   - Operating Agitator Speed (N):       218 RPM (3.63 rev/s) [VFD: 50 to 250 RPM]
   - Impeller Reynolds Number (Re):      4.36 × 10⁵ (Fully Turbulent)
   - Net Shaft Power Demand (P):         4.00 kW (5.36 HP)
   - Power per Unit Volume (P/V):        1.00 kW/m³
   - Impeller Tip Speed (v_tip):         6.39 m/s
   - Nominal Operating Torque (Tq):      175.4 N·m
   - 95% Blending Time (t_95):           22.1 seconds
   - Zwietering Njs for Pd/C:            88 RPM (Operating 218 RPM >> Njs = Full Suspension)

5. DRIVE & TRANSMISSION SPECIFICATION:
   - Electric Motor Rating:              7.5 kW (10 HP), 4-Pole, 415V/50Hz, Flameproof Ex-d
   - Gearbox Type & Ratio:               Helical-Bevel Unit, Ratio = 1:6.6
   - Solid Shaft Diameter (ds):          65 mm (SS316L Solid Forging)
   - First Critical Speed (N_crit):      410 RPM (Operating speed 218 RPM <= 0.53 N_crit)
   - Mechanical Seal:                    Double Balanced Cartridge Seal, Pressurized Barrier Pot
================================================================================

# 32. Vendor Data Verification Checklist

When reviewing bids from agitator vendors, process engineers must reject generic catalog statements ("This agitator fits a 5 KL reactor") and demand the following verified engineering deliverables:

  1. Certified Impeller Power & Pumping Curves (NpN_p and NQN_Q vs. ReRe): Verified across the vessel's specific D/TD/T, C/TC/T, and baffle geometry.
  2. Finite Element Shaft Deflection & Stress Analysis: Showing combined von Mises stresses \le allowable fatigue limits under peak starting torque.
  3. Modal Critical Speed Calculation Report: Proving first lateral natural frequency Ncrit1.40×NmaxN_{crit} \ge 1.40 \times N_{max}.
  4. Mechanical Seal Thermal Balance & Barrier Fluid Plan: Certifying barrier fluid flow rate and seal face temperature rise under full vacuum and maximum steam jacket temperatures.

# 33. The 11-Step Agitator Design Hierarchy

The systematic process engineering design workflow must always proceed sequentially through the following hierarchy:

[1. Process Objective Defined] (Blending, Solids, Heat Transfer, Gas-Liquid, Crystallization)
              |
              v
[2. Fluid Physical Properties] (Density, Viscosity Rheology, Solids Density/Size, Gas Flow)
              |
              v
[3. Vessel Geometry Established] (Tank Diameter T, Liquid Height H, Head Profile)
              |
              v
[4. Impeller Family Selected] (Axial Hydrofoil, Pitched Blade, Radial Rushton, High-Visc Ribbon)
              |
              v
[5. Geometric Ratios Fixed] (D/T = 0.33, C/T = 0.30, S/D = 1.25, Baffles B = T/10)
              |
              v
[6. Flow Regime & Reynolds Number Calculated] (Re = rho * N * D^2 / mu)
              |
              v
[7. Power Number & Power Demand Computed] (P = Np * rho * N^3 * D^5)
              |
              v
[8. Hydrodynamic Validation] (P/V, Tip Speed v_tip, Blend Time t_mix, Zwietering Njs)
              |
              v
[9. Scale-Up Protocol Selected] (Constant P/V, Tip Speed, or Njs)
              |
              v
[10. Shaft Mechanical & Dynamic Sizing] (Combined Torque/Bending, Diameter ds, Critical Speed Ncrit)
              |
              v
[11. Motor, Gearbox & Seal System Integrated] (Service Factor Fs >= 1.3, Double Mechanical Seal)

# 34. Core Engineering Conclusion

Agitator design is not about selecting a motor size or arbitrarily assigning an RPM from previous plant recipes.

It is about precisely tailoring three-dimensional fluid motion, shear distribution, and turbulent energy dissipation to satisfy the chemical and physical transport requirements of the process.

The most efficient agitator is never the one with the highest horsepower motor. It is the system that delivers the exact required process performance with the optimum balance of bulk volumetric circulation, localized shear, thermal boundary layer disruption, and multi-decade mechanical reliability.


# 35. Interactive Batch Reactor Scale-Up Calculator

Need to calculate power demand (PP), Reynolds number (ReRe), torque (TqT_q), tip speed (vtipv_{tip}), and scale-up parameters for your chemical or API reactor?

Launch the Interactive Batch Reactor Scale-Up Calculator →

Compute impeller power numbers, motor kilowatts, blend times, and multi-scale scale-up trajectories across laboratory, pilot, and commercial batch reactors.


# Applicable Engineering Standards & Codes Used

The engineering methodologies, design correlations, and safety criteria detailed in this article adhere to the following international standards and industry codes:

  • ASME Boiler and Pressure Vessel Code (BPVC) Section VIII Division 1 & 2: ASME Boiler and Pressure Vessel Code (BPVC) Section VIII Division 1 & 2
  • API 620 & API 650: Welded Tanks for Oil, Chemical and Liquid Storage
  • TEMA Class R, C & B: Tubular Exchanger Manufacturers Association Standards
  • DIN EN 13445: Unfired Pressure Vessels European Standard
  • IS 2825: Code for Unfired Pressure Vessels (Bureau of Indian Standards)
Reactor DesignAgitationScale-UpImpeller SelectionFluid DynamicsMixingPower NumberZwietering NjsAPI ManufacturingChemical Engineering
Comments (0)

Discussion

Please Log In to participate in the technical discussion.

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