Overall Mass Balance

The overall mass balance ensures mass conservation throughout the packed bed.

ϵt(1PPt1TTt)+TPz(uPT)+RTPρbi=1Nqit=0\epsilon_t(\frac{1}{P}\frac{\partial P}{\partial t} - \frac{1}{T}\frac{\partial T}{\partial t}) + \frac{T}{P}\frac{\partial}{\partial z}(\frac{u P}{T}) + \frac{RT}{P}\rho_b\sum_{i=1}^{N}\frac{\partial{q_i}}{\partial{t}} = 0
SymbolDescriptionUnit
PPPressurePa
TTTemperatureK
ϵt\epsilon_tTotal void fraction-
uuSuperficial velocitym s⁻¹
qiq_iComponent i mass sourcemol kg⁻¹ s⁻¹
ttTimes
zzAxial coordinatem
RRIdeal gas constantJ mol⁻¹ K⁻¹

Component Mass Balance

The component mass balance accounts for axial dispersion, convection, and adsorption/desorption accordnig to the rate model selected. Axial dispersion can be constant or calculated from different empirical equations.

ϵt(yit+yiPPtyiTTt)+TPz(uyiPT)ϵbTPz(DaxPTyiz)+RTPρbqit=0\epsilon_t (\frac{\partial y_i}{\partial t} + \frac{y_i}{P}\frac{\partial P}{\partial t} - \frac{y_i}{T}\frac{\partial T}{\partial t}) + \frac{T}{P} \frac{\partial}{\partial z}(\frac{u y_i P}{T}) - \epsilon_b \frac{T}{P}\frac{\partial}{\partial z}(\frac{D_{ax} P}{T} \frac{\partial y_i}{\partial z}) + \frac{RT}{P} \rho_b \frac{\partial q_i}{\partial t} = 0
SymbolDescriptionUnit
yiy_iMole fraction of component ii in gas phase-
PPPressurePa
TTTemperatureK
ϵb\epsilon_bBed void fraction-
uuSuperficial velocitym s⁻¹
DaxD_{ax}Axial disperion coefficientm² s⁻¹
qiq_iComponent i mass sourcemol kg⁻¹ s⁻¹
ttTimes
zzAxial coordinatem
RRIdeal gas constantJ mol⁻¹ K⁻¹