Distribuzione di tempi di residenza in reattori chimici

1. Distribuzione di tempi di residenza in reattori chimici

2. Introduzione • Finora abbiamo considerato solo reattori ideali come plug flow e mixed flow. • Quando si progettano reattori si cerca di avvicinarci ai reattori ideali. • Nei reattori reali i flussi deviano dall’idealità • Obiettivi di questa parte: simulare flussi non ideali • Residence Time Distribution (RTD) • Misura di RTD: metodi di input a Pulse & Step • Caratterisitche di RTD • RTD nei Reattori Ideali • Modellazione di reattori con RTD • Modello a dispersione

3. Comportamento di flusso non ideale

4. Caratterizzazione di reattori Non-Ideali Perchè ? Come ? • per modellare il comportamento di reattori non ideali • per capire come la non idealità impatta sulle performance • Residence Time Distribution (RTD)

5. Residence Time distribution - 1 • Qual’è la caratterisitca di un flusso ideale • Ogni atomo o molecola ha lo stesso tempo di residenza • Non è il caso nei reattori reali. La distribuzione di questi tempi è detta residence time distribution (RTD) • La RTD è rappresentata dalla funzione E [time-1] • La RTD può essere vista come carattristica del mixing che avviene dentro al reattore • La RTD viene normalizzata:

6. Residence Time Distribution - 2 • Frazione più giovane di t1: • Frazione più vecchia di t1:

7. Residence Time Distribution - 3 • Ragioni della deviazione dal flusso ideale: - Channeling del fluido - Creazione di zone stagnanti • La RTD può essere misurata sperimentalmente iniettando un tracciante (atomo inerte o molecola con proprietà fisiche simili a quelle del fluido) - Input ad impulso - Pulse Input - Input a gradino - Step Input

8. Effluent Feed Detector Tracer Come si Misura RTD ?? REACTOR • Caratteristiche del tracciante • facilmente rilevabile • non-reattivo • non-adsorbente • proprietà fisiche simili a quelle del fluido

9. Metodo del Pulse Input - 1 • Un ammontare noto di tracciante viene immesso nel reattore • Si misura la concenrtazione del tracciante all’uscita (la curva C) • La RTD è data da: • Notare che

10. Metodo del Pulse Input - 2 • Possibili problemi nel pulse input: • L’impulso deve avvenire in un tempo molto breve a confronto con il tempo di residenza medio. • La dispersione deve essere trascurabile tra il punto di iniezione e l’ingreso del reattore. • Se la curva C ha una lunga coda, ci può essere una certa inaccuratezza. • Infatti, una integrazione sulla curva C può essere errata anche di molto. • A volte è meglio estrapolare la curva e contiunare analiticamente.

11. Integrale di Convoluzione • La RTD è responsabile per ilcambiamento di forma traCin e Cout • Il tracciantecheentra t’ secondi prima di t è rappresentato dal rettangolo A, quindi

12. Integrale di Convoluzione • Matematicamente (stringendo il rettangolo): Che è equivalente a

13. Effluent Feed Detector Tracer C(t) C(t) C(t) C(t) t t t t Come si Misura RTD ?? REATTORE Input Response Pulse Input Method Response Input Step-Input Method

14. Input Response Pulse Input Method C(t) t t + D t C(t) t Residence-time distribution function Determinazione di RTD con Pulse Input • L’ammontare di traccianteDN che lascia il reattore tra il tempo “t” e “t+Dt” è • DividendoDN per No, si ottiene la frazione di materiale che ha un tempo di residenza tra “t” e “t+Dt”.

15.  RTD con Pulse Input (cont.) • Se Non è noto precisamente, può essere calcolato come segue: • Si può definirela funzione RTD, E(t)come: • La frazione di materiale che lascia il reattore che è rimasta nel reattore per tempi tra t1 e t2 è

16. C(t) C(t) t t Determinazione di RTD con Step Input Input Response Co Step-Input Method • Il tracciante è alimentato a concentrazione Co costante, ad iniziare da t=0 • Trascuranedo I dettagli, e sostituendo l’integrale di convoluzione, la concentrazione in uscita Cout(t) è:

17. RTD con Step Input (Cont.) • Dividendo l’equazione di sopra per Co, si ottiene la curva F(t) • La curva E(t) si ottiene differenziando F(t) rispetto al t

18. Problemi nel metodo Step Input • E’ a volte difficile da mantenere la concentrazione di tracciante ocstante in carica • Prevede derivazioni che possono portare a grossi errori (depende dalla “qualità” del segnale in uscita). Può essere necessario lisciare la curva. • E’ necessario una grosso quantitativo di tracciante. Probleme se tracciante è caro.

19. RTD da dati di Pulse Input Le concentrazioni riportate nella tabella seguente sono un responso continuo ad un impulso di ingresso in un reattore chimico. Tabulare e graficare la curva E(t).

20. RTD da dati di Pulse Input

21. Quasi plug-flow Packed Bed Reactor con channeling Caratteristiche del RTD - Esempi

22. CSTR CSTR with dead-zone & channeling Caratteristiche del RTD - Esempi

23. Relazioni integrali • F(t) è definito da Danckwerts come la funzione di distribuzione cumulativa

24. 1.0 0.8 80% del fluido spende20s o meno nel reattore F(t) 0.0 t=20 s t Caratteristiche del RTD • La frazione di materiale che risiede nel reattore per un tempo inferiore a t può essere espressa da • F(t) viene anche chiamata funzione di distribuzione cumulativa

25. Momenti di RTD • Tempo di residenza medio (Mean residence time): il valore medio è uguale al primo momento della funzione RTD, E(t): • varianza: (secondo moment): indica lo “spread” della distribuzione.

26. dove, e Caratteristiche del RTD • Mean Residence Time (tm) • Funzione RTD Normalizzata E(Q) Q è il numero di volumi di reattore basato sulle condizioni di entrata del fluido che hanno percorso il reattore nel tempo “t”

27. RTD per reattori Plug-Flow • Tutti gli atomi e le molecole spendono lo stesso tempo nel reattore. La funzione di distribuzione è quindi un picco di altezza infinita e larghezza zero. • Matematicamente, questo picco viene rappresentato dall funzione del delta di Dirac:

28. Out In 1.0 E(t) t=0 t=t t=t F(t) t t RTD nei reattori ideali • RTDs in Batch e PFRs

29. v CSTR V RTD nei reattori ideali (Cont.) • RTDs in CSTR Si può derivare RTD per un CSTR dai bilanci di materia su un tracciente inerte introdotto con un pulse input. Al tempo t=0, la concentrazione del tracciante nel CSTR è Co Input - Output + Gen = Accumulo Per t>0 Integrando l’equazione di sopra con le condizioni iniziali C=Co at t=0

30. RTD in reattori ideali (Cont.) • RTDs in CSTR Conoscendo la concentrazione del tracciante in uscita Si ottiene E(t) Quindi,

31. RTDs in reattori con flusso laminare Profilo di velocità per flusso laminare in un tubo circolare Il tempo (t) che serve ad un elemento di fluido al raggio r per attraversare un tubo di lunghezza L RTD in reattori ideali

32. r dr Portata Volumetrica di fluido che flusce tra r, er+dr  Frazione di fluido totale che scorre tra r, e r + dr RTD in reattori ideali • RTDs in reattori a flusso laminare

33.  E(t) RTD in reattori ideali • RTDs in reattori a flusso laminare (cont.) Dobbiamo riformulare eq. di sopra in termini di “t”

34. RTD in reattori ideali • RTDs Function per reattori laminari a flusso (cont.) La distribuzione cumulativa, F(t), si esprime con Il tempo di residenza medio, tm, si può derivare

35. Come possiamo utilizzare la funzione RTD per modellare un reattore ?

36. Reattore reale come serie CSTR + PFR (1) La zona agitata vicino alle eliche è modellata da un CSTR perfetto. (2) L’ingresso e uscita da un PFR. Si combinano i due modelli per il reattore reale. PFR (p) CSTR (s) E(t) Ritardo di p Inpulso tracciante t PFR (p) CSTR (s) Ritardo di p Inpulso tracciante Non ci sono differenze in RTD tra (1) CSTR + PFR e (2) PFR + CSTR !

37. Reazione del secondo ordine Una reazione del secondo ordine (CA0 = 1 kmol/m3) è condotta in un CSTR reale che vuole essere modellato da due diversi sistemi reagenti: (1) CSTR + PFR (2) PFR + CSTR Trovare la conversione in ciascun sistema nel caso che s = p = 1 min (1) CSTR + PFR CA0 CAi CA CSTR PFR bilancio massa bilancio massa

38. (1) PFR + CSTR CA0 CAi CA PFR CSTR bilancio massa bilancio massa C’è una differenza !

39. Considerazione • Reattori in serie (PFR e CSTR) per modellare zone morte di un CSTR • Importa l’ordine dei reattori? • No, in termini di RTD. • Basta che la somma dei tempi di residenza sia la stessa. • Ma… • Se la reazione è del II ordine (es 13-4 pag 834) … • … c’è una differenza • RTD non è una descrizione completa della struttura di un reattore o sistemi di reattori. • RTD è unica per un reattore • Ma un reattore non è unico per un data RTD !!!

40. Modellazione di reattori con RTD • La RTD ci dice quanto tempo molecole ed atomi sono rimaste nel reattore, ma non ci dice nulla sul miscelamento nel reattore. • Per reazioni del primo ordine, la conoscenza del tempo di stazionamento di ciascuna molecola nel reattore è sufficiente, poichè la conversione è indipendente dalla concentrazione [e.g. dX/dt = k(1-X)] • Per ordini di reazione diversi da uno, la conoscenza della RTD non è sufficiente per la previsione della conversione. Quindi è necessario sviluppare modelli che tengano conto dei fenomeni di miscelazione nel reattore.

41. Modellazione di reattori con RTD • ZERO parameteri aggiustabili • Modello a flusso segregato • Modello a massima miscelazione • UN parametero aggiustabile • Modello dei tank in serie • Modello a dispersione • DUE parameteri aggiustabili • (reattori reali come combinazione di reattori ideali)

42. Macromixing e Micromixing • Macromixing: produce una RTD ma non specifica come le molecole di età diversa si miscelano. • Quello che abbiamo fatto finora!! • La RTD è sufficiente per descrivere la non idealità • Micromixing: specifica come molecole di età diversa si miscelano. Due casi estremi di micromixing: • Tutte le molecole della stessa età rimangono assieme nel loro viaggio attraverso il reattore e non vengono miscelate fino all’uscita del reattore • Segregazione completa • Si parla di macrofluido • Molecole di diversi gruppi di età sono mescolate subito all’ingresso del reattore • Micromixing completo • Si parla di microfluido

43. Modellazione di reattori con RTD Modello a flusso segregato Si modella il reattore reale come un numero di piccoli reattori batch, ciascuno dei quali con un tempo diverso. Tutte le molecole che stazionano nel reattore per lo stesso tempo (i.e., che hanno la stessa età) rimangono assieme nello stesso globulo (i.e., reattore batch). La miscelazione dei globuli con diversa età avviene solo all’uscita el reattore reale..

44. Modello a flusso segregato

45. Modello a flusso segregato • Gli elementi di fluido di diversa età non si mescolano tra loro. • Caso del PFR: lo vediamo come un flusso di tanti reattori batch che passsano nel reattore. • Conoscere la RTD e la velocità di reazione in un modello a flusso segregato è sufficiente per calcolare la conversione Da dati sperimentali Dall’eq. di progetto di un “reattore batch” (i.e., ciascun globulo) Reazione del I ordine

46. Derivare le equazioni per una reazione del primo ordine utilizzando il modello a flusso segregato se la RTD equivale ad (a) PFR ideale, e (b) CSTR ideale. (1) PFR ideale RTD: Identiche! Reazione I ordine in PFR: (2) CSTR ideale RTD: Identiche! Reazione I ordine in CSTR:

47. Conversione media usando la RTD Per una reazione del I ordine, assumendo miscelazione completa (CSTR ideale) o segregazione completa (globuli) in un CSTR, otteniamo gli stessi risultati in termine di conversione. Questo perchè la reazione è del primo ordine; il micromixing non influenza la reazione. Per reazioni del I ordine la RTD è sufficiente per calcolare la conversione in qualsiasi tipo di reattore. Conversione media con il modello a segregazione XE t

48. Modello a massimo miscelamento - 1 • Situazionespecularerispetto a prima: entratesfalsate • Il mescolamentoavvieneall’ingresso •  è il tempo cheimpiegailfluido per muovereda un puntoparticolareall’uscita del reattore • Portatavolumetricacheentradallato del reattoretra +  e  • Portatavolumetricaalimentata al reattorefino a  • Volume di fluido con tempo di stazionamentotra  +  e  :

49. Modello a massimo miscelamento - 2 Bilancio di massa su A tra  + e: Condizioni al contorno: a  = , CA= CA0 (X = 0). Per ottenere la soluzione, le equazioni devono essere integrate indietro numericamente, iniziando da un valore di  molto grande e finendo con  = 0

50. Esempio: reazione del secondo ordine Una reazione in fase liquida del secondo ordine: La reazione è condotta a 320K e la carica è di A puro con CA0=8 mol/dm3. Il reattore è non ideale e viene modellato da due CSTRs with interscambio. Il reattore è di 1000 dm3 e la portata in carica è di 25 dm3/min. Si dispone di un test per la misura della RTD. Determinare I valori minimi e massimi di conversione per diversi gradi di micromixing del reattore. E(t) 1-F(t) C(t) I limiti sulla conversione sono dati dai due modelli di completa segregazione e di massimo miscelamento. (1) Completa segregazione E(t)