← All calculators
Beedy Academy · Calculators

Free calculators for process engineering

17 chemical process calculations: reactors, kinetics, mass balances, heat, vapour pressure and distillation. Formula shown, with practice problems.

Reactions and reactors

Conversion, yield and selectivity

FormulaX = (NA0 − NA) ÷ NA0  ·  Y = ν × NP ÷ NA0  ·  S = ν × NP ÷ (NA0 − NA)
kmol/h
kmol/h
kmol/h
Conversion of A (X)—%
Product yield (Y)—%
Selectivity (S)—%

Yield = conversion × selectivity. Falling selectivity often points to a temperature that is too high or an ageing catalyst.

Space time and conversion: CSTR and PFR

Formulaτ = V ÷ Q  ·  XCSTR = kτ ÷ (1 + kτ)  ·  XPFR = 1 − e−kτ
m³
m³/h
h⁻¹
Space time (τ)—h
Conversion in a CSTR—%
Conversion in a PFR—%

First-order reaction, constant density, isothermal. For the same volume a PFR always converts more than a CSTR.

Reactor volume for a target conversion

FormulaVCSTR = Q × X ÷ (k × (1 − X))  ·  VPFR = −Q × ln(1 − X) ÷ k
m³/h
h⁻¹
%
Volume as a CSTR—m³
Volume as a PFR—m³

First order, isothermal, constant density. Above 90 % conversion a single CSTR becomes very large: several tanks in series are used instead.

Arrhenius: effect of temperature on rate

Formulak₂ ÷ k₁ = exp[(Eₐ ÷ R) × (1 ÷ T₁ − 1 ÷ T₂)]  (T [K])
kJ/mol
°C
°C
Ratio k₂ ÷ k₁—×
Rate or constant at T₂ (k₂)—

Rule of thumb: with Eₐ ≈ 50 to 55 kJ/mol the rate roughly doubles every 10 °C. Careful with exothermic reactions: heat release follows the same curve.

Activation energy from two runs

FormulaEₐ = R × ln(k₂ ÷ k₁) ÷ (1 ÷ T₁ − 1 ÷ T₂)  (T [K])
min⁻¹
°C
min⁻¹
°C
Activation energy (Eₐ)—kJ/mol
Eₐ ÷ R—K

k₁ and k₂ in the same unit. Only two points: for a robust value, plot ln k against 1/T over at least three temperatures.

Batch cycle time and productivity

Formulatcycle = tcharge + treaction + tdischarge + tcleaning  ·  P = m × 24 ÷ tcycle
h
h
h
h
kg
Cycle time—h
Batches per day—
Productivity (P)—kg/d

Dead times (charging, discharge, cleaning) often weigh as much as the reaction itself: that is where capacity is gained.

Mass balances and blending

Mixing two streams

Formulam₃ = m₁ + m₂  ·  w₃ = (m₁ × w₁ + m₂ × w₂) ÷ m₃
kg/h
% w/w
kg/h
% w/w
Mixed flow (m₃)—kg/h
Mixed concentration (w₃)—% w/w
Solute carried—kg/h

Balance on mass, not volume: volumes do not always add up. Also works for blending two tanks (kg instead of kg/h).

Diluting a strong acid or base

Formulam₁ = m₂ × w₂ ÷ w₁  ·  V₁ = m₁ ÷ ρ₁  ·  mH₂O = m₂ − m₁
% w/w
kg/L
% w/w
kg
Concentrate needed (m₁)—kg
Concentrate needed (V₁)—L
Water needed (mH₂O)—kg

Always add acid to water, slowly and with mixing (strong heat of dilution). Work in mass: volumes do not add up.

Evaporator: water removal

FormulaP = F × xF ÷ xP  ·  W = F − P  ·  S = W ÷ e
kg/h
% w/w
% w/w
kg water/kg steam
Concentrate produced (P)—kg/h
Water evaporated (W)—kg/h
Heating steam (S)—kg/h

Solids do not evaporate: the balance sets the water to remove. Typical economy: 0.8 to 0.9 single effect, about 0.8 × n for n effects.

Process gas: density and actual flow

Formulaρ = P × M ÷ (R × T)  ·  Q = QN × (PN ÷ P) × (T ÷ TN)
kPa abs
°C
g/mol
Nm³/h
Density (ρ)—kg/m³
Actual flow—m³/h
Mass flow—kg/h

Ideal gas (Z = 1): good at low pressure; above about 10 bar correct with the compressibility factor Z. Air: M = 28.96; methane: 16.04; nitrogen: 28.01.

Heat and energy

Reaction heat and adiabatic temperature rise

FormulaQ = −ṅ × ΔHr  ·  ΔTad = −ṅ × ΔHr ÷ (ṁ × cₚ)
kmol/h
kJ/mol
kg/h
kJ/kg·K
Heat released (Q)—kW
Adiabatic rise (ΔTad)—K

ΔTad = temperature rise if all cooling is lost. For a batch use the batch kmol and kg. Above about 50 K, carry out a runaway study.

Jacketed reactor: cooling water and area

Formulaṁc = Q ÷ (cₚ × (Tc2 − Tc1))  ·  A = Q ÷ (U × ΔTlm)
kW
°C
°C
°C
kJ/kg·K
W/m²·K
Cooling water flow—kg/h
Log mean temperature difference (ΔTlm)—K
Heat transfer area needed (A)—m²

Well-mixed reactor (uniform temperature). Typical U for a water-cooled jacket: 300 to 600 W/m²·K, lower with fouling or a glass-lined vessel.

Heating a batch: energy and steam

FormulaQ = m × cₚ × (T₂ − T₁)  ·  P = Q ÷ t  ·  mv = Q ÷ Lv
kg
kJ/kg·K
°C
°C
h
kJ/kg
Energy needed (Q)—MJ
Average duty (P)—kW
Steam used (mv)—kg

Sensible heat only; losses and vessel mass ignored (add 10 to 20 %). Lv ≈ 2133 kJ/kg at 3 bar abs, 2015 kJ/kg at 10 bar abs.

Vapour-liquid equilibrium and distillation

Vapour pressure (Antoine equation)

Formulalog₁₀ Psat [mmHg] = A − B ÷ (C + T [°C])  ·  Tb = B ÷ (A − log₁₀ 760) − C
°C
Vapour pressure (Psat)—kPa
Vapour pressure (Psat)—mmHg
Normal boiling point (101.325 kPa)—°C

Constants in mmHg and °C (standard chemical engineering tables); do not extrapolate outside the range shown. Where Psat equals the vessel pressure, the liquid boils.

Relative volatility and K-values (Raoult)

FormulaKi = Pisat ÷ P  ·  α = KL ÷ KH = PLsat ÷ PHsat
°C
kPa abs
Relative volatility (α)—
K-value of light (KL)—
K-value of heavy (KH)—

Ideal mixture (Raoult), fair for similar compounds (benzene-toluene, hexane-heptane). Methanol-water and ethanol-water are non-ideal: real α differs.

Minimum number of stages (Fenske)

FormulaNmin = ln[(xD ÷ (1 − xD)) × ((1 − xB) ÷ xB)] ÷ ln α
mole fraction
mole fraction
Minimum theoretical stages (Nmin)—

Total reflux, constant α, binary mixture; Nmin includes the reboiler. In real operation about 2 to 2.5 times Nmin are needed.

Minimum reflux (Underwood, binary)

FormulaRmin = [xD ÷ xF − α × (1 − xD) ÷ (1 − xF)] ÷ (α − 1)  ·  R = f × Rmin
mole fraction
mole fraction
Minimum reflux ratio (Rmin)—
Operating reflux ratio (R)—

Saturated-liquid feed (q = 1), constant α, constant molar overflow. Usual factor f: 1.2 to 1.5.

Learning and checking tools. The ranges shown are general guidance: always follow your plant's procedures, permits and set points. Nothing is sent anywhere: the calculations run in your browser.