Sesion 2

FENOMENOS DE TRANSPORTE
Luis Carrasco Venegas

Unidad de Posgrado de la Facultad de Ingeniería Química, Universidad Nacional del Callao
Unidad de Investigación de la Facultad de Ingeniería Química, Facultad de Ingeniería Química, Universidad Nacional del Callao
Centro de Investigación en Bioenergía.
Correo: lacarrascov@unac.edu.pe
ORCID: https://orcid.org/0000-0002-7832-3366
Web: https://www.luiscarrascovenegas.com/
Docente Renacyt
Colab. Rocio Vasquez Edquen
Estudiante del X Ciclo de la carrera profesional de Ingenieria Química.
V TRANSPORTE MACROSCOPICO DE MATERIA
5.1. Introducción
El balance macroscópico de materia es una de las herramientas fundamentales de los ingenieros
químicos y sus áreas afines como la metalurgia, alimentos, ambiental, etc. Se aplica en los
procesos de transferencia de masa donde no se toma en cuenta los aspectos difusionales, es decir
la transferencia de masa por diferencia de concentración. En ausencia de reacciones químicas, se
conoce también como la ecuacion de continuidad en su forma integral. En la mayoría de los
procesos ordinarios, se aplica este principio de conservación de materia en su forma
macroscópica. Este tipo de balance de materia también se aplica a sistemas de mezclado, sistemas
reaccionantes tanto en régimen estacionario como en régimen no estacionario, tanto en
condiciones isotérmicas y no isotérmicas; en este último caso se requerirá también del balance
macroscópico de energía.
5.2 Balance macroscópico de sustancias puras
El balance macroscópico de materia es la ecuacion de conservación de materia, en cuya
deducción no se considera el termino difusional; es decir, la transferencia de masa
por diferencia de concentración; se
considera los términos de transporte de masa por flujo o transporte advectivo.
El flujo másico 𝑚 es la cantidad de masa que fluye a través de una sección transversal por unidad
de tiempo, tal como se observa en la Figura 4.1. En a) se observa el esquema genérico donde
ndA es el vector diferencial de área y la componente del vector velocidad debe ser paralela al
vector diferencial de área; en b) se observa el mismo caso aplicado a un sistema de flujo
considerando un flujo normal al área superficial y un flujo no normal respecto al punto de entrada
en la superficie de control
La ecuacion general de transferencia de masa para un sistema no reaccionante, está dada por:
velocidad  velocidad  velocidad 
 de entrada   de salida   acumulacion 
     
 de materia  −  de materia  =  de materia  (4.1)
     
 al volumen   del volumen  en volumen el 
 de control   de control   de control 
Para un elemento pequeño de área, dA = ndA (diferencial de área como vector es igual vector
direccional por el diferencial de área en su forma escalar) en la superficie de control , entonces
la velocidad neta de transferencia de masa en la dirección del flujo está dada por
( Vx dA) = ( V cos  dA) = ( V .ndA)
Donde la velocidad neta de transferencia de masa debida al flujo es:

A
( Vx dA) =  ( V cos  dA) =  ( V .ndA)
A A
(4.2)
La velocidad acumulación de materia viene dada por:

t 
 dv (4.3)
v
Reemplazando (4.2) y (4.3) en (4.1)

-  ( V .ndA) = t v  dv
A
(4.4)
Reordenando (4.4)

t 
 dv +  ( V .ndA) =0 (4.5)
v A
La ecuacion (4.5) puede ser obtenida a partir del Teorema de Reynolds, yes conocida como la
ecuacion de conservación de masa total o ecuacion de continuidad.
Ha de decirse que existe considerable relativa. flexibilidad en la selección de un volumen de

control cuando se resuelve un problema. Sin embargo, suele ser conveniente seleccionar un
volumen de control cuya superficie de control sea normal al flujo en todas las localizaciones en
donde cruza el flujo. Ya que de esta manera el producto punto 𝑉 ∙ 𝑛 es simplemente la magnitud
de la velocidad.
Ahondando un poco más sobre el volumen de control, si consideramos la figura xxx en la cual se
tiene un volumen y una superficie de control
En tanto que, la componente normal de la velocidad Vn estaría dada por:
Vn = V cos  = V .n (4.6)
Y el flujo másico diferencial dm podría reescribirse entonces como:
dm = Vn dA = V .n dA (4.7)
Sí la expresión anterior se integra sobre la superficie de control anterior se tendrá el flujo másico
neto
mneto = mentra − msale (4.8)
que interactúa con el volumen de control.
mneto =   (V .n)dA
SC
(4.9)
Por lo tanto el principio de conservación de masa podría reescribirse como:
mentra − msale = mVC o mVC + msale − mentra = 0 (4.10)
La ecuacion (4.10) puede ser escrita mediante:


  dv + SC  (V .n)dA = 0
t VC
(4.11)
Las ecuaciones (4.5) y (4.11) son las mismas ecuaciones con diferente notación
Sí la integral de superficie se divide por partes se tendrá que:
  (V .n)dA =   V
SC salida
n A− 
entrada
 Vn A (4.12)
Donde 𝐴 representa el área de una entrada o salida, y los símbolos de sumatoria son usados para
enfatizar que todas las entradas y salidas deben ser consideradas.
La ecuacion de continuidad o conservación de masa total, si bien tiene bastante uso; sin embargo,
normalmente se utiliza como una ecuacion auxiliar en los cálculos de ingeniería
Algunas aplicaciones.
Régimen estacionario.- Supongamos el esquema siguiente:
La ecuacion de conservación de materia para este caso seria:
  (V .n)dA +   (V .n)dA −   (V .n)dA = 0

S1 S2 S3
(4.13)
Lo que implica que
1  V1  S1 + 2  V2  S2 − 3  V3  S3 = 0 (4.14)
La ecuacion (4.14) es :
m1 + m2 − m3 = 0 (4.15)
Es decir, la suma de todas las masas de entrada al volumen de control es igual a la suma de todas
las masas de salida del volumen de control.
Este principio se puede ampliar a múltiples entradas y salidas de flujo masico.
Cuando el sistema es móvil, el fluido se desplaza a una velocidad diferente al sistema

La velocidad relativa = Velocidad de flujo – Velocidad del sistema: Vr = V f − Vs
  (V
SC
r . n)dA =   (V
SC
r − Vs ) . ndA (4.16)
A partir de la ecuacion (4.11), al aplicar el teorema de Gauss al termino de flujo neto se tiene:
  
  t + .( V ) dv = 0
VC
(4.17)
Luego, al retirar la integral se tiene

+ .( V ) (4.18)
t
La ecuacion (4.18) también es una forma de la ecuacion de continuidad expresada en su forma
diferencial. Para un fluido incompresible, se transforma en:
.V = 0 (4.20)
Problema 4.1 Conservación de materia en una tobera divergente-convergente.
Por una tobera divergente-convergente circula agua, tal como se observa en la figura P_4.1. El
area de la sección de entrada es de 4 cm2 y su velocidad de 5 m/s. La salida tiene una sección de
2.5 cm2. Determine la velocidad media a la salida.
Solucion
El esquema respectivo es el siguiente:
La ecuacion de conservación de masa total a partir de la ecuacion (4.11) viene dada por:

t VC
 dv +   (V .n)dA = 0 (1)
SC
Dado que en la tobera no se produce acumulación de materia, (1) se transforma en:
  (V .n)dA = 0
SC
(2)
La ecuacion (2) Puede ser escrita para la entrada y salida, es decir

SCS
 (V .n)dA −
SCE
  (V .n)dA = 0 (3)
La ecuacion (3) es equivalente a:

2  V2  S2 − 1  V1  S1 = 0 (4)
Dado que el agua es un fluido incompresible, se puede considerar que la densidad es la misma,
es decir: 1 = 2 , con lo cual, la ecuacion (4) se tanto, a partir de la ecuacion (4) se obtiene:
S1
 V2 = V1  (5)
S2
Reemplazando datos en (5)
La velocidad media a la salida de la tobera es 8m/s. También reemplazando en (4) se puede

obtener el flujo de masa.
m = 2  V2  S2 = 1  V1  S1 (6)
Este flujo de masa puede ser encontrada a partir de cualquiera de las igualdades
m = 2  V2  S2 = 1000kg / m3  8m / s  2.5 10−4 m2 = 2kg / s
Problema 4.2
Perfil de velocidad en una rendija
Se hace circular aire entre dos placas planas y lisas; el flujo en la entrada es uniforme y presenta
las siguientes características: P1=3.5 kgf/cm2 abs, T1=27 oC, y <V1>=1.2 m/s. Suponiendo que el
flujo es isotérmico y la presion de descarga P 2 es igual a 1.75 kgf/cm2 a la misma temperatura,
¿Cuál es la velocidad en la línea central (velocidad máxima)?. Considere que la ecuacion del perfil
de velocidad viene dada por:
  y 2 
Vx = Vx max 1 −   
  H  
solución
El esquema representativo del problema es el siguiente:
La ecuacion de conservacion de masa a partir de la ecuacion (4.11) es:


  dv + SC  (V .n)dA = 0
t VC
(1)
Dado que en el sistema no hay acumulación de materia:

  dv = 0
t VC
(2)
Entonces la ecuacion (1) se transforma en:
  (V .n)dA = 0
SC
(3)
La ecuacion (3) desarrollada para la condición de entrada y salida, del siguiente modo:
−  1 (V .n)dA +  2 (V .n)dA = 0 (4)

A1 A2
La primera integral puede ser simplificada teniendo en cuenta que la velocidad de entrada <V1>
es constante, mientras que la segunda ecuacion depende del perfil de velocidades.
− 1  V1  S1 + 2  (V .n)dA = 0 (5)
S2
Reemplazamos la ecuacion del perfil de velocidad

W H
− 1  V1  S1 + 2  Vx max
1 − ( y / H ) 2  dzdy = 0 (6)
0 −H
Integrando (6)
H
 y3 
− 1  V1  S1 + 2WVx max  y − =0
3H 2  − H
(7)

Reemplazando los limites
4
− 1  V1  S1 + 2WVx max  H = 0 (8)
3
Por otro lado, S1 = 2 HW , con lo cual
3 
Vx max =  1  V1  (9)
2 2
Es decir que la velocidad máxima, la misma que se da en la línea central depende de la velocidad
media de entrada y de la relación de densidades, que a su vez dependen de la presion y
temperatura.
Problema 4.3
Tasa de acumulación de materia dentro de un volumen de control

Suponiendo que la configuración mostrada en la Figura xxx sea bidimensional, calcule la rapidez
de variación de la masa dentro del sistema.
Solucion
La ecuacion de conservación de materia considerando un ingreso, una salida y la acumulación

viene dada por:

t v
 dv +   2 (V .n)dA −  1 (V .n)dA = 0 (1)
A2 A1
Reemplazando los perfiles de velocidad para las condiciones de entrada y salida mostradas en la
Figura xx, al despejar el termino de acumulación se tiene:

W h2 W h1
  dv = −  2   Vmax 1 − y / h2  dzdy + 1   Vmax

 1 − ( y / h1 ) 2  dzdy = 0 (2)
t v 0 0 0 0
Simplificando la ecuacion (2)

h
2 1 h
 
t v
max  max  1 − ( y / h1 )  dy = 0
 dv = −  2WV 1 − y / h2 dy + 1WV   2
 (3)
0 0
Integrando y simplificando se tiene:
 1 2

t v 2 3
 h1 = 0
 dv = − 2WVmax h2 + 1WVmax (4)
Problema
Descarga de un fluido que se encuentra contenido en un tanque cilíndrico
Un tanque cilíndrico de radio R y altura H inicialmente llena de agua, esta se descarga por un
orificio situado en la parte inferior cuyo radio es R 0. Encuentre una ecuacion que permita obtener
la altura instantánea del agua en el tanque.
Datos
H = 4 m , R = 4 m , R0 = 0.04 m , Cd = 0.98 , t0 = 0 , t f = 2000 s
Solución
El esquema respectivo es el siguiente
La ecuacion de conservación de materia macroscópica es la siguiente

t v
 dv +   2 (V .n)dA −  1 (V .n)dA = 0 (1)
A2 A1
Dado que no se tiene entrada de materia al tanque, la ecuacion (1) se reduce a:

t v
 dv +  2 (V .n)dA = 0 (2)
A2
Desarrollando cada una de las integrales
dv
 = −   V  S2 = − Q2 (3)
dt
Dado que el fluido es incompresible, la densidad del fluido puede ser simplificado, luego (3) se
transforma en:
dv
= −Q2 (4)
dt
El flujo volumétrico a la salida se obtiene mediante la ecucion de Torricelli
Q2 = Cd  V  A0 = Cd 2 gh   R02 (5)
Donde
Cd : Coeficiente de descarga, cuyo valor depende del diseño y es cercano a la unidad.
R0 : Es el radio del orificio

Los demás términos tienen la misma nomenclatura usada en este capitulo
Por otro lado el volumen puede ser expresado del siguiente modo: v =  R h
2
Derivando el volumen:
dv dh
=  R2 (6)
dt dt
Reemplazando (5) y (6) en (4)
dh
 R2 = −Cd 2 gh   R02 (7)
dt
0 t
dh R2

H h
= −Cd 2 g  02  dt
R 0
(8)
La integración resulta en:

2 1/2
R H  1
t = 2     (9)
 R0   2 g  Cd
La ecuacion (9) permite obtener el tiempo total de descarga del fluido
La altura instantánea se obtiene a partir de (7)

2
 Cd  R0 2 
h =  −    (2 g )1/2 t + H 1/2  (10)
 2  R  
El programa respectivo es el siguiente

# Descarga de un fluido desde un tanque cilíndrico a través de un orificio
d(h)/d(t) = (-Cd * (R0 / R) ^ 2 * (2 * g * h) ^ 0.5)
h(0) = 4
g = 9.81
Cd = 0.98
R=2
R0 = 0.04
t(0) = 0
t(f) = 2000
# Solucion analitica
H=4
t1 = 2 * (R / R0) ^ 2 * (H / (2 * g)) ^ 0.5 / Cd
# Altura instantanea (para comprobar)

h1 = (-0.5 * Cd * (R0 / R) ^ 2 * (2 * g) ^ 0.5 * t2 + H ^ 0.5) ^ 2
t2 = 0
Para la altura de fluido dada, se obtiene un tiempo de descarga de 2303.69 s. equivalente a 38.4
min.
Por otro lado para un tiempo de 2000 s., la altura instantánea expresada en forma grafica es la
siguiente.
Problema
Descarga de un tanque cilíndrico con ingreso continuo
Encuentre la altura instantánea y el tiempo de descarga de agua desde un tanque cilíndrico de

radio R. La altura inicial del fluido es H; en el fondo del tanque hay un orificio de radio R 0 , tal
como se muestra en la Figura P_5.xxx; además por la parte superior se alimenta un flujo Q1
Los datos son los siguientes
H = 5 m , R = 1 m , R0 = 0.02 m , Cd = 0.98 , Q1 = 0.001m3 / s , t0 = 0 , t f = 1000 s
Solución:
El desarrollo de este problema es en forma similar al anterior problema.
La ecuacion general de conservación de masa total viene dada por:

t v
 dv +   2 (V .n)dA −  1 (V .n)dA = 0 (1)
A2 A1
Desarrollando la ecuacion (1)
dv
 = Q1 − Q2 (2)
dt
Considerando la densidad constante y reemplazando la derivada del volumen respecto al tiempo
por la derivada de la altura respecto al tiempo, además, usando la ecuacion de Torricelli en el
punto de salida, se tiene:
dh
 R2 = Q1 − Cd 2 gh   R02 (3)
dt
Reordenando la ecuacion (3)
dh Q R2
= 1 2 − Cd 2 gh  02
dt  R R
# FLUJO NO ESTACIONARIO DE UN FLUIDO EN UN TANQUE

d(h1)/d(t) = If (h1 > 0) Then ((V0 / (pi * R ^ 2)) - Cd * (2 * g) ^ 0.5 * h1 ^ 0.5 * (R0 / R) ^ 2) Else (0)
V0 = 0.001
h1(0) = 5
pi = 3.1416
Cd = 0.98
R0 = 0.02
g = 9.81
R=1
t(0) = 0
t(f) = 1000
La altura instantánea para 1000 segundos se obtiene de manera gráfica y es la siguiente
El tiempo de descarga se obtiene de forma iterativa, asignando diferentes valores de tiempo final;
es asi como se obtuvo un tiempo de 3000 segundos, donde la altura del fluido tiene un valor muy
próximo a cero, tal como se muestra en la figura. De hecho no se llegara nunca a cero debido a
que hay un ingreso permanente
Problema
Descarga de agua desde un tanque cónico invertido
Obtenga la altura instantánea del agua que se encuentra en el interior de un tanque cónico
invertido de radio R y un orificio en la parte inferior R 0. El tanque tiene una altura H y el agua se
encuentra inicialmente lleno en el tanque, tal como se muestra en la figura
Los datos para utilizar son los siguientes:
H = 4 m , R = 2.5 m , R0 = 0.02 m , Cd = 0.98 , t = 2879 s
Solución
La ecuacion de balance de masa total viene dada por:

t v
 dv +   2 (V .n)dA −  1 (V .n)dA = 0 (1)
A2 A1
Dado que no hay ingreso de agua al tanque, se tiene

t v
 dv = −   2 (V .n)dA = 0 (2)
A2
A partir de (2) se obtiene:
dv
= −Q2 (3)
dt
El volumen del tanque conico depende de la altura y el radio
1 1
v = A h =   r 2  h (4)
3 3
A : Area de la sección transversal a una determinada altura h.
Ahora por semejanza de triángulos se debe poner el radio en función a la altura, esto es:
r R R
= o r= h (5)
h H H
Reemplazando (5) en (4)
2 2
1 R  1 R
v =   h   h =      h3 (6)
3 H  3 H 
Por tanto su derivada es:
2
dv R dh
=      h2 (7)
dt H  dt
Reemplazando en (3) y utilizando la ecuacion de Torricelli
dh R02 −1.5
= −Cd H 2 g  2 h
2
(8)
dt R
La ecuacion (8) es de variables separables
R02
 h dh = −Cd H 2g  R2 
3/2 2
dt (9)
La integral es la siguiente
2 5/2 R2
h = −Cd H 2 2 g  02 t + C (10)
5 R
2 5/2
La condición de frontera es la siguiente: t = 0 h = H , con lo cual C = H y luego
5
reemplazando en (10)
2 5/2 R2 2
h = −Cd H 2 2 g  02 t + H 5/2 (11)
5 R 5
A partir de la ecuacion (11) se obtiene
2/5
 5  R2 2  
h =    −Cd H 2 2 g  02 t + H 5/2   (12)
 2  R 5  
La ecuacion (12) permite calcular la altura en función del tiempo. Para obtener el tiempo de
descarga se puede obtener a partir de (12) o (9).
Utilizando (9)
0 t
R02
h dh = −Cd H 2 g  2  dt
3/2 2
(13)
H
R 0
Lo cual permite obtener

1/2
2 1 R2  H 
t =   2   (14)
5 Cd R0  2 g 

# Descarga de un tanque cónico invertido
d(h)/d(t) = -Cd * (R0 / R) ^ 2 * H ^ 2 * (2 * g) ^ 0.5 * h ^ (-1.5)
h(0) = 3.99
t(0) = 0
t(f) = 2879
R0 = 0.02
R = 2.5
Cd = 0.98
g = 9.81
H=4
# Tiempo de descarga
td = (2/5) * (1 / Cd) * (R / R0) ^ 2 * (H / (2 * g)) ^ 0.5
La solución numérica conduce a la siguiente figura:
De la ejecución del programa se obtiene el tiempo de descarga del agua igual a 2879 s.,
equivalente a 47.98 min.
Problema
Descarga de agua desde un tanque cónico
Obtenga la altura instantánea del agua que se encuentra en el interior de un tanque cónico de
radio R y un orificio en la parte inferior R0. El tanque tiene una altura H y el agua se encuentra
inicialmente lleno en el tanque, tal como se muestra en la figura.
Los datos para utilizar son los siguientes:
H = 4 m , R = 2.5 m , R0 = 0.02 m , Cd = 0.98 , t = 7678 s
Solución
El esquema utilizado es el siguiente

La ecuacion de balance de masa total viene dada por:

t v
 dv +   2 (V .n)dA −  1 (V .n)dA = 0 (1)
A2 A1
Dado que no hay ingreso de agua al tanque, se tiene

t v
 dv = −   2 (V .n)dA = 0 (2)
A2
A partir de (2) se obtiene:
dv
= −Q2 (3)
dt
El volumen del fluido que queda en el tronco de cono se obtiene mediante:
v = vcono mayor − vcono menor (4)
1 1 1
v =   R2  H −   r 2  ( H − h) =   ( R 2  H − r 2  (H − h)) (5)
3 3 3
De la relación de radios
r R−h R2
= o r2 = 2
( H − h) 2 (6)
h H H
Luego el volumen será
1 R2
v =   ( R 2  H − 2 ( H − h)2  ( H − h)) (7)
3 H
1  ( H − h) 3 
v =   R2 H −
H 2 
(8)
3 
Derivando la ecuacion (8)
dv 1 2  ( H − h) 2   dh 
=   R  −3 − 
H 2   dt 
(9)
dt 3 
El flujo volumétrico a la salida viene dada por la ecuacion de Torricelli
Q2 = Cd 2 gh   R02 (10)
1  ( H − h) 2   dh 
  R 2  −3   − dt  = −Cd 2 gh   R0
2
(11)
 
2
3  H
A partir de (11) se obtiene la derivada de la altura instantánea respecto al tiempo, es decir:

2
dh R  H2
= −  0  Cd 2 gh  (12)
dt  R ( H − h) 2
2
 R0 
Hacemos el siguiente cambio de variable A =   Cd 2 g  H , con lo cual (12) se
2
 R
transforma en:
dh h1/2
= − A (13)
dt ( H − h) 2
La ecuacion (13) es una ecuacion diferencial de variables separables
( H − h) 2
 h1/2 dh = − A dt (14)
Hacemos el siguiente cambio de variables: h = y dh = 2 ydy

2
Reemplazando en (14)
( H − y 2 )2
 y 2 ydy = − A dt (15)
Desarrollando la ecuacion (15)
2 ( H 2 − 2 Hy 2 + y 4 )dy = − A dt (16)
Integrando
4 2
2H 2 y − Hy 3 + y 5 = − A  t + C (17)
3 5
Retornando a la variable original
4 2
2H 2 h1/2 − Hh3/2 + h5/2 = − A  t + C (18)
3 5
La condición inicial es: t = 0 h = H ; luego
4 2 16
C = 2H 5/2 − H 5/2 + H 5/2 = H 5/2 (19)
3 5 15
4 2 16
2H 2 h1/2 − Hh3/2 + h5/2 = − A  t + H 5/2 (20)
3 5 15
Para hallar el tiempo de descarga hacemos h = 0
2 1/2
16 H 5/2 16  R   H  1
t= 2
=      (21)
15  R0  15  R0   2 g  Cd
  Cd 2 g  H
2
 R
El programa es el siguiente
# Descarga de un tanque cónico
d(h)/d(t) = -A * h ^ (0.5) / (H - h) ^ 2
A = (R0/R) ^ 2 * Cd * (2 * g) ^ 0.5 * H ^ 2
H=4
h(0) = 3.9999
t(0) = 0
t(f) = 7600
R0 = 0.02
R = 2.5
Cd = 0.98
g = 9.81
# tiempo de descarga
td = 16 / 15 * (H / (2 * g)) ^ (0.5) * (R / R0) ^ 2 / Cd
El tiempo de descarga obtenido es: 7678.973 equivalente a: 127.98 s. Asimismo la altura

instantánea es la siguiente
Problema
Relación de tiempos de descarga de dos tanques cónicos

Suponiendo que los recipientes cónicos llenos de agua cuyas dimensiones son exactamente
iguales, se abren simultáneamente la llave de descarga. Cuál de los recipientes se descargara
primero?
solución
El tiempo de descarga del cono invertido es:

2 1/2
2 1  R  H 
tCI =       (1)
5 Cd  R0   2 g 
El tiempo de descarga del tanque cónico normal es
2 1/2
16 1 R H 
tCN ==      (2)
15 Cd  R0   2 g 
tCN 2/5 3
= = (3)
tCI 16 /15 8
Es decir que el tanque del cono invertido se descarga más rápido respecto al tanque cónico
normal.
La relación de tiempos de descarga puede ser observado en la Figura xxxx
Problema
Descarga de un tanque cilíndrico horizontal
Consideremos un cilindro de radio R y longitud L lleno de fluido se encuentra posicionado

horizontalmente. En su parte inferior (zona de contacto con el piso de apoyo), se encuentra un
orificio de descarga de radio R0. Encuentre la variación instantánea de la altura de fluido en
función al tiempo para las condiciones iniciales siguientes:
t = 0 h1 (0) = 2 R , h2 (0) = 1.8R , h3 (0) = 1.5R , h4 (0) = R , h5 (0) = 0.5R
Considere el caso especifico :
R = 1 m L = 2 m R0 = 0.02 m Cd = 0.98
Solución
El esquema grafico es el siguiente
Aplicamos la ecuacion de balance de masa total

t v
 dv = −   2 (V .n)dA = 0 (1)
A2
El desarrollo de la ecuacion (1) conduce a:
dv
= −Q2 (2)
dt
De la ecuacion de Torricelli
Q2 = Cd 2 gh   R02 (3)
Del grafico
dv = 2 x Ldh (4)
Utilizamos la ecuacion del circulo
( x − 0)2 + (h − R)2 = R 2 (5)
De la ecuacion (5) se obtiene:
x = 2hR − h 2 (6)
dv = 2 2hR − h 2 Ldh (7)
A partir de (7) se obtiene
dv dh
= 2 2hR − h2 L (8)
dt dt
dh C
2R − h = − d 2 g   R02 (9)
dt 2L
dh Cd 2 g   R02
=− (10)
dt 2L 2R − h
La ecuacion (10) de hecho tiene solución analítica y numérica, cuyo resultado final dependerá
de las condiciones iniciales del problema, es decir en el tiempo cero, cual es la altura del fluido?
(esta demás explicar que tanto el resultado analítico como el numérico deben dar los mismos
valores)
La solución numérica es inmediata; se escribe directamente la ecuacion diferencial, se da el valor

inicial y se obtiene la altura instantánea.
Para la solución analítica, se debe integrar la ecuacion (10) del siguiente modo
Cd
 2R − hdh = −
2L
2 g   R02  dt (11)
Cd
Hacemos el siguiente cambio de variable: A = − 2 g   R02 , luego la ecuacion (11) resulta
2L
en:
 2R − hdh = A  t (12)
Para integrar (12) utilizamos el siguiente cambio de variable
u = 2R − h u 2 = 2R − h dh = −2udu
Luego la ecuacion (12) se transforma en:
A
− u 2 du = t (13)
2
A partir de (13)
u3 A
− = t +C (14)
3 2
Volviendo a la variable original
2(2 R − h)3/2
− = At + C (15)
3
Reemplazando el valor de la constante A
4(2 R − h)3/2 Cd
= 2 g   R02t + C (10)
3 L
En esta etapa las condiciones de frontera determinan el valor de la constante, donde el valor de
h debe fluctuar entre 0  h  2R
Si por ejemplo, el cilindro se encuentra completamente lleno, entonces h = 2R , por lo tanto, la
constante será igual a cero.
Si por ejemplo, h = 1.8R cuando t=0 entonces, reemplazamos en (10)
4(2 R − 1.8R)3/2 4(0.2 R)3/2

C= =
3 3
Del mismo modo si h = 1.5R en t=0 entonces:
4(2 R − 1.5R)3/2 4(0.5R)3/2

C= =
3 3
Del mismo modo: h = R en t=0 , entonces
4(2 R − R)3/2 4 R3/2

C= =
3 3
Y se puede proseguir para cualquier valor de h comprendido en el intervalo mencionado
El tiempo de descarga en cada caso se obtiene despejando el valor de t de la ecuacion (10) una
vez conocido el valor de C y haciendo h=0
El programa desarrollado es el siguiente

# Descarga de un tanque cilíndrico horizontal
# solución numérica y analítica
d(h1)/d(t) = If (h1 > 0) Then (-Q1 / A1) Else (0)
Q1 = Cd * pi * R0 ^ 2 * (2 * g) ^ 0.5
A1 = 2 * (2 * R - h1) ^ 0.5 * L
h1_A = 2 * R - ((3 / (4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5 * t) ^ (2 / 3)
h1_AA = If (h1_A > 0) Then (2 * R - ((3 / (4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5 * t) ^ (2 / 3)) Else (0)
F1 = ((3 / (4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5) ^ (2 / 3)
td1 = (2 * R / F1) ^ (3 / 2)

Q2 = Cd * pi * R0 ^ 2 * (2 * g) ^ 0.5
A2 = 2 * (2 * R - h2) ^ 0.5 * L
h2_A = 2 * R - ((3 / (4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5 * t + (0.2 * R) ^ (3 / 2)) ^ (2 / 3)
h2_AA = If (h2_A > 0) Then (2 * R - ((3 / (4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5 * t + (0.2 * R) ^ (3 / 2)) ^ (2
/ 3)) Else (0)
B2 = (0.2 * R) ^ (3 / 2)
F2 = (3/(4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5
td2 = ((2 * R) ^ (3 / 2) - B2) / F2

Q3 = Cd * pi * R0 ^ 2 * (2 * g) ^ 0.5
A3 = 2 * (2 * R - h3) ^ 0.5 * L
h3_A = 2 * R - ((3 / (4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5 * t + (0.5 * R) ^ (3 / 2)) ^ (2 / 3)
/ 3)) Else (0)
B3 = (0.5 * R) ^ (3 / 2)
F3 = (3/(4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5
td3 = ((2 * R) ^ (3 / 2) - B3) / F3

Q4 = Cd * pi * R0 ^ 2 * (2 * g) ^ 0.5
A4 = 2 * (2 * R - h4) ^ 0.5 * L
h4_A = 2 * R - ((3 / (4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5 * t + (R) ^ (3 / 2)) ^ (2 / 3)
h4_AA = If (h4_A > 0) Then (2 * R - ((3 / (4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5 * t + (R) ^ (3 / 2)) ^ (2 / 3))
Else (0)
B4 = (R) ^ (3 / 2)
F4 = (3/(4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5
td4 = ((2 * R) ^ (3 / 2) - B4) / F4

Q5 = Cd * pi * R0 ^ 2 * (2 * g) ^ 0.5
A5 = 2 * (2 * R - h5) ^ 0.5 * L
h5_A = 2 * R - ((3 / (4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5 * t + (1.5 * R) ^ (3 / 2)) ^ (2 / 3)
/ 3)) Else (0)
B5 = (1.5 * R) ^ (3 / 2)
F5 = (3/(4 * L)) * pi * R0 ^ 2 * Cd * (2 * g) ^ 0.5
td5 = ((2 * R) ^ (3 / 2) - B5) / F5
L=2
Cd = 0.98
g = 9.81
pi = 3.1416
R0 = 0.02
R=1
h1(0) = 1.999
h2(0) = 1.8
h3(0) = 1.5
h4(0) = 1
h5(0) = 0.5
t(0) = 0
t(f) = 500
En el grafico se observa que para las tres primeras condiciones iniciales, el tanque aun no se ha
descargado, mientras que para los dos últimos casos la descarga ya concluyo.
Se verifica los siguientes tiempos de descarga
h= 2R m, td =1382.698 s
h=1.8R m, td = 1338.973 s
h=1.5R m, td = 1209.861 s.
h=R m, td = 893.8403 s.
h=0.5R, td = 484.6092 s
Problema
Flujo no estacionario en un tanque con recirculación
En el esquema que se muestra en la figura, obtenga el volumen instantáneo en el tanque,

considerando los siguientes datos:
Densidad del fluido : 1000 kg/m3
Flujo de masa a la entrada: 0.1 kg/s
Altura del tanque: 5 m
Altura inicial del fluido: 4 m.
Diámetro del tanque: 2m
Diámetro del orificio: 0.04 m.
Considere además un flujo de masa de reciclado
Solución
El balance de masa total o ecuacion de continuidad, esta dada por:

t v
 dv = −   (V .n)dA = 0 (1)
A2
La integración de la ecuacion (1) conduce a:
dv
= −Q2 (2)
dt
El volumen instantáneo viene dado por:
v =  D2h (3)

dv dh
=  D2 (4)
dt dt
De acuerdo al esquema, se requiere expresar las ecuaciones en términos del flujo de masa, por
tanto, los términos de (4) se requiere multiplicar por la densidad.
dv dh
 =  D2 (5)
dt dt
El termino de flujo neto de masa viene dado por:
  (V .n)dA = m − m
A2
1 2 + m2 − Q2 (6)
Igualando (5) y (6)
 dh
 D2 = m1 − m2 + m2 − Q2 (7)
4 dt
 dh 
 D2 = m1 −  d 2 2 gh (8)
4 dt 4
2
dh 4m1  d 
= −  2 g h1/2 (9)
dt  D 2  D 
Hacemos el siguiente cambio de variables
2
4m1 d
A= B=  2g
 D2 D
La ecuacion (9) se transforma en:
dh
= A − B h1/2 (10)
dt
dh
 A− B h 1/2
=  dt (11)
Hacemos: h = y 2 dh = 2 ydy y luego reemplazando en (11)
2ydy
 A − By = t + C (12)
1 ( A − u)
Hacemos: u = A − By , du = − Bdy , dy = − du , y =
B B
Reemplazando estos cambios, la ecuacion (12) se transforma en:
2 ( A − u)du
B2 
− = t +C (13)
u
La integración de (13) conduce a:
2 2
− 2
A ln u + 2 u = t + C (14)
B B
En términos de la variable y la ecuacion (14) se convierte en:
2 2
− 2
A ln( A − By) + 2 ( A − By) = t + C (15)
B B
Ahora en términos de la variable h
2 2
− 2
A ln( A − Bh1/2 ) + 2 ( A − Bh1/2 ) = t + C (16)
B B
Ahora aplicamos la condición inicial: t = 0 h = H
2 2
C=− 2
A ln( A − BH 1/2 ) + 2 ( A − BH 1/2 ) (17)
B B
2 2 2 2
− 2
A ln( A − Bh1/2 ) + 2 ( A − Bh1/2 ) = t − 2 A ln( A − BH 1/2 ) + 2 ( A − BH 1/2 ) (18)
B B B B
Simplificando (18)
2 2 2 2
t =− 2
A ln( A − Bh1/2 ) + 2 ( A − Bh1/2 ) + 2 A ln( A − BH 1/2 ) − 2 ( A − BH 1/2 )
B B B B
2 ( A − BH 1/2 ) 2 1/2 1/2

t = 2 A ln + (H − h )
B ( A − Bh1/2 ) B
# Flujo no estacionario en un tanque con recirculación
d(h)/d(t) = If (h < 0) Then (0) Else (A - B * h ^ (1 / 2))
h(0) = 5
A = 4 * m1 / (RHO * pi * D ^ 2)
B = (d/D) ^ 2 * (2 * g) ^ 0.5
m1 = 0.1
RHO = 1000
d = 0.04
D=2
pi = 3.1416
g = 9.81
td = 2 * (A / B ^ 2) * ln((A - B * H ^ (1 / 2)) / (A - B * h1 ^ (1 / 2))) + (2 / B) * (H ^ 0.5 - h1 ^ 0.5)
h1 = 0.001
H=5
t(0) = 0
t(f) = 3000
Resultados
El perfil de velocidad es el siguiente
El tiempo de descarga tanto por el método analítico como por el método numérico conducen a
t=2591.624 s.
Problema
Flujo no estacionario en un recipiente cónico
Un tanque cónico tal como se muestra en la figura se encuentra inicialmente lleno de agua cuya
altura es H. Por la parte superior tiene un flujo de entrada y por la parte inferior un flujo de salida
la cual esta dada por la ecuacion de Torricelli. Encuentre la altura instantánea del agua en el
tanque.
Solución
El esquema es el siguiente
Aplicamos la ecuacion de continuidad o ecuacion de balance de masa total


  dv + SC  (V .n)dA = 0
t VC
(1)
En base a los flujos en el recipiente cónico se tiene:
v
− Q1 + Q2 +  =0 (2)
t
Simplificando la ecuacion (2) y con la ayuda de la ecuacion de Torricelli
dv
= Q1 − Cd A0 2 g  h1/2 (3)
dt
El volumen instantáneo en el tanque cónico esta dado por:
1
v= A h (4)
3
El área instantánea de la sección transversal del cono es:
A =  r2 (5)
De la relación de triángulos que se forma entre el tanque cómico (real) y el cono que se forma
con el volumen de agua contenido en el cono (imaginario) se tiene:
r R
= (6)
h H
1 R2
v =  2  h3 (7)
3 H
dv R2 dh
=  2  h2 (8)
dt H dt
Igualando las ecuaciones (8) y (3)
R2 dh
 2  h2 = Q1 − Cd A0 2 g  h1/2 (9)
H dt
 Q  H 2 2

dh 1 2  R0 
=    − Cd H   2 g  h1/2 
1
(10)
dt h 2    R   R 
El programa desarrollado es el siguiente:

d(h)/d(t) = (1 / h ^ 2) * ((Q1 / pi) * (H / R) ^ 2 - Cd * H ^ 2 * (R0 / R) ^ 2 * (2 * g) ^ 0.5 * h ^ 0.5)
h(0) = 5
H=5
R=3
R0 = 0.02
Cd = 0.98
Q1 = 0.005
pi = 3.1416
g = 9.81
t(0) = 0
t(f) = 1000
El resultado grafico es el siguiente
Nota: Para los datos del problema, el nivel del agua en el tanque tiende a disminuir, dado que el
flujo volumétrico de entrada es menor al flujo volumétrico de salida.
Problema
Perdida de un gas ideal desde un tanque esférico
Un tanque esférico de 2 m de radio contiene gas propano a 200 atm de presion. El tanque tiene
una válvula de descarga de 2 cm de diámetro. Estando el gas a la presion indicada, la válvula sufre
un desperfecto, lo cual genera una fuga de gas a la atmosfera; encuentre la variación instantánea
de presion en el tanque.
Solución
Balance de masa.

t  d = −   ( V.n ) dA ...(1)
Integrando la ecuación (1)

( ) = − V  A ...(2)
t
Dado que el volumen del tanque permanece constante, la ecuación (1) se transforma en:
1 d V
=−  A ...(3)
 dt 
Por la ley de los gases ideales, se tiene:  = P  M/RT
Diferenciado se tiene:
M
d = .dP ...(4)
RT
Reemplazando las ecuaciones (3) y (4)
RT M dP P 
  =−  A ...(5)
P M RT dt 
: Volumen específico, : Volumen del recipiente
RT M dP P  m / M
  =−  A ...(5)
P  M RT dt 
m : Volumen molar del gas
Simplificando la ecuación (5)
1 dP P  RT A RT A
=−  =−  (6)
P dt PM  M 
Integrando la ecuación (6), en la cual se introduce el factor de corrección debido a la

transformación de atmosferas a N/m2
P t RT(101325) R20
 
dP
= −  dt (7)
P0 P 0 M 
La integración de la ecuación (7) se tiene:
 RT(101325) R20 
P = P0 exp  −  t  ...(8)
 M  
La ecuación (8) permite obtener la presión instantánea en el tanque y la ecuación (6) se integra
numéricamente y se obtiene los mismos valores.
En base a la expresión (8), se puede indicar que el vaciado de gas del tanque sigue una cinética
de primer orden.
El programa es el siguiente:
# Pérdida de gas de un tanque esférico
d(P)/d(t) = -B * P
B = (R * T * 101325 / M) ^ 0.5 * pi * R0 ^ 2 / vol
R = 0.082
T = 298
R0 = 0.02
vol = (4/3) * pi * RT ^ 3
M = 44
pi = 3.1416
RT = 2 # Radio del tanque
P(0) = 200
P1 = P0 * exp(-B * t)
P0 = 200
t(0) = 0
t(f) = 500
5.3 Balance de materia de mezclas no reaccionantes que contienen solutos
Las ecuaciones de conservación de materia también pueden ser aplicadas a sistemas que
contienen uno o mas solutos disueltos en un determinado solvente. Este tipo de aplicaciones
por lo general se realiza a mezclas de sales con agua en diversas proporciones, las cuales a su vez
pueden ser mezcladas con otras soluciones ya sea con la finalidad de incrementar o disminuir las
concentraciones de las soluciones para ser aplicados en un determinado proceso. En estos casos
se suele ignorar los cambios de volumen debido al mezclado pues en mezclas de sustancias
solubles entre si, 1 litro de agua mas 1 litro de etanol no produce 2 litros de mezcla; se puede
producir un incremento o una disminución del volumen dependiendo de la naturaleza de los
componentes; sin embargo, estos cambios son mínimos y para las aplicaciones dentro del campo
de la ingeniería quimica, resulta irrelevante.
Problema
Mezclado de un soluto en un tanque (continuo para la mezcla y semicontinuo para el soluto)
Suponer que se tiene un tanque de mezclado en el cual se tiene un volumen de soluto con una
concentración inicial de soluto A igual a 0 . En el instante t = 0, se deja que ingrese al tanque
•
el mismo soluto A con una concentración AE y un flujo, volumétrico . A la salida, la solución
tiene una concentración instantánea A , considerando que el flujo volumétrico de salida es el
•
mismo, es decir .
La concentración de A esta expresada en masa de soluto/ volumen de soluciones que no es igual
a la densidad, pues está es la masa de solución/volumen de solución.
Pueden ser utilizados los siguientes datos:

•
A0 = 4 g/L; AE = 2 g/L; v =500 L; . =0.8 L/min; t : 900 min
Solucion
Consideremos el esquema respectivo.
Figura 5.18 Esquema para el balance de materia en un tanque de mezclado sin

variación de volumen.
Se aplica el balance macroscópico de masa.

t  d = −   ( V .n ) dA (1)
El termino de acumulación de (1) puede ser escrita del siguiente modo:
 
t  d =
t
( A  ) (2)
En la ecuación (2) el volumen es constante, puesto que el flujo volumétrico de entrada y salida,
es el mismo, de tal forma que el termino de acumulación es independiente del volumen
 
( A  ) =  A (3)
t t
El término de flujo de la ecuación (1) está dado por

• •
 V . n dA =  −  AE  (4)
A • •
 = AE − A  (5)
t
Usan la siguiente definición.
•

S= (6)

•
: Flujo volumétrico  : Volumen del tanque.
S: Tiempo espacial (inversa del tiempo de residencia) que representa al número de volúmenes
que pasa en un determinado tiempo. Con estos cambios, la ecuación (5) se transforma en:
A
= AE S − A  S (7)
t
dA
 AE −A

= sdt (8)
Cuyo valor integrado resulta.
−Ln( AE − A ) = st + C (9)
Usamos la condición de frontera t = 0, A = AO ; con la cual la constante resulta.
C = −Ln ( AE − AO ) (10)
Reemplazando (10) en (9) y reordenado se obtiene.
A = AE − ( AE −AO ) e−St (11)
En el caso particular AE = 0, es decir cuando entra agua pura el tanque que se tiene.
A = AO e−st (12)
En el caso particular AO sea igual a cero, se obtiene
A = AE (1 − e−st ) (13)
El programa desarrollado describe las ecuaciones (11), (12) y (13)

d(Y)/d(t) = 0
Y(0) = 0
RHOA0 = 4 # g/L
RHOAE = 2 # g/L
VOL = 500 # L
FVOL = 1 # L/min
RHOA = RHOAE - (RHOAE - RHOA0) * exp(-S * t)
S = FVOL / VOL
# Si RHOAE=0
RHOA1 = RHOA0 * exp(-S * t)
# Si RHOA0=0
RHOA2 = RHOAE * (1 - exp(-S * t))
t(0) = 0
t(f) = 900 # minutos
Problema
Mezclado en un tanque con flujos de masa de entrada y salida variable (sistema semicontinuo
tanto para la masa total como para la masa de soluto)
Consideremos que en el sistema mostrado en la Figura xxx el flujo de masa y la concentración

varía con el tiempo, se debe establecer al balance de masa total y el balance de masa del soluto
respectivo. Al tanque que contiene inicialmente una masa m0 con una fracción de soluto x0, se

dejan que ingrese un flujo de masa m 1 con una fracción de soluto x1. A la salida se tiene un flujo

de masa m 2 con una fracción de soluto x 2 . Encuentre la ecuacion que permite el calculo de la
masa instantánea de soluto contenida en la solución.
Figura 5.19 Esquema para el balance de materia en términos de flujo de masa.
Datos
Masa inicial en el tanque: m0 = 100 kg


Flujo de masa a la entrada: m1 =10 kg/min

Flujo de masa a la salida del tanque: m 2 =5 kg/min
Fracción másica inicial de soluto: x0 =0.20
Fracción másica a la entrada del tanque: x1 =0.15
Tiempo de proceso: t =60 min
Solución
La ecuacion de transporte de masa total viene dada por:

t v
 dv = −   2 (V .n)dA (1)
A2
A partir de la ecuacion (1) se obtiene
dm  
= m1 − m2 (2)
dt

m= m ( 
1 − m2 t + m0 ) (3)
Ahora aplicamos el balance de masa de soluto

t v
 A dv = −   A (V .n)dA (4)
A2
El termino de flujo (transporte advectivo) de la ecuacion (4) viene dado por
  ( V.n ) d
A A

=m 
2 x2 − m1 x1 (5)
El termino de acumulación de la ecuacion (4) viene dado por
 d mdx2 dm
t   dV = dt (mx ) =
A 2
dt
+ x2
dt
(6)
Remplazando (5) y (6) en (4)
dx2 dm  
+ x2 + m2 x 2 − m 1 x1 = 0 (7)
dt dt
m
dx2 
+ x2 m (  )
1 − m2 + m2 x 2 − m1 x1 = 0 (8)
dt
dx 2 dt
= (9)

m 1 ( x1 − x 2 ) (m − m ) t + m

1

2 0
Integrando la ecuación (5.46)
 m (x  (m − m ) t + m
dx2 dt
= (10)
1 − x2 )
  
1 1 2 0
El resultado conduce a:
−
1 
Lnm 1 ( x1 − x2 ) =
1 (
 
Ln  m  ) 
1 − m2 t + m0  + C (11)
  
m1 m1 − m2
Usando la condición t = 0, x2 = x0, se obtiene la constante C, es decir.

m
C = Ln ( x1 − x 0 ) +
1
Lnm0 (12)
 
m1 − m2
( )
( )
 / m
−m  −m

  
1 1 2
m1 − m2 t + m0 
x 2 = x1 + ( x 0 − x 1 )  (13)
 m0 
La ecuación (13) permite calcular la concentración del soluto en cualquier instante.
El programa adjunto permite obtener la concentración másica instantánea en función al tiempo

d(X2)/d(t) = (-X2 * (m1p - m2p) - m2p * X2 + m1p * X1) / m
m = (m1p - m2p) * t + m0
X2(0) = 0.2
X2A = X1 + (X0 - X1) * (((m1p - m2p) * t + m0) / m0) ^ a

a = -m1p/(m1p - m2p)
X1 = 0.15
X0 = 0.2
m1p = 10 # kg/min
m2p = 8 # kg/min
m0 = 100 # kg
t(0) = 0
t(f) = 60 # min
Problema
Balance de materia en dos tanques de mezclado sin variación de volumen
Consideremos dos tanques de mezclado cuyos volúmenes iniciales de solución son V 1 y V2 cuya
concentración inicial es 0 en ambos tanques. Al tanque 1 se admite un flujo V1 de agua pura El
flujo mezclado V1 pasa al segundo tanque y de este se descarga fuera del tanque 2. Encuentre
la concentración instantánea en cada tanque.
Los datos son los siguientes

•
Flujo volumétrico:  =5L/min
Concentracion inicial: 0 = 0.2 kg/L
Volumen inicial en el tanque 1: v1 =200 L
Volumen inicial en el tanque 2: v2 =100 L
Figura 5.21 Esquema para el balance de materia de dos tanques de mezclado sin
variación de volumen.
Solucion
El balance de masa del soluto en el tanque 1, está dada por:

t   d +   ( V.n ) dA −   ( V. n ) dA = 0

A
S
A
E
A (1)
Como a la entrada solo hay agua para hacer el balance de masa del soluto, se tiene
 E
A ( V.n ) dA = 0
La ecuación (1) se transforma en:

t  dm +  

1
S
A ( V.n ) dA = 0 (2)
Dado que la masa se expresa como: m1 = 1 1 la ecuación (2) permite obtener.
d1 •
1 + 1 1 = 0 (3)
dt
Usar el siguiente cambio de variable.

1 1
= (4)
1 
Donde S1 = 1 / 1 , es el tiempo espacial. Con este cambio de variables la ecuación (3), se

transforma en:
d1
= −1 S1 (5)
dt
Integrando la ecuación (5) se tiene.
1 = 0e−S1t (6)
Por aplicación de la ecuación (1), el balance de masa en el tanque 2, conduce a:
d2 • •
2 = 1 1 − 2 1 (7)
dt
De la ecuación (7) se tiene:
• •
d2 1 1 1
= − 2 (8)
dt 2 2

• •
d2 1 1
= 0 e− S1 t − 2 (9)
dt 2 2
Expresando la ecuación (9) en términos de la velocidad espacial
d2
+ 2S2 = 0 S2e− S1 t (10)
dt
La ecuación (10) es una ecuación diferencial lineal, cuya forma es la siguiente.
dy ( )
+ P t .y(t) = Q(t) (11)
dt
Cuya solución conocida tiene la forma.
ye = e

P(t)dt P(t)dt
Q(t)dt + C (12)
Siguiendo el procedimiento dado por la ecuación (12), la solución de la ecuación (100), está
dada por:
2e = e

S2dt S2dt
.S20e−Setdt + C (13)
S20 ( S2 −S1 )t
eeS2t = e +C (14)
S2 − S1
Usamos la condición inicial siguiente:
Si t = 0, 2 = 0
Reemplazando en la ecuación (14)
0 S2 S
C = 0 − =− 0 1 (15)
S2 − S1 S2 − S1
Despejando (15) en (14)
0 S2 ( S2 −S1 )t 0 S1
2eS2t = e − (16)
S2 − S1 S2 − S1
La simplificación de la ecuación (16), permite obtener.
0  −S1t
S2e − S1e 
− S2t
2 = (17)
S2 − S1
Esta ecuación permite calcular la concentración de soluciones a la salida de dos tanques de
tamaño diferente, pero de flujo constante.
Ahora supongamos que ambos volúmenes 1 y 2 son iguales y dado que el flujo volumétrico
•
1 , es constante, entonces S1 = S2, luego, la ecuación (17) se transforma en:
0
2 = lim S2e− S1 t − S1e− S2t  (18)
S2 →S1 S2 − S1
La aplicación del límite conduce a:
2 = 0e− S1t + S1 0 t e−S1t (19)
La ecuación (19) permite calcular la concentración a la salida del segundo tanque con volúmenes
y flujo volumétrico iguales en ambos tanques.
El programa desarrollado es el siguiente
d(RHOA)/d(t) = -RHOA * S1
RHOA(0) = 0.2
d(RHOB)/d(t) = S2 * (RHOA - RHOB)
RHOB(0) = 0.2
RHO1 = RHO0 * exp(-S1 * t)

RHO2 = (RHO0/(S2 - S1)) * (S2 * exp(-S1 * t) - S1 * exp(-S2 * t))
S1 = FV1 / V1
S2 = FV1 / V2
FV1 = 5 # L/min
V1 = 200 # L
V2 = 100 # L
RHO0 = 0.2
# Cuando el tiempo espacial es el mismo en ambos tanques
RHO2A = RHO0 * exp(-S1A * t) + S1A * RHO0 * t * exp(-S1A * t)
S1A = FV1 / V1
t(0) = 0
t(f) = 120 # min
Problema
Balance de masa de una mezcla en una configuración de tres tanques con flujo constante.
Realizar el balance de materia de soluto que fluye en una solución a través de tres tanques de
volúmenes diferentes, tal como se muestra en la figura.
Utilice los siguientes datos
E = 0.15kg/L, 10 = 0.4kg/L, 20 = 0.2kg/L, 30 = 0.1kg/L

• L
0 = 2.5 , 1 = 400 L , 2 = 100 L , 3 = 200 L
min
t = 150 min
Figura 5.21 Esquema del balance de masa en 3 tanques de volúmenes diferentes, con
flujo constante.
SOLUCIÓN
De acuerdo con el esquema mostrado en la figura, este sistema es estacionario respecto al flujo
de masa total, sin embargo, respecto al soluto, el balance de masa es un régimen no estacionario.
En este caso, en base al mismo esquema, hay un flujo de entrada de soluto, E .
El balance de masa en cada tanque está dado por:
Tanque 1:
d1 • •
1 = E 0 − 1 0 (1)
dt
Tanque 2:
d2 • •
2 = 1 0 − 2 0 (2)
dt
Tanque 3:
d3 • •
3 = 2 0 − 3 0 (3)
dt
Las condiciones iniciales son:
t = 0, 1 = 10 , 2 = 20 , 3 = 30
Reordenando las ecuaciones (1), (2) y (3)
d1
= S1 ( E − 1 ) (4)
dt
d2
= S2 ( 1 − 2 ) (5)
dt
d3
= S3 ( 2 − 3 ) (6)
dt
Las ecuaciones (4), (5) y (6) son resueltas en forma simultanea utilizando un algoritmo numérico;
sin embargo, también puede ser resuelta haciendo uso de un método analítico.
# BALANCE DE MASA DE SOLUTO EN TRES TANQUES CON FLUJO CONSTANTE
d(RHO1)/d(t) = S1 * (RHOE - RHO1)
d(RHO2)/d(t) = S2 * (RHO1 - RHO2)
d(RHO3)/d(t) = S2 * (RHO2 - RHO3)
RHO1(0) = 0.4 # kg/L

RHO2(0) = 0.2 # kg/L
RHO3(0) = 0.1 # kg/L
RHOE = 0.15 # kg/L
S1 = V0 / V1 # 1/min
S2 = V0 / V2 # 1/min
S3 = V0 / V3 # 1/min
V0 = 2.5 # L/min
V1 = 400 # L
V2 = 100 # L
V3 = 200 # L
t(0) = 0
t(f) = 150 # min
Problema
Balance de materia de una mezcla en tres tanques con flujo volumétrico variable
Realizar el balance de materia en un sistema de tres tanques de mezclado donde los flujos
volumétricos de entrada y salida son variables tal como se muestra en el esquema.
Datos:
10 = 1.2 g/L , 20 = 0.8 g/L , 30 = 0.6g/L, E = 0.5 g / L

• L • L L • L •
0 = 2 , 1 = 0.8 , 2 = 1.2 , 3 = 0.5
min min min min
10 = 500 L , 20 = 800 L , 30 = 1000 L
Figura 5.22 Esquema del balance de masa en tres tanques con flujo variable.
SOLUCIÓN
La ecuación de balance de masa en el tanque (1), está dada por:

t   d +   ( V.n ) dA = 0
1 1 1 (1)
 • •
( 1 1 ) = E 0 − 1 1 (2)
t
Siguiendo el mismo procedimiento se obtiene las siguientes ecuaciones diferenciales en los
tanques (2) y (3)
 • •
( 2 2 ) = 1 1 − 2 2 (3)
t
 • •
( 3 3 ) = 2 2 − 3 3 (4)
t
Dado que el volumen es variable en cada tanque, además de las concentraciones en masa,
derivando la ecuacion (2) se tiene:
dv1 d • •
1 + 1 1 = E 0 − 1 1 (5)
dt dt
Por otro lado por aplicación del balance de masa global
d1 • •
=  0 − 1( ) (6)
dt
Reemplazando (6) en (5) y simplificando
•
d1 0
= ( E − 1 ) (7)
dt 1
Por otro lado, al integrar la ecuación (6), se tiene
( •
1 = 0 − 1 t + 10
•
) (8)

•
d1 0 ( E − 1 )
= (9)
dt ( −  )t + 
•
0
•
1 10
De forma similar se obtiene las ecuaciones respectivas para los tanques (2) y (3)
•
d2 1 ( 1 − 2 )
= (10)
dt ( −  )t + 
•
1
•
2 20
•
d3 2 ( 2 − 3 )
= (11)
dt ( −  )t + 
•
2
•
3 30
Se resuelve simultáneamente las ecuaciones (9), (10), (11) con las indicaciones iniciales
1 (0) = 10 , 2 (0) = 20 y 3 (0) = 30
# Balance de masa de soluto en tres tanques con flujo variable

d(RHO1)/d(t) = V0 * (RHOE - RHO1) / ((V0 - V1) * t + V10)
d(RHO2)/d(t) = V2 * (RHO1 - RHO2) / ((V1 - V2) * t + V20)
d(RHO3)/d(t) = V3 * (RHO2 - RHO3) / ((V2 - V3) * t + V30)
RHO1(0) = 1.2 # g/L
RHO2(0) = 0.8 # g/L
RHO3(0) = 0.6 # g/L
RHOE = 0.5 # g/L
V10 = 500 # L
V20 = 800 # L
V30 = 1000 # L
V0 = 2 # L/s
V1 = 0.8 # L/s
V2 = 1.2 # L/s
V3 = 0.5 # L/s
t(0) = 0
t(f) = 1800
5.3 Balance de materia de mezclas reaccionantes
Una serie de reacciones químicas se llevan a cabo en reactores químicos discontinuos (batch),
semicontinuos y continuos (CSTR: Continuous stirred tank reactor), donde los aspectos de
transferencia de masa por difusión molecular no son tomados en cuenta, pues se considera que
la mezcla tiene agitación perfecta. Para las aplicaciones se utiliza el balance de masa global o
ecuacion de continuidad en su forma integral y el balance de materia por cada uno de los
componentes de la mezcla. Para el caso de los reactores batch y semicontinuos, se obtienen
sistemas de ecuaciones diferenciales ordinarias y para el caso de los reactores CSTR se obtiene
un sistema de ecuaciones algebraicas lineales o no lineales os cuales dependen del modelo
cinetico.
El balance de materia aplicado a un determinado componente (por lo general al reactivo

limitante) es el siguiente
Tasa de   Flujo de   Flujo de  Tasa de 

 acumulacion  = entrada  +  salida  +  generacion  (5.21)
       
 de materia   de materia   de materia   de materia 
El primero y segundo termino de (5.21) representa el flujo neto de materia. Por tanto, la
expresión matemática puede ser escrita del siguiente modo.
d  n

 i
dt 
C d = −   Ci ( v  n ) dA +  ri d
i =1 
(5.22)
Reactores y sistema de reactores de mezcla completa
Este tipo de reactores se utiliza fundamentalmente para obtener grandes volúmenes de

producción, pues se elimina las etapas de carga y descarga de una operación discontinua. Como
consecuencia de la aplicación de los balances de materia se obtienen sistemas algebraicos, los
cuales son resueltos utilizando un software apropiado.
Problema
Consideremos un reactor químico continuo (CSTR), donde se lleva a cabo una reacción quimica
de primer orden. Encuentre la concentración del reactante a la salida del reactor, asi como su
conversión. Repita el problema considerando una cinética de segundo orden irreversible. La
reacción es del tipo: nA ⎯⎯
k
→ B (n=1: cinética de primer orden; n=2: cinética de segundo
orden)
Datos:
•
C AE = 2 mol/L; k =0.5 min-1,  =0.5 L/min;  =100 L.
Solucion
El balance de materia para el componente A viene dado por:
   
 C A d +  C A ( v  n ) dA =  rA d (1)
t   
Dado que un CSTR opera en régimen estacionario, el termino de acumulación se omite,
obteniéndose la siguiente ecuacion.
 
 C A ( v  n ) dA =  rA d (2)
 
La integración de (2) para una cinética de primer orden para un volumen constante conduce a:
• •
C A  − C AE  = rA (3)
Para una cinética de primer orden, la velocidad de reacción viene dada por: rA = −kC A y
reemplazando en (3) se tiene:
• •
C A  − C AE  = −kC A (4)
Despejando se tiene:
•
C AE  C AE
CA = = (5)
•
 + k 1 + k
Reemplazando datos
2
CA = = 0.0198 mol/L
1 + 0.5  (100 / 0.5)
Ahora la conversión puede ser obtenida a partir de:
C A = C AE (1 − x) (6)
CA 0.0198
x = 1− = 1− = 0.9901
C AE 2
Ahora si consideramos una cinética de segundo orden, a partir de (3) se obtiene

• •
C A  − C AE  = −kC A2 (7)
Ahora esta ecuacion se puede resolver utilizando el procedimiento de una ecuacion de segundo
grado o utilizando un software. Utilizaremos el software.
f(CA) = CA * FV - CAE * FV + k * CA ^ 2 * V
FV = 0.5
k = 0.5
V = 100
CAE = 2
X = 1 - CA / CAE
CA(max) = 2
CA(min) = 0.01
Se obtiene los siguientes resultados
C A =0.1365; x =0.9317
Problema
Reactor tipo tanque continuo con una cinética reversible
Supongamos que en un reactor tipo tanque (CSTR) se lleva a cabo una reacción química
⎯⎯
k1
reversible del tipo A ⎯ → B . El reactor tiene un volumen constante  . La concentración a la
⎯
k 2
entrada es CAE y el flujo de entrada es  . Encuentre la concentración de A y B a la salida del CSTR.

El esquema se muestra en la figura xxx
Datos:
k1 = 2 min–1 k2 = 0.5 min–1

CAE = 2 mol /L  = 2 L / min  = 1000 L
Solución
A partir de la ecuación (5.22) al eliminar el termino de acumulación se obtiene
 
 C A ( v  n ) dA =  rA d (1)
 
• •
 C A ( v  n ) dA = CA − CAE  (2)
 r d = ( −k C
A 1 A + k2CB )  (3)

• •
CA − CAE  = ( −k1CA + k2CB )  (4)
Dado que la ecuación (4) contiene dos variables (CA y CB), este inconveniente se puede resolver
de dos modos:
a) Transformando CA y CB en términos de conversión.

b) Realizando el balance de materia de la especie B, con lo cual se consigue dos ecuaciones
con dos variables.
Usando el primer procedimiento, con la ayuda de la estequiometría
C A = C AE − C AE x (5)
CB = CAE x (6)
( CAE − CAE x ) − CAE = ( −k1 (CAE − CAE x ) + k2CAE x )  (7)
Simplificando la ecuación (6), se obtiene
k1
X= •
(8)
(k1 + k2 ) + / 
Luego se reemplaza en (5) y (6) y se obtiene las concentraciones respectivas de A y B
Usando el segundo procedimiento, haciendo el balance de B en régimen estacionario:

 
 CB

( v  n ) d=  r d
 B
(9)
 
• •
CB − CBE  = RB  (10)
De acuerdo con la estequiometria, RB = –RA. Resumiendo, se tiene:

• •
CA  = CAE − ( −k1CA + k2CB )  = 0 (11)
• •
CB − CBE − (k1CA − k2CB )  = 0 (12)
Se resuelve en forma simultánea las ecuaciones (11) y (12) y se obtiene las concentraciones de
A y B y a partir de estas se obtiene la conversión de A.
Utilizaremos un programa para resolver en forma simultanea ambos métodos y luego hacer las
comparaciones.

# Método 1
X1 = k1 / ((k1 + k2) + (FV / V))
k1 = 2.5
k2 = 0.5
CAE = 2
FV = 2
V = 1000
CA1 = CAE * (1 - X1)
CB1 = CAE * X1
# Método 2
f(CA2) = CA2 * FV - CAE * FV - rA2 * V
CA2(0) = 0.5
f(CB2) = CB2 * FV - rB2 * V
CB2(0) = 0.2
rA2 = -k1 * CA1 + k2 * CB2
rB2 = -rA2
X2 = 1 - CA2 / CAE
Ambos métodos reportan los mismos valores de concentración de A y B y de la conversión
C A = 0.3344; CB = 1.6655; x = 0.8327
Problema
Reacciones tipo tanque continuo en serie

En una batería de tres reactores en serie de diferentes volúmenes como se muestra en la figura
xxx se lleva a cabo una reacción del tipo A ⎯k⎯
→B . Encuentre la concentración del reactante a la
salida de cada uno de los reactores.
Datos
C AE = 1.5 mol / L , 1 = 500 L , 2 = 800 L , 3 = 1200 L , k = 0.002 min −1 ,

•
 = 0.5 L / min
Solución
El balance de masa esta dado por:
El balance de masa para los tres reactores está dado por:
 
 C A ( v  n ) dA +  rA d = 0 (1)
 
La forma integrada de la ecuación para cada reactor está dada por:
• •
Reactor 1: CAE  − CA1  − kCA11 = 0 (2)
• •
Reactor 2: CA1  − CA2  − kCA22 = 0 (3)
• •
Reactor 3: CA2  − CA3  − kCA33 = 0 (4)
De la ecuación (2) se obtiene el valor de CA1

•
CAE  CAE C
CA1 = •
= = AE (5)
kv1 1 + 1
+ kv1 1+ •

Donde 1, representa el tiempo de residencia del reactante A en el primer reactor.
A partir de la ecuación (3) se obtiene el valor de CA2, es decir:
•
C A1  C A1 C A1
C A2 = •
= = (6)
1 + 1• 2 1 + k2
kv
+ k1 v2

CAE
CA2 = (7)
(1 + k1 )(1 + k22 )
A partir de la ecuación (4) se despeja el valor de CA3
•
C A2  C A2 C A2
C A3 = •
= = (8)
1 + • 3 1 + k3
kv
+ k3v3

Reemplazando (7) en (8), se tiene
CAE
CA3 = (9)
(1 + k1 )(1 + k2 )(1 + k3 )
De esta manera se obtiene la concentración a la salida de los tres reactores.
En el caso límite, cuando los volúmenes y los flujos volumétricos son constantes e iguales
entonces  1 =  2 =  3 =  3 , se reemplaza en las ecuaciones respectivas y se obtiene:
CAE
CA1 = (10)
1 + k
C AE
C A2 = (11)
(1 + k )2
C AE
C A3 = (12)
(1 + k )3
Con los datos dados se resuelve simultáneamente las ecuaciones (2),(3), y (4) o mediante las
ecuaciones (10),(11) y (12). Se obtiene:
CA1 = 0.5 mol / L , CA2 = 0.1190476mol / L , CA3 = 0.017006 mol / L
También puede generalizarse para “n” reactores del mismo volumen y flujo volumétrico
constante, es decir:
C AE
C An = (13)
(1 + k )n
Se puede desarrollar un programa que resuelva directamente las ecuaciones (2), (3) y (4)
f(CA1) = CAE * FV - CA1 * FV + RA1 * V1
f(CA2) = CA1 * FV - CA2 * FV + RA2 * V2
f(CA3) = CA2 * FV - CA3 * FV + RA3 * V3
RA1 = -k * CA1
RA2 = -k * CA2
RA3 = -k * CA3
CA1(0) = 0.2
CA2(0) = 0.2
CA3(0) = 0.2
FV = 0.5 # Flujo volumétrico (L/min)

V1 = 500 # Volumen
V2 = 800 # Volumen
V3 = 1500 # Volumen
k = 0.002 # 1/min
CAE = 1.5
X1 = 1 - CA1 / CAE # conversión
# comprobación por el método analítico

C1 = CAE / (1 + k * TAU1)
C2 = C1 / (1 + k * TAU2)
C3 = C2 / (1 + k * TAU3)
TAU1 = V1 / FV
TAU2 = V2 / FV
TAU3 = V3 / FV
Al ejecutar el programa desarrollado se obtiene los mismos valores anteriores, es decir:

CA1 = 0.5 mol / L , CA2 = 0.1190476mol / L , CA3 = 0.017006 mol / L . También se
obtiene la conversión en cada tanque respecto a la concentración de entrada, x1 = 0.6667 ;
x2 = 0.9206 ; x3 = 0.9886
Problema
Batería de reactores continuos de volumen constante y cinética compleja

Se tiene una batería de tres reactores en serie tal como se muestra en la figura xxx. Encuentre
la concentración a la salida del reactor así como su conversión respecto de la concentración de
entrada
La cinética es del tipo:
⎯⎯
→ 3R + 1.25
2A ⎯
⎯
k1
k2
Datos:
CAE = 2 mol /L 3 = 500L

•
 L  −1
 = 0.8 L / min k1 = 0.0002   min
 mol 
−3.2
 L 
1 = 200 L k2 = 0.00005   min−1
 mol 
2 = 1000 L
Solución
Dado que el sistema se encuentra en régimen estacionario, el balance de masa de cada
componente en el primer reactor 1 está dado por:
 
Balance de masa de A:  CA

( v  n ) dA +  r d = 0
 A
(1)
 
 
Balance de masa de B:  CR

( v  n ) dA +  r d = 0
 R
(2)
 
 
Balance de masa de C:  CS

( v  n ) dA +  r d = 0
 S
(3)
 
La velocidad de reacción puede ser expresada en cada reactor del siguiente modo:
rAi = −k1C2Ai + k23CRi C1.5

Si (4)
Donde i es el número de reactor
En base a la estequiometría se tiene la relación de velocidades de reacción.
1 1 1
− rAi = rRi = rS (5)
2 3 1.2 i
A partir de la ecuación (1) se tiene:
• •
CAE − CA1 + rA1 1 = 0 (6)
De la ecuación (2), se obtiene:

•
0 − CR1 + rR1 1 = 0 (7)
De la ecuación (3), se obtiene:

•
0 − CS1 + rS1 1 = 0 (8)
Se procede del mismo modo para el reactor 2 y se obtiene:

• •
Balance de materia de A: CA1 − CA2 + rA2 2 = 0 (9)
• •
Balance de materia de R: CR1 − CR2 + rR2 2 = 0 (10)
• •
Balance de materia de S: CS1 − CS2 + rS2 2 = 0 (11)
Del mismo modo, se aplica el balance de materia para el reactor 3

• •
Balance de masa de A: CA2 − CA3 + rA3 3 = 0 (12)
• •
Balance de masa de R: CR2 − CR3 + rR3 3 = 0 (13)
• •
Balance de masa de S: CS2 − CS3 + rS3 3 = 0 (14)
Se resuelve simultáneamente las ecuaciones del (6) al (14) lo cual permite obtener las
concentraciones de salida de cada uno de los reactores.
Asimismo, la conversión en cada etapa del reactor se obtiene mediante:
x1 = 1 − C A1 / C AE , x2 = 1 − C A2 / C AE , x3 = 1 − C A3 / C AE
El programa desarrollado es el siguiente:

# Batería de reactores continuos de volumen constante y cinética compleja
f(CA1) = CAE * FV - CA1 * FV + RA1 * V1
f(CR1) = 0 - CR1 * FV + RR1 * V1
f(CS1) = 0 - CS1 * FV + RS1 * V1
f(CA2) = CA1 * FV - CA2 * FV + RA2 * V2

f(CR2) = CR1 * FV - CR2 * FV + RR2 * V2
f(CS2) = CS1 * FV - CS2 * FV + RS2 * V2
f(CA3) = CA2 * FV - CA3 * FV + RA3 * V3

f(CR3) = CR2 * FV - CR3 * FV + RR3 * V3
f(CS3) = CS2 * FV - CS3 * FV + RS3 * V3
RA1 = -k1 * CA1 ^ 2 + k2 * CR1 ^ 3 * CS1 ^ 1.25

RA2 = -k1 * CA2 ^ 2 + k2 * CR2 ^ 3 * CS2 ^ 1.25
RA3 = -k1 * CA3 ^ 2 + k2 * CR3 ^ 3 * CS3 ^ 1.25
RR1 = (-3/2) * RA1

RR2 = (-3/2) * RA2
RR3 = (-3/2) * RA3
RS1 = -(1.2 / 2) * RA1

RS2 = -(1.2 / 2) * RA2
RS3 = -(1.2 / 2) * RA3
CAE = 2
FV = 0.8
V1 = 200
V2 = 1000
V3 = 500
k1 = 0.0002
k2 = 0.00005
CA1(0) = 1.2
CA2(0) = 1
CA3(0) = 0.8
CR1(0) = 0.4
CR2(0) = 0.4
CR3(0) = 0.3
CS1(0) = 0.2
CS2(0) = 0.1
CS3(0) = 0.1
X1 = 1 - CA1 / CAE
X2 = 1 - CA2 / CAE
X3 = 1 - CA3 / CAE
Se obtiene los siguientes resultados:
CA1 = 1.8321 CA2 = 1.3747 CA3 = 1.2115

CR1 = 0.2517 CR2 = 0.9378 CR3 = 1.1826
CS1 = 0.1006 CS2 = 0.3751 CS3 = 0.4730
Y las conversiones respectivas
x1 = 0.08391 , x2 = 0.3126 , x3 = 0.3942
Problema
Batería de tres reactores con recirculación
En una batería de tres reactores tipo tanque continuo se lleva a cabo una reacción del tipo
⎯⎯
→ B. Encuentre la concentración de A y B a la salida de los reactores, considerando el
2A ⎯
⎯
k1
k2
reciclo en cada uno de los reactores, tal como se muestra en la figura. Encuentre las
concentraciones de A y B a la salida de cada uno de los reactores
Datos:
1 k2 0.0004
2 k1 0.002
7 CAE 2.
8 V1 1000.
9 V2 800.
10 V3 500.
11 FV0 2.
12 FV4 1.5
13 FV5 0.5
14 FV6 1.
15 FV1 3.
16 FV2 3.5
17 FV3 3.5
Solución
Balance de materia en el reactor 1
Balance de materia de A
   
 CA d = − C A ( v  n ) dA +  rAd (1)

t   
 
 
En régimen estacionario:  C A d = 0
t 

El término de flujo se descompone en:

   
 CA

( v  n ) dA =  CA ( v  n ) dA −  CAE ( v  n ) dA −  CA ( v  n ) dA
1 3
(2)
   
Reemplazando la ecuación (2) en (1) se tiene:

• • •
CAE 0 + CA3 6 − CA1 1 + rA1 1 = 0 (3)
rA1 = −k1C2A1 + k2CB1
Siguiendo el mismo procedimiento, se continua el desarrollo para el componente B
Balance de materia de B
• •
CB3 6 − CB1 1 + rB1 1 = 0 (4)
1
rB1 = − rA1
2
• • •
CA1 1 + CA3 5 − CA2 2 + rA2 2 = 0 (5)
rA2 = −k1C2A2 + k2CB2
• • •
CB1 1 + CB3 5 − CB2 2 + rB2 2 = 0 (6)
1
rB2 = − rA2
2
• •
CA2 2 − CA3 3 + rA3 3 = 0 (7)
rA3 = −k1C2A3 + k2CB3
• •
CB2 2 − CB3 3 + rB3 3 = 0 (8)
1
rB3 = − rA3
2
Se resuelve simultáneamente las ecuaciones (3), (4), (5), (6), (7) y (8).
# Balance de materia de A y B en el reactor 1

f(CA1) = CAE * FV0 + CA3 * FV6 - CA1 * FV1 + RA1 * V1
f(CB1) = CB3 * FV6 - CB1 * FV1 + RB1 * V1
RA1 = -k1 * CA1 ^ 2 + k2 * CB1
RB1 = -(1 / 2) * RA1
CA1(0) = 0.1
CB1(0) = 0.2
f(CA2) = CA1 * FV1 + CA3 * FV5 - CA2 * FV2 + RA2 * V2
f(CB2) = CB1 * FV1 + CB3 * FV5 - CB2 * FV2 + RB2 * V2
RA2 = -k1 * CA2 ^ 2 + k2 * CB2
RB2 = -(1 / 2) * RA2
CA2(0) = 0.1
CB2(0) = 0.2
f(CA3) = CA2 * FV2 - CA3 * FV3 + RA3 * V3
f(CB3) = CB2 * FV2 - CB3 * FV3 + RB3 * V3
RA3 = -k1 * CA3 ^ 2 - k2 * CB3
RB3 = -(1 / 2) * RA3
CA3(0) = 0.1
CB3(0) = 0.2
# Datos
CAE = 2
V1 = 1000
V2 = 800
V3 = 500
FV0 = 2
FV4 = 1.5
FV5 = 0.5
FV6 = FV4 - FV5
FV1 = FV0 + FV6
FV2 = FV1 + FV5
FV3 = FV2
k1 = 0.002 # (L/mol)/min
k2 = 0.0004 # (1/min)
X1 = 1 - CA1 / CAE
X2 = 1 - CA2 / CAE
X3 = 1 - CA3 / CAE
Resultados
1 CA1 0.9706722627
2 CB1 0.5146638687
3 CA2 0.7304764321
4 CB2 0.634761784
5 CA3 0.5905605236
6 CB3 0.7047197382
20 X1 0.5146638687
21 X2 0.634761784
22 X3 0.7047197382
Reactores discontinuos (batch) y reactores semicontinuos
Los reactores batch se utilizan principalmente para aplicaciones a baja escala para la manufactura
de productos dentro del área de quimica fina, donde se carga el reactante o reactantes en el reactor
y se espera que se todos los productos de reacción queden dentro del reactor para su posterior
separación por diversos medios dependiendo de los productos formados; si alguno de los
productos es un gas y son evacuados por algún punto del reactor, este deja de ser un batch y se
denomina reactor semicontinuo.
Hay una serie de configuraciones y modos de operaciones de los reactores semi continuos, dentro
de los cuales podemos considerar que los flujos de entrada y salida son diferentes pero constantes
ingresos de un solo reactivo al reactor, mientras el otro se encuentra cargado; los flujos de entrada
y salida variable, múltiples reactores con flujo de entradas y salidas constantes o variables
operando en condiciones isotérmicas, adiabáticas o con transferencia de calor al medio, etc.
Cabe mencionar que un reactor batch, es un caso especial de un reactor semicontinuo y debido a
ello podemos hacer el tratamiento en forma indistinta.
Problema
Batería de reactores semi continuos.

En una batería de tres reactores semi continuos en serie se lleva a cabo la reacción del tipo
⎯⎯
→ B. Encuentre la concentración y volumen instantáneo en cada uno de los reactores
2A ⎯
⎯
k1
k1
mostrados.
Solución
Dado que el sistema se encuentra en régimen no estacionario, se aplica el balance de materia
en cada uno de los reactores.
   
 Cid +  Ci ( v  n ) dA =  rd (1 )

t    i
 
La ecuación (1), puede ser escrita del siguiente modo:
   
 Cid = − Ci ( v  n ) dA +  rd (2)
t    i
  
Que al integrar permite obtener
d ( Ci i ) • •
= CEi 0 − Ci 0 + rii (3)
dt
Aplicando la ecuación (3) para la conservación de A en el reactor 1
d( CA1 1 ) • •
= CAE 0 − CA1 1 + rA1 1 (4)
dt
La ley de velocidad se expresa mediante: rA1 = −k1C2A1 + k2CB1
A partir de la ecuación (4) se tiene
dCA1 d1 • •
1 + CA1 = CAE 0 − CA1 1 + rA1 1 (5)
dt dt
Por otro lado del balance de masa total, con la consideración que la densidad es constante, se
tiene:
d1 •
(
•
= 0 − 1 ) (6)
dt
( •
1 = 0 − 1 t + 10
•
) (7)

•
dCA1 0 ( CAE − CA1 )
= + rA1 (8)
dt 1
Usando la ecuación (3), se aplica el balance de materia de B
dnB1 (9)
= −nB1 + ng
dt
d •
dt
( CB1 1 ) = −CB1 1 + rB1 1 (10 )

•
dCB1 CB1 0
=− + rB1 (11 )
dt 1
1
Con rB1 = − rA1
2
Se prosigue del mismo modo para los demás reactores y se obtiene:
Para el componente A
dCA2
=
( CA 1
− CA2 ) •
1 + rA2 (12 )
dt 2
rA2 = −k1C2A2 + k2CB
( •
2 = 1 − 2 t + 20
•
) (14)
Para el componente B:
•
dCB 1 ( CB1 − CB2 )
= + rB2 (15)
dt 2
1
rB2 = − rA2
2
Se continua el procedimiento del mismo modo para el reactor 3
Para el componente A
•
dCA3 2 ( CA2 − CA3 )
= + rA3 (16)
dt 3
rA3 = −k1C2A3 + k2CA3 (17)
( • •
3 = 2 − 3 t + 30 ) (18)
Para el componente B
•
dCB3 2 ( CB2 − CB3 )
= + rB3 (19)
dt 3
1
rB3 = − rA3
2
Se resuelve simultáneamente las ecuaciones (8), (11), (12), (15), (16) y (19), con las condiciones
iniciales respectivas.
d(CA1)/d(t) = FV0 * (CAE - CA1) / V1 + RA1
d(CB1)/d(t) = FV0 * (CBE - CB1) / V1 + RB1
RA1 = -k1 * CA1 ^ 2 + k2 * CB1
RB1 = -(1 / 2) * RA1
CA1(0) = 0.2
CB1(0) = 0
d(CA2)/d(t) = FV1 * (CA1 - CA2) / V2 + RA2
d(CB2)/d(t) = FV1 * (CB1 - CB2) / V2 + RB2
RA2 = -k1 * CA2 ^ 2 + k2 * CB2
RB2 = -(1 / 2) * RA2
CA2(0) = 0.1
CB2(0) = 0
d(CA3)/d(t) = FV2 * (CA2 - CA3) / V3 + RA3
d(CB3)/d(t) = FV2 * (CB2 - CB3) / V3 + RB3
RA3 = -k1 * CA3 ^ 2 + k2 * CB3
RB3 = -(1 / 2) * RA3
CA3(0) = 0.05
CB3(0) = 0
# Datos generales
FV0 = 1.5 # L/min
FV1 = 1.8 # L/min
FV2 = 2.2 # L/min
FV3 = 1.2 # L/min
CAE = 2 # mol/L
CBE = 0
V10 = 800 # L
V20 = 600 # L
V30 = 500 # L
k1 = 0.15 # (L/mol) min^-1
k2 = 0.005 # min^-1
V1 = (FV0 - FV1) * t + V10
V2 = (FV1 - FV2) * t + V20
V3 = (FV2 - FV3) * t + V30
X1 = 1 - CA1 / CAE
X2 = 1 - CA2 / CAE
X3 = 1 - CA3 / CAE
t(0) = 0
t(f) = 180 # min
Resultados
Concentracion de A y B en el reactor 1
Concentracion de A y B en el segundo reactor

Concentracion de A y B en el tercer reactor
Volumen instantáneo en cada uno de los reactores

Perfil de conversión de A en cada uno de los reactores

Sesion 2

Cargado por

Copyright:

Formatos disponibles

Sesion 2

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Sesion 2

Cargado por

Copyright:

Formatos disponibles

FENOMENOS DE TRANSPORTE

Luis Carrasco Venegas

V TRANSPORTE MACROSCOPICO DE MATERIA

Donde la velocidad neta de transferencia de masa debida al flujo es:

La velocidad acumulación de materia viene dada por:

Reemplazando (4.2) y (4.3) en (4.1)

Ha de decirse que existe considerable relativa. flexibilidad en la selección de un volumen de

En tanto que, la componente normal de la velocidad Vn estaría dada por:

Y el flujo másico diferencial dm podría reescribirse entonces como:

mneto = mentra − msale (4.8)

que interactúa con el volumen de control.

Por lo tanto el principio de conservación de masa podría reescribirse como:

mentra − msale = mVC o mVC + msale − mentra = 0 (4.10)

La ecuacion (4.10) puede ser escrita mediante:

La ecuacion de conservación de materia para este caso seria:

  (V .n)dA +   (V .n)dA −   (V .n)dA = 0

Lo que implica que

Este principio se puede ampliar a múltiples entradas y salidas de flujo masico.

Cuando el sistema es móvil, el fluido se desplaza a una velocidad diferente al sistema

Luego, al retirar la integral se tiene

Problema 4.1 Conservación de materia en una tobera divergente-convergente.

El esquema respectivo es el siguiente:

Dado que en la tobera no se produce acumulación de materia, (1) se transforma en:

La ecuacion (2) Puede ser escrita para la entrada y salida, es decir

La ecuacion (3) es equivalente a:

Reemplazando datos en (5)

La velocidad media a la salida de la tobera es 8m/s. También reemplazando en (4) se puede

m = 2  V2  S2 = 1000kg / m3  8m / s  2.5 10−4 m2 = 2kg / s

Perfil de velocidad en una rendija

El esquema representativo del problema es el siguiente:

La ecuacion de conservacion de masa a partir de la ecuacion (4.11) es:

Dado que en el sistema no hay acumulación de materia:

Entonces la ecuacion (1) se transforma en:

−  1 (V .n)dA +  2 (V .n)dA = 0 (4)

Reemplazamos la ecuacion del perfil de velocidad

Tasa de acumulación de materia dentro de un volumen de control

La ecuacion de conservación de materia considerando un ingreso, una salida y la acumulación

  dv = −  2   Vmax 1 − y / h2  dzdy + 1   Vmax

Simplificando la ecuacion (2)

Integrando y simplificando se tiene:

H = 4 m , R = 4 m , R0 = 0.04 m , Cd = 0.98 , t0 = 0 , t f = 2000 s

La ecuacion de conservación de materia macroscópica es la siguiente

Dado que no se tiene entrada de materia al tanque, la ecuacion (1) se reduce a:

Desarrollando cada una de las integrales

Cd : Coeficiente de descarga, cuyo valor depende del diseño y es cercano a la unidad.

R0 : Es el radio del orificio

La integración resulta en:

La ecuacion (9) permite obtener el tiempo total de descarga del fluido

La altura instantánea se obtiene a partir de (7)

El programa respectivo es el siguiente

# Altura instantanea (para comprobar)

Descarga de un tanque cilíndrico con ingreso continuo

Encuentre la altura instantánea y el tiempo de descarga de agua desde un tanque cilíndrico de

Los datos son los siguientes

H = 5 m , R = 1 m , R0 = 0.02 m , Cd = 0.98 , Q1 = 0.001m3 / s , t0 = 0 , t f = 1000 s

El desarrollo de este problema es en forma similar al anterior problema.

La ecuacion general de conservación de masa total viene dada por:

Desarrollando la ecuacion (1)

# FLUJO NO ESTACIONARIO DE UN FLUIDO EN UN TANQUE