Effect of Tube Diameter on Plug Flow Reactor
A firstorder, 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 smalldiameter 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 smallerdiameter 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/m^{2}), z is distance down the reactor (m), ρ is fluid density (kg/m^{3}), f is friction factor, D_{h} is the hydraulic diameter, which is equal to reactor diameter in absence of heating or cooling (m), v is volumetric flow rate (m^{3}/s), A_{c} is crosssectional area of the reactor (m^{2}), F_{A} 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:
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 N_{Re} is the Reynolds number, ϵ is wall roughness (m), μ is dynamic viscosity (kg/[m s]), and D_{r} is reactor diameter (m).
Mass balance:
Mass balance:
where C_{A} is the reactant concentration (mol/m^{3}), and k is the rate law constant (1/s).
Energy balance:
Energy balance:
where U_{a} is the heat transfer coefficient time area/volume (J/[m^{3} s K]), T_{a} is the temperature of the coolant (K), ΔH is heat of reaction (J/mol), and c_{P} is heat capacity of reactants (J/[mol K]).