Effect of Tube Diameter on Plug Flow Reactor

A first-order, exothermic reaction takes place in a plug flow reactor (PFR) that has pressure drop and heat transfer through the walls. Use the slider to observe the effects of changing reactor diameter; the total feed flow rate is kept constant by changing the number of parallel reactors (number of equivalent reactors) so that the total reactor cross section does not change. Use buttons to observe how reactant molar flow rate, temperature, and pressure depend on reactor length. For small-diameter reactors, the pressure drop is higher, which increases the volumetric flow rate and reduces the residence time; this lowers conversion. Since heat transfer is more efficient for smaller-diameter reactors because the surface area per volume is larger, the temperature increases less; this also lowers conversion. 

Download the CDF file to view the simulation offline using the free Wolfram CDF player





Details
Ergun equation for pressure drop in a plug flow reactor:
where P is pressure (N/m2), z is distance down the reactor (m), ρ is fluid density (kg/m3), f is friction factor, Dh is the hydraulic diameter, which is equal to reactor diameter in absence of heating or cooling (m), v is volumetric flow rate (m3/s), Ac is cross-sectional area of the reactor (m2), FA is molar flow rate of A (mol/s), T is temperature (K), R is ideal gas constant (J/[mol K]), and the subscript 0 refers to the inlet.

The pipe is assumed to be commercial steel. For laminar flow, the following equation is used to calculate the Darcy friction factor. For turbulent and transitional flow, Serghide's explicit solution to the Colebrook equation is used:
where NRe is the Reynolds number,  ϵ is wall roughness (m),  μ is dynamic viscosity (kg/[m s]), and Dr is reactor diameter (m).

Mass balance:
where CA is the reactant concentration (mol/m3), and k is the rate law constant  (1/s).

Energy balance:
where Ua is the heat transfer coefficient time area/volume (J/[m3 s K]), Ta is the temperature of the coolant (K),  ΔH is heat of reaction (J/mol), and cP is heat capacity of reactants (J/[mol K]).