1.3.
REACTORES TUBULARES PARA REACCIONES
HOMOGÉNEAS EN FLUJO LAMINAR Y EN FLUJO
TURBULENTO.
REACTORES TUBULARES EN FLUJO LAMINAR
Perfil de velocidad: líneas de corriente de elevada
velocidad en el centro y velocidad casi nula en las
paredes
Gradiente de composición con el radio compensado
por la difusión molecular
D: 10-9 m2/s – 10-10 m2/s
𝐷𝐷.𝑡𝑡̅
• Difusión molecular despreciable < 3. 10−3
𝑅𝑅 2
Para una línea de velocidad,
𝑑𝑑𝑑𝑑
𝑣𝑣𝑧𝑧 (𝑟𝑟) = 𝑟𝑟𝐴𝐴 , C(r, 0) = Cin
𝑑𝑑𝑑𝑑
C(r, L) = Cout(r)
1
Integrando,
−𝑘𝑘. 𝑧𝑧
𝐶𝐶 (𝑟𝑟, 𝑧𝑧) = 𝐶𝐶𝑖𝑖𝑖𝑖 . exp � �
𝑣𝑣𝑧𝑧 (𝑟𝑟)
A la salida del reactor
𝐶𝐶 (𝑟𝑟, 𝐿𝐿) = 𝐶𝐶𝑜𝑜𝑜𝑜𝑜𝑜 (𝑟𝑟) =
−𝑘𝑘. 𝐿𝐿
= 𝐶𝐶𝑖𝑖𝑖𝑖 . exp � � = 𝐶𝐶𝑖𝑖𝑖𝑖 . exp(−𝑘𝑘. 𝑡𝑡)
𝑣𝑣𝑧𝑧 (𝑟𝑟)
𝐿𝐿
𝑡𝑡 = 𝑣𝑣 , tiempo de residencia para una línea de corriente
𝑧𝑧 (𝑟𝑟)
Propiedades físicas constantes,
𝑟𝑟 𝑛𝑛+1
Ec. de Kouette/Ostwald-de-Waele 𝑣𝑣𝑧𝑧 (𝑟𝑟) = 𝑢𝑢�. �1 − �𝑅𝑅� �
𝑟𝑟 2 𝑟𝑟 2
𝑣𝑣𝑧𝑧 (𝑟𝑟) = 2𝑢𝑢�. �1 − 𝑅𝑅2 � = 𝑢𝑢𝑚𝑚𝑚𝑚𝑚𝑚 �1 − 𝑅𝑅2 �
Introduciendo en C(r,L),
−𝑘𝑘. 𝐿𝐿
𝐶𝐶𝑜𝑜𝑜𝑜𝑜𝑜 (𝑟𝑟) = 𝐶𝐶𝑖𝑖𝑖𝑖 . exp � �=
𝑟𝑟 2
𝑣𝑣𝑚𝑚𝑚𝑚𝑚𝑚 �1 − � 2 ��
𝑅𝑅
Concentración promedio a la salida “mixing cup”:
1 𝑅𝑅
̅
𝐶𝐶𝑜𝑜𝑜𝑜𝑜𝑜 = � 𝑐𝑐𝑜𝑜𝑜𝑜𝑜𝑜 (𝑟𝑟). 𝑣𝑣𝑧𝑧 (𝑟𝑟). 2. 𝜋𝜋. 𝑟𝑟. 𝑑𝑑𝑑𝑑
𝑄𝑄 0
𝑅𝑅
𝑄𝑄 = � 𝑣𝑣𝑧𝑧 (𝑟𝑟). 2. 𝜋𝜋. 𝑟𝑟. 𝑑𝑑𝑑𝑑 = 𝜋𝜋. 𝑅𝑅 2 . 𝑢𝑢�
0
2
• Difusión molecular y convección:
𝜕𝜕𝜕𝜕 1 𝜕𝜕𝜕𝜕 𝜕𝜕 2 𝐶𝐶 𝜕𝜕 2 𝐶𝐶
𝑣𝑣𝑧𝑧 (𝑟𝑟) = 𝐷𝐷𝐴𝐴 � + + � + 𝑟𝑟𝐴𝐴
𝜕𝜕𝜕𝜕 𝑟𝑟 𝜕𝜕𝜕𝜕 𝜕𝜕𝑟𝑟 2 𝜕𝜕𝑧𝑧 2
Difusión de Fick en dirección axial y radial
𝑧𝑧 𝑟𝑟
Variables adimensionales 𝑧𝑧 = 𝐿𝐿 𝑟𝑟 = 𝑅𝑅
𝑣𝑣𝑧𝑧 (𝑟𝑟) 𝜕𝜕𝜕𝜕 𝐷𝐷𝐴𝐴 . 𝑡𝑡̅ 1 𝜕𝜕𝜕𝜕 𝜕𝜕 2 𝐶𝐶 𝑅𝑅 2 𝜕𝜕 2 𝐶𝐶
= 2 � + +� � � 𝑟𝑟𝐴𝐴
� + 𝑡𝑡.
𝑢𝑢� 𝜕𝜕𝜕𝜕 𝑅𝑅 𝑟𝑟 𝜕𝜕𝜕𝜕 𝜕𝜕𝑟𝑟 2 𝐿𝐿 𝜕𝜕𝑧𝑧 2
𝐷𝐷𝐴𝐴 .𝑡𝑡̅
𝑅𝑅2
importancia de la difusión molecular
𝑅𝑅 2
� 𝐿𝐿 � importancia de la difusión radial comparada con la
axial
(R/L)2 < 10-2
Reactores normales
(R/L)2 < 10-4
Condiciones de contorno,
C = Cin z=0
𝜕𝜕𝜕𝜕
𝜕𝜕𝜕𝜕
=0 r=0
𝜕𝜕𝜕𝜕
𝜕𝜕𝜕𝜕
=0 r=1
Solución del balance de materia C(r,z) para cada
componente
Concentración a la salida, C(r,1) = Cout(r) concentración
promedio
3
Perfiles de temperatura en régimen laminar:
Balance de energía: convección, conducción y liberación de
calor
𝑣𝑣𝑧𝑧 𝜕𝜕𝜕𝜕 𝛼𝛼 𝑇𝑇 . 𝑡𝑡̅ 1 𝜕𝜕𝜕𝜕 𝜕𝜕 2 𝑇𝑇 ∆𝐻𝐻𝑅𝑅 . 𝑟𝑟𝐴𝐴 . 𝑡𝑡̅
= 2 � + �−
𝑢𝑢� 𝜕𝜕𝜕𝜕 𝑅𝑅 𝑟𝑟 𝜕𝜕𝜕𝜕 𝜕𝜕𝑟𝑟 2 𝜌𝜌. 𝑐𝑐𝑝𝑝
𝛼𝛼 𝑇𝑇 difusividad térmica
𝑅𝑅 2
� 𝐿𝐿 � ↓↓ conducción axial despreciable
Condiciones de contorno:
1) T = Tin z=0 , condición de entrada
𝜕𝜕𝜕𝜕
2) 𝜕𝜕𝜕𝜕 = 0 r=0 , simetría radial
𝜕𝜕𝜕𝜕
3) 𝜕𝜕𝜕𝜕 = 0 r=1 , operación adiabática
T = Twall r=1 , Tª constante en la pared
T = Twall (z) r = 1 , perfil específico
Solución del balance de energía T(r,z)
Temperatura a la salida, T(r,1) = Tout(r)
1
𝑣𝑣𝑧𝑧 (𝑟𝑟)
𝑇𝑇 = 2 � 𝑇𝑇𝑜𝑜𝑜𝑜𝑜𝑜 (𝑟𝑟) 𝑟𝑟. 𝑑𝑑𝑑𝑑
0 𝑢𝑢�
4
Coordenadas rectangulares:
Platos paralelos muy anchos, balance de materia,
𝑣𝑣𝑧𝑧 𝜕𝜕𝜕𝜕 𝐷𝐷𝐴𝐴 . 𝑡𝑡̅ 𝜕𝜕 2 𝐶𝐶
= 2 � 2 � + 𝑡𝑡. � 𝑟𝑟𝐴𝐴
𝑢𝑢� 𝜕𝜕𝜕𝜕 𝑅𝑅 𝜕𝜕𝑦𝑦
y=R y=1
2R y=0 y=0
y = -R y = -1
• Perfil de velocidad para fluido newtoniano, con
viscosidad constante,
𝑦𝑦 2
𝑣𝑣𝑧𝑧 (𝑦𝑦) = 1,5 . 𝑢𝑢� �1 − 2 � = 1,5 . 𝑢𝑢� . (1 − 𝑦𝑦 2 )
𝑅𝑅
Condiciones de contorno,
C = Cin z=0
𝜕𝜕𝜕𝜕
𝜕𝜕𝜕𝜕
=0 y=0
𝜕𝜕𝜕𝜕
𝜕𝜕𝜕𝜕
=0 y=1
Solución del balance de materia C(z,y)
5
• Perfil lineal de velocidad,
𝑦𝑦
𝑣𝑣𝑧𝑧 (𝑦𝑦) = 𝑣𝑣0 .
𝐻𝐻
𝑣𝑣𝑧𝑧 𝜕𝜕𝜕𝜕 𝐷𝐷𝐴𝐴 . 𝑡𝑡̅ 𝜕𝜕 2 𝐶𝐶
= 2 � 2 � + 𝑡𝑡. � 𝑟𝑟𝐴𝐴
𝑢𝑢� 𝜕𝜕𝜕𝜕 𝑅𝑅 𝜕𝜕𝑦𝑦
Condiciones de contorno aproximadas,
𝜕𝜕𝜕𝜕
𝜕𝜕𝜕𝜕
=0 y = 0, H
Transmisión de calor,
T = Tupper y=H
T = Tlower y=0
Viscosidad variable,
µ = µ (T, C)
Proximidad a la pared: si hay calentamiento µ ↓↓ perfil
de velocidad se aplana se aproxima a PFR
enfriamiento µ ↑↑ perfil v(r) se elonga
6
Para Re < 100 en reacciones con polímero
i) El momento del fluido es despreciable comparado con
las fuerzas viscosas
ii) 𝑣𝑣𝑟𝑟 es despreciable comparado con 𝑣𝑣𝑧𝑧
𝑑𝑑𝑑𝑑 1 𝑑𝑑 𝑑𝑑𝑣𝑣𝑧𝑧
0= − + �𝜇𝜇 . 𝑟𝑟 �
𝑑𝑑𝑑𝑑 𝑟𝑟 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝑑𝑑𝑣𝑣𝑧𝑧
=0 𝑟𝑟 = 0 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆í𝑎𝑎 𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟
𝑑𝑑𝑑𝑑
𝑣𝑣𝑧𝑧 = 0 𝑟𝑟 = 𝑅𝑅
Solución,
𝑅𝑅 𝑟𝑟
1 −𝑑𝑑𝑑𝑑 𝑟𝑟 𝑑𝑑𝑑𝑑 𝑅𝑅 ∫𝑟𝑟 𝜇𝜇
𝑑𝑑𝑑𝑑
𝑣𝑣𝑧𝑧 = � �� = 𝑢𝑢0 𝑅𝑅 𝑟𝑟
2 𝑑𝑑𝑑𝑑 𝑟𝑟 𝜇𝜇 ∫0 𝜇𝜇 𝑑𝑑𝑑𝑑
𝑑𝑑𝑑𝑑
𝑢𝑢0 𝑦𝑦 𝑑𝑑𝑑𝑑
constantes locales que se determinan a partir de
los caudales másicos:
𝑅𝑅
� 2. 𝜋𝜋. 𝑟𝑟. 𝑣𝑣𝑧𝑧 . 𝜌𝜌. 𝑑𝑑𝑑𝑑 = 𝑄𝑄. 𝜌𝜌̅ = 𝜋𝜋. 𝑅𝑅 2 . 𝑢𝑢�. 𝜌𝜌̅
0
(restricción de flujo másico cte.)
𝑅𝑅 𝑟𝑟
𝑣𝑣𝑧𝑧 (𝑟𝑟) 𝑅𝑅 2 ∫ 𝑟𝑟 𝜇𝜇
𝑑𝑑𝑑𝑑
=
𝑢𝑢� 2 ∫𝑅𝑅 𝑟𝑟 𝜌𝜌 ∫𝑅𝑅 𝑟𝑟 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
0 𝜌𝜌̅ 𝑟𝑟 𝜇𝜇
Solución 𝑣𝑣𝑧𝑧 (r)
7
𝜇𝜇 𝑦𝑦 𝑣𝑣𝑧𝑧 varían con z también
Si 𝑣𝑣𝑟𝑟 es ↓↓ y no influye en 𝑣𝑣𝑧𝑧
1 𝜕𝜕 𝜕𝜕(𝑣𝑣𝑧𝑧 . 𝜌𝜌)
(𝑟𝑟. 𝑣𝑣𝑟𝑟 . 𝜌𝜌) + =0
𝑟𝑟 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
1 𝑟𝑟 𝜕𝜕𝑣𝑣𝑧𝑧
𝑣𝑣𝑟𝑟 = − � 𝑟𝑟. 𝜌𝜌. 𝑑𝑑𝑑𝑑
𝑟𝑟. 𝜌𝜌 0 𝜕𝜕𝜕𝜕
Líneas de corriente,
𝑟𝑟0 𝑟𝑟
� 𝑟𝑟. 𝑣𝑣𝑧𝑧 (𝑟𝑟, 0)𝜌𝜌. 𝑑𝑑𝑑𝑑 = � 𝑟𝑟 ′ . 𝑣𝑣𝑧𝑧 (𝑟𝑟 ′ , 𝑧𝑧)𝜌𝜌𝜌𝜌𝜌𝜌′
0 0
8
REACTORES TUBULARES EN FLUJO TURBULENTO
Flujo turbulento velocidad ≠ 0 en las 3 direcciones
𝑉𝑉 = 𝑣𝑣 +Ψ
Ψ, vel que fluctúa 𝑣𝑣, promedio en el tiempo
1 𝑡𝑡
𝑣𝑣 = lim � 𝑉𝑉 𝑑𝑑𝑑𝑑
𝑡𝑡→∞ 𝑡𝑡 0
Flujo turbulento en tubos largos y vacíos
𝑣𝑣𝑟𝑟 = 𝑣𝑣𝜃𝜃 = 0
𝑣𝑣𝑧𝑧 (𝑟𝑟) tiene perfil promedio con el tiempo que varía muy
rápidamente cerca de la pared del tubo “capa límite
fluidodinámica”
Régimen laminar
Turbulencia desarrollada
z
• Modelo de dispersión axial
𝑑𝑑𝑑𝑑
𝑢𝑢� = 𝑟𝑟𝐴𝐴 𝐹𝐹𝐹𝐹
𝑑𝑑𝑑𝑑
𝑑𝑑𝑑𝑑 𝑑𝑑 2 𝐶𝐶
𝑢𝑢� = 𝐷𝐷 + 𝑟𝑟𝐴𝐴
𝑑𝑑𝑑𝑑 𝑑𝑑𝑧𝑧 2
Zona laminar,
𝑢𝑢�2 𝑅𝑅2
𝐷𝐷 = 𝐷𝐷𝐴𝐴 +
48 𝐷𝐷𝐴𝐴
9
𝐿𝐿 �𝑅𝑅
𝑢𝑢
Reactores laminares si 𝑅𝑅 > 8 𝐷𝐷
𝐴𝐴
Condiciones de contorno
𝑢𝑢�. 𝐶𝐶𝑖𝑖𝑖𝑖 = 𝑢𝑢�. 𝐶𝐶 (0− ) 𝑧𝑧 = 0
𝑑𝑑𝑑𝑑
𝑢𝑢�. 𝐶𝐶𝑖𝑖𝑖𝑖 = 𝑢𝑢�. 𝐶𝐶 (0+ ) − 𝐷𝐷 𝑧𝑧 = 0
𝑑𝑑𝑑𝑑 0+
𝐶𝐶 (𝐿𝐿− ) = 𝐶𝐶 (𝐿𝐿+ ) = 𝐶𝐶𝑜𝑜𝑜𝑜𝑜𝑜 𝑧𝑧 = 𝐿𝐿
“Condiciones de
𝑑𝑑𝑑𝑑 Danckwerts”
=0 𝑧𝑧 = 𝐿𝐿
𝑑𝑑𝑑𝑑
D = 0 a la entrada y a la salida del reactor
Reactor isotermo y reacción de 1er orden, solución
𝑃𝑃𝑃𝑃 𝑧𝑧 + 𝑃𝑃𝑃𝑃 𝑧𝑧 − 𝑃𝑃𝑃𝑃 𝑧𝑧
𝐶𝐶 (𝑧𝑧) = 𝑒𝑒𝑒𝑒𝑒𝑒 � � �C1 exp � s � + C2 exp � s ��
2 𝐿𝐿 2 𝐿𝐿 2 𝐿𝐿
C1, C2 constantes
z=L
𝑃𝑃𝑃𝑃
𝐶𝐶𝑜𝑜𝑜𝑜𝑜𝑜 4 𝑠𝑠 exp � 2 �
=
𝐶𝐶𝑖𝑖𝑖𝑖 𝑠𝑠 𝑃𝑃𝑃𝑃 −𝑠𝑠 𝑃𝑃𝑃𝑃
(1 + 𝑠𝑠)2 exp � � − (1 − 𝑠𝑠)2 exp �
2 2 �
�.𝐿𝐿
𝑢𝑢
𝑃𝑃𝑃𝑃 = 𝐷𝐷
, nº de Peclet
4 𝑘𝑘 𝑡𝑡̅
𝑠𝑠 = �1 +
𝑃𝑃𝑃𝑃
10
• PFR
𝐶𝐶𝑜𝑜𝑜𝑜𝑜𝑜
lim = 𝑒𝑒 −𝑘𝑘.𝑡𝑡̅
𝑃𝑃𝑃𝑃→∞ 𝐶𝐶𝑖𝑖𝑖𝑖
• CSTR
𝐶𝐶𝑜𝑜𝑜𝑜𝑜𝑜 1
lim =
𝑃𝑃𝑃𝑃→0 𝐶𝐶𝑖𝑖𝑖𝑖 1 + 𝑘𝑘. 𝑡𝑡̅
• Validez para Pe > 8
𝑃𝑃𝑃𝑃
𝐶𝐶𝑜𝑜𝑜𝑜𝑜𝑜 4. 𝑠𝑠. 𝑒𝑒𝑒𝑒𝑒𝑒 � 2 (1 − 𝑠𝑠)�
=
𝐶𝐶𝑖𝑖𝑖𝑖 (1 + 𝑠𝑠)2
� 𝑅𝑅
𝑢𝑢
R. turbulento 0,8 ≤ ≤ 2,5
𝐷𝐷
𝐿𝐿 � 𝐿𝐿
𝑢𝑢 𝐿𝐿
R. laminar 0,8 ≤ ≤ 2,5 𝑅𝑅
𝑅𝑅 𝐷𝐷
Dispersión axial en condiciones no isotérmicas:
Perfil de conc. plano Perfil de Tª plano
T(r,z) = T(z)
Balance de materia
𝑑𝑑𝑑𝑑 𝑑𝑑 2 𝐶𝐶
𝑢𝑢� = 𝐷𝐷 + 𝑟𝑟𝐴𝐴
𝑑𝑑𝑑𝑑 𝑑𝑑𝑧𝑧 2
Balance de energía
𝑑𝑑𝑑𝑑 𝑑𝑑 2 𝑇𝑇
𝑢𝑢� = 𝐸𝐸 + 𝐽𝐽
𝑑𝑑𝑑𝑑 𝑑𝑑𝑧𝑧 2
11
E, coef. de dispersión axial de calor
J, suma de todas las fuentes (W+R)
Régimen turbulento desarrollado,
𝑢𝑢�. 𝐿𝐿 𝑢𝑢�. 𝐿𝐿
(𝑃𝑃𝑃𝑃) 𝑇𝑇ª = ≈ 𝑃𝑃𝑃𝑃 =
𝐸𝐸 𝐷𝐷
𝑑𝑑𝑑𝑑 𝑑𝑑 2 𝑇𝑇 2 𝑈𝑈 (𝑇𝑇 − 𝑇𝑇𝑒𝑒𝑒𝑒 ) ∆𝐻𝐻𝑅𝑅 . 𝑟𝑟
𝑢𝑢� = 𝐸𝐸 − −
𝑑𝑑𝑑𝑑 𝑑𝑑𝑧𝑧 2 𝜌𝜌. 𝑐𝑐𝑝𝑝 𝑅𝑅 𝜌𝜌. 𝑐𝑐𝑝𝑝
𝑑𝑑𝑑𝑑
𝑢𝑢�. 𝑇𝑇𝑖𝑖𝑖𝑖 = 𝑢𝑢�. 𝑇𝑇(0+ ) − 𝐸𝐸 𝑧𝑧 = 0
𝑑𝑑𝑑𝑑 0+
𝑑𝑑𝑑𝑑
=0 𝑧𝑧 = 𝐿𝐿
𝑑𝑑𝑑𝑑 L
Dispersión axial de materia despreciable
Re > 10.000
Dispersión axial de calor despreciable
PFR
12
