9.problemas Aplicados A La Ingenieria Quimica Pag 39

Tema 9: Problemas aplicados a la Ingeniera Qumica
TEMA 9.
PROBLEMAS APLICADOS A LA
INGENIERIA QUIMICA
Clculo numrico en Ingeniera
RegresionesRLM.1. Transmisinlinealesmltiplesdecalor. en una tubera

La transmisin de calor a travs de tuberas se sabe que el nmero de Nusselt
(Nu) funcin del nmero de Reynolds y del nmero de Prandlt:
Nu f (Re, Pr)
hD Dv C
f ,
k k
Una correlacin ampliamente utilizada es la ecuacin sugerida por Sieder y
Tate:
Nu 0.023Re 0.8 Pr1/ 3
Los coeficientes de la anterior ecuacin se han obtenido por correlacin de
muchos datos experimentales. Aunque en general la anterior ecuacin se
puede escribir como
Nu a Re b Prc
Donde a, b, c son parmetros ajustables.
Para los siguientes datos experimentales obtenidos en una tubera de de
pulgada de dimetro externo
a. compruebe la ecuacin original de Sieder y Tate, calculando el error
para cada punto y representndolo grficamente)
b. Utilice una correlacin lineal mltiple para hallar los parmetros a, b, c. Y
calcule y represente nuevamente los errores obtenidos en esta
correlacin.
Punto Re Pr Nu
1 49000 2.3 277
2 68600 2.28 348
3 84800 2.27 421
4 34200 2.32 223
5 22900 2.36 177
6 1321 246 114.8
7 931 247 95.9
8 518 251 68.3
9 346 273 49.1
10 122.9 1518 56
11 54.0 1549 39.9
12 84.6 1521 47
13 1249 107.4 94.2
14 1021 186 99.9
15 465 414 83.1
16 54.8 1302 35.9
Clculo Solucin: numrico en Ingeniera
Para linealizar la ecuacin basta tomar logaritmos. As:
log( Nu ) = log( ) + b log(Re) + clog(Pr)
Haciendo un ajuste lineal mltiple, se obtiene que
log( a ) = -0.1684 a = 10 -0.1684 = 0.6786
b = 0.5374; c = 0.2435
La representacin de los errores utilizando la ecuacin original y la obtenida por
el ajuste es la siguiente:
20
20
0 10
-20 0
-40
-10
-60
Nu
-20
-80
-30
-100
-40
-120
-140 -50
-1600 2 4 6 8 10 12 14 16 -600 2 4 6 8 10 12 14 16
puntos
En azul aparecen los errores de la correlacin de Sieder y Tate, y en rojo los

errores cuando se ajustan los datos experimentales.
RLM.2. Transferencia de calor en un reactor de lecho fluidizado

Dow y Jacob [Dow, W.M. y Jacob, M. Chem. Eng. Prog. 47, 637 (1951)]
propusieron las siguientes ecuaciones adimensionales para modelar la
transmisin de calor en un tubo vertical con un slido fluidizado con aire:
a2 a3 a a5
D D 1 C s 4 DG
Nu = a t t s t
1
L p
D g g Cg
Determine los parmetros a1-a5 y los errores cometidos.

h D
m t D Dt 1 C Dt G
t
kg L Dp g Cg g
469 0.636 309 833 256
913 0.636 309 868 555
1120 0.641 309 800 786
234 0.285 309 800 255
487 0.285 309 800 555
709 0.283 309 767 850
581 0.518 683 795 254
650 0.521 683 795 300
885 0.524 683 795 440
672 0.455 1012 867 338
986 0.451 1012 867 565
1310 0.455 1012 867 811
1190 0.944 1130 1608 343
1890 0.974 1130 1608 573
2460 0.985 1130 1608 814
915 0.602 1130 1673 343
1260 0.602 1130 1673 485
1690 0.617 1130 1673 700
h D
Nu = m t Nmero de Nusselt
kg
hm Coeficiente individual de transmisin de calor
Dt Dimetro de tubo
Dp Dimetro de partcula
L Longitud del lecho fluidizado caliente
G Velocidad msica del gas
k g , g , Cg , g Propiedades de la fase gas
Cs , s Propiedades de la fase slida
Solucin
Se debe en primer lugar linealizar la ecuacin. Basta para ello tomar logaritmos
en ambos lados de la ecuacin:
D D 1 C s DG
t t s t
log( Nu ) = log( a1 ) + a 2 log + a3 log + a 4 log + a5 log
D C
L p g g g
Utilizando ahora una rutina para regresin lineal mltiple:
Coeficientes Intervalo de confianza del 95%
log( a1 ) 0.0714 -0.3753, 0.5181
a
2 0.7402 0.6221, 0.8583
a3 0.3454 0.2648, 0.4260
a4 -5.704110-4 -0.1613, 0.1602
a5 0.7865 0.7082, 0.8648
Clculo
El resultado numrico
muestra quelos enIngeniera
intervalos de confianza para el primer parmetro y para el cuarto cortan el valor cero, lo que
indica que, posiblemente a1 se pueda fijar a uno y que la dependencia del cuarto trmino se
pueda despreciar, al menos en el intervalo de valores experimentales testado.
0.04
0.03
0.02
0.01
error [log(Nu)]
-0.01
-0.02
-0.03
-0.04
-0.050 2 4 6 8 10 12 14 16 18
puntos
RLM.3. Correlacin de los parmetros de la ecuacin de Antoine para el

calculo de la presin de vapor de un componente puro (propano)
La ecuacin de Antoine permite calcular la presin de vapor de un componente
puro en funcin de la temperatura.
B
A
P =10 T C
v
Linealice la ecuacin anterior y determine, utilizando regresin lineal
mltiple, los parmetros A, B y C para los siguientes datos de presin de
vapor del propano
Punto T (F) P (psia) Punto T (F) P (psia)
1 -70 7.37 11 30 66.3
2 -60 9.72 12 40 78
3 -50 12.6 13 50 91.8
4 -40 16.2 14 60 107.1
5 -30 20.3 15 70 124
6 -20 25.4 16 80 142.8
7 -10 31.4 17 90 164
8 0 38.2 18 100 187
9 10 46 19 110 213
10 20 55.5 20 120 240
Represente los valores calculados y experimentales en una grfica P vs T.
Represente los errores cometidos en P vs T para cada dato.
Solucin:
Para poder obtener los parmetros el primer paso es obtener una
expresin lineal. Para ello se debe reescribir la ecuacin
B
A
P = 10 T C B ;
v
; log( P ) = A +
v
log( Pv )(T + C ) = A(T + C ) + B
T+C
log( Pv )T + C log( Pv ) = AT + AC + B; log( Pv )T = ( AC + B ) + AT C log( Pv )
Recordando que en una regresin lineal mltiple
y = a0 + a1 x1 + a2 x2 + ... identificando trminos
y = log( Pv ) T ; x1 = T ; x2 = log( Pv )
a0 = AC + B ; a1 = A; a 2 = C
Haciendo el ajuste lineal mltiple se obtiene que:
a0 = 677.3026; a1 = 5.2258; a2 = -428.1183
Y por lo tanto:
A = 5.2258 ; B = -1560.0 ; C = 428.1183
Las representaciones de los puntos experimentales y los calculados a
travs de la regresin, as de los errores cometidos en cada punto son:
250 1
Calculados
Experimentales
0.8
200
0.6
Presin(psi
error(psia)
150 0.2
a)
100
0.4
0
50
-0.2
0 -0.4
-80 -60 -40 -20 0 20 40 60 80 100 120 0 2 4 6 8 10 12 14 16 18 20
Temperatura (F) punto
RLM.4. Capacidad calorfica del propano.

La capacidad calorfica del propano gas se puede expresar como un polinomio
de la temperatura:
C p = a0 + a1T + a 2T 2 + a3T3 + ...
Ajuste los datos de la tabla siguiente y a polinomios de diferentes grados y
comprelos entre s
Cp Cp
-1 -1
Punto T (K) (kJ kmol K ) Punto T (K) (kJ kmol-1 K-1)
1 50 34.06 11 700 142.67
2 100 41.3 12 800 154.77
3 150 48.79 13 900 163.35
4 200 56.07 14 1000 174.6
5 273.16 68.74 15 1100 182.67
6 298.15 73.6 16 1200 189.74
7 300 73.93 17 1300 195.85
8 400 94.01 18 1400 201.21
9 500 112.59 19 1500 205.89
10 600 128.7
Represente calores de Cp calculados y experimentales vs T en la misma
grfica
Nota: en este caso es necesario escalar los valores de T, para ello defina una
T
nueva variable: T* = (donde Tm es el valor medio de las temperaturas) de
Tm
tal forma que
Clculo numrico en Ingeniera + b T * + b T * 2 + b T * 3 + ...
1 2
C=
p 0 3
Y por lo tanto : = b1 ;a = b22 ; a = b3 ;...... a = bii
1 T 2 T 3 T3 i T
m m m m
Solucin:
De acuerdo la sugerencia del enunciado calculamos Tm =672.1742 K
Utilizaremos directamente T* para el resto del problema.
Los resultados los ajustamos a polinomios de diversos grados (entre 2 y
) los resultados se dan en la siguiente tabla, se incluye el valor de los
parmetros para un intervalo de confianza del 95%.
Polinomio grado 2 Intervalo de confianza del 95%
a0 17.7427 [14.3303 ; 21.1552]
a1 146.3905 [138.6433 ; 154.1376]
a2 -27.8513 [-31.2312 ; -24.4715]
Polinomio de grado 3 Intervalo de confianza del 95%

a0 20.5235 15.8850 25.1620
a1 132.0625 113.3747 150.7504
a2 -12.1304 -31.2679 7.0070
a
3 -4.5947 -10.1096 0.9203
Los resultados del ajuste se pueden ver a continuacin.

220 6
error polin. grado 2

200 Cp experimental 5 error polin. grado 3
Polinomio grado 2
180 Polinomio grado 3 4
160 3
140 2
error de Cp
120 1
Cp
100 0
80 -1
60 -2
40 -3
200 0.5 1 1.5 2 2.5 T* -40 0.5 1 1.5 2 2.5

T*
RLM.5. Correlacin de los coeficientes de actividad de una mezcla binaria

utilizando la ecuacin de Margules.
La ecuacin de Margules para la correlacin de de coeficientes de actividad en
mezclas binarias es:
1 = exp x22 (2 B A) + 2 x23( A B)
2 = exp x12 (2 A B ) + 2 x13( B A)
Donde x1, x2 son las fracciones molares de los componentes 1 y 2

respectivamente y 1 , 2 son los coeficientes de actividad. Los parmetros A y
B son constantes para cada mezcla binaria.
Las dos ecuaciones anteriores se pueden combinar para calcular la energa
libre de exceso de Gibbs.
g = GE = x ln( ) + x ln( ) = x x ( Ax + Bx )
RT 1 1 2 2 1 2 2 1
Para la mezcla Benceno (1) y n-heptano (2) estn disponibles los siguientes
datos:
Punto x1 1 2
1 0.0464 1.2968 0.9985
2 0.0861 1.2798 0.9998
3 0.2004 1.2358 1.0068
4 0.2792 1.1988 1.0159
5 0.3842 1.1598 1.0359
6 0.4857 1.1196 1.0676
7 0.5824 1.0838 1.1096
8 0.6904 1.0538 1.1664
9 0.7842 1.0311 1.2401
10 0.8972 1.0078 1.4038
Calcule los valores de los parmetros A y B
Solucin:
En este caso podemos reordenar la ecuacin que nos da la energa libre de
exceso de la siguiente manera:
2 2
g = GE = x ln( ) + x ln( ) = A x x + B x x = Az + Bz
RT 1 1 2 2 1 2 1 2 1 2
2 2
donde z = x x ; z = x x La dificultad en este caso es que el ajuste
1 12 2 1 2
lineal mltiple exige que la recta de correlacin pase por el origen.
Teniendo esto en cuenta se obtiene:
A = 0.2511; B = 0.4609
0.1 1.5 x 10-3
experimental
0.09 calculado
0.08 1
0.07
error(GE/RT
GE/(RT
0.5
0.06 0
)
0.05
0.04
0.03
-0.5
0.02
0.01 -1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
fraccion molar de benceno fraccion molar de benceno
RLM.6. Regresin de velocidades de reaccin cataltica heterognea.

La Tabla siguiente presenta la velocidad de reaccin para la siguiente reaccin
cataltica heterognea, representada por A B
PA( atm) PB ( atm) r 105
1 0.0010 5.10
0.9 0.1016 5.40
0.8 0.2021 5.55
Clculo numrico en Ingeniera0.7 0.3011 5.85
0.6 0.4004 6.00
0.5 0.5043 6.15
0.4 0.6040 6.30
0.3 0.7020 6.45
Se ha sugerido la siguiente ecuacin para correlacionar los datos:
k1 P A
r= donde k1, K A y KB son coeficientes que deben
1+ K A PA + K B PB 2
determinarse por regresin.
Utilizando regresin lineal, determine dichos coeficientes.
Solucin
La ecuacin anterior se puede lienalizar reordenando los trminos, invirtiendo y
tomando la raz cuadrada.
P 12 1 KA KB
A
= + PA + PB = a1 + a2 PA + a3 PB
r k 1 k 1 k 1
El resultado es: a1 = 177.8217; a 2 = 317.5511; a3 = 215.8415 y por lo tanto
k1 = 1 2
= 3.1625 105
177.8217
K A = a 2 k1 = 1.7858
K B = a3 k1 = 1.2138
1
0.8
0.6
0.4
0.2
(PA/r)0.5
-0.2
-0.4
-0.6
-0.8
-11 2 3 4 5 6 7 8
Puntos
ENLEcuacione.1.s noDimetrolinealesde una tubera a travs del flujo

Utilizando la ecuacin de Chen, un alumno de Ingeniera qumica obtuvo un
factor de friccin ( f ) de 0.01 para una tubera con rugosidad relativa ( ) 0.006.
Si por la tubera circulaba agua con un caudal de 0.01 m 3/s. Cul era el
dimetro de la tubera?. (Utilice diferentes mtodo numricos para calcularlo)
Ecuacin de Chen:
1.1098
1 = 4 log 5.0452 log + 5.8506
0.898
f 3.7065 Re 2.8257 Re
Solucin:
Se debe reescribir la ecuacin igualndola a cero, por ejemplo:
1 5.0452 1.1098 5.8506
+ 4 log log + =0
0.898
f 3.7065 Re 2.8257 Re
Y aplicar cualquier mtodo numrico (biseccin Newton, secante, etc)
D = 0.3869 m.
ENL.2. Ecuacin de Van der Waals
La ecuacind e Van der Waals para un gas es:
P+ A
V B = R T donde la P es la presin (atm). V es el volumen
2
V
molar (l/mol). T es la temperatura absoluta (K) y R=0.082054 atm l /mol K.
Hallar el volumnen molar del CO2 a 500 K y 10 atm, utilizando todos los
mtodos estudiados.
Nota: A= 3.592 l2 atm / mol2 B= 0.04267 l / mol. Solucin
Para los mtodos de punto fijo es necesario despejar el volumen molar del gas
de alguna manera, por ejemplo:
RT
V= A +B
P +
V2
(Tenga en cuenta que no siempre converger, depende de cmo
despeje el volumne). Resolviendo por diferentes mtodos se obtiene: V=
4.055 l
Para mtodos como el de Newton la ecuacin hay que reorganizarla de forma
que quede igualada a cero
P + A V B R T = 0
V2
ENL.3. Evaporador Flash
Una mezcla de composicin (fraccin molar) como la de la tabla se introduce
en un evaporador flash a la temperatura de 160 C y 1 atm. Determinar que
fraccin de mezcla de ha evaporado. Suponer mezcla ideal.
comp j zj Aj Bj
1 0.1 -3848 8.358
2 0.2 -4002 7.833
3 0.2 -4997 8.665
4 0.2 -5595 8.847
Clculo numrico en Ingeniera5 9.261 0.2 -6301
6
10.078 0.1 -7296
*
Ln K j =
Aj
T + 460

+ Bj Kj=
Pj
P
(T C)
Nota: la forma ms fcil de resolver el problema es resolver
nx j y j = 0 donde xj, yj son funciones de la fraccin evaporada.

j1
Calcule tambin las temeraturas e roco y burbuja de la mezcla.
Solucin.
Este problema se puede resolver planteando el conjunto completo de
ecuaciones o bien reducirlo a una nica ecuacin que dependa slo de la
fraccin de vapor. Vemos como se puede llegar hasta este punto.
V
Haciendo un balance de material.

F=L+V
Fz j = Lx j + Vyj
V
Defiendo la fraccin de vapor como =
F
z j = (1 )x j + yj
Como y j = K j xj , sustituyendo en la ecuacin anterior se tiene que:

z
j
xj =
(1 ) + K j
Como sugiere el enunciado se puede resolver una sola ecuacin de la
forma
n
n
j
n
j
zj
zj
z j (1 K j ) j
y =0= K
(1
x =
j1 j1 )+Kj (1 ) + K j j1 (1 ) + K j
El resultado es:
= 0.6876
comp j zj xj yj
1 0.1 0.0161 0.1381
2 0.2 0.0658 0.2610
3 0.2 0.1272 0.2331
4 0.2 0.2251 0.1886
5 0.2 0.3382 0.1372
6 0.1 0.2276 0.0420
ENL.4. Reacciones de equilibrio
Sea el siguiente conjunto de reacciones:
A+ B C+D
C+B X+Y
A+ X Z
Calcule las concentraciones de equilibrio de cada una de las especies si
inicialmente se introduce una concentracin de:
CAo = 1.5 M CBo = 1.5 M
Datos:
Las constantes de equilibrio para el sistema son:
K 1 =CC CD = 1.06 K 2 = CX CY = 2.63 K 3= CZ = 5 (l/mol)
CA C B CC CB C A CX
Resuelva el sistema de ecuaciones por el mtodo que considere oportuno.
Justifique su eleccin. Tome valores iniciales razonables para cada una de las
variables. (No use el valor cero como valor inicial).
Solucin:
Quiz la mejor forma de resolver el problema sea utilzar grados de avance de
cada ecuacin. As, si i es el grado de avance de cada una de las reacciones
se tiene que:
C A = CA0 1 3
C B = CB0 1 2
CC = 1 2
CD = 1
CX = 2 3
CY = 2
CZ = 3
Para evitar divisiones por cero conviene pasar los trminos que estn
dividiendo multiplicando y dejar el sistema de ecuaciones de la siguiente
manera:
CC C D 1.06 C A CB = 0
C X CY 2.63 CC CB = 0
C Z 5 C A CX = 0
El resultado, aplicando el mtodo de Newton es:
1 = 0.7053; 2 = 0.5518; 3 = 0.3740
NOTA: Los valores iniciales son muy importantes en este caso. Conviene
comprobar que el resultado final de concentraciones lleva a soluciones con
sentido fsico, porque desde el punto de vista puramente matemtico el
sistema tiene diferentes soluciones.
ENL.5. Circulacin de fluidos en una red de tuberas

Agua a 25 C est circulando por la red de tuberas que se muestra en la figura.
La presin a la salida de la bomba es de 15 bares (15 105 Pa) sobre la
atmosfrica, y el agua se descarga a presin atmosfrica al final de la lnea.
Todas las tuberas son de 6 catlogo 40 (0.154 m de dimetro interno), con
ClculorugosidadnumricorelativaenIngeniera0.006. Las longitudes
equivalentes de todas las tuberas conectando los diferentes nodos son
L01 = 100 m; L12 = L23 = L45 = 300 m; L13 = L24 = L34 =1200m
Calular los flujos y presiones en los nodos 1, 2, 3 y 4 para la red. Asuma que
el factor de friccion es constante e igual a 0.005 para todas las tuberas.
Informacin adicional y notas:
Para la solucin del problema es conveniente expresar la prdida de
presin en cada nodo desde i hasta j como:
P
i,j=
k
i,j qi, j 2
Donde qi, j es el flujo volumtrico y Pi, j es la cadida de presin. El trmino
ki, j est relacionado con los factores de friccin de fanning y la velocidad
media de los fluidos a travs de la ecuacin
L
i, j
2f
k = D = 32 f Li, j
D / 4 2 D5
i, j
2 2
0 q01 1
q12
q13
5 4 q24
2
q45
q23
q34 3
Repita el problema pero ahora usando una ecuaci adecuada para calcular
el factor de friccin.
Solucin:
Hay dos relaciones que gobiernan el proceso.1. La suma de los flujos en
cada nodo dbe ser cero. 2. La suma de todas las caidas de presion en un
bucle cerrado debe ser cero.
q01 q12 q13 = 0
q12 q24 q23 = 0
q23 + q13 q34 = 0
q24 + q34 q45 = 0
P +P =P
P01 + P12 + 24 45 bomba
P13 P23 P12 = 0

P + P
23 34 P24 = 0
Que junto con la ecaucin que relaciona los flujos y las caidas de presin
define completamente el sistema de ecuaciones.
P
i,j=
k
i,j qi, j 2
Dado que los flujos en m3/s pueden ser muy pequeos en comparacin con
las presiones combiene escalar el problema, por ejemplo trabajando en kPa
y en litros/s
Caso a:
Tramo Q (l/s) P( kPa)
0-1 0.3118 1.8194
1-2 0.2059 2.3812
1-3 0.1059 2.5165
2-3 0.0491 0.1352
2-4 0.1569 5.5261
3-4 0.1549 5.3909
4-5 0.3118 5.4583
Caso b. Debemos ademas utilizar una ecuacin adecuada para los
factores de friccin. Por ejemplo la ecaucin de Coolebrock y White.
Ahora los factores de friccin son tambin incognitas de nuestro sistema.
Nuevamente para facilitar la resolucin numrica es conveniente escalar
las variables. En este caso se puede escalar por un factor de 100 el
factor de friccion.
1
D 1.256
= 4log +
f 3.7 Re f
Tramo Q (l/s) P( kPa) f.103

0-1 0.3915 0.2809 0.2809
1-2 0.2682 0.3159 0.3159
1-3 0.1232 0.4160 0.4160
2-3 0.0697 0.5248 0.5248
2-4 0.1985 0.3494 0.3494
3-4 0.1929 0.3530 0.3530
4-5 0.3915 0.2809 0.2809
ENL.6. Velocidad terminal

Determinar la velocidad terminal de una partcula slida cayendo en un fluido.
Considere una partcula esfrica uniforme de hierro cayendo en aire.
Datos para la partcula:
p = 7860 kg / m3
D p = 0.5 mm
Datos para el aire
= 1.23 kg / m3
= 1.79 10 5 Pa s
Notas: La velocidad terminal de una partcula esfrica cayendo en un medio
fluido se puede determinar por la siguiente expresin
4g
p Dp
vt =
3 CD
ClculoElcoeficientenumricodeen rozamiento Ingeniera depende del
nmero de Reynolds de la sigueinte manera:
24
D = Si Re < 0.1
Re
24 0.7
) Si 0.1 Re 1000
D = (1 + 0.14 Re
Re
vD
Donde Re = t p

Solucin. Realmente basta con resolver una ecuacin con una incognita que es
la velocidad terminal.Suponiendo un nmero de Reynolds maor que 0.1, se
obtiene vt = 7.586 m / s . Debemos adems comprobar el nmero de Reynolds
que es a la velocidad terminal 260.7
ENL.7. Grfica de Moody

Dibuje la grfica de Moody utilizando la ecuacin de Coolebrok para el factor de
friccin en rgimen turbulento
Solucin: La grfica de moody representa el Nmero de Reynolds frente a 4f .
Reordar que
f = 16 Si Re < 2100
Re

1 = 4log D + 1.256 Si Re 5000
f 3.7 Re f
100
10-1
4f
10-2
102 103 104 105 106 107 108
Re
ENL.8. Temperatura de llama

Cuando se quema gas con aire, la temperatura mxima terica que se puede
alcanzar en la llama, es lo que se conoce como temperatura de llama
adiabtica. Esta temperatura depende de la composicin del gas as como de
la cantidad de aire con la que este se queme. Determine la temperatura de
llama adiabtica cuando se quema gas natural. Considere los casos los que
la relacin aire/gas natural vara entre 0.5 y 20.
Composicin del gas natural:
CH4: 90% en volumen, C2H6 : 6% en volumen; N2: 4% en volumen.
Entalpa combustin: CH4: -212798 cal/mol; C2H6: -372820 cal/mol.
La capacidad calorfica de los diferentes gases imp se puede expresar
por:
2
Cp = a + bT + cT ( T en K y Cp en cal/mol K).
a b c
CH4 3.381 18.044 -4.3
C2H6 2.247 32.201 -11.049
CO2 6.214 10.396 -3.545
H2O 7.256 2.298 0.283
O2 6.148 3.102 -0.923
N2 6.524 1.25 -0.001
La temperatura ambiente y de todas las entradas al quemador es de 25C.
Solucin:
Para determinar la temperatura de llama adiabtica debemos realizar un
balance de materia en funcin de conversin y un balance de energa.
Las reacciones de combustin son las siguientes:
CH 4 + 2O2 CO2 + 2 H 2O
7
C 2 H 6 + 2O2 2 CO2 + 3 H 2O
Nota: Si no hay oxigeno suficiente para la combustin se forman productos de
de combustin incompleta. Sin embargo vamos a suponer que el metano y el
etano no reaccionaran mas y simplemente se calentaran, adems no hay
preferencia por ninguna de las reacciones y ambas se dan en la misma
extensin (lo que nos permite tratarla como si fuera una sola reaccin)
La composicin del gas natural (fraccin molar = fraccin en volumen)
x
CH = 0.75;xC H = 0.23; xN = 0.02
4 2 6 2
Y la composicin del aire
xO2 = 0.21; xN2 = 0.79
aire aire
Llamando r a la relacin aire/gas natural y tomando como base de clculo 1
mol de gas natural las cantidades iniciales de cada componente en la mezcla
sern:
El oxgeno necesario para la combustin completa ser:

7
O = 2 CH + C H
2 combustion completa 4 0 2 2 6 0
El oxigeno consumido por cada reaccin (suponiendo que la reaccin fuese
completa) sera:
(O2 ) 0 (CH4 ) 0
O2 reac1 = 2
2 (CH 4 ) 0 + 3.5(C 2 H6 ) 0
(O2 ) 0 (C 2 H6 ) 0
O
2 reac2 = 3.5 2 (CH 4 ) 0 + 3.5(C 2 H6 ) 0
LaClculoconversinnumricodependerenIngenierade la cantidad de oxigeno disponible. As:
0.5O2 reac1
Si O2 combustion completa < =
CH4 0
En otro caso =1
Donde es la conversin.
Los moles finales se pueden calcular en funcin de la conversin de la
siguiente manera:
CH 4 = CH4 0 1
C 2 H 6 = C 2 H6 0 1
CO2 = CO2 0 + CH 4 0 + 2 C 2 H6 0
H 2 O = ( H 2 O ) 0 + 2 CH 4 0 + 3.5 C 2 H6 0
O2 = (O2 ) 0 2 CH 4 0 3.5 C 2 H6 0
N2=(N2)0
Finalmente de un balance de energa
H + H
C , CH CH 4 CH 4 C,CH C2H6 C 2 H6 +
4 0 2 6 0
( moles finales ) i TT C p (T ) i = 0
i 0
Si resolvemos utilizando un mtodo de Newton la ecuacin anterior para

diferentes valores de r se obtiene. Como esa de esperar la temperatura
mxima se obtiene justo cuando se aade la cantidad estequiomtrica de aire.
(primer punto de conversin total).
2400 1
2200 0.9
2000 0.8
conversindemetan
1800 0.7
Temperatura
1600 0.6
o
1400 0.5
K
1200 0.4
1000 0.3
800 0.2
600 0.1
400 0 2 4 6 8 10 12 14 16 18 20 0 0 2 4 6 8 10 12 14 16 18 20
relacion aire/gas natural relacion aire/gas natural
EcuacionesLIN.1.Balanceslinealesde materia en estado estacionario ara una

secuencia de separacin.
El p-xileno, estireno, tolueno y benceno se separan en la siguiente secuencia
de columnas de destilacin.
D1 7% Xileno
4% Estireno
54% Tolueno
35% Benceno
15% Xileno 18 % Xileno

25% Estireno 24 % Estireno
40% Tolueno B1 42 % Tolueno
20% Benceno 16 % Benceno
15 % Xileno
D2 10 % Estireno
F = 70 kmol/min 54 % Tolueno
21 % Benceno
24 % Xileno
B2 65 % Estireno
10 % Tolueno
1 % Benceno
a) Calcula el flujo molar de las corrientes D1, D2, B1 y B2.

b) Reduce el flujo original de la corriente de alimento en un 1% y
luego en un 2% y repite el apartado (a). Explica los resultados.
c) Determina el flujo molar y las composiciones de las corrientes B y
D del apartado (a).
LIN.2. Recuperacin de Cromato potsico cristalino
A continuacin se muestra el diagrama de flujo de un proceso en rgimen
estacionario para la recuperacin de cromato potsico cristalino, a partir de una
disolucin acuosa de esta sal.
Agua
Slido
4500 kg/h
Evaporador Cristalizador
Disolucin
R=36.36% peso cromato potasico
Se combinan 4500 kg/h de una disolucin que contiene 1/3 (fraccin msica)
de cromato potsico con un flujo de recirculacin que contiene 36.36% en peso
de cromato. El flujo combinado se alimenta a un evaporador. El flujo
concentrado que sale del evaporador contiene un 49.9% en peso de cromato.
EsteClculo numrico flujo se
enIngeniera
alimenta un cristalizador donde se enfra (produciendo la
precipitacin de cristales de cromato potsico). Una parte se filtra y se recircula. La otra parte
que es una mezcla de slido y disolucin al 36.36% en peso se extrae del sistema. Se sabe
que el 95% en masa de toda la parte extrada del sistema es cromato slido. Plantee un
sistema de ecuaciones lineal y calcule todas las variables desconocidas del sistema
EcuacionesEDO.1. diferencialesCinticadereaccinordinariasreaccin
Calcule, utilizando los mtodos numricos que considere oportunos, la
variacin con el tiempo de las especies A, B y C en la siguiente reaccin.
1
k k k = 0.5min 1 A(0) = 1; B (0) = 0; C (0) = 0 t 0,10
A B C
1 2
1
k =1.0min
2
Solucin: Las ecuaciones que hay que integrar son las siguientes:
dC
A = k C
dt 1 A
dC
B =kC k C Con las condiciones de contorno dadas en el
dt 1 A 2 B
dC
C= kC
dt 2 B
enunciado.
El resultado es:
1
A
0.9 B
C
0.8
0.7
Concentracin
0.6
0.5
0.4
0.3
0.2
0.1
00 1 2 3 4 5 6 7 8 9 10
tiempo (s)
EDO.2. Reactor de lecho fluidizado

En un reactor de lecho fluidizado de 0.8 m de longitud y 1 m de dimetro se
lleva a cabo una reaccin cataltica. El modelo de flujo aplicable es el de
partculas de tamao intermedio, con lo que las concentraciones de reactivo A
en la fase burbuja (CAb) y en la emulsin (CAe) siguen las ecuaciones:
d CAb
- ub = b k CAb + Kbe( CAb - CAe )
dz
d CAe
-(1- )umf = (1- )(1- mf )k CAe - Kbe ( CAb - CAe )
dz
En estas ecuaciones:
z es la distancia hasta el fondo del lecho (altura)
k es la constante cintica de la reaccin (igual a 3 seg-1 en nuestro
caso) B es la fraccin de slidos en la estela (igual a 10-3 en nuestro caso)
Clculomfes lanumricovelocidadenIngenieramnima de fluidizacin (igual a 0.21 m/s en nuestro caso)
mf es la porosidad el punto de mnima fluidizacin (igual a 0.45 en nuestro

caso)
b es la velocidad de ascenso de la burbuja, que es una funcin de su
dimetro:
ub = 1.6 u o umf + 113.d b0 . 5 d t1. 35 + ubr , m / s
u0 es la velocidad de alimentacin del gas (igual a 0.4 m/s en nuestro

caso) ubr es la velocidad de ascenso de una burbuja individual:
ubr = 0,711(g db )0,5
dt es el dimetro del reactor (igual a 1 m en nuestro caso)
db es el dimetro de una burbuja, que es una funcin de la posicin:
1,11 0.81
db = 2,3(u0 - umf ) z [SI]
es la fraccin de burbujas, que puede calcularse con la expresin:
-
= u0 umf
ub
Kbe es el coeficiente de transferencia de materia desde la burbuja a la emulsin
que debe calcularse mediante el uso de K bc (burbuja-nube) y K ce (nube-
emulsin):
1/Kbe = (1/Kbc)+(1/Kce)
umf D0,5 g1/ 4
K bc = 4,5 +5,85
[s -1 ]
db d b5 / 4
0,5
D mf ubr
K = 6,77 [ s-1 ]
ce db 3
D es la difusividad del reactivo (9 10-5 m2/s en nuestro caso).
Se desea conocer cmo varan a lo largo de la altura del reactor:
a) las concentraciones CAb y CAe, teniendo en cuenta que a
la entrada son ambas iguales a 0.1 mol/L
b) los coeficientes de transferencia Kbc y Kce
c) el dimetro y la velocidad de las burbujas
Solucin: Basta resolver las ecuaciones diferenciales del enunciado, y luego
sustituir para encontrar cada una de los parmetros que se piden:
0.1
Concentracion en fase burbuja
0.09 Concentracion en fase emulsion
0.08
(mol/L)
0.07
0.06
ntracio
Conce
0.05
0.04
0.03
0.02
0.01
00 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

Altura en el reactor
0.35 3
0.3
2.5
0.25
2
velocidad
diametro
0.2
1.5
0.15
1
0.1
0.5
0.05
00
00
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
altura altura
105
Kbc
Kce
104
103
102
101
100
10-1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

altura
EDO.3. Reactor encamisado exotrmico

EnClculo numrico un enIngeniera
reactorencamisado
se lleva a cabo una reaccin exotrmica. Para este tipo de sistema se comprueba mediante
balances de materia y energa que se cumple la siguiente ecuacin:
dT 4UT Tc
d = Tad dt CP k1
donde es la conversin en el reactor (cero para = 0 ) y Tc es la temperatura
del encamisado. A la entrada del sistema ( = 0 ) la temperatura es de 340 K.
En el sistema, de primer orden, se cumple la ecuacin de diseo:
d
= k1
d
en la que la constante de velocidad sigue la ley de Arrhenius,
E
k = k0 exp
RT
y se dispone de los siguientes datos:
T (K) k (s-1)
300 0.00012239
320 0.0013828
330 0.003864
340 0.0107844 En
el caso particular que se trata:
- la temperatura inicial del sistema de reaccin es 340 K,
- se midi el valor de Tad =146 K ,
- se conoce 4U = 0.2 s1
dt CP
En este ejercicio:
1) Se desea seguir la evolucin T- para tres temperaturas de encamisado
distintas, 330, 335 y 337 K. Realiza los programas necesarios para
visualizar las diferencias en una misma grfica.
2) Tambin se quiere seguir la evolucin T- en una misma grfica para
los tres casos presentados.
NOTAS:
1. Tenga en cuenta que en la ecuacin 1 la variable independiente es la
conversin, y en la 2 es el tiempo. Cuando integre ambas deben
depender del tiempo
2. Integre para los primeros 200 segundos.
Solucin:
La primera parte requiere estimar los valores de la Energa de activacin y el
factor pre-exponencial de la ecuacin cintica. Para ello se puede llevar a cabo
1
un ajuste lineal mltiple representando ln( k ) vs.
T
k 0 = 3.99551012 s1
Haciendo esto se obtiene que: E
= 11403.42 K
R
-4.5
Experimental
-5 Regresion
-5.5
-6
-6.5
Ln(k)
-7
-7.5
-8
-8.5
-9
-9.5
2.9 2.95 3 3.05 3.1 3.15 3.2 3.25 3.3 3.35
1/T x 10-3
El siguiente paso es resolver el sistema de ecuaciones diferenciales ordinarias,

simplemente poniendo ambas ecuaciones en funcin del tiempo:
dT 4U T Tc
= k 1 Tad
d d t C Pk 1
d
= k1
d
El resultado final
355 1
0.9
350 0.8
0.7
Temperatura K
345 0.6
conversion
0.5
340 0.4
0.3
335 0.2
0.1
3300 20 40 60 80 100 120 140 160 180 200 0 0 20 40 60 80 100 120 140 160 180 200
tiempo (s)
tiempo (s)
355
350
Temperatura K
345
340
335
3300 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
conversion
EDOClculo.4numrico.DescomposicinenIngeniera de SH2 en una partcula esfrica
Se pretende descomponer el sulfuro de hidrgeno en sus elementos. Para ello
se dispone de un sistema en el que se difunde SH 2 dentro de las partculas
esfricas de un catalizador poroso de slica/almina. El sistema se caracteriza
por la ley de Fick, que aplicada a un sistema con una reaccin de primer orden,
como es el caso, tiene la forma:
d2C A
De kCA
2
= 0 dx
donde De representa la difusividad efectiva, CA la concentracin de SH2, k la
constante de velocidad de reaccin y x la distancia desde el centro de la
partcula.
CA= conocido
CA presenta un mnimo
x=0 x=L
Las condiciones de contorno de la ecuacin anterior son:
1. La concentracin de SH2 en la superficie del catalizador es conocida, e
igual a 10 mol/m3.
2. La concentracin de SH2 en el centro de la partcula presenta un mnimo,
por lo que dCA/dx=0 en el centro.
Se desea determinar el perfil de concentracin de SH2 a lo largo de la partcula.
Datos:
De= 510-6
m2/s k=1 s-1
Radio de la partcula= 2 mm.
Solucin: El primer paso es transformar la ecuacin diferencial de segundo
orden en un sistema de ecuaciones diferenciales de primer orden:
dCA x = 0 CA = 10
=V
dx con condiciones de contorno: dC
dV kCA x=L V= A
=0
= dx
dx De
El problema es que ahora las condiciones de contorno vienen dadas en dos
puntos diferentes. Para solucionar el problema se puede utilizar la tcnica
conocida como Shooting method que consiste en suponer uno de los valores
iniciales, integrar la ecuacin diferencial, y comprobar si se satisface el otro
valor inicial. Si no es as el valor supuesto se modifica hasta que el segundo
valor inicial se cumpla. Para hacer esto de forma sistemtica se puede
introducir la ecuacin diferencial dentro de un mtodo para resolver
ecuaciones.
En este caso es conveniente suponer la concentracin en x=L. La razn es que
el valor de la concentracin en dicho punto es fcilmente acotable (estar
siempre entre 0 y 10), lo que nos permite utilizar un mtodo como el de la

biseccin que podemos asegurar que va a converger iempre. Sin embargo de
la otra alternativa, suponer el valor de la derivada de la concentracin en el
punto inicial, slo sabemos que debera ser negativ . El resultado es el
siguiente:
Valor de la concentracin en x=L, C A( x = L) = 7.0058
10.5 3500
10 3000
9.5 2500
Concentracin
aC
9 2000
Derivad
8.5 1500
a
8 1000
7.5 500
70
0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Radio -3 Radio -3
x 10 x 10
EDO.5. Cintica de descomposicin formacin

La concentracin de un reactivo A vara con el tiempo segn la siguiente
ecuacin:
dC A 3500
= 9420exp CA t en minutos, CA en mol/L
dt 298
que no depende de la temperatura. Se realiza un experimento y al cabo de 25
minutos se encuentran en el sistema 3.4 mol/L de A. Cual era la
concentracin inicial? Dibuje la variacin de la concentracin con el tiempo
para este reactivo.
b) La concentracin de un reactivo B vara con el tiempo segn la siguiente
ecuacin:
dC B 15000 2
= 51015 exp CB t en minutos, CB en mol/L, T en
dt T
K
y la temperatura vara con el tiempo segn la ecuacin:
T = 298+ 5t t en minutos, T en K
En el sistema se parte de una concentracin de B de 7 mol/L. Cul es la
concentracin al cabo de 25 minutos? Dibuje la variacin CB-t.
c) Un tercer reactivo C est presente en el sistema. La concentracin de C es
una funcin de las ecuaciones que describen la variacin de A y B, de forma
que:
dCC = 0.05(C + C
A B
) dt
La concentracin inicial de C es 20 mol/L. Dibuje la variacin de CC con el
tiempo para los 25 primeros minutos
Solucin:
En este caso tenemos primero una ecuacin diferencial, cuya nica dificultad
es que se debe integrar hacia atrs (incrementos de tiempo negativos). El
resultado es
25
20
15
m
C
A
L
/
l
10
00 5 10 15 20 25
tiempo (min)
El valor de la concentracin de A a tiempo cero es de : 22.00 mol/L.

El segundo apartado slo requiere tener en cuenta que la temperatura vara de
forma lineal con el tiempo. El resultado es:
7
4
CB mol/L
00 5 10 15 20 25
tiempo (min)
Por ltimo, en el tercer caso se tiene un sistema de ecuaciones diferenciales

ordinarias donde se debe utilizar el valor inicial de concentracin de A,
calculado en el primer apartado. El resultado es
25
CA
CB
CC
20
15
mol/L
C
10
0 0 5 10 15 20 25
tiempo (min)
EDO.6. Reactor cataltico de membrana

La deshidrogenacin de un compuesto A tiene lugar en un reactor con
membrana selectiva bajo condiciones isotermas.
A B1C
2
La reaccin es reversible con constante de equilibrio KC .
La ventaja de utilizar un reactor de membrana est en que uno de los
productos atraviesa de forma selectiva la membrana, de tal manera que si el
componente que atraviesa la membrana es uno de los productos, la reaccin
se puede llevar hasta conversin completa.
Permeado H2
H2
Particulas de
catalizador
Alimento
Efluente
La velocidad de transferencia hasta la membrana depende no slo de la

resistencia ofrecida por la membrana, si no tambin de las capas de
concentracin en cada lado de sta. El coeficiente de transferencia de materia
a travs de la membrana se puede correlacionar a travs de la siguiente
ecuacin:
1 1
v 2 v 2
kC ( v ) kC ( v0 ) kC0
v v
0 0
Si el paso limitante es el flujo a travs de la membrana entonces el coeficiente

se puede considerar independiente de la velocidad, pero este no es el caso.
El componente C es transportado a travs de la membrana, se puede calcular
como:
rC ( membrana ) kC a CC para valores de perneado C bajos.
CalculeClculonumrico losflujos en Ingeniera
molare, en
funcin de la masa de catalizador, as como la prdida de presin para los siguientes valores:
dm 3 dm 3 kmol 0.5
k = 0.5 ; v0 = 50 ; KC = 0.5
kg min min dm3
P = 10 atm; F = 10 kmol ; = 0.002 atm
0 A0 min kg
2
dm dm
k = 0.1 ; a=2 ; W = 200 kg
C0 min kg
Solucin: Las ecuaciones para este sistema seran las
siguintes. Balances de Materia:
dF
A = rA'
dW
dFB
=
rA' dW
dFC = 1 r ' k
A C
aCC dW 2
Ecuaciones de velocidad
r' C C1/2
=kC B C
A A
KC
Relacin flujo concentracin y estequiometra
C = FA ; C = F B ; C = FC
A v B v C v
F+F+F P F+F+F 1 P
v = v0
A
F
B C 0
= v0
A B
F
C
con y =
A0 P A0 y P0
Cada de presin
F + F + FC
A F B
dy = A0
dW 2y
Coeficiente de transferencia de materia
v 12
k =k
C C0
v
0
La integracin de las cuatro ecuaciones diferenciales en funcin de la masa de
catalizador (entre 0 y 200 kg) lleva a:
10 1
FA
9 FB
FC 0.95
8
perdidadepresionP/P0
0.9
flujomolarkmo
l/min
5 0.85
4
0.8
3
2
0.75
1
0 0 20 40 60 80 100 120 140 160 180 200 0.7 0 20 40 60 80 100 120 140 160 180 200
masa de catalizador masa de catalizador
EDO.7. Asociacin de RCTA isotermos

Considere que la siguiente reaccin elemental e irreversible que tiene lugar en
fase lquida.
A+BC
Se lleva a cabo en una asociacin de 3 reactores en serie de tipo tanque
agitado. Las especies A y B se adicionan de forma individual al primer tanque
con un flujo volumtrico (cada una de ellas) de 6 dm3/min. El volumen de cada
reactor es de 200 dm3. Considere que inicialmente cada reactor est lleno con
un disolvente lquido inerte. La concentracin inicial de cada reactante es
3
C A 0 = C B0 = 0.5 mol3 y la constante de reaccin es de k = 0.5 dm
dm mol min
a) Calcule las concentraciones en estado estacionario a la salida de cada
reactor
b) Calcule el perfil de concentraciones hasta que se alcanza el estado
estacionario
c) Repita el apartado b suponiendo que ahora el componente B se divide
en tres partes iguales que se alimentan a cada reactor.
C C
A1 A3
C
A2
C C C
A0 B1 B3
C
B2
C
B0
Solucin: La primera parte consiste en resolver un sistema de ecuaciones no

lineales:
vC vC kCC V
0AA0 1 A1A1 B1
=0 1
V1
v C v C k C
0B B0 1 B1
C V V
A1 B1 1 = 0 1
vC v C k C C V
1 A1 2 A2 A2 B2 2
= 0 V2
vC v C k C C
1 B1 2 B2 A2 B2
V
2 = 0 V2
v C v C k C C
2 A2 3 A3 A3 B3
V
3 = 0 V3
v C v C k C C V
2 B2 3 B3 A3 B3 3 =0
V3
C A 0 = C B0 = 2.0; k = 0.5; V1 = V2 = V3 = 200
Con
v0 A = v0 B = 6; v1 = v2 = v3 = 12
C A1 = CB1 = 0.29167
Se obtienen los siguientes valores: C A 2 = CB2 = 0.13644
C A 3 = CB3 = 0.08132
En la segunda parte debemos integrar las siguientes ecuaciones:
dC v C v C k C C V
A1 = 0A A0 1 A1 A1 B1 1
dt V1
dC v C v C k C C V
B1 = 0B B0 1 B1 A1 B1 1
dt V1
dC v C v C k C C V
A2 = 1 A1 2 A2 A2 B2 2
dt V2
dC B 2 v1C B1 v2 C B 2 k C A 2 C
=
V
B2 2 dt V 2
dC v C v C k C
A3= 2 A2 3 A3
C V
A3 B3 3 dt V3
dC v C v C k C
B3= 2 B2 3 B3 A
C V
3 B3 3 dt V3
0.35
CA1
CA2
0.3
CA3
0.25
concentracion
0.2
0.15
0.1
0.05
0 0 5 10 15 20 25 30 35 40 45 50
tiempo
En el tercer caso debemos modificar las ecuaciones referentes a CB2 y CB3 de

la siguiente manera;
v C +v C v C k C C V
dC B2 = 0 B B 0 1 B1 2 B 2 A 2 B2 2
dt V2
v C +v C v C k C C V
dC 0B B0 2 B2 3 B3 A 3 B3 3
B3 =
dt V3
Pero ahora con v0 B = 2; v1 = 8; v2 = 10; v3 = 12
1
CA1
0.9 CA2
CA3
0.8
CB1
CB2
0.7
CB3
0.6
concentracion
0.5
0.4
0.3
0.2
0.1
0 0 5 10 15 20 25 30 35 40 45 50
tiempo
EDO.8. Tiempo de vaciado de un depsito

Eldepsito de la figura,inicialme nte lleno de agua, se descarga (vaca) a la
atmsfera a travs del tubo inferior. Calcule y

represente grficamente, para las dos primeras
D
horas de vaciado: 1
(1)
Variacin de la altura (H) del agua con el

tiempo dentro del depsito.
Ho
Variacin de la velocidad de salida con el
tiempo (punto 2).
Tenga en cuenta que la velocidad en cada instante
L
se puede calcular a partir de un balance de energa
mecnica:
2 (2)
V2 2L D
2 g H + 2 f V 2D = 0
donde: V2 = velocidad en la seccin de salida (2). g =

constante de gravitacin universal
f = factor de friccin
L = longitud de tubera D
= Dimetro de tubera
H = altura del agua contando desde el punto (2)
el factor de friccin se puede calcular a partir de la ecuacin de Coolebrook y

White:

1 = 4log
D +
1.256
f 3.7 Re f
donde = rugosidad
Re = nmero de Reynolds.
Datos:
= 10-3 Pas (Viscosidad del agua)
= 1000 kg/m3 (Densidad del agua)
= 5 10-4 m (Rugosidad)
L=1m Longitud de tubo
D = 0.5 cm dimetro del tubo en el punto (2)
Ho = 10 m Altura inicial del agua medida desde el punto
(2).
D1 = 0.2 m Dimetro del depsito.
Adems, utilizando la ecuacin de continuidad se tiene que la variacin de

altura con el tiempo viene dada por la expresin:
dH
V S =S
2 2 1 dt
donde S1 = seccin en el punto (1)
S2 = seccin en el punto (2)
Solucin: En este problema se debe resolver una ecuacin diferencial ordinaria.
Sin embargo, en cada paso de la integracin se debe adems resolver un
sistema de ecuaciones no lineales que permitan determinar el valor de la
velocidad y el factor de friccin.
10 3.5
9
3
8
7 2.5
6
2
altura(m)
Velocid
ad(m/s
)
4 1.5
3
1
2
1 0 1000 2000 3000 4000 5000 6000 7000 8000 0.5 0 1000 2000 3000 4000 5000 6000 7000 8000
tiempo (s) tiempo (s)
0.0263
0.0262
0.0261
factor de friccin
0.026
0.0259
0.0258
0.0257
0.0256 0 1000 2000 3000 4000 5000 6000 7000 8000

tiempo (s)
EDO.9. Simulacin del enfriamiento del agua dentro de un botijo.

Un botijo se encuentra expuesto al sol en un da muy caluroso (39 C). el
objetivo es simular como vara la temperatura y la masa de agua contenida en
el botijo con el tiempo.
Para ello haga las siguientes suposiciones:
El botijo es perfectamente esfrico con un radio de 0.1 m.
La densidad del agua es constante 1000 kg/m3
No hay prdidas de agua por goteo o exudacin
El coeficiente de transferencia de materia en la superficie exterior y en la
interior es el mismo.
La superficie del lquido en la interfase con el aire esta a Temperatura
constante y en equilibrio con el aire
La pared seca alrededor del lquido se mantiene constante a 39C e
irradia a la superficie interior del lquido a (S 24.2 C )
Clculo numricoElcoeficienteenIngenieraglobal de transmisin de
calor U incluye todas las resistencias de conveccin en el lquido y conduccin en el
lquido y a travs de las paredes.
El lquido est perfectamente mezclado.
La capacidad calorfica de la arcilla es mucho menor en comparacin
con la capacidad calorfica del agua, de tal forma que el contenido de
calor del botijo se puede despreciar.
El factor de forma de radiacin f es constante.
La renovacin de aire en el exterior es constante.
La humedad y la temperatura del aire permanecen constantes.
Las ecuaciones que definen el sistema son entonces:
dV = k ' a H H
s
dt
d L 4 4
VC
p = hC a G S + f 273 + G 273 + s 4 R 2 S
dt
dV
Ua L S W
dt
La primera ecuacin expresa la velocidad de evaporacin como una funcin del
coeficiente de transferencia de materia del agua (k) de la superficie total (a = A
+ S) y la diferencia entre la humedad de saturacin del aire (Hs) y la humedad
ambiental (H). donde:
H s = 0.018 kg agua / kg aire seco
H = 0.011 kg agua / kg aire seco
La segunda ecuacin es un anlisis global que expresa la variacin trmica del
lquido y corresponde a un balance entre la transferencia de calor del aire al
agua (conveccin trmica sobre la superficie A+S y radiacin trmica del aire
en la cmara sin agua). Las prdidas de calor en la interfase lquido vapor y la
evaporacin del lquido.
Datos y Notas:
M 0 = 3.161kg Masa inicial de agua
Coeficiente individual de
k ' = 80 kg agua / ( h m2 ) transferencia de materia
a=A+S rea total
H s = 0.018 kg agua / kg aire seco Humedad de saturacin
H = 0.011 kg agua / kg aire seco Humedad ambiental
Cp = 1 kcal / ( kg K) Capacidad calorfica del agua
G = 39 C Temperatura del Aire
Temperatura del vapor justo encima
de la superficie del agua (coincide
S = 24.2 C
con la temperatura de saturacin
adiabtica)
Temperatura del agua (vara con
L
tiempo)
W = 583 kcal / kg Calor de vaporizacin del agua
f = 3.0510 8 kcal / ( h m 2 K4 )
2 Coeficiente global de transmisin de
U = 22 kcal / ( h m k)
calor
2
hc 3
De acuerdo con Sherwood Pigford : k ' = s Le
s = 0.24 kcal / ( kg K) Calor hmedo
Le = 1.15 Nmero de Lewis
G
Superficie interior del Agua A
h
G
S r
R 0.10 m
Superficie exterior mojada S

Hay que tener en cuenta que el volumen ocupado por el agua es

V= R3
4
3Rh 2 h3
3 3
El rea interior del agua:
A = 2 R h h2
Y la superficie exterior mojada
S = 2 R 2R h
Solucin: Se deben integrar las ecuaciones diferenciales anteriores con el
clculo simultneo de V, A, S y h. Quizs lo ms sencillo sea escribir las
ecuaciones en funcin de h.
3.5 40
3 38
36
2.5
Temperatura del agua (C)
masa de agua (kg)
34
2
32
1.5
30
1
28
0.5 26
00 10 20 30 40 50 60 70 240 10 20 30 40 50 60 70
tiempo (h) tiempo (h)
EDO Clculo.10numrico.ReactorenIngenieranoisotermo
En un reactor qumico tiene lugar la reaccin A R, que tiene lugar en fase
lquida. Se conocen los siguientes datos sobre la capacidad calorfica de la
mezcla, la cintica y el calor de reaccin:
C ps = 250cal / molK
Hr = 15000cal / mol
E / R =12000 K
k0 = 81013 min1
Las condiciones de operacin suponen un tiempo espacial =100 min y una
temperatura de entrada de la alimentacin T0 = 29.5 C
Un balance de materia en estado no estacionario, indica que para el arranque

del reactor se cumple la siguiente ecuacin:
dx A E/R xA
= k 0 exp 1 xA
dt T
Asimismo, un balance de energa en rgimen no estacionario proporciona:
dT H r k E/R T T
= 0 exp 1 xA
0
C
dt ps T
En las ecuaciones anteriores, xA representa la conversin de A y T la
temperatura.
Represente las variaciones de la conversin y la temperatura si la temperatura
inicial es de 300 K.
Solucin:
1 370
0.9 365
0.8 360
0.7 355
Temperatura
0.6 350
conversion
0.5 345
(K)
0.4 340
0.3 335
0.2 330
0.1 325
0 0 20 40 60 80 100 120 140 160 180 200 3200 20 40 60 80 100 120 140 160 180 200
tiempo (s) tiempo (s)
Ecuaciones diferenciales enderivadas parciales

EDP.1. Per fi ld etem p er atur a s e nl a par ed de un h or no
Calcule como vara con el tiempo el perfil de temperaturas de la pared de un

horno cuyo interior se mantiene temperatura constante de 1000 C. Para los
siguientes casos.
El lado exterior se mantiene 25C.
El lado exterior se puede considerar una superficie perfectamente
aislante.
El lado exterior pierde calor al ambiente.
Datos y Notas.
L = 0.1m Grosor de la pared
Cp = 400 J/(kg K) Capacidad calorfica
3
= 2000 kg/m Densidad de la pared
To = 25 C Temperatura ambiente
Tint =1000 Temperatura de la pared interior
h = 20 W/(m2 K) Coeficiente individual de transmisin
de calor
k = 0.55 W/(m K) Conductividad calorfica
Para los puntos interiores se debe cumplir que:
2
T = T = k
t x2 Cp
Se deben cumplir las siguientes condiciones iniciales y de contorno
En todos los casos
T (0, x ) = T0 x L; T (0, L ) = Tint
Para el primer caso:
T (t ,0) = T0 ; T (t , L ) = Tint
Para el segundo caso
q = k T =0 T (0, L ) = T
int
x x 0
Para el tercer caso
q = k T = h (T T ) T (0, L ) = T
x0 0 int
x x 0
Solucin: La solucin en estos casos supone la representacin con el tiempo
de todos los perfiles de temperatura, se darn slo algunos valores.
1000 1000
900 900
800 800
t = 1000 s t = 2000 s
700 700
C
Temperatura C
600 600
500 500
Clculo400 numrico en Ingeniera

m
400
T
e
p
e
a
r
r
t
300 300
200 200
100 100
0 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
x (m) x (m)
1000
900
t = 4000 s
800
700
Temperatura C
600
500
400
300
Superficie aislante
200
Perdidas al exterior
100 Paredes a T constante
0 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
x (m)
EDP.2. Difusin ms reaccin de un gas

Sobre una pelcula lquida se absorbe un gas (A) que reacciona para dar un
compuesto B. dicho compuesto a su vez reacciona para dar otro C.
k
1
A+L B rA = k1 CA
k2 r =kC k C
B C B 1 A 2 B
Un balance en estado no estacionario viene dado por
2
CA CA
= DA + rA
2
t x
2
C B = D CB + r
2
t B x B
Con condiciones iniciales y de contorno.
t = 0 CA = 0 x 0; CB = 0
x = 0 CA = CAo ; CB = 0 t
dC A d CB
x=L = 0; =0
dx dx
Utilizando un mtodo de diferencias finitas realice un programa en Matlab que
permita visualizar la variacin de la concentracin de A y B con el tiempo para
los 25 primeros segundos. (Utilice 50 puntos para discretizar el espacio x).
Indique claramente el resultado de aplicar el criterio de estabilidad.
Para el compuesto B realice un grafico en 3 dimensiones que represente la
concentracin frente al tiempo y la posicin (considere los 25 primeros
segundos).
Datos:
L = 0.5 cm;
DA = 0.001 cm2 / s
DB = 0.0001 cm2 / s
k1 = 0.05 s-1
k2 = 0.04 s-1
CAo = 1 mol / L
CBo = 0 mol / L
1 1
CA CA
0.9 t = 10 s CB 0.9 t = 40 s CB
0.8 0.8
0.7 0.7
0.6 0.6
Concentracin
Concentracin
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0 5 10 15 20 25 30 35 40 45 50 0 0 5 10 15 20 25 30 35 40 45 50
posicin (cm) posicin (cm)
1 1
CA CA
0.9
t = 80 s CB
0.9
t = 200 s CB
0.8 0.8
0.7 0.7
0.6 0.6
Concentracin
Concentracin
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0 5 10 15 20 25 30 35 40 45 50 0 0 5 10 15 20 25 30 35 40 45 50
posicin (cm) posicin (cm)
EDP.3. Flujo ms dispersin en una tubera

A lo largo de un tubo est circulando una corriente de agua. En un momento
dado se comienza a introducir en uno de sus extremos una cierta sustancia, de
tal manera que, manteniendo el caudal constante, la concentracin vara con el
tiempo adimensional de acuerdo a la siguiente expresin:
2
(t 0.4)
Co = exp
0.005
donde t hace referencia a unidades de tiempo adimensional y Co la
concentracin que entra a la tubera en un instante dado. La variacin de la
concentracin en funcin del tiempo a lo largo de la direccin axial dentro del
tubo se puede calcular mediante la expresin
2
C=M CC
x
t 2x
donde C es la concentracin, t es tiempo adimensional, M es el modulo de

difusin y x es una longitud adimensional.
C
Con condicin de contorno en el punto de salida = 0 en x=1
x
Utilizando el mtodo de Crank-Nicholson de diferencias finitas calcule el perfil
de concentraciones que tendra el tubo al cabo de una unidad de tiempo
adimensional para los siguientes tres valores de M (0.01; 0.001; 0.0001).
represente en una misma figura con distintos colores los perfiles obtenidos.
Notas:
Utilice 200 nodos para la posicin
Utilice un incremento de tiempo adimensional de 0.01 unidades.
Para las derivadas temporales utilice una diferencia de primer orden
hacia delante
Para las derivadas respecto a la posicin utilice un esquema de
diferencias centrales de segundo orden, excepto para la condicin de
contorno, donde tiene libertad para elegir el esquema que considere
oportuno.
1 1
M = 0.01 M = 0.01
0.9 t=0.6 M = 0.001 0.9 t=0.8 M = 0.001
M = 0.0001 M = 0.0001
0.8 0.8
0.7 0.7
0.6
concentracin
0.6
concentracin
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
posicin posicin
1 1
M = 0.01 M = 0.01
0.9 t=1.0 M = 0.001
0.9 t=1.2 M = 0.001
M = 0.0001 M = 0.0001
0.8 0.8
0.7 0.7
0.6 0.6
concentracin
concentracin
0.5 0.5
0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
posicin posicin
En este caso es interesante tambin ver cual es el perfil a la salida del tubo con
el tiempo (remarcar que las graficas anteriores estn todas referidas a la
posicin):
1
M = 0.01
0.9 M = 0.001
M = 0.0001
Concentracin a la salida
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0.8 1 1.2 1.4 1.6 1.8 2
tiempo
EDP.4. Enfriamiento de una esfera

Una esfera de dimetro D = 0.1 m. con una temperatura de 90 C. inicialmente
igual en todos sus puntos, se introduce en un bao de agua a 20 C.
Datos y Nomenclatura
R = radio de la esfera
k = 10 W/(m C) Conductividad trmica
Cp = 410 J/(kg K) Capacidad calorfica
h = 220 W/(m2 K) Coeficiente individual con el exterior
r = 2000 kg/m3 Densidad de la esfera

Ta = 20 C Temperatura del agua
To = 90 Temperatura inicial de la esfera.
a) Haciendo un anlisis global, es decir, suponiendo que todos los puntos

de la esfera estn siempre a la misma temperatura, si se plantea un
balance de energa se obtiene que:
Entrada + Generacin = Salida + Acumulacin
2 4 R3
0+ 0 = h4R T Ta t+ Cp Tt + t Tt
3
Dividiendo por t y tomando lmites se llega a:

Clculo numdTricoh en4IngenieraRTTa
3hT Ta
2
dt = 4 R 3 = R Cp
Cp
3
a.1) Calcule cmo vara la temperatura de la esfera durante los primeros
100 segundos y represente T vs t.
a.2) Calcule cul sera la temperatura de la esfera a los 40 segundos.
Admita en todos los casos que la temperatura del agua no vara
b) Haciendo un anlisis local o sea, si no se admite que la temperatura de

todos los puntos de la esfera varan a la misma velocidad, sino que la
temperatura depende del radio (r ), siguiendo un procedimiento anlogo
al anterior se llega a la siguiente expresin:
T 1 2 T
= 2 kr
t r Cp r r
Con condiciones de contorno:
T = hT Ta
k
r r R
k T =0
r r 0
b.1) Calcule utilizando un mtodo en diferencias finitas explcito, el perfil

de temperaturas a los 40 segundos. (Divida el radio en 30 puntos). Y
deduzca las condiciones de estabilidad.
b.2) Represente el perfil obtenido en el apartado anterior junto con el

perfil a 40 segundos del apartado a (nota: para el apartado a ser una
lnea horizontal porque en el caso a todos los puntos de la esfera estn a
la misma temperatura).
b.3) Calcule la temperatura media de la esfera a los 40 s, comente la
validez de la aproximacin hecha en el apartado a).
_
0
R
T(r)dr
Nota : T =
0R dr
Admita en todos los casos que la temperatura del agua no vara
Solucin:
Apartado a) el anlisis global lleva una ecuacin diferencial ordinaria, cuyo

resultado se muestra grficamente continuacin.
90
80
70
Temperatura (C)
60
50
40
30
200 20 40 60 80 100 120 140 160 180 200

tiempo (s)
Apartado b1) La temperatura vara tanto con la posicin como con el tiempo. El
perfil de temperaturas a diferentes tiempos es el siguiente.
100
90
t= 20 s
80
70 t= 40 s
Temperatura (C)
60
t= 60 s
50
40 t= 80 s
30
20
10
0 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05
Radio m
Apartado b.2)
La temperatura a los 40 s resultante del anlisis global es de 56.79 C. Si
integramos para conseguir la temperatura media a los 40 s en el anlisis local
la temperatura que se obtiene es de 66.47 C. Se puede observar una
desviacin importante.
Una modificacin interesante es calcular la temperatura media ponderada en
funcin de la masa (o el volumen a densidad constante). En este caso la
temperatura media es de 59.55 C
T_ 0
R
V ( r )T ( r )dr
0RV ( r )dr
Los resultados grficamente son los siguientes (las lneas rectas indican
temperaturas medias).
C

(
)
75
perfil de T a los 40 s
T media a los 40 s analisis local
T media a los 40 s analisis global
70
T media a los 40 s ponderada volumen
m
T
e
p
e
a
r
r
t
65
60
55
500 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05
Radio m
EDP.5. Simulacin de una columna cromatogrfica

En una columna de cromatografa de un metro de longitud, se introduce por
uno de sus extremos una mezcla de dos sustancias A y B. La concentracin en
el extremo de entrada para cada una de las sustancias viene dada por la
siguiente expresin:
C j
2
1 ( t 1.5)
= exp donde vA = 0.2 y vB = 0.4
x0
v
j 0.2
Donde t = tiempo (min); Cj = concentracin (mol/m3)
La ecuacin que describe el avance a travs de la columna para cada
componente j es la siguiente:
2
Cj Cj Cj
= Dj uj j = A, B
2
t x x
Donde t = tiempo (min)
Dj = difusividad efectiva
(m2/min) vj = velocidad (m/min)
Inicialmente en la columna la concentracin de cada uno de los componentes

es cero.
Como condicin contorno en el extremo de salida use:
C C
j = u j
j
t x x L
DATOS:
Componente A Componente
B
u (m/min) 0.2 0.4
D (m2/min) 210-3 10-4

Longitud Columna = 1 m.
Utilizando un mtodo en diferencias finitas totalmente implcito, desarrolle
en Matlab un programa que permita:
a) Visualizar simultneamente la variacin de concentracin de cada

uno de los componentes A y B a lo largo de la columna con el tiempo.
b) Representar, en otra figura, la variacin con el tiempo de la

concentracin de A y B (en una misma grfica) en el punto de salida
de la columna.
c) Para el componente A: Calcule el nmero de moles totales por

unidad de rea que entra a la columna; el nmero de moles totales
por unidad de rea que sale de la columna y el nmero total de moles
por unidad de rea transversal que hay dentro de la columna a los 4
minutos. Estos valores deben ser iguales, y son una medida de cmo
se cumple el balance de materia. Qu porcentaje de error est
cometiendo?. Para la integracin utilice el mtodo de los trapecios.
NOTA: Para calcular los moles por unidad de rea de un componente

que han atravesado una determinada seccin en un intervalo de tiempo
dado debemos resolver:
S t1 Q j C j dt S t 1 S u j C j dt t 1 u j C j dt
1 t2 1 t 2 t 2
Para calcular el nmero de moles de un componente dentro de la

columna en un momento dado basta resolver:
L
0 C j dx
Utilice un esquema en diferencias finitas con 80 puntos en el espacio y
utilice un incremento de tiempo de 0.01 min. Integre durante los 10
primeros minutos.
Solucin: Los perfiles a diferentes tiempos son los siguientes:

5 5
Componente A Componente A
4.5 t = 2 min Componente B 4.5 t = 3 min Componente B
4 4
Concen
Concentracio
tracion
3.5 3.5
3 3
n
2.5 2.5
2 2
1.5 1.5
1 1
0.5 0.5
0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
z/L z/L
5
Componente A
4.5 t = 3.5 min Componente B
4
Concent
racion
3.5
2.5
1.5
0.5
00 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

z/L
El perfil de concentraciones en funcin del tiempo a la salida de la columna

cromatogrfica es el siguiente.
2.5
Componente A
Componente B
2
Concentracion
1.5
0.5
00 1 2 3 4 5 6 7 8 9 10
tiempo (min)
Las integrales numricas de los moles iniciales, finales y a los 4 minutos son:
0.79266; 0.79258; 0.79248. Lo que lleva a errores del 0.011%, 0.0009%
0.012% respectivamente. Se comprueba que el error es suficientemente
pequeo, al menos en lo que a balance de materia se refiere

9.problemas Aplicados A La Ingenieria Quimica Pag 39

Cargado por

Copyright:

Formatos disponibles

9.problemas Aplicados A La Ingenieria Quimica Pag 39

Cargado por

Información del documento

Descripción original:

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

9.problemas Aplicados A La Ingenieria Quimica Pag 39

Cargado por

Copyright:

Formatos disponibles

Tema 9: Problemas aplicados a la Ingeniera Qumica

RegresionesRLM.1. Transmisinlinealesmltiplesdecalor. en una tubera

En azul aparecen los errores de la correlacin de Sieder y Tate, y en rojo los

RLM.2. Transferencia de calor en un reactor de lecho fluidizado

Determine los parmetros a1-a5 y los errores cometidos.

RLM.3. Correlacin de los parmetros de la ecuacin de Antoine para el

RLM.4. Capacidad calorfica del propano.

Polinomio de grado 3 Intervalo de confianza del 95%

Los resultados del ajuste se pueden ver a continuacin.

error polin. grado 2

200 0.5 1 1.5 2 2.5 T* -40 0.5 1 1.5 2 2.5

RLM.5. Correlacin de los coeficientes de actividad de una mezcla binaria

2 = exp x12 (2 A B ) + 2 x13( B A)

Donde x1, x2 son las fracciones molares de los componentes 1 y 2

RLM.6. Regresin de velocidades de reaccin cataltica heterognea.

ENLEcuacione.1.s noDimetrolinealesde una tubera a travs del flujo

nx j y j = 0 donde xj, yj son funciones de la fraccin evaporada.

Haciendo un balance de material.

Como y j = K j xj , sustituyendo en la ecuacin anterior se tiene que:

ENL.5. Circulacin de fluidos en una red de tuberas

P13 P23 P12 = 0

Tramo Q (l/s) P( kPa) f.103

ENL.6. Velocidad terminal

ENL.7. Grfica de Moody

ENL.8. Temperatura de llama

El oxgeno necesario para la combustin completa ser:

Si resolvemos utilizando un mtodo de Newton la ecuacin anterior para

EcuacionesLIN.1.Balanceslinealesde materia en estado estacionario ara una

15% Xileno 18 % Xileno

a) Calcula el flujo molar de las corrientes D1, D2, B1 y B2.

R=36.36% peso cromato potasico

EDO.2. Reactor de lecho fluidizado

mf es la porosidad el punto de mnima fluidizacin (igual a 0.45 en nuestro

ub = 1.6 u o umf + 113.d b0 . 5 d t1. 35 + ubr , m / s

u0 es la velocidad de alimentacin del gas (igual a 0.4 m/s en nuestro

00 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

10-1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

EDO.3. Reactor encamisado exotrmico

El siguiente paso es resolver el sistema de ecuaciones diferenciales ordinarias,

siempre entre 0 y 10), lo que nos permite utilizar un mtodo como el de la

EDO.5. Cintica de descomposicin formacin

El valor de la concentracin de A a tiempo cero es de : 22.00 mol/L.

Por ltimo, en el tercer caso se tiene un sistema de ecuaciones diferenciales

EDO.6. Reactor cataltico de membrana

La velocidad de transferencia hasta la membrana depende no slo de la

Si el paso limitante es el flujo a travs de la membrana entonces el coeficiente

EDO.7. Asociacin de RCTA isotermos

Solucin: La primera parte consiste en resolver un sistema de ecuaciones no

En el tercer caso debemos modificar las ecuaciones referentes a CB2 y CB3 de

EDO.8. Tiempo de vaciado de un depsito

atmsfera a travs del tubo inferior. Calcule y

Variacin de la altura (H) del agua con el

Tenga en cuenta que la velocidad en cada instante

donde: V2 = velocidad en la seccin de salida (2). g =

el factor de friccin se puede calcular a partir de la ecuacin de Coolebrook y

Adems, utilizando la ecuacin de continuidad se tiene que la variacin de

0.0256 0 1000 2000 3000 4000 5000 6000 7000 8000

EDO.9. Simulacin del enfriamiento del agua dentro de un botijo.