TP1 Fenomenos de Transporte Final

Departamento de Ingeniera Qumica Segundo
Cuatrimestre 2016
76.47
Fenmenos de
Transporte
Ingeniera
Qumica
Trabajo Prctico N1
Integrantes
Compte Alvaredo Juan Ignacio
Padrn
96216
Klenk Camila
95638
Louzau Denise
95796
Romero Juan Sebastian

Pacchy Matas
/09/2016
96042
76.47 Fenmenos de Transporte Ingeniera Qumica

Trabajo Prctico N1
Alvaredo Klenk - Louzau Pacchy - Romero
Ejercicio 1
a) En este problema se pide obtener una expresin de la variacin de la
presin en el lquido. Para ellos se establecen las siguientes
suposiciones:
El sistema se encuentra esttico y en equilibrio.
La base del tanque es un cuadrado de 5m x 5m.
Se establece el sistema de coordenadas como se indica a
continuacin
Al encontrarse el sistema esttico y en equilibrio se aplica la ecuacin de

hidrosttica
P= . g
La orientacin del sistema de coordenadas es tal que la gravedad tiene
una nica componente, vertical, hacia abajo.
Como el aire posee muy baja densidad, se puede inferir que la variacin
de la presin con la altura de la columna de aire es despreciable, por ende,
sobre la superficie del lquido la presin ser la atmosfrica. Esta misma
apreciacin no se puede realizar en el seno del lquido, donde la densidad
del agua no es desestimable, y a medida que se desciende en el mismo, la
altura de la columna de agua modifica la presin. A partir de estas
consideraciones se puede reescribir la ecuacin como:
dP
=g
dy
P
Patm
dP= g dy
0

Trabajo Prctico N1
P=Patm+ g y
P=101325 Pa+ 1000
kg
m
9,8 2 y
3
m
s
De este modo las isobaras en el tanque se disponen de la siguiente

manera:
Y=0m -> P=101325 Pa Pa
Y=4,5m -> P=145425 Pa

Y= 5m -> P=150325 Pa
b) Para calcular la fuerza de presin resultante en el plano inclinado, se integra

la expresin de la presin en funcin de la altura.
dF=P n dS= P n dz dL
El signo menos indica que la fuerza es contraria a la normal saliente de la
superficie. En el siguiente esquema se representa a la normal en el mismo
sentido que la fuerza y su descomposicin en sus correspondientes
coordenadas x e y.
n
L

Trabajo Prctico N1
Por trigonometra se descompone la altura del tringulo en L= y / (sen),
y la normal en
n =cos i +sin j
; de esta manera se llega a la siguiente
expresin:
cos i+ sin
dy
( Patm+ gy ) ( j )dz
sin
5
dF=
4,5
F=739375 N (cos 30 i+sin 30 j)

F=640317,5 N i +369687,5 N j
El punto de aplicacin se calcula igualando dos fuerzas, una que va

desde y=4,5 m, hasta el punto de aplicacin (L) y otro que va desde y=5m
hasta el mismo punto.
cos i+ sin
cos i+sin
dy
( Patm+ gy ) ( j )dz
sin
dy
( Patm+ gy ) ( j ) dz
sin
0
L
0
5
4.5
y
g4,5
g5
L
Patmy + g Patm4,5
=Patm5+
PatmL g
2
2
2
2
Los resultados de dicha ecuacin cuadrtica son:
y= -25,43 (este valor se descarta porque una altura no puede ser

negativa)
y=4,752 m (a 0,252 m del comienzo del plano inclinado)
Para calcular el correspondiente punto en x, se obtiene

trigonomtricamente como:

Trabajo Prctico N1
x=
0,252m
=0,436 m
tan30
5 m0,436 m=4,564 m
Entonces, el punto de aplicacin se encuentra sobre el plano inclinado en

x=4,564m e y=4,752m.
Ejercicio 2
a) Se elige como modelo de barrido el modelo de Mezclado Perfecto por
sobre el de Flujo-Pistn debido a la existencia de un agitador dentro del
tanque, lo cual nos permite suponer que el contenido del mismo se
encuentra perfectamente mezclado. Asimismo, resulta difcil pensar en la
presencia de fluidos inmiscibles dentro del recipiente, con lo cual es
descartado el modelo de Flujo-Pistn. Esta suposicin de mezclado perfecto
garantiza la homogeneidad de propiedades punto a punto, dentro del
volumen de control. Adems se supone que el fluido se mezcla
instantneamente ni bien ingresa al mismo.
Para poder calcular las concentraciones de las tres especies A, B y C en
funcin del tiempo es necesario separar el proceso en dos etapas. La etapa
inicial es aquella en la cual el tanque comienza a llenarse y la especie A
comienza inmediatamente a descomponerse de acuerdo a la reaccin
A B C
; la etapa final del problema es aquella en el cual el tanque ya
se encuentra lleno, comienza a salir fluido del mismo y las reacciones

qumicas continan su curso. Tambin se consideran necesarias para la
resolucin del problema las siguientes suposiciones:

Trabajo Prctico N1
Modelo de Mezclado perfecto: garantiza propiedades como la densidad ()

constantes en todo el volumen de control.
Solucin diluida: de esta forma las densidades parciales (

expresarse como
pueden
i=x i i
b)
Etapa inicial
Para poder calcular las concentraciones de A, B y C hasta el momento de
t llenado ) se toma como volumen de control el fluido ingresando al
llenado (
tanque.
Datos de los que se dispone:
me=2
kg
s
x entrada
=0,05
a
V =3 m
k 1=5
k 2=0,02
kg
m3 s
=1000
kg
m3
kg
3
m s
Se conocen adems las expresiones de las velocidades de las reacciones

de descomposicin en cadena
r 1=k 1 x a
r 2=k 2
A B C :
(1)
(2)
1. Se comienza planteando el correspondiente Balance de masa total:
dV = (v n) d
t Vc
Sc

Trabajo Prctico N1
De acuerdo con las suposiciones adoptadas para la resolucin del
ejercicio se plantea constante en todo el volumen de control, de manera
que el balance de masa adopta la forma:
dV (t) me
=
dt
dV (t )=
me
dt
(3)
Se integra la ecuacin (3) entre las condiciones iniciales y finales para

esta etapa, es decir
t inicial =0 ;
V (t
=0
=0 )
inicial
t final =t
V (t
final
=V (t)
=t )
respectivamente,
V (t )
dV (t)=
0
me
dt
Se obtiene la siguiente expresin para el volumen de fluido en el tanque

en funcin del tiempo:
V (t)=
me
(4)
Despejando el tiempo se obtiene:
t=V (t)
me
Y reemplazando en la ecuacin (5) con

obtener
(5)
V (t)=V (t
llenado
=3 m3
se puede
t llenado :
kg
3
m
=1500 s
kg
2
s
1000
3
t llenado =3 m
2. Se plantea a continuacin un balance de masa parcial para el componente

A:

Trabajo Prctico N1
x dV = x a ( v n ) ds r a dV
t Vc a
Sc
Vc
Reemplazando en la ecuacin anterior con la expresin (1) de consumo
de componente A se obtiene:
d
entrada
x a dV = x a
< v >S k 1 x a (t) dV
dt Vc
Vc
d ( xa ( t)V ( t) )
x entrada
me
k 1 x a (t)V (t)
a
=dt
(6)
Se integra la ecuacin (6) entre las condiciones iniciales y finales para

esta
etapa,
x ta V (t
=t )
x a (t )V (t )
inicial
inicial
=0
=0)
t final =t
respectivamente,
d ( x a (t) V (t) )
x entrada
me
k 1 x a (t)V (t)
a
x ta V (t
t inicial =0 ;
decir
=x a (t)V (t)
final
final
es
= dt
0
x entrada
me
k 1 x a (t)V (t) t
1
a
ln
=
k1
x entrada
me
a
Al despejar
x a (t )
de la ecuacin anterior se obtiene:
k
me
x entrada
a
x a (t)=
1e
k 1 V (t )
Reemplazando la expresin (4) hallada para la dependencia del volumen

con el tiempo:
k
x entrada
x a (t)= a
1e
k1 t
(7)

Trabajo Prctico N1
Del cual se es posible obtener
x a ( t llenado ) =0,0067

B:
x dV = x b ( v n ) ds+ r b dV
t Vc b
Sc
Vc
Reemplazando con las expresiones (1) y (2) de produccin y consumo de

B, respectivamente; y las expresiones (4) para V (t) y (7) para
x a (t) , se
obtiene:
d ( x b (t)V (t))=
entrada me
entrada k me
me
xa
xa
e
2 t k dt
2
1
Se integra la ecuacin entre las condiciones iniciales y finales para esta

etapa,
x tb V (t
es
=x b (t)V (t)
final
final
tinicial
t inicial =0 ;
decir
=t )
x b V (t
inicial
=0
=0)
t final =t
respectivamente y se obtiene:
2
entrada
t
k
k 2t
me xa
me
entrada
x b (t) V (t)=
xa
t
+
(e 1)
Reemplazando nuevamente con la expresin (4) para V(t),
x b ( t )= x
entrada
a
t
k
k2 t
x entrada
a
(1e
)
2
k1 t
De la cual es posible obtener
x b ( t llenado ) =0,02834 .

C:

Trabajo Prctico N1
x dV = x c ( v n ) ds + r c dV
t Vc c
Sc
Vc
d (x c (t)V (t))
=k 2 V (t )
dt
De este balance se obtiene la siguiente expresin para la fraccin msica

de C en funcin del tiempo:
x c (t)=k 2
t
2
De la cual es posible obtener
x c ( t llenado )=0,015 .
Tabla de Resultados Obtenidos en la Etapa inicial del proceso
t llenado
x a ( t llenado )
1500 s
0,0067
x b ( t llenado )
0,02834
x c ( t llenado )
0,015
c)
Etapa final
Para poder calcular las concentraciones de A, B y C a partir el momento
de llenado (
t llenado ) se toma como volumen de control el fluido que se
encuentra dentro del tanque.

1. Se plantea el correspondiente Balance de Masa General para la Etapa Final
del problema. En esta fase del proceso se considera que la masa que
ingresa al volumen de control es la misma que lo abandona, es decir que se
supone estado estacionario para la misma.
(v n)ds=0
Sc
ms=
me=
2. Se plantea un balance de masa parcial para el componente A, donde esta

vez el volumen no vara con el tiempo:
10

Trabajo Prctico N1
d ( x a (t) )
=m
x entrada
m
x a (t)k 1 x a (t )V
a
dt
d ( x a (t) )
x entrada
m
( k V + m)
a
x a (t ) 1
V
V
=dt
Integrando la expresin entre las condiciones iniciales (
x ta =x a (t llenado)
inicial
t final=t
) y finales (
x ta =x a (t )
final
t inicial =t llenado ;
) para esta etapa del
proceso se obtiene:
( m+ k V )
x entrada
x entrada
(t t
m
m
a
a
x a (t)=
+ x a ( t llenado )
e V
1V )
1V )
( m+k
( m+k
llenado
(8)
De la expresin (8) es posible obtener el valor de
x a ( 50 min=3000 s )=0,00588
3. Se plantea un balance de masa parcial para el componente B:
d ( x b (t))
=m
x b (t)+ k 1 x a (t )V k 2 V
dt
Reemplazando con la expresin (8) para
x a (t)
diferencial lineal no homognea de la forma
se obtiene una ecuacin
d xb
+ P ( t ) x b =f (t ) . Se halla
dt
una solucin de la misma mediante el mtodo del factor integrante. Dicho

mtodo consiste en multiplicar trmino a trmino la ecuacin por el factor
(t )=e
P ( t ) dt
=e
ms
dt
V
ms
=e V
t+C
11

Trabajo Prctico N1
m
m
(
t +C m
t +C
x b (t) Vm t+C V
d ( x b (t)) V
1
e
+
e
=e
( A ) +( B ) e
dt
V
xentrada
m
k1
a
A=
k 2
m+
k1V
Donde
m
+k 1V
( ttllenado )
V
x entrada
k1 m
k1
a
B=x a ( t llenado )
( m+
k1 V )
Reagrupando, es posible obtener:
( x (t ) e
b
m
t+C
V
)' =e
m
t+ C
V
(
1
( A ) +( B ) e
m+k
1V
(t t llenado )
V
Integrando la expresin entre las condiciones iniciales (
t inicial
b
t final
b
m
t inicial +C
V
m
t final +C
V
=x b (t llenado ) e
=x b (t) e
m
t llenado +C
V
m
t +C
V
finales
t final =t
) para esta etapa del problema y despejando
x b (t) , se cancelan las constantes de integracin C y se obtiene la

siguiente expresin:
(t
V
x b ( t )= A + e V
llenado
t )
[(
)(
B
1e
k1
k 1
( tt llenado)
)+ x ( t
b
llenado
De esta expresin es posible obtener el valor de
x b ( 50 min=3000 s )=0,01964
) A m
12

Trabajo Prctico N1
4. Se plantea por ltimo, el correspondiente balance de masa parcial para el
componente C:
d ( x c (t) )
=dt
x c (t)
k2 m
V
Integrando esta expresin entre las condiciones iniciales (
x tc =x c (t llenado )
inicial
) y finales (
t final =t
x tc =x c (t )
final
) y despejando
x c (t)
se obtiene:
m
x c ( t llenado ) ) V
(t
k 2 V ( k 2 V m
xc=
e
m
llenado
t )
De la cual puede obtenerse el valor de
x c ( 50 min=3000 s ) =0,0245
Tabla de Resultados Obtenidos en la Etapa final del proceso

t
3000 s
d)
x a ( 50 min=3000 s )
0,00588
x b ( 50 min=3000 s )
0,01964
x c ( 50 min=3000 s )
0,0245
13

Trabajo Prctico N1
Xa, Xb y Xc en funcion del tiempo

0.05
0.05
0.04
0.04
Fraccionesmsicas
0.03
Xb
0.03
Xa
0.02
Xc
0.02
0.01
0.01
0
0
500 1000 1500 2000 2500 3000 3500

Tiempo [s]
Para poder evaluar las expresiones de
x a ( t inicial =0 )
x b ( t inicial =0 )
es
preciso calcular el lmite con t tendiendo a cero. Es posible observar que a

medida que el tiempo avanza, las concentraciones de cada componente
tienden a los valores que se obtendran en el estado estacionario. Dichos
valores salen de plantear los siguientes balances:
Balance de Masa Parcial para el componente A en estado estacionario
para la etapa final del proceso:
x entrada
m
a
x a ( t )=
=0,005882
m+k
1V
Balance de Masa Parcial para el componente B en estado estacionario
k 1 x entrada
k V
a
x b ( t )=
2 =0,00125
m+
k1 V
m
Balance de Masa Parcial para el componente B en estado estacionario

14

Trabajo Prctico N1
x c ( t )=
k2 V
=0,03
Ma, Mb y Mc en funcin del tiempo

90
80
70
60
ma
50
Msas de cada componente [kg]
mb
40
mc
30
20
10
0
0
1000 2000 3000 4000

Tempo [s]
Es posible observar que a medida que el tiempo aumenta, los valores de

Ma, Mb y Mc tienden a los valores que resultan de multiplicar las fracciones
msicas del estado estacionario por V(t). Entre t=0 y
t llenado
las curvas
difieren de las encontradas para las fracciones msicas para este espacio de
tiempo pues, el volumen vara con el tiempo en esta etapa del proceso.
e)
Se muestra en los siguientes grficos las formas que adoptan las curvas
de fracciones msicas al hacer tender las constantes cinticas a infinito y a
cero respectivamente:
15

Trabajo Prctico N1
Fracciones msicas variando K1

0.06
Fracciones msicas
0.05
Xa K1
0.04
Xb K1
Xc K1
0.03
Xa K10,889
0.02
Xb K10,889
0.01
Xc K10,889
0
0
2000
4000
6000
Tiempo [s]
En este grfico puede observarse que al hacer tender a cero los valores
de K1, las concentraciones del componente B se vuelven negativas. Esto
puede explicarse ya que al mantenerse constante el valor de K2, se
consume ms B de lo que se produce. Adems, las concentraciones del
componente C no varan respecto a lo hallado para los valores originales de
las constates cinticas. Como la expresin hallada para
x c (t)
no depende
de K1, las curvas para los valores de Xc cuando la constante tiende a cero o
a infinito son los mismos y se muestran superpuestos en el grfico.

0.06
Fracciones Msicas
0.05
Xa K20,03
0.04
Xb K20,03
Xc K20,03
0.03
Xa K20
0.02
Xb K20
0.01
Xc K20
0
0
2000
4000
Tiempo [s]
6000
16

Trabajo Prctico N1
En este grfico puede observarse que al hacer tender a cero el valor de
K2, las concentraciones Xc y Xb tambin son nulas. Lo mismo ocurre para
Xa pero tanto en el caso en el que la constante cintica tiende a cero como
a infinito. Cuando K2 toma valores muy grandes la concentracin del
componente B se vuelve negativa, pues se consume ms B de lo que se
produce(recordar que el valor de K1 sse mantiene constante); por el
contrario la concentracin de C aumenta.
Dado que las fracciones msicas nunca pueden ser negativas se buscan
los valores lmites para K1 y K2 de manera que no existan incongruencias.
Los grficos hallados se muestran a continuacin:

0.06
Fracciones msicas
0.05
Xa K1
0.04
Xb K1
Xc K1
0.03
Xa K10,889
0.02
Xb K10,889
0.01
Xc K10,889
0
0
2000
4000
6000
Tiempo [s]
De acuerdo a lo expuesto anteriormente K1 no puede ser menor a 0,889

kg/m3s.
17

Trabajo Prctico N1

0.06
Fracciones Msicas
0.05
Xa K20,03
0.04
Xb K20,03
Xc K20,03
0.03
Xa K20
0.02
Xb K20
0.01
Xc K20
0
0 100020003000400050006000
Tiempo [s]
De acuerdo a lo expuesto anteriormente K2 no puede ser mayor a 0,03

kg/m3s.
Ejercicio 3
a) En la siguiente figura, se tiene un codo que se encuentra unido
a un tanque mediante una brida. Por el mencionado tanque, circula
aire.
Datos
Caudal msico de aire(m): 5 kg s-1

Temperatura ambiente: 25C
Dimetro de la caera (D): 0,45 m
Presin de entrada (Pent): 10 atm
Presin de salida (Psal): 1 atm
Altura del codo (R): 0,75 m
Densidad del aire (, a 25C): 1,1769 kg m -3
Tensin de rotura admisible de la caera (adm): 300 MPa
Aceleracin gravitatoria (g): 10 m s-2
18

Trabajo Prctico N1
Presin atmosfrica (Patm): 1 atm

Dimetro de los bulones: 0.01m
Viscosidad del aire (, a 25C): 1.8464E-5 Pa s
En primer lugar, se tomara como volumen de control al sistema: fluidocodo. Una vez elegido este, se proceder a realizar un diagrama de cuerpo
libre del definido sistema.
presion
F
P
Donde
es la fuerza peso y
presion
F
es la fuerza ejercida por la
presin del fluido.
Suposiciones
Estado estacionario
Fluido con rgimen turbulento
Propiedades constantes a lo largo del tiempo
19

Trabajo Prctico N1
Peso del codo despreciable
Planteadas las suposiciones, se plantea el balance de masa
d dV
Vc
dt
= ( v . n ) dS
Sc
Por estado estacionario, el trmino de la izquierda es nulo. Por lo tanto, y

teniendo en cuenta que se tiene un rgimen turbulento, la igualdad se
reduce a:
v> 2 dS
v >2 dS=
S2
S1
En ambos trminos tenemos la definicin de caudal msico, entonces:
me =ms=m
Esto permite calcular la velocidad media a partir de la expresin:
S vm
Siendo S la seccin de la caera.
Sabiendo la velocidad media, se puede hallar el nmero de Reynolds:
D v>
=
Obtenindose el valor de Re=766588,05. Este mismo es mucho mayor
que Re=10E4, por lo tanto, es correcto suponer un rgimen turbulento.
A continuacin, se plantea el balance macroscpico
movimiento:
de cantidad de
20

Trabajo Prctico N1
d v dS
Vc
dt
= v ( v . n ) dS+ F ext
Sc
Como se supone estado estacionario, el trmino de la izquierda de anula.

El segundo trmino de la izquierda se refiere a la sumatoria de fuerzas
externas aplicadas al volumen de control elegido. En esto caso se
consideraran la fuerza peso, las fuerzas de presin que comprimen al
volumen de control, las fuerzas ejercidas por la presin atmosfrica y la
fuerza que ejerce el sistema sobre al brida, siendo esta ltima la de inters.
Entonces se puede reescribir el balance as:
Vc
Sc
Sc
Sc
+ g dV p
F
n dS + patm
n dS= v ( v . n ) dS
Despejando la fuerza sobre la brida:
Sc
Vc
Sc
Sc
sist /brida= v ( v . n ) dS g dV + pn dS patm

F
n dS
Se proceder a desarrollar cada una de las integrales por separado y
agrupar los trminos resultantes de acuerdo a la direccin correspondiente,
ya que se trata de un balance vectorial.
v >
|n|dS
v >
|n|dS
v >2S j
v >
A)
Sc
s2
s1
s2
v ( v . n ) dS= v( v . n ) dS+ v( v . n ) dS=
B)
g dV =
Vc
21

Trabajo Prctico N1
C)
Sc
s2
s1
pn dS= p|n|dS p|n|dS=PsalS i PentS j
D) Para calcular la ltima integral, se va

expresin:
considerar la siguiente
Patm n dS=0
S
2
S1
Entonces, teniendo en cuenta la superficie cerrada A=S+S 1 +S2
S1
S2
Patm n dS= Patm n dS+ Patm n dS + Patm n dS=0

La primera integral es la que se despeja, obteniendo:
Patm n dS=PatmS jPatmS i

S
Por ltimo, se suman las expresiones obtenidas en A, B, C y D, quedando

la fuerza de sostn expresada como:
v> 2S+ PsalSPatmS
2
v > S + gRSPentS+ PatmS
j
F sosten=
Dando como resultado:
22

Trabajo Prctico N1
sist /brida=133,63 N i145094,11 N j

F
Pero la fuerza calculada es la que ejerce el sistema sobre la brida.
Entonces, por principio de accin reaccin, la fuerza buscada es opuesta en
signo pero de igual modulo a la obtenida.
brida/ sist=133,63 N i +145094,11 N j

F
b) Como segundo punto de problema, se debe calcular la cantidad de
bulones requeridos para sostener el codo. Para ello, se utilizara la siguiente
formula:
adm=
2Fsosten
Ab N b
Donde el Ab es el rea de los bulones, Nb es el nmero de bulones y
adm
es tensin de rotura admisible de la caera. Dada esta, se despeja
el nmero de bulones, dando como resultado:
N b=12,32 13
c) Como ltimo punto, se pide calcular el momento que experimenta la
unin tanque-codo. Para ello, se elegir como volumen de control al sistema
fluido+ codo. Se planteara el balance de cantidad de movimiento angular,
respecto del centro de la base del codo:
d r ( v ) dV
Vc
dt
ext
= r ( v )( v . n ) dS+ M
Sc
Por estado estacionario, el trmino de la izquierda es cero. En el caso del

momento externo, solo se tendr en cuenta el ejercido por las fuerzas de
presin a la salida del codo, ya que tanto la fuerza peso, las fuerzas de
presin a la entrada y las fuerzas debido a la presin atmosfrica son
paralelas al vector posicin, por lo tanto, sus momentos respectivos sern
nulos. El balance se reduce a la siguiente expresin:
S2
S2
sist /brida= r ( v )( v . n ) dS+ r ( p n ) dS

M
Se analiza cada integral por separado:
23

Trabajo Prctico N1
E) Primero se procede a resolver el producto vectorial dentro de la
integral.
0,0,R< v >
r =( 0, R , 0 ) v =( v > ,0, 0 ) =
Entonces la integral se reduce a:
v > S k
2
r ( v )( v . n ) dS=R
S2
F)
r =( 0, R , 0 ) n =( 1,0,0 )=(0,0,R)
Por lo tanto:
r ( p n ) dS=RPsalS k
S2
Entonces, al sumar ambas expresiones obtenidas se tiene:
v > 2S + RPsalS
R
sist /brida=
M
Ejercicio 4
Suposiciones:
a)
La temperatura del vapor es constante en la caera.

Fluido incompresible. Por lo tanto, Cp = Cv = C.
La densidad y la capacidad calorfica son constantes.
Mezclado perfecto: la temperatura en el fluido es uniforme.
La energa interna es independiente de la posicin.
Las paredes del recipiente son adiabticas.
24

Trabajo Prctico N1
Para proceder a la resolucin del ejercicio y simplificacin de cuentas, se
separar el problema en dos procesos: uno con balance de masa no
estacionario y otro, estacionario.
La eleccin arbitraria del volumen de control se tomar de forma tal que
ste sea el fluido.
1er. Intervalo de tiempo:

Se realiza un balance macroscpico de masa sin reaccin qumica no
estacionario, resultando:
dV = v dS
t Vc
Sc
Luego,
M = V = mt=
V (t)=
m
t
Procediendo anlogamente con la expresin de balance de energa total,

se llega a:
d
^ Q
U =m
he+
dt tot
Para el trmino de energa interna y entalpia, necesito elegir una Tref. En

particular, para la simplificacin de los clculos voy a elegir Tref = Te, Te:
temperatura de entrada a mi Volumen de Control (VC). De esta forma,
href =he hehref =0
Luego,
U =Q
dt tot
Por defincin:
25

Trabajo Prctico N1
U tot = V C ( T Te )
Q=U
o A ( TsT ) ,Ts :temperatura del vapor en el IQ
Reemplazando:
d
( V (t) ( T Te ) )=U o A ( t )( TsT )
dt
Expresando al volumen de fluido y al rea de intercambio de calor como

funcin del tiempo, se obtiene:
t
dV m
V ( t )= m

=
dt
A (t)=
V
Ao , Vo : volumen de flui do cuando el tanque estlleno
Vo
A ( t ) =m
t Ao
Vo
Reemplazando en el balance V(t), A(t) y dV/dt, y aplicando la regla de la

derivada del producto,
C ( T Te )
d ( T Te )
dV
t Ao
+ C V
=U o m
( TsT )
dt
dt
Vo
C ( T Te ) m+
C V
d ( T Te )
t Ao
=U o m
( TsT )
dt
Vo
Para facilitar las cuentas, se proceder a adimencionalizar la ecuacin:
d (T Te )
t Ao
Uom
(TsT )
C ( T Te ) m
dt
Vo
+
=
C ( TsTe ) m
C ( TsTe ) m
C ( TsTe ) m
C V
26

Trabajo Prctico N1
Simplificando los trminos redundantes:
T Te V d T Te
Uo Ao t TsT
+
=
TsTe m
dt TsTe
C Vo TsTe
Recordando que del balance de masa resultaba:
t=
V
m
Luego,
T Te
d T Te
Uo Ao t TsT
+t
=
TsTe dt TsTe
C Vo TsTe
Para facilitar la integracin se va a recurrir a un cambio de variables:
T Te
, :temperatura adimensional
TsTe
Uo Ao t
, :tiempo adimensional
C Vo
Siendo que,
=
dt d
TsT
=1
TsTe
Reagrupando,
d
= ( 1 )
d
Y luego,
27

Trabajo Prctico N1
d
1
+ 1+ =1
d
( )
sta es una ecuacin diferencial de 1er. Orden, cuya posible solucin es:
=1
1k e
Para hallar el valor de la constante de integracin, debo analizar las

condiciones de contorno y reescribirlas con el cambio de variables
correspondiente:
A t=0,T =Te
A =0,=0
Planteando los lmites correspondientes:
1k e
lim =0=1lim
0
0
1k e L H
1=lim
lim k e =k
0
0
'
Luego,
1e
=1
Reescribiendo la ecuacin con las variables originales, obtenemos la

funcin buscada:
Uo Ao
T Te
1e CVo
=1
TsTe
Uo Ao
t
C Vo
28

Trabajo Prctico N1
Ya modelizada la funcin de la temperatura, podramos averiguar la
temperatura cuando el tanque se llen, entonces:
t llenado =
Vo
m
Que reemplazada en la ecuacin modelizadora de la temperatura como

funcin del tiempo, se simplifica a:
Uo Ao
T llenadoTe
1e C m
=1
TsTe
Uo Ao
Cm
Luego, con los datos del problema:
t llenado =5026,5 seg y T llenado =369,6 K=
2do. Intervalo de tiempo:
Realizando el balance de masa macroscpico para el volumen de control,

que ahora pasa a ser estacionario en masa:
m
e= m
s= m
Anlogamente, para el balance de energa total, resulta:
^
dM U
=m
( h^ eh^ ) + Q
dt
Como el sistema es estacionario en masa:
^
dU
=m
( h^ eh^ ) + Q
dt
29

Trabajo Prctico N1
Reemplazando con la expresin para la transferencia de calor, con la
definicin de entalpa y energa interna para el fluido, y la masa total como
funcin del volumen, resulta:
Vo C
dT
=m
C ( TeT ) +Uo Ao (TsT )
dt
Despejando e integrando:
m C+Uo Ao
T
mCTe+Uo Ao Ts
dT
Que resulta en:
T=
mC
Te+Uo Ao Ts
m
C Te+Uo Ao Ts
+
e
mC+Uo Ao
m C+Uo Ao
m C +Uo Ao
t
Vo C
Luego calculando t para T =Teb=373 K , se obtiene t=718,9 seg
Por ltimo , t eb =718,9+5026,5=5745,4 seg 96 min 1,6 hr
b)
Graficando la Temperatura con respecto al tiempo para dos caudales
diferentes, se obtiene:
30

Trabajo Prctico N1
T(K) vs. t(s)

400
350
300
250
200
150
100
50
0
1000
2000
3000
T(m=5)
4000
5000
6000
7000
T(m=50)
Grfico 4.1. Variacin de la temperatura con respecto al tiempo a dos caudales

diferentes: uno a 5kg/s (azul) y otro a 50 kg/s (gris).
Observando la ecuacin:
Uo Ao
t
CVo
T Te
1e
=1
TsTe
Uo Ao
t
C Vo
Se puede ver fcilmente que la temperatura no va a depender

directamente del caudal msico de entrada. Lo que va a variar va a ser: los
intervalos de tiempo con balance de masa estacionario y no estacionario,
que estn relacionados directamente en la evolucin de la temperatura con
el tiempo y las condiciones en las que empez el sistema estacionario en
masa; y la relacin del rea del fluido que entra en contacto con la caera
de vapor, es decir, con el intercambio de calor del sistema.
Por ejemplo: para un caudal de 0,5 kg/s, el tiempo que va a tardar en
llenarse el tanque va a ser de 50265 segundos. En este caso, el lquido
entrara en ebullicin cerca de los 5700 segundos, mucho antes de llenarse
el tanque, es decir, un volumen menor. El hecho de que el tiempo se
encuentre en el mismo orden que para el caudal del ejercicio (5kg/s) se
31

Trabajo Prctico N1
debe a que al llenarse ms lento, el rea de intercambio de calor es menor
y, por consiguiente, se intercambia menos calor en ese tiempo. Con un
caudal de 5kg/s, como se vio, se llega a ebullicin en 5745 s, pero un
volumen de fluido 10 veces mayor. Mientras que para un caudal de 50 kg/s,
el tanque se va a llenar en, aproximadamente, 503 segundos, pero el
sistema, una vez que es estacionario en masa y en energa, se puede ver
por el grfico, que nunca llega a ebullicin, sino que su cota para la
temperatura es de, aproximadamente, 318K.
c)
El caudal mximo que este sistema puede llevar a ebullicin va a ser el
caudal que en el lmite con infinito tienda a esa temperatura. Planteando
esto matemticamente, resulta:
mC
Te+Uo Ao Ts
m
C Te+Uo Ao Ts
+
e
mC+Uo Ao
m C+Uo Ao
[]=
lim T =lim
lim T =lim
t
mx =
m
mC
Te+Uo Ao Ts
mC +Uo Ao
=Teb, siendo que
<0
mC+Uo
Ao
Vo C
Uo Ao(TsTeb)
Kg
=6,11
C (TebTe)
s
Ejercicio 5
m C +Uo Ao
t
Vo C
32

Trabajo Prctico N1
Introduccin
El ejercicio n5 presenta una situacin que involucra un chorro de agua
que sale de la boca de una caera con una cierta velocidad, en rgimen
turbulento. La misma choca con un disco de determinada masa que slo
puede moverse en el eje vertical. Dicha situacin se representa en el
siguiente diagrama.
z
Figura N1: Diagrama del
problema n5.
v2
h
R1
v1
h: altura de la columna de
agua
V1: velocidad de salida
V2: velocidad despus del
choque
y R1: radio del disco de salida
del chorro
A continuacin el problema plantea tres situaciones distintas:

1 El chorro impacta con el disco y forma un ngulo de 90
respecto de la velocidad inicial, v 1 distribuyndose por toda la
superficie.
2 El chorro impacta y vuelve en direccin contraria a la que
vena.
3 El chorro impacta y forma un ngulo de 110 respecto de su
velocidad inicial.
El problema pide hallar la altura de la columna de agua en estos tres
casos y analizar sus diferencias.
Los datos aportados son los siguientes:
Masa del disco: 0.5 kg

Velocidad media de salida del chorro: 2.5 m/s
Radio de la superficie de salida del chorro: 2 in5.08 cm
Rgimen de flujo turbulento
A estos datos se sumaron ciertas suposiciones:
El agua se encuentra a temperatura ambiente: 25 C

La densidad del agua se aproxima a 1000 kg/m 3
El sistema se encuentra en estado estacionario
33

Trabajo Prctico N1
Se desprecia el rozamiento del aire a medida que el chorro
asciende
El agua forma un columna de agua perfecta
Flujo incompresible
El choque con el disco es completamente elstico
Desarrollo
Para la resolucin de este problema se plantearon dos volmenes de
control. El primero corresponde a la columna de agua desde la base hasta
antes de que el lquido impactara con el disco. El segundo es el diferencial
de volumen que resta de la columna para que el agua impacte contra el
disco. A continuacin se muestra un diagrama de dichos volmenes de
control.
Figura N2:
diagrama de
Disc
o
Volumen de control
N2
Volumen de control
N1
En el primer volumen de control se realiz un balance de energa

mecnica, paso que se realiz para las tres situaciones y cuya ecuacin
permanece invariable en todas las situaciones.
En el segundo volumen de control se realiz un balance de cantidad de
movimiento. Ser este balance el que contenga la informacin de cmo
interacciona el chorro con el disco, dato que depende de la situacin que
estemos tratando.
A continuacin se presenta la resolucin de la primera situacin.
Resolucin caso N1: ngulo de 90
Primero realizamos un balance de energa mecnica en el volumen de
control N1.
BEM)
34

Trabajo Prctico N1
d
v
v
p
+ gh dV =
+ gh+ ( v n^ ) dS+W eje +W friccion +W SM +W comp/ exp
dt
2
2
Al no haber sistemas de superficies mviles ni de ejes, ni tampoco

expansin y compresin y se desprecia el rozamiento con el aire nuestra
ecuacin se reduce a:
d
v
v
p
+ g h dV =
+ g h+ ( v n^ ) d S
dt
2
2
Por estado estacionario:
v2
p
+ g h+ ( v n^ ) dS=0
2
Las superficies involucradas en esta integral son la superficie de la boca

de salida del cao, que llamaremos S 1, y esa misma superficie elevada
cierta altura h, que llamaremos S 2. De esta manera nuestro balance
queda:
v1
p
v
p
+ g h1+ 1 ( v1 n^ ) d S 1+ 2 + g h2 + 1 ( v2 n^ ) d S2 )=0
2
La altura h1 es cero debido a nuestro sistema de referencia presentado

en la figura N1. En ambas superficies la presin es la misma, es decir el
fluido sale expulsado debido a una diferencia de presin pero una vez fuera
de la caera la presin es la atmosfrica. De esta forma los trminos de
presin se anulan.
Por otro lado las velocidades van en la misma direccin que las normales
salientes de cada superficie, pero en S 1 va en sentido contrario mientras que
en S2 va paralela a la normal.
Entonces nuestra ecuacin queda:
v 13
v 23
(
S +
+ g h v 2 S)=0
2
2
v1
v
S= 2 +g h v 2 S
2
2
35

Trabajo Prctico N1
De esta expresin despejamos h, la altura de la columna de agua.
v 13 v 23
= + g h v 2
2
2
h=
(1
)
1
( v 3v23 )
2 gv2 1
Esta ecuacin se utilizar en todos los casos planteados por el problema.

A continuacin se realiz un balance de cantidad de movimiento en el
volumen de control N2.
BCM)
d
v dV = v ( v n^ ) dS+ F ext
dt
Existen dos fuerzas externas en este volumen de control, la del peso del
fluido contenido en el mismo y la del peso del disco. La primera es
despreciable y la segunda es posible calcularla ya que conocemos la masa
del disco. A esto se el estado estacionario del sistema, entonces el balance
queda:
0= v ( v n^ ) dSMg k^
Donde M es la masa del disco y g la aceleracin de la gravedad.
En este caso las superficies de control sern discos de misma rea que la
del cao, distanciadas por un diferencial de altura. Es por esto que podemos
despreciar el peso del fluido ya que el lquido contenido en el volumen de
control es despreciable. Adems al ser un volumen de altura muy pequea
las velocidades en ambas superficies son iguales en mdulo (o en este caso
velocidad media en la superficie), pero no en sentido ni direccin.
En este caso:
v2= v 2 > k^
v3= v 2 > ^ , siendo
disco.
v3
la velocidad luego del choque con el
36

Trabajo Prctico N1
d S 2=r dr d , siendo
choque.
d S 3=R dz d , siendo
S2
el disco de radio R antes del
S 3 , la superficie lateral del diferencial
de cilindro que tomamos como volumen de control.
^ sin ^j
^=cos i+
Reemplazando:
2
v 2 > ( cos i^ + sin ^j ) R dz d
2 ^
v > S k
Mg k^
0=
2
^ sin ^j ) d
v 2 >2 hR ( cos i+
0
v2 > S k^
Mg k^
0=
2
Debido a la periodicidad de las funciones coseno y seno, la integral

anterior se anula y la expresin queda como sigue:
2 ^
v 2 > S kMg
k^
0=
De donde:
2
v 2 > =
Mg
S
Reemplazando:
v 2 >2 =
0.5 kg10 m/s

kg
1000 3 8.107. 103 m2
m
v 2 > 0.785 m/s
37

Trabajo Prctico N1
Reemplazando en (1):
m 3
m
2.5
0.785
s
s
h=
m
m
210 20.785
s
s
) (
h 0.964 m

Se realiz el balance de cantidad de movimiento como en el caso
anterior, pero en este caso como el chorro choca y se devuelve en sentido
contrario la superficie que presentan flujo es solo una sola, el disco de radio
R.
0= v ( v n^ ) dSMg k^
En este caso se sabe que:
v2= v 2 > k^
v3= v 2> k^ , siendo
v3
la velocidad despus del choque.
S2 y S3 coinciden pero sus normales son opuestas
Reemplazando en el balance de cantidad de movimiento y recodando

que el fluido esta en rgimen turbulento:
v 2 >2 S k^ + ( v 2 > k^ ) <v 2 >dS

0=
v 2 >2 S k^
^
v 2 >2 S k
Mg k^
0=
Operando queda:
38

Trabajo Prctico N1
^
v 2 >2 S kMg
k^
0=2
v 2 >2 =
Mg
2 S
v 2 > 0.555 m/s

Reemplazando en (1) (resultado de balance de energa mecnica para
todos los casos):
3
m
m
2.5
0.555
s
s
h=
m
m
210 20.555
s
s
) (
h 1.39 m
Es de importancia remarcar que en los dos casos anteriores el flujo de
agua, luego del choque, se daba solo por una superficie. En el primer caso el
flujo se daba por la superficie lateral del cilindro, en el segundo por la
misma superficie por donde entraba y en este caso ser una combinacin
de los dos. Es decir habr dos superficies por donde el chorro de agua fluir
luego del choque. El balance de cantidad de movimiento es:
0= v ( v n^ ) dSMg k^
En esta situacin se sabe:
v2= v 2 > k^
Podemos observar en la figura N3 que la velocidad despus

Figura N3:
del choque tendr una componente radial, cuya direccin depender
Diagrama de
del ngulo , y una componente en el eje vertical negativa.
ngulos entres
velocidades antes
y despus del
choque
39

Trabajo Prctico N1
Entonces,
vector:
descomponiendo
dicho
v3= v 2 >cos(20 ) ^ v 2 >sin 20 k^
^ sin ^j
^=cos i+
De
esta
manera
el
balance
queda:
v 2 > S k^ + v ( v n^ ) d Slat + v ( v n^ ) d S3

0=
2
v 2 > S k^ + A+ B

0=
2
Desarrollo A como la integral en la superficie lateral
A= v ( v n^ ) R d dz
n^ =^
v = v 2 >cos (20 ) ^ v 2 >sin 20 k^

v n^ = v 2> cos20
Reemplazando
integrando
en
^
v 2 >2 hR ( cos 20 2 cos i+cos
20 2 sin ^jsin 20 cos 20 k^ ) d
0
A=
2
v 2 > hRsin 20 cos 20 k^
A=2
z:
40

Trabajo Prctico N1
En esta expresin h es la altura del volumen de control N2, por lo
tanto:
lim A=0
h0
Desarrollo B como la integral en el disco de radio R

En esta integral se realiza un cambio de notacin para la densidad que
pasa a llamarse .
B= v ( v n^ ) d d
n^ =k^
v = v 2 >cos (20 ) ^ v 2 >sin 20 k^
v n^ = v 2> sin 20
Reemplazando en la expresin:
2
^
v 2 > sin 20 ( cos 20 cos i+cos
20 sin ^jsin20 k^ ) d
2
B=
2
2
^
v 2 > sin 20 (k)
R2
B= 2
2
Volviendo al balance original y sustituyendo, adems ahora volvemos a

la notacin original para la densidad:
v 2 >2 S sin 20 2 k^ Mg k^
v 2 > 2 S k^ +
0=
Acomodando trminos la expresin para la velocidad media a la altura
h es:
41

Trabajo Prctico N1
2
v 2 > =
Mg
S(1+sin2 20 )
Es importante aclarar que el ngulo dentro de la expresin es el ngulo

que forma la velocidad despus del choque con respecto al plano xy. A
esta variable la llamamos .
Reemplazando valores en la expresin anterior llegamos a que:
v 2 > 0.743 m/s

Sustituyendo en la expresin (1) hallamos que h es:
h 1.02m
Resultados y anlisis
A continuacin se muestran los resultados de los tres casos estudiados y
se proceder a analizar sus diferencias.
ngulo
ngulo
respecto de la []
velocidad inicial
[]
90
0
Velocidad a la
Altura
de
la
altura
mxima columna de agua
[m/s]
[m]
0.785
0.964
110
20
0.743
1.02
180
90
0.555
1.39
Se puede observar que a medida que el ngulo aumenta la velocidad

disminuye y la altura de la columna aumenta. Esto puede deberse a la
influencia que tiene el ngulo de salida en la transferencia de cantidad de
movimiento hacia el disco.
Si observamos la expresin (2) podemos ver que en
=90
el
denominador alcanza un mximo, que implica un mnimo en la velocidad a

la altura h. Podemos deducir, entonces, que la transferencia de cantidad
de movimiento alcanza un mximo cuando el chorro impacta contra el disco
y se devuelve en la direccin contraria.
Es importante aclarar un aspecto de la simetra del problema que
permite comprender por qu se da el mximo en este valor.
42

Trabajo Prctico N1
Para cada el chorro choca contra el disco y se devuelve con este
ngulo formando un tronco de cono en las situaciones menores a 90 y
mayores a 0. En el caso de 90 el tronco de cono se fusion con la columna
de agua y en el caso 0 el cono se abre tanto que se transforma en un
plano.

TP1 Fenomenos de Transporte Final

Cargado por

Copyright:

Formatos disponibles

TP1 Fenomenos de Transporte Final

Cargado por

Información del documento

Descripción 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

TP1 Fenomenos de Transporte Final

Cargado por

Copyright:

Formatos disponibles

Departamento de Ingeniera Qumica Segundo

Romero Juan Sebastian

76.47 Fenmenos de Transporte Ingeniera Qumica

Al encontrarse el sistema esttico y en equilibrio se aplica la ecuacin de

76.47 Fenmenos de Transporte Ingeniera Qumica

De este modo las isobaras en el tanque se disponen de la siguiente

Y=0m -> P=101325 Pa Pa

Y=4,5m -> P=145425 Pa

b) Para calcular la fuerza de presin resultante en el plano inclinado, se integra

76.47 Fenmenos de Transporte Ingeniera Qumica

; de esta manera se llega a la siguiente

F=739375 N (cos 30 i+sin 30 j)

El punto de aplicacin se calcula igualando dos fuerzas, una que va

y= -25,43 (este valor se descarta porque una altura no puede ser

Para calcular el correspondiente punto en x, se obtiene

76.47 Fenmenos de Transporte Ingeniera Qumica

Entonces, el punto de aplicacin se encuentra sobre el plano inclinado en

; la etapa final del problema es aquella en el cual el tanque ya

se encuentra lleno, comienza a salir fluido del mismo y las reacciones

76.47 Fenmenos de Transporte Ingeniera Qumica

Modelo de Mezclado perfecto: garantiza propiedades como la densidad ()

Solucin diluida: de esta forma las densidades parciales (

t llenado ) se toma como volumen de control el fluido ingresando al

Datos de los que se dispone:

Se conocen adems las expresiones de las velocidades de las reacciones

1. Se comienza planteando el correspondiente Balance de masa total:

76.47 Fenmenos de Transporte Ingeniera Qumica

Se integra la ecuacin (3) entre las condiciones iniciales y finales para

Se obtiene la siguiente expresin para el volumen de fluido en el tanque

Despejando el tiempo se obtiene:

Y reemplazando en la ecuacin (5) con

2. Se plantea a continuacin un balance de masa parcial para el componente

76.47 Fenmenos de Transporte Ingeniera Qumica

Se integra la ecuacin (6) entre las condiciones iniciales y finales para

de la ecuacin anterior se obtiene:

Reemplazando la expresin (4) hallada para la dependencia del volumen

76.47 Fenmenos de Transporte Ingeniera Qumica

Del cual se es posible obtener

3. Se plantea a continuacin un balance de masa parcial para el componente

Reemplazando con las expresiones (1) y (2) de produccin y consumo de

Se integra la ecuacin entre las condiciones iniciales y finales para esta

Reemplazando nuevamente con la expresin (4) para V(t),

De la cual es posible obtener

4. Se plantea a continuacin un balance de masa parcial para el componente

76.47 Fenmenos de Transporte Ingeniera Qumica

De este balance se obtiene la siguiente expresin para la fraccin msica

De la cual es posible obtener

Tabla de Resultados Obtenidos en la Etapa inicial del proceso

t llenado ) se toma como volumen de control el fluido que se

encuentra dentro del tanque.

2. Se plantea un balance de masa parcial para el componente A, donde esta

76.47 Fenmenos de Transporte Ingeniera Qumica

Integrando la expresin entre las condiciones iniciales (

) para esta etapa del

De la expresin (8) es posible obtener el valor de

3. Se plantea un balance de masa parcial para el componente B:

Reemplazando con la expresin (8) para

diferencial lineal no homognea de la forma

se obtiene una ecuacin