Progetto di reattori Non isotermi…

1. Progetto di reattori Non isotermi… Trasferimento di calore (Q) a CSTR e PFR/PBR Reattori continui a flusso in stato stazionario non isotermi (CSTR, PFR)

2. Trasferimento di calore refrigerante - reattore (da un bilancio di energia sul refrigerante) Trasferimento di calore in un CSTR Assumendo che la temperatura,T sia uniforme nel CSTR Per alte portate di refrigernate (Ta1Ta2=Ta):

3. Ta T Trasferimento di calore in un PFR/PBR Ricordando che in PFR/PBR la concentrazione e le velocità di reazione variano nel reattore. Q varierà anche. Calore totale trasferito al reattore a = area di scambio/volume Portata di calore ad una posizione nel PFR Portata di calore ad una posizione nelPBR

4. Equazioni di bilancio di energia in reattori a flusso

5. Reattori a flusso non isotermi Ora che sappiamo esprimere “Q” , torniamo ad esaminare l’equazione generale di bilancio energetico Allo stato stazionario, bilancio energia per un sistema reagente singolo Medio su T0 e T Medio su Tref e T

6. XEB T Reattori a flusso non isotermi Questa è l’equazione di bilancio energetico allo stato stazionario. Applicazioni Applicazione-1: caso speciale; reattore adiabatico senza lavoro Ws If Cp term<<DHRxn

7. Reattori a flusso non isotermi Applicazione-2: CSTR con scambio termico; no Ws Vediamo di applicare questi concetti alla soluzione di un CSTR.

9. CSTR Non-isotermo: esempio • Produzione di ossido di propilene • Reazione già vista: FA0 T = 58F T = 58F FB0 FM0 T = 75F Adiabatico V = 300 gal

10. Bilancio Moli e Cinetica Stechiometria Bilancio Energia (adiabatico)

11. Risolvendo le due equazioni simultaneamente si ottiene X = f (T)

13. Simile al problema precedente. Adesso si utilizza un serpentino di raffreddamento di 40 ft2 di superficie e con portata di acqua sufficientemente grande da considerare costante la T del fluido a 85 F. Il coefficiente globale di scambio (tipico) è di 100 Btu/h ft2 F.

14. Risolvendo le due equazioni simultaneamente si ottiene X = f (T)

15. Soluzione

16. FA0 T0 FAe Te V+DV V Reattori a flusso non isotermi Applicazione-3: PFR adiabatico Ta T Dalla 1 si ottiene T(X), quindi k(T) in funzione di X e quindi –rA in funzione di X. Bilancio energia Bilancio moli

17. Algoritmo per un reattore adiabatico • Scegli X  Calcola T  Calcola k  (se gas) calcola T°/T  calcola Kc  Calcola (-ra) • Incrementa X e ripeti i calcoli • Al termine riporta tutto in un diagramma di Levenspiel … Levenspiel Plot per una reazione esotermica, adiabatica

18. Normal butano viene isomerizzato ad isobutano in un PFR. Reazione adiabatica all’equilibrio in fase liquida sotto pressione usando tracce di catalizzatore liquido che da una velocità di reazione specifica di 3.1 h-1 a 360 K. Calcolare il volume del PFR necessario a processare 100000 gal/giorno di una miscela al 90 mol % di n-butane e 10 mol % i-pentano (inerte). La carica entra a 330 K.

19. Stechiometria Bilancio moli (PFR) Cinetica Bilancio Energia: = 0

21. Equazioni da risolvere simultaneamente X Xe V T V All’equilibrio: -rA = 0 -rA Max. Vel. Reazione nel tubo V

22. Soluzione • Calcolo a mano • Integrando: V= 2.6 m3 • Soluzione polymath V (X=.4) = 1.15 m3 • CSTR: V (X=.4) = 1.00 m3 X T -rA Xe Max. Vel. Reazione nel tubo V V V

23. FA0 T0 FAe Te V+DV V Reattori a flusso non isotermi Applicazione-4: PFR con scambio termico Ta T Differenziando l’eq. di bilancio di energia rispetto a V Valutiamo ciascun termine per ottenere un’equazione in dT/dV

24. Reattori a flusso non isotermi Applicazione-4: PFR con scambio termico (cont.) Sostituendo i termini differenziali nel bilancio Rearrangiando,

25. Reattori a flusso non isotermi Applicazione-4: PFR con scambio termico (cont.) Come si risolve il problema di un PFR non-isotermo? o Si DEVE risolvere le due equazioni differenziali, g(X,T) e f(X,T), simultanemente. Serve un Ordinary Differential Equation (ODE) solver -- Polymath

26. Torniamo indietro alla equazione generale: X non è sempre il parametro più conveniente…... Differenziando rispetto a V

27. PFR/PBR CSTR (integrando rispetto V) Usando (-rA) invece che X

28. Si può notare che per reazioni esotermiche, questo termine porta ad un aumento di T Forme alterantive dell’equazione di EB per PFR

29. Cracking in fase vapore dell’acetone per dare ketene e metano (passo principale nella produzione di anidride acetica): Reazione del I ordine rispetto all’acetone e [unità : k (s-1) eT(K)] Si alimenta 8000 kg di acetone per ora ad un reattore tubolare. Il reattore consiste di un banco di 1000 tubi da 1-inch schedule 40. Considera due casi: 1. Reattore operato adiabaticamente. 2. Reattore circonfdato da uno scambiatore di calore con coefficiente di scambio di 110 J/m2 s K, e temperatura ambiente di 1150 K. T e P di ingresso sono le stesse nei due casi: 1035 K e 162 kPa. Graficare conversione e temperatura lungo il reattore.

30. Stechiometria Bilancio moli (PFR) Cinetica E poi, bilancio di energia… per trovare X = f (T)

31. Caso 1, operazione Adiabatica Capacità termiche molto sensibili alla temperatura Dove in generale solo A a V = 0 (t = 0)

32. Tutto da risolvere simultaneamente.

33. T X Vc V Reazione endotermica adiabatica. La reazione virtualmente si spegne dopo Vc.= 3.5 m3

34. Caso 1, Operazione con scambio termico Capacità termiche molto sensibili alla temperatura Dove in generale solo A a V = 0 (t = 0)

35. a: area di scambio per unità di volume di tubo = v0 era la carica “totale” iniziale (caso adiabatico). Quando si considera lo scambio termico deve essere basata per ciascun tubo: ...

36. T X V La reazione è adiabatica endotermica ed il calore fornito dallo scambiatore fa aumentare la conversione.

37. PFR • PFR: area ombreggiata è il volume

38. CSTR • CSTR: area ombreggiata è il volume

39. CSTR + PFR • area ombreggiata è il volume