Energy Balance

Skarstrom offers two thermal models, selected in Column Properties → Energy Balance via the Thermal equilibrium model drop-down.

Local Thermal Equilibrium (LTE)

This is the default. It assumes heat transfer between the gas and solid is fast enough that the two share a single temperature TT, so one lumped energy equation is solved per cell.

(ϵbCp,gρg+Cp,sρb)Ttλ2Tz2ϵbPt+uCp,gρgTz+4hidi(TTw)RTPρbi=1N(ΔHi)qit=0\left(\epsilon_b C_{p,g} \rho_g + C_{p,s} \rho_b\right) \frac{\partial T}{\partial t} - \lambda \frac{\partial^2 T}{\partial z^2} - \epsilon_b \frac{\partial P}{\partial t} + u C_{p,g} \rho_g \frac{\partial T}{\partial z} + \frac{4 h_i}{d_i} (T - T_w) - \frac{RT}{P}\rho_b\sum_{i=1}^{N}(-\Delta H_i)\frac{\partial{q_i}}{\partial{t}} = 0
SymbolDescriptionUnit
Cp,gC_{p,g}Gas heat capacityJ kg⁻¹ K⁻¹
Cp,sC_{p,s}Adsorbent heat capacityJ kg⁻¹ K⁻¹
TTGas/Adsorbent temperatureK
TwT_wWall temperatureK
PPPressurePa
uuSuperficial velocitym s⁻¹
ϵb\epsilon_bBed porosity (gas phase only)-
ρg\rho_gGas densitykg m⁻³
ρb\rho_bAdsorbent bulk densitykg m⁻³
λ\lambdaEffective axial thermal conductivityW m⁻¹ K⁻¹
hih_iInternal heat transfer coefficient (gas–wall)W m⁻² K⁻¹
did_iWall internal diameterm
qiq_iComponent i mass sourcemol kg⁻¹ s⁻¹
ΔHi\Delta H_iComponent i heat of adsorptionJ mol⁻¹

The final term is the heat released (or absorbed) by adsorption, using the per-component isosteric heat of adsorption ΔHi\Delta H_i entered on the Adsorbent Material → Heat of Adsorption tab. The wall-coupling term 4hidi(TTw)\tfrac{4 h_i}{d_i}(T - T_w) links this equation to the wall energy balance.

The effective axial thermal conductivity λ\lambda and the internal coefficient hih_i are either constants or computed from the transport correlations.

Local Thermal Non-Equilibrium (LTNE)

When gas–solid heat transfer is not fast (for example large particles or fast cycles), select Local thermal non-equilibrium (LTNE). The solver then carries separate fluid and solid temperature equations coupled by a gas–particle heat transfer term. This requires a Particle heat transfer model and coefficient, described in the transport correlations. With LTNE, the solid temperature is available as its own plotted variable.

Adsorbed-phase heat capacity

The Adsorbed phase heat capacity selector (same tab) controls how the adsorbed molecules' heat capacity is treated:

  • Constant — enter a molar heat capacity per component.
  • Gas phase — the adsorbed-phase heat capacity is taken equal to the gas-phase value.

Energy source (heating / cooling)

Any step can apply a heat duty to a layer — the basis of temperature swing (TSA) operation. On the Step Data → Energy Sources tab, tick a layer and enter a Power [W]; positive values heat the layer, negative values cool it. The duty is distributed over the layer and added as a source term to the energy balance for the duration of that step.