Table of Contents
Cover
Title Page
List of contributors
Preface
PART 1: Principles of catalytic reaction engineering
CHAPTER 1: Catalytic reactor types and their industrial significance
1.1 Introduction
1.2 Reactors with fixed bed of catalysts
1.3 Reactors with moving bed of catalysts
1.4 Reactors without a catalyst bed
1.5 Summary
References
CHAPTER 2: Microkinetic analysis of heterogeneous catalytic systems
2.1 Heterogeneous catalytic systems
2.2 Intrinsic kinetics of heterogeneous reactions
2.3 External (interphase) transport processes
2.4 Internal (intraparticle) transport processes
2.5 Combination of external and internal transport effects
2.6 Summary
Nomenclature
Greek letters
References
PART 2: Two‐phase catalytic reactors
CHAPTER 3: Fixed‐bed gas–solid catalytic reactors
3.1 Introduction and outline
3.2 Modeling of fixed‐bed reactors
3.3 Averaging over the catalyst particle
3.4 Dominant fluid–solid mass transfer
3.5 Dominant fluid–solid mass and heat transfer
3.6 Negligible mass and thermal dispersion
3.7 Conclusions
Nomenclature
Greek letters
References
CHAPTER 4: Fluidized‐bed catalytic reactors
4.1 Introduction
4.2 Key hydrodynamic features of gas‐fluidized beds
4.3 Key properties affecting reactor performance
4.4 Reactor modeling
4.5 Scale‐up, pilot testing, and practical issues
4.6 Concluding remarks
Nomenclature
Greek letters
References
PART 3: Three‐phase catalytic reactors
CHAPTER 5: Three‐phase fixed‐bed reactors
5.1 Introduction
5.2 Hydrodynamic aspects of three‐phase fixed‐bed reactors
5.3 Mass and heat transfer in three‐phase fixed‐bed reactors
5.4 Scale‐up and scale‐down of trickle‐bed reactors
5.5 Trickle‐bed reactor/bioreactor modeling
Nomenclature
Greek letters
Subscripts
Superscripts
Abbreviations
References
CHAPTER 6: Three‐phase slurry reactors
6.1 Introduction
6.2 Reactor design, scale‐up methodology, and reactor selection
6.3 Reactor models for design and scale‐up
6.4 Estimation of transport and hydrodynamic parameters
6.5 Advanced computational fluid dynamics (CFD)‐based models
6.6 Summary and closing remarks
Acknowledgments
Nomenclature
Greek letters
Subscripts
References
CHAPTER 7: Bioreactors
7.1 Introduction
7.2 Basic concepts, configurations, and modes of operation
7.3 Mass balances and reactor equations
7.4 Immobilized enzymes and cells
7.5 Aeration
7.6 Mixing
7.7 Heat transfer
7.8 Scale‐up
7.9 Bioreactors for animal cell cultures
7.10 Monitoring and control of bioreactors
Nomenclature
Greek letters
Subscripts
References
PART 4: Structured reactors
CHAPTER 8: Monolith reactors
8.1 Introduction
8.2 Design of wall‐coated monolith channels
8.3 Mapping and evaluation of operating regimes
8.4 Three‐phase processes
8.5 Conclusions
Nomenclature
Greek letters
Superscripts
Subscripts
References
CHAPTER 9: Microreactors for catalytic reactions
9.1 Introduction
9.2 Single‐phase catalytic microreactors
9.3 Multiphase microreactors
9.4 Conclusions and outlook
Nomenclature
Greek letters
Subscripts
References
PART 5: Essential tools of reactor modeling and design
CHAPTER 10: Experimental methods for the determination of parameters
10.1 Introduction
10.2 Consideration of kinetic objectives
10.3 Criteria for collecting kinetic data
10.4 Experimental methods
10.5 Microkinetic approach to kinetic analysis
10.6 TAP approach to kinetic analysis
10.7 Conclusions
References
CHAPTER 11: Numerical solution techniques
11.1 Techniques for the numerical solution of ordinary differential equations
11.2 Techniques for the numerical solution of partial differential equations
11.3 Computational fluid dynamics techniques
11.4 Case studies
11.5 Summary
Nomenclature
Greek letters
Subscripts/superscripts
References
PART 6: Industrial applications of multiphase reactors
CHAPTER 12: Reactor approaches for Fischer–Tropsch synthesis
12.1 Introduction
12.2 Reactors to 1950
12.3 1950–1985 period
12.4 1985 to present
12.5 The future?
References
CHAPTER 13: Hydrotreating of oil fractions
13.1 Introduction
13.2 The HDT process
13.3 Fundamentals of HDT
13.4 Process aspects of HDT
13.5 Reactor modeling and simulation
Nomenclature
Greek letters
Subscripts
Non‐SI units
References
CHAPTER 14: Catalytic reactors for fuel processing
14.1 Introduction—The basic reactions of fuel processing
14.2 Theoretical aspects, advantages, and drawbacks of fixed beds versus monoliths, microreactors, and membrane reactors
14.3 Reactor design and fabrication
14.4 Reformers
14.5 Water-gas shift reactors
14.6 Carbon monoxide fine cleanup: Preferential oxidation and selective methanation
14.7 Examples of complete fuel processors
Nomenclature
References
CHAPTER 15: Modeling of the catalytic deoxygenation of fatty acids in a packed bed reactor
15.1 Introduction
15.2 Experimental data for stearic acid deoxygenation
15.3 Assumptions
15.4 Model equations
15.5 Evaluation of the adsorption parameters
15.6 Particle diffusion study
15.7 Parameter sensitivity studies
15.8 Parameter identification studies
15.9 Studies concerning the deviation from ideal plug flow conditions
15.10 Parameter estimation results
15.11 Scale‐up considerations
15.12 Conclusions
Acknowledgments
Nomenclature
Greek letters
References
Index
End User License Agreement
List of Tables
Chapter 02
Table 2.1 Individual terms of LHHW rate equations for the surface reaction‐controlling cases of various catalytic reactions.
Chapter 03
Table 3.1 Fixed‐bed reactor mathematical model (definition of timescales
τ
in Table 3.2).
Table 3.2 Timescales of the different processes in the model for a fixed bed.
Chapter 04
Table 4.1 Some typical key characteristics of catalytic fluidized‐bed reactors compared with those of alternative types of reactor.
Table 4.2 Solid catalyzed gas‐phase reactions which have been carried out in commercial fluidized‐bed reactors.
Table 4.3 Typical operating ranges and features of catalytic fluidized‐bed reactors.
Table 4.4 Comparison of typical properties of catalytic and gas–solid fluidized‐bed reactors.
Table 4.5 Common instrumentation required in fluidized‐bed reactors.
Table 4.6 Complementary components of fluidized‐bed reactor system.
Chapter 05
Table 5.1 Model parameters.
Table 5.2 Values of kinetic parameters.
Table 5.3 Parameters used in simulations.
Table 5.4 Two‐bed reactor operating conditions.
Chapter 06
Table 6.1 Three‐phase catalytic reactors.
Table 6.2 Comparison of multiphase reactors (qualitative rating: more stars mean superior performance on the pertinent metric).
Table 6.3 Some illustrative applications of three‐phase slurry reactors.
Table 6.4 Overview of vessel designs and performance attributes of three‐phase slurry reactors.
Table 6.5 Idealized flow and axial dispersion models (steady state).
Table 6.6 Mixing cell model.
Table 6.7 Correlations for estimation of
k
l
a
in three‐phase slurry and fluidized beds.
Table 6.8 Liquid holdup, mass transfer coefficients, and effective interfacial area in gas–liquid reactors.
Table 6.9 Correlations for liquid dispersion coefficients in three‐phase slurry and three‐phase fluidized beds.
Table 6.10 Governing equations for Euler–Euler formulation.
Table 6.11 Turbulence closures.
Table 6.12 Closures for interphase momentum exchange.
Table 6.13 Closures for solid phase.
Chapter 08
Table 8.1 Monolith reactors classified according to flow, materials, and operation features.
Table 8.2 Typical dimensions and pressure drop in some monolith reactor processes [46].
Table 8.3 Friction factors and transfer coefficients for some common cross‐sectional shapes in monolith channels [39].
Table 8.4 Solution strategies for the Graetz–Lévêque problem with wall reaction.
Table 8.5 Bridging the gap between convection and diffusion regimes.
Table 8.6 Calculation of the effectiveness factor in uniform and nonuniform washcoats.
Table 8.7 Accuracy of effectiveness factor calculation methods for nonuniform geometry (circle‐in‐square shape with a first‐order reaction).
Table 8.8 Experimental ranges used in the development of empirical mass transfer correlations.
Table 8.9 Comparison between models for catalytic combustion in a monolith.
Table 8.10 Vertices in the Damköhler–Graetz plot with reaction–transport and profile development regimes.
Chapter 09
Table 9.1 Asymptotic Sherwood number for constant reactant concentration.
Table 9.2 Two‐phase laminar–laminar frictional pressure drop correlations.
Table 9.3 Comparison of mass transfer parameters in different gas–liquid contactors [81].
Chapter 11
Table 11.1 Linear system solvers available in COMSOL Multiphysics.
Table 11.2 Model equations used to simulate heat‐exchange integrated microchannel reactor operation shown in Figure 11.8.
Table 11.3 Boundary conditions related to the mathematical model given in Table 11.2.
Chapter 12
Table 12.1 Comparison of reaction engineering models for the Fischer–Tropsch synthesis in slurry bubble column reactor.
Table 12.2 Kinetic models for Table 12.1.
Table 12.3 Existing operating Fischer–Tropsch plants.
Chapter 13
Table 13.1 Properties of various Mexican crude oils (obtained at IMP analytical lab).
Table 13.2 Standard enthalpies and equilibrium constants of representative HDT reactions.
Table 13.3 Typical process conditions of various hydrotreating and hydrocracking processes.
Table 13.4 Properties of the feedstocks.
Chapter 15
Table 15.1 The values of input variables in laboratory‐scale unit [15, 20].
Table 15.2 Model equations.
Table 15.3 Physical properties and variables of the model system.
Table 15.4 Parameter values of the sensitivity study.
List of Illustrations
Chapter 01
Figure 1.1 Outline of a chemical process.
Figure 1.2 Schematic presentation of a packed‐bed reactor.
Figure 1.3 Furnace configurations for multitubular packed‐bed reformers.
Figure 1.4 Side‐fired tubular reformer design by Haldor‐Topsøe.
Figure 1.5 Heat transfer strategies in multitubular packed‐bed reactors. (a) Cross‐flow, (b) parallel flow, and (c) boiling‐water cooling.
Figure 1.6 Various configurations of vessel‐type packed‐bed reactors. (a) Single‐bed adiabatic packed‐bed reactor, (b) adiabatic reactor with interstage gas injection, and (c) multiple adiabatic beds with interstage heat exchange.
Figure 1.7 Packed‐bed reactor with multiple adiabatic beds for ammonia synthesis.
Figure 1.8 Packed‐bed reactor configuration for autothermal reforming of methane to synthesis gas.
Figure 1.9 Schematic presentation of a monolith reactor.
Figure 1.10 Comparison of pressure drop in various configurations of monoliths and packing structures.
Figure 1.11 Radial flow reactor concept.
Figure 1.12 Radial flow ammonia synthesis converter by Haldor‐Topsøe.
Figure 1.13 Disk reactor concept.
Figure 1.14 Schematic presentation of a microchannel reactor. (a) Machined plates with microchannels, (b) microchannel reactor block obtained after bonding the plates, and (c) characteristic section of the multichannel reactor.
Figure 1.15 Riser cracking process by UOP. (a) Reactor, (b) stripper, (c) riser, (d) slide valve, (e) air grid, and (f) regenerator.
Figure 1.16 High‐temperature Fischer–Tropsch synthesis reactors. (a) Sasol Synthol circulating fluidized‐bed reactor.(b) Sasol Advanced Synthol turbulent fluidized‐bed reactor.
Figure 1.17 Slurry bubble column reactor for low‐temperature Fischer–Tropsch synthesis.
Figure 1.18 Stirred‐tank reactor with typical dimensions.
Figure 1.19 Heat transfer strategies in stirred‐tank reactors. (a) Jacket, (b) internal coils, (c) internal tubes, (d) external heat exchanger, (e) external reflux condenser, and (f) fired heater.
Chapter 02
Figure 2.1 Potential energy curves representing the action of a solid catalyst.
Figure 2.2 Effect of temperature on amount of gas adsorbed for simultaneous physical adsorption and activated chemisorption.
Figure 2.3 Reactant concentration profiles in different global rate regimes: I, external mass transfer limitation; II, pore diffusion limitation; III, both external and internal mass transfer limitations; IV, no mass transfer limitations on the intrinsic rate.
Figure 2.4 Total pressure dependence of initial rates for the reaction
A
→
B
+
C
.
Figure 2.5 Types of laboratory catalytic reactors mostly used: (a) plug‐flow packed‐bed reactor and (b) spinning basket reactor mimicking CSTR performance.
Figure 2.6 The effect of external mass transport limitations on conversion rates.
Figure 2.7 Correlation of experimental data on external mass transport between bulk fluid and particle surface in packed beds by various researchers.
Figure 2.8 Steady‐state multiplicity caused by external temperature gradients in simple nonporous surface‐catalyzed exothermic reactions when external mass transfer resistance is either considerable or negligible.
Figure 2.9 Representation of a catalyst particle having micro‐ and macropore structure.
Figure 2.10 Internal effectiveness factor as a function of Thiele modulus for porous particles of various shapes.
Figure 2.11 Reactant concentration and temperature profiles inside a porous spherical catalyst particle with exothermic reaction occurring at steady state.
Figure 2.12 Nonisothermal internal effectiveness factors for first‐order reactions in porous spherical particles.
Chapter 03
Figure 3.1 Fixed‐bed reactors. Schematic representation of several operating configurations and arrangements: (a) Cylindrical with axial flow (a.1‐single, a.2‐multistage, a.3‐multitubular); (b) cylindrical with radial flow; (c) spherical reactor (c.1‐axial and c.2‐radial flow).
Figure 3.2 Transport–reaction mechanisms in a fixed‐bed reactor at different scales.
Figure 3.3 Effectiveness factor curve (
η
vs.
φ
) for a spherical catalyst particle with Dirichlet boundary condition and first‐order exothermic reaction. Numerical results from Ref. [99]. A curve for
is depicted (with
), while the others were calculated for
and several values of
β
.
Figure 3.4 Asymptotic forms of the effectiveness factor (
) as a function of the Thiele modulus (
φ
) for a slab catalyst with a first‐order exothermic reaction (
and
) for
. Numerical values in the intermediate region are given in Ref. [80].
Figure 3.5 Relative error (
δ
%) associated with the calculation of conversion from a pseudo‐homogeneous model for a first‐order isothermal reaction as a function of
, with
or 0.5.
Figure 3.6 Relative error (
δ
%) associated with the calculation of conversion from a pseudo‐homogeneous model for an
m
th‐order isothermal reaction as a function of
, with
or 0.5.
Figure 3.7
regime diagram for a second‐order reaction occurring in an isothermal fixed bed. Iso‐
δ
lines given by Equation 3.85.
Figure 3.8 Cross‐sectional average concentration (a) and temperature (b) profiles according to a pseudohomogeneous model (full lines: two‐dimensional; dashed lines: one‐dimensional). After Ref. [133] with
°C and
°C.
Figure 3.9 Deviation in the reaction rate between pseudo‐homogeneous and heterogeneous models in the partial oxidation of methanol to formaldehyde, as calculated in Ref. [135].
Figure 3.10 Comparison of the axial temperature (a) and methanol concentration (b) according to several models (cross‐sectional averages for 2D models) for a fixed‐bed reactor to produce formaldehyde from methanol [135].
Chapter 04
Figure 4.1 Schematic of principal features of a modern fluid catalytic cracking (FCC) reactor system.
Figure 4.2 Schematic of main features of fluidized‐bed reactor for production of acrylonitrile showing immersed serpentine heat transfer tubing, separate distributors for air and ammonia + propane or propene, and one representative cluster of internal cyclones in series.
Figure 4.3 Diagrammatic representation of flow regimes of gas fluidization with increasing superficial gas velocity, with transition velocities indicated on the arrows between adjacent flow regimes. Dashed rectangular box encloses the five flow regimes of primary interest with respect to catalytic fluidized‐bed reactors.
Figure 4.4 Pictorial representation of particles and void regions in flow regimes of principal interest for catalytic fluidized‐bed reactors.
Figure 4.5 Simple single‐phase model predictions for first‐order irreversible catalytic ozone decomposition reaction in comparison with experimental fluidized‐bed reactor data of Sun and Grace [44].
Figure 4.6 Schematic of two‐phase model representation of bubbling or slugging fluidized‐bed reactor.
Chapter 05
Figure 5.1 Fixed‐bed reactors for gas–liquid–solid catalyzed systems. (a) Two‐phase downflow fixed‐bed reactor (trickle‐bed reactor). (b) Fixed‐bed reactor with countercurrent flow. (c) Two‐phase upflow fixed‐bed reactor (packed‐bed bubble flow reactor).
Figure 5.2 Hydrodynamic regimes for two‐phase downflow in fixed beds—nonfoaming (a) and foaming (b) liquids.)
Figure 5.3 Hydrodynamic regimes for two‐phase upflow in fixed beds: transition zones (/////) from Ref. [30] and transition lines from Ref. [28].)
Figure 5.4 Representation of connectivity between subregions in the representative elementary volume:
V
g
,
V
ℓ
, and
V
s
for the bulk phases and
A
gℓ
and
A
ℓs
for the interfaces (a), interaction of interface
A
gℓ
with
V
g
and
V
ℓ
, bulk phases in creating and destroying gas–liquid interfacial area (b).)
Figure 5.5 Mass transfer resistances in three‐phase fixed‐bed reactors.
Figure 5.6 Liquid and gas loads in typical trickle‐bed reactors for gas oil hydroprocessing.)
Figure 5.7 Sketch of trickle‐bed reactor (a) and local‐scale representation of the two‐phase flow (b). Sketch showing the resistances involved in the G/L reaction on the catalyst particle (c, d).)
Figure 5.8 Two‐phase pressure drop ratio (a) and bed volume average specific deposit (b) versus time at different values of fine particle concentration (
d
f
= 5 µm,
T
in
= 598 K,
P
in
= 10 MPa).
Figure 5.9 Effect of fine particle deposition process on dibenzothiophene conversion at
t
= 485 min (
T
in
= 598 K,
P
in
= 10 MPa): (a)
d
f
= 5 µm, (b)
d
f
= 1 µm.
Figure 5.10 Biofilm thickness versus time at different values of inlet cell concentration (only cell detachment was considered): (a)
; (b)
(
,
,
,
,
α
= 0.1,
D
f
= 2.0).
Figure 5.11 Number of aggregates versus axial distance (
mg/l,
mg/l,
mg/l,
g/l,
α
= 0.1,
D
f
= 2.0,
t
= 2000 min, surface 1—
, surface 2—
).)
Figure 5.12 Two‐phase pressure drop ratio. (a) Cell + aggregate volume average specific deposit and biofilm thickness at
z/H
= 0.11 (b) versus time in the presence or absence of biomass detachment (
,
,
,
,
,
α
= 0.5,
D
f
= 1.5).
Figure 5.13 Schematic diagram of two‐bed reactor system for DME production from glycerol via an integrated process involving aqueous‐phase glycerol reforming coupled with DME synthesis process.
Figure 5.14 Steady‐state axial profiles of the volumetric gas flow rate (a) and temperature (b) for aqueous‐phase glycerol reforming and DME synthesis without
in situ
H
2
O removal.
Figure 5.15 Steady‐state axial profiles of the volumetric gas flow rate for aqueous‐phase glycerol reforming and DME synthesis with
in situ
H
2
O removal. (a)
. (b)
.
Figure 5.16 Steady‐state axial profiles of the temperature for aqueous‐phase glycerol reforming and DME synthesis with/without
in situ
H
2
O removal. (a)
. (b)
.
Chapter 06
Figure 6.1 Schematic diagrams of industrial three‐phase slurry reactors. (a) Slurry bubble column. (b) Three‐phase fluidized bed. (c) Three‐phase agitated vessel.
Figure 6.2 Schematic of the overall methodology for relating performance of the reactor with the hardware and the operating protocol.
Figure 6.3 Possible ideal contacting patterns in three‐phase slurry reactors. (a) Countercurrent (gas and liquid in plug flow). (b) Co‐current (gas and liquid in plug flow). (c) Mixed (gas and liquid in mixed flow).
Figure 6.4 Reactant profiles in neighborhood of gas (bubble) and catalyst particle.
Figure 6.5 Hydrodynamic simulation of methanol synthesis in gas–liquid slurry bubble column reactors. (a) Reactor operating conditions and simulation grid. (b) Simulated instantaneous distribution of gas and solid holdup. (c) Simulated instantaneous distribution of temperature.
Figure 6.6 Predicted instantaneous distribution of various species in (a) gas and (b) slurry phases at catalyst concentration of 20%.
Chapter 07
Figure 7.1 An estimate on the current use in industry of bioreactors with different basic characteristics regarding aeration and stirring.
Figure 7.2 Typical configuration of commonly used bioreactors.
Figure 7.3 Modes of bioreactor operation: (a) batch, (b) continuous, and (c) fed batch.
Figure 7.4 Flow through a tubular reactor of total length
L
, where a differential element of volume,
dV
=
Adz
, is highlighted.
Figure 7.5 Different stages of cell growth in batch processes: (1) lag phase, (2) exponential growth phase, (3) stationary phase, and (4) death phase.
Figure 7.6 CSTR with recycle.
Figure 7.7 Variation of substrate concentration,
C
A
, from bulk liquid medium to center of support particle.
Chapter 08
Figure 8.1 Representation of the several scales in a catalytic monolithic reactor. (a) Monolith honeycomb [1] . (b) Single washcoated monolith channel. (c) Catalytic washcoat.
Figure 8.2 Conversion calculated according to Graetz series with different number of terms. Equation 8.26 is plotted with 1, 5, and 12 terms as
X
Graetz
for intermediate mass transfer control (
). The solution with only one term corresponds to the fully developed concentration profile. A uniformly valid solution for conversion
X
(plotted as a solid line) is given by Equation 8.37 and is in good agreement with the numerical solution.
Figure 8.3 Normalized dependence of the first eigenvalue on the Damköhler number:
versus
curve. The analytical result given by Equation 8.34b was derived in Lopes et al. [40]. Numerical results for laminar flows are also shown (dashed lines).
Figure 8.4 Fully developed concentration profile given by Equation 8.34. Numerical (full lines) and analytical (dashed lines) refer to laminar flow in a circular channel. The approximate solution uses eigenvalues and coefficients calculated from Equations 8.34b and 8.34c for several values of Damköhler number.
Figure 8.5 Reactant concentration (mixing cup) as a function of the rescaled Damköhler number (
Da
*
) given in Equation 8.36b for laminar flow between parallel plates. The dashed lines are plotted according to Equation 8.35b for several values of the Graetz parameter (
αPe
m
/
z
) in the ranges of
Da
*
leading to less than 5% relative error.
Figure 8.6 Conversion as a function of the dimensionless axial distance for laminar flow in a circular channel with instantaneous wall reaction. Three regions can be identified in different ranges of the Graetz parameter: Graetz and Lévêque descriptions, separated by a transition regime.
Figure 8.7 Convection–diffusion limiting cases for channel mass transfer with wall reaction in a Damköhler–Graetz parametric map (
). Four regions can be depicted: fully developed profile (for
αPe
m
/
z
below the prediction shown by the full lines for given values of
ε
fd
), developing profile (for Graetz numbers above the ones estimated by dashed lines for given
e
dev
), transition region (separating the two previous limits at high
Da
), and overlapping region (where similar predictions from both models are obtained at low
Da
). All boundaries were obtained analytically and the respective expressions can be found in Ref. [43].
Figure 8.8 Calculation approaches for the effectiveness factor in thin catalytic coatings. The distinction between the conventional method and the one proposed in Lopes et al. [94] is illustrated. The Thiele modulus for a first‐order reaction is given by
. The ratio between the catalyst and channel transverse length scales is given by
and
σ
is the catalyst shape factor.
Figure 8.9 Damköhler–Graetz plot for linear kinetics with mapping of reaction–transport and convective–diffusive regimes in a washcoated monolith. Two values of the diffusion ratio are examined, and laminar flow inside a circular channel with
is considered. Regime boundaries are plotted for given values of the criteria (internal control,
; external control,
; no internal limitation,
; no external limitation,
; developed profile,
; developing profile,
). Vertices (
V
) delimiting the intermediate regime are also shown.
Figure 8.10
operating map for a circular channel with laminar flow coated with an annular catalyst layer (
) and a first‐order reaction. External and internal mass transfer control observed for
and
, respectively. Negligible resistances for
and
. Full lines associated with conditions in channel are plotted for values of the Graetz parameter. The area for overall kinetic control (Ov. KC) is depicted as well as the ones for mass transfer control: overall (Ov. MTC) and internal (Intraph. MTC).
Figure 8.11 Damköhler–Graetz diagram for nonlinear kinetics. A second‐order reaction occurs in the washcoat, but the remaining conditions are identical to previous representations.
Figure 8.12 Damköhler–diffusion ratio diagram for nonlinear kinetics (second‐order reaction). Same conditions from previous representations apply.
Figure 8.13 Design in the presence of a constraint on conversion. Two distinct problems are depicted: (i) increment in pressure drop due to the existence of a maximum reaction rate above which operation is not feasible and (ii) reduction in energetic requirements by increasing the allowed pressure drop.
Figure 8.14 Taylor flow in washcoated monoliths. Gas–liquid–solid mass transfer routes are represented schematically. The relevant geometrical dimensions can be also depicted.
Chapter 09
Figure 9.1 (a) Influence of channel diameter variations on the reactor conversion; ideal reactor (thick line), reactor with
= 0.1 (thin line). (b) Influence of the parameter
γ
on the conversion in microreactor with
= 0.01;
γ
= 15, 30, and 60.
Figure 9.2 Typical ranges of force ratios in multiphase flow for packed beds, micropacked beds, and fine porous media.
is the gravity/surface tension ratio and
. The term
u
is the superficial velocity,
ρ
is the density,
l
is the characteristic hydrodynamic length, that is, wetted area over wetted perimeter,
g
is gravity,
μ
is the viscosity, and
σ
is the surface tension. The subscripts
L
and
G
refer to liquid and gas, respectively.
Figure 9.3 Flow transition for cyclohexene/H
2
cocurrent flow in a micro trickle bed reactor (■) versus superficial mass velocity [23]. Charpentier flow map (♦) is presented for comparison [24].
G
is the superficial mass velocity.
2
2
2
Figure 9.8 Mechanisms of the formation of gas bubbles in the microchannels: (a, b) squeezing and (c) dripping.
Figure 9.9 Schematic of Taylor flow showing the definitions of the unit cell, gas bubble length
L
b
and the liquid slug length
L
s
. The lengths of the nose
L
nose
and tail
L
tail
sections of the gas bubble are also indicated. Mass transfer contributions in the Taylor flow are indicated with the numbers: (1) bubble to wall through film, (2) bubble to slug, and (3) slug to wall.
Chapter 10
Figure 10.1 Conceptual diagram of a molecular beam scattering (MBS) experiment comprised of pulsed beam source, sample target, mass spectrometer detector, and a differentially pumped ultrahigh‐vacuum system.
Figure 10.2 Plot showing the scattered CO intensity as a function of time from an initially clean Pd(111) surface. The initial sticking coefficient can be determined from the mass spectrometer signal
I
(
t
) using the relationship
S
(
t
) = 1 −
I
(
t
)/
I
(∞).
Figure 10.3 S
0
for O
2
on Pd (111) as a function of substrate temperature.
Figure 10.4 Adaptation of the “pressure/materials” correlation proposed by Bonzel [37] indicating pressure regimes and types of materials found in various types of catalyst characterization experiments.
Figure 10.5 Key components of a TAP experiment.
Figure 10.6 (a) Schematic of TAP single particle microreactor configuration. The 400 µm diameter Pt particle is packed within a bed of inert quartz particles with diameters between 210 and 250 µm. (b) Image comparing a 400 µm Pt particle to a pencil point. (c) SEM image showing the complex surface structure of a polycrystalline Pt particle. (d) Higher magnification (15 000×) of the particle shown in (c), which shows the surface is nonporous [86].
Figure 10.7 Illustration of a TAP pump‐probe experiment in which O
2
/Ar and CO
2
/Ar are pulsed in an alternating sequence and the CO
2
transient response is measured during each pulse.
Figure 10.8 Comparison of CO
2
produced during TAP vacuum pump‐probe experiments and atmospheric flow experiments for CO oxidation over single Pt particle with the same composition of reactants. (a) A typical set of pump‐probe CO
2
responses (
m
/
e
= 44) for reaction at 140, 170, and 350°C. There is a shift in the amount of CO
2
produced during both CO and oxygen pulses as temperature increases. (b) CO
2
production observed from atmospheric flow experiment. The CO
2
produced while increasing reactor temperature is less than the CO
2
produced during reactor temperature decrease as shown by the counterclockwise hysteresis loop. (c) CO
2
production observed from vacuum pump‐probe experiment. The black line represents the total CO
2
yield. The circle and diamond points represent the CO
2
yield on the oxygen pulse and CO pulse, respectively.
Figure 10.9 (a) Oxygen uptake over reactor‐equilibrated VPO at 800 torr O
2
. (b) O
2
desorption spectrum under vacuum from
18
O
2
‐treated (VO)
2
P
2
O
7
.
Figure 10.10 (a)
n
‐Butane pulse response curves over oxygen‐treated VPO at various temperatures. Each curve is obtained for the same initial VPO oxidation state. (b) Arrhenius plot obtained from the temperature dependence of the
n
‐butane conversion giving an activation energy of 12 kcal/mol.
Figure 10.11 (a) Arrhenius plots for
n
‐butane oxidation over a single VPO sample reduced by a long series of
n
‐butane pulses. As the surface is reduced, the activation energy increases. (b) Apparent equilibrium constant for various products as a function of VPO oxidation state.)
Figure 10.12 Pump‐probe data showing CO
2
production as a function of temperature and pump‐probe interval. At 150°C the CO
2
production is essentially independent of the pump‐probe interval out 9 s separation. At 350°C the CO
2
production drops as the pump‐probe interval increases, indicating a drop in the active oxygen concentration with time.
Figure 10.13 Normalized CO
2
production on the CO pulse calculated from the zeroth moment of the pulse response curve. CO
2
production is constant at 150°C and drops to ca. 0.5 times its value in 9 s at 350°C.
Chapter 11
Figure 11.1 Structure of a finite element grid.
Figure 11.2 One‐dimensional finite elements. (a) A linear element, (b) a quadratic element, and (c) linear and (d) quadratic variation of field variable
ϕ
over an element. Black circles represent the nodes, while the white, numbered circles represent the particular element.
Figure 11.3 Schematic presentation of a monolith reactor with porous catalyst‐coated channels
Figure 11.4 Steady‐state reactor temperature profiles obtained at different values of molar CH
4
/O
2
and H
2
O/CH
4
ratios at the inlet.
Figure 11.5 Transient reactor temperature profile obtained at CH
4
/O
2
= 2.24 and H
2
O/CH
4
= 1.17.
Figure 11.6 Transient reactor temperature profile obtained upon changing CH
4
/O
2
from 2.24 to 1.89 (H
2
O/CH
4
= 1.17).
Figure 11.7 Transient reactor temperature profile obtained upon changing H
2
O/CH
4
from 1.17 to 1.56 (CH
4
/O
2
= 1.89).
Figure 11.8 Frontal view of the heat‐exchange integrated microchannel reactor concept and the two‐dimensional unit cell (shown at the lower inset).
Figure 11.9 The first one‐tenth of the computational grid used in the 2D simulation of the microchannel reactor for coupling catalytic methane combustion and
i
‐octane steam reforming.
Figure 11.10 Temperature profiles obtained from the solution of 3D model (side length of the microchannel = 5.6 × 10
−4
m, thickness of the wall separating the channels = 2 × 10
−4
m, reactor material: AISI‐steel).
Figure 11.11 Profiles of average bulk temperatures along the steam reforming channel obtained at different wall thickness values by the solution of the 2D model.
Figure 11.12 Hydrogen yields and
i
‐octane conversions (in parentheses) obtained for different (a) wall thickness values, (b) wall materials, and (c) channel side lengths by the solution of the 2D model.
Figure 11.13 Two‐dimensional unit cell with (a) straight‐through (ST) and (b) microbaffled (BF) channel wall configurations.
Figure 11.14 The first one‐third of the computational grid used in the 2D simulation of the microchannel reactors with ST and BF wall configurations for coupling catalytic methane combustion and
i
‐octane steam reforming.
Chapter 12
Figure 12.1 Scale of catalytic materials and reactors [2].
Figure 12.2 Schematic of some possible forms of catalyst pellets, the smallest chemical reactor [9].
Figure 12.3 Schematic of German reactors used in the pre‐1955 period according to Tramm [10].
Figure 12.4 Schematic of the slurry reactor with multilayer beds [11].
Figure 12.5 Possible reactors for Fischer–Tropsch synthesis. (a) Slurry bubble column reactor, (b) multitubular trickle bed reactor, (c) circulating fluidized‐bed reactor, and (d) fixed fluidized‐bed reactor [15].
Figure 12.6 Relative activity profile of an precipitated iron catalyst after different periods of time on stream. Times are in relative units (
: 1 tu;
, 50 tu;
, 270 tu; x, 1000 tu (tu = time units)).
Figure 12.7 Flow regime maps for trickle‐bed reactors [17].
Figure 12.8 Examples of monolith reactors illustrating cell size and shape possibilities [41].
Figure 12.9 The Velocys commercial‐scale Fischer–Tropsch reactor.
Figure 12.10 Comparison of catalyst productivity for Shell, Sasol, and Velocys Oxford Catalysts Group (OCG) reactors (courtesy of Velocys).
Figure 12.11 Productivity of commercial Velocys microchannel, fixed‐bed, and slurry reactors compared to that of the reactor scale (courtesy of Velocys).
Figure 12.12 Profile of reactor temperature at the top, middle, and bottom of a small channel reactor during Fischer‐Tropsch synthesis [48].
Figure 12.13 CO conversion, hydrocarbon productivity, and methane selectivity for Fischer–Tropsch synthesis at 214°C with a Pt–Co–alumina catalyst [48].
Figure 12.14 Model for structured bed of catalyst in a straight flow channel reactor [50].
Figure 12.15 Steady‐state performance of a structured catalyst bed (25 bar) [50].
Figure 12.16 Illustration of possible types of slurry bubble column reactors. (a) Simple bubble column, (b) cascade bubble column with sieve trays, (c) packed bubble column, (d) multishaft bubble column, and (e) bubble column with static mixers [61].
Figure 12.17 Schematic of an Exxon slurry reactor (left) [68]. Schematic of Sasol slurry bubble column reactor (right) [69].
Figure 12.18 Illustration of some concepts for the sparger for the bubble column reactor. Left, from [61]: (a) dip tube, (b) perforated plate, (c) perforated ring, and (d) porous plate. Right: closeup view of a sparger from an Exxon patent [68].
Figure 12.19 Typical particle trajectories within three different flow regions around a rising bubble in a slurry bubble column reactor [70].
Figure 12.20 Model of the flow scheme in the slurry bubble column reactor [73].
Figure 12.21 Hydrodynamic model of a slurry bubble column reactor in the heterogeneous flow regime [82].
Figure 12.22 Experimental data for gas holdup in a 0.1 m diameter bubble column operating with air–water system spanning both the homogeneous and heterogeneous flow regimes [113].
Figure 12.23 Examples of structured packing: (a) open cross flow structure (OCFS), (b) closed cross flow structure (CCFS), (c) knitted wire, and (d) aluminum foam [117].
Figure 12.24 Schematic representation of structured catalytic packing made by BUCT: (a) cross section, (b) side view, and (c) actual photography [118].
Figure 12.25 Membrane reactors for FT synthesis from the literature: (a) distributed feeding of reactants
A
and
B
, (b)
in situ
water removal by selective membrane (F, feed; S, sweep), (c1) plug‐through contactor membrane (PCM) with wide transport pores, (c2) forced‐through flow membrane contactor, product and heat removal by circulated liquid product, (d) zeolite encapsulated FT catalyst,
P
′, modified product [123].
Figure 12.26 Density versus pressure of the mixture of 55% hexane and 45% pentane at 220°C [131].
Figure 12.27 Schematic of a laboratory‐scale supercritical fixed‐bed reactor [131]. T/C: thermocouple.
Figure 12.28 CO conversion versus time on stream on a 25% Co/γ‐Al
2
O
3
slurry‐phase impregnation catalyst in a fixed‐bed reactor with varying partial pressure of supercritical fluid (SCF).
Figure 12.29 Olefin selectivity as an function of carbon number for a variety of reactors.
Figure 12.30 Cost of plant versus the plant capacity [135].
Figure 12.31 Schematic of a FT reactor suitable for a mobile biorefinery [136].
Chapter 13
Figure 13.1 Flow diagram of a typical hydrotreating unit.
Figure 13.2 Hydrotreating reactor technologies.
Figure 13.3 Process scheme of a typical petroleum refinery.
Figure 13.4 World consumption of oil refined products from 1990 to 2010.
Figure 13.5 Crude oil production in Mexico.
Figure 13.6 Reaction pathways of dibenzothiophene HDS.
Figure 13.7 Reaction scheme of quinoline HDN.
Figure 13.8 Simplified representation of the HDM of metalloporphyrins.
Figure 13.9 Examples of saturation reactions.
Figure 13.10 Reaction mechanism of paraffin hydrocracking.
Figure 13.11 First‐order rate coefficient values for the HDS of various sulfur compounds.
Figure 13.12 Multibed hydrotreating reactor with quenching.
Figure 13.13 Hydrogen circuit in a hydroprocessing unit.
Figure 13.14 Representation of the trickle‐flow regime in a hydroprocessing reactor.
Figure 13.15 Concentration profiles in a TBR reactor.
Figure 13.16 General representation of liquid quenching‐based processes. (a) Multiple feeding. (b) Product recycling.
Figure 13.17 Hydroprocessing reactor internals. (—) Axial delta‐
T
. (‐‐‐) Radial delta‐
T
.
Figure 13.18 Distributor tray design parameters. (a) Effect of tray spacing on liquid distribution, (b) comparison of tray spacing, and (c) discharge pattern of several tray designs.
Figure 13.19 Onstream catalyst replacement reactor.
Figure 13.20 H‐Oil ebullated‐bed reactor.
Figure 13.21 Slurry‐phase hydroprocessing reactor.
Figure 13.22 Sulfiding conditions of an HDS catalyst.
Figure 13.23 IMP heavy oil upgrading process.
Figure 13.24 Details of the experimental program.
Figure 13.25 Sulfur, nitrogen, and metals molar concentration profiles. Simulated: (—) liquid phase, (‐‐‐) solid phase; (■) experimental.
Figure 13.26 Hydrogen and hydrogen sulfide molar concentration profiles. Simulated: (—) liquid phase, (‐‐‐) solid phase.
Figure 13.27 Evolution of the gas‐phase composition. (—) Simulated. (■) Experimental.
Figure 13.28 H
2
/oil ratio and superficial gas velocity profiles. Simulated: (—) 380°C, (‐‐‐) 400°C; (■) experimental.
Figure 13.29 Process performance during time‐on‐stream: (symbols) experimental, (—) simulated.
Figure 13.30 HDM performance and MOC: (symbols) experimental, (—) simulated. The experimental values (■) of MOC at the end of the run were estimated based on the metals balance between feed and products.
Figure 13.31 Simulated axial MOC profiles. The MOC units are referred to the total amount of fresh catalyst in the reactors.
Figure 13.32 Catalyst wetting efficiency as function of superficial liquid mass velocity: (■) bench scale, (
) semi‐industrial scale, (—) simulated.
Figure 13.33 Comparison of semi‐industrial and bench‐scale performances.
Figure 13.34 Simulation of the proposed industrial reactor configuration. Evolution of reactor temperature, H
2
/oil ratio, and concentration of the chemical lumps. (‐‐‐) Average temperature and H
2
/oil ratio.
Figure 13.35 Average H
2
/oil ratio as function of the number of hydrogen quenches.
Figure 13.36 Evolution of the axial temperature profiles.
Chapter 14
Figure 14.1 Ethylene formation versus time on stream; (a) autothermal reforming of synthetic diesel (
O
/
C
= 1.12,
S
/
C
= 1.25 [36]).(b) autothermal reforming of commercial diesel (
S
/
C
=3,
O
/
C
=1.6 [35]).
Figure 14.2 Different feed injection systems; (a) ATR‐5, (b) ATR‐8 [37].
Figure 14.3 Different feed injection systems; (a) first generation [42].(b) improved injection system without recirculation zone [43].
Figure 14.4 Methane conversion and reactor exit temperature versus
GHSV
for a fixed bed compared to a metallic monolith (
S
/
C
= 3) [63].
Figure 14.5 (a) Jellyroll monolith prepared from platinum/ceria catalyst, (b) comparison of propane conversion, hydrogen yield and H
2
/CO ratio for fixed bed, jelly roll, and monolithic reactors as determined for different
GHSV
[64].
Figure 14.6 Hot spot formation and methane and oxygen concentration profiles in two ceramic monoliths of different cell density for partial oxidation of methane at
O
/
C
= 0.54 and at a feed inlet temperature of 362°C; (a) temperature profile for 400 cpsi, (b) temperature profile for 115 cpsi, (c) concentration profiles for 400 cpsi, (d) concentration profiles for 115 cpsi [65].
Figure 14.7 Ethylene formation versus temperature at different S/C ratio for autothermal reforming of low sulfur diesel (
O
/
C
= 0.45) [34].
Figure 14.8 Effect of feed inlet temperature and
O
/
C
ratio (here shown as
O
2
/
C
ratio) on reforming efficiency (
O
/
C
= 2.5,
GHSV
= 9 900–10 800 h
−1
) [44].
Figure 14.9 (a) System start‐up concept, (b) reformer start‐up time [43].
Figure 14.10 Species concentration versus
O
/
C
ratio (here expressed as
O
2
/
C
) for an autothermal diesel reformer at
S
/
C
= 2.3. Symbols represent experimental data, the lines represent simulation results [76].
Figure 14.11 Simulation of a cocurrently operated plate heat exchanger for diesel steam reforming coupled to a catalytic diesel burner. Left: flow arrangement, second left: fuel mass fraction in the reformer, third left: fuel mass fraction in the burner, right: temperature in the plate.
Figure 14.12 Scheme of the membrane reactor for methanol reforming [109].
Figure 14.13 Conversion in a membrane reactor applying different membrane types. For comparison, thermodynamic equilibrium conversion and conversion as achieved in a fixed catalyst bed are added [113].
Figure 14.14 Scheme of the membrane reactor for methanol oxidative steam reforming [115].
Figure 14.15 Scheme of the fluidized bed membrane reactor [117].
Figure 14.16 Scheme of the electrochemical membrane reactor for partial oxidation of methane [119].
Figure 14.17 Glycerol conversion at different
WHSV
in a membrane reactor (MR) and in a traditional reactor (TR) [125].
Figure 14.18 Schematic of the staged‐separation membrane reactor [126].
Figure 14.19 (a) Integrated water–gas shift reactor/heat exchanger designed for 5 kW fuel cell system, (b) internal temperature profiles determined during operation of this reactor; the different profiles correspond to different load (100% load is the 5 kW equivalent). The values shown for 100% load were determined at the fresh reactor (0 h) after 24 h under operation (24 h) and after an accidental exposure to reformate of too low temperature for the water–gas shift reaction (after LT) [141].
Figure 14.20 Thermal conductive catalyst plates integrated into a water–gas shift reactor. The dotted line indicates the flow paths of reformate and air [143].
Figure 14.21 (a) Axial, (b) centrifugal, and (c) centripetal flow arrangements of the water–gas shift reactors [144].
Figure 14.22 CO conversion as a function of reaction temperature for a membrane reactor (MR) and the traditional process at different
GHSV
for 350°C feed temperature and a reaction pressure of 15 bar [145].
Figure 14.23 Volume reduction by application of membrane technology for water–gas shift for different feed pressure and reaction temperature [145].
Figure 14.24 Cocurrently liquid‐cooled folded‐plate reactor for preferential oxidation [166].
Figure 14.25 (a) Cocurrently operated microreactor for preferential oxidation cooled by water evaporation, (b) CO concentration in the off‐gas of the reactor at different total reformate flow rates, CO inlet concentrations and
O
/
CO
values [167].
Figure 14.26 Temperature profile of a cocurrently cooled preferential oxidation reactor with 5 kW
el,net
power equivalent.
Figure 14.27 Flow scheme of a natural gas fuel processor/fuel cell system [168].
Figure 14.28 Gas temperature, carbon monoxide, and hydrogen content in the gas phase at different components of the natural gas fuel processor/fuel cell system [168].
Figure 14.29 Start‐up energy demand of the different components of the natural gas fuel processor/fuel cell system [168].
Figure 14.30 Sankey diagram of the 2 kW
el
CHP PEM fuel cell/methane fuel processor system [171].
Figure 14.31 Integrated methanol fuel processor with 100 W power equivalent [16, 86].
Figure 14.32 (a) Fuel processor of the VeGA system [16].(b) 250 W
el
fuel cell/fuel processor system VeGA developed by a cooperation of TRUMA and IMM [106].
Chapter 15
Figure 15.1 Simplified reaction scheme of stearic acid deoxygenation.
Figure 15.2 The experimental, the first‐order power law, and the Langmuir batch reaction curves for stearic acid [20]. Parameters:
p
2
= 0.011 l/mol,
p
3
= 25 l/mol.
Figure 15.3 Dimensionless initial concentration of stearic acid and heptadecane inside spherical catalyst particles located near the reactor inlet.
Figure 15.4 Effect of the rate constant of the main reaction on stearic acid concentration.
Figure 15.5 Effect of the rate constant of coking on stearic acid concentration.
Figure 15.6 Effect of the final activity on stearic acid concentration.
Figure 15.7 Stearic acid profiles inside the reactor.
Figure 15.8 Correlation between the main reaction rate and the final activity (
Pe
= 3).
Figure 15.9 Correlation between the deactivation rate and the final activity (
Pe
= 3).
Figure 15.10 Correlation between the main and coking reaction rate parameters (
Pe
= 3).
Figure 15.11 Identification between the rate parameter of the main reaction and the final activity.
Figure 15.12 Identification between reaction rate parameters of the main reaction and coking.
Figure 15.13 The effect of Peclet number on stearic acid consumption (catalyst: (Pd/C (Aldrich))).
Figure 15.14 The dependence of parameter estimation results on the
Pe
number.
Figure 15.15 Estimated and experimentally observed concentrations of stearic acid (catalyst: (Pd/C (Aldrich))).
Figure 15.16 Estimated and experimentally observed concentrations of stearic acid (catalyst: (Pd/C (Sibunit))).
Figure 15.17 Simulated stearic acid concentration change with time in the pilot reactor (
T
= 633°C,
p
= 20 bar, residence time, and active metal concentration same as in the laboratory scale).
Figure 15.18 Simulated stearic acid concentration change with time in the pilot reactor (
T
= 633°C,
p
= 10 bar, active metal loading same as in the laboratory‐scale unit).
Figure 15.19 Cumulative production of the pilot reactor.
Figure 15.20 Stearic acid conversion (%) as a function of time on stream.
Guide
Cover
Table of Contents
Begin Reading
Pages
iii
iv
x
xi
xii
xiii
1
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
269
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381