Comprehensive Chemical Engineering Unit Operations: Distillation, Fluid-Solid Separations, Agitation, Heat Transfer, Reaction Engineering, and Process Intensification

Core Paradigms and Scales of Chemical Engineering

• Concept of Unit Operations:

  • Rationalizes and simplifies the study of chemical engineering equipment.
  • Equipment performing the same physical or chemical transformation is grouped under a single term (e.g., distillation, reaction, filtration, drying).
  • A complete production facility is conceptualized as a succession of unit operations, each analyzed with unified tools and principles.
  • Facilitates the modular design of new manufacturing processes.

• Core Paradigms of the Discipline:

  • Material and energy balances.
  • Coupling of mechanical, thermal, and chemical phenomena.
  • Multi-scale approach to process analysis.

• Spatial Scales in Chemical Engineering:

  • Macroscopic Scale: Human-level scale; phenomena observable with the naked eye or macroscopic imaging.
  • Mesoscopic Scale: Intermediate internal apparatus dimensions, without descending to micronic levels.
  • Microscopic Scale: Scale of thin boundary films and individual particles, typically on the order of $1\,\mu\text{m}$.
  • Molecular Scale: Nanometer level ($\text{nm}$) or smaller, governing molecular interactions and quantum phenomena.
  • Megascale / Metascale: Scale of manufacturing units, complete production sites, and corporate networks.

Distillation Principles and Vapor-Liquid Equilibria

• Fundamental Mechanism:

  • Separation of homogeneous mixtures (typically liquid) based on the formation of a vapor phase enriched in the most volatile component through coupled heat and mass transfer.
  • While energy-intensive, it presents minimal operational risk, making it the most widely deployed industrial separation operation.

• Laboratory Distillation Configuration:

  • Heat Source (1): Bunsen burner or heating mantle supplying vaporization energy.
  • Distillation Flask (2): Round-bottom vessel containing the boiling mixture.
  • Distillation Head (3): Directs ascending vapors to the condenser.
  • Thermometer (4): Positioned at the vapor takeoff to monitor boiling temperatures.
  • Refrigerant / Condenser (5): Double-walled cooling jacket operating with counter-current cold water (inlet 6, outlet 7).
  • Receiving Flask (8): Collects the condensed distillate.
  • Vacuum Adapter (9, 10): Connects the closed system to a vacuum pump to lower boiling temperatures when needed.

• Industrial Applications:

  • Petroleum Refining: Atmospheric and vacuum distillation producing gasoline, kerosene, diesel, and petrochemical feedstocks.
  • Cryogenic Air Separation: Liquefaction and fractionation of air to produce pure $\text{O}_2$ and $\text{N}_2$.
  • Alcoholic Spirits Production: Separation of ethanol from fermented agricultural broths.

• Process Classifications:

  • Flash Distillation: Single-stage equilibrium separation, analogous to simple batch distillation.
  • Rectification: Multi-stage continuous or batch distillation featuring internal vapor-liquid counter-current contact.
  • Binary Distillation: Separation restricted to two distinct chemical components.

• Operational Criteria for Binary Distillation:

  • Relative volatility $\alpha_{1,2}$ must be sufficiently large:     $\alpha_{1,2} \approx \frac{P_1^{\text{sat}}}{P_2^{\text{sat}}} \ge 1.2$     where $P_1^{\text{sat}}$ is the saturation vapor pressure of the light (more volatile) component and $P_2^{\text{sat}}$ is that of the heavy (less volatile) component.
  • Industrial continuous flow rates typically exceed $100\,\text{kg\,h}^{-1}$.
  • Components must be chemically and thermally stable at process temperatures with no explosion hazards.

• Thermodynamic Vapor-Liquid Equilibrium (VLE):

  • Isobaric Diagram ($T$ vs. $x, y$):
  • Bubble Point Curve: Lower boundary indicating the temperature at which the first vapor bubble forms during heating.
  • Dew Point Curve: Upper boundary indicating the temperature at which the first liquid drop condenses during cooling.
  • Two-Phase Region: Intermediate domain between the bubble and dew curves where liquid and vapor coexist.
  • Equilibrium Curve ($y$ vs. $x$):
  • Relates vapor mole fraction $y$ to liquid mole fraction $x$ at thermodynamic equilibrium.
  • Anchored at coordinates $(0,0)$ and $(1,1)$.
  • For ideal mixtures, the curve lies entirely above the first diagonal ($y = x$).
  • Distillation of a Zeotropic Liquid of Composition $x_1$:
  • Initial heating reaches the bubble point at temperature $T_1$, generating the first vapor bubble of composition $y_1 > x_1$.
  • Continued heating into the two-phase zone at temperature $T_2$ establishes an equilibrium between liquid of composition $x_2$ and vapor of composition $y_2$, where $x_2 < x_1 < y_2 < y_1$.
  • Reaching the dew point temperature $T_3$ vaporizes the final drop of liquid, yielding total vapor of composition $y_3 = x_1$.
  • Azeotropic Behavior:
  • Occurs when bubble and dew point curves contact tangentially at a pinch point.
  • The equilibrium curve intersects the diagonal $y = x$, preventing further enrichment beyond the azeotropic composition by conventional distillation (e.g., ethanol-water azeotrope at $96\%$ by mass).

Distillation Column Balances and the McCabe-Thiele Method

• Column Architecture and Mass Exchange:

  • Liquid flows downward by gravity from top to bottom.
  • Vapor generated in the reboiler (boiler) ascends counter-currently from bottom to top.
  • Rectifying (Enriching) Section: Located above the feed tray; liquid reflux scrubs out the heavy component, enriching vapor in the light component.
  • Stripping Section: Located below the feed tray; ascending vapor strips the light component from descending liquid, concentrating the heavy component at the bottom.
  • Reflux Ratio ($R$):$R = \frac{L_0}{D}$     where $L_0$ is the liquid reflux rate returned to the column top and $D$ is the distillate product rate.

• Macroscopic Conservation Balances:

  • Total Material Balance:$A = B + D$
  • Light Component Material Balance:$A \cdot z_A = B \cdot x_B + D \cdot x_D$
  • Enthalpy Balance:$A \cdot h_A + \phi_R = D \cdot h_D + B \cdot h_B + \phi_C$
  • Nomenclature:
  • $A, B, D$: Total mass or molar flow rates of feed, bottoms residue, and distillate.
  • $z_A, x_B, x_D$: Mass or mole fractions of the light component in feed, bottoms, and distillate.
  • $h_A, h_B, h_D$: Specific or molar enthalpies of feed, bottoms, and distillate.
  • $\phi_R, \phi_C$: Thermal power duties supplied to the reboiler and removed by the condenser.

• McCabe-Thiele Method Assumptions (Lewis Hypotheses):

  • The column operates adiabatically (perfectly insulated, heat loss through walls is zero).
  • Heat of mixing of the two components is negligible.
  • Latent heats of vaporization of both components are nearly identical, and sensible heat changes are negligible compared to latent heat.
  • Constant Molar Overflow (CMO): Liquid and vapor molar flow rates remain constant within each section of the column, yielding straight operating lines on a McCabe-Thiele plot.

• Step-by-Step McCabe-Thiele Construction:

  • Step 1: Plot Reference Points: On the $y = x$ diagonal line, plot the distillate point $(x_D, x_D)$, feed point $(z_A, z_A)$, and residue point $(x_B, x_B)$.
  • Step 2: Rectifying Operating Line: Draw the line passing through $(x_D, x_D)$ with slope:$\text{Slope} = \frac{R}{R + 1}$
  • Step 3: Feed Line ($q$-line): Draw the feed condition line through $(z_A, z_A)$. The slope depends on the thermal quality of the feed:
  • Subcooled liquid: slope $> 1$
  • Saturated (boiling) liquid: vertical line ($\text{slope} = \infty$)
  • Partially vaporized mixture: negative slope
  • Saturated vapor: horizontal line ($\text{slope} = 0$)
  • Superheated vapor: positive slope between $0$ and $1$
  • Step 4: Stripping Operating Line: Connect the residue point $(x_B, x_B)$ to the intersection point of the rectifying operating line and the feed line.
  • Step 5: Step Off Stages: Construct horizontal and vertical steps between the operating lines and the equilibrium curve, beginning at $(x_D, x_D)$. The total number of steps represents the Number of Equilibrium Stages ($\text{NES}$). The step crossing the intersection of the two operating lines identifies the optimal feed tray location.

• Limiting Reflux Conditions:

  • Total Reflux ($R \rightarrow \infty$):
  • No distillate is withdrawn ($D = 0$); all condensed overhead vapor is returned as reflux ($L_0 = V$).
  • The operating lines merge with the $y = x$ diagonal.
  • Yields the minimum number of theoretical stages ($\text{NES}_{\min}$).
  • Minimum Reflux ($R_{\min}$):
  • The rectifying operating line intersects the equilibrium curve at the feed line junction (pinch point).
  • Requires an infinite number of stages ($\text{NES} \rightarrow \infty$) to achieve the separation.
  • Industrial Operating Range: Typical economic optimum is set to $R = 1.2 \cdot R_{\min}$ to $1.5 \cdot R_{\min}$.

• Alternative Stage Determination Methods:

  • Ponchon-Savarit Method: Rigorous graphical method using enthalpy-concentration coordinates, avoiding the constant molar overflow assumption.
  • Commercial Simulators: Multi-component rigorous numerical tray-by-tray modeling based on MESH equations and thermodynamic activity/fugacity databases.

Column Internals and Sizing Methodologies

• Packed Columns:

  • Vertical cylindrical shells filled with packing material to promote contact between downflowing liquid and rising vapor.
  • Random (Dumped) Packings: Disordered small units (Raschig rings, Pall rings, saddles) providing low cost and high specific area; susceptible to compaction, fouling, and unpredictable channeling.
  • Structured Packings: Arranged corrugated sheets of metal gauze or perforated sheet; provide exceptionally low pressure drops, high void fractions, and predictable mass transfer; significantly higher capital cost.
  • Column Internals: Require liquid distributors at the top, support grids, and liquid redistributors every few meters to counter wall channeling.
  • Height Equivalent to a Theoretical Plate (HETP):$\text{Total Height of Packing} = \text{NES} \cdot \text{HETP}$     where HETP is determined empirically from packing geometry and system fluid properties.

• Tray (Plate) Columns:

  • Staged vessels with horizontal trays spaced regularly along the column height.
  • Sieve Trays: Perforated flat plates; vapor bubbles through holes while liquid crosses the tray horizontally.
  • Valve Trays: Perforations covered by liftable caps/valves that open with rising vapor flow, offering wide turndown ratios.
  • Bubble-Cap Trays: Fixed risers capped with slotted hoods; maintain liquid seal at all flow rates, eliminating weeping at very low vapor velocities.
  • Downcomers: Channels that convey liquid from an upper tray to the tray below.
  • Actual Stages: Calculated by dividing $\text{NES}$ by overall tray efficiency ($E_0$).

• Preliminary Column Design Protocol:

  • 1) Graphically establish $\text{NES}_{\min}$ at total reflux using McCabe-Thiele.
  • 2) Graphically determine $R_{\min}$ at the pinch point.
  • 3) Select operating reflux: $R = 1.2 \cdot R_{\min}$ to $1.5 \cdot R_{\min}$.
  • 4) Step off the operational $\text{NES}$ on the McCabe-Thiele diagram.
  • 5) Calculate column height via HETP (for packed columns) or tray efficiency and tray spacing (for tray columns), followed by hydraulic diameter sizing to avoid flooding.

Fluid-Solid Mechanical Separations: Fundamentals and Hydrodynamics

• Operational Regimes by Particle Size:

  • Clarification/Decantation: Low solids concentrations with fine particles ($1$ to $10\,\mu\text{m}$).
  • Centrifugation: Fine to medium particles ($10$ to $500\,\mu\text{m}$) in medium-concentration suspensions.
  • Filtration: Concentrated slurries containing broad particle distributions ($5$ to $1,000\,\mu\text{m}$).

• Hydrodynamics of Fluid Flow Around Isolated Particles:

  • Particle Reynolds Number ($\text{Re}_p$):$\text{Re}_p = \frac{\rho_f \cdot u \cdot d_p}{\mu}$     where $\rho_f$ is fluid density, $\mu$ is fluid dynamic viscosity, $d_p$ is particle diameter, and $u$ is relative fluid-particle velocity.
  • Total Drag Force ($T$):$T = \frac{1}{2} C_x \cdot \rho_f \cdot u^2 \cdot \Omega_c$     where $C_x$ is the drag coefficient and $\Omega_c$ is the projected area of the particle perpendicular to flow.
  • Archimedes Number ($\text{Ar}$):$\text{Ar} = \frac{d_p^3 \cdot g \cdot \rho_f \cdot (\rho_p - \rho_f)}{\mu^2}$     where $\rho_p$ is solid particle density and $g$ is gravitational acceleration.

• Terminal Fall Velocities of an Isolated Particle ($u_{t0}$):

  • Laminar Regime (Stokes's Law):
  • Criterion: $\text{Ar} < 27.6$
  • Reynolds relationship: $\text{Re}_{t0} = \frac{\text{Ar}}{18}$
  • Settling velocity:       $u_{t0} = \frac{d_p^2 \cdot g \cdot (\rho_p - \rho_f)}{18 \cdot \mu}$
  • Intermediate Regime (Van Allen's Law):
  • Criterion: $27.6 < \text{Ar} < 4.4 \times 10^4$
  • Reynolds relationship: $\text{Re}_{t0} = 0.153 \cdot \text{Ar}^{0.714}$
  • Turbulent Regime (Newton's Law):
  • Criterion: $4.4 \times 10^4 < \text{Ar} < 1.1 \times 10^9$
  • Reynolds relationship: $\text{Re}_{t0} = \sqrt{3 \cdot \text{Ar}}$
  • Settling velocity:       $u_{t0} = \sqrt{\frac{3 \cdot d_p \cdot g \cdot (\rho_p - \rho_f)}{\rho_f}}$

• Hindered Settling of Particle Swarms:

  • In concentrated suspensions, the fall velocity $u_t$ decreases relative to $u_{t0}$:     $u_t = u_{t0} \cdot (1 - \varphi)^n$     where $\varphi$ is the solid volume fraction in suspension.
  • Richardson-Zaki Exponent ($n$):
  • $n = 4.65$ for $\text{Re}_{t0} < 0.2$
  • $n = 4.4 \cdot \text{Re}_{t0}^{-0.03}$ for $0.2 < \text{Re}_{t0} < 1$
  • $n = 4.4 \cdot \text{Re}_{t0}^{-0.1}$ for $1 < \text{Re}_{t0} < 500$
  • $n = 2.4$ for $\text{Re}_{t0} > 500$
  • Suspension Bulk Properties:
  • Density: $\rho = \varphi \cdot \rho_p + (1 - \varphi) \cdot \rho_f$
  • Viscosity: $\mu = \mu_f \cdot 10^{1.82 \cdot \varphi}$

• Centrifugation Velocity:

  • Centrifugal acceleration replaces gravity: $\omega^2 \cdot r = K \cdot g$, where $K = \frac{\omega^2 \cdot r}{g}$.
  • Laminar regime: $u_c = K \cdot u_t$
  • Intermediate regime: $u_c = K^{2/3} \cdot u_t$
  • Turbulent regime: $u_c = \sqrt{K} \cdot u_t$

• Flow Through Porous Media:

  • Darcy's Law:$\Delta P = \frac{\mu \cdot L}{B \cdot \Omega} \cdot Q_v$     where $B$ is bed permeability (measured in Darcy, where $1\,\text{Darcy} = 10^{-12}\,\text{m}^2$; sand $\approx 100\,\text{Darcy}$, cigarette $\approx 1,000\,\text{Darcy}$, soil $\approx 1\,\text{Darcy}$, limestone $\approx 10^{-2}\,\text{Darcy}$).
  • Ergun's Equation:$\frac{\Delta P}{L} = K_h \cdot \frac{(1 - \varepsilon)^2}{\varepsilon^3} \cdot a_p^2 \cdot \mu \cdot u_s + h_B \cdot \frac{1 - \varepsilon}{\varepsilon^3} \cdot a_p \cdot \rho_f \cdot u_s^2$
  • First term represents the Carman-Kozeny viscous loss equation ($K_h \approx 4.2$).
  • Second term represents the Burke-Plummer inertial loss equation ($h_B \approx 0.3$).
  • $\varepsilon$: Bed porosity (void fraction).
  • $a_p$: Specific surface area per unit volume of particle ($a_p = \frac{6}{d_p}$ for spheres).
  • Superficial/surface velocity: $u_s = \frac{Q_v}{\Omega}$.
  • Interstitial average pore velocity: $u_p = \frac{u_s}{\varepsilon}$.

Decantation, Cyclones, and Filtration Technologies

• Gravity Decantation (Settling):

  • Uses settling basins or clarifiers (rectangular, cylindrical, or cylindroconical) equipped with bottom mechanical rakes.
  • Single-Particle Criterion: Separation occurs if particle fall time is shorter than tank residence time:     $t_{\text{fall}} = \frac{H}{u_{t0}} < \tau = \frac{V}{Q_v}$
  • Hydraulic capacity of a conventional clarifier equals ground surface area multiplied by particle terminal settling velocity.
  • Lamellar Decanters: Incorporate $n$ inclined plates of dimensions $L \times l$ tilted at angle $\alpha$ to the horizontal:     $Q_{v,\max} = n \cdot L \cdot l \cdot \cos(\alpha) \cdot u_{t0}$
  • Kynch theory is applied for concentrated, zone-settling sludges.

• Centrifugal Decantation:

  • Preferred when density differences are small ($\Delta \rho / \rho_f > 10\%$), kinematic viscosity is low ($< 1.5 \times 10^{-4}\,\text{m}^2\,\text{s}^{-1}$), or particles are small ($> 0.5\,\mu\text{m}$).
  • Tubular Centrifuges: High-speed units with centrifugal acceleration factors $K > 3,000$.
  • Disk-Stack Centrifuges: Bowls measuring $0.15$ to $1.0\,\text{m}$ in diameter rotating up to $12,000\,\text{rpm}$; handle suspensions up to $15\%$ solids by volume with flow rates up to $100\,\text{m}^3\,\text{h}^{-1}$.
  • Sigma Theory for Decanter Centrifuges:$Q_v = \Sigma \cdot u_t$$\Sigma = \frac{\omega^2 \cdot V}{g \cdot \ln\left(\frac{R}{R_0}\right)}$     where $V$ is the liquid volume, $R_0$ is the inner liquid surface radius, and $R$ is the outer bowl wall radius.

• Cyclone and Hydrocyclone Separators:

  • Use tangential fluid feed entries (velocities of tens of meters per second) converting kinetic energy to a vortex, forcing solids outward to the walls while clarified fluid discharges axially.
  • Remove particles from $1$ to $100\,\mu\text{m}$.
  • Have no moving parts; operate up to $1,000^\circ\text{C}$ and pressures up to $500\,\text{bar}$.

• Cake Filtration Theory:

  • Filtration Flux Equation:$\frac{dV}{dt} = u \cdot \Omega$
  • Pressure Drop Model (Darcy Basis):$\Delta P = R \cdot \mu \cdot u = (R_s + R_g) \cdot \mu \cdot u$     where $R_s$ is the clean filter medium resistance ($\text{m}^{-1}$) and $R_g$ is the accumulated cake resistance ($\text{m}^{-1}$).
  • Specific Cake Resistance ($\alpha$ in $\text{m\,kg}^{-1}$):$R_g = \frac{\alpha \cdot M}{\Omega}$     where $M$ is the mass of dry cake solids.     $\alpha = \frac{1 - \varepsilon}{\rho_s \cdot B_0}$
  • Cake Thickness ($Z$):$Z = \frac{w \cdot V}{\Omega \cdot (1 - \varepsilon) \cdot \rho_s}$     where $w = \frac{M}{V}$.
  • Suspension Moisture Relationships:$m = 1 + \frac{\varepsilon \cdot \rho_{\text{liq}}}{(1 - \varepsilon) \cdot \rho_s}$$w = \frac{s \cdot \rho_{\text{liq}}}{1 - m \cdot s}$     where $m$ is wet cake mass divided by dry cake mass, and $s$ is suspension solid mass fraction.
  • General Differential Equation of Filtration:$\Delta P = \mu \cdot \left(\frac{R_s}{\Omega} + \frac{\alpha \cdot w \cdot V}{\Omega^2}\right) \cdot \frac{dV}{dt}$

• Analytical Integrations of the Filtration Equation:

  • Constant-Flow Filtration ($u = \frac{dV}{\Omega \cdot dt} = \text{const}$):$\Delta P = R_s \cdot \mu \cdot u + \alpha \cdot \mu \cdot w \cdot u^2 \cdot t$
  • Constant-Pressure Filtration ($\Delta P = \text{const}$):$\frac{t}{V} = \frac{R_s \cdot \mu}{\Omega \cdot \Delta P} + \frac{\alpha \cdot \mu \cdot w}{2 \cdot \Omega^2 \cdot \Delta P} \cdot V$     Plotting $t/V$ versus $V$ yields a straight line where the slope provides $\alpha$ and the intercept provides $R_s$.
  • Compressible Cakes: Specific resistance increases with pressure drop across the cake:     $\alpha = \alpha_0 \cdot (\Delta P_g)^n$     with compressibility index $0 < n < 0.8$ (typically $n \approx 0.1$).

Agitation and Mixing (Stirring)

• Objectives of Stirring in Chemical Engineering:

  • Miscible fluid homogenization (blending).
  • Suspension of solid particles (dissolution, catalytic slurry reactors).
  • Liquid-liquid dispersion (emulsions, liquid extraction).
  • Gas-liquid dispersion (gas absorption, bioreactors).

• Flow Regimes and Patterns Generated by Impellers:

  • Axial Flow: Pumping parallel to the shaft axis; creates a single circulation loop across the vessel height (e.g., marine propellers, hydrofoil impellers).
  • Radial Flow: Pumping perpendicular to the shaft toward vessel walls; splits into two loops (upper and lower) with high shear near the blades (e.g., Rushton flat-blade turbines).
  • Tangential Flow: Liquid rotates in circles along the wall; creates poor bulk vertical turnover unless scraped (e.g., anchors, helical ribbons for high-viscosity media).

• Vortex Dynamics and Baffling:

  • Central unbaffled rotation pulls down a deep surface vortex, leading to motor overload, mechanical instability, and undesirable gas ingestion.
  • Remedies: Off-centering the shaft, tilting the shaft, or installing vertical wall baffles.
  • Standard Baffling: Four vertical plates set perpendicular to the wall, each with width:$b = \frac{T}{10}$

• Standard Tank Geometry Proportions:

  • Liquid height: $L_z \approx T$ (tank diameter).
  • Impeller diameter: $d_A = \frac{T}{3}$.
  • Impeller off-bottom clearance: $z_A = \frac{L_z}{3}$.

• Dimensionless Characterization in Stirring:

  • Impeller Reynolds Number ($\text{Re}$):$\text{Re} = \frac{\rho \cdot N \cdot d_A^2}{\mu}$     where $N$ is rotational frequency ($\text{s}^{-1}$) and $d_A$ is impeller diameter ($\text{m}$).
  • Power Number ($N_p$):$N_p = \frac{P}{\rho \cdot N^3 \cdot d_A^5}$     where $P$ is power dissipated into fluid ($\text{W}$).
  • Pumping Number ($N_{QP}$):$N_{QP} = \frac{Q_p}{N \cdot d_A^3}$     where $Q_p$ is the volumetric flow rate pumped by the impeller ($\text{m}^3\,\text{s}^{-1}$).

• Power Number Regimes on Log-Log Coordinates:

  • Laminar Regime ($\text{Re} < 10$): Viscous forces dominate; power number is inversely proportional to Reynolds number:     $N_p \propto \frac{1}{\text{Re}}$
  • Turbulent Regime ($\text{Re} > 10^4$ in baffled vessels): Inertial forces dominate; $N_p$ reaches a constant plateau value $N_{p0}$.

• Impeller Performance Parameters in the Turbulent Regime:

• Impeller Selection Rules:

  • Turbulent Homogenization: High pumping, low shear, low power requirement $\rightarrow$ marine propellers, hydrofoil blades.
  • Viscous Homogenization ($\mu > 10\,\text{Pa\,s}$): High wall clearance, bulk positive displacement $\rightarrow$ helical ribbons, anchors, screw impellers.
  • Solid Suspension: High axial flow, moderate shear $\rightarrow$ hydrofoil, inclined blade propellers; shift to turbines if mass transfer is limiting.
  • Gas Dispersion: High shear to break bubbles, high power $\rightarrow$ Rushton turbines, rotor-stator devices.

• Stirred Tank Scale-Up Rules:

  • Rule 1: Constant Power per Unit Volume ($P/V = \text{const}$):     Keeps fluid micro-mixing uniform across scales:     $N_1^3 \cdot d_{A1}^2 = N_2^3 \cdot d_{A2}^2$     (in the turbulent regime with identical fluids and geometrical similarity).
  • Rule 2: Constant Blade Tip Velocity ($u_{\text{tip}} = \pi \cdot N \cdot d_A = \text{const}$):     Preserves shear stress distribution and pump efficiency:     $N_1 \cdot d_{A1} = N_2 \cdot d_{A2}$

Heat Exchangers and Thermal Design

• Industrial Roles: Feed preheating, reboiling/evaporation, reactor effluent cooling, vapor condensation, process heat recovery.

• Major Industrial Technologies:

  • Tubular Heat Exchangers:
  • Coaxial Double-Pipe: Simple inner pipe surrounded by concentric outer pipe; runs counter-currently or co-currently.
  • Shell and Tube (TEMA Standards): Tube bundle enclosed within an outer shell. TEMA nomenclature standardizes front heads (stationary channels), shell styles (E, F, G, H, J, K, X), and rear heads (fixed, floating, U-tube).
  • Plate Heat Exchangers:
  • Series of corrugated metal plates clamped together; corrugated channels generate turbulence at low Reynolds numbers, yielding high heat transfer coefficients and compact footprints.
  • Types: gasketed (limited to moderate temperatures and pressures; susceptible to chemical attack) vs. welded/brazed (handle higher pressures/temperatures; difficult to mechanically clean).

• Overall Thermal Energy Balances (Without Ambient Heat Losses):$\phi = Q_{mc} \cdot Cp_c \cdot (T_{ce} - T_{cs}) = Q_{mf} \cdot Cp_f \cdot (T_{fs} - T_{fe})$   where $Q_{mc}, Q_{mf}$ are mass flow rates of hot and cold fluids ($\text{kg\,s}^{-1}$), $Cp_c, Cp_f$ are specific heat capacities ($\text{J\,kg}^{-1}\,\text{K}^{-1}$), and indices $e$ and $s$ indicate inlet and outlet.

• General Heat Transfer Rate Equation:$\phi = H \cdot S \cdot \Delta T_{ml}$

  • $H$: Overall heat transfer coefficient ($\text{W\,m}^{-2}\,\text{K}^{-1}$, accounting for internal/external convection and wall conduction).
  • $S$: Heat transfer surface area ($\text{m}^2$).
  • $\Delta T_{ml}$: Logarithmic Mean Temperature Difference (LMTD).

• Derivation of the Logarithmic Mean Temperature Difference (LMTD):

  • Assumptions: Steady-state operation; perfect thermal insulation ($\phi_{\text{loss}} = 0$); negligible axial heat conduction; constant fluid physical properties ($Cp, \mu$); uniform heat transfer coefficient $H$.
  • Energy Balance on Differential Area $dS$:$d\phi = H \cdot (T_c - T_f) \cdot dS$
  • For Co-Current Flow:
  • $d\phi = -Q_{mc} \cdot Cp_c \cdot dT_c = Q_{mf} \cdot Cp_f \cdot dT_f$
  • $d(T_c - T_f) = dT_c - dT_f = -d\phi \cdot \left(\frac{1}{Q_{mc} \cdot Cp_c} + \frac{1}{Q_{mf} \cdot Cp_f}\right)$
  • Substituting $d\phi$:       $\frac{d(T_c - T_f)}{T_c - T_f} = -H \cdot \left(\frac{1}{Q_{mc} \cdot Cp_c} + \frac{1}{Q_{mf} \cdot Cp_f}\right) \cdot dS$
  • Integrating from $S = 0$ to $S$:       $\ln\left(\frac{T_{cs} - T_{fs}}{T_{ce} - T_{fe}}\right) = -H \cdot S \cdot \left(\frac{1}{Q_{mc} \cdot Cp_c} + \frac{1}{Q_{mf} \cdot Cp_f}\right)$
  • Since $\frac{1}{Q_{mc} \cdot Cp_c} = \frac{T_{ce} - T_{cs}}{\phi}$ and $\frac{1}{Q_{mf} \cdot Cp_f} = \frac{T_{fs} - T_{fe}}{\phi}$:       $\ln\left(\frac{T_{cs} - T_{fs}}{T_{ce} - T_{fe}}\right) = -\frac{H \cdot S}{\phi} \cdot [(T_{ce} - T_{fe}) - (T_{cs} - T_{fs})]$
  • Rearranging:       $\phi = H \cdot S \cdot \frac{(T_{ce} - T_{fe}) - (T_{cs} - T_{fs})}{\ln\left(\frac{T_{ce} - T_{fe}}{T_{cs} - T_{fs}}\right)}$
  • For Counter-Current Flow:
  • $d\phi = -Q_{mc} \cdot Cp_c \cdot dT_c = -Q_{mf} \cdot Cp_f \cdot dT_f$
  • Yields the terminal driving force expression:       $\phi = H \cdot S \cdot \frac{(T_{ce} - T_{fs}) - (T_{cs} - T_{fe})}{\ln\left(\frac{T_{ce} - T_{fs}}{T_{cs} - T_{fe}}\right)}$
  • General Formula:$\Delta T_{ml} = \frac{\Delta T_1 - \Delta T_2}{\ln\left(\frac{\Delta T_1}{\Delta T_2}\right)}$     where $\Delta T_1$ and $\Delta T_2$ are temperature differences between streams at the two ends of the exchanger.

• Worked Digital Sizing Application:

  • Parameters:
  • Hot product stream: $Q_{mc} = 29.6\,\text{kg\,s}^{-1}$, $T_{ce} = 420\,\text{K}$, $T_{cs} = 380\,\text{K}$, $Cp_c = 2,200\,\text{J\,kg}^{-1}\,\text{K}^{-1}$.
  • Cold crude oil stream: $T_{fe} = 295\,\text{K}$, $T_{fs} = 330\,\text{K}$, $Cp_f = 1,986\,\text{J\,kg}^{-1}\,\text{K}^{-1}$.
  • Heat transfer coefficient: $H = 433.8\,\text{W\,m}^{-2}\,\text{K}^{-1}$.
  • Step 1: Total Heat Exchanged ($\phi$):$\phi = Q_{mc} \cdot Cp_c \cdot (T_{ce} - T_{cs}) = 29.6 \times 2200 \times (420 - 380) = 2,604,800\,\text{W}$
  • Step 2: Cold Fluid Flow Rate ($Q_{mf}$):$Q_{mf} = \frac{\phi}{Cp_f \cdot (T_{fs} - T_{fe})} = \frac{2,604,800}{1986 \times (330 - 295)} = 37.47\,\text{kg\,s}^{-1}$
  • Step 3: Counter-Current LMTD ($\Delta T_{ml}$):
  • $\Delta T_1 = T_{ce} - T_{fs} = 420 - 330 = 90\,\text{K}$
  • $\Delta T_2 = T_{cs} - T_{fe} = 380 - 295 = 85\,\text{K}$$\Delta T_{ml} = \frac{90 - 85}{\ln(90/85)} = 87.48\,\text{K}$
  • Step 4: Required Heat Transfer Surface Area ($S$):$S = \frac{\phi}{H \cdot \Delta T_{ml}} = \frac{2,604,800}{433.8 \times 87.48} = 68.6\,\text{m}^2$
  • Co-Current Comparison:
  • $\Delta T_1 = 420 - 295 = 125\,\text{K}$
  • $\Delta T_2 = 380 - 330 = 50\,\text{K}$$\Delta T_{ml,\text{co}} = \frac{125 - 50}{\ln(125/50)} = 81.85\,\text{K}$     Counter-current operation delivers a higher mean temperature driving force ($87.48\,\text{K}$ vs. $81.85\,\text{K}$), reducing required surface area.

• Multipass and Crossflow Calculation Methods:

  • Correction Factor Method:$\phi = F \cdot H \cdot S \cdot \Delta T_{ml,\text{counter}}$     where $F$ ($\le 1$) is read from chart nomograms based on temperature efficiency $eff$ and capacity ratio $R$:     $eff = \frac{T_{fs} - T_{fe}}{T_{ce} - T_{fe}}$$R = \frac{T_{ce} - T_{cs}}{T_{fs} - T_{fe}}$
  • Number of Transfer Units ($\text{eff-NUT}$) Method:$\phi = eff \cdot (Q_m \cdot Cp)_{\min} \cdot (T_{ce} - T_{fe})$$\text{NUT} = \frac{H \cdot S}{(Q_m \cdot Cp)_{\min}}$$R = \frac{(Q_m \cdot Cp)_{\min}}{(Q_m \cdot Cp)_{\max}}$

Chemical Reaction Engineering: Ideal and Non-Ideal Reactors

• Classification of Chemical Reactions and Reactors:

  • Phase Status: Homogeneous (single fluid phase) vs. Heterogeneous (multi-phase: solid-catalyzed, gas-liquid, fluid-solid consumable).
  • Thermodynamic Attributes: Complete vs. Reversible; Endothermic vs. Exothermic vs. Athermic.
  • Regime Modes: Transient/unsteady vs. Stationary/steady-state.
  • Operational Modes: Batch (closed), Semi-batch (semi-open), Continuous (open).

• Reaction Metrics and Progress Quantities:

  • Fractional Conversion of Key Reactant A ($X_A$):
  • Closed Reactor: $X_A = 1 - \frac{n_A}{n_{A0}}$
  • Open Reactor: $X_A = 1 - \frac{F_A}{F_{A0}}$
  • Moles / Flux of Species $j$ Related to Conversion ($\nu_A = -1$):
  • Closed: $n_j = n_{j0} - \frac{\nu_j}{\nu_A} \cdot n_{A0} \cdot X_A$
  • Open: $F_j = F_{j0} - \frac{\nu_j}{\nu_A} \cdot F_{A0} \cdot X_A$
  • Generalized Extent of Reaction ($X_i$) for Multiple Reactions:
  • Closed: $n_j = n_{j0} + n_0 \cdot \sum_i \nu_{ij} \cdot X_i$
  • Open: $F_j = F_{j0} + F_0 \cdot \sum_i \nu_{ij} \cdot X_i$
  • Total Moles (Closed): $n_{\text{total}} = n_I + n_0 \cdot \left(1 + \sum_i \Delta\nu_i \cdot X_i\right)$
  • Total Molar Flux (Open): $F_{\text{total}} = F_I + F_0 \cdot \left(1 + \sum_i \Delta\nu_i \cdot X_i\right)$
  • Stoichiometric Dilation: $\Delta\nu_i = \sum_j \nu_{ij}$

• Material Balances for Ideal Reactor Archetypes:

  • General Material Balance on Species $j$:$F_j^e + \left(\sum_i \nu_{ij} \cdot r_i\right) \cdot V = F_j^s + \frac{dn_j}{dt}$
  • Closed Ideally Stirred Batch Reactor:$\nu_j \cdot r \cdot V = \frac{dn_j}{dt}$
  • Continuous Stirred-Tank Reactor (CSTR / Open ISTR at Steady State):$F_j^e + \nu_j \cdot r \cdot V = F_j^s$     Since mixing is instantaneous and complete, exit composition and temperature match the uniform interior conditions.
  • Plug Flow Reactor (PFR at Steady State):$dF_j = \nu_j \cdot r \cdot dV$     Fluid elements pass through with uniform velocity profile and no axial mixing.

• Reactor Volume Optimization via the Levenspiel Plot ($C_{A0}/r$ vs. $X_A$):

  • PFR Volume: Proportional to the integral area under the curve between inlet and outlet conversion:     $V_{\text{PFR}} = F_{A0} \cdot \int_0^{X_A} \frac{dX_A}{r}$
  • CSTR Volume: Proportional to the rectangle area evaluated at the outlet conversion:     $V_{\text{CSTR}} = F_{A0} \cdot \frac{X_A}{r(X_A)}$
  • Comparison Rules:
  • Monotonically increasing curve (positive order kinetics) $\rightarrow$ PFR requires smaller volume than CSTR.
  • Monotonically decreasing curve (autocatalytic kinetics) $\rightarrow$ CSTR requires smaller volume than PFR.
  • Curve displaying a local minimum $\rightarrow$ Optimal volume is achieved by a CSTR operating up to the conversion of the minimum ($X_{A,\min}$), followed by a PFR up to the final target conversion ($X_{As}$).

• Selectivity and Yield in Complex Networks:

  • Reaction Types: Main reaction (yields product of interest P); Parallel/competing reactions (consume reactant without producing P); Serial/consecutive reactions (consume product P).
  • Overall Yield ($Y_{P/A}$):$Y_{P/A} = \frac{F_P}{\nu_{P/A} \cdot F_{A0}}$
  • Local Relative Efficiency ($\eta'_{P/A}$):$\eta'_{P/A} = \frac{R_P}{-\nu_{P/A} \cdot R_A}$     where $R_j = \sum_i \nu_{ij} \cdot r_i$.
  • Overall yield represents the integral of local relative efficiency over conversion:     $Y_{P/A} = \int_0^{X_{As}} \eta'_{P/A}\,dX_A$
  • If $\eta'_{P/A}$ decreases with $X_A$, a PFR maximizes yield; if it increases with $X_A$, a CSTR maximizes yield.

• Thermal Balance and Operation of Reactors:

  • Steady-State Energy Balance:$\phi = Q_m^e \cdot Cp \cdot (T_s - T_e) + F_0 \cdot \sum_i \Delta_r H_i(T_s) \cdot X_i$
  • Endothermic or Irreversible Exothermic Reactions: Operate isothermally at the maximum safe allowable temperature to maximize reaction rate and conversion.
  • Reversible Exothermic Reactions: Exhibit an Optimum Temperature Progression (OTP). High initial temperatures accelerate kinetics; progressively lower temperatures preserve favorable thermodynamic equilibrium at higher conversions. Industrially achieved using a cascade of adiabatic beds with intermediate gas coolers.

• Residence Time Distribution (RTD) in Non-Ideal Reactors:

  • Used to diagnose macroscopic flow anomalies (short-circuiting, bypassing, recirculation, stagnant dead zones).
  • Experimental Protocol: Inject an inert tracer (Dirac pulse or Heaviside step) into the feed stream and monitor its exit concentration $C(t)$.
  • Residence Time Distribution Function ($E(t)$):$E(t) = \frac{C(t)}{\int_0^\infty C(t)\,dt}$
  • Cumulative Distribution Function ($F(t)$):$F(t) = \int_0^t E(t)\,dt$
  • Two-Point Convolution for Inaccessible Inlets:$C_s(t) = \int_0^t C_e(t - t_s) \cdot E(t_s)\,dt_s$
  • Hydrodynamic Diagnostic Rules:
  • Ideal PFR exhibits a single Dirac delta peak shifted by passage time $\tau$.
  • Ideal CSTR exhibits an immediate jump followed by exponential decay: $E(t) = \frac{1}{\tau} \cdot e^{-t/\tau}$.
  • Mean residence time $\bar{t}_s < \tau$ indicates stagnant dead zones or channeling.
  • Mean residence time $\bar{t}_s > \tau$ indicates internal recirculation loops.

Heterogeneous Catalysis and Multi-Phase Systems

• Seven Consecutive Steps of Heterogeneous Solid Catalysis:

  • 1) External mass transfer of reactants from bulk fluid to the external pellet surface across the boundary layer.
  • 2) Internal diffusion of reactants through the porous network to active sites.
  • 3) Adsorption of reactants onto active catalytic surface sites.
  • 4) Chemical surface reaction transforming adsorbed reactants to adsorbed products.
  • 5) Desorption of products from active sites.
  • 6) Internal diffusion of product molecules out through the pellet pores.
  • 7) External mass transfer of products from the external surface into the bulk fluid.

• Apparent Reaction Rate ($r$):$r = k \cdot C_A^n$   combining steps 3, 4, and 5 into an apparent rate constant $k$ that may deviate from classical Arrhenius temperature dependence.

• Effectiveness Factors and Thiele Modulus:

  • Rate Relationships:$\bar{r} = \eta_s \cdot r_s = \eta_e \cdot r_e$
  • External Effectiveness Factor ($\eta_e$):$\eta_e = \frac{k_D}{k_D + k \cdot L}$     where $k_D$ is external mass transfer conductance and $L$ is characteristic pellet dimension.
  • Thiele Modulus ($\varphi_s$):$\varphi_s = L \cdot \sqrt{\frac{k \cdot C_s^{n-1}}{D_e}}$     where $D_e$ is the internal effective diffusion coefficient.
  • Internal Effectiveness Factor ($\eta_s$):
  • Spherical pellet:       $\eta_s = \frac{3}{\varphi_s} \cdot \left(\frac{1}{\tanh(\varphi_s)} - \frac{1}{\varphi_s}\right)$
  • Flat wafer pellet (half-thickness $L$):       $\eta_s = \frac{\tanh(\varphi_s)}{\varphi_s}$

• Regime Identification:

  • Chemical Kinetic Regime: $\varphi_s^2 \ll 1$; $\eta_s \approx 1$. Reactant concentration is uniform throughout pellet pores; chemical surface kinetics controls the rate.
  • Diffusion-Controlled (Physical) Regime: $\varphi_s^2 \gg 1$; $\eta_s \approx \frac{1}{\varphi_s}$. Fast reaction consumes reactant near the pellet exterior; internal diffusion controls the overall rate.
  • Mixed Regime: $0.3 < \varphi_s < 3$.

• Prater Thermicity Criterion ($\beta_P$):$\beta_P = \frac{(-\Delta_r H) \cdot D_e \cdot C_s}{\lambda_e \cdot T_s}$   evaluates the maximum intraparticle temperature gradient ($\Delta T_{\max} = \beta_P \cdot T_s$) where $\lambda_e$ is effective thermal conductivity of the pellet.

• Catalytic Reactor Selection Criteria:

  • Low reaction rate $\rightarrow$ Fixed bed with large catalyst pellets to minimize pressure drop.
  • Rapid catalyst deactivation $\rightarrow$ Moving bed or fluidized bed allowing continuous catalyst regeneration.
  • Slow internal diffusion $\rightarrow$ Eggshell (shell-impregnated) catalysts rather than mass-impregnated pellets.
  • High thermal runaway risk $\rightarrow$ Multitubular packed bed exchangers or fluidized beds providing high heat transfer.

Process Intensification and Multi-Scale Simulation

• Process Intensification (PI):
  • Pioneered by Imperial Chemical Industries (ICI) in the late 1970s via HiGee technology (centrifugal fields replacing a $30\,\text{m}$ distillation column with a $1.5\,\text{m}$ rotating packed bed).
  • Formal Definition: Substantial reduction in the physical size-to-capacity ratio of equipment, accelerated transport phenomena, and functional integration of unit operations.
  • Intensification Technologies:
  • Equipment: Microreactors, static mixers, monolithic reactors, rotating packed beds, compact plate heat-exchanger reactors.
  • Methods: Multifunctional reactors (reactive distillation, reactive extraction), alternative energy forms (microwaves, ultrasound, plasma).
  • Phenomenological Strategy Matrix: 17 potential equipment-level transport/thermodynamic limitations mapped against 17 intensification strategies to overcome mass/heat transfer bottlenecks and thermodynamic equilibrium limits.
  • Commercial Deployment Status: Over 2,000 industrial reverse-flow reactors, 300 hydrogen separation membrane units, 200 reactive distillation units, and 100 dividing-wall distillation columns.

• Hierarchical Categories of Process Models (Meyer Classification):

  • Pure Knowledge Models: Formulated exclusively from fundamental physical and chemical principles without empirical fittings.
  • Phenomenological ("Gray Box") Models: Conservation balances (mass, momentum, energy) coupled with constitutive thermodynamic equations and empirical transport correlations.
  • Hybrid Models: First-principle macroscopic balances integrated with machine learning (neural networks), commonly used for complex combustion or bioprocesses.
  • Empirical ("Black Box") Models: Derived purely from experimental data (polynomial regressions, transfer functions, neural networks, dimensional analysis).

• Multi-Scale Simulation Architecture:

  • Plant / Flowsheet Level: Modular simulators (Aspen Plus, ProSim, Pro/II) using thermodynamic databases to resolve macroscopic heat and material balances across interconnected units.
  • Apparatus / Equipment Level: 1D/2D phenomenological models resolving internal state variables.
  • Mesoscopic / Microscopic Level: Computational Fluid Dynamics (CFD via Navier-Stokes discretization) and Discrete Element Method (DEM for particle-particle contact mechanics).
  • Molecular Level: Molecular dynamics and Monte Carlo simulations predicting transport and phase equilibrium properties.
  • CAPE-OPEN Standard: International open interface specification managed by COLaN ensuring modular interoperability between different simulation platforms.

• Pilot-Scale Validation Facilities:

  • Academic and industrial pilot units bridge the gap between benchtop R&D and commercial implementation.
  • Essential for evaluating non-ideal hydrodynamics, heat and mass transfer limitations, and operational stability.
  • Standard pilot equipment: multipurpose stirred reactors, pilot rotary kilns, automated continuous distillation columns, multitubular heat exchangers, single-effect evaporators, and automated filter presses.