La rilevanza del processo di desolforazione in situ in combustori a letto fluidizzato è stata già sottolineata in una lezione del Corso di Processi e Impianti di Trattamento Reflui, presente su federic@, e dedicata alle applicazioni industriali dei reattori a letto fluidizzato.
In questa lezione, si vuole presentare il modello di solfatazione “a grani”, che rappresenta un’estensione del modello nucleo-mantello visto nelle lezioni precedenti.
La reazione di interesse è la solfatazione del sorbente calcareo, che prevede il coinvolgimento, in un processo eterogeneo, di CaO (derivante dalla veloce calcinazione del calcare) ed SO2 (prodotto dalla combustione), in presenza di O2, per dare CaSO4.
La particella di CaO si assume sferica, di raggio R costante. Essa ha una porosità ε, mentre il grado di pieno 1-ε tiene conto della presenza di grani (sub-particelle) di CaO.
La concentrazione di SO2 vale CSO2,g all’esterno della particella (bulk gassoso, valori tipici intorno a 2000 ppm, come prodotto dalla combustione), mentre all’interno della particella essa è funzione della coordinata radiale r e del tempo t.
La Figura mostra, in sezione, la particella di CaO considerata.
Come mostrato in Figura, ogni grano si solfata seguendo lo Shrinking Core Model. Esso cioè ha raggio costante rg, ed è costituito da un mantello di prodotto (CaSO4) che circonda un nucleo di reagente da convertire (CaO), di raggio rc variabile nel tempo.
Il grano, al proprio esterno, vede la concentrazione CSO2(r,t), ossia quella interna alla particella complessiva. La concentrazione di SO2 è una funzione, detta CSO2,shell, nel mantello del grano, ed è una funzione, detta CSO2,i, all’interfacies core-shell.
Quindi, mentre nello SCM classico si assume che la reazione chimica avvenga all’interfacies core-shell dell’intera particella, qui la trattazione è più raffinata, poiché si considera la reazione chimica all’interfacies core-shell di ogni singolo grano costituente una particella.
Il grado di conversione locale, ossia riferito al singolo grano di una particella di CaO, è definito (in accordo allo SCM che caratterizza il grano) come:
E’ possibile, mediante media volumetrica, definire ora il grado di conversione medio (“average”) per l’intera particella:
Si nota quindi, per usare la (2), che è necessario conoscere la funzione rc(r), cosa che verrà ottenuta mediante l’applicazione di questo modello. Ovviamente, ci si aspetta che quando r diminuisce (verso il centro della particella), rc aumenti (nuclei meno convertiti, perché più lontani dall’atmosfera gassosa).
La velocità di solfatazione, espressa come moli reagenti per unità di tempo e di superficie reattiva (quella di interfaccia), è:
dove k è un’opportuna costante cinetica, espressa in lunghezza/tempo. Tenendo presente il rapporto tra superficie reattiva e volume del singolo grano:
si può scrivere la velocità di solfatazione anche come moli reagenti per unità di tempo e di volume del grano:
Considerando il sistema pseudo-stazionario (aspetto già presentato a proposito dello SCM), scegliendo come volume di controllo quello compreso tra coordinata radiale r ed r+dr, ricordando (nel termine di Fick) che SO2 si muove per diffusione in verso contrario alla definizione di r, e tenendo presente che la frazione di particella occupata da grani vale 1-ε, si può scrivere il seguente bilancio di materia:
nel quale è stata introdotta l’area specifica corretta a:
Questo bilancio tiene conto di termini di ingresso ed uscita per diffusione di SO2 all’interno della particella, e di scomparsa per reazione chimica all’interfacies nucleo-mantello di ogni singolo grano.
La (5) può essere matematicamente riarrangiata in:
ossia, sviluppando la derivata del prodotto:
La (8) assume la formulazione finale:
che rappresenta l’equazione di bilancio di materia cercata. Essa è differenziale, alle derivate parziali, e di secondo ordine. Necessita quindi di due condizioni ai limiti. La prima, semplice, è che CSO2=CSO2,g per r=R.
Per la seconda, si considera che il decorso della reazione sia trascurabile per coordinate radiali molto vicine al centro della particella. Pertanto, per r=0 si può considerare nullo il gradiente, ed imporre la seguente condizione (detta “alla Neumann”, perché imposta non sulla funzione, ma sulla sua derivata):
La soluzione dell’equazione (9) restituisce CSO2(r,rc,CSO2,i).
Il bilancio di materia su SO2 circoscritto a valori di r appartenenti all’intervallo [rc,rg], ovvero al solo mantello di un singolo grano, si scrive in maniera del tutto analoga a quanto fatto per l’intera particella, tenendo però in conto dei soli fenomeni diffusivi, poiché nel mantello di CaSO4 non può esserci reazione chimica.
Quindi, si riscrive la (9) tenendo conto di quanto appena detto, e considerando naturalmente che adesso la variabile è CSO2,shell:
La (10) necessita di due condizioni ai limiti. La prima, semplice, è che CSO2,shell=CSO2 (la concentrazione esterna al grano ed interna alla particella, ottenuta come soluzione dal bilancio precedente (9)) per r=rg.
Per la seconda, si impone che il flusso che arriva all’interfacies prenda parte alla reazione chimica, ovvero, per r=rc:
La soluzione della (10) restituisce CSO2,shell(r,rc,CSO2,CSO2,i). Una soluzione particolare, e di interesse, può essere cercata per r=rc, dove CSO2,shell=CSO2,i: è una soluzione, quindi, del tipo CSO2,i(rc,CSO2).
Poiché la soluzione della (9) era, come visto, CSO2(r,rc,CSO2,i), la combinazione di queste due ultime soluzioni rende CSO2(r,rc). E’ chiaro quindi che, una volta trovata la relazione rc(t), sarà possibile esprimere il profilo di concentrazione di SO2 in funzione della coordinata radiale della particella e del tempo, CSO2(r,t).
Il bilancio su CaO (per n numero di moli) si scrive:
dove:
per CCaO concentrazione di CaO nel nucleo (a tutti gli effetti, è la densità molare di CaO, poiché il nucleo è qui esclusivamente costituito da CaO).
Combinando, è:
dove si è ricordata la funzionalità, già trovata, per CSO2,i.
Imponendo la condizione iniziale rc=rg (t=0), la (13) restituisce la funzionalità rc(t,CSO2). Essa, combinata con la già trovata CSO2(r,rc), consente di determinare il desiderato profilo CSO2(r,t).
La procedura operativa può adesso essere schematizzata nel seguente modo:
1. considerata una particella di determinato raggio R, si fissa un tempo t;
2. si fissa una coordinata radiale r all’interno della particella;
3. dalla soluzione CSO2(r,t), si determina CSO2;
4. dalla soluzione rc(t,CSO2), si determina rc.
Iterando la procedura per ogni r, si ottiene il profilo rc(r). Esso è quello necessario, come detto, per risolvere l’integrale presente nella (2), ossia per ottenere il grado di conversione medio della particella.
L’approccio mostrato vale per un particolare istante t. Estendendo nel tempo (cioè, iterando la procedura per ogni istante), e considerando (come visto altrove) la distribuzione dei tempi di residenza E(t) delle particelle di sorbente aventi lo stesso raggio R, si può ottenere il grado di conversione medio per una determinata classe granulometrica:
In conclusione, è possibile (come visto altrove) ottenere il grado di solfatazione dell’intero campione di sorbente, mediandolo sulle diverse classi granulometriche mediante y(Ri), la frazione volumetrica di particelle aventi raggio Ri:
1. Carbone e biomasse: generalità
2. Processi reattivi del char: aspetti cinetici e diffusivi - parte prima
3. Processi reattivi del char: aspetti cinetici e diffusivi - parte seconda
4. Processi reattivi del char: aspetti reattoristici
5. Combustione: aspetti cinetici e bilanci di materia
6. Desolforazione in situ in combustori a letto fluidizzato - parte prima (modello a grani)
7. Desolforazione in situ in combustori a letto fluidizzato - parte seconda (bilancio di popolazione)
8. Gassificazione - parte prima
9. Gassificazione - parte seconda
10. Rimozione di particolato da reflui gassosi: applicazioni esercitative - parte prima
11. Rimozione di particolato da reflui gassosi: applicazioni esercitative - parte seconda
12. Assorbimento con reazione chimica: applicazioni esercitative
13. Adsorbimento: modellazione frattale di dati cinetici
14. Evoluzione dinamica di processi chimici: trasformate di Laplace per funzioni di interesse
15. Evoluzione dinamica di processi chimici: applicazioni esercitative
16. Reattori chimici: integrazioni esercitative
17. Analisi granulometrica e porosimetrica: applicazioni esercitative