MaximaLabs
All unit operations

Kinetic reactor simulation

1 inlet -> 1 outlet (CSTR/PFR)

Governing equations

The exact equations the solver works for a kinetic reactor — the same math shown in the app's "Theory" panel, not a black box.

r=kf ⁣react ⁣Cnkr ⁣prod ⁣Cn    (reversible),k=k0eEa/RTr = k_f\!\prod_{react}\! C^{n} - k_r\!\prod_{prod}\! C^{n}\;\;(\text{reversible}),\quad k=k_0 e^{-E_a/RT}
rLHHW=r/(1+iKiCi)m    (adsorption inhibition)r_{LHHW} = r\,/\,(1 + \textstyle\sum_i K_i C_i)^m\;\;(\text{adsorption inhibition})
CSTR: ξ=r(ξ)VPFR: dn˙idV=νir\text{CSTR: } \xi = r(\xi)\,V \qquad \text{PFR: } \frac{d\dot n_i}{dV} = \nu_i\,r
dTdV=r(ΔHr)rr(UA/V)(TTc)in˙icp,i    (non-isothermal)\frac{dT}{dV} = \frac{\sum_r(-\Delta H_r)\,r_r - (UA/V)(T-T_c)}{\sum_i \dot n_i\,c_{p,i}}\;\;(\text{non-isothermal})
rr
reaction rate [mol/m³/s]
k0k_0
pre-exponential factor
EaE_a
activation energy [J/mol] — the temperature sensitivity of the reaction
RR
gas constant
TT
temperature [K]
CjC_j
concentration of species j [mol/m³]
njn_j
reaction order in j
ξ\xi
reaction extent [mol/s]
VV
reactor volume [m³]
n˙i\dot n_i
component molar flow [mol/s]
νi\nu_i
stoichiometric coefficient (negative for reactants)
kf,krk_f, k_r
forward / reverse Arrhenius rate constants (reversible reaction)
KiK_i
LHHW adsorption constant of species i
mm
adsorption-term exponent
ΔHr\Delta H_r
heat of reaction r [J/mol] (− exothermic)
UAUA
overall heat-transfer coefficient × area to the coolant [W/K]
TcT_c
coolant temperature [K]
cp,ic_{p,i}
molar heat capacity [J/mol/K]

Parameters

mode 'cstr'|'pfr', volume [m^3], k0 [Arrhenius pre-exp], Ea [J/mol], stoichiometry, key_component, optional orders {component: order}, temperature [K], heat_of_reaction [J/mol]; reversible: k0_rev/Ea_rev (+ orders_rev); LHHW: adsorption {component: K_i} + lhhw_exponent; energy balance: adiabatic=true OR ua [W/K] + coolant_temperature [K]

Example flowsheets that use it

Stop fighting legacy software. Build your first flowsheet in 60 seconds.