CAPITULO 5.

Las temperaturas de burbuja y rocío de la mezcla se pueden evaluar utilizando

modelos matemáticos que correlacionan las propiedades del sistema como
temperatura, presión y composiciones. Se calculó el coeficiente de reparto
expresado como la relación de concentraciones en ambas fases, dependiente de
la temperatura y teniendo en cuenta la presión ambiental, se optó por usar la
ecuación de Antoine que describe la relación entre la temperatura y la presión de
saturación del vapor, teniendo en cuenta que se conocen las composiciones en la
alimentación dichas temperaturas se pueden calcular a partir de la ecuación de
Rachford-Rice ya que este es la única variable en la ecuación.

B Ecuación 1
PS = A+
T +C

PS Ecuación 2
Ki=
P

Un valor K  alto, indica que el componente tiende a concentrarse en el vapor; si es

bajo, tiende a concentrarse en el líquido. Si el valor K es la unidad, el componente
se dividirá igualmente entre el vapor y el líquido.

Calculando las temperaturas de burbuja ( = 0) y rocío ( = 1) de la mezcla.

( Ki−1 ) Zi Ecuación 3
f ( φ )=
1+(1−φ)( Ki−1)

A partir de esto se obtuvo que la temperatura de burbuja y rocío es 78,16°C y

82,58°C respectivamente,

Con los valores del coeficiente de distribución se calculó la volatilidad relativa,

definida como la razón de los valores K del componente más volátil a uno de los
menos volátiles. Si la volatilidad relativa es alta, un componente tiene una
tendencia mucho mayor a vaporizar (es más volátil) que el otro, y será fácil
separarlos uno del otro mediante la vaporización de uno; del otro lado cuando la
tendencia a vaporizar de los componentes es similar, la volatilidad relativa se
acerca a la unidad. Para efectos de los cálculos se seleccionó el propanol como el
compuesto clave con volatilidad intermedia. A partir de dichos cálculos se
determinó que los componentes concentrados en la fase vapor son acetato de
etilo y etanol, el n-propil acetato y el metoxipropanol que son los componentes
menos volátiles se encuentran en fase líquida.

Prueba de destilación.

La prueba de destilación se realizó para obtener la respectiva curva destilación

conforme al método establecido por NTC 6191, donde se tomó la temperatura y el
volumen del condensado para que el resultado se refleje en forma de la
temperatura de ebullición en función del volumen del destilado.

La tabla X muestra los resultados obtenidos en la destilación y de acuerdo a estos

resultados, se realizaron la curva de destilación.

Tabla X. Resultado destilación de la muestra de residuo.

Prueba 1 Prueba 2
Volumen Temperatura Tiempo Volumen Temperatura Tiempo
(mL) observada (°C) (min) (mL) observada (°C) (min)
0 74 0,00 0 74 0,00
10 74 8,08 10 75 4,24
20 74 11,09 20 76 6,22
30 78 17,32 30 80 9,55
40 80 20,26 40 81 12,41
50 81 23,05 50 82 14,59
60 82 26,24 60 84 26,05
70 84 29,40 70 86 28,12
80 86 33,11 80 88 32,62
90 88 36,42 90 89 39,18

1
NTC 619- Método de ensayo para el intervalo de destilación de líquidos orgánicos volátiles.
Figura X. Curva de destilación.

90

85
Temperatura (°C)

80
Prueba 1
Prueba 2
75

70
0 10 20 30 40 50 60 70 80 90 100
Porcentaje volumen recuperado (%)

De acuerdo a los resultados mostrados en la tabla x y figura x , se establece la

temperatura inicial y final de ebullición de la mezcla es de 74 y 88,5 °C
respectivamente, durante la experimentación se encontraron variables no
controlables como la temperatura, angulo de inclinación del condensador, agua de
refrigeración entre otros, esto porque se disponia con un equipo de destilación
simple lo cual dificultaba el manejo de dichas variables, según esto se establece
que el control de la temperatura juega un papel muy importante en la operación.
Ademas al comparar la temperatura de burbuja calculada con la temperatura
inicial de destilación observada en el laboratorio se verifico que son similares.

Como residuo de la destilación se obtiene un lodo, que dadas las características

de las materias primas involucradas está compuesto por resinas sintéticas y
pigmentos. En la tabla x se muestra la relación

Tabla X. Resultado destilación residuo sólido.

Prueba 1 Prueba 2
Peso inicial residuo solvente(g) 86,674 86,674
Peso final residuo sólido(g) 8,947 9,456
Porcentaje (%) 10,323 10,910

Se realizó una prueba de destilación con un rotaevaporador para observar el

efecto de la presión sobre el sistema, para ello se tomaron 200 mL de residuo,
ajustando las condiciones de destilación según el manual del equipo de la
siguiente manera: la temperatura del baño de 60°C, temperatura de vapor de 40°C
y la temperatura de agua de refrigeración máximo 20°C. En la tabla x se muestran
los resultados obtenidos.

Tabla X. Resultado destilación al vacío.

Prueba 1 Prueba 2
Presión (mbar) 130 175
Temperatura destilación (°C) 78,4 60,0
Volumen inicial residuo (mL) 200 200
Volumen recuperado (mL) 173 182
Tiempo destilación (min) 131 190
Recuperado (%) 86 91

La presión de vacio se selecciono de acuerdo a las especificaciones del equipo,

donde se debe identificar el componente en mayor proporcion en la mezcla, es
decir, el etanol, disolvente encontrado en la tabla de datos del manual, el cual
tiene establecido las condiciones de temperatura y presión para trabajar al vacio.

RESULTADOS Y COMPARACIÓN CROMATOGRAFIA DE GASES.

Operación destilación por lotes.

La operación de una columna de destilación discontinua convencional se puede

describir convenientemente en tres partes: el período de puesta en marcha, el
período de producción, y el período de parada. Antes de la puesta en marcha, la
alimentación se carga normalmente al rehervidor y luego se calienta.
Durante el período de puesta en marcha, el sistema se hace funcionar bajo reflujo
total, es decir, sin ningún producto retirado de la parte superior, hasta que se
alcanza un estado estacionario o el producto de la composición de destilación
alcanza una pureza deseada. El período de producción es normalmente el que
consume más tiempo, pero para las separaciones difíciles tales como
separaciones de alta pureza o azeotrópicas, el tiempo de inicio también puede ser
significativo.
Durante el período de producción, se elimina continuamente una fracción del
condensado de cabeza de la columna de acuerdo con el reflujo elegido. A medida
que transcurre el tiempo, el material que se destila se vuelve menos rico en los
componentes más volátiles y la destilación de un corte se detiene cuando la
composición del destilado acumulado alcanza el valor deseado.

La destilación por lotes o discontinua utiliza dos modos de operación:

 Reflujo constante: La relación de reflujo se ajusta a un valor predefinido,

que se mantiene durante la duración del corte dado. Por lo tanto, la
 composición destilada instantánea cambia continuamente. La destilación se
continúa hasta que el destilado medio ha alcanzado el valor deseado, y
después se cambia el corte.
 Reflujo variable: Cuando se desea mantener una composición de
destilación constante, se debe aumentar constantemente el reflujo a lo
largo del corte o de la operación. A medida que avanza el tiempo, el
componente más volátil se está continuamente agotado. La relación de
reflujo aumentara progresivamente.

Aunque las operaciones para cada modo de operación pueden dar la misma
cantidad de producto con la misma pureza, el tiempo total requerido para cada
operación difiere generalmente. En el caso de una operación de reflujo constante,
la composición del producto se deja variar, pero el reflujo se mantiene constante.
Esto significa que la relación de reflujo utilizada puede ser mayor que la necesaria
al comienzo del corte, y por lo tanto el producto destilado tiene una pureza muy
por encima de la especificación. Por otra parte, hacia el final del corte, la relación
de reflujo es demasiado baja y el producto destilado está entonces por debajo de
la especificación.

Un argumento similar puede hacerse para el modo de operación de composición

de destilado constante, donde la composición es fija pero se permite que la
relación de reflujo varíe. Ambas tienen sus limitaciones. Generalmente se
encuentra que el modo de operación de composición de destilado constante
produce un mejor rendimiento que el de reflujo constante, pero es más compleja
de implementar.

El diseño de una columna de destilación por lotes es similar al diseño de una

columna de destilación continua en la que deben determinarse los siguientes
parámetros principales de diseño:
1. Número de etapas o placas en la sección de columna (y por lo tanto la altura de
la columna)
2. Diámetro de la columna
3. Tamaño del rehervidor (tamaño del lote y área de transferencia de calor)
4. Tamaño del condensador (área de transferencia de calor)
5. Tamaño del tambor de reflujo

DESTILACIÓN INSTANTANEA.

La destilación instantánea o destilación flash, es uno de los procesos de

separación más sencillos y comunes, en este proceso parte de la corriente de
alimentación que se evapora para separar vapor y liquido en equilibrio, donde el
componente más volátil estará en mayor proporción en el vapor. 2

Equilibrio vapor líquido sistema multicomponente.

El proceso de destilación involucra un equilibrio entre la fase líquido-vapor, el

modelo para el equilibrio se representa mediante la ecuación 4.

y i P=x i Pisat Ecuación 4

Dónde:

 xi: Fracción mol en fase líquida

 yi: Fracción mol en fase vapor
 Psat : Presión de vapor de las especies puras i.

Teniendo en cuanta las desviaciones de la idealidad de las soluciones en la fase

liquida, la ley de Raoult modificada se obtiene cuando i, coeficiente de actividad
se adiciona.

y i P=γ i x i P sat
i
Ecuación 5

γ i x i P sat
i
Ecuación 6
y i=
P

yi P Ecuación 7
x i= sat
γi P i

xi Ecuación 8
x ¿=
∑ xi
i

La presión de vapor se determina mediante la ecuación de Antoine, la cual

describe la relación entre la temperatura y la presión de vapor de las sustancias
puras.

2
WANKAT, Ingenieria de procesos de separación. 41 p.
 B Ecuación 9
ln Psat
i =A + + D∗LnT + E T F
C +T

B Ecuación 10
+ D∗ln T isat + E T sat F
LnP= A+ sat i
C +T i

Tabla . Constantes de Antoine.

Sustancia A B C D E F
Etanol 86,486 -7931,1 0 -10,2498 6,39E-06 2
Metoxipropanol 55,425 -6886,7 0 -5,5896 2,46E-17 6
Propanol 79,4625 -8294,91 0 -8,90961 1,82E-06 2
N-propil acetate 98,6228 -8038,41 0 -12,4522 8,86E-06 2
Acetato de etilo 88,376 -7147,86 0 -10,9917 8,55E-06 2

La determinación del coeficiente de actividad se establecio mediante el arbol de

selección de Carlson3, conociendo de antemano la composición de la mezcla, la
presión del sistema y la naturaleza de las sustancias involucradas, se definio que
el modelo termodinamico que más se ajusta a dichas condiciones y por
disponibilidad de parametros es el NRTL.

El modelo NRTL es una extensión del concepto de Wilson aplicable a sistemas

multicomponentes del tipo vapor-líquido, líquido-líquido, y vapor-líquido-líquido. Al
requerir solo constantes de pares binarios, retiene hace que el método sea simple
además de basarse en el parámetro α ij, se basa en una fracción molar de la
especie considerada. Esta se expresa:

x exp ⁡(−α ji τ ji ) Ecuación 11

x ji = NC j
∑ x k exp ⁡(−α ki τ ki )
k=1

La expresión que permite calcular los coeficientes de actividad para la fase líquida
de mezclas multicomponentes en el modelo NRTL está dada por:

3
CARLSON, Eric C. Don’t gamble with physical properties for simulations. Chemical Engineering Progress.
1996. p 36-46
 NC NC
Ecuación 12
∑ τ ji G ji x j NC
G ij x j ∑ τ kj Gkj x k
ln γ iL= j=1NC +∑ NC
∗¿(τ ¿ ¿ ji− k=1NC )¿¿
j=1
∑ Gki x k ∑ Gkj x k ∑ Gkj x k
k=1 k=1 k=1

Donde:

G ji =exp ⁡(−α ji τ ji ) Ecuación 13

Los coeficientes  vienen dados por las siguientes expresiones:

gij −g jj Ecuación 14
τ ij =
RT

g ji −gii Ecuación 15
τ ji =
RT

Donde gij, gjj son las energías libre molares de Gibbs de interacción entre pares de
moléculas. En las ecuaciones anteriores se debe considerar Gji ≠ Gij, ij ≠ ji,
Gii=Gjj=1, ii=jj=0.

Las constantes de la ecuación NRTL para las sustancias involucradas se tomaron

de la base de datos de Hysys Aspen.

Tabla. Constantes ecuación para NRTL

Etanol Propanol NPA Metoxi Acetato 

Etanol 0 3122,736 947,937 -561,891 1240,399 0,3
Propano
l -1909,146 0 298,514 -1922,062 392,428 0,3
NPA -475,594 -319,105 0 -1467,698 535,626 0,3
Metoxi 1035,64 5487,552 2449,609 0 2709,21 0,3
Acetato -716,228 -425,852 -463,923 -1642,714 0 0,3

Determinación punto de burbuja y roció.

Los cálculos de la temperatura de burbuja y rocío requieren procesos iterativos

dado que la resolución de un sistema de ecuaciones llega a ser complejo.
Los coeficientes de actividad son función de la temperatura y de la composición de
la fase líquida. Ya que ∑i yi = 1 y ∑i xi = 1 la ecuación 5 es la suma sobre todas las
especies y se obtiene:

P=∑ xi Psat
i γi
Ecuación 16
i

De manera alternativa se tiene para xi:

yi P Ecuación 17
1=∑ sat
i γi P i

La primera estimación para saber si una mezcla multicomponente producirá o no

una mezcla de dos fases en equilibrio cuando se so mete a una operación de
flash, a una temperatura y presión dadas, se puede obtener mediante los valores
K.4

yi Ecuación 18
k i=
xi

Con respecto a la ecuación 5 se indica que el valor de K para la ley de Raoult

modificada es:

γ i Psat
i
Ecuación 19
k i=
P

Conforme a la relación de equilibrio (ecuación 6), yi=Kixi. La suma con ∑i yi = 1

produce:

∑ K i x i=1 Ecuación 20
i

En estos términos, para los cálculos del punto de burbuja, donde se conoce xi,
definiendo los valores de K que satisfagan la ecuación 11.

y Ecuación 21
∑ Ki =1
i i

4
Operación de separación por etapas. 332 p
De igual manera, conforme a la relación de equilibrio yi/Ki=xi. La suma con ∑i xi =1,
para los cálculos del punto de rocío, donde se conoce y i, definiendo los valores de
K que satisfagan la ecuación 21.

La presión de vapor para cada especie PiSat son funciones débiles de la

temperatura, por lo cual se establece una nueva relación:

sat y i P sat
j
Ecuación 22
P j =P ∑ ( ¿ )¿
i γ i Pisat

La sumatoria se debe hacer sobre todas las especies incluyendo a j, la cual es

elegida del grupo de sustancias involucradas en la mezcla. La temperatura
correspondiente a la presión de saturación se determina mediante la ecuación 6
en función de T, para encontrar la T inicial para dar comienzo al proceso iterativo.

Para una estimación inicial se establece:

T =∑ x i T sat
i
Ecuación 23
i

La figura x muestra el esquema iterativo para el cálculo de la temperatura de

burbuja.

Figura Algoritmo calculo temperatura de burbuja.

En la figura x se muestra el esquema para la determinación del punto de rocio y
como en el caso anterior el modelo iterativo es controlado por la temperatura
mediante:

T =∑ y i T isat Ecuación 24
i

Figura Algoritmo calculo temperatura de rocío

Muestra de cálculo temperatura de burbuja.

Peso Molecular Fracción Fracción

Compuesto (g/mol) másica molar
Etanol 46,06844 0,47973 0,65502
Propanol 60,09502 0,34217 0,027643
NPA 102,13 0,10633 0,210734
Metoxipropanol 90,12 0,02641 0,074215
Acetato de etilo 88,1 0,04536 0,032385

A partir de la ecuación 12 se determinó el coeficiente de actividad para cada

componente.
 Mediante la ecuación 9, se calcula la presión de vapor de cada componente

Psat
etanol =exp ¿

Psat
etanol∗γ etanol
k etanol=
75000 Pa

A partir de esto y mediante el uso de solver de Excel se obtuvo que la temperatura

de burbuja 356,346K (83,196°C).

Tabla . Resultados determinación punto de burbuja.

Sustanci
a Ln   xi Ps (Pa) Ki Ki*xi yi
Etanol 0,105533 1,111304 0,47973 98534,26 1,460019 0,700415 0,700415
Propanol -2,046099 0,129238 0,02641 45631,80 0,078631 0,002077 0,002077
NPA 0,052707 1,054121 0,34217 45586,19 0,640712 0,219232 0,219232
Metoxi -0,326819 0,721214 0,10633 21663,87 0,208323 0,022151 0,022151
102729,7
Acetato -0,101668 0,903329 0,04536 8 1,237318 0,056124 0,056125

Muestra de cálculo punto de rocío.

Inicialmente se supone que el coeficiente de actividad tiene un valor igual a 1.

Mediante la ecuación 6, se determinó la temperatura de saturación para cada

compuesto:

75 kPa=exp ¿

T sat
etanol=343,91 K

Tabla . Temperatura de saturación.

Sustancia Tsat (K)

Etanol 343,911
Propanol 384,287
NPA 362,573
Metoxi propanol 365,119
Acetato de etilo 341,328
A partir de estos datos y por medio de la ecuación 15 se obtiene la temperatura de
saturación de la mezcla, con un valor de 350,45 K, con este valor se estima la
presión de saturación para cada componente.

Psat
etanol =exp ¿

Psat
etanol =97881,10 Pa

Tabla . Presión de saturación a 350,45 K

Sustancia Psat (Pa)

Etanol 97881,10
Propanol 45302,21
NPA 45311,86
Metoxi propanol 21516,33
Acetato de etilo 102167,93

Se calculan las fracciones molares de líquido para cada componente, mediante la

ecuación 5 y con estos valores se estiman la fracción molar normalizada para
cada componente.

0,65502∗75000 Pa
x etanol= =0,50189
97881,10 Pa

0,50189
x etanol= =0,425717
0,50189+0,045764 +0,348816+0,258694 +0,023774

Tabla . Fraccion molar líquida

Sustancia Xi XNi
Etanol 0,501896 0,425717
Propanol 0,045764 0,038818
NPA 0,348816 0,295871
Metoxi propanol 0,258694 0,219428
Acetato de etilo 0,023774 0,020165
Sumatoria 1,178945 1

La relación de equilibrio se determina a partir de la ecuación 10.

1∗97881,10 Pa
k etanol= =1,30508
75000 Pa
De acuerdo al algoritmo para la determinación del punto de roció, se debe estimar
el coeficiente de actividad y definir que la diferencia de i calculado con el i
supuesto (1) sea menor.

Tabla . Relación de equilibrio y coeficiente de actividad.

Sustancia ki i 
Etanol 1,3050814 1,0574936 0,057493601
Propanol 0,60402943 0,09011223 -0,909887766
NPA 0,60415812 1,14036225 0,140362250
Metoxi propanol 0,28688436 0,83598303 -0,164016973
Acetato de etilo 1,36223914 0,97167714 -0,028322857

Como se cumple con la relación se prosigue con el procedimiento. Ahora se debe

calcular la PjSat, donde se selecciono al etanol como especie j, para ello se calculo
el valor de yi para cada especie (ecuación 5).

1,0574936∗0,501896∗97881,10 Pa
y etanol = =0,69267
75000 Pa

Psat
j =115,39 kPa

Tabla . Parametros para PjSat

yi yi/i Pjsat/Pisat
0,6926749 0,6550157
4 3 1
0,0024909
9 0,0276432 2,160625537
0,2403199 0,2107399
3 9 2,160165297
0,0620426 0,0742152
7 2 4,549154984
0,0323858
0,0314686 6 0,958041335

Con el valor de PjSat se determina un nuevo valor de temperatura para cada

especie y se establece la diferencia entre las temperaturas calculadas hasta que
se cumpla T< A partir del modelo iterativo se obtuvo que la temperatura de
rocío converge cuando es igual a 361,922 K (88,773°C).
 Tabla . Resultados determinación punto de rocio.

Sustancia  Ps (Pa) Xi XNi Ki

Etanol 1,057544 151891,33 0,323429 0,421637 2,025218
Propanol 0,110005 73075,66 0,028371 0,036985 0,974342
NPA 1,138421 67487,73 0,234198 0,305311 0,899836
Metoxipropanol 0,860553 33830,05 0,164532 0,214492 0,451067
Acetato de etilo 0,982521 146777,17 0,016548 0,021573 1,957029

Balance de materia.

F=L+V Ecuación 25

z i F=x i L+ y i V Ecuación 26

La relación de equilibrio (ecuación 18), es una medida que indica la tendencia de

una especie dada a la partición de sí misma, entre las fases líquida y de vapor. El
valor de K es una medida de la tendencia de un componente para favorecer la
fase vapor, cuando Ki es mayor a la unidad, se dice que la especie i exhibe una
mayor concentración en la fase vapor, pero cuando es menor, presenta mayor
concentración en la fase líquida, y se considera como un componente pesado.

yi Ecuación 27
x i=
ki

1=L+V Ecuación 28

z i=x i (1−V )+ y i V Ecuación 29

Entonces despejando Xi:

zi ki Ecuación 30
y i=
1+V (k i−1)

zi Ecuación 31
x i=
1+ V (k i−1)
Los dos conjuntos de fracciones mol deben sumar la unidad ∑i xi =∑i yi =1. Se suma
la yi para todas las especies y se le resta la unidad, la diferencia Fy debe ser cero:

zi ki Ecuación 32
F y =∑ −1=0
i 1+ V (k i−1)

zi Ecuación 33
F x =∑ −1=0
i 1+V (k i−1)

El valor de V se determina con la solución de Fy - Fx=F:

z i (k i−1) Ecuación 34
F=∑ =0
i 1+ V (k i −1)

Derivando:

dF z i (k i −1)2 Ecuación 35
=−∑ 2
dV i [1+V (k i −1)]

Para solucionar este sistema de ecuaciones se utilizó el proceso iterativo de

Newton, con el fin de obtener el valor de V. Este cálculo se realiza si la
temperatura asignada se encuentra entre la temperatura de burbuja y de rocío.

Inicialmente se hará uso de los resultados de temperatura de rocío y temperatura

de burbuja, las cuales brindan el valor inicial para la determinación del coeficiente
de actividad y V. Teniendo en cuenta que para el punto de rocío V=1 y para el
punto de burbuja V=0.
El calor absorbido por un líquido para pasar a vapor se denomina calor de
vaporización. En un proceso de ebullición a presión constante el calor absorbido
por la sustancia que cambia de fase sin cambiar su temperatura viene dado por la
ecuación 1. El calor latente de vaporización se determinó empleando la entalpia
de vaporización en el punto de ebullición normal como lo establece Watson 5
ecuación 2.

Q=m∗λ Ecuación 36

Dónde:

 Q: Calor
 m: masa
 : Calor latente de vaporización.

∆ Hv 1−T r 0,38
Ecuación 37
∆ Hv
2
=
1
( )
1−T r
2

Dónde:

 Hv : Calor de vaporización de líquido a T2 (kJ/mol)

2

 Hv : Calor de vaporización de líquido a T1 (kJ/mol)

1

 T v : Temperatura
2

 T v : Temperatura
1

Q=m∗c∗ΔT Ecuación 38

Dónde:

 C: Calor especifico (kJ/mol)

 T: Calor de vaporización de líquido a T 1 (kJ/mol)

Dimensionamiento del tanque.

El tanque se dimensiona teniendo en cuenta el flujo de alimentación y la densidad

del residuo.

5
HIMMELBLAU, D. Principios básicos y cálculos en ingeniería química. Balances de energía. 6 ed. New York:
Bernard Goodwin, 1997, 420 p.
 m Ecuación 39
V=
ρ

46,79 kg /d
V= 3
=0,0575m 3 /d
813,93 kg /m

Teniendo en cuenta que el recipiente por seguridad debe de estar lleno máximo al
80%, se tiene entonces un factor de seguridad de 20%.

V nuevo =0,0575 m3∗1,2=0,0689 m 3

El diámetro y altura del tanque se calculó con base a la relación L/D y la geometría
del recipiente:

H
=2
D

π Ecuación 40
V = ∗D 2∗H
4

Dónde:

 V: Capacidad (m3)
 D: Diámetro del recipiente. (m)
 H: Altura del recipiente (m).

0,0689 m3∗2
D=
√
3

π
=0,3528 m

H=0,3528 m∗2=0,7056 m

Balance de Energia.

Para determinar la potencia del equipo es importante tener en cuenta las

ganancias y perdidas de calor en el sistema, teniendo como punto de partida dos
etapas, la primera consiste en el calor requerido para iniciar el proceso, este
involucra el calor que se necesario para calentar el equipo, el calor que ha de
recibir la mezcla para aumentar su temperatura e iniciar la separación y el calor
que se pierde, en la segunda etapa ya se tiene en cuenta el calor latente o energia
necesaria para que se de el cambio de fase e igualmente se tiene en cuanta el
calor perdido durante esta etapa. Esto se simplifica en la ecuación 41

Q=Q Equipo +QResiduo +Q Perdido+ QLatente +Q Perdido Ecuación 41

Q=m∗Cp∗∆T Ecuación 42

Dónde:

 m: masa de la sustancia (kg)

 Cp: Capacidad calorífica (J/kg K)
 T: Diferencia de temperatura (K).

El calor absorbido por el recipiente se determinó teniendo en cuenta el material del

equipo que para este caso es acero inoxidable AISI 304, las propiedades se
obtuvieron a partir de la ficha técnica de este.

Tabla 23. Propiedades Acero inoxidable AISI 304

Propiedad Valor
Densidad (kg/m3) 7.930
Capacidad calorifica (J/kg K) 510

m=V ∗ρ

kg
m=0,06899 m 3∗7930 =547,06 kg
m3

J
Q equipo =547,06 kg∗510 ∗( 350,61−288,15 ) K=1,74 x 107 J =1,21 kJ /s
kg K
El calor del residuo o calor sensible, corresponde al calor absorbido por la mezcla
para que esta aumente su temperatura desde la temperatura ambiente hasta el
punto de ebullición de la mezcla.

0,003249 kg J
Q residuo = ∗8370 ∗( 350,61−293,15 ) K =1,56 kJ /s
s kg K

0,002785 kg
Q latente = ∗8,80 kJ /kg=0,024523 kJ /s
s

Se determino la perdida de calor asociada con el cambio de fase.

0,002785 kg J
Q perdido = ∗8710 ∗( 353,15−350,61 ) K =0,0616 kJ / s
s kg K

Q=Q Equipo +Q Residuo +Q Perdido + Q Latente +Q Perdido

kJ kJ kJ kJ kJ
Q=1,21 +1,56 +1,56 +0,024523 +0,0616 =4,42 kW
s s s s s

CONDENSADOR

Inicialmente se calculó el promedio de cada uno de los parámetros obtenidos que

se muestran en el capítulo anterior, primero para determinar la media térmica
logarítmica balanceada, se tuvo en cuenta un balance para cada etapa:

En la etapa de enfriamiento, el flujo de agua se determinó con la siguiente

ecuación:

ṁvapor c p ( T h 1−T h 2 )=ṁ agua c p ( T c 2−T c1 )

vapor agua
Ecuación 43

ṁvapor c p ( T h 1−T h 2 ) Ecuación 44

ṁagua = vapor

cp agua
( T c 2−T c1 )
 Donde:

ṁvapor : Flujo másico del vapor (kg/s).

ṁ agua: Flujo másico del refrigerante (kg/s).
c p : Coeficiente de calor calorífica de la mezcla vapor (J/kg*K), a
vapor

la temperatura media (351,88 K).

c p : Coeficiente de calor calorífica del agua 6 (J/kg*K), a la
agua

temperatura media (288,15 K).

T h 1: Temperatura de roció a la entrada del condensador.
T h 2: Temperatura de burbuja a la salida del condensador.
T c1 : Temperatura de entrada del refrigerante.
T c2 : Temperatura de salida del refrigerante.

0,002785417∗8,71E+03∗( 353,15−350,61 )
ṁ agua =
4186∗( 293,15−283,15 K )

ṁagua =0,001472845 kg/ seg

El coeficiente de calor calorífico de la mezcla líquida a T h 1, se determinó a partir de

la siguiente expresión7:
NC
Ecuación 45
Cpm=∑ X i∗Cpi
i=1

Cpm= ( ( 1,54E+02∗0,34217 ) + ( 2,19E+02∗0,10633 )+ ( 1,15E+02∗0,47973 )+ ( 1,84E+02∗0,04536 ) + ( 2,24E+02∗0

Cp m=8,73E+03 J /mol k

Para el etanol:

Cpliq =A + ( B∗T ) + ( C∗T 2 )+ ( D∗T 3 ) Ecuación 46

Cpliq =59,342+ ( 0,3636∗353,15 )+ ( −1,22E-03∗353,152 ) + ( 1,80E-06∗353,153 )

6
Calor específico del agua. Available from Internet <
[Link]
7
RODRIGUEZ, Jorge. Estimación de propiedades termodinámica. Modelado, simulación y optimización de
procesos químicos. [Link] from Internet
<[Link]
Cpliq =1,15E+02( J /mol K)

La viscosidad de la mezcla vapor a T h 1, se determinó a partir de la siguiente

ecuación8:
n
Ecuación 47
ln μm =∑ X i ln μi
i=1

ln μm =( ( 0,34217∗ln ( 0,637621986 ) )+ ( 0,10633∗ln ( 0,284271055 ) ) + ( 0,47973∗ln ( 0,43003087 ) ) + ( 0,04536∗ln ( 0

ln μm =0,00046364 Pa. seg

La densidad de la mezcla a T h 1, se determinó a partir de la siguiente expresión 9:

1 Ecuación 48
ρ=
Xi
∑ ρi

1
ρ=
0,34217 0,10633 0,47973 0,04536 0,02641
+ + + +
0,749631994 0,820119515 0,732356437 0,824585037 0,863750015
ρ=753,7284198 kg/m3

La conductividad térmica de la mezcla vapor a T h 1, se halló con la siguiente

ecuación10:
n n
Ecuación 49
k m=∑ ∑ ∝i ∝ j k ij
i=1 j=1

k m=( 0,377225063∗0,352241641∗0,074067758 ) + ( 0,068283728∗0,06112043∗0,377225063 ) + ( 0,078989169∗

k m=0,02490725 W /m K

En cuanto al balance de energía para el vapor, se utiliza la siguiente ecuación:

8
Algunos Métodos de Estimación para Viscosidad. Available from Internet <
[Link]
9
RODRIGUEZ, Jorge. Estimación de propiedades termodinámica. Modelado, simulación y optimización de
procesos químicos. [Link] from Internet <
[Link] >.
10
Algunos Métodos de Estimación para conductividad térmica. Available from Internet <
[Link]
 q h=ṁvapor c p ( T h 1−T h 2 )
vapor
Ecuación 50

q h=0,002785417∗8,71E+03∗(353,15−350,61)
q h=61,65330843 W
Para la media térmica logarítmica, considerando que el intercambiador opere en
contracorriente, se usa la siguiente ecuación:
∆ T 1−∆T 2 Ecuación 51
MLDT =
∆T 1
ln ( )
∆T 2

60−67,46
MLDT =
60
ln(67,46 )
MLDT =63,65716339 K

Para estimar el área a partir de un supuesto de coeficiente global, para el caso se

toma un valor de U D = 350 W/m2 k usado para líquidos orgánicos, mediante la
siguiente ecuación:

q 61,65330843W
A= = =0,002770102 m2
U D∗MLDT 350∗63,65716339

Ahora para determinar el factor de corrección ( F T ), se deben hallar los parámetros,

mostrados a continuación:
∆T 2 353,15−350,61
R= = =0,254
∆T 1 293,15−283,15
T c 2−T c1 293,15−283,15
S= = =0,142857143
T h 1−T c1 353,15−283,15
1−S
√ R2 +1∗ln 1−RS
F T= =0,998953715
2−S ( R+1− √ R2 +1 )
( R−1 )∗ln
2−S ( R+1+ √ R2 +1 )
 En cuanto a los tubos, se seleccionó OD (In) de 1/2" con un diámetro externo de
0,0127m y diámetro interno de 0,0066m del documento Boiler Tube Expanders 11.
Por conveniencia se seleccionó arreglo triangular, para el cual se determinó los
siguientes parámetros:

Numero de tubos12:

Nt=
0,875∗ ( CTP
CL )
∗DC 2

2
Pt

0,9
Nt=
0,875∗ ( 0,87 )∗8,875108761 =8,875108761
2

2
0,04

Donde:

CTP: Se refiere a la constante que toma en cuenta los espacios entre la coraza,
los tubos y los pasos del lado del tubo.

DC: Diámetro equivalente de la coraza en metros.

CL: Constante de distribución.

Pt: Espaciado de los tubos, metros.

También se tuvo en cuenta la distancia imaginaria o claro, hallada de la siguiente

manera13:

1,6+(0,004∗D is )
C=
1000

1,6+ ( 0,004∗203 )
C= =0,002412 m
1000

En cuanto a la coraza, conociendo que es un paso por el tubo y dos por la coraza,
se determinó el espacio entre deflectores (B) con la quinta parte del diámetro
11
INC,Thomas [Link]. Boiler Tube [Link] York.
12
ONTIVEROS, Jesús F. Operaciones unitarias. Dimensionamiento de Equipos de Transferencia de Calor.
Available from Internet < [Link]
13
OSORIO AGUILAR, Rita y GONZÁLEZ FLORES, Rafael. Metodología para el diseño de un intercambiador
de calor de tubo de coraza. En: Revista de la Asociación Mexicana de Metodología de la Ciencia y de la
Investigación, A.C. Junio, 2010. vol.2, no.1, p.56-75.
interno de la coraza, dando como resultado 0,0406m. Luego se halla el área de
flujo por el lado de la coraza, así:

Ds∗C∗B 0,203∗0,002412∗0,0406
a s= = =0,000496981m2
N PS∗Pt 1∗0,04

Mientras que el diámetro equivalente de la coraza para arreglo triangular, se

determina con la siguiente expresión14:

((
2

De =
1
4∗ ∗Pt∗0,86∗Pt −
2
2
)4 (
∗π∗( DE ) ∗1
)
1
∗π∗( DE )
2

((
2

De =
1
4∗ ∗1,575∗0,86∗1,575 −
2
2
4 ) (
∗π∗( 0,5 ) ∗1
)
1
∗π∗(0,5)
2

D e =4,932499335 pulg=0,125250993 m

Donde:

DE: Es el diámetro exterior del tubo, pulg.

El área de flujo de los tubos (a t ¿ , se determinó con la siguiente ecuación:

N t∗π
∗Di 2
4
a t=
N Pt

8,875108761∗π
∗0,00662
4
a t=
2

a t=0,000151817 m 2

14
ARANDA DELGADO, Carlos Eduardo. Diseño de un intercambiador de calor de tubos y coraza para fluidos
líquido –liquido. Perú: Universidad Nacional del Callao; 2014.
Además, se determinan las propiedades a la temperatura media, junto con los
siguientes modelos:

La velocidad másica de la mezcla, se halla con la siguiente expresión:

m∙
G=
a

0,002785417
G=
0,000496981

G=5,604679545 kg /s m2

Donde:

m ∙ : Es el flujo másico del vapor, kg/s.

a: Área transversal para el lado de la corza,m 2

El número de Reinolds ( ℜ¿, para la mezcla, se determina con la siguiente

ecuación:

DG
ℜ=
μ

0,125250993∗5,604679545
ℜ=
0,000472431

ℜ=1485,913301

Donde:

D: Diámetro equivalente de la coraza.

G: velocidad másica de la mezcla.

μ: Viscosidad de la mezcla.

El número de Prandtl (Pr), para la mezcla, se determina con la siguiente

expresión:

Cp∗μ
Pr=
k

8714,375254∗8714,375254
Pr=
0,024972955

Pr=164,8560225
El número de Nusselt (Nu), para la mezcla, se halla de la siguiente forma:

Nu=0,36∗ℜ0,55∗Pr 1/ 3

Nu=0,36∗1485,9133010,55∗164,85602251 /3

Nu=¿ 109,6311194

Mientras que los coeficientes de convección para la mezcla y el agua, se

determinan mediante las siguientes ecuaciones:

Nu∗ρ
h=
G

109,6311194∗755,0990214
h=
5,604679545

h=14770,22019W /m2 ° C

h∗G
hie =
Di

117,1315397∗9,701429275
hie =
0,0066

hie=¿ 172173,2346 W /m 2° C

Donde:

ρ :Es la densidad a temperatura promedio

El coeficiente global limpio ( Uc). Se calcula con la siguiente ecuación:

h ie∗h h
U c=
hie + hh
172173,2346∗14770,22019
U c=
172173,2346+14770,22019
U c =13603,23948 W /m 2° C

Finalmente, como tanto el RD calculado como el Ud son mayores que los

propuestos se acepta el diseño térmico.
REBOILER

Para el diseño del rehervidor, se seleccionó el aceite Mobiltherm 605, para

calentar más rápido el sistema, por sus excelentes propiedades en cuanto a
conductividad y calor especifico15.
El flujo del aceite diátermico se determino con la siguiente expresión 16:

( F∗Cp )∗( T 2−T 1 )

M=
H 3−H 4

( 46.79∗8.37E+03 )∗( 350.61−293.15 )

M=
224591.4−82797

M =158.70 kg /dia

Donde:

F: Flujo masico de la mezca, kg/d.

Cp: Coeficiente de calor calorífico de la mezcla líquida a temperatura medía,

J/kg*K.

T 2: Temperatura de burbuja de la mezcla líquida, K.

T 1: Temperatura ambiente, K.

H 3:Entalpia del aceite diatermico a temperatura 127.46 °C.

H 4 :Entalpia del aceite diatermico a temperatura 62.46 °C.

El calor ganado se determino con la siguiente ecuación:

Q g=F∗Cp∗( T 2−T 1 )

15
Mobiltherm 605. ExxonMobil Corporation. Esso Petrolera Argentina S.R.L.
16
BARDERAS VALIENTE, Antonio. Capítulo VII rehervidores. Ingeniería del calor. Available from Internet
<[Link]
 46.79∗8.37E+03∗( 350.61−293.15 )
Q g=
1∗3600

Q g=6,250.89W

El calor perdido se halla con la siguiente expresión:

Q p=m.∗( H 3 −H 4 )

Q p=158.70∗( 224591.4−82797 )

Q p=6,250.89 W

Los resultados del calor denotan que se cumple con el balance de energia.

Para la media térmica logarítmica, se utiliza la misma ecuación que para el

condensador:

∆ T 1−∆T 2
MLDT =
∆T 1
ln ( )
∆T 2
50−42.46
MLDT =
50
ln(42.46 )
MLDT =46 ° C
Luego, se realizó un proceso iterativo, en la primera iteración, se halló el área del
rehervidor, para lo cual se supuso el coeficiente global de transferencia de calor
(250 W/m 2 °C), con la siguiente ecuación:
Qp
A=
U∗MLDT
6,250.89
A=
250∗46

A=0,543555361m2

Donde:

Q p: Calor perdido, W.

U: Coeficiente global de transferencia de calor, W/m2°C.

MLDT : La media térmica logarítmica, ºC.

A partir de los resultados anteriores, se determinó el diámetro interno, con la

función Solver, utilizando la siguiente función objetivo:
fobjetivo=( 2∗π∗r )∗( r +4∗r )

fobjetivo= ( 2∗π∗Dint
2 )∗(
Dint 4∗Dint
2
+
2 )
Donde:
R: Es el radio interno del rehervidor, m.
Dint: Es el diametro interno del rehervidor, m.

Para la segunda iteración, se realizan los siguientes cálculos:

 El coeficiente de transferencia de calor para el fluido interno (aceite) se

calcula con la siguiente expresión17:
1 0.14
jk ( μ
hic =
Dt ( )
NPr ) 3∗
μw

1
260∗0,1315889 ( 0.14
hic = 254,2560765 ) 3∗( 1 )
0,263073335

hic =823,8969361 W /m2 º C

Donde:

hic:Coeficiente convectivo o de pelicula interno, W /m2 º C .

j:Factor J18.

k: Conductividad térmica del aceite diátermico, W/m°C.

Dt:Diametro interno, m.

NPr: Nùmero de Prandtl.

μ : Viscosidad dinámica a la temperatura calórica,kg/m s.

μw : Viscosidad dinámica a la temperatura de la pared del tubo, kg/m s.

17
GODOY LEMUS, Pablo Arturo. Evaluación de alternativas para la recuperación de calor en un reactor de
surfactante en la industria de agroquímicos. Guatemala: Universidad de San Carlos de Guatemala; 2003.
18
KERN, Donald Q. Procesos de transferencia de calor. Coeficientes de transferencia de calor para chaquetas
y serpentines. p. 816.
El número de Prandtl (NPr), se determina con la siguiente ecuacion:

Cp∗ϑ∗ρ
NPr=
k

1928∗0,000020397∗850,78
NPr=
0,1315889

NPr=254,2560765

Donde:

Cp: Coeficiente de calor calorífico del aceite a la temperatura medía, J/kgºC.

ϑ: Viscosidad cinematica a la temperatura media, m 2 /s .

ρ : Densidad del aceite diátermico a la temperatura media,kg /m3.

k: Conductividad térmica del aceite diátermico a la temperatura media, W/m°C.

Con el anterior coeficiente de transferencia de calor ( h ic ), se determina el

coeficiente de transferencia de calor para el fluido exterior del tubo (hio ) , con la
siguiente expresión:

hic∗DI
hio =
DE

823,8969361∗0,263073335
hio =
0,269173335

hio =805,2258025W /m2 º C .

Donde:

W
hic : Coeficiente convectivo o de pelicula interno, ºC.
m2

DI: Diámetro interno,m.

DE: Diámetro externo se tuvo en cuenta un Calibre 11-AISI 304 ,m.

 El coeficiente global limpio ( Uc). Se calcula con la siguiente ecuación:

h io∗hic
U c=
hio +hic
805,2258025∗823,8969361
U c=
805,2258025+823,8969361
U c =407,2271879W /m 2 ° C

 Para el área, se empleó la siguiente ecuación:

Qp
A=
Uc∗MLDT
6,250.89
A=
407,2271879∗46
A=0,332771762m2
 El número de Reynolds ( ℜ¿ , se determina con la siguiente ecuación:
D i∗V
ℜ=
ϑ
0,263073335∗1,5
ℜ=
0,000020397
ℜ=19346,47264

Donde:

Di: Diámetro interno, m.

V:Velocidad del aceite diátermico,m/s.

ϑ :Viscosidad cinematica a la temperatura media 19, m2 /s .

 A partir de los resultados anteriores, se volvió a determinar el diámetro

interno, con la función Solver, utilizando la siguiente función objetivo:
fobjetivo=( 2∗π∗r )∗( r +4∗r )

fobjetivo= ( 2∗π∗Dint
2 )∗(
Dint 4∗Dint
2
+
2 )
Con la misma forma anteriormente mostrada, se repitio el proceso iterativo 14
veces más, hasta que convergío el área, dando como resultado el área del reboiler
de 0,20665 m 2.

Respecto a las resistencias, se realizan los siguientes calculos 20:

 La potencia, se halla con la siguiente igualdad:

P=Q
19
Centro de tecnología de transferencia de calor. Ingeniería Industrial Aeronáutica de Terrassa. p.D 2.
20
RINALDO, Paul. Guía internacional del [Link]. Marcombo [Link]ña.1995. Available from
Internet< [Link]
id=frr2fA0T3ssC&pg=PA49&dq=unidad+voltios+sobre+amperios&hl=es-
419&sa=X&ved=0ahUKEwje0YOUlYrUAhXKKyYKHfPEAIwQ6AEIIjAA#v=onepage&q=unidad%20voltios
%20sobre%20amperios&f=false >
 6,250.89
P=
24∗2600
P=0.07234822517Wattios
 La corriente se determina con la siguiente expresión:
P
I=
V
0.07234822517
I=
24

I =0.003014509382 Amperios

 La resistencia, se halla con la siguente ecuación:

V
R=
I

24
R=
0.003014509382
R=7961.494544 Ohmios

Donde:

P: Potencia del rehervidor, W.

V:Voltaje del rehervidor,V

CHAQUETA DE CALENTAMIENTO

La altura de la chaqueta se determinó con la siguiente expresión:

ht
h=
1+ f

0,7056
h=
1+0,10

h=0,64 m

Donde:

ht : Altura del tanque, m.

f : Factor de seguridad.
El número de Reynolds ( ℜ¿ , se determina con la siguiente ecuación:

D i∗V
ℜ=
ϑ
0,35962∗1,50000
ℜ=
5,30E-06

ℜ=1,02E+05

Donde:

Di: Diámetro interno de la chaqueta, m.

V:Velocidad del aceite diátermico,m/s.

ϑ :Viscosidad cinematica a la temperatura media, m2/s.

Al ser el regimen turbulento, se halla f con la siguiente ecuación:

f =( 0,79∗ln ( ℜ ) −1,64 )−2

f =0,02

Mientras que el número de Nusselt (Nu), se determino con la siguiente ecuación:

f
∗( ℜ−1000∗Pr)
8
Nu=
f 0,5 ( 0,666666666666667 )−1
1+12,7∗
8 () ∗ Pr

0,02
∗(1,02E+05−1000∗14,00)
8
Nu=
0,02 0,5 (
1+12,7∗
8( ) ∗ 14,000,666666666666667 )
−1

Nu=812,49

Donde:

f : Expresión dependiendo el regimén.

Re: Número de Reinolds.

Pr: Número de Prandtl.

AISLANTE
Se seleccionó como aislante de la chaqueta de calentamiento, la fibra de vidrio
para disminuir las perdidas de calor, ademas de ser un material que se usa
frecuentemente para dicho objetivo. Mediante la siguiente ecuación se halló la
resistencia termica del aislante ( Raislante ¿:

x a x F .V x ac
Raislante = + +
k a k F . V ka

0,0068 0,0254 0,0037

Raislante = + +
15,00 0,04 15,00

Raislante =0,6357000 m2 ° C /W

Donde:

x a:Espesor de la lámina de acero inoxidable AISI-304,m.

x F .V : Espesor de la fibra de vidrio, m.

x ac : Espesor de la chaqueta de calentamiento Calibre 15 AISI-304,m

k a : Conductividad térmica del acero inoxidable AISI-304, W/m°C.

k F . V : Conductividad térmica de la fibra de vidrio, W/m°C.

El calor con aislante, se hallo con la siguiente ecuación:

2∗π∗h∗(T 1−T 2 )
q aislante =
r r
ln 2 ln 3
1 r1 r2 1
+ + +
h1∗r 1 k a k F . V h2∗r 3

2∗π∗0,7056∗(127,46000−20,00000)
q aislante =
0,18161 0,20701
ln ln
1 0,17981 0,18161 1
+ + +
289,19550∗0,17981 15,00 0,04 5,00000∗0,20701

q aislante =111,87W /m

Donde:

h:Altura del tanque, m.

T 1: Temperatura del aceite, °C.

T 2: Temperatura del medio ambiente, °C.

h1 : Coeficiente de transferencia de calor por convección del aceite, W /m2 ° C .

h2 : Coeficiente de transferencia de calor por convección del aire, W /m2 ° C .

r 1: Radio externo del tanque de destilado, m.

r 2: Radio externo de la chaqueta de calentamiento, m.

r 3: Radio externo de la chaqueta de calentamiento y espesor de la fibra de

vidrio ,m.

k a : Conductividad térmica del acero inoxidable AISI-304, W/m°C.

k F . V : Conductividad térmica de la fibra de vidrio, W/m°C.

El calor sin aislante, se hallo con la siguiente ecuación:

2∗π∗h∗(T 1−T 2)
q sinaislante =
r
ln 2
1 r1 1
+ +
h1∗r 1 k a h2∗r 2

2∗π∗0,7056∗(127,46000−20,00000)
q sinaislante =
0,18161
ln
1 0,17981 1
+ +
289,19550∗0,17981 15,00 5,00000∗0,18161

q sinaislante =424,94 W /m

Tanque

Se utilizo la siguiente relación:

h
=2
D

A
D=
√ π∗1,25
 0.0001509876004
D=
√ π∗1,25

D=0.006200699695m

Mientras que la altura se hallo con la anterior ecuación:

h=2∗D

h=2∗0.006200699695

h=0.01240139939 m

4∗w∗c 13 μ 0.14

hi =
( (
1.86∗
π∗K∗L ) ( ) )∗K
∗
μw
D
1 0.14
4∗350∗0.5044 26.80

hi =
( (
1.86∗
π∗0.0729∗12.41620763
3
∗ ) (
25.295049 ) )∗0.0729
0.15

hi =454.0770792 BTU /h ft 2 ° F

Donde:

w: Flujo másico del aceite diátermico, lb/h.

c: Calor específico del fluido21, BTU/lb°F.

K: Conductividad térmica, BTU/h ft °F.

L:Longitud del rehervidor, ft.

μ : Viscosidad a la temperatura calórica22 , lb/ft h.

μw : Viscosidad a la temperatura de la pared del tubo,lb/ft h.

Se determina el coeficiente de diseño (Ud), por medio de la siguiente expresión:

1
Ud=
1
+ Rd
U

21
Regenerador de disolventes. Manual de instrucciones, uso y mantenimiento. Italia sistema tecnologici.
22
WRIGHT, Wes. Simple equations to approximate changes to teh properties of crude oil with changing
[Link] from Internet< [Link]
temperature?page=3#.WScJWus1-aG >.
 1
Ud=
1
+0.002
250

U d =166.6666667 W / m 2° C

Donde:

Rd : Es el factor de ensuciamiento, adimensional.

El coeficientes de transferencia de calor para el fluido exterior (hio) ,se halla con la
siguiente expresión23:

hi∗Di
hio =
De

200∗0.15
hio =
0.006200699695

hio =9831.14858W /m2 ° C

La temperatura de pared para los fluidos caliente y frio respectivamente,se

determinan con los siguientes modelos:

T hP=T h 2+ F c∗(T h 1−T h 2)

T hP=62.46+ 0.48492∗(127.46−62.46)

T hP=93.9798 ° C

T cP=T c 1+ F c∗(T c2−T c 1)

T cP=20+0.48492∗(77.46−20)

T cP=47.8635032° C

Donde:

Fc: Fracción calórica, adimensional24.

23
RAMIREZ ARBOLEDA, Lina María. Obtención de oseina, fasfato dicalcico y grasa a partir de hueso de res en
la empresa productos alimenticios Don rafa LTDA. Manizales: Universidad Nacional De Colombia; 2003.
24
GONZÁLEZ, Mahuli [Link] de calor. Universidad nacional experimental Francisco de Miranda
Available from Interne < [Link]
[Link] >.
Ahora , la temperatura calorica mencionada anteriormente, se determino con la
siguiente expresión:

T c =T 2 + F c∗(T 1−T 2)

T c =62.46+0.48492∗(127.46−62.46)

Sistema de Tubería.

El sistema consiste en xx m de tuberia de acero inoxidable para el tramo 1 con un

diametro interior de 1” y el tramo 2 con un diametro interno de ½”.

Inicialmente se tuvo en cuenta

v∗Di∗ρ
ℜ=
μ

Dónde:

 v: Velocidad del fluido (m/s)

 Di: Diámetro de la tubería (m)
 : Densidad de la mezcla (kg/m3).
 : viscosidad de la mezcla (Pa*s)

A flujo =π∗r 2

A flujo =π∗( 0,01215 m )2=0,0004638 m 2

Q
v=
A flujo

Dónde:

 v: Velocidad del fluido (m/s)

 Q: Caudal del fluido (m3/s)
 A: Área del flujo (m2)

3,682 x 10−6 m3 / s
v= =0,0079396 m/ s
0,0004638 m2
 m kg
0,0079396 ∗0,0243∗827,76 3
s m
ℜ=
0,00177 Pa∗s
Recomendación.

Caracterizar y realizar un estudio del residuo solido proveniente de la destilación

con el fin de buscar si este residuo puede llegar a tener algun uso.