Isotherm models

Single Site Langmuir (Partial Pressure)

qi=m1b1exp(q1RT)Pi1+b1exp(q1RT)Piq_i^* = \frac { m_1 b_1 \exp( \frac {q_1}{R T} ) P_i }{ 1 + b_1 \exp( \frac {q_1}{R T} ) P_i }
SymbolDescriptionUnit
qiq_i^*Adsorbed amount of component ii at equilibriummol kg⁻¹
m1m_1Isotherm parameter 1mol kg⁻¹
b1b_1Isotherm parameter 2bar⁻¹
q1q_1Isotherm parameter 3J mol⁻¹
PiP_iPartial pressure of component iibar
TTTemperatureK
RRIdeal gas constantJ mol⁻¹ K⁻¹

Single Site Langmuir (Molar Concentration)

qi=m1b1exp(q1RT)Ci1+b1exp(q1RT)Ciq_i^* = \frac { m_1 b_1 \exp( \frac {q_1}{R T} ) C_i }{ 1 + b_1 \exp( \frac {q_1}{R T} ) C_i }
SymbolDescriptionUnit
qiq_i^*Adsorbed amount of component ii at equilibriummol kg⁻¹
m1m_1Isotherm parameter 1mol kg⁻¹
b1b_1Isotherm parameter 2mol mol⁻¹
q1q_1Isotherm parameter 3J mol⁻¹
CiC_iMolar concentration of component iimol m⁻³
TTTemperatureK
RRIdeal gas constantJ mol⁻¹ K⁻¹

Dual Site Langmuir (Partial Pressure)

qi=m1b1exp(q1RT)Pi1+b1exp(q1RT)Pi+m2b2exp(q2RT)Pi1+b2exp(q2RT)Piq_i^* = \frac { m_1 b_1 \exp( \frac {q_1}{R T} ) P_i }{ 1 + b_1 \exp( \frac {q_1}{R T} ) P_i } + \frac { m_2 b_2 \exp( \frac {q_2}{R T} ) P_i }{ 1 + b_2 \exp( \frac {q_2}{R T} ) P_i }
SymbolDescriptionUnit
qiq_i^*Adsorbed amount of component ii at equilibriummol kg⁻¹
m1m_1Isotherm parameter 1mol kg⁻¹
b1b_1Isotherm parameter 2bar⁻¹
q1q_1Isotherm parameter 3J mol⁻¹
m2m_2Isotherm parameter 4mol kg⁻¹
b2b_2Isotherm parameter 5bar⁻¹
q2q_2Isotherm parameter 6J mol⁻¹
PiP_iPartial pressure of component iibar
TTTemperatureK
RRIdeal gas constantJ mol⁻¹ K⁻¹

Dual Site Langmuir (Molar Concentration)

qi=m1b1exp(q1RT)Ci1+b1exp(q1RT)Ci+m2b2exp(q2RT)Ci1+b2exp(q2RT)Ciq_i^* = \frac { m_1 b_1 \exp( \frac {q_1}{R T} ) C_i }{ 1 + b_1 \exp( \frac {q_1}{R T} ) C_i } + \frac { m_2 b_2 \exp( \frac {q_2}{R T} ) C_i }{ 1 + b_2 \exp( \frac {q_2}{R T} ) C_i }
SymbolDescriptionUnit
qiq_i^*Adsorbed amount of component ii at equilibriummol kg⁻¹
m1m_1Isotherm parameter 1mol kg⁻¹
b1b_1Isotherm parameter 2mol mol⁻¹
q1q_1Isotherm parameter 3J mol⁻¹
m2m_2Isotherm parameter 4mol kg⁻¹
b2b_2Isotherm parameter 5mol mol⁻¹
q2q_2Isotherm parameter 6J mol⁻¹
CiC_iMolar concentration of component iimol m⁻³
TTTemperatureK
RRIdeal gas constantJ mol⁻¹ K⁻¹

Henry's Law (Partial Pressure)

qi=kexp(QRT)Piq_i^* = k \exp( \frac {Q}{R T} ) P_i
SymbolDescriptionUnit
qiq_i^*Adsorbed amount of component ii at equilibriummol kg⁻¹
kkIsotherm parameter 1mol kg⁻¹
QQIsotherm parameter 2J mol⁻¹
PiP_iPartial pressure of component iibar
TTTemperatureK
RRIdeal gas constantJ mol⁻¹ K⁻¹

Rate models

Linear driving force

qit=k(qiqi)\frac{\partial q_i}{\partial t} = k (q_i^* - q_i)
SymbolDescriptionUnit
kkMass transfer coefficients⁻¹
qiq_i^*Equilibrium adsorbed concentrationmol kg⁻¹

Macropore Diffusion

The macropore mass transfer coefficient is calculated using the following equation:

k=60Deffdp2k = 60 \frac{D_{\text{eff}}}{d_p^2}

where the effective diffusivity DeffD_{\text{eff}} combines Knudsen diffusivity DkD_k and molecular diffusivity DmD_m through:

1Deff=1Dk,eff+1Dm,eff\frac{1}{D_{\text{eff}}} = \frac{1}{D_{k,\text{eff}}} + \frac{1}{D_{m,\text{eff}}}

The effective Knudsen and molecular diffusivities account for the porosity and tortuosity of the medium:

Dk,eff=DkϵτD_{k,\text{eff}} = D_k \frac{\epsilon}{\tau} Dm,eff=DmϵτD_{m,\text{eff}} = D_m \frac{\epsilon}{\tau}

The Knudsen diffusivity DkD_k is estimated as:

Dk=97rpTMD_k = 97 \, r_p \sqrt{\frac{T}{M}}
SymbolDescriptionUnit
kkMass transfer coefficients⁻¹
DeffD_{\text{eff}}Effective diffusivitym²/s
Dk,effD_{k,\text{eff}}Effective Knudsen diffusivitym²/s
Dm,effD_{m,\text{eff}}Effective molecular diffusivitym²/s
DkD_kKnudsen diffusivitym²/s
DmD_mMolecular diffusivity (input)m²/s
ϵ\epsilonPorosity of the medium
τ\tauTortuosity of the medium
dpd_pParticle diameterm
rpr_pPore radiusm
TTTemperatureK
MMMolecular weightg/mol