Reaction engineering, a core subfield of chemical engineering, is concerned with the analysis, design, and optimization of chemical reactors. Its central questions revolve around understanding how reaction kinetics, thermodynamics, and transport phenomena (mass, heat, and momentum transfer) interact within a reacting system to determine reactor performance, selectivity, and yield. The historical evolution of the field is marked by successive formalizations of theory and the emergence of durable, often competing, methodological schools for modeling and design.
The field's pre-scientific roots lie in Empirical and Heuristic Design, where reactor operation was based on craft knowledge, rules of thumb, and simple batch processes. The foundational formalization began in the early 20th century with the development of Chemical Kinetics and the Ideal Reactor Models. The work of Damköhler and others established the continuum description of reacting systems, leading to the canonical ideal models: the Batch Reactor, Continuous Stirred-Tank Reactor (CSTR), and Plug Flow Reactor (PFR). These models, based on simplifying assumptions of perfect mixing or no mixing, became the cornerstone of reactor analysis and introduced the core methodology of Mole Balances.
A major historical transition was the systematic incorporation of transport limitations. The Resistance-in-Series Model and the development of Heterogeneous Catalysis theory forced the field to move beyond homogeneous kinetics. This led to the paradigm of Macro-Kinetics, which treats the observed rate as a function of intrinsic kinetics and transport steps. The Langmuir-Hinshelwood-Hougen-Watson (LHHW) Kinetics framework provided a dominant mechanistic model for surface-catalyzed reactions. To account for internal diffusion in catalyst pellets, the Thiele Modulus and Effectiveness Factor concepts were formalized, creating a lasting school for analyzing pore diffusion effects.
The mid-20th century saw the rise of sophisticated analytical and semi-analytical schools for modeling non-ideal flow. The Residence Time Distribution (RTD) theory, pioneered by Danckwerts, provided a diagnostic tool to characterize flow patterns without knowing the exact flow field. This supported the Tanks-in-Series and Dispersion Models as competing paradigms for representing flow non-ideality, often used in conjunction with the Segregated Flow Model and Maximum Mixedness Model to predict conversion.
The advent of digital computation catalyzed a new era but did not erase prior schools; it created new, coexisting methodological families. The Finite Difference Method and Finite Element Method (FEM) became dominant computational schools for solving the complex partial differential equations (PDEs) describing multidimensional reacting flows with coupled transport. This enabled the direct numerical solution of detailed Computational Fluid Dynamics (CFD) models coupled with reaction, a paradigm that now rivals and complements simpler engineering models.
A key modern methodological divide exists between continuum-level and molecular-level approaches. The continuum Population Balance Modeling (PBM) school emerged as the canonical framework for systems where particle size, age, or other distributed properties evolve (e.g., polymerization, crystallization). Concurrently, the Microkinetic Modeling school, which builds reaction networks from elementary steps on catalytic surfaces, represents a rival, more fundamental paradigm to empirical LHHW kinetics. The current landscape is defined by the challenge of Multiscale Modeling, which seeks to formally integrate models across scales—from molecular (Density Functional Theory (DFT)) to pellet (PBM, microkinetics) to reactor (CFD, ideal models)—creating a new methodological frontier rather than a single unified school.
**