REATORES NÃO-ISOTÉRMICOS

ESTADO ESTACIONÁRIO

OPERAÇÃO ADIABÁTICA

REATORES II – AULA 06

Prof. Dra Vivian Stumpf Madeira

 Projeto de Reatores Adiabáticos
• •
Q− We = FAo *  i * C pi * (T − Tio ) + FAo * X * ( H rçp (TR ) + (T − TR ) * C p ,rç )
n

i =1

 O balanço de energia no reator, em estado estacionário, relaciona a

temperatura dos gases (ou da fase fluida) (T) com a conversão (X), o calor
e o trabalho. Se nenhum calor ou trabalho são trocados entre o sistema e a
vizinhança, a relação entre a temperatura e a conversão é:

− σ𝑛𝑖=1 𝜃𝑖 ∗ 𝐶𝑝𝑖 ∗ 𝑇 − 𝑇𝑖𝑜

𝑋(𝑇) = 𝑝 = 𝑋𝐸𝐵
Δ𝐻𝑟ç (𝑇𝑅 ) + 𝑇 − 𝑇𝑅 ∗ Δ𝐶𝑝,𝑟ç

𝑝
−𝑋 ∗ Δ𝐻𝑟ç (𝑇𝑅 ) + σ𝑛𝑖=1 𝜃𝑖 ∗ 𝐶𝑝𝑖 ∗ 𝑇𝑖𝑜 + 𝑋 ∗ Δ𝐶𝑝,𝑟ç ∗ 𝑇𝑅
𝑇(𝑋) =
σ𝑛𝑖=1 𝜃𝑖 ∗ 𝐶𝑝𝑖 + 𝑋 ∗ Δ𝐶𝑝,𝑟ç
 Relação temperatura-conversão para
Processo adiabático
 Linhas de operação
adiabática
 Para o escoamento
XBE pistonado (PFR ou PBR),
o fluido no reator se
move progressivamente
ao longo da curva [T(V)
e X(V)].
 Para o CSTR, o fluido
T imediatamente pula
para o seu valor final
σ 𝜃𝑖 𝐶𝑝𝑖 𝑇𝑜 σ 𝜃𝑖 𝐶𝑝𝑖 𝑇
𝑋𝐵𝐸 = 𝑝 − 𝑝
na curva.
Δ𝐻𝑟ç (𝑇𝑅 ) Δ𝐻𝑟ç (𝑇𝑅 )
Reator tubular – Estado estacionário - Adiabático
(FPR) – Reação única reversível

⎯⎯⎯
k1 (T )
n − bu tan o ⎯⎯ → i − bu tan o
⎯
k −1 ( T )
1A 1B

Exemplo 8-3 – Fogler: Rç de isomerização, reversível, fase líquida

(v=vo) , PFR adiabático.
 Dados:
 k1(360K)=31,1 (1/h);
 FTo=163.000 mol/h; To=330K; Hprç=-6.900,0 J/moln-butano
 Ea=65.700,0 J/mol; Kc (333K)=3,03; CAo= 9.300 mol/m³
 yn-butano,o=0,9 mol/mol; yi-pentano,o=0,1 mol/mol (inerte);
 Cpn-B = Cpi-B = 141 J/mol.K; Cpi-P = 161 J/mol.K
 VPFR para X=70% = ???

FAo (mol/s) FA (mol/s)

FIo (mol/s) FB (mol/s)
To rA FI (mol/s)
FA V FA V +V
ΔV T
Da dependência da constante de equilíbrio com a
temperatura

d ln K p H rç (T ) Kp
=  Eq. de Van ' t Hoff Com K c =
( RT )
2 
dT RT
d c b 
Se  =  + − − 1 = 0, 0; Cp = 0, 0
a a a 

C p  0,0  H rç0  1 1  
Kc (T ) = Kc (T1 )*exp  −
 R  T T  
*
  1 

ln  =
( )
 K p (T )  H rç (TR ) − TR * C p  1 1  C p
 *  −  +
T 
* ln 
 K (T ) 
 p 1  R  T1 T  R  T1 
 rç única ⎯⎯⎯
k1 (T )
n − bu tan o ⎯⎯⎯ → i − bu tan o
1) Balanço Molar : k −1 ( T )
PFR
1A 1B
dFA dX −rA
= rA ; FA = FAo *(1 − X ) → = → −rA ( X ) = ??
dV dV FAo
2) Lei de Velocidade :
 CB 
−rA = k1 (T ) *  C A −  Reação reversível
 K c (T ) 
 Ea  1 1    H rço  1 1  
k1 (T ) = k1 (T1 ) *exp  *  −   ; kc (T ) = kc (T2 ) *exp  *  − 
 R  T1 T    R  T2 T  
3) Estequiometria :
F (mol / s ) FAo *(1 − X )
CA = A = ( fase líquida ) → C A = C Ao (1 − X )
v(l / s ) vo
FB (mol / s ) FAo *( B + X )
CB = = → CB = C Ao X
v(l / s ) vo
 X 
−rA ( X , T ) = k1 (T ) * C Ao  (1 − X ) − 
 K c (T ) 
4) Combinação :
dX k1 (T )  X 
= *  (1 − X ) − 
dV vo  K c (T ) 
fç ( X ,T )

 Ea  1 1    H rço  1 1  
k1 (T ) = k1 (T1 ) *exp  *  −   ; kc (T ) = kc (T2 ) *exp  *  − 
 R  T1 T    R  T2 T  
T ( X ) = ?? ou T (V ) = ??
PFR
Adiabático
5) Balanço de Energia :
= 0,0
  i =n

 X * ( rç r ) ( r
−H (T ) + X * T * Cp rç )  o 
+

T *  i Cpi 

( )
=
est . estac . i =1
T(X ) i =n

(
 i Cpi + ( X * Cprç ) )
• •
Q =Ws =0,0

i =1

X * ( −H rç (Tr ) )
T ( X ) = To + i =n
= To + 43, 43* X
 ( Cp )
i =1
i i

Equação algébrica
 PFR adiabático - Solução Numérica EDO -
Método: Runge-Kutta
⎯⎯⎯
k1 (T )
n − bu tan o ⎯⎯⎯ → i − bu tan o
k −1 ( T ) Reação em fase líquida; Reversível; Exotérmica
1A 1B

𝑑𝑋 𝑘1 (𝑇) 𝑋
= ∗ (1 − 𝑋) − → (1)
𝑑𝑉 𝑣𝑜 𝐾𝑐 (𝑇)
𝑓ç(𝑋,𝑇)

𝐸𝑎 1 1
𝑘1 (𝑇) = 𝑘1 (𝑇1 ) ∗ exp ∗ − → (1.1)
𝑅 𝑇1 𝑇
𝑜
Δ𝐻𝑟ç 1 1
𝑘𝑐 (𝑇) = 𝑘𝑐 (𝑇2 ) ∗ exp ∗ − → (1.2)
𝑅 𝑇2 𝑇

𝑇(𝑋) = 𝑇𝑜 + 43,43 ∗ 𝑋 → (2) PFR adiabático

(eq. algébrica)
Ponto inicial:
V=0,0 − T=T0 − X=0,0
 Balanço Equações a serem resolvidas acopladas:
diferencial PFR
𝑑𝑋 −𝑟𝐴 adiabático
para o PFR = (eq.
𝑑𝑉 𝐹𝐴𝑜
adiabático diferencial)
𝑋 Uma outra
−𝑟𝐴 = 𝑘1 (𝑇)𝐶𝐴𝑜 (1 − 𝑋) −
𝐾𝑐 (𝑇) maneira de
𝐸𝑎 1 1 resolver.
Solução 𝑘1 (𝑇) = 𝑘1 (𝑇1 ) ∗ exp
𝑅
∗ −
𝑇1 𝑇
𝑜
numérica 𝑘𝑐 𝑇 = 𝑘𝑐 𝑇2 ∗ exp
Δ𝐻𝑟ç
∗
1 1
−
𝑅 𝑇2 𝑇
EDO
𝑑𝑇 𝑟𝐴 ∆𝐻𝑟ç (𝑇)
=
𝑑𝑉 𝐹𝐴𝑜 σ 𝜃𝑖 𝑐𝑝𝑖 + 𝑥∆𝑐𝑝𝑟ç
Método:
Runge-Kutta Ponto inicial:
V=0,0 − T=T0 − X=0,0
Conversão de A (n-butano) ao
longo do reator (PFR)
σ𝑛𝑖=1 𝜃𝑖 𝐶𝑝𝑖 𝑇 − 𝑇𝑜
𝑋𝐵𝐸 = 𝑝
−Δ𝐻𝑟ç (𝑇𝑅 )
X X equilíbrio
0,80
X Xeq
0,70 0,80

0,60 0,70

0,50 0,60
Conversão, X

0,40 0,50

Conversão (X)
0,30 0,40
A introdução de
0,20 0,30
inertes poderia
0,10 0,20
aumentar a
conversão!
- 0,10
0,0 1,0 2,0 3,0 4,0
Volume, m³ -
320 330 340 350 360 370

temepratura ao longo da reação/reator, T

Velocidade da reação ao longo do reator
(PFR)  X 
−rA ( X , T ) = k1 (T )* C Ao  (1 − X ) − 
70.000
 K c (T ) 
Velocidade
20
60.000 máxima
com
X35%
18 𝑑𝑋 −𝑟𝐴
50.000
=
16 𝑑𝑉 𝐹𝐴𝑜
14
-ra, mol/m³.h

40.000 𝐹𝐴𝑜
12 𝑉𝑃𝐹𝑅 = න 𝑑𝑋
−𝑟𝐴

FA0/-rA
30.000
10

20.000 8

6
10.000
4

- 2 𝐹𝐴𝑜 𝑥𝑠
- 1,00 2,00 3,00 4,00
VCSTR até 35% 𝑉𝐶𝑆𝑇𝑅 =
conversão −𝑟𝐴,𝑠
Volume, m³ -
- 0,10 0,20 0,30 0,40 0,50 0,60 0,70
Conversão, X
Perfil de temperatura (reator adiabático –
reação exotérmica)
365

360

355 X X equilíbrio
0,80
350 0,70
Temperatura, K

0,60

Conversão, X
0,50
345
0,40
0,30
340
0,20
0,10
335 -
-0,5 0,5 1,5 2,5 3,5
330 Volume, m³

325
- 0,50 1,00 1,50 2,00 2,50 3,00 3,50
Volume, m³
 Constante de equilíbrio da reação
Constante de equilíbio da rç
- Kc(T) - rç exotérmica
3,50

3,00

2,50

2,00

1,50 370

360

Temperatura, K
1,00
350

0,50 340

330
-
- 0,50 1,00 1,50 2,00 2,50 3,00 320
- 1,00 2,00 3,00
Volume, m³ Volume, m³
 Constante cinética X Constante de
equilíbrio
Constante cinética da rç - ka(T)
cte cinética cte equilíbrio

26,00

 − Ea 
k A (T ) = A *exp  
21,00
 RT 
16,00

 H rç0  1 1  
Kc (T ) = Kc (T1 )*exp  −
 R  T T  
11,00 *
  1 
6,00

1,00
- 0,50 1,00 1,50 2,00 2,50 3,00
Volume, m³
 Calor gerado na reação - Porque a
curva tem este formato?
Calor gerado na reação, (J/s)/m 3reator
450.000.000

400.000.000

350.000.000 𝑑𝑇 𝑟𝐴 ∆𝐻𝑟ç (𝑇)

=
300.000.000
𝑑𝑉 𝐹𝐴𝑜 𝑐𝑝𝐴

250.000.000

200.000.000

150.000.000

100.000.000

50.000.000

-
- 0,50 1,00 1,50 2,00 2,50 3,00 3,50 4,00 4,50
Volume PFR, m3
Reator de Mistura Perfeita (CSTR) - Adiabático - Estado
Estacionário – Exemplo 8.3 aplicado ao CSTR

𝑟ç ú𝑛𝑖𝑐𝑎
1) 𝐵𝑎𝑙𝑎𝑛ç𝑜 𝑀𝑜𝑙𝑎𝑟ቚ : FAo (mol/s)
𝐶𝑆𝑇𝑅
𝐹𝐴𝑜 − 𝐹𝐴 𝐹𝐴𝑜 − 𝐹𝐴 𝐹𝐴𝑜 ∗ 𝑋𝑠𝑎í𝑑𝑎
𝑉𝐶𝑆𝑇𝑅 = ; 𝑋= → 𝑉𝐶𝑆𝑇𝑅 =
−𝑟𝐴 𝐹𝐴𝑜 −𝑟𝐴,𝑠𝑎í𝑑𝑎
CSTR
2)𝐿𝑒𝑖 𝑑𝑒 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑𝑒:
𝑋𝑠𝑎í𝑑𝑎 -rA, mol/L.s
−𝑟𝐴,𝑠𝑎í𝑑𝑎 = 𝑘1 (𝑇𝑠𝑎í𝑑𝑎 ) ∗ 𝐶𝐴𝑜 (1 − 𝑋𝑠𝑎í𝑑𝑎 ) −
𝐾𝑐 (𝑇𝑠𝑎í𝑑𝑎 )
FA (mol/s)
3)𝐵𝑎𝑙𝑎𝑛ç𝑜 𝑑𝑒 𝐸𝑛𝑒𝑟𝑔𝑖𝑎:

𝑝
−𝑋𝑠𝑎í𝑑𝑎 ∗ Δ𝐻𝑟ç (𝑇𝑅 ) + σ𝑛𝑖=1 𝜃𝑖 ∗ 𝐶𝑝𝑖 ∗ 𝑇𝑖𝑜 + 𝑋𝑠𝑎í𝑑𝑎 ∗ Δ𝐶𝑝,𝑟ç ∗ 𝑇𝑅
𝑇𝑠𝑎í𝑑𝑎 (𝑋) =
σ𝑛𝑖=1 𝜃𝑖 ∗ 𝐶𝑝𝑖 + 𝑋𝑠𝑎í𝑑𝑎 ∗ Δ𝐶𝑝,𝑟ç
Escoamento
𝑽𝑪𝑺𝑻𝑹 𝟏
pistonado
endotérmico CSTR ou PFR 𝑭𝑨𝒐
= 𝑿
−𝒓𝑨,𝒔 𝒔
PFR
adiabáticos? 𝑽𝑷𝑭𝑹
=න
𝒙
𝟏
𝒅𝑿
𝑭𝑨𝒐 𝒐 −𝒓𝑨

Escoamento
pistonado
exotérmico

CSTR
exotérmico

Nestes casos, a associação de reatores em série

(CSTR+PFR), iria minimizar o volume total de
reação para uma mesma conversão final.
EXERCÍCIO 01: Reator tubular; Estado estacionário; Reator
adiabático; Reação única; Fase gasosa.

Suponha agora que a mesma reação de isomerização

seja realizada em fase gasosa, em um PBR, com os
parâmetros apresentados na Tabela abaixo.
Considerando-se que há perda de carga no reator (P≠Po)
e que ocorre uma variação na temperatura da fase
gasosa, (T≠To) pergunta-se:
a) Qual a conversão alcançada para 50Kg de
catalisador?
b) O que você, como eng. responsável pela planta, faria,
para aumentar a conversão para 40%, (supondo-se a
mesma massa de catalisador)?
 Dados para o k1(360K), cte cinética da rç = 31,1 1/h
Exercício 01: Energia de ativação = 65.700,0 J/mol

Kc (333,15) = cte de equilíbrio da rç = 3,03 -

-
ΔH0rç = 6.900,0 J/mol n-butano

Fto, vazão molar total da fase gasosa = 1.630,0 mol/h

Temperatura na entrada, To = 350,0 K

Pressão na entrada, Po = 10,0 atm

mol n-
Yn-butano,o = YAo = 1,0 butano/mol

Cpn-Butano = Cpi-Butano = 141,0 J/mol.K

f = porosidade do leito catalítico (permanece cte durante a rç) = 0,450
rsólido (massa específica) = 1923,00 kg/m³
diâmetro do reator (PBR) = 0,041 m
ro gás (massa específica do gás na entrada - To, Po) = 0,413 lbm/ft³
Mo (viscosidade do gás) = 0,067 lbm/ft.h
Diâmetro das partículas de catalisador, dp = 0,006 m
Conversão obtida para as condições
dadas no exercício

Conversão obtida para as condições do problema

0,80

0,70

0,60

0,50
Coinversão, X

0,40

0,30
𝐶𝑜𝑛𝑡𝑟𝑜𝑙𝑒 𝑐𝑖𝑛é𝑡𝑖𝑐𝑜
0,20

0,10

-
- 10,0 20,0 30,0 40,0 50,0 60,0
Massa de catalisador

x; dp=6mm; To=350K X eq
Queda de pressão e Vazão volumétrica do gás →
Para as condições dadas no exercício.
Perda de carga ao longo do reator
1,20

1,00

0,80 𝐸𝑙𝑒𝑣𝑎𝑑𝑎 𝑝𝑒𝑟𝑑𝑎 𝑑𝑒 𝑐𝑎𝑟𝑔𝑎

y=p/po

0,60

0,40

0,20

- Vazão volumértica da fase gasosa ao longo

- 10,00 20,00 30,00 40,00 50,00 60,00 do reator
massa catalisador, kg 25,00

20,00

15,00
𝐸𝑥𝑝𝑎𝑛𝑠ã𝑜 𝑣𝑜𝑙𝑢𝑚é𝑡𝑟𝑖𝑐𝑎 𝑑𝑜 𝑔á𝑠
10,00

5,00

-
- 10,00 20,00 30,00 40,00 50,00 60,00
Massa Kg
 Como reduzir a perda de carga no
reator?

dY  T P Reduzir o

=− * ; Y
valor de
alfa!
dW 2Y T0 Po
2 * o
  0,00;  =
 cat. * (1 -  ) * A t * Po
1,75 * G  1 −  
2
o = *  
o * Dp   ³ 
 Sugestões para aumentar a conversão:
1) Utilização de dois reatores em paralelo

Um único reator - dp=6mm - To=350K Dois reatores em paralelo – dp=6mm –

0,80 To=350K
0,80
0,70
0,70
0,60
0,60

0,50 G (cte) = 23.823,41 kg/m².h G (cte) = 11.911,70 kg/m².h

Conversão, X

0,50

Conversão, X
bo= 12,03 kPa/m bo= 3,10 kPa/m
0,40
X 0,40
a= 0,017 1/kg a= 0,004 1/kg X
X eq
0,30 X eq.
0,30

0,20 0,20
𝑋 = 20%
0,10 0,10

- -
- 10,0 20,0 30,0 40,0 50,0 - 5,0 10,0 15,0 20,0 25,0
Massa de catalisador massa catalisador, kg
2) Dois reatores em paralelo e aumento no dp (hipótese:
rç limitada pela rç qmca intrínseca)

0,80

0,70

0,60

0,50
y=p/po

0,40 X, dp=12mm
X eq.
0,30 X, dp=6mm

0,20 𝑋 = 21%

0,10

-
- 5,00 10,00 15,00 20,00 25,00
massa catalisador, kg
 2) Dois reatores em paralelo, aumento no dp,
aumento na temperatura de entrada
0,80

385
0,70

380
0,60

375
0,50
𝑋 = 40% 370

Temperatura (K)
0,40

365
0,30

360
0,20

355
0,10

350
-
0 5 10 15 20 25
massa catalisador, kg 345
- 10 20 30 40 50 60
X, dp=12mm, 360K X eq. X, dp=12mm, 350K massa catalisador, kg

2 reatores; dp 12mm 1 reator; dp 6mm

 n espécies ^

dT
−  H i * ri
= i= A
n espécies
dV
 F * cp
i= A
i i
  
dP 1 G  1 −   150* (1 −  ) *  
=− * *   *  + 1,75 *G 
dZ  Dp   ³   Dp 
   termo turbulento
termo la min ar

 Solução_Prof. Vivian
Temperatura fase fluida, K vazão volumétrica da fase fluida, v (m3/s)
510 1,40

500 1,20

490 1,00

480 0,80

0,60
470

0,40
460

0,20
450

-
440 0 10000 20000 30000 40000 50000 60000 70000 80000
0 10000 20000 30000 40000 50000 60000 70000 80000

Fração molar X Massa catalisador Perda de carga X Massa catalisador

50% 60,0

45%
50,0
40%

35%
40,0
30%

25%
YC 30,0
20%
YA
15%
20,0
10%

5% 10,0
0%
0 20000 40000 60000 80000
massa catalisador, kg -
0 10000 20000 30000 40000 50000 60000 70000 80000
 Temperatura fase fluida, K Aquecimento do reator
600
com Ta=cte=490K
500

Conversão, X(%)
400 35%

300 30%

25%
200

20%
100
15%
0
0 20000 40000 60000 80000 10%
T com aquecimento T limite T adiabático
5%

0%
0 20000 40000 60000 80000
X com aquecimento X adiabático
 Perda de carga X Massa catalisador
60,0 Adiabático com Redução na
50,0 perda de carga (dp=6mm;
40,0 porosidade=0,5)
30,0

20,0
Conversão, X(%)
40%
10,0 35%

- 30%
0 20000 40000 60000 80000 25%
dp=6mm dp=2mm 20%
15%
10%
Temperatura fase fluida, K
5%
600
0%
500 0 20000 40000 60000 80000

400 dp=6mm; alfa=5*10^-7

dp=2mm; alfa=4*10^-6
300

200

100

0
0 20000 40000 60000 80000
T_dp=6mm T limite T_dp=2mm
 Temperatura fase fluida, K
600
500 Adiabático – Aumento
400 no valor de To
300
200
100
0
0 20000 40000 60000 80000
To=450; adiabático; dp=2mm
Conversão, X
To=446; adiabático; dp=2mm 0,40

T_limite 0,35
0,30
0,25
0,20
0,15
0,10
0,05
-
0 20000 40000 60000 80000
X- To=450; adiabático; dp=2mm
X-To=446; adiabático; dp=2mm
 Reatores adiabáticos
com Estágios

Conversão de
equilíbrio

Atenção!
Reatores adiabáticos com Estágios
Reatores adiabáticos com Estágios
 Exemplo 8-6 e 8-7 Reator Adiabático 2
1,00
Xequilíbio E Xbalanço energia 0,90
0,80
Reator Adiabático 1 0,70

X (conversão)
0,60
1,00
0,50
0,90
0,40
0,80
0,30
0,70
X (conversão)

0,20
0,60
0,50 0,10
0,40 -
0,30 300 320 340 360 380 400 420 440 460
0,20 Temperatura, K
0,10
-
300 320 340 360 380 400 420 440 460 480 500 Reator Adiabático 3
Temperatura, K 1,00

0,80
σ 𝜃𝑖 𝑐𝑝𝑖 Balanço de
𝑟𝑒𝑎𝑡𝑜𝑟 𝑗
𝑋𝑠,𝑗 = 𝑇𝑠,𝑗 − 𝑇0,𝑗 0 energia, XBE 0,60
X (conversão)
∆𝐻𝑟ç (𝑇𝑟)
0,40
𝐹𝐴0,𝑗 − 𝐹𝐴,𝑗 Definição de conversão
𝑋𝑠,𝑗 = 0,20
𝐹𝐴0,𝑗 para cada reator
-
𝐶𝐵,𝑗 Constante de equilíbrio (para o 300 320 340 360 380 400 420 440 460
𝐾𝑒𝑞,𝑗 = อ (0,20)
𝐶𝐴,𝑗 cálculo de Xeq.) em cada reator Temperatura, K
𝑒𝑞.
 Exemplo 8-6 e 8-7 - Xequilíbio E Xbalanço energia
1,00

0,90 𝐹𝐴,1 = 𝐹𝐴0 1 − 𝑋1 = 𝐹𝐴0,1

𝐹𝐴,2 = 𝐹𝐴0,1 1 − 𝑋2 = 𝐹𝐴0,2
0,80 78% conversão final
𝐹𝐴,3 = 𝐹𝐴0,2 1 − 𝑋3
0,70
𝐹𝐴0 − 𝐹𝐴𝑗
0,60 𝑋𝑎𝑡é 𝑜 𝑟𝑒𝑎𝑡𝑜𝑟𝑗 =
X (conversão)

resfriamento
𝐹𝐴0
0,50
Qual será a
0,40
resfriamento conversão final
0,30 obtida se for
adicionado um
0,20
inerte no primeiro
0,10 Reator?
-
300 320 340 360 380 400 420 440 460 480 500
Temperatura, K
FIo = 10,00 molInerte/s
Xeq,R1 X BE,1 X BE,2 X BE,3 cpInerte = 25,00 cal/mol K