Teoria della penetrazione

Nello studio dei fenomeni di trasporto la teoria della penetrazione stima i coefficienti di scambio lato liquido nell'operazione unitaria di assorbimento gas-liquido, e può essere applicata anche ai casi di diffusione da gocce o bolle in quiete. Fu sviluppata nel 1935 da Ralph Higbie[1] Più in generale, la teoria della penetrazione può essere applicata ad una qualunque interfaccia gas-liquido, con un piccolo gradiente di velocità.

Secondo questa teoria, un differenziale di massa del liquido esposta per un certo tempo (chiamato appunto "tempo di esposizione") alla corrente gassosa modificherà il proprio profilo di concentrazione (ovvero la distribuzione della concentrazione di soluto), seguendo una legge del tipo:

N ¯ A = Δ c ⋅ 4 D π t e x p {\displaystyle {\bar {N}}_{A}=\Delta c\cdot {\sqrt {\frac {4{\mathcal {D}}}{\pi t_{exp}}}}}

dove:

  • N ¯ A {\displaystyle {\bar {N}}_{A}} il flusso di materia del soluto (chiamato componente A) attraverso l'interfaccia gas-liquido mediato nell'intervallo di tempo che corrisponde al tempo di esposizione
  • Δ c {\displaystyle \Delta c} la differenza di concentrazione di soluto tra la concentrazione C A i {\displaystyle C_{Ai}} all'interfaccia e la concentrazione C A 0 {\displaystyle C_{A0}} nel bulk del liquido.
  • D {\displaystyle {\mathcal {D}}} il coefficiente di diffusione di materia
  • π {\displaystyle \pi } è il cosiddetto "pi greco" (una costante pari a circa 3,14)
  • t e x p {\displaystyle t_{exp}} è il "tempo di esposizione".

Il termine 4 D π t e x p {\displaystyle {\sqrt {\frac {4{\mathcal {D}}}{\pi t_{exp}}}}} è il coefficiente di scambio di materia lato liquido KL.

Si noti che KL è proporzionale alla radice del coefficiente di diffusione di materia, mentre secondo la teoria del film KL è direttamente proporzionale al coefficiente di diffusione di materia. Confrontando questi due risultati teorici con i risultati sperimentali, si ottiene come la teoria della penetrazione predica con più accuratezza il valore di KL rispetto alla teoria del film (o del doppio film).

Una generalizzazione della teoria della penetrazione è la teoria del rinnovamento superficiale (Danckwerts, 1951).[2]

Determinazione del flusso di materia mediato nel tempo

L'equazione vista sopra è solamente il risultato finale di una teoria molto più complessa. Il procedimento analitico che ci permette di comprendere a fondo la relazione di Higbie parte dalla scrittura di un bilancio di materia per il liquido:

( S d z ) ( C A | t + d t − C A | t ) = ( S d t ) ( − D ∂ C A ∂ z | z ) − ( S d t ) ( − D ∂ C A ∂ z | z + d z ) {\displaystyle (Sdz)(C_{A}|_{t+dt}-C_{A}|_{t})=(Sdt)\left(-{\mathcal {D}}{\frac {\partial {C_{A}}}{\partial {z}}}|_{z}\right)-(Sdt)\left(-{\mathcal {D}}{\frac {\partial {C_{A}}}{\partial {z}}}|_{z+dz}\right)}

in cui:

  • z è la distanza dall'interfaccia gas-liquido
  • S è l'area della superficie di scambio di materia

Si ricava da questo bilancio la seconda legge di Fick:

∂ C A ∂ t = D ∂ 2 C A ∂ z 2 {\displaystyle {\frac {\partial {C_{A}}}{\partial {t}}}={\mathcal {D}}{\frac {\partial {^{2}C_{A}}}{\partial {z^{2}}}}}

La seconda legge di Fick va quindi particolarizzata, imponendo le condizioni al contorno del caso, che sono:

{ C A = C A 0         ( t = 0 ) C A = C A i         ( z = 0 ) C A = C A 0         ( z → ∞ ) {\displaystyle \left\{{\begin{aligned}&C_{A}=C_{A0}\ \ \ \ (t=0)\\&C_{A}=C_{Ai}\ \ \ \ (z=0)\\&C_{A}=C_{A0}\ \ \ \ (z\rightarrow \infty )\end{aligned}}\right.}

Risolvendo il sistema alle derivate parziali dato dall'equazione alle derivate parziali rappresentata dalla seconda legge di Fick e dalle sue condizioni al contorno associate, si perviene alla determinazione del flusso di materia istantaneo per z = 0 (cioè il flusso in corrispondenza dell'interfaccia gas-liquido):

N A ( t ) = Δ c ⋅ D π t {\displaystyle N_{A}(t)=\Delta c\cdot {\sqrt {\frac {\mathcal {D}}{\pi t}}}}

A questo punto bisogna mediare lungo l'intervallo di tempo tra t = 0 e t = texp per ricavare il flusso di materia mediato nel tempo:

N ¯ A = ∫ 0 t e x p N A ( t ) d t t e x p = Δ c ⋅ 4 D π t e x p {\displaystyle {\bar {N}}_{A}={\frac {\int _{0}^{t_{exp}}N_{A}(t)dt}{t_{exp}}}=\Delta c\cdot {\sqrt {\frac {4{\mathcal {D}}}{\pi t_{exp}}}}}

Note

  1. ^ Ralph Higbie, The rate of absorption of a pure gas into still liquid during short periods of exposure, in American Institute of Chemical Engineers, vol. 31, 13-15 maggio 1935, pp. 365-389.
  2. ^ In inglese "surface-renewal theory".

Bibliografia

  • Alan S. Foust, Leonard A.Wenzel; Curtis W. Clump; Luis Maus; L. Bryce Andersen, I principi delle operazioni unitarie, Ambrosiana, 1967, pp.223-236, ISBN 88-408-0117-0.
  • (EN) J. F. Richardson, J. H. Harker, Coulson & Richardson's Chemical Engineering, volume 2, 1ª ed., Butterworth-Heinemann, 2002, ISBN 0-7506-4445-1.
  • (EN) James J. Carberry, Chemical and Catalytic Reaction Engineering, Dover Publications, 2001, pp.259-260, ISBN 0-486-41736-0.
  • (EN) Warren McCabe, Julian Smith, Peter Harriott, Unit Operations In Chemical Engineering, 6ª ed., Tata Mcgraw Hill Publishers, 2005, pp. 526-528, ISBN 0-07-060082-1.

Voci correlate

Collegamenti esterni

  • Teorie per il calcolo dei coefficienti di trasporto di materia (PDF), su polymertechnology.it.
  • (EN) Penetration Theory (PDF), su crsim.utah.edu. URL consultato il 6 marzo 2009 (archiviato dall'url originale l'11 giugno 2010).
  • (EN) Correction of the penetration theory based on mass-transfer data from bubble columns operated in the homogeneous regime under high pressure, DOI:10.1016/j.ces.2007.07.030.
  • (EN) Gas-liquid reactions (PDF), su vanelk.nl.
  Portale Chimica
  Portale Fisica
  Portale Ingegneria