PRACTICAS DE SIMULACION DE PROCESOS
PRACTICA N° 7
SIMULACIÓN DE FENOMENOS DE TRANSPORTE DE ENERGÍA
1. INTRODUCCION
Mediante Excel es factible poder desarrollar simuladores para verificar una serie
de pruebas y poder optimizar un proceso, así mismo se puede hacer con cualquier
lenguaje de programación para esto es importante, de antemano, conocer las
operaciones y los proceso que se desarrollan en un proceso real y modelar las
variables importantes que nosotros quisiéramos evaluar y luego hacer uso de
soluciones analíticas y numéricas de acuerdo a la dificultad y necesidad de
solución del problema.
2. 1RA SIMULACION
CALCULO DEL FACTOR DE COMPRESIBILIDAD
a) Utilice Excel para calcular el volumen del vapor de agua (la temperatura critica
es Tc=647.4 K y la presión critica es Pc=218.3 Atm) a Tr = 1y Pr = 1.2. Compare
los resultados obtenidos con los valores que figuran en una base de datos de
propiedades físicas (V=0.052456 L/g-mol). Complete los cálculos para Tr=3.0 y
Pr=10 (V=0.0837 L/g-mol). Realice solo los cálculos si C>0.
b) Calcule el factor de compresibilidad y el volumen molar del vapor de aguas
utilizando Excel para las temperaturas y presiones reducidas recogidas en la
tabla siguiente.
Pr Pr Pr Pr Pr Tr
0.10 2.00 4.00 6.00 8.00 1.00
0.20 2.20 4.20 6.20 8.20 1.20
0.40 2.40 4.40 6.40 8.40 1.50
0.60 2.60 4.60 6.60 8.60 2.00
0.80 2.80 4.80 6.80 8.80 3.00
1.00 3.00 5.00 7.00 9.00
1.20 3.20 5.20 7.20 9.20
1.40 3.40 5.40 7.40 9.40
1.60 3.60 5.60 7.60 9.60
1.80 3.80 5.80 7.80 9.80
10.00
Confecciones una tabla y una gráfica del factor de compresibilidad frente a Pr y
Tr, así como una tabla y un gráfico del volumen molar frente a la presión y Tr, La
presión y el volumen deberían estar en una escala logarítmica en la segunda
grafica.
La ecuación R-K generalmente se escribe (Shacham y colaboradores) como:
𝑅𝑇 𝑎
𝑃= − (4.5.1)
𝑉−𝑏 𝑉(𝑉+𝑏)√𝑇
5/2
𝑅 2𝑇𝑐
Donde: 𝑎 = 0.42747 ( ) (4.5.2)
𝑃𝑐
1
DR. ROLANDO S. BASURCO CARPIO
𝑅𝑇𝑐
𝑏 = 0.08664 ( ) (4.5.3)
𝑃𝑐
y P= presión en Atm
V= volumen molar en l/g-mol
T= temperatura en K
R= constante de los gases (R=0.08206 atm.l/g-mol.K)
Tc= temperatura critica en K
Pc= presión critica en atm
El factor de compresibilidad viene dado por
𝑃𝑉
𝑧= (4.5.4)
𝑅𝑇
La ecuación 4.5.1 también se puede escribir, de forma algebraica, en función del
factor de compresibilidad como una ecuación cubica
𝑓(𝑧) = 𝑧 3 − 𝑧 2 − 𝑞𝑧 − 𝑟 = 0 (4.5.5)
donde:
𝑟 = 𝐴2 𝐵 (4.5.6)
𝑞 = 𝐵2 + 𝐵 − 𝐴2 (4.5.7)
𝑃𝑅
𝐴2 = 0.42747 ( 5/2 ) (4.5.8)
𝑇𝑅
𝑃
𝐵 = 0.08664 ( 𝑅 ) (4.5.9)
𝑇 𝑅
En la que Pr es la presión reducida (P/Pc) y Tr, es la temperatura reducida(T/Tc).
La ecuación 4.5.5 se puede resolver analíticamente para raíces cubicas. Algunas
de estas raíces son complejas.
Considerando solo las raíces reales, la secuencia de cálculos implica lo
siguientes pasos:
𝑓 3 𝑔 2
𝐶 =( ) +( ) (4.5.10)
3 2
donde:
−3𝑞−𝑞
𝑓= (4.5.11)
3
−27𝑟−9𝑞−2
𝑔= (4.5.12)
27
Si C>0, hay una sola solución real para z dada por:
1
𝑧 =𝐷+𝐸+ (4.5.13)
3
donde
𝑔
𝐷 = (− + √𝐶)1/3 (4.5.14)
2
2
PRACTICAS DE SIMULACION DE PROCESOS
𝑔
𝐸 = (− − √𝐶)1/3 (4.5.15)
2
Si C<0, hay tres soluciones reales:
−𝑓 𝜙 2𝜋(𝑘−1) 1
𝑧𝑘 = 2√ 𝑐𝑜𝑠 [( ) + ]+ 𝑘 = 1, 2, 3 (4.5.16)
3 3 3 3
donde
𝑔2 /4
𝜙 = 𝑎 𝑐𝑜𝑠√ (4.5.17)
(−𝑓3 )/27
En la región supercrítica, cuando Tr >= 10, dos de las soluciones son negativas,
por lo que se selecciona el valor máximo zk, como el factor de compresibilidad
verdadero.
3. 2DA SIMULACION
DISTRIBUCIÓN DE AGUA
En la Figura se muestra una tubería que distribuye agua a una temperatura
constante, T= 60 °F, desde el punto 1 donde la presión es p 1 = 150 psi y la
elevación es z1 = 0 ft al punto 2 donde la presión es la presión atmosférica y la
elevación es z2 = 300 ft. La densidad y la viscosidad del agua se pueden calcular
a partir de las ecuaciones siguientes:
= 62,122 + 0,0122T- 1,54xl0-4xT2 + 2,65 x 10-7 T3 - 2,24x10-10T4
1057,51
Ln = - 11,0.318 +
𝑇+214,624
Donde T se expresa en °F, en Ibm/ft3 y en Ibm/ft • s.
a) Calcule la velocidad de flujo q (en gal/min) para una tubería con una longitud
efectiva L = 1.000 ft y hecha de un tubo de acero comercial tipo 40 de 8 pulgadas
de diámetro nominal. (Solución: u = 11,61 ft/s, gpm = 1.811 gal/min).
b) Calcule las velocidades del fluido en ft/s y en gal/min para tuberías a 60 °F con
longitudes efectivas de L= 500,1.000,... 10.000 ft y hechas de tubos de acero
comercial tipo 40 de 4, 5, 6 y 8 pulgadas de diámetro nominal. Utilice el método
de sustituciones sucesivas para resolver las ecuaciones de los diferentes casos y
presente los datos en forma tabular. Prepare las gráficas de la velocidad del fluido
frente a D y L y la velocidad de flujo q frente a D y L.
c) Repita el apartado a) a las temperaturas T = 40, 60 y 100 °F y exponga los
resultados en una tabla indicando la temperatura, densidad, viscosidad y
velocidad de flujo.
Ecuaciones y Datos
3
DR. ROLANDO S. BASURCO CARPIO
El balance de energía mecánica general de un líquido incompresible aplicado a
este caso da lugar a:
1 2 𝑔𝑐 ∆𝑃 𝑓𝐹 𝐿𝑣 2
− 𝑣 + 𝑔∆𝑧 + +2 =0
2 𝜌 2
donde es la velocidad del fluido en ft/s, g es la aceleración de la gravedad, g =
32,174 ft/s2, z = z2 – z1 es la diferencia de elevación (ft), gc es un factor de
conversión (en unidades inglesas gc = 32,174 ft • lbm/lbf • s2), P = P2 – P1 es la
diferencia en presión Ibm/ft2),fF es el factor de fricción de Fanning, L es la longitud
del tubo (ft) y D es el diámetro interior del tubo (ft). La utilización del método de
sustituciones sucesivas requiere que la Ecuación anterior sea resuelta para :
𝑔𝑐 ∆𝑃 𝑓𝐹 𝐿
𝑣 = √(𝑔∆𝑧 + ) / (0,5 − 2 )
𝜌 𝐷
La ecuación para el cálculo del factor de fricción de Fanning depende del número
de Reynold, Re = D/, donde es la viscosidad en Ibm/ft • s. Para régimen
laminar (Re < 2.100), el factor de fricción de Fanning se puede calcular a partir de
la ecuación:
16
𝑓𝐹 =
𝑅𝑒
Para régimen turbulento (Re > 2.100) se puede utilizar la ecuación de Shacham
𝜀/𝐷 5,02 𝜀/𝐷 14,5 2
𝑓𝐹 = 1/16 {𝑙𝑜𝑔 [ − 𝑙𝑜𝑔 ( + )]}
3,7 𝑅𝑒 3,7 𝑅𝑒
Donde /D es la rugosidad de la superficie del tubo ( =0,00015 ft para tubos de
acero comercial).
La velocidad del fluido en las tuberías se puede transformar en velocidad de flujo
(caudal) multiplicando aquélla por la sección transversal del tubo, por la densidad
del agua (7,481 gal/ft) y por el factor (60 s/min). De este modo q tiene las unidades
de (gal/min). Los diámetros interiores (D) de los tubos de acero comercial tipo 40
de 4, 5, 6 y 8 pulgadas, son 4,026, 5,047, 6,065 y 7,981 respectivamente.
