UNIVERSIDAD NACIONAL DE SAN CRISTÒBAL DE

HUAMANGA
FACULTAD DE INGENIERÌA QUÌMICA Y METALURGIA
ESCUELA DE FORMACIÒN PROFESIONAL DE
INGENIERÌA QUÌMICA

TRABAJO Nº 2
DISEÑO DE UN CONDENSADOR

DOCENTE : Ing. MENDOZA ROJAS, Cipriano

ALUMNOS : Ronald Yamil BONZANO HUAMANÍ

CURSO : DISEÑO Y SELECCIÓN DE EQUIPOS

SERIE : 500

AYACUCHO –PERÙ
2013
TRBAJO N°2 Diseño de un condensador de metanol

Diseñar un condensador horizontal de metanol para producir un destilado de 68 kg/h de metanol de

concentración XD = 0.98 molar y el resto agua. El metanol sale por el tope de la columna de
destilación como vapor saturado e ingresa al condensador a la presión de 570 mmHg( ciudad de
Arequipa). La relación de reflujo que opera la columna de destilación es 3.0. El agua de refrigeración
ingresa al condensador a 18°C en verano y sale a 30°C. Para fines de esquematizar el cálculo se
representa el siguiente esquema.

Se desea determinar:

a) Flujos de reflujo, y vapor a partir de balance de materia alrededor del condensador.

b) La carga térmica del condensador.
c) El número de tubos y el diámetro de la carcasa.
d) La caída de presión para el lado del tubo.
Resolución

Suponiendo que el metanol entra por la columna a 58°C

Para determinar la presión de saturación, determinamos con las constantes de Antoine

Metanol(1) Agua(2)
A0 B0 C0 A0 B0 C0
7.1973 1574.99 238.87 7.07406 1657.459 227.02
M = 32.042 kg/kmol M = 18.015 kg/mol
P = 570 mmHg = 76 kPa T 1= 58°C = 331.15 k

Parámetros de NRTL para metanol (1) – Agua (2)

b12 = -253.88
b21 = 845.21
ᾳ = 0.3
Para el caso de sistemas binarios, la expresión generalizada de NRTL se resume las
siguientes ecuaciones:

Lnϒ1 = x22 [ 𝜏21.G212 + 𝜏12.G12 ]

2
(x1 + x2.G21) (x2 + x1.G12)2

Lnϒ2 = x12 [ 𝜏12.G122 + 𝜏21.G21 ]

(x2 + x1.G12)2 (x1 + x2.G21)2

𝜏12 = b12 𝜏21 = b21

RT RT
G12 = Exp(- ᾳ 12𝜏12 ) G21 = Exp (- ᾳ 21𝜏21 )

R(cal/mol.K) = 1.98721
 Remplazamos los datos y obtenemos como vemos en la tabla

𝜏12 = -0.385797
𝜏21 = 1,284387
G12 = 1,122703
G21 = 0,680236

ϒ1 =Exp ( x22 [ 𝜏21.G212 + 𝜏12.G12 ] )

2
(x1 + x2.G21) (x2 + x1.G12)2

ϒ2 = Exp ( x12 [ 𝜏12.G122 + 𝜏21.G21 ])

2
(x2 + x1.G12) (x1 + x2.G21)2

Donde las únicas incógnitas seria x1 y x2 y esos valores los tendremos dando valores para
x1 y así determinaríamos los valores para los coeficientes de actividad

Por la ley de Raoult modificada  yiP = xi.ϒi.Pisat

Como la presión total es constante vamos a calcular la temperatura para cada iteración de xi

. y1 Pt = x1ϒ1P1sat ………………………….. (1)

. y2 Pt = x2ϒ2P2sat ………………………….. (2)

Sumamos las ecuaciones (1) y (2) y obtenemos:

Pt = x1ϒ1P1sat + x2ϒ2P2sat

76 = x1ϒ1 .10(A1 – B1/(C1+T)) + x2ϒ2.10( A2 - B2/( C2+ T)) ……………………. (4)

En la ecuación (4) x1 se obtiene dando valores y de acuerdo a x1 se obtiene los
coeficientes de fugacidad y nuestra única incógnita seria T, pero se obtendría
resolviendo la ecuación (4) en la que para un valor de T debe cumplir la condición. Con la
ayuda de Excel se encuentra con el método de Newton Raphson para determinar la
temperatura de equilibrio.

Se muestra la tabla y así poder graficar la T vs x,y :

X1 X2 ϒ1 ϒ2 T* P1 satura P2 satura Y1
0 1 2.58499032 1 92.136 274.8579516 75.99821829 0
0.05 0.95 2.22548509 1.00378904 83.927 208.0192265 55.42445741 0.30456821
 0.1 0.9 1.95400325 1.0143767 79.11 175.4646933 45.69066557 0.45112971
0.15 0.85 1.74577702 1.03078169 75.906 156.2333276 40.0462644 0.53831951
0.2 0.8 1.58400121 1.05220977 73.578 143.3815298 36.32349693 0.59767504
0.25 0.75 1.45701143 1.07799403 71.766 133.9951893 33.63167878 0.64221225
0.3 0.7 1.35654074 1.10755204 70.274 126.6540857 31.54342176 0.678203
0.35 0.65 1.27661554 1.1403547 68.989 120.6018761 29.83366453 0.70903658
0.4 0.6 1.21284036 1.17590333 67.84 115.3954425 28.37182289 0.73661184
0.45 0.55 1.1619247 1.21371274 66.784 110.7764461 27.08216851 0.76212171
0.5 0.5 1.12136391 1.25329891 65.793 106.5826125 25.91735576 0.78630194
0.55 0.45 1.08922012 1.29417015 64.848 102.7073264 24.84639703 0.80959194
0.6 0.4 1.06396922 1.33582121 63.938 99.08715043 23.85077093 0.83230798
0.65 0.35 1.04439245 1.37772979 63.054 95.67281666 22.91617722 0.85457867
0.7 0.3 1.02949849 1.41935524 62.192 92.43845048 22.03493472 0.876522
0.75 0.25 1.01846665 1.46013923 61.349 89.36408507 21.20109395 0.89816783
0.8 0.2 1.0106051 1.49950807 60.522 86.43144892 20.4092618 0.9194533
0.85 0.15 1.00531985 1.5368768 59.711 83.63402067 19.65728521 0.94035658
0.9 0.1 1.00209146 1.57165473 58.916 80.965596 18.94312688 0.96080841
0.95 0.05 1.00045762 1.60325228 58.138 78.42349713 18.26571646 0.98074232
1 0 1 1.63108918 57.376 75.99871955 17.62232159 0.99998315

Grafica de Tvs x,y Para Metanol(1) y Agua(2) a P = 570 mmHg = 76 kPa

 95

90

85

80

75
T

T Vs x

70 T vs Y

65

60

TV
TL

55

50
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x,y
xD
Calculo de la energía neta total retirada por el
condensador 100
T Vs x
To = 90°C = 363.15 K 90

xD = 0.98 --------> Tv = 58.3°C = 331.45 K 80

TL = 57.7 °C = 330.85 K 70

60 TV
Las lecturas de Tv y TL se realizan sobre el diagrama de TL
temperatura 50
0 0.5 1
HL = [ xDCpAMA + (1-xD)CpBMB] (TL- To) + ∆Hsol xD

TmL = (TL + To)/2 = (57.7 + 90)/2 = 73.85 °C = 347 K a esta temperatura hallamos Cp

Calculo de CpL
CP
A B C CP = A +B*T+ C*T^2 (J/mol.K) (Kcal/kg.°C)
Metanol 74.86274 -0.102315 0.000407 88.365898 0.658704036
Agua 103.8 -0.1748 0.0002748 76.2327932 1.010444099

Donde : -----------> HL = - 679.8579871 kcal/kmol

Ahora determinamos Hv

Hv = xDMA [CpA(Tv-To) + λA] + (1-xD)MB[CpB (Tv- To) + λB ]

TmV = (TV + To)/2 = (58.3 + 90)/2 = 74.15 °C = 347.3 K a esta temperatura hallamos Cp

Cp
CP = A +B*T+ C*T^2 CP
(J/mol.K) (Kcal/kg.°C) M(Kg/mol)
88.41997753 0.65910716 32.042
76.23759129 1.0105077 18.02

Determinamos el calor latente de vaporización a la temperatura de Tv = 58.3 °C= 331.45K

. λ = λb [ ( Tc – Tv ) ]0.38
(Tc – Tb)
λb normal(kJ/mol) Tb (K) Tc (K) λ(kJ/mol) λ(kcal/Kg)
Metanol 35.254 337.6 513 35.71870727 266.2572009
Agua 40.669 373.15 647.3 42.91713216 568.8544405
Dónde : HV =7942.75755 kcal/kmol

Ahora determinamos la carga de energía neta total retirada por el condensador

Qc = DT ( 1+R).(Hv-HL)

Donde el destilado total es: DT = 68 kg x 1 kmol = 2.1222 kmol

h 32.042 kg h

Qc = (2.122 kmol).(1+3).( 7942.75755 +-679.8579871 ) kcal/kmol

h
Qc = 73196.1621 kcal = 73196.1621kcal x 1h x 1 J = 85.126 kW
h h 3600s 0.23885cal

Ahora construimos nuestra verdadera temperatura

T1 = 90°C

T2 = 57.7 °C ∆T1 = T1 – t2 = 60°C

t2 = 30 ∆T2 = T2 – t1 = 39.7°C

t1 = 18°C

Hallamos la temperatura media logarítmica

∆TmL = (∆T1 - ∆T2)/ln(∆T1/∆T2) = 49.15 °C

Hallamos la capacidad calorífica del agua a Tm = ( 30 + 18 )/2 = 24 °C

CP H2O = 4.23 kJ/kg°C

Hallamos el flujo del agua refrigerante

Flujo de agua refrigerante = 85.126/(Cpx∆t) = /(4.23x(30-18)) = 1.68 kg/s

Usamos un paso de la carcasa y un pasos de tubos

R = (T1 – T2) /(t2-t1) = (90-57.7)/(30-18) = 2.69

S = (t2-t1)/(T1-t1) = (30-18)/(90-18) =0.17

De la figura 12.9 --------------> Ft = 0.98

Y luego hallamos el fator de corrección de la temperatura

∆Tm = Ftx∆TmL = 0.975x49.15 = 47.92°C

De la tabla 12.1 ------------------> U = 250 W/m2.°C

Área de provisional------------> A = QCH3OH /(U.∆Tm) = 85.126x103/(250x47.92) = 7.106m2

De acuerdo a las dimensiones de los tubos estándares escogemos de acuerdo a nuestros

problemas o los más usados:

20 mm de OD , 16 mm ID y la longitud del tubo de: L = 1.85 m---------> L = 1.85 m

Hallamos el área de un tubo

At = L xDox∏= 1.85x20x10-3x3.1416 = 0.12 m2

Hallamos los números de tubos

Nt = A/At = 7.106/0.12 = 59.22≡ 60

Tomamos el fluido del lado de la carcasa usamos 1.25 triangular pitch

Db = do(Nt/k1)1/n1 = 20(60/0.319)1/2.142 = 230.6mm

De la figura 12.10 usamos Split-ring floating head ----> 50 mm

Ds = 230.6+50 = 280.6 mm

Hallamos los coeficientes al lado del tubo

Temperatura media del agua-------> = (30 + 18) / 2 = 24°C

𝜋
Área de la sección atravesada del tubo = x162 =201 mm2
4

Flujo del área del tubo = 60x201x10-6 = 0.01206m2

Velocidad de la masa de agua por área = 1.68/(0.01206) = 139.3 kg/s.m2

Densidad del agua = 998.204 kg/m3

viscosidad del agua = 1.005 mNs/m2

Velocidad lineal del agua = 139.3/998.204 = 0.1396 m/s

Re = d.v.di/u = 998.204x0.139x16x10-3/0.8x10-3 = 2775.04

hi = 4200(1.35+ 0.02xtm)v0.8 = 4200x(1.35+0.02x24)0.13960.8 = 913.66 W/m2°C

di0.2 160.2
Hallando los coeficientes al lado de la carcasa

Seleccionamos los espacios de los bafles = Ds/5 = 280.6/5 =56.12 mm

Tubo pitch = 1.25x20 = 25 mm
Calculo del área de cruce --> As = (Pt-do).DsLb/Pt = (25-20)x280.6xx10-6/25 = 0.00315 m2
Velocidad de la masa------> Gs = (272/3600)x1/As = 23.99kg/sm2
Diámetro equivalente: de= 1.10 (Pt2 – 0.917.do2)/do = 1.10x(252-0.917x202)/20 = 14.2 mm
Temperatura media al lado de la caracasa-- > = (90+57.7)/2 = 73.85 °C = 347.5 K
Densidad del metanol = 719.74 kg/m3
Viscosidad ----> u = 0.2873 cP = 0.2873 mNs/m2
Capacidad calorífica = 2.758 kJ/kg°C
Conductividad térmica ----------> Kf = 0.19 W/m°C
Re = Gs de/u = 23.99 x 14.2x10-3 / 0.2873 x 10-3 = 1186.14

Pr = Cp.u/Kf = 2.758 x 0.2873/0.19 = 4.17

De la figura 12.29 obtenemos jn= 1.8x 10-2

hs = kf x jn x Re x Pr1/3 = 0.19 x 1.8x 10-2 x 1186.14x 4.171/3 = 459.82 w/m2°C

de 14.2x10-3

Hallamos el coeficiente global

1 = 1 + 1 + do.Ln(do/di) + do x 1 + do x 1
Uo ho hod 2Kw di hid di hi

ho = 459.82W/m2°C
hod = 5000 de la table 12.2
Kw = 50 W/m°C de la tabla 12.6
do = 20 mm del factor estándar de diseño de los tubos
di = 16 mm del factor estándar de diseño de los tubos
hid = 3000 W /m2°C de la tabla 12.2
hi = 913.66 W/m2°C

Al reemplazar los dato obtenemos :

Uo = 238 W/m2°C como vemos hey probado con los demás valores de U supuesto y es el que
más se acerca.
Determinamos la caída de presión a lado del tubo

De la figura 12.24 para Re = 2775.04 obtenemos en el gráfico.------> jf = 5.5X 10-3

Donde :
Np = 1
jf = 5.5x10-3
L = 1.85 m
di = 16 mm = 16 x 10-3
u = uw : descuidando el factor de corrección de la viscosidad
d = 998.204 kg/m3
vt = 0.1396 m/s

Al reemplazar los datos obtenemos la caída de presión al lado del tubo

∆Pt = 73.80 Pa