UNIVER9IDA II AUTONOMA METROPOLIT.

4 NA
UNIDAD IZTAPALAPA
DIVISION CBI

DEPARTAMENTO DE I.P.H.
AREA D E RECURSOS EVEGETICOS
I,iCENCIA TURA INGENIERIA EN ENERCIA

SIMULA CION DINAMICA

DE
TRENDE EK4YORADORES

Enrique Z(wmri Pano.

Aiurnno e la licenciritmn (le Ingeniern en Energin

Mrtrcuia: 92323 783

Mxico, D. F., Noviembre de 1 9 3 .

Serriintirio de proyectos.

hnrique Zamora Pano.

Universidad Autnoma Metropolitana, Iztapalapa.
Alumno de la licenciatura de Ingeniera en Energa.
Alejandro Torres Aldaco.
Universidad Autnoma Met copolitana, Iztapalapa.
Asesor, Profesor Investigador.

UAM-Mapa 1ap a
Divisin de Ciencias Bsicas e Ingenierta
Are a de Ingenie ra en Recursos Energticos
Departamento de Ingenierta de Procesos e Hidrulica
Mxico, D.F.
1998.

--

Indice+
~

Tema:
Prologo.
Capitulo 1. Evaporacin
..
1.1. Introduccion ................................................................
.
.
1.2 Objetivo .......................................................................
1.3. Evaporacin ................................................................

Capitulo 11. Tipos de evaporadores y accesorios.

11.1, Generalidades de los evaporadores .........................
11.2. Tipos de evaporadores ....................................................
11.2.1. Evaporador de tubos horizontale
II.2.2. Evaporador de tubos verticales .....................
11.2.3. Evaporador de pelcula agitada ............
.................
11.3. Accesorios de los evaporadores .....
11.3.I . Condensadores ...................
II.3.2. Purgadores ...............

Pag.

5
6
8
12
13
13
I4

Capitulo 111. Capacidad y ecom ina de un evaporador.

ILL 1. Capacidad de un evapordor .............
TII.l.l. Diferencia de temperatura ...........
HI.1.2. Elevacin del punto de ebullicin .........
................ 17
TII. 1.3. Coeficiente de transmisin de calor
III.1.4. Coeficiente de la pelcula del vapor
IIL1.5. Coeficientes del lado del lquido .
HI.1.6. Coeficientes globales ........................
................ 20
111.2. Economa de un evaporador .............................................
20
,,

Capitulo 1V. Evaporador de simple efecto.

IV. 1. Evaporador de simple efecto ............................................
21
lV.2. Balances de entalpa de un evaporador de simple efecto ............ 23

Capitulo V. Evaporador con mltiples efectos.

V. 1. Evaporador de mltiples efectos .....................
V.l.l. Funcionamiento de un mltiple efecto ...................
V. 1.2. Capacidad de evaporacin y superficie de
V.1.3. Mtodos de alimentacin ...
V.2. Capacidad y economa de evapo
V.3. Clculo de un mltiple efecto ...........
V.4. Nmero ptimo de efectos .....
................

Capituio VI. Elevacin de la temperatura de ebullicin.

VI.1, Efecto de la carga del lquido y de la elevacin de la tempertura
de ebullicion ..............................................................

31

Capitulo VIL Recompresin del vapor.

VI1.1.Recompresin del vapor ...............................................
VII.l.l.Recompresin mecnica ..........................................
., trmica .............................................
VII.1.2.Recompresion

33
33
34

. I

Capitulo VIII. Metodologa de Simulacin.

., ..............................................................
VIII.1. Introduccion
VIII.2.Mtodo de Newton ....................................................
VIII.3,Implementacin al cdigo de simula.cin ...........................
VIII.4. Aplicacin y ejemplos ................................................
Conclusiones.
Bibliografa.
Anexo A. Cdigo de simulacin de tren de evaporadores.
Anexo B. Propiedades termodinmicas del va.por de agua.

37
37
40

42

A la memoria de:
Jos Luis Zcrmorn Hernandez.
Dicen que a travs de los hijos es coiiio se obtiene
la inmortalidad, no s realmente si esto sea cierto, pero
en lo que a m respecta vives en mi corazn y
te llevo siempre en mis recuerdos.

Y por ms que deseo encontrar,

no existe palabra o frase que pueda expresar por s misma
lo que realmente siento, sin embargo, quiero expresarte
mi gratitud por todo l o que lograste hacer de m y por lo que nunca
me hizo falta an en la disiancia, un Padre y un Amigo.

Este libro pretende introducir al tema de ewprndoores a aquellas personas que por primera
vez abordan el tema, y para estudiantes di: ingeniera dar un breve resumen de la
importancia de los evaporadores como parte i i:egral de las operaciones unitarias.
Un objetivo ms especifico en esta obra es presentar un cdigo de simulacin de trenes de
evaporadores para facilitar todo clculo matemtico relacionado con las variables del
sistema.
Esta obra abarca sohrnente en su contexto a lcls evaporadores, lo que facilita su et:tudio con
el marco tebrico y aplicacin con el cdigo de simulacin. En principio se presentan un
concepto fundamental de la evaporacin como proceso fsico pasando a la cisscripcin de
equipos de evaporacin existentes y sus accesorios como conjunto y sistema,
posteriormente se describen los arreglos posibles y nmero de etapas; as como capacidad y
economa de los evaporadores, y por ltimo el cdigo de simulacin del tren de
evaporadores.
Se ha logrado en este texto darle un seguimiento segn las necesidades 3. un tema dentro de
un curso de operaciones unitarias o de fenmenos de transporte, los requerimientos

mninios para la compresin total es derivado de cursos de ecuaciones diferenciales que

pueden llevar no solo estudiantes de ingeniera en energa, sino tambi6n estudiantes de
ingeniera qumica.

Los ejemplos que se anexan en el captulo VIII, se dan corno apoyo para la introduccin en
el uso de el simulador; en ste mismo capitulo se comparan mtodos algebraicos con
mtodos nurnericos implementados al programa de simulacin
En las ltimas hojas se anexan tablas termodinmicas de las propiedades fsicas del vapor y
un diagrama de Mollier, que nos servirn para un manejo mucho ms cmodo del programa

de simulacin.
Quiero agradecer sinceramente a compaeros estudiantes que por sus valiosos consejos y
sugerencias aportaron ideas oportunas para la realizacion de esta obra, tambin agradezco a
mi profesor asesor por sus comentarios precisos que lograron darle profundidad y seriedad
a esta obra.

E. Zamora Pano.

_I_-.

Cnptrtlo I

Evap oracin.

O bj et ivo.

Preseiitar i i n marco tedrico del concepo de eixporacili y lit) pchornmn general de

los problemas tpicos que se preseir~nridebido n la evaporacin.

Capitulo I

I.1. INTRUDUCCION.

La evaporacin de lquidos, como proceso fisico, es solo un aspecto particular de la

transmisin de calor. Se utiliza este proceso para la separacin, por ebullicin, de una parte
del lquido contenido en una disolucin o suspensin.
El calor necesario para ello puede proceder de: cualquier medio de calefaccin. Cuando este
medio es el vapor de ag::a condensante, el aparato se llama evaporudor,y la evaporacin se
estudia separadamente como operacin bsico,. La clasificacin de la evaporacin dentro de
la ingeniera qumica responde al empleo de aparatos especiales y mtodos particulares.
En una operacin de evaporacin se concentra una solucin al eliminar disolvente por
ebullicin. Por lo general, el producto deseacio es la solucin concentrada, pero en algunns

ocasiones, el producto principal es el disalvente evaporado. Us posible continuar la

concemtracin hasta la saturacin de la solucion y ms aun, hasta que el soluto precipita
como un slido cristalino. En este caso, I operacin se conoce con frecuencia como
cristnlizrrciri.

En realidad, el proceso de evaporacin es equivalente a una separacin en una sola etapa.

Sin embargo, los componentes se separan en forma tan definida que con frecuencia se
olvida la composicin del equilibrio. En aquellas aplicaciones en las que la fase vapor
constituye un producto valioso, este equilibrio puede ser importante. A pesar de la multitud
de procesos de velocidad que ocurren, el ingeniero interesado en estas operaciones casi
siempre puede considerar el proceso comp'ieto en trminos de transferencia de calor del
calentador a la solucin. Este paso controla la velocidad total y se analiza con mayor
facilidad. L a s contribuciones de los otros procesos de transferencia tienen un efecto menor
sobre la operacin total. Adems, se conoc(: tan poco, fuera de las escasas compaas que
coristniyen evaporadores, de la mayor parte de los detalles en el diseiio, que no existe razn
para tratar aqu ms que los balances de caior y material, y la formulacin simple de la
velocidad de transferencia de calor a travs de la superficie de calentamiento.
En el captulo ltimo de ste trabajo se describe como se logra adaptar la simulacin a la
evaporacin en una y mltiples etapas mediante un cdigo desarrollado para encontrar
variables propias del sistema.

Copittilo I. Evaporacin.

I. 2. OIIJETIVO.

El objetivo de la evaporacin es concentrar una disolucin consistente en un soluto no

voltil y un disolvente voltil. E11 la mayor parte de las evaporaciones el disolvente es agua.
La evaporacin se realiza vaporizando una paite del disolvente para prodi zir una solucin
concentrada. La evaporacin difiere del secado en que el residuo es un lquido (a veces
altamente viscoso) en vez de un slido; difere de la destilacin en que el vapor es
generalmente un solo componente y, aun cuando el vapor sea una mezcla, en la
evaporacin no se intenta separar ei vapor en fracciones; difiere de la cristalizacin en que
su inters reside en concentrar una disolucin y no en formar y obtener cristales.

Normalmente, en evaporacin el producto valioso es el :quido concentrado mientras que el

vapor se condensa y se desecha. Sin embargo, en algn caso concreto puede ocurrir lo
contrario. El agua salubre se evapora con frecuencia para obtener un producto exento de
slido para la alimentacin de calderas, para procesos con requerimientos especial,:s o para
el consumo humano. Esta tcnica se conoce con frecuencia con el nombre de desticrcin de
clgrrn, pero se trata en realidad de evaporacin. Se han desarrollado procesos de
evaporacin a gran escala utilizndose para la recuperacin de agua potable a partir de agua
de mar. En este caso el agua condensada es el producto deseado. Solamente se recupera una
fraccin del agua contenida en la alimentacin, mientras que el resto se devuelve al mar.
La evaporacin por cualquier otro procedimiento puede estudiarse h t r o de la tranmisin
de calor, siempre que no presente aspectos fundamentales que justifiquen un tratamiento
independiente.
I. 3. E VAP O W CION.

Como se analiz con anterioridad, la evaporacin es la operacin de concentrar una

solucin mediante la eliminacin de di:jolvente por ebullicin. Por lo general la
concentracin se detiene antes que el soluto comience a precipitarse de la solucin. En
consecuencia, un evaporador debe consistir bsicamente de un intercambiador de calor
capaz de hervir la solucin y un dispositivo para separar la fase vapor del lquido en
ebullicin. En su forma ms simple, puede ser una charola de lquido colocada sobre una
placa caliente. La superficie de la placa caliente es un intercambiador de calor simple y el
vapor se desprende en la gran rea para flu)o de vapor y su consecuente baja velocidad de
flujo. En la operacin industrial, por lo ger,eral el equipo se construye para una operacin
continua, la superficie de intercambio de calor se incrementa de un modo notable, la
ebullicin es sensiblemente ms violenta y la evolucin del vapor es rpida. Se encuentran
problemas tales como formacin de espuma, incrustacin, sensibilidad al calor, corrosin y
limitaciones de espacio, Estos problemas han resultado en variaciones y refinamientos del
diseo de evaporadores para lograr diferentes combinaciones de propiedades de la solucin
y condiciones econmicas.

Capi[irlo I. Evariorncion.

Para poder entender los problemas que se presentan dentro de la evaporacin estudiemos de
manera general aquellas limitaciones que posteriormente servirn como antecedente para la
correcta eleccin de cualquier tipo de evaporador y SLI correcto funcionamiento.
Cnrncterslicns del lquido. La solucin prctica a un problema de evaporacibn est
profundamente afectada por el carcter del lquido que se quiere concentrar. Precisamente
es la gran varier! id de caracteristicds de dicho:; lquidos (que demanda criterio y experiencia
en el diseo y operacin de evaporadores) lo que ampla esta operacin desde una sencilla
transmisin de calor hahta algo mucho ms; complicado. Ms adelante mencionaremos
algunas propiedades de los lquidos que se eviipn'an.
Concentracin. Aunque la disolucin que entra como alimentacin de un evaporador
puede ser suficientemente diluida teniendo riuchas de las propiedades fsicas del agua, a
medida que aumenta la concentracin de la cisolucin adquiere cada vez un carcter ms
particular. La densidad y la viscosidad aumentan con el contenido de slido hasta que la
disolucin o bien se transforma en saturada o resulta inadecuada p<ra una transmisin de
calor adecuada. La ebullicin continua de una dsulucin saturada da lugar a la formacin
de cristales, que es preciso separar pues de lo contrario obstruyen los tubos. La temperatura
de ebullicin de la disolucin puede tambin aumentar considerablemente al aumentar el
contenido del slido, de forma que la temperatura de ebullicin de una disolucin
concentrada puede ser mayor que la del agua a la misma presin.
Formacin de esplima. Algunos materiales, especialmente sustancias orgnicas, forman

espuma durante la vaporizacin. Una espuma estable acompaa al vapor que sale del
evaporador dando lugar a un importante arrastre. En casos extremos toda la masa de lquido
puede salir con el vapor y perderse.
Sensihilidad a fa temperdim.

Muchos productos qumicos finos, farmacuticos y

alimenticios se daan cuando se calientan a temperaturas moderadas durante tiempos
relativamente cortos. En la concentracin de estos materiales se' necesitan tcnicas
especiales para reducir tanto la temperatura del lquido como el tiempo del calentamiento.

Formacin de costras. Algunas disoluciones depositan costras sobre las superficies de

calefaccin. En estos casos el coeficiente global diminuye progresivamente hasta que llega
un momento en que es preciso interrumpir la operacin y limpiar los tubos. Cuando las
costras son duras e insolubles, la limpieza resulta dificil y costosa.
Materiales de construccin. Siempre que c:s posible, los evaporadores se construyen con
algn tipo de acero. Sin embargo, muchas disoluciones atacan a los metales frreos y se

produce contaminacin. En estos casos se utilizan materiales especiales tales como cobre,
nquel, acero inoxidable, grafito y plomo. Puesto que estos materiales son caros, resulta
especialmente deseable obtener elevadas velocidades de transmisin de calor con el fin de
minimizar el costo del equipo.
Estos problemas y propiedades de los lquidos son de los ms importantes a considerar,
entre otros estn la corrosin, incrustacin y cristalizacin de solidos.

Iplwos de evmorndores v accesorios.

Objetivo.
Coriocer los difrentes tipos de evoporadores ms comunes y comerciales que
existen, as como las cnrclctersiicas principales que ellos poseen. Adems de conocer as
mismo el uso y fitnciones de los accesorios que cotifienen.

IT.1. GENERALIDADES DE LOS EVAPOKADORES.

Los componentes esenciales de un evaporador son: la cmara de calefaccin y la cmara de

evaporacin, separadas por una superficie de calefaccin. La forma y disposicin de ambas

cmaras, diseadas para lograr un funcionamiento eficaz y un valor mximo del coeficiente
de calos, varan de unos a otros tipos de evaporadores.

El evaporador ms sencillo est formado por una cmara de calefaccin (camisa de vapor)
que rodea el recipiente donde se efecta la evaporacin. La superficie de transmisin de
calor tiene aqu un Sirea muy limitada, y el di;positivo slo sirve para evaporaciones en
pequea escala. !Zn caso contrario, hemos de recurrir a la superficie de calefaccin tubular,
que permite incluir un rea de transmisin de calor muy extensa en un aparato de
dimensiones mnimas.
Los evaporadores pueden operar bien como unidades con un paso a travs o con
circulacin. En la operacin de un paso a trivs, el lquido de alimentacin pasa una sola
vez d travs de los tubos, desprende el vapor y sale de la unidad como disolucin
concentrada. Toda la evaporacin tiene lugar en un solo paso. La relacin de evaporacin a
alimentacin est limitada en una unidad de un solo paso, por tanto, estos evaporadores se:
adaptan bien a la operacin en mzltiple efecto, don de la concentracin total puede
conseguirse en varios efectos. Los evaporadores de pelcula agitada operan siempre con un
paso a travs; los evaporadores de pelcula ascendente y de pelcula descendente pueden
tambin operar de esta forma.
Los evaporadores con un solo paso a travs son especialmente tiles para materiales
sensibles al calor. Operando a vaco elevado se puede mantener baja la temperatura del
lquido. Con un solo paso rpido a travs de los tubos de la disolucin concentrada est a la
temperatura de evaporacin, pero solamente en cuanto sale del evaporador.
En los evaporadores con circulacin se ma.ntiene una masa lquida dentro del equipo. La
alimentacin que entra se mezcla con la masa del lquido y despus pasa a travs de los
tubos. El lquido no evaporado descarga de los tubos y retorna al equipo, de forma que en
cada paso solamente ocurre una parte de la evaporacin total. Todos los evaporadores de

Capitulo II. Epos de Evaporadores y accesorios.

circulacin forzada operan en esta forma; los evaporadores de pelcula ascendente son
generalmente unidades de circulacin.
L a disolucin concentrada procedente de un evaporador con circulacin se retira de la masa
del lquido. Por tanto, sta ha de tener la mxima concentracin, Puesto que el lquido que

entra en los tubos puede contener varias partes de disolucin concentrada por cada parte de
alimentacin, su concentracin, densidad, viscosidad y temperatura de ebullicin son
prximas al valor mximo. En consecuencia, el coeficietice de transmisin de calor tiende a
ser bajo.
Los evaporadores de circulacin no son muy adecuados para concentrar lquidos sensibles
al calor. Con un vaco razonablemente bueno la temperatura de la masa de lquido puede
ser no destructiva, pero el lquido est repetidamente expuesto al contacto con los tubos
calientes; por consiguiente, algo del lquido puede calentarse a una temperatura
excesivamente elevada. Aunque el tiempo de residencia del lquido en la zona de
calentamiento puede ser corto, parte del lquido est retenido en el evaporador un tienipo
considerable, y el calentamiento prolongado de aun una pequea parte de un material
sensible al calor, tal como lo son los productos alimenticios, puede atropear todo el
producto.
Sin embargo, los evaporadores de circulacion pueden operar en un elevado intervalo de
concentracin entre alimentacin y la disolucin concentrada utilizando solo una unidad,
adaptndose bien a la evaporacin de un solo efecto. Pueden operar indistintamente con
circulacibn natural, con flujo a travs de los lubos inducido por la diferencia de densidades,
o con circulacin forzada, promoviendo el flujo mediante una bomba.
II.2. TIPOS DE E VA POR4 DORES.

Los principales tipos de evaporadores tubuliires calentados con vapor de agua actualmente
utilizados son:
1. Evaporadores de tubos horizontales
2. Evaporadores de tubos verticales.
a) Flujo ascendente (pelcula ascendente).
b) Flujo descendente (pelcula descendente).
c) Circulacin forzada.
3. Evaporadores de pelcula agitada.
En el desarrollo histrico de los evaporadores de tubo han surgido varios tipos diferentes,
que siguen emplendose en la actualidad. La r z n de la existencia de los diversos tipos hay
que buscarla ms en la costumbre de cada industria que en un criterio de utilizacin
racional.

Capitulo 11. Tipos de Evopormiores y accesorios.

II.2.1.Eiwporatlor de tubos horizontnles.

El primer tipo construido, el evaporador de hibos horizontales con vapor por el interiv,
responda a la necesidad de la superficie di: calefaccin contenido el un aparato de
volumen limitado. Posteriormente, y siempre por la tendencia a aumentar la capacidad de
evaporacin, se ha tratado de mejorar la-, condiciones de transmisin de calor desde el
vapor condensante hasta el lquido en ebullicin, consiguindose a igualdad de superficie
una mayor capacidad de evaroracin. Para ello es necesario aumentar el coeficiente de
coiweccin en el lado del lquido hirviente, haciendo que sea ste, y no el vapor, el que
circule por el interior de los tubos, mientras que el vapor condensa en el exterior. Los
distintos modelos de evaporador difieren en el modo en que se ha activado esta circulacin;
si en unos de ellos es superior la capacidad dl2 evaporacin por unidad de superficie, para
lograrlo ha sido necesario introducir complicaciones mecnicas que pueden no ser
compensadas por aquel aumento.
Lus tubos horizontales interfieren con la circu'lacin natural del lquido en ebullicin y, por
tanto, minimizan la agitacin del lquido. Como resultado, el coeficiente total de
transferencia de calor es menor que en otra; formas de evaporadores, en especial si la
solucin es viscosa. No tiene aditamentos para romper la espuma que se presenta debido a
la accin de la ebullicin. Adems, la incrustacin originada por la evaporacin de la
solucin se acumula en el exterior de los tubos de donde no puede eliminarse con facilidad
como sera en el interior de los tubos.

Fig.IT. 1. Evaporador de tubos horizontales.

Capitulo II. Tipos de Evaporadoresy accesorios.

El tipo normal de evaporadores de tubos horizontales es el representado en la Fig.II.1. La

cirniii-a de evaporacin est formada por un cuerpo cilndrico vertical, cerrado en sus dos
bases, con salidas para el vapor en la parte superior y para la solucin concentrada por el
fondo. Esta ltima puede ser cnica o tener h r m a apropiada para la recogida de cristales,
cuando la evaporavijn d lugar a una cristalizacin. La parte inferior atravesada por un haz
de tubos insertados en dos chapas colectoras que forman las dos caras internas de las cajas
de vapor. El vapor entra por una trubera de conexin y se condensa a lo largo de los tubos,
arrastrando las pequeas cantidades de gases ciondensables que puede contener hacia el otro
extremo, donde se ha colocado el grifo de purga de aire. Este evaporador se construye de
chapa de acero o hierro fundido. Se emplea con preferencia para soluciones no viscosas que
no depositen cristales o impurezas en la cristalizacin, ya que la operacin de limpieza de
las superficies es relativamente dificil.

El evaporador Yaryan responde a una conisepcin bastante diferente, ya que el vapor

condensa en el exterior de los tubos. Consta de dos cuerpos cilndricos horizontales: el
primero contiene u n haz de tubos que es atravesado por la disolucin a concentrar, y el
segundo haz (:s una cmara de separacin con tabiques deflectores contra los que choca la
mezcla del lquido concentrado y el vapor que sale de los tubos. El funcionamiento de este
evaporador se explica fcilmente considerando la Fig.II.2.

Fig.II.2. Evaporador de Yaryan.

Mediante las placas de distribucin C, y Cr, se disponen los tubos en series de 3 5, que
reciben la alimentacin del colector L y expulsan la mezcla de concentrado y vapor sobre
los tabiques deflectores. El vapor entra por uno de los extremos del cuerpo de evaporacin,
y en el extremo opuesto estn la salida del ondensado y la purga de aire.
La mayor ventaja del evaporador de Yaryan es la completa destruccin de la espuma por el
choque rpido contra los tabiques deflectores. La circulacin del iquido, activada por la
formacin de burbujas de vapor, hace aumentar la eficiencia de la superficie de calefaccin.
El evaporador se desmonta fcilmente, io que facilita la limpieza de la superficie interna de
los tubos.

Capitulo II. Tipos de Evaporadores y accesorios.

II.2.2. K~lrpornrlorde tubos verticdes.

Los evaporadores de tubos verticales posteriores a los modelos ya descritos, han desplazado
a stos en muchos casos por sus ventajas indudables.
El tipo llamado Esttidar, uno de los ms conocidos, consta de un cuerpo cilndrico vertical
anlogo al descrito para el evaporador de tubos honzonta!es. En su parte inferior tiene una
cmara anular de vapor, atravesada por haces de tubos abiertos en sus dos extremos, una
mejor ilustracin est dada en la Fig.IT.3.

Fig.TI.3, Evaporador Estndar de tubos verticales.

El lquido a evaporar que entra en los haces de tubos adquiere en ellos un movimiento
ascendente al producirse las primeras burbujas de vapor. Las porciones de lquido
arrastradas del fondo han de ser reemplazadas por otras que descienden por el conducto
central B, producindose as una circulacin activa que facilita la evaporacin.
La mayor desventaja del evaporador Estndar es el arrastre de espuma o gotas del lquido
con el vapor, que slo se atena mediante la introduccin de cortaespumas y tabiques
deflectores. Este dispositivo es ms fcil de aplicar en el evaporador de cestu, representado
en la Fig.II.4. Este evaporador tiene un conducto de descenso perifrico en sustitucin del
conducto central. El elemento de calefaccin forma un cuerpo compacto, que puede sacarse
del evaporador para su limpieza. Las restantes caractersticas son anlogas a las del
evaporador Estndar, y las diferencias entre ambos son ms bien de ripo constructivo.
8

__.

Cnpitulo I,? Tips de Bwporndores.v ncceson'os.

El evaporador de irtos Icirgos iwticciles, conocido como evaporador de Kestner (Fig. 11.4.),
constituye un notable avance sobre los anteriores, ya que en l se logra mejorar el
coeficiente integral de transmisin de calor por el aumento de la velocidad de circulacin
del lquido en los tubos. El nivel del lquido se mantier:e un poco por encima del extremo
inferior. L a formacin de burbujas de vapor origina una corriente ascendente muy rpida
d \ I lquido a evaporar, que pasa una sola vez p o i los tubos. El deflector M, separa del vapor
l a mayor parte de las partculas de lquido, y el resto queda retenido en el separador de
cicln E, volviendo al tubo de salida del concentrado. Mediante deflectores situados entre
el haz de tubos se consigue que el vapor circule por 1 cmara de calefaccin de tal modo
que los gases no condensables se arrastran hacia la parte inferior, donde se dispone
fcilmente el grifo de purga correspondiente.

Fig.11.4. Evaporador de cesta y de Kestner, respectivamente.

En ambos tipos. los tubos se insertan dentro de lminas de tubos rolando o soldando, lo
cual reduce en forma considerable el costo en comparacin con los castillos empacados que
se utilizan tradicionalmente en el evaporador de tubos horizontales. Estos evaporadores
superan la mayor parte de las desventajas operacionales del evaporador de tubos
horizontales.

--

Cnpitulo II. Ts 'ins de Evaporadores y nccesorios.

Es posible manejar lquidos viscosos pero la circulacin es lenta y se obtienen cieficientes

bajos. En consecuexia, los evaporadores de tubos verticales son por completo
satisfactorios para la mayor parte de las demandas de evaporacin y slo son imprcticos
cuando el lquido que se est evaporanda es naiy viscoso, forma espuma de mmera
importante o slo puede estar sujeto a las temper )turas del evaporador durante periodos
muy cortos.
En la Fig.iI.5, se muestra un evaporador de circiikicit~fortnc.lo. En estos evaporadores, el

lquido que se est evaporando se bombea a travs de un intercambiador de calor, en donde

el medio de calentamiento rodea a los tubos que conducen la solucin. La cada de presin
y la carga hidrosttica combinadas son, con frecuencia, los sutictenteniente grandes para
evitar que la solucin hierva en los tubos del intercambiador, de manera que el vapor
generado es s e p r a d o de manera instantnea a medida que el lquido entra ai espacio de
separacin.

Fig.II.5. Diagrama de seccin transversa.1 de un evaporador de tubos verticales con

circulacin forzada.
Dado que l a velocidad de la mezcla que se separa al instante es alta, los desviadores son
importantes para minimizar el arrastre. Un desviador de diseo adecuado promover la
fusin de pequeas burbujas adems de cambiar sio la direccin de flujo. Por lo general,
los evaporadores actuales de circulacin forzada estn equipados con calentadores externos,
como el que se muestra en la Fig.11.6.
10

Copiirlo U. Tipos de Evoporodoresy occesorios.

El equipamiento externo hace que la limpieza y la sustitucin de los tubos corrodos o

erosionados sea ms simple que con el elemento interno de calentamiento. Tambin
permite la construccin de una unidad ms compacta, de manera que pueda instalarse en
espacios de poca altura. AI evaporar alcunas soluciones, es importante evitar la ebullicin
en el intc.rior de los tubos, a fin de reducir 12, deposi,.:in de slidos. En el evaporador con
calentamiento externo, es fcil evitar la ebullicin hajando el caleiitador en relacin con el
espacio de separacin. E:::o no puede evitarse tan fcilmente cuando la superficie de
intercambio de calor se enci; :ntra en el interior del cuerpo del evaporadtx.

Fig.TI.6. Diagrama de seccin transversal cie un evaporador de circulacin forzada con un

calentador horizontal externo.

En los evaporadores de circulacin forzada, el coeficiente de transferencia de calor

depender de la velocidad de circulacin, as como de la diferencia lograda de
temperaturas total, y las propiedades del sistema. Con velocidades de circulacin bajas la
ebullicin ocurre en buena medida a travs de la longitud del tubo. La ebullicin
incrementa la turbulencia y puede hacer que el coeficiente del lado de la ebullicin sea
hasta del doble de valor que tendra sin ebu'ilicin.
La fraccin del licor evaporado a su paso a travs del tubo es muy pequea; por tanto, la
circulacin total a travs de los tubos s:1 muchas veces mayor que la velocidad de
alimentacin.

Copittilo II. Tipos de Evoporodores y accesorios.

II,2.3. Eiviporntlor de pelculn cigitarla

La principal resistencia a la transmisin de calor desde el va,)or de agua que condensa hasta
el lquido que hierve en un evaporador reside del lado del lquido. Por tanto, cualquier
mtodo para disminuir esta resistencia tirodiice un considerable aumento del coeficiente
global de transmisin de calor. En evaporadores de tubos larp3s, especialmente con
circulacin forzada, la velocidad del lquido a travs de los tubos es elevada. El flujo del
lquido es altamente turbulento y la velocidad de transmisin de calor es elevada. Otra
forma de aumentar la turbulencia es mediante agitacin mecmica de la pelcula del liquido,
tal como lo muestra la Fig 11.7, que es un evaporador de pelcula descendente modificado,
con un solo tubo encamisado que contiene un agitador interno. La alimentacin entra por la
parte superi<Jrde la seccin encamisada, mientras que el vapor asciende desde la zona de
vaporizacin hasta un separador no encamisado cuyo dimetro es algo mayor que el tubo de
evaporacin. En el separador las palas del agitador proyectan hacia fuera el lquido
arrasirado que choca contra unas placas verticales estacionarias. Las gotas coalescen sobre
estas placas y retornan a la seccibn de vaporizacin. El val >r exento de lquido sale a travs
de orificios situados en la parte wperior de la unidad

Fig.II.7. Evaporador de pelcula agitada.

La principal ventaja de un evaporador de pelcula agitada es su capacidad para conseguir
elevadas velocidades de transmisin de calor con lquidos viscosos. Es particularmente
eficaz con materiales viscosos sensibles al calor tales como gelatina, ltex de caucho,
antibiticos y zumos de mitas. Sus desventajas son el elevado costo, las partes internas
mviles que pueden requerir un importante mantenimiento.

12

Capituio11. Tipos de Evaprodores y accesorios.

11.3.ACCESORIOS DE LOS ET/APOWDORES.

Describiremos algunos accesorios brevemente, que son fundamentales para el

funcionamiento de los evaporadores y que encuentran un importante lugar en el m p l e o
hacia otras operaciones tecnolgicas.
I13.1. Conrlensariores.

Cuando un evaporador 'ia de trabajar a presin i i xior a la atmosfrica, el procedimiento

ms eficaz es conectar la cniara de evaporacin ,on una bomba de vaco, a traves de un
condensador enfriado por agua. La presin sobre la cmara de evaporacin nunca puede
llegar a ser menor que la tensin del v a p x del igua a la temperatura del condensador.
Aunque pueden emplearse condensadores de superficie, casi siempre se emplean los de
contacto directo, por ser ms econmicos.

El condensador de superficie es un intercambiador de calor, y puede servir co~iiotal un

calentador tubular de vapor. Las dimensiones del calentador se calculan de acuerdo con las
leyes de trmsmisin de calor, ya conocidas.

En el condensador de contacto, la condensacin del vapor se efecta por contacto con el

agua de refrigeracin. El contacto puede hacerse eri corriente paralela o en contra corriente;
en el primer caso, los gases no condensab1i:s se enfran hasta la temperatura de salida del
agua; en el ltimo, hasta la temperatura de entrada. El agua y los gases no condensables
pueden salir del condensador juntos o separados, segn el diseo que se adopte.

Fig.II.8.Condensador en contracomente (cascada).

13

--

Cap ti110 II. Tipos de Evoporodoresy nccesorios.

En 11 Fig 11.8, se representa un condensador en contracomente, con salidas separadas del

aguz y de los gases incondensables. El agua cae .n dos cascadas, ponindose en contacto
con la mezcla gaseosa procedente del evaporador. Este condensador requiere una bomba de
vaco para extraer el aire, conectada al extremo inferior. Esta ltima puede sustituirse por el
dispositivo llamado columna baromtnca.

El sistema ms elemental para la extraccin del condensado consiste en un grifo rt gulahle

que se gi ~diiade modo que el caudal de salida sea igual a la cantidad de agua formada por
condensa. in; pero este sistema requiere de stenciii continua de un operario, y, a pesar de
ello, es fcil que se produzcan prdidas de vapor o anegamientos de la cmara de
calefaccin.
IZ. 3.2. Purgadores.

Los aparatos automticos empleados reciben el nombre d:: purgadores. Su misin es

permitir la descarga del condensado sin dejar salir el vapor. Exi::ten dos tipos principales de
purgadores: los termostticos y los mecniccls
Los purgadores termosidticos basan su distincin entre el vapor y el condensado en la
diferencia de temperaturas que tienen ambos junto a la parte activa del purgador. El
funcionamiento de uno de estos aparatos se comprende estudiando la Fig.II.9, que
corresponde a un purgador termosttico de dilatacin de lquido.

Fig.II.9. Purgador de condensado, de tipo termosttico.

En el interior del cartucho A hay un tubo de paredes onduladas a cuyo extremo (de la
izquierda) se une el vstago de la vlvula C. El espacio entre el tubo A y el fuelle est lleno
de aceite. Ai calentarse el cartucho metlico por contacto con el vapor, se dilata el aceite,
obligando a contraerse el fuelle B, y con ello la vlvula C se cierra sobre su asiento. Ai
enfriarse nuevamente el elemento se abre la vlvula, dando salida al condensado.

14

Capitulo II. Tipos de Ev'voporodores.v nccesorios

Los purgadores mechtricos se diferenci,.n de los termostticos que la diferenciacin de

temperatura se realiza mecnicamente, por medio de una cubeta o flotador. Este purgador
se ilustra en 1 Fig.IT.10.

l3g.H. 10. Purgador de condensado, de tipo mecnico

La cubeta A mantiene su posicin horizontal hasta que se llena de condensado; en este
momento gira sobre el eje de la izquierda, abriendo, al caer, la vlvula de salida C, por
desplazamiento del vstago B.Parte del condensado sale por la vlvula, y al descargarse la
cubzta adopta nuevamente la posicin horizontal.

Los purgadores mecnicos pueden emplearse para la extraccin de grandes volmenes de

condensado, y se adaptan bien a las variaciones de caudal. En muchos casos puede
acoplarse un dispositivo de rel, con un purgador pequeo dispuesto en derivacin con la
corriente principal de condensado, que acta sobre la vlvula de salida. Los purgadores
descritos descargan el condensado a una presin siempre inferior a la existente en la cmara
de condensacin. Otros diseos mecnicos ms complicados pueden efectuar la descarga
contra un espacio que se encuentre a presin a la del vapor condensante.
El vapor de agua siempre arrastra aire y otros gases no condensable, en mayor o menor
proporcin. Estos gases se acumularn en la cmara de calefaccin, haciendo disminuir
progresivamente la temperatura de condensacin y el coeficiente de transmisin. Por esto,
ha de proveerse un medio para la evacuacibn de aire. Cuando la cantidad es pequea, basta
con un pequeo grifo de purga, que se maneja a mano; pero, en instalaciones grandes, es
conveniente el empleo de dispositivos de separacin automtica del aire sin que arrastre
vapor, que, en lneas generales, estn basados en los mismos principios que el purgador de
condensado.

15

Caplitdo III

Capacidad y ecolzomia de un evaporador.

Objetivo.

Coriocer los efectos q:*eestablecen r'cr capncidad de

saber como obmer la economci del mismo.

z~iievaporador, y

de esa forma

Capitulo III

III.1. CAPACIDAD D E UN EVAPORADOR

Las principales caractersticas de funcionamiento de un evaporador tubular calentado con

vapor de ajua son la capacidad y la economa. La capacidad se define como el nmero de
libras de agua evaporada por hora.

El funcionaniiento de un evaporador est sometido a las leyes de transmisin del calor. De

acuerdo con ellas, la cantidad de calor transmitida es proporcional a tres factores:
1. Superficie de evaporacin.
2. Diferencias de temperaturas entre el vapor condensante y el lquido hirviente.
3. Y el coeficiente global de transmisin de calor.
La superficie de evaporacin tiene un valor fijo en un evaporador ya construido. No
podemos decir lo mismo de la diferencia de temperaturas y del coeficiente de transmisin,
que vara con las condiciones del vapor condensante y del lquido a evaporar.

ITA 1.1. D$?rencia de tempernturns.

Este factor es funcin directa de las temperaturas de la cmara de condensacin y la cmara
de evaporacin.
La temperatura de la cmara de condensac:in depende de la presin del vapor, que se
mantiene constante en una instalacin normal. Slo cuando se acumulan gases no
condensables puede descender la temperatura de condensacin pera una misma presin
total, ya que sta ser la suma de la pr'esin del vapor y la que ejercen los gases
permanentes.
Para una proporcin relativamente pequea tie gases puede cesar en absoluto la transmisin
de calor, por anularse la diferencia de temperaturas.

16

-.

Capirillo 111. Capacidad y economa de un evaporador.

.
I

En un aparato que est evaporando agua, la temperatura de la cmara de evaporacin es la

de ebullicin del agua a l a presin dada; determinando sta con un manmetro, hallaremos
la temperatura en las tablas de vapor. En la prctica, la disolucin que se evapora contiene
mayor o menor proporcin de slidos disueltos, y, en consecuencia, el punto de ebullicin
es superior al del agua a la misma presin.
Los evaporadores funcionan, en general, en marcha continua, y as, la concentracin del
lquido puede suponerse constante e igual a la concentracin de salida. Para los fines de la
transmisin de calor supondremos que el lquido se calienta hasta la temperatura de
ebullicin, antes de llegar a la superficie de calefaccin.

IIL 1.2. Eleiwcicin del punto de ebullicin.

El aumento en el punto de ebullicin, aunque dado por las leyes de la fsicoqumica, se

determina ordinariamente mediante la reghi de Diihring: la temperatura de ebullicin de

una disolucin acuosa a distintas presiones es funcin lineal de la temperatura de ebullicin
del agua a la misma presin. Basndonos en esta regla, si conocemos la temperatura de
ebullicin a dos diierentes presiones, podemos construir un diagrama (como el de la
Fig.ITI.1 para el hidrxido sdico) que nos permite obtener el punto de ebullicin a
cualquier presin.
La elevacin en el punto de ebullicin vara poco con la presin, y es prcticamente
constante dentro del intervalo que puede interesamos. Una vez conocido el valor de la
elevacin podemos sumarlo a la temperatura deducida de las condiciones del manmetro y
de las tablas de vapor, y as hallar la verdadera temperatura de ebullicin. Para las
disoluciones concentradas, este trmino de correccin puede ser del mismo orden de
magnitud que la diferencia real de temperaturas.
La temperatura as determinada corresponde a la ebullicin en la superficie libre del lquido
en el evaporadoi. Por debajo de este nivel, la temperatura ser superior, como corresponde
al aumento de la presin total con la carga hidrosttica del lquido. Aunque la elevacin por
este concepto es corrientemente de poca iniportancia, representa un factor digno de ser
tenido en cuenta en un evaporador que opere con diferencias de temperaturas del orden de
IOOC.
Por otra parte, al estudiar la transferencia de calor se observa cmo la existencia de cargas

hidrostticas altera el esquema ideal supuesto para la evaporacin, de modo que un parte
muy considerable de la superficie de calefaccin el proceso real puede ser el calentamiento
del lquido y no la evaporacin.
Esto explica que el coeficiente de transmisin aparente vare con el nivel del lquido en el
evaporador, presentando, a veces, un mximo muy marcado para una altura ptima.

17

Cwitulo !I. Cnwcidod Y econ m i n de un evnnoroclor.

Punto de ebullii:in del agua,

"F

Fig.III.1. Lneas de Dhring para el sistema NaOH - HzO.

Cnpitolo lil. *mopocidnd.v

econonia de un evaporotlor.

III. 1.3. Coeficiente de trtinsniisitn de (.olor.

En el coeficiente integral de transmisin se han combinado los coeficientes superficiales del

vapor condensante y del lquido hirviente con las resistencias de conduccin de la pared de
separacin y los depsitos slidos que pueden formarse en ambas caras. A excepcin de la
resistencia de la pared, y quiz tambin del coeficiente superficial de condensacin, todos
los factores que determinan el coeficiente integral son las variables. Aunque supongamos
que las superficies del evaporador se mantienen perfectamente limpias, queda como
variable el coeficiente superficial de evaporacin.

El coeficiente de tr::nsmisin vara de uno:r puntos a otros de la superficie, y el que se

emplea en los clculos es un valor medio determinado experimentalmente para la
evaporacin del lquido dado. Como este coeficiente se ha despejado de la igualdad
q

UAAt

ec. 1

su valor depende del incremento de temperatura considerado. Normalniente hemos de

tomar como tal la diferencia entre la temperatura de condensacin, deducida de la presin,
y la temperatura real de ebullicin en la superficie del lquido. No es aconsejable el empleo
del coeficiente nparenfe, en cuya deduccin se hace caso omiso de la elevacin del punto
de ebullicin, porque podra llevar a resultados totalmente falsos.

III. I.4. Coeficiente de in pelcula de vapor.

El coeficiente de la pelcula del vapor de agua es elevado, aun para la condensacin en la
pelcula. A veces se adicionan promotores al vapor de agua para provocar la condensacin

en gotas y aumentar todava ms el coeficiente.

Puesto que la presencia de gases no condensables reduce sustancialmente el coeficiente de
pelcula del vapor, es preciso tomar precauci'ones para purgar los no condencables de la caja
de vapor y prevenir la entrada de aire cuando el vapor de agua est a una presin inferior a
la atmosfkrica.

111.1.5. Coeficientes del Indo del lquido.

El coeficiente del lado del lquido depende en gran medida de la velocidad del lquido sobre
la superficie de calefaccin. En la mayora de los evaporadores, y especialmente en los que
tratan materiales viscosos, la resistencia del 'lado del lquido controla la velocidad global de
transmisin de calor hacia el lquido en ebullicin.
La formacin de costras sobre los tubos de un evaporador aade una resistencia termica
equivalente a un factor de ensuciamiento.

19

Cnpitiito III. Copncirlnc.'y economa de un evnporodor.

.

III.1.6. CocJcientes globales.

Debido a la dificultad de medir los elevados coeficientes de pelcula individu.iles en un

evapurador, los resultado experimentales generalmente se expresan en funcin de
coeficientes globales, que estn basados en la cada neta de temperatura corregida para
tener en cuenta la elevacin del U P 3 de ebullicihn Los coeficientes globales estn, por
supuesto, influenciados por los mismos factores que los coeficientes individualey pero si
una resistencia (tal como la de pelcula del lquido) es controlante, variaciones importantes
de las dems resistencias apenas tendran efecto sobre el efecto global.
I(..2. ECONOMIA DE UN EVAPORADOR

La economa es el nmero de libras vaporizadas por libra de vapor vivo que entra a la
unidad. En un evaporador de simple efecto la economa es siempre algo menor que la
unidad, pero en los evaporadores de mltiple efecto puede ser considerablemente mayor.
El principal factor que influye sobre la economa de un evaporador es el nmero de efectos.
Mediante un diseo adecuado, la entalpa de vaporizacin del vapor de agua que entra en el
primer efecto puede utilizarse una o ms veces dependiendo del nmero de efectos. La
economa tambin est influenciada por la temperatura de la alimentacin.

Si la temperatura es inferior a la de ebullicin en el primer efecto, para el calentamiento de

la carga se utiliza una parte de la entalpa de vaporizacin del vapor de agua y solamente
una parte queda disponible para la ocupaci6n.
Si la alimentacin est a una temperatura superior a la de ebullicin, la vaporizacin sbita
que se produce contribuye a generar una evaporacin adicional a la producida por la
coriderisacin del vapor de agua. Desde el punto de vista cuantitativo la economa de un
evaporador es totalmente una cuestin de balances de entalpa.

20

Evnporador de simpe efecto.

Objetivo.

Preseninr el concepto de evnporndor de simple efecto.

Capitulo IV

i K 1 . EVAPORADOR DE SMYLE EFECTO.

El clculo de un evaporador se efecta pos aplicacin de varios conceptos bsicos de la

ingeniera qumica, como son el balance de materias, el balance calorfico, y la ecuacin
para la intensidad del paso de calor.
Todos los clculos se refieren a una unidad de tiempo determinada, generalmente la hora.
El balance de materias nos permite conocer los caudales de todas las corrientes de entrada y
salida al evaporador en funcin de una de ellas y las concentraciones inicial y final.
Para conocer la cantidad de calor que ha tie pasar de una a otra cmara a trav. de la
superficie de calefaccin planteamos un balance calorfico en el que intervienen los calores
especficos de la alimentacin y del producto, el vapor latente de vaporizacin del agua, los
calores latentes de dilucin, de cristalizacin, etc. El sumado ms importante es el calor de
vaporizacin del agua, que suele tomarse a temperatura de ebullicin del agua a la misma
presin. Los datos sobre calores de dilucin y cristalizacin son poco comentes, y como su
orden de magnitud es pequeo suelen despreciarse frente al calor de vaporizacin. No
puede decirse lo mismo del calor sensible de la disolucin, que representa una fraccin
considerable del calor necesario cuando la alimentacin entra a temperatura baja. La suma
del calor sensible y el calor de dilucin se puede determinar exactamente por diferencias de
entalpas cuando se dispone de los datos entalpa concentracin. (Fig.iV. 1, para el NaOH)
De las consideraciones antes presentadas, los clculos necesarios para fijar los
requerimientos de rea de calentamiento de un evaporador simple, a partir del conocimiento
de su carga, resultan directos. Slo necesitan resolverse de manera simultnea la ecuacin
de velocidad de transferencia de calor con el balartce de calor y los balances de masa total y
soluto. Como se mencion con anterioridad!, el coeficiente global de transferencia de calor
se obtiene de la experiencia previa.
Igualando el calor total, q, ganado por la disolucin al perdido por el vapor de calefaccin
obtendremos la cantidad de vapor que ha de condensarse. La superficie de calefaccin se
despeja de la ecuacin general q = UAAt; la diferencia de temperaturas se deduce de las
condiciones fijadas.

21

Copitiilo I% Evnporodor de siiitple e j c t o .

Fig.IV. 1. Diagrama entalpa-concentracicjn para soluciones acuosas de NaOH a una

presin de 1 atm. El estado de referencia peua el agua lquida a 32OF a su presin de vapor.
Para NaOH, es estado de referencia es el de una solucin infinitamente diluida a 68F.

22

Copit )lo IV. Evoporodor de simple efecto.

--

IK2. BAIANCESDE ENTALPA EN UN EVAPORADOR DESIMPLE EFECTO.

En un evaporador de simple efecto el calo: latente de condensacin del vapor de agua es
transmitido a travs de una superficie de calefaccin para vaporizar agua de una disolucin
a ebullicin. Se necesitan dos balances de entalpa, uno para el vapor de agua y otro para el
lado de la disolucin del vapor.

La velocidad de flujo del vapor de agua y del condensado es m,, la del lquido de
alimentacin es mfi y la del lquido concentrado es m.La velocidad de flujo de vapor hacia
el condensador, suponiendo que no se precipitan slidos del lquido, es ni/- m.Por otra
parte, sea T, la temperatura de condensacin del vapor de agua, T la temperatura de
el~iillicindel lquido en el evaporador, y T3la temperatura de la alimentacin como lo
muestra la Fig.TV.2.

Fig.IV.2. Balances de materia y energa en un evaporador.

Se supone que no hay fugas o arrastre, y que no es preciso tener en cuenta las prdidas de
calor en el evaporador. La comente de vapor de agua que entra a la cmara de
condensacin puede estar sobrecalentada, y el condensado gener-lmente abandona la
cmara de condensacin algo subenfriado por debajo de su temperatura de ebullicin. Sin
embargo, tanto el sobrecalentamiento del vapor de agua como el subenfriamiento del
condensado son pequeos y resulta aceptable despreciarlos al aplicar un balance de
entalpa. Los pequeos errores que SI: cometen ai despreciarlos se campensan
aproximadamente al no tener en cuenta las perdidas de calor desde el evaporador.
Con estas suposiciones la diferencia entre la. entalpa del vapor de agua y la del condensado
es simplemente A,, el calor latente de condensacin del vapor de agua.

23

Criprtulo IV. Evoporador de simple eJicto.

.~

..

El balance de entalpa para el lado del vapor de agua es

qr = nis (H, - HJ

donde

m, A,

ec. 2

velocidad de transmisin de calor a travs de la superficie de calefaccin

desde el vapor de agua.
U; = entalpa especfica del w p o r tie agua.
H, = entalpa especfica del conderisado.
A, = calor latente de condensacin del vapor de agua.
M, = velocidad de flujo del vapor de agua.
(I. =

El balance de eiitalpa para el lado de la disolucin es

q = (mf - m) H, - m#lf
donde q

mH

ec. 3

Velocidad de transmisin de calor desde la superficie de calefaccin

hacia el lquido.
H,= entalpa especifica del vapor.
Hf = entalpa especfica de la disolucin diluida.
H = entalpa especfica de la disolucin concentrada.
En ausencia de prdidas de calor, el calor transmitido desde el vapor de agua hacia los
tubos es igual al transmitido desde los tubos hacia la disolucin y, por tanto, qr = q. As,
igualando las ecuaciones 2 y 3.
=

q=mS~3=(m~-m)H,-m/l;r/+mH

ec. 4

Las entalpas del lado de la disolucin H, H/; y H dependen de las caractersticas de la

disolucin que se concentra. La mayor parte de las disoiuciones cuando se concentran o se
diluyen a temperatura onstante no producen un g a n efecto trmico. Esto se cumple para
disoluciones orgnicas, as como para disoluciones moderadamente concentradas de la
mayor parte de las sustancias inorgnicas. As, el azcar, la sal comn y las disoluciones de
fbricas de papel no poseen calores de dilucin o de mezcla apreciables.
Por otra parte, el cido sulfrico, el hidrxido sdico y el cloruro clcico, especialmente en

disoluciones concentradas, desarrollan una cantidad de calor considerable cuando se

diluyen y, por tanto, poseen importantes calores de dilucin. Adems del calor latente de
vaporizacin, cuando las disoluciones de estas sustancias se concentran hasta densidades
elevadas se requiere una cantidad equivalente de calor.

24

CaptuLo v

Evuporudor con mltiples efectos.

Objetivo.
Presentar el concepo de mrlfipies Igectos,los mtodos de alimentacin y mmero
optimo de efectos.

Capitulo v

K 1. EVAPORADOR D E MULTIPLES EFECTOS.

El calor latente del vapor puede aprovecharse para evaporar una disolucin que hierve a
temperatura inferior a la de condensacin c!e aqul. La disolucin puede ser la misma de
que aqul que procede, si se lleva a un segundo evaporador donde se mantiene una presin
inferior a la del primero. Del mismo modo, podemos emplear el vapor formado en el
segundo evaporador para hacer hervir el lquido, a presin an mas baja, en un tercer
evaporador; y as sucesivamente.

La instalacin as formada se llama evaporacin en mltiples efectos o simplemente

miltiple efec::o, y los distintos evaporadore:; se denominan primer efecto, segundo efecto,
etc. El nmero posible de efectos tiene L a lmite en la prctica, como veremos ms
adelante. La Fig.V. 1, muestra tres evaporadares que forman un sistema de triple efecto.

m,efectos primero, segundo y tercero.

Fig.V.1. Evaporador de triple efecto: I, Ii,

Capituio

K 1.1.

t;:flp#rador con mlittiptes ejktos.

FUNCIONAMIENTO DE UN MULTIPLE EFECTO.

Sea un triple efecto, como el que se muestra en la Fig.V.l. El lquido a evaporar entra al
primer efecto, y pasa sucesivamente por 10s otros dos. Suponemos que los tres
evaporadores, en un principio a la temperatura ambiente, estn llenos hasta el nivel normal,
y que mediante un dispositivo se mantendr este nivel en ;ualquier momento.
Para iniciar el funcionamiento abrimos la llave del vapor y ponemos en marcha la bomba
de vaco, que mantiene un vaco constante sobre la cmara de ebullicin del ltimo efecto.

AI entrar el vapor en el prinier efecto desplaza el aire all contenido, que pasa por el
purgador VI a la cmara de ebullicin. Comienza la condensacin, y el calor cedido se
emplea primero en aumentar el lquido hasta la ebullicin, y despus en la produccin de
vapor. Este vapor pasa al segundo efecto, desaloja el aire y se condeisa, calentando y
evaporando el lquido.
El proceso se repite para el tercer efecto, donde la temperatura de ebullicin est fijada por
el vaco de la bomba. AI empezar a hervir el lquido disminuye la diferencia de
temperaturas entre las dos cmaras, y con ello el consumo de vapor, aumentando la presin
en el espacio de evaporacin del segundo efecto.
En consecuencia, disminuye la diferencia de: temperaturas en ste y aumente la presin en
el primer efecto. Se llega as aun equilibrio tie funcionamiento con temperaturas constantes
en cada efecto.
Las vlvulas F,,F2,y F, han ido abrindose para dar paso ai lquido que compensar la
prdida de nivel por evaporacin. Cuando el lquido alcanza en el tercer efecto la
concentracin deseada, se abre la vlvula T, y se regula la velocidad de salida de modo que
se mantenga aquella concentracin.
Para mayor simplicidad, supongamos que en la dilucin no hay aumento de la temperatura
de ebullicin respecto del agua pura, y que el vapor no pierde presin al pasar de uno a otro
efecto. Las presiones y temperaturas en cada. cmara sern las indicadas en la figura, y las
diferencias de temperatura en los tres efectos resultan ser:

La suma de las tres diferencias, que es

At! i-At-

At3 = to - f 3

= At,

equivale al incremento de temperatura en un simple efecto que trabaja con el mismo vapor
de calefaccin y la bomba de vaco conectada a su cmara de evaporacin.

26

__-."
..

Capiid o

V. Evaporador con niiiipies efectos.

K i . 2 . CAPACIDAD DE EVAPORACION Y SUPERFICIE DE CALEFACCION.

Para hacer un calculo aproximado de la capacidad de evaporacin, es costumbre despreciar

las variaciones de calor sensible del lquido en los tres efectos, as como las prdidas de
calor con el exterior.
Supongamos una instalacin formada por n efectos, todos ellos de igual superficie de
evaporacin, A. A partir de W Kg. de vapor de calefaccin, se produce en el primer efecto
una cantidad aproximadamente igual de vapor a presin p~ (ya que los calores latentes
varan poco con la temp, j.atura, y el calor de condensacin del uno se emplea en la
evaporacin del otro).
Aplicando el mismo razonamiento a los restantes efectos, deducimos que las cantidades de
calor transmitidas y las evaporaciones son iguales en todos ellos. La primera consecuencia
que obtenemos es: la evaporacin total en el mltiple efecto es n-veces superior a la
obtenida en el primer efecto. El consumo de vapor es el mismo, y la superficie de
calefaccin, n-veces superior.
Para comparar el sistema con un solo efecto que trabaje entre las temperaturas to y t,,,
tendremos en cuenta que, al ser mayor la diferencia de temperaturas, tambin los es, y en el
mismo grado, la cantidad de calor que puede transmitirse. Una serie de suposiciones
aproximadas nos permite deducir que aquel incremento es n veces superior al de cada uno
de los efectos. En consecuencia, podemos afirmar que: la evaporacin total en el mltiple
efecto es igual a la obtenida en un solo efecto que funcione con la diferencia de
temperaturas lo - t,,. El consumo de vapor es n-veces menor, y la superficie de calefaccin
n-veces mayor.
Resumiendo las anteriores conclusiones, que son slo aproximadas, establecemos que si,
para efectuar una determinada evaporacin, sustituimos el evaporador sencillo por un
sistema de n efectos, dividiremos por el n el consumo de vapor, y multiplicaremos p o r n la
superficie de evaporacin necesaria.

V.1.3. METODOS DE ALIMENTACION.

El mtodo habitual de alimentar un evaporador de mltiple efecto consiste en introducir
mediante una bomba la disolucin diluida en el primer efecto y hacerla circular despus a
travs de los dems efectos, tal como se muestra en la Fig.V.i.a. Este mtodo recibe el
nombre de alimentacin directa. La concentracin de la disolucin aumenta desde el
primer efecto hasta el ltimo. Este modelo de flujo del lquido es el ms sencillo. Requiere
una bomba para introducir la alimentacin en el primer efecto, ya que con frecuencia este
efecto esta a una presin superior a la atmosfrica, y una bomba para extraer la disolucin
concentrada del ltimo efecto. Sin embargo, el paso de un efecto a otro se realiza sin
bombas puesto que el flujo es en el sentido de presiones decrecientes, y todo lo que se
requiere son vlvulas de control en las lneas de unin.

27

Copitulo K Evriporodor con nitittipies efectos.

Otro mtodo comn es el de nlimetrlncitr i17versc1,en el que la alimentacin diluida se

introduce en el ltimo efecto y se bombea despus a travs de los sucesivos efectos hasta el
primero, tal como se muestra en la Fig.V.2.b. este mtodo requiere una bomba entre cada
pareja de efectos, adems de la bomba de disolucin concentrada, ya que el flujo tiene lugar
en el sen;ido de presiones crecientes. La alimentacin inversa conduce con frecuencia a una
mayor capacidad que la alimentacin directa cuando la disolucin es viscosa, pero puede
producir una menor economa que la alimentacin directa cuando la alimentacin est fra.

,. . ..

,---.

Vapw d
agua

,
ondwih
conccntradi

AiimrhtmiOn

(a)

Cona
sad

Otsaluci6n
concentrada

- r;--

['.

Di'

Aliment acin

Aiimen-

N&.

Fig.V.2. Modelos de flujo del Iquido en evaporadores de mltiple efecto: (a) alimentacin
directa; (b) alimentacin inversa; (c) alimentacin mixta; (d) alimentacin paralela.

28

Capiiulo Y Evaporodor con nritiples efectos.

A veces se utilizan otros mtodos de alimentacin. En la alimenlacin mixfa la disolucin

diluida entra en un efecto intermedio, circula con alimentacin direi:ta hasta el extremo de
la serie, y despus se bombea hacia atrR; a los primeros efectos para conseguir la
concentracin final, tal como se muestra ?n la Fig.V.2.c. Esta forma de operar elimina
algunas de las bombas que se requieren eu la alimentacin inversa y permite realizar la
evaporacin final a temperatura ms elevada. En los evaporadrires con cristalizacin se
introduce directamente en cada efecto para dar lugar a lo que se llama alimentacin
paralela, tal como se representa en la Fig.V.2.d. En la alimentacin paralela no hay
transporte de lquido entre efectos.

V.2. CAPAUDAD Y ECONOMIA DE EVAPORADORES DE MULTIPLE EFECTO.

El aumento de economa mediante la evaporacin en mltiple efecto se consigue a costa de

un; capacidad reducida. Es claro que aumeritarido varias veces la superficie de calefaccin
auni-iitar la capacidad de evaporacin, pera ste no es el caso.
La capacidad total de un evaporador de mliiple efecto generalmente no es superior a la de
uno de simple efecto que tiene igual superficie de calefaccin en cada uno de los efectos y
opera con las mismas condiciones extremas, pero cuando hay u n apreciable aumento del
punto de ebullicin dicha capacidad con frecuencia es considerablemente menor.
Cuando la elevacin del punto de ebullicin es despreciable el ATefectivo global es igual a
la suma de los AT de cada efecto, y la cantidad de agua evaporada por unidad de Oren de
superficie en un evaporador d e N efectos es aproximadamente igual a (IIN) de la de un
evaporador de simple efecto.
Sea cual h e r e el nmero de efectos, si los coeficientes globales son los misnios, la
capacidad no ser superior a la de un simple efecto que tenga igual iiea que cada efecto de
un sistema mltiple. La elevacin del punto de ebullicin tiende a que la capacidad de un
evaporador de mltiple efecto sea menor que la del correspondiente a uno de efecto simple.
Sin embargo se modificarn los coeficientes globales de un evaporador de simple efecto.

K3. CALCULO D E UNMLTLTIPLE EFECTO.

En el diseo de un evaporador de rnliltiple efixto los resultados generalmente deseados son

la cantidad de vapor de agua consumido, el rea de la superficie de calentamiento que se

requiere, las temperaturas aproximadas en los distintos efectos y la cantidad de vapor que
abandona el ltimo efecto.

Lo mismo que en un evaporador de simple efecto, estos valores se obtienen a partir de

balances de materia, balances de entalpa, y la ecuacin de capacidad. Sin embargo, en un
evaporador de mltiple efecto se utiliza un mtodo de tanteo en vez de una solucin
algebraica directa.

29

Capitulo I< Evnporndor con niltipies efectos.

Consideremos, por ejemplo, un evaporador de triple efecto. Se pueden escribir siete

eciiaciones: un balance de entalpa para ca' a efecto, una ecuacin de capacidad para cada
efecto y la evaporacin total, que es conocida, o bien la diferencia entre los flujos de las
disoluciones concentrada diluida.
Si se supone que la superficie de calentarniento en cada efecto es la misma, hay siete
incgnitas en estas ecuaciones: (1) el flujo de vapor vivo que entra en el primer efecto, ( 2 ) a
(4) la velocidad de flujo de cada efecto, (5) l a temperatura de ebullicin en el primer efecto,
(6; la temperatura de ebullicin en el segundo efecto, y (7) la superficie de calefaccin por
efecto. Es posible obtener las siete incgnitas a partir de estas ecuaciones, pero el mtodo
resulta tedioso.
Otro mtodo de clculo es el siguiente:

I . Se suponen temperaturas de ebullicin en el primer y segundo efecto.

2. A partir de balances de entalpa se obtienen las velocidades de flujo de vapor de agua de
disolucin de un efecto a otro.
3. Se calcula la superficie de calefaccin que se requiere en cada efecto a partir de las
ecuaciones de capacidad.
4. Si las reas de calefaccin as encontradas no son aproximadamente iguales, se estiman
nuevos valores para las ternperaturas de ebullicin y se repiten las partes 2 y 3 hasta que
las superficies de calefaccin sean iguales.

En la prctica estos clculos se realizan con un ordenador digital.

V: 4. NUMERO OPTIRfO BE EFECTOS.
El coste de cada efecto, por unidad de superficie, es una funcihn de su rea total y
disminuye con el rea, tendiendo a un valor asinttico para grandes instalaciones. Por tanto,
a inversin necesaria para la adquisicin de un evaporador de N efectos es
aproximadamente N veces la de un evaporador de simple efecto de la misma capacidad. El
nmero ptimo de efectos se obtiene a partir de un balance econmico teniendo en cuenta
el ahorro de vapor de agua que se obtiene con la operacin de mltiple efecto y la inversin
adicional que de requiere.

30

Elevacin de la ianperatura de ebullicin.

Objetivo.
Dnr n conocer los efectos pie cmsa /a carga en el lquido dentro de los
evnporndores.

Vl EFECTO DE LA CARGA DE LIQUIDO Y B E LA ELEVACION DEL PUNTO DE

El3 ULLICION.
La carga de lquido y la elevacin del punto de ebullicin influyen sobre la capacidad de un
evaporador de mltiple efecto todava ms que en el caso de un simple efecto. La reduccin
de capacidad oricinada por la carga del lquido, lo mismo que antes, puede estimarse
cuantitativamente. Comll se recordar, la carga de lquido reduce la cada de temperatura
disponible en el evaporador. La elevacin del punto de ebullicin tambin reduce la cada
de temperatura disponible en cada efecto en la forma que se expone a continuacin.
Consideremos un evaporador que est concentrando una disolucin con una gran elevacin
del punto de ebullicin. El vapor procedente de esta disolucin a ebullicin est a la
temperatura de la disolucin y, por tanto, ast sobrecalentado en una cantidad igual a la
elevacin del punto de ebullicin.
Por tanto, la cada de temperatura en cada efecto se calcula a partir de la temperatura del

vapor saturado a la presin de la cmara de condensacin y no a partir de la temperatura de

ebullicin del lquido en el evaporador precederite. Esto quiere decir que de la cada total de
temperatura disponible se pierde la elevacin del punto de ebullicin en cada efecto, de
forma que la prdida de capacidad que resulta es con frecuencia importante.
En la Fig.Vi. 1, se muestra la influencia de estas prdidas en la cada de temperatura para un
evaporador de mltiple efecto. Los tres diagramas se esta figura representan la cadas de
temperatura en un evaporador de simple, doble y triple efecto.
Los tres diagramas de esta figura representan las cadas de temperatura en un evaporador de

simple, doble y triple efecto. Las condiciones extremas son las mismas en los tres casos; es
decir, la presin del vapor de agua en ei primer efecto y la temperatura de saturacion del
vapor que sale del ltimo efecto son idriticas en los tres evaporadores. Cada efecto
contiene un lquido con elevacin del punto de ebullicin. La altura total de cada columna
representa el intervalo de temperatura desde la del vapor del vapor vivo hasta la de
saturacin del vapor procedente del ltimo efecto.

31

Cnpilrtlo 1 4 Elevccin de la eniperniuro de eb ullicidn.

:VAPOR DE AGUA AL PRIMER EFECTO

CAIDA DE

L.__l TEMPERAlURA
ELEVACIONES D E L PUNTO

DE E6ULL:CION

TEMPERATURA DE SATURACION DEL VAPOR

SIMPLE

DOBLE

PROCEDENTE DEL ULTIMO EFECTO

Fig.Vi.1. Efecto de la elevacin de la temperatura de ebullicin sobre !a capacidad de

evaporadores.
La Fig.VI.1, muestra que en los casos extremos de un gran nmero de efectos o elevaciones
del punto de ebullicin e r u:' evaporador propuesto podran ser mayores que la cada total
de temperatura disponible. La operacin en estas condiciones es imposible. En este caso
sera preciso revisar las condiciones de operacin del evaporador para reducir el nmero de
efectos o aumentar la cada total de temperatura.
La economa de un evaoorador de efecto mltiple no est influenciada por las elevaciones
del punto de ebullicin si se desprecian factores menores, tales como la temperatura de la
alimentacin y variaciones de los calores de vaporizacin. Un kilogramo de vapor de agua
que condensa en el primer ::fecto genera aproximadamente un kilogramo de vapor, que a su
vez condensa en el segundo efecto generando un kilogramo de agua, y as sucesivamente.
L a economa de un evaporador de mltiple efecto depende de consideraciones del balance
de calor y no de la velocidad de transmisin de calor. La capacidad, por el contrario, se
reduce a causa de la elevacin del punto de ebullicin. La capacidad de un evaporador de
doble efecto que concentra una disolucin 'ron elevacin en el punto de ebullicin es
generalmente menor que la mitad de la capacidad de dos efectos simples que operan con la
misma cada global de temperatura. La capacidad de un triple efecto tambin es
generalmente menor de un tercio que tres efectos simples con las mismas temperaturas
extremas.

32

Recompresin de vapor.

Objetivo.
Implementar y dar a' conocer otros mlodos para incrementar la economa de un
evnporador.

Capitulo WI

VIL I. RECOMPRESION DEI. VAPOR.

Otro mtodo para incrementar la economa del vapor es mediante recompresin de los
vapores que salen del evaporador y su alimentacin a la chaqueta del vapor de un efecto
apropiado. Se puede usar compresin mecnica o inyeccin con chorros de vapor. En el
caso de la inyeccin del vapor, sc utiliza un eyector de vapor de alta presin, para arrastrar
y comprimir una gran parte de los vapores que salen del evaporador. Una parte del vapor se
condensa por separado para compensar el vapor de agua aadido.
La energa del vapor formado en la ebullicin de una disolucin puede utilizarse para
vaporizar ms agua siempre que exista una cada de temperatura en la dircccin deseada
para la transmisin de calor. En un evaporador de mltiple efecto esta cada de temperatura
se crea debido a una progresiva disminucin del punto de ebullicin de la disolucin en una
serie de evaporadores que operan con presiones absolutas que van disminuyendo. La
deseada fuerza impulsora se puede obtener tambin aumentando la presin (y, por
consiguiente, la temperatura de condensacin) del vapor formado, utilizando una
recompresibn mecnica o trmica.

VII.1.1. RECOiWPRESXON MECANICA.

En la Fig.VII.1, se ilustra el fundamento de la recompresin mecnica del vapor. La

alimentacin fia se precalienta hasta una temperatura prxima a la de ebullicin mediante

intercambio de calor con la disolucin concmtrada y se bombea a travs de un calentador
como en el caso de un evaporador convencional de conveccin forzada. Sin embargo, el
vapor formado no se condensa directamente sino que se comprime hasta una presin algo
ms elevada por medio de un compresclr centrfugo o de desplazamiento positivo,
transformndose as en vapor vivo que entra como alimentacin en el calentador.
Puesto que la temperatura de saturacin del vapor comprimido es superior a la de ebullicin
de la alimentacin, el calor fluye desde el vapor hacia la disolucin generando ms vapor.
Se requiere una pequea reposicin d e vapor de agua.

33

Capitulo VIL Reconipresih del iwpor.

La cada de temperatura ptima para un sistema tpico es del orden de 10F. La utilizacin
de energa en un sistema de este tipo es muy buena: la economa, basada en el vapor de
agua equivalente a la energa que se requiere para accionar el compresor es ia
correspondiente a un evaporador de 1O a 15 efectos.

Fig.VII. 1. Diagrama esquemtico de un evaporador con recompresin mecnica.

Las aplicaciones ms importantes d e la eva.poracin con recompresin mecnica son la

concentracin de disoluciones radiactivas muy diluidas y la produccin de agua destilada.

VI11.2. RECOMPRESION T E M I C A .
En un sistema de recompresin trmica el vapor se recomprime utilizando vapor de alta
presin en un eyector, tal como lo muestra la Fig.VII.2.
Esto da lugar a ms vapor del que se requiere para la ebullicin de la disolucin, de forma

que el exceso de vapor se purga o se condensa. La relacin entre el vapor de agua motriz y
el vapor procedente de la disolucin depende de la presin de evapoiacin; para muchas
operaciones a baja temperatura, con vapor de agua de 8 a 10 atm de presin, la relacin
entre el vapor de agua que se requiere y la masa de agua evaporada es del orden de 0.5.

34

Capiiulo VI[. Recompresin del vapor.

Fig.ViI.2. Diasrama de flujo para un sistema evaporador Buflovak de pelcula turbulenta

ascedente-descendente, de doble efecto con recompresin trmica.

Puesto que los chorros de vapor pueden tratar grandes volmenes de vapor de baja
densidad, la recompresin trmica resulta ms adecuada que la recompresin mecnica para
la evaporacin a vaco. Los eyectores son m i s baratos y de ms fcil mantenimiento que los
compresores y las soplantes. Las principales desventajas de la recompresin trmica son la
baja eficiencia mecnica y la falta de flexibilidad del sistema frente a variaciones de las
condiciones de operacin.
En los Estados Unidos, estos sistemas se utilizan slo en aplicaciones especiales, debido a

diversas desventajas econmicas y mecnicas

1. El compresor es costoso y tiene altos cnstos de mantenimiento, mayores que el resto
del sistema de evaporacin.
2. En intervalos de operacin razonables, la diferencia de temperaturas obtenida es
pequea, de alrededor de 10F. Como resultado, el equipo evaporador debe ser grande
para cualquier velocidad de produccin razonable. Se podran obtener diferencias de

35

Cnpitirlo WI. Reconipresin del iwpor.

temperaturas mayores con una relacin de compresin ms alta, pero la eficiencia

disminuira proporcionalmente.
3. Un evaporador con recompresin no puede arrancarse sin vapor o suministro de calor
auxiliar.
4. Si la solucin tiene un IPE grande, el costo de recompresin aumenta con rapidez
debido a que el vapor se debe comprimir a una presi tal que su temperatura de
saturacin est por encima de la temperatura a la cual hierve la solucin.

En una situacin en donde los costos de combustible son altos, pero la energ'i elctrica es
barata, la evaporacin con recompresin resulta econmica y se recomienda utilizarla. En
los Estados Unidos, los costos de combustible han sido bajos en relacin con los costos de
energa elctrica, por tanto, la evaporacin con recompresin por lo general no es atractiva.
En Europa, el combustible es relativamente caro y escaso, mientras que la energa por
abastecimiento de hidroelctricas es barata en algunas reas. Como resultado, existen varias
unidades grandes de evaporacin con recompresin operando en plantas qumicas europeas.

Metodologa de simulacin.

Objetivo.

Organizar y ndecirar el marco terico implementndolo a un cdigo de simzilacin

y resolver alglinos ejemplos tpicos.

Capitulo U I I

VIXI.I. INTROD UCCIOh!

La simulacin dinmica (de trenes de evaporadores) es una buena opcin en la

comprensin de las variaciones de un sistema por mnimas que estas sean.
A veces no es posible obtener la solucin de una ecuacin (sea esta una ecuacin

diferencial), pero s se puede encontrar una satisfactoria aproximacin. Estas

aproximaciones se hallan usando mtodos numricos.
Existen dos factores que contribuyen a la importancia de los mtodos de aproximacin
numricos. El ]rimero es el hecho de que mucnos de los conceptos matemticos que se
presentan en las aplicaciones no pueden ser calculados mediante mtodos exactos. El
segundo es la invencin de las computadoras y calculadoras de alta velocidad que han
convertido en prcticos los mtodos numricos. En realidad es ms fcil calcular algo en
forma aproximada usando una calculadora (y obtener la respuesta con la exactitud deseada)
que usar mtodos exactos, an cuando stos estn disponibles.
Luego entonces, necesitamos la herramienta que nos ayude a resolver problemas de las
diversas industrias y, en general, de la fsica aplicada, que contienen dos o ms variables
independientes. Generalmente estas variables son: en tiempo, distancia, temperatura, flujo,
presin, etc. Por tanto no podemos prescindir de algoritmos computacionales para obtener
soluciones de forma vertiginosamente rpida.

MII.2.METODO DE NEWTON
Hay u n mtodo general muy conocido para resolver problemas de ingeniera, algunas
personas lo conocen como Saliendo deipaso y otras ms lo conocen como Ensayoy
error .
Para la simulacin dinmica de este proyecto se usa ese mtodo (ensayo y error), se utiliza
el mtodo de solucin numrica de Newton para varias variables no lineales.

77

Capi(u1o Vll Metodologa de siniiilocrn.

Consideraciones generales en matemticas y en ingeniera toman en cuenta la necesidad de

encontrar las races (soluciones) de una ecuztcinf(x) = O. Sif(x) es un polinomio y de ste
necesitamos las races para manejar y entender el sistema, el mejor y eficaz mtodo de
solucin es el de Newton.
En especfico ste mtodo expone como c.slcular aproximadamente un cero real de una
funcin derivable, es decir, un nmero real x tal quef,(x) = O.
Las soluciones de una ecuacin de la forma

se llaman races de la ecuacin o ceros de la ncinf: S i f e s una funcin polinomial de

grado menor que cinco, existen frmulas para obtlxqer los ceros Se conocen ya bien las
frmulas si f tiene grado uno (funcin lineal) o dos (funcin cuadrtica). Para una funcin
polinomial de grado tres o cuatro, el metodo genernl de obtencin de los ceros es
complicado. Adems, ->aralos ceros de una funcin polinomial de grado cinso o mayor
existe un teorema, el c d se acredita <t Niels Abel, que expresa que no puede haber frmula
general en trminos de un nmero finito de operaciones sobre los coeficientes. Sin
embargo, existen procesos numricos para aproximar soluciones de ecuaciones y ahora son
ms importantes que antes debido a uso creciente de cornputadoras prograrnables.

Uno de los procesos implica una aplicacin (de la derivada y fue desarrollado por Sir Isaac
Newton en el siglo'CvI1. A ste se le conoce como el mtodo de Newton.
Consideramos todava el problema de resolvlrr la ecuacin f(x) = O para obtener una raz r.
Supngase que f es derivable, de modo que la grfica de y = f(x) tiene una tangente para
cada punto. Si podemos encontrar una primera aproximacin XI para r mediante la grfica o
algn otro recurso, deber haber una mejor aproximacin x2, en la interseccin con el eje de
las x de la tangente en el punto (XI,fix,)). Usando x2 como nueva aproximacin se obtendr
otra an mejor, x3, etctera.

En forma ms general, expone un algoritmo llamado frmula de recurrencia, que va

haciendo aproximaciones numricas hasta encontrar la solucin a cero, que posteriormente
servir para realizar el cdigo de simulacin. La extensin del mtodo de Newton-Raphson
a funciones de varias variables se denomina mtodo de Newton, Este mtodo se deriva de
una nianera anloga a la demostracin anterior.
Asuma que se desea encontrar los valores de ;c y y que hacen simultneamente
y 5 (x, y) = O. Si cada funcin se representa por una expansin de la serie de
Taylor alrededor del conjunto de valores (xk, y$, se obtiene el siguiente resultado cuando se
desprecian las derivadas de orden mayor a 1.

I; (x,y) = O

38

Copituio VIIi. bletodolonia de sinrulncidn.

Para el conjunto de valores & + I , yk+J para los cuales

Las ecuaciones VIII.1 y VII1.2 se reducen a

ec.VIiI.3

ec.VI11.4

Tambin, se sobrentiende que las funciones y sus derivadas son evaluadas a x = x k y y

El problema se reduce entonces a dos ecuaciones con dos incgnitas, d.Ek+l y dyk+i.

= yk.

Los valores calculados X k + / y y k + / del ensaya k se utilizan como valores a asumir para el
ensayo @c + I). Este procedimiento se repite hasta que se alcance a exactitud deseada de
las laces.
Para el caso general de n ecuaciones independientes con n incgnitas, se obtienen n
ecuaciones de la forma de las ecuaciones VI11.3 y VnI.4 al aplicar el mtodo de Newton.
Las ecuaciones obtenidas pueden representarse por la siguiente ecuacin matricial.
Jk d x k =

-fk

Donde J k es la matriz cuadrada jacobiana de orden n y

ec.ViII.5

d x k

y f k son los vectores columna

Captulo VIl. blerodologa de sinririncin.

En Is aplicacin de ste mtodo, se recomiend!a que las caractersticas de convergencia sean

probadas resolviendo una gran variedad de ejemplos. Tambin deben ser investigados
diferentes grupos de valores iniciales.
Si nicamente se desea encontrar las races positivas de las funciones, se deben tener
mtodos alternos para la seleccin de variables a utilizar en el ensayo siguiente cuando uno
o ms valores negativos aparecen en los clculos de los ensayos intermedios.
A continuacin se presenta un conjunto de condiciones suficientes para la convergencia del
mtodo de Newton.

Considere el vector columna a, como la solucin de las variables independientes,

a = [alCLZ ... a,lT
Como su nombre lo implica para este conjunto de valores significa que

fi(a)=O
O

f(a)= o

(1 < i < n )

ec.VIiI.6

Asuma que los elementos de f y J son continuos para todos los vectores iniciales X,-, y para
todos los vectores Xk generados al utilizar el mtodo de Newton.
Adems asuma que el determinante de J es 1dift:rente de cero (14 to) para todo XOy Xr.
Entonces se puede demostrar de que para cualquier vector inicial XO,suficientemente
cercano al vector solucin a, el procedimiento de Newton converge al vector solucin a.
vIII.3. IMPLEMENTACIONAL CODICO DE SIMULACION.

El problema bsico de impternentar un cdigo de simulacin difiere uno de otro en lo

consistente del marco terico, es decir, mientras mayor informacin de alimentacin se
tenga mayor ser la informacin que arrojar.
En nuestro caso tomamos fragmentos de un algoritmo presentado en el libro de Recetas
mimricas y codificado en un compilador de turbo C U , adaptado a nuestras necesidades.

Este algoritmo esta definido con anterioridad, y es el nico algoritmo que puede resolver
ecuaciones simultaneas de orden superior no llneales con la mayor exactitud requerida.

in

Copitulo iWi.

hetodoloain de sinrulncin.

Para lograr el objetivo de la simulacin dinmica de trenes de evaporadores se hace

necesario una presentacin adecuada dado qiJe este cdigo no solo resuelve las ecuaciones
que se requieren para un evaporador de simple efecto, sino que tambin resuelve de dos y
tres efectos con alimentacin directa, paralelo y mixta.
No difiere la simulacin de un sistema a otro, bsicamente es capturar los datos de
alimentacin al sistema; en realidad la nica diferencia es el mtodo de alimentacin del
sistema a simular, pero esto esta contemplado dentro del cdigo !e si1:iulacin de tal
manera que para el usuario no tenga mayor problema en realizar la simulacin
No es dificil implementar el marco terico a cdi(.;o de simulacin, la nica dificultad de
realizar una siniulacin dinmica radica en como pt esentar los resultados de tal manera que
sean comprensibles y de utilidad para los fines requeridos (estos son: observar las
variaciones del sistema).
La simplicidad que implica usar el programa de simulacin no implica de igual manera l a
simplicidad de elaboracin del cdigo de simulacin.
El algoritmo elaborado de manera estructurada y elaborado en un compilador de Turbo
C t t es el fin cent 11 de este texto, y de wanera que el lector pueda aplicar y usar con

amplitud el material didctico se ha logrado implementar en el cdigo de simulacin una

presentacin en forma de tablas, ya sea para capturar la entrada de las condiciones del
sistema, o bien, para el anlisis de los resultados obtenidos.
Antes de interaccionar con el cdigo, se recomienda leer y analizar el algoritmo de solucin
pwsentado en el anexo A.
Este anexo lo hemos implementado con la inica finalidad de que el lector interesado en
interaccionar con el programa simulador logre tener la facilidad de codificar su propio
programa a partir del algoritmo de solucin (en el caso de no contar con l).
Para finalizar, no debemos olvidar que la simulacin da facilidad y rapidez para encontrar
ia solucin a un problema dado.
Pero la simulacin depende mucho de la informacin de alimentacin del sistema a
simular.
Como ayuda a lograr una mayor comodidad hemos anexado tab!as termodinmicas de
vapor de agua saturado con entradas por temperaturas y por presiones, as como tablas
termodinmicas del vapor de agua recalentado y el diagrama de Miillier para entalpia y
entropa que sern de gran utilidad al momento de realizar la simulacin (stas ltimas se
presentan en el anexo B).

41

Copiulo UII. ii,[eiolio/o.en ie sinrirlncin.

VIII.4. API,ICACION Y EJEMPLOS.

En esta seccin comenzaremos a introducirlos directamente a la aplicacin del cdigo de

simulacin y operacin del programa.

Este programa puede ser utilizado bsjo ambiente Windows 3, 3.11 (trabajo en grupos) o 95,
as como de igual manera, sino se cuenta coi1 la ambientacin adecuada tambin podr ser
utilizado por uciiarios que solo cuenten ccm el sistema operativo bsico (D.O.S. Disk
Operative System). Figura VIII. 1.

Fig. VIII. 1. Apertura del cdigo de simulaciijn.

Como hemos mencionado, la ventaja de urilizar un sirnulador se basa principalmente en
que podemos monitorear las variaciones de un sistema determinado, de esta forma
proseguiremos analizando los trenes de evaporacin.
En lo subsecuente se analizarn dos sistemas propios de los trenes de evaporacin haciendo

una comparacin de entre resolverlo mediante mtodos algebraicos o resolverlo con el

programa simulador.
Estos dos sistemas que analizaremos son: un evaporador de simple efecto y un tren de

evaporacin con alimentacin en paralelo.

Procedamos entonces a plantear algunos problemas.
A?

-~

Cnptulo VIii. Aletodo sigin de sitrrulocin.

Supongamos que deseamos ver como se comportan los sistemas de evapor:idor de simple
efecto y triple efecto con alimentacin en paralelo respectivamente al ejemplo 1 y 2.
Planteado algebraicamente necesitamos solamente la formulacin correspondiente a cada
uno de los sistemas, mientras tanto en el program?, de simulacin nos proporcionan un
men de inicio en el cual podemos elegir de entre los siguientes sistemas a simular segn
muestra la figura VITi.2.

Fig. VIII.2. Men de inicio para simulaciri del tren de evaporizacin.

Una vez que elijas la opcin a simular solamente lo que resta es capturar los datos de
entrada y de ser necesario segn el caso tambin se tendrn que calcular los datos de salida.
Comencemos a resolver el ejemplo 1 de manera algebraica y posteriormente de manera
numrica mediante el programa de simulacin.
Ejemplo 1:

Un evaporador de simple efecto ha de concentrar 20,000 k g h de una disoluci6n de NaOH

al 20 % hasta un 50 YOde slidos. L a presi.n manomtrica del vapor de agua ser 20 atm y
la presin absoluta en el espacio de vapor 100 mm Hg. El coeficiente global se estima que
ser 250 W/mzh0C. La temperatura de la ;alimentacin es de 100 "C. Calclese la cantidad
de vapor de agua consumido. La economa y la superficie de calefaccin que se requiere.

C:r,pfuloViIi. ibetodologo de siniulocin.

Evaporador de simple efecto.

Evaporado E

--+

4.

Conden. V

Balance de energa en la cmara de calentamiento.

VHv q = Vhv
q = V (hv -Hv)
ec. 1 q = V hcond.
Sabemos tambin que
q=UAAT
ec.2 A = q / U A T ,
VHV+ Fhr= Chc + Vhv i- EHE
V (Hv- hv) = Chc + EHE - Fhf

ec.3 V = Chc + EHg,Fhf

hcond.
Balance global.
ec.4 F = C + E
Balance de soluto.
ec.5 FXf = CX, + EXE
Nota: XE = O por ser un soluto no volatl.
Diferencia de temperatura.
ec.6 AT = Tvap - Teb

44

--

Coptulo Vlll. Metodologa de siniulocin.

Solirciri.

La cantidad de agua se obtiene a partir de un balance de materia. La alimentacin contiene

SO120 = 4 kg de agua por kg de slido; la disolucin concentrada contiene 50/50kg de agua
por kg de slido. La cantidad evaporada es 4 .-I = 3 kg de agua por kg de slido, o sea
3 x 20,000 x 0.20 = 12,000 k g h
La velocidad de flujo de disolucin concentrada m es 20,000 - 12,000 = 8,000 kg/h.
Co?isiimode vapor de opa.Puesto que con disoluciones concentradas de NaOH el calor de
dilucin no es despreciable, la velocidad de transmisin de calor se obtiene a partir de la
ec.VI.4 y la fig. IV.l. La temperatura de vaporizacin de la disolucin de 50 % a una
presin de 100 mm Hg se obtiene como sigue, a partir de la ec.6.

Temperatura de ebullicin del agua a 100 mm Hg = 124 "C

Temperatura de ebullicin de la disolucin = 197 "C
Elevacin del punto de ebilllicin = 197 - 124 = 73 "C
Las entalpias de la alimentacin y de la disolucin concentrada se obtienen a partir de la
fig. IV. 1.
Alimentacin, 20 % de slidos, 1 O0 "C:
H/= 55 kcaiJkg
Disolucin concentrada, 50 % de sidos, 197 "C: H = 221 kcakg
La entalpa del vapor que abandona el evapclrador se obtiene a partir de las tablas del vapor
de agua. La entalpa del vapor de agua sobrecalentado a 197 "C y 1.93 atrn es 1149 kcalkg,
que es H,en la ec. 4 del captulo IV.
El calor de vaporizacin del vapor de agu.; ;1, a una presi manomtrica de 20 atm, de
acuerdo a tablas es de 939 kcakg.
La velocidad de transmisin de calor y el consumo de vapor de agua pueden obtenerse
ahora a partir de las ecs. 1 y 3.

+ 8,000 :< 221 - 20,000 x 55 = 14,456,000 kcaVh

m, = 14,456,000 1939 = 15,400 kg/h

q = (20,000 - 8,000)(1149)

Economa. La economa es 12,000 I 15,400 = 0.78.

Szcpe~ciede cnlefaccion. La temperatura de condensacin del vapor de agua es 259 "C.El
rea de calefaccin requerida es, ec.2:
A = 14,456,000 / 2513(259 - 197) = 930 mz

Este tratamiento se realizo de manera algebraica y enseguida realizaremos el mismo

ejemplo pero con el programa de simulacin y compararemos los resultados.

45

Conitulo TII. h fetodolowa de siitrulocin.

Fig. VIII.3.Condiciones bajo la simulacin de un simple efecto.

Como se comento con anterioridad, l o nico que requisita el programa es que el usuario
capture los datos de entrada.
En la anterior figura VTII.3, podemos observar las condiciones en las cuales se va a simular
el sistema de un evaporador de simple efecto. Para el caso de un evaporador de simple
efecto no se realiza ninguna consideracin como se harn para dos y tres efectos; esto es
porque realmente lo que se simula en ei simple efecto son las condiciones reales de
operacin.

En la siguiente figura VII1.4, se muestran los resultados de la simulacin.

Comparemos los resultados algebraicos y los obtenidos mediante la simulaciri y
observemos que no difieren grandemente Lino de otro, lo que significa que no importa cual
mtodo se utilice siempre obtendremos el mismo resultado, con la nica diferencia en
tiempo de los mtodos utilizados.
Tambin cabe sealar que la figura VII1.4,tiene la opcin de regresar y editar nuevamente
las condiciones de entrada lo que nos facilita demasiado su entendimiento y nos ayuda a

comprender como se comporta el sistema con la mnima variacin que se realice en las
condiciones de entrada.

46

Fig . VIII.4. Resultados de la simulacin de cin simple efecto

Como se comento con anterioridad realizaremos dos ejemplos, enseguida realizaremos el
ejemplo 2, que consiste en obtener los resultados de operacin de un tren de evaporadores
en de doble efecto en corriente directa.

Observemos como se complica el lgebra a medida de que se agregan etapas al tren de

evaporadores, y tambin que tipos de suposiciones se realizan para poder simular estos
sistemas.

Ejemplo 2.
Se desea concentrar 10,000 kgih de una solucin de NaOH del 12% al 40% en un tren de
evaporadores de doble efecto con alimentacin en corriente directa a 50 "C. Despus de
calentarlo en un intercambiador y disponiendo de vapor a 2.5 aim para el calentamiento en
el primer efecto, en el segundo efecto se mantiene un vaco a 545 mm Hg manomtricos y
el vapor procedente del evaporador se condensa a 25 "C.Si los coeficientes globales son
1900 y 1200 kcal/m2h"C para el primero y segundo efecto respectivamente calcule la
superficie de calentamiento para cada evaporador y la economa del tren de evaporadores.

Primero observemos grficamente como podremos resolver el ejemplo.

47

Copirulo P7.'l.

i\de/odologode siriiiilnci9n.

Tren de evaporadores de doble efecto.

con alimentacin en corriente directa.
Ei
l-----l

E2

+
V

Ualnnce [le soluto pnrtl

Efecto I
1) F = E i + C i
2) X r F = C i Xci
Efecto n
3) C i = E i + C z
4) x c i ci = c2 xc2

C A ~ efecto
A

Bnlrrnce de ninstl globnl

F =El + E2 + C2 = ETot
5) X F F = C ~ X C ~
6) F =Etot + C2

Bulnnces de Energa
8) V = Et hi - F Cp (Tf - Tebl) / hi
9) E2 = E2 h2 - Ct Cp (Tebl - Teb2) hi
Velociritld de trtlnsmisin de cnlor
10) qi = Ui Ai ATi
11) qi = Vi hcondi

12) AT = (Tebl - Teb2)

13) AT1 = AT tot I [ 1 + Ui/U2 ]

Sometido a las siguientes supocisiones y restricciones:

o Las entalpas de condensacin son iguales hi = hz
o Las condiciones de diseo requieren que Al = A2
o La vel de transmisin de calor en cada efecto debe ser igual qi = (12
Esto anterior es bsicamente y en esencia el criterio de convergencia.
Ahora resolvamoslo algebraicamente.

Soliicin.
Con los datos de alimentacin podemos hacer el balance de soluto
ec.5. C2 = FXr/ Xc2 = (10,000)(0.12) 0.40 = 3,000 k g h
supongamos que el evaporado del efecto I El = O, por l o que en la ec.6 nos da
F = E, + E2 + C2 ; Etot = F - C2 = 10,000 - 3,000 = 7,000 kg/h
Por lo tanto EI = E2 = 3,500 kgh
ec.3 Ci = F -El = 10,000 - 3,500 = 6,500 kg/h
calculando las fracciones en las ecs. 2 y 4
&i=XrF/
=(0.12)(10,000)/6,500=C~.184= 18.4%

Coprulo I TU. Afetocidogio

siitiidnoon.

Ahora bien, para calcular la temperatura de la entrada del vapor

P = 545 mm Hg ; entonces (760 - 545 )mm Hg - 1 atm / 760 mm Hg = 0.282 atm
Por lo que buscando en tablas de vapor corresponde a una temperatura de 128 "C
Buscando entonces las diferencias de temp:rs:tura ec.12 y 13.
Ar1= ATtot / [ 1 + UiAJz ] = (128 - 95) / 2.583 = 12.77 "C
AT2=ATtot -AT1 = 128 - 12.77= 115.23 "C
Referidos a las tempt'ratura referidas y el diagrama de entalpa cmcentracin
HE1 = hg + CpAT = 643.4 kcalkg HV = 649.85 kcal/kg
hfgl= 520.8 Hf=50 kcallkg
hcl = 95 kcal/kg
Por l o que en las ecs.8 y 9
y
El = 3,660 k g h
V = 4,549.54 kgfh

Finalmente calculando para la transmisin de calor y las reas en las ecs. 10 y 11.
(11 = V h= (4,549.54)(520.8) = 2.37 x lo6 kciih
q 2 = E1 h = (3,660)(520.8) = 1.93 x lo6 kcaV1i
por lo tanto
A I = (11 / UI ATl = 97.6 m2
Az= q2 / UzATz = 79.6 m2
Este mismo ejemplo resolvamoslo con el programa de simulacin, dando las condiciones
de entrada tenemos que en la Fig. VIIT.5.

Fig. VIII.5. Condiciones de entrada para el sistema de tren de evaporadores de doble efec$0
con alimentacin en comente directa.

49

C'npilu!o Wll. iletodologio de sinrulncibn.

..

Lo que nos dara como resultado los siguientes datos en la Fig. ViII.6.

Fig. VIII.6. Resultados arrojados por el programa al simular el ejemplo 2.

Este ejemplo se puede de i c .ual manera generalizar a un sistema de tren de evaporacin con
tres efectos. Resolvamos ihora un tercer ejemplo en que el sistema sea un poco ms
complejo como lo sera un sistema de tren de evaporacin de tres efectos con alimentacin
en paralelo y analicemos el sistema.

Ejemplo 3.
Se desea disear un sistema de evaporacin triple efecto para concentrar el soluto de una
solucin del 10 % a una solucin del 50 %'O en peso con un flujo de alimentacin de 50,000
k g h y entra al primer efecto como lquido a 100 "C debe usarse alimentacin en paralelo.
Para cumplir con los requisitos de ca1entar;iiento del primer efecto, se utiliza vapor saturado
a 250 "C y el tercer efecto debe ser operado a una presin correspondiente ai punto de
ebullicin del solvente puro de 125 "C. Considere despreciable la elevacin en los puntos
de ebullicin. Considere tambin el Cp para alimentacin y otras corriente lquidas de 1
kcaVkg"C. Los coeficientes globales son 500, 300 y 200 kcal/hrm2"C respectivamente para
cada etapa. Y suponga adems os calores de vaporizacin iguales a h = 1,000 k c a n g .

Veamos el mtodo de resolucin algebraica y sus modelos matemticos.

Copitulo Vlll. hfetotiologo de siniulocin.

.-

Tren de evaporadores de triple efecto.

con alimentaci8n en paralelo.

Bnlrnce de soluto pnrn codla efecto

Efecto I
1) F = E i + C i
2) XfF=CiXci
Efecto TI
3) C i = E 2 + C 2
4) x c i ci = c2 XCZ

Rnlnnces de Energa
9) V = E l hi - F Cp (Tf - Tebi) / hi
10) E2 = E2 1 2 - CI Cp (Tebi - Teb?) / hi
11) E3 = E3 h3 - CZCp (Tebz - ATtot) / hi
Veloci~lnrlde transmisin de calor
1;) qi=UiAiATi
13) qi = Vi hcondi
14) AT = (Tebi - Teb,)
15) AT1 = AT tot / [ 1 + UiAJ2 + Uim3 ]

Blalnnce de masa gfobnl

F = E , +E2 +E3 + C3 =ETot
7) XfF=C3Xc3
8) F = E t o t + C3

Sometido a las siguientes supocisiones y rest.ricciones:

o Las entalpas de condensacin son iguales hi = h2 = A3
Las condiciones de diseo requieren que Ai = AZ= A3
o La vel de transmisin de calor en cada efecto debe ser igual qi = q 2 = q 3
Esto anterior es bsicamente y en esencia el criterio de convergencia.

51

Ctipitulo W I . ibetociologa de siniulocin.

Solrlcitr.

Primero realicemos los balances de masa (de solvente y soluto mediante las siguientes
ecuaciones.
Concentrado en el efecto 3.
ec.8) C3 = xF * F ixC3 = 0.10 * 50,000 10.50 = 10,000 kglh
La cantidad de evaporacin total.
ec.7) Etot = F - C3 = 50,000 - 10,000 = 40,000 k g h
Realicemos la siguiente suposicin, cada etaisa evapora la misma cantidad de agua.
E! = E2 = E3 = 1/3 Etot = 40,000 / 3 = 13,333.33 kc/h
Ahora realicemos el balance por efecto individual
Balance en el efecto I
Ci = F -El = 50,000 - 13,333.33 = 36,667 kgh
ec.1) F = EI + CI;
xCi = xF * F /Ci = 0.136
ec.2) F * xF =xC1 * Ci;
Balance en el efecto I1
C2 = Ci -E2 = 36,667 -J3,333.33
ec.3) C I = E2 + C2;
ec.4) xCi * Ci = XCZ* C2; xC2 = xCi * Ci / CZ= 0.2137
Balance en el efecto III
ec.5) CZ= E3 + C3;
ec.6) xC2 * C2 = xC) * C,;

= 36,667

kfl

C3 = 10,000 k g h
XCS= 0.50 coridicin de diseo

Ahora encontremos las diferencias de temperaturas para aplicar posteriormente un balance

de energa.
ec.14) ATtot = (TvapI - Tebr~r)= (250 - 125) = 125 "C
ec.15) AT1 = ATtot / ( 1 + 1J1/U;?+ Ui/U3 ) = 24.193 "C
AT2 = ATi ( Ui/LJ2 ) = 40.322 "C
AT3 = AT! ( Ur/U3 ) = 60.483 "C
AT1 = Tv - Teb~;

Tebr = Tv - AT1 = '225.81 "C

TebiI = Tebr - AT2 = 185.39 "C
TebiiI = 125 "C condicin de diseo

Procedamos ahora a realizar el balance de energa por efecto

Balance de energa efecto I
ec.9) V = E l l - FCp (Tref - Tebr) / h = 19,623.5 k g h

57

Coorulo JfIf.

I\

eiodoloeio de anitrlocin.

Balance de energa efecto Ti

ec.10) E2 = E2h - ClCp (Tebi - Tebll)/ h = 12,255.12 kg/h
Balance de energa efecto III
ec.11) E3= E3h - C2Cp (Tebii - ATtot) / h = 1 1,923.92 k g h
Ahora con estos datos obtenidos podemos calcular la velocidad de transmisin de calor y a
su vez la superficie de calentamiento. Ecs. 1 3 y 12 respectivamente.
q l = V h = 19.64 x lo6 kcal/h
q z = E i h = 12.55 x IO6 kcallh
q 3 = E& = 1 1.92 x 1O6 kcallh
Para el rea de calentamiento
AI = qi / UI AT1 = 1,622.44 mz
-42 = q 2 / UzAT2 = 1,013.14 mz
A3 = q3 / U3 AT3 = 986.91 mz
Este clculo de igual manera lo haremos con el programa, lo primero que tenemos que
hacer es capturar los datos de entrada como lo muestra la siguiente figura VIII.7.

Fig. ViII.7. Condiciones a simular de un tren de evaporadores de triple efecto con

alimentacin en paralelo

53

Un2 vez capturado los datos, mostremos los resultados de la programacin mediante el
simulador y comparemos los resultados de aiiibos mtodos. Figura VIII.8.

Fig.. VIII.8. Resultados de la simulacin de un tren de triple efecto con alimentacin en

paralelo.
Haciendo una comparacin del mtodo algebraico al mtodo numrico encontramos que no
hay ninguna diferencia y mucho mejor an, que la rapidez de solucin con el simulador es
ms corto que el algebraico.
Por otro lado podemos hacer pequeas variaciones sobre las condiciones de entrada y
observar como repercuten stas en los resuitados. As mismo podemos caracterizar un
diseo de tren de evaporadores y tener otra:; alternativas mediante la simulacin.
Para demostrar que si se puede interaccionar con el simulador la siguiente figura VIII.7, nos
muestra como se pueden variar las condiciones de entrada, esto puede hacerse para cada
tren individualmente, aqu para una mejor y rpida demostracin hemos elegido a un
evaporador de simple efecto.
Supongamos que en la figura VIII.3 del ejemplo 1, no equivocamos al dar una condicin de
entrada o simplemente que queremos ver como se comporta el sistema con una variacin.

54

Cnpitulo Vili.

Figi
efec

I\ etodologin de

siiiiulocirt.

nple

La figura VII.9. Muestra como hubo un error o se desea corregir las condiciones de
entrada para el inciso e), que es la presin bajo la cual opera el evaporador que
anteriormente era de 1.93 atm y el nuevo valor ser 2.53.
Con el simulador se puede realizar esto de manera inmediata y observar los cambios y

variaciones que presenta el sistema por mnimos que estos sean, mientras que en el mtodo
algebraic0 la nica solucin para esto es repetir nuevamente todo el proceso.

55

Como conclusin puedo decir que el niarco terico presentado es niuy bsico y conciso
dirigido a estudiantes de ingeniera que comienzan con el estudio de operaciones unitarias,
esto logra mayor integracin e interpretacin de los sistemas de evaporadores.
Adicionalniente en el desarrollo de este trabajo fue de gran ayuda u n algoritmo y cdigo de
simulaciri de solucin a las eciiaciones del sistema de evaporadores. Estas ecuaciones por
mtodos convencionales matemticos o algebraicos no hubieran podido resolverse tan
fcilmente, adicionalmente si se llevaran a cabo tratamientos matemticos dificultara
demasiado su comprensin.
El objetivo especfico de este trabajo es el de lograr un cdigo de simulacin, este
algoritmo de solacin (anexo A), presenta caractersticas que por su facilidad de uso y
compresin mediante el cdigo de simulacin unifican al usuario, marco terico y
entendimiento de los sistemas de evaporadores dando cabida a que este objetivo se
satisficiera de forma por ms que especifica.

Otra conclusin que encuentro apropiada, es la facilidad con que se pueden resolver los
sistetnas de evaporacin con el programa de simulacin, adems de que si se resuelven
alternativamente por ambos mtodos (algebraico y programa) el tiempo de resolucin es
mucho menor con el simulador que por el mtodo algebraico.
Por otro lado veo que los resultados de simulacin son idnticos a los obtenidos con el

netodo algebraico que se realiza con gran exactitud matemtica. Otra ventaja que es
evidente es que con el programa se puede corregir o modificar las condiciones de entrada
observando no solo un sistema especfico, sino que se pueden estudiar sistemas alternativos
y hacer comparaciones para el diseo adecuado de un tren de evaporacin.

No dudo que pudieran existir ventajas intrnsecas en el uso de la siniulacin para sistemas
de dificil acceso como lo son los trenes de evaporadores, pero confo en que cuando se
difunda con mayor seriedad y empeo lo importante de sus resultados as como su utilidad
stas sern ms visibles.
S i n embargo considero que la evolucin natural de la industria y el avance de la tecnologa
cornputacional que se presenta en nuestros das lograr que se susciten mayores y mejores
programas de simulacin que dar como resultado una excelente preparacin a aquellos
nuestros futuros ingenieros de una manera intrnseca.

Max S . Peters & Klaus D. Timmerhaus, Plant design and economics for chemical
engineering, McGraw Hill, 1980 Tokio, J.spon.

T.F. Edgar & D. M. Himmelblau, Optimization of chemical processes, McGraw Hill,

1988 USA.

Coulson & Richardsons, Chemical engineering, Vol. 6 (Chemical engineering

design), Pergamon Press, 1994 Great Britain.
o

Angel Vian y Joaqun Ocon, Elementos cle ingeniera qumica (Operaciones bsicas),
Aguilar, 1979 Madrid.

Joaqun Ocon y Gabriel Tojo, Problemas de ingeniera qumica (Operaciones bsicas)

Tomo I, Aguilar, 1974 Madrid.

Gael D. Ulrich, Diseo y economa de los procesos de ingeniera qumica,

Interamericana, 1986 Mxico.

Warren L. McCabe y Julian C. Smith, Operaciones bsicas de ingeniera qumica,

McGraw Hill, 1991 Espaa.

Alan S. Foust y Leonard A. Wenzel, Principios de operaciones unitarias, Continental,

1996 Mxico.

Charles D. Holland, Fundamentos y modelos de procesos de separacin, Prentice

Hail, 1981 Madrid, Espaa.

Alejandro Torres Aldaco, Notas de la UEA, Transferencia de Calor y Masa, 97-0

AKneAxo
ACdigo de Simulacin de Tren de evaporadores.

.-tnexo i l . C(jdi~p
de siinir!ocin de tren de evoporodores

/'Proqrsma

Con la
#include
#include
#include
#include
#include

re?li.-?do

asesocia d e
<stdin.h>
<math.h>
<conin.h>
<ctype.h>
<process.h >

p o r rnrj. Enriqiie ?.wiora !?ario.

M.

en

C.

Alejandrn Torre!; Aldaco.*/

char vaior,casos, fin, resp;

int E;
flnst ms:
doiible F.::F.:~:C,pv,pe,Ta,Te,Td,T~,Tp,Hf,Hc,H~~,1,U,rl,r?,C,a,b,c,d,e,q,ee,A;
double T$.-,'Pv,h C 1 , t i t ~ l , i J ~ , : ~ F F , ~ : [ : ~ ' ,S~I~- '. ,~E, l , C l ,xC1,xCCl,dTl,dTZ. hf ,ql,qZ,hEl,Al,A?;
d<?iit>leE,ET,C3,XC?, XC3,4Tt,dTj,Tebl, T4,2,'ret>3,U3,Tvl,Cp,q3,A3,V,
E2, E3;
void datos (veid)
do( 8 = 1;
cl rscr I ) ;
printf i"\n\t\t\t**' Tabla de datos **'\n"J;
print f I " ',na I F l u j o de alimentaciCn.(kq/h)
print f ( " \ nb ) Poncentaje en peso de la disoluciCn.
D r i n t f I " \ nc i Porcentale en peso del concentrado.
print f ( " \ndI PresiCn del vapor de aqua.(atmi
printf ("\ne) PresiCn del evaporador.(atm)
printf I"\nf) Temperatura de alimentaciCn.(SCl
princf ("\nrJ)Temperatura de ebulliciCn del aqua. (5C)
print f ( " \ rih ) Temperatura de ebulliciCri de le disoluciCn.(C)
printf ("\ni) Temperatura de condensdcicn del. aqua.(SC)
printf ("\nj) Entalpia de alimentaciCn.(kCal,'kq)
princf ("\nti Entalpia de la disoluciCn concentrada.(kCal/kg)
printf ("\nl) Entalpia del vapor de agua. (kCal/kq)
print f ( " \ nm i Calor de vaporizaciCn del vapor de agua.IkCal/kg)
print f ( " \ nn ) Coeficiente gl@b.~:.(W/m2hSC)
printf ("\n\nLos datos de e:itrada son correctos.. . ( S / N I " 1 ;
resp = getche0;
if iresp == 'n' I 1 resp == I N ' )

%f",Fl;
%f",xF);
% f",XC);

%f",pv);
% f",pe) ;
%f",Tal ;
% f'' ,Te ;
%f",TdJ;
% f" ,Tc 1 ;
%f",Hf ) ;
% f" ,HC ;
% f",Hv);
% f",1 I ;
6 f",u I ;

i
printf("\nCu 1 dato no es
valor=getcheO;
switchlvalor)
I
case ' 3 l :
printf("\n\nCu 1 es
scanf("9lf",hF);
break;
case 'b' :
print("\n\nCu 1 es
scanf ( " % 1 f" , S x F ) ;
break;
case IC':
printf(',\n\nCu 1 es
scan f i " 8 1 f" ,6xC1 ;
break;
case 'd' :
print("\n\nCu 1 es
scan f ( " 8 1 f spv 1 ;
break;
case # e ' :
printf("\n\nCu 1 es
scanf ( " % 1 f",spe J ;
break;
case I f ' :
printf("\n\nCu 1 es
scanf ("%lf",hTal;
break;
case '-2'
:
printf("\n\nCu 1 es
scan f I ** 8 1 f" ,6Te) ;
break;

correcto...

"1;

el flujo d e ilimentaciCn. (kg/hl " I ;

el poncentaje en peso de la disoluciCn.(%) " I ;

el porcentaj,? en peso del concentrado.(%) " I ;

la pcesiCn dtel vapor de aqua.(atml

"1;

'I,

la presien d,?l evaporador. (atml " I ;

la temperatura de alimentaciCn.(C) " I :

la temperatura de ebullicicn del aqua.(C) " I ;

I: ite.ro .:!

C35'

Ctkli,po de sirriidncin de tren de evnporotiores

'3':

printf["\n\nCu 1 es la remperatura de ebullicicn de la disoluciCn.(C) " I ;

scan f ( " * 1 f" , 6Td I ;
break;

case 'i':
printf("\n\nCu 1 es la temperatura de condensacien del aqua.lC) ' * I ;
scanf l"SlC",&TcI;
bredt ;

caso 'j':
printf("\n\nCu 1 es 1.3 erit3lp;a de alimentacien. (kCal/kgI " I ;
scan f ( " $ 1 f", &H f I ;
hceat;
case I t * :
prinrf ("?n\nCii 1 es la entdlpjo de l a disnluciCn concentrada. (kCal/kql " ) ;
scqnf ( " 1 1 f",&lcI;
break;
case ' 1 ' :
printfl"\n\nCu 1 es la ent3lp;a del vapor de aqua.(kCal/kq) " I ;
scanf [ "9 1 f" ,6Hv I ;
break;
case ' i n ' :
pcintfl"\n\nCu 1 es el calor d e vaporizaciCn del vapor de agua.(kCal/kql "1;
scan f ( " % 1 f" ,6 1 1 ;
break;
case In' :
printf("\n\nCu 1 es el roeficienie global. (W/mZhCI " I ;
scanf ( "9 1 f" , s l l I ;
break;
default :
ClKSCr( 1;
printf ( " \ n \ n \ n \ n \ n \ n \ n \ n \ n \ n \ n \ : \ t T l e n e s
que seleccionar alguna opcicn.. . " ) ;
qetche ( I ;
break;
I

1
else I
if (resp == ' s ' Ilresp == 'C'I@=O;
else (
pcintf("\n\nColo responde si los idatos son correctos
getchell;)
1
1 whilo(B1;
rl=100 .O-xF;
rZ=10O.O-xC;
a=r1 / :< F;
b=r 3/ XC ;
c=a-b;
d=c* F*xF/l OO. O;
e=$-d;
Tp=Td-Te;
q= (F-el-(Hvl+e*Hc-P'Hf;
ms=q/ 1 ;
ee= (d/ms1'100. O;
A=q/ (U'(T.C-T~I I ;
ClrSCK( 1;
printf("\t\t\t\tTabia de resultados...");
printf ("\n\ri\nEl consumo de vapor de aqua e.s.(k3/h):
printfl"\nLa cantidad de aqua evaporada es.lkg/hl:
prinrl"\nLa velocidad de flujo de diso1uci:n es.lkg/h):
printf("\nLa velocidad do transmisiCn de calor es.(kCal/hl:
printf("\nLa elevacicn da1 punto de etullicicn es.(SCI:
printt("\nLa economja del evaporadoc es.:
princf("\nLa superficie de calentamiento es.lm21:
printf("\n\n\n\n\n\n\n\n\n\n\n\n\nPresiona r para volver...");
pcintf("\no presiona q para terminar..."];
fin = getche0;
switch (fin)

o no.

"I;

% .Zf",msI i
%.Zf",d);
% . 2 f " ,e I ;
%.2f",ql;
%.Zf",Tpl;
%.2f",ee/100.0):
%.Z",R) ;

case ' q ' : break;

case 'r': d a t o s o ;

Arteso .4. C6digo de sirtiulncin de tren de evoporndores

____

default:

break;
cl rscr ( 1 ;
printl"\n\n\n\n\n\n\n\n\n\n\n\tTornaste

una mala opcitn, tal vez quieras

modificar la tabla.");
qetche I I ;
ddtos ( 1 ;
break;
1
)

void 3datosds (void)

do( 8 = 1;
cl rscr(1;
printfl"\n\t\t\t'** Tabla de datos *"\n");
printf("\nal Filijo de alimentacien. (tq/h)
printf("\nb) Poncentaje en peso de la disoluciCn.
printfI"\nc) Porcentaje en peso del rcmcentrado.
printfi"\nd) PresiCn del vapor de aqua. (atm)
printf("\ne) PresiCn del evaporador. (atm)
printfi"\nf) Temperatura de alimentacicn. ( J C )
printfi"\nq) Temperatura de ebulliciCn a la salida.(SC)
printfi"\nh) Temperatura de VaporizaciCn de agua.(%)
printf("\nil Temperatura de condensaciCn del aqua.15CI
printf("\njI Entalpia de alimentaciCn.lkCal/kql
printf("\nt) Entalpia de vaporizacicn del efecto Z.ikCal/kg)
printf("\nl) Entalpia del vapor de aqua. (kCal/kgl
printf("\nml Entalpia de la disoluciCn concmtrada.lkCal/kgl
printf("\nn) Calor de vaporizaciCn del vapor de aqua. (kCal/kq)
pcintf("\nol Coeficiente global del efecto i.(W/mZhSC)
printf("\npl Coeficiente global del efecto ;!.iW/mZhC)
printf("\n\nLos aatos de entrada son carrectos... lS/N) " I ;
resp = getche0;
if iresp == In' I I resp == 'N'I
L

print("\nCu 1 dato no es
valor=getche();
switch(valor1
I
case 'a':
printfl"\n\nCu 1 es
scanf 1 " % 1 f" ,6 L);
break;
case 'b' :
pcintf("\n\nCu 1 es
scan f 1 " % 1 C " , 6 x F 1 ;
break:
case 'c':
printfl"\n\nCu 1 es
scan f ( " % 1 f '' ,6 xC 1 ;
break;
case Idt:
printf("\n\nCu 1 es
scani"%lf",spvl;
break;
case 'e' :
pcintf ("\n\nCu 1 es
scanfl"%lf", 6pel;
break;
case 'E':
printfi"\n\nCu 1 es
scanf("%lf",cTa);
break;
case '3':
pcintfl"\n\nCu 1 es
scanf i "%1E", sTe2);
break;
case 'h' :
printE("\n\nCu 1 es
scan f i " % 1 f " ,6Tv ) ;

correcto..

%f", Fl;
% f",x F ) ;
% f " ,xc 1 ;
% " , PVI ;

%E", pe i

%",Tal ;
% ",Te2I ;
9 ",Tv);
% f",Tc) ;
% f " ,Hf 1 ;
%f",hCll;
% f" ,Hv ;
%f",HEll;
%f",l);
% f" ,u 1 ;
%",U21 ;

"I;

el flujo d e alimentacitn. lkq/hl " I ;

el poncentajo en peso de la disoluciCn.(%) " I ;

el porcentaje en peso del concentrado.(%) " I ;

la presien dol vapor de aqua.latm1 " I ;

la presiCn del evaporador.(atm) " I ;

la temperatura de alimentaciCn.lC: " I ;

la temperatura de ehulliciCn a la salida.iCI " I ;

la temperatura de evaporizaciCn del aqua.(%)

"I;

break;
case '1':
printEl"\n\nCu 1 es la temperatura de condensacien del aqua. (SC] " ) ;
scanf I " % l E", 6Tcl ;
break;
case ' J ' :
priniEl"\n\nCu 1 es la entalpia de alimentaciCn.lkCal/kg) ' O ;
scan E 1 " % 1 f " , S Hf 1 ;
break;
case ' i : ' :
prirttfl"\n\nCu 1 es la ent31pia de vapnri:aci+n del efecto 2. lkCal/kq) " ) ;
scan f 1 '"t 1 E", S hi31 1 ;
break;
case '1':
printEl"\n\nCu 1 es la -ntilp;a (del vapor de agua.lkCal/kg) "1;
scanf l"$lf",SHY);
break;
case * m ':
printfl"\n\nCu 1 es la entalpia ie disolucien concentrada.(kCal/kq) " I ;
scanf l " % l f " , SHE1) ;
break;
case In' :
printfl"\n\nCu 1 es el calor de vaporizaciCn del vapor de aqua.lkCal/kql " ) ;
scanf I " % 1 E", 6 1 I ;
break;
case 'o':
printEl"\n\nCu 1 es el coeficiente global del efecto l.(W/mZhSC) " I ;
scan E I" a 1 t " ,SU ) ;
break;
case 'p':
printE("\n\nCu 1 es el coeficierte global del efecto Z.(W/rnZhC) " I ;
scanfl"6lf",6IJ21 ;
break;
default :
ClKSCK[ ;
printf ( " \ n \ n \ n \ n \ n \ n \ n \ n \ n \ n \ n \ t \ t T i e n e s
que seleccionar alguna opciCn
1;
getchel I ;
break;

..."

1
else (
if (resp == ' 5 ' Il'resp == 'S'lB=O;
else (
printf("\n\nSclo respnde si los datos son ~ ~ c r e c t oos no.
qetchel ) ; I
1
} while(B1;
xPF=xF/100.0;
xcc=xc/100. o;
C2=Pcx FF/xCC;
El 2=F-C2;
El-E12/2.0;
Cl-F-E1 ;
xCl=xFF'F/Cl;
xcc1=xc1*100.0;
dTl=llT~-Te~l/I1.0+lU/~J~l~~;
dT2=ITv-Te?l-dTl;
V-lE1*HE1+Cl*hC1-PChf]/l;
ql=V'l;
q2=E1'1HEl-hE11;
Al=ql/ (U'dT1 I ;
A2=qZ/(U?'dTZI;
clrscr I I ;
printfl"\t\t\t\tTabla de resultados...");
printf("\n\n\nEl f l u j o de la disoluciCn del efecto 2 es.lkg/h):
printf("\nEl flujo de la disoluciCn del efecto 1 es.lkq/hl:
printfl"\nEl concentrado del efecto 1 es.lkq/hl:
printfl"\nLa fracciCn que concentra el efezto 1 es.
printfl"\nLa delta de temperatura en el efecto 1 es.1921:
prrntfl"\nLa delta de temperatura en el efecto 2 es.ICl:

"1;

%.2:",C2);
%.2f",Cl);

% . 2 f " ,El I ;

XCCl 1 ;
%.Zf",dTlI;
%.Z",dTZ);
% Zf",

___-

..liie.ro . I. Cci%qo
de .~iniuincin
de tren de evopormiores

?.2f",qll;
printfI'*\nl..svel. de t r a n s m i s i C n de r i l n r e f e c t o 1 es.(i:C3!/h):
printf("\riI.i vel. do transmi::iCn ci? c:il,Jr efecto 2 es. (kCai/h): %.2f",q2);
printf l"\nl,a canrid3d tot.11 de evipora,?itn es. lkCal/h):
% .2E", v I ;
printfl"\nl.s superficie de calentani~erirodel efecto 1 es. (m2): %.2f",Al);
pr'intE("\nLa superficie de raienr3iiiienr~del efecto 2 es. lm21: %.2f",A2);
print f ( " \ n \ n \ n \ n \ n \ n \ n \ n \ n \ n P r e s i o n a
r para volver.. . " I ;
printfl"\no presiona q para terminar.. . " I ;
fin = qetchell;
switch ( f i n 1
I
case ',-,': hCe3k;
case ' K' : i?it^Srl.; I ) ;
break;
d e f a u I t : clrscr ( I ;
printf ( " \ n \ n \ n \ n \ n \ n \ n \ n \ n \ n \ n \ f T o n a s t euna mala opcicn. tal vez quieras
modi fica r 1,3 t ab13 " ) ;
qerche ( 1 ;
darosds I I ;
break;
1

i
void datosdp (void)
(

do( 5
1;
clrscr( ;
pcintf("\n\t\t\t*** Tabla de datos +*+\n"l;
ori ri c f ( " \ na I Flitjo ,de alimentxitn. Ik,y/h)
p r i n t E I '' \ nb 1 Poncentaje en peso de la disolucien.
prlntf("\nc) Porcentaje en peso del concentrado.
printf ("\n4) Presien del vapor de aqua.latm1
p r int f ( "\ne 1 PresiCn del evaporador. l a t n i )
print f I " \ n E I Tempe:atura de a 1 inientaci C n . 15C)
p r i n t f ( " \ nq I Temperatura de ebullicitn a la salida.(SC)
printf l"\nhl Temperatura de vaporisacicn de agua.15C)
printf ("\niI Temperatura de condensacitn del agua.(5C)
printf("\nj I Entalpia de alimencacien.(kCal/kql
pr i n t f I \ nk I Entalpia de vaporizacien del efecto 2.1kCal/kg)
print f 1 " \ nl 1 Entalpja del vapor de aqua.IkCai/kq)
printE("\nml Entalpia de la disoluciCn concentrada.ikCal/kql
p r in t E ( *' \ nn ) Clor de vaporizacicn del vapor de agua.(kCal/kg)
printf ("\no) Coeficiente global del efecto l.(W/m?hSC)
p r i n t E ( I' \ np ) Coeficiente t3loba1 del efecto I.(W/m?hSC)
printf("\n\nLos datos de entrada son ccrrectos... ( S / P I )
resp = qetche0;
if (resp == ' n ' I 1 resp == 'NI)
printE("\nCu 1 dato no es correcto...
vaIor=qetche();
switch (valorI

%f",F);

%f",xF);
% f XCI ;
% f" ,pv1;
% f",pe);
% E " ,Ta 1 ;
% f " ,Te 2 ) ;
%f",Tvl;
% f",Tc) ;
Bf", Hf 1 ;
%E", hC1 I ;
% E " , Hv) ;
", HE1 I ;
%f",l1;
% f" ,U ) ;
'1,

%f",U2);

");

'I);

case lai:
printf("\n\nCu 1 es
scanf ("%lf",sFI;
break;
case 'b' :
printf l"\n\nCu 1 es
scanf("%lE",&xF);
break;
case
print("\n\nCu 1 es
scanf ("Ylf", SxCI;
break;
case Id' :
printEl'*\n\nCu A es
scan I " % 1 fI * , Spv I ;
break;
case ' e ' :
printfl"\n\nCu 1 es
scan f I " % 1f" ,spe I ;

el f l u j o de alimentacien. lkg/hl " I ;

el poncentaje en peso de la disoluciCn.i%) " I ;

el porcentaje en peso del concentrado.(%)

");

la presicn del vapor de agua.(atm) " I ;

la presiCn del evaporador. (atml

");

Anexo . I. C M g o de siitiiriocidn de Iren

de evoporndores

break;
case ' f :
printf("\n\nCu 1 es la temperatura de alimentaciCn.(5Cl " I ;
scan f ! " % 1 f " , 6Ta I ;
break;
case '1':
print("\n\nCu 1 es la temperatura de ebulliciCn a la salida.(C) " I ;
scan E i "R 1 E", 6Te2I ;
break;

case 'h':
printf("\n\nCu 1 es
scan f 1 '<% 1 E", 6Tv I ;
break;
case *i':
printE("\n\nCu 1 es
scanf i " B 1 f" , G Tc I ;
break;
case '3':
printf("\n\nCu 1 es
scanf ( " % 1 f" , 6 H f 1 ;
break;
case 'k':
printE("\n\nCu 1 es
scanf ( "9 If", 6hCL 1 ;
break;
case ' 1 :
prrntf("\n\nCu 1 es
scan("lf", 6Hvl ;
break;
case 'ml:
printf("\n\nCu 1 es
scanf ( " 3 1f " , 6 HE 1 1 ;
break;
case 'n' :
printf("\n\nCu 1 e s
scan( " S 1f", 6 11 ;
break;
case 'o':
pcintEl"\n\nCu 1 es
scanf i " % 1 f ",SUI ;
break:
case ' p ' :
pcintf i"\n\nCu 1 es
scan E ( " % 1 E" , 6 U2 I ;
break;
default :

la tpmperatura de evaporizaciCn ael agua.(%)

13

"I;

temperatura de condensacien del agua.(SC) " I ;

la entalpia de alimentaciCn.[kCal/kql " I ;

la entalpia de vaporizacien del efecto Z.(kCal/kgl " I ;

la entalpia del vapor de aqua.LkCal/kg) " I ;

la entalpia de disoluciCn concentrada.(kCal/kgl " I :

el calor de vaporizacicn del vapor de agua.(kCal/kgl " I ;

el coeficiente global del efecto l.(W/m2hCI " I ;

I

el coeficiente global del efecto 2 . iW/rn2hC) " I ;

clrscr ( ) ;

print ("\n\n\n\n\n\n\n\n\n\n\n\t\tTienes que seleccionar alguna opcicn

getche( I ;
break:

I
1
else I
if iresp =- ' 5 ' Ilresp == 'S'IB-O;
else [
printE("\n\nSolo responde si los datos son correctos o no.
getcheil;)

..."I :

"I;

while(B1;
xFF-xF/ 100.O;
xcc=xc/loo. o;
C2= F'xFF/ xCC;
El 2= E-C2;
El=E12/2.0;
Cl=F-El;
xCl=xFF* F/Cl;
xcc1=xc1+100.0;
dT1=~Tv-Te21/11.0+1U/U2));
dTZ'=(Tv-Tell-dTl;
V=iE1*HE1+ClChCl-FChEl /1;

.'!ne.roit. Cdigo de siriririncin

de tren de evoporodores

ql=V*l;
q2=El*IHEI-hEll;
A l = (V'l I / (U'dT1 I ;
A2=q2/(Ul'dT21;
c1rsrro;
printf("\t\t\t\tTabla de resultadm.. . " I ;
%.2f",Cl)i
printL("\n\n\nEl flujo de 1 3 aisnlucitn del efecto 2 es. (kg/hl:
printf("\nEl flujo de 13 disnliiciCri del efecto 1 es. (kg/h):
% . Y " , c2 1 ;
printfl"\,nEl cnnc?ntcado del e f e r t o 2 e s . ikq/hI:
% .2f",El I ;
printfl"\nLa EracriCn que concentra el efecto 2 es.
% .2f xcc1 I ;
% .2f",dT2 I ;
prir#tfi"\ni.adelta de tempecatiira en el efecto 1 es. ISCI:
princf ("\nLa delta de teniperatiira on el efecto 2 es. (4Cl:
% . 2 f ,dT1 I ;
priritC("\nLa vel. de :ransinisiCri de c a l o r efecto 1 es.(kCal/hl: % .2f" ,q2 I ;
printf("\nI.a vel. de transniisicn de calor afecto 2 es. (kCal/hl: % . 2 " , q l I ;
printf("\nLa cantidad total de evaporaciCn es.(kCzl/hl:
% . 2",V1;
% . 2 f " ,A2 I ;
printf l"\nLa superficie de calentaniiento dr?l efecto 1 es. (m!):
p r i n t f i"\nI.a superficie de calentaniiento d:l efec:o
2 es. ( m 2 1 :
% .2f",Al I ;
p r i n t ( " \ n \ n \ n \ i i \ n \ n \ n ? n \ n \ n \ n e r e s i o n a r para volver.. . " I ;
printf("\no presiona q para terminar.. . " I ;
fin = qetcheil;
swi:ch
ifinl
Z',

case ' , : J ' : break;

case Ir': datosdpll;
break;
default: clrscril;
peintf("\n\n\n\n\n\n\n\n\n\n\n\tTomaste
una

mala opcien, tal vez quieras

modificar la tabla.");
getche i I ;
datosdpl);
break;
1

void datostp (void1

do( B = 1;

ClKSCrI 1 ;
c>rintfl"\n\t\t\t+** Tabla de datos "*\n"I;
print f I * * na
\ 1 Flujo de alimentaciCn.(kq/hl
% f " . Fl ;
print E ("\nbl Poncentaje en peso de la disoluciCn.
%f",xFI ;
printf l"\ncl Porcentaje en peso del concentrado.
% f" ,xc31 ;
pci nt f ( "\nd I Temperatura de alimencaciCn.(Cl
%f",Tal ;
pr i nt f i" \ne I Temperatura de vapor de aqua. (5C)
% E", Tv I ;
% f" ,Teb3);
print f I " \ n f I Temperatura de ebullicitn a l a salida.lSC1
Calor
especifico
de
la
sustarcia.(kcal/kqCl
printf l"\nql
%f",Cpl;
printf ("\nh I Calor de vaporizacien del vapor de aqua. lkCal/kql % f", 11 ;
printf ["\ni) Coeficiente global del efectc, l.iW/mZhSCI
% f'*,u 1 ;
printf ("\njI Coeficiente qlobal del efecto 2.(W/mZhCI
%f",U2l i
print E ( " \nk I Coeficiente qlobal del efecto 3.(W/m2hSCI
% ",U31;
printfl"\n\n\n\n\n\nLos datos d e entrada son correctos... (S/N1 " I ;
resp = qetche ( 1 ;
if lresp == 'n' I I resp -= ' N ' I

printf("\nCu 1 dato no es correcto...

valor-qetche();
switch (valor)

"I;

case 'al:
printf("\n\nCu 1 es el flujo de alimentaciCn.(kq/hl " I :
scanf l"%lf", 6 F l ;
break;
case 'b' :
printfl"\n\nCu 1 es el poncentaje en peso de la disoluciCn.i%I "I;
scanf ("%lf", sxF) ;
break;
case IC' :
printf("\n\nCu 1 es el porcentaje en peso del concentrado.(%) " I ;
scanf ("%If", 6xC3);
break;

W':
prlrttf("\n\rLu 1 es
scarif 1 "P 1 E " , CT.3I ;
break:
case I r : , :
pcinrf("\n\nCu 1 e s
s c a n F I " 9 1 E " , 6Tv 1 ;
b re 3 :i ;
case 'f,:
prrntfl"\n\nrii 1 es
sc3 n f ( " 9 1 f " , 6 T rb3 1 ;
b re ~i
L: :
cast? '87' :
priritf("\n\nCii 1 es
scan f I " i1 E " , 6Cp t ;
braat;

case

cas-

la

rem?-raturs de alimentacicn. (SCI

"ti

la temperatura de vapor de agua a la entrada.(C) " I ;

la temeperatiira de ebiilliciCn a la salida.(5Cl " I ;

-1

calor espt?~;
fico de la solucicn. (kcal/kqSCI " I ;

'h':

printf("\n\nC~1 es P I c i l o r de 'iaporizacitn del v a p o r de agua.(kCal/kcg) ' * I ;

scan f I "Y 1 f " , 6 1 ) ;
break;
case ' i t :
Frintf("\n\nCu 1 es el coeficiente global del efecto l.IW/m2hSCl " 1 ;
scan[ " P l f " , 6 1 1 ) ;
break;

case ' j l :
printi("\n\nCu 1 es el coeficiente qlobal d e l efecto 2. iW/rnZhCl " I ;
s c3 n E ( " 9 1 E " , s U2 t ;
break:
case 'Ir' :
pcincff"\n\nCu 1 es el cneficiente global del efecto 3.(W/mEhSCI " I ;
scan E ( ' ' 9 ; 1f " ,6 U3 t ;
break:
default :
Cl KSCC ( ) ;

pcntf ("\n\n\n\n\n\n\n\n\n\n\n\t\tTienes que seleccionar alguna opciCn

getche ( 1 ;
bresk;

..."I ;

t
1

else [
i f (resp == ' s ' llresp == 'S')=C;
else I
printE("\n\nSolo responde si los datos son correctos o no.
getchp0;l

"I;

t whilelet:
C3-.xF+F/.?C3;
EC=F-C3;
E=Et/3;
Cl=F-E;
xC1=xP*F/Cl;

CZ=Cl-E:
xC2=xCl+Cl/C?;
dTt-Tv-Teb3:
dTl=irfTtI/il.O+(U/lJZI+lU/U3tI ;
aT2=U*dTl/U?:
Tebl=Tv-dT1;
TebZ=TebI-dT!;
dT3-TebZ-Teb3;
V=l(E'l)-l~'Cp*iTa-Tebtltl/1;
E2=(~E*1t-(C1*CptlTebl-Teb211)/1;
E3=((E*lt-(C2*Cp*lTeb2-dTttIt/l;

q1-v* 1 ;
q2=E?+l;
q3=E3*1;
Al=ql/ (U'dTlI ;
A2=qC/(U2*dTZl;
A3=q3/(U3'1T31;
ClCSCCI I :
printf("\t\t\t\tTabla de resultados...");

printC("\n\nEl flujo d e la disolucitn del ef

.lkg/hl:
printf("\nEi fii.ijo de la disnluciCn d e l efec
tq/hl:
printf ("\nEl flujo de l a disoluciCn del efec
kq/hl:
printf ("\nEl concentrado ?,?tal del sistema e
pcintf("\nEl concentrado de1 efecto 1 os.lkq/hl:
printfl"\nEl concentrado del efecto 7 es.I%q/h):
printf("\nEl concentrado del efectc 3 es.(kq/h):
printt("\nI.a fraccien que concentra el e f e c m 1 es.
printf('*\nLa fraccitn quo concentra el efect.3 2 es.
printf("\ni,a fracciCn que concentra el e f e c t ~3 es.
printf("\nLa dn1t.a de temperatura en el efecto 1 es.(SCl:
printf("?nL.n rir'lt3 de temperariir3 en el efecto 2 es.(4Cl:
printf ("\nt.a deltp de ternperatiira en e l efecto 3 es. (SCI:
priritfl"\nLa vel. de transmisiCn de ?alar efecto 1 es. (kCal/hl:
printfl"\nLa vel. d e transmisiCn de m l n r efecto 2 es.lkCal/hl:
printf("\nLa vel. d e transmisitn d e c a l o r efecto 3 es.(kCal/hl:
printf l"?nLa superficie de calentaniiento del efecto 1 es. (m21:
priritf("\nLa superficie de calentamiento del efecto 2 es.(mZI:
print f("\nLa superfimzie de calentamiento rfel efecto 3 es. (m21:
printf("\n\nPresiona r para volver.. . " I ;
printfl"\no presiona q para terminar.. . " I ;
fin = getche(1;
switch (fin)

%.2f",Cll;
%.2E",C21;
%.2f",C31;
% .2f",Et I ;
% .2f " ,VI ;
% .Zf", E2 I ;
% .2f",E31 ;
% .2 f" ,xc1 I ;
% .2f",xc2 I ;
% .2 f" ,xc3 I ;
%.2f",dT11;
% .2f" ,dT2 I ;
% .2 f " ,dT3 I ;
%.Zf",qlI;
%.Zf",q21;
%.2E",q3);
%.Zf",Al);
%.2f",A21;
%.Zf",A31;

case ' q ' : break;

case 'r': datostpll;
break;

default:

ClrSCrO;
printfl"\n\n\n\n\n\n\n\n\n\n\n\tTomaste
una m a l a opcien, tal vez quieras

modificar la tabla.");
qetche I I ;
datostp0;
break;
I
I
void datsstcf (void1

do( B = 1 ;
clrscrll;
printfI"\n\t\t\t*** Tabla de datos ***\n"l;
printfl"\nal Flujo de alimentaciCn. 1%3/h)
printfi"\nbl Poncentaje en peso de la disolucicn.
printf("\nc) ?orcentaje en peso del concentrado.
print f I"\nd I Tempera t u ca de al imentaci en. I SCI
printE("\nel Temperatura de vapor d e agua.(SCl
printfl"\nil Temperatura de ebulliciCn a l a salida.(ECI
pcintfl"?nql Calor especifico de la sustancia.lkcal/kqSCI
printfl"\nhl Calor de vaporizacicn del vapc,r de aqua. (kCal/kgl
printE("\nii Coeficiente global del efecto l.IW/mZhSCl
printfl"\njl Coeficiente global del efecto 2.lW/m2hSCl
printf("\n%l Coeficiente global del efecto 3.(W/rnZhSCl
printfl"\n\n\n\n\n\nLos datos de entrada son correctos.
resp = qetche0;
if lresp == In' / I resp == ' N ' I
printf("\nCu 1 dato no es correcto...
valnr=getche( I ;
switchlvalorl

..

%E",FI;

%f",xF);
%E", XCLI i
% f" ,Ta 1 ;
% "

,T v ) ;

%E",Teb31 ;
8 f",cpt i
%f",ll;

%f",Ul;
%",U21;
% " ,U31 ;

(S/N1 " I ;

ml;

case *a':
pcintf("\n\nCu 1 e s el f l u j o de alimentaciCn.lkg/hl " I ;
scanf ( " % 1f", h Fi;
break;
case 'b' :
printf("\n\nCu 1 es el poncentaje en peso de la disoluciCn.(%I " I ;
scanf l"%lf",S x F ) ;
break;
case k':
printf("\n\nCu 1 es el porcentaje en peso d e l concentrado.(%) " I ;

scan E I " % 1 I",6xCl I ;

break;
case 'd' :
printfl"\n\nCu 1 es la temperitL.ra de alimentaciCn.(CI " I ;
scenf ("E 1 f", 6Ta I ;
break;
case l e ' :
pcintf("\n\nCLi 1 es 1 3 temperatura de vapor de agua a la entrada.lCI " ) :
scanI " e 1 f", STY);
break;
case 6 ' :
prin: ("\n\nCii 1 es la temp-ratura do ebirlliciCn a la salida.15C) " ) ;
scan: ( " Y L E " , ST-t.31;
break;
case ' 3 ' :
printf("\n\nCu 1 es el c a l o r espec;fico de la solucien.(kcal/kqSCI " I ;
scan f ( " 9 1 f " ,cCp 1 ;
bKeak;
case 'h' :
printf("\n\nCu 1 es e l cabe de vaporizacien del vapor de aqua. lkCaI/i:q) "1;
sca n f ( " P 1 I",6 1 I ;
break;
case Ii':
printfI"\n\nCil 1 es el coeficiente global del efecto l.IW/mZhC) " ) :
ssanf ( " % 1 " ,S U I ;
break;
case j :
printf("\n\nCu 1 es el coeficiente global del efecto Z.IW/mZhCI " ) ;
scan E 1 " % 1 " ,6112I ;
break;
case ' R I :
printfl"\n\nCu 1 es el coeEiciente global del efecto 3.(W/mZhC) " I ;
scanf I " % 1f" ,su3 ;
break;
default :
ClrSCK( 1 ;
printf I"\n\n\n\n\n\n\n\n\n\n\r\t\tTienes que seleccionar alguna opciCn ..."I ;
qetche 1 I ;
break;
1

I
else [
if (resp == ' s ' I iresp == ' S ' ) = = O ;
else I
pcintf("\n\nCnlo responde si l o s datos son correctos o no.
qetche0;l

");

whileIB1;
C l = x F + F/x C 1 ;

Et=F-C1;
E=Et/3;
C3=F-E;
xC3=xF'F/C3;
CP=C3-E;
xc2=sc3+c3/c2;
dTt-Tv-Teb3;
dTl=(dTtl/ (l.O+(U/UZIt(U/U3)I ;
dT?=U*dTl/U?;
Tebl=Tv-dT1;
TebZ=Tebl-dT?;
dT3=Teb:-Teb3;
V=((E+ll-(C?*Cp*(TebZ-Tebll1)/1;
E2=(lE*11-lC3+Cp*ldTt-TebZ)~l/1;
E3-llE~1I-lF-Cp*lTa-dTtll~/1;
ql = v i 1 ;

qZ=E?*l;
q3=E3* 1;
A1=ql/lUCdT1~;
AZ=q?/(U?+dT?I;
A3=q3/ (U3'dT3);

10

~~

clrscr(I ;
"i ;
p r i n t f ( " \ r \ i \ t \ tT.zk-1 a d e resii 1 t ados
p r i n t f ( " \ n \ n E l floijo de l a d i s o l u c i c n d e l e f e c t o 1 e s . ( k g / h l :
p r i n r i l " \ n E l f l i i j o de l a d i s o l u c i t n .le1 e f e c t o 2 e s . ( k q / h ) :
p r i i i t E ( " \ n E l f l u j n d e l a d i s o l i i c i C n .del efecto 3 es.(kq/hl:
p r i n t f ( ' * \ r , E l concenr.ra,io t o r 3 1 .del s i s t e m a ES. ( k q / h I :
p r i r i t f ( " \ i i E l concenrrado d e l e f o c t n 1 c s .
pririr f ( " \ n E l cqnrentradci .del + - , - t o 2 e s .
p r i n t f ( " ' \ r i E l concenr-cado ,ir1 r f e c t r . 3 e s .
p r i r , t ( < ' \ r i l . . ? f r a c c i c n ?!ir CiJncentra P I P I
p r i n t f ( " \ i i l i 3 f r s - c i t n qii- c o n c e n r r a c.1 e f e c t o ? e s .
p r i n t f ( " \ n l . . i frar.-iCn q'ie c n n c - n r r a e l e f n c t . o 3 e s .
p r i n t f ( " \ r i L o , i t ? l r n ,de temperitlira - n r l e f e c . t o 1 e s . ( W ) :
p r i n t f ( " \ n l . a d e l t . 3 ,de teniperatiira en e l ? e
p r i n t ( " \ n . a 1 ~ 1 r ad e ceiiiperatiira en P I e f o
p r i n t f ( " \ n l . a v e l . ,de t r a n m i
p r i n t ( " \ n L a v e l . ,de t r a n s m i s i c o d e c a l o r e
p r I n r f ( " \ n l . a vol. ,is r r a n s m i s i c n d- c , l o r e
Frri nr f I " t r < l . a s i i p e r f i r i P rir c.a 1 entanii-rito .de
p r i n t f ("\riLa s u p e r f i c i e d e c a l e i i t s n i i e n t o de
p r i n t f ( " \ t f i L 3 s u p e r f i c i e a e c 1 r l w t d m l r n t o ,del e f e c t o 3 e s . (m2):
printf("?ri\neresi^na r p a r a vnlver.. . " I ;
p r i n t ( " \ n o p r e s i o n a q para t e r m i n a r . . . " I ;
fin = qetche(1;
switch ( f i n 1

. ..

%.?f",CIl;
%.2f",C2);
%.2f",C31;
8 . 2 f 'I, E t 1 ;
% .2 f " ,v I ;
% .2f", E2 1 ;
% 2f " ,E3 1 ;
%.2f",xC3);
% Zf",xc2 I ;
% . 2 f 8' ,xC11 ;
S . ? f ",dT3) ;
'8.2 f",dT? 1 ;

.
.

%.?f",A3);

case 'q': b ~ p i k ;
c a s e Ir?: d a t n s t c f o ;
break;
default:
clrscr(1;
prin? f i " \ n \ n \ n \ n \ n \ n \ n \ n \ n \ n \ n \ t T w n a s t e
m o d i f i c a r la t d o l a . " ) ;
getcho ( I ;
datostcfo;
hrejk;

una mala o p c i o n , t a l vez q u i e r a s

main ( v o i d )

ClrSCK( 1 ;
i n t i=8. 0 ;
texccolor(i+71:
textbacCqcnund(iA1 ;
C l KSCE (

1;

p r i n t f ( " \ n \ t \ t \ t U N I V E R C i D A D AUTrYlONA ME'."ROPOLITANA. " 1 ;

printf ("\n\t\t
UH1DP.D IZATAPALAPA." I ;
print f ("\n\t\t
DIVISION C E 1 " J ;
p r i rit f ( " \ n \ n \ t \ t \ t
LICENCIATURA E N E R G I A . " ) ;
p r i n t f I " \ n \ n \ n \ n \t
Programa d e SimulaciCn d i n mica de Tren de E v a p o r a d o r e s " ) ;
p r i n r f ( " \ n \ n \ n \ n \ n \ n E s t < . e s un proqrama r e a l i z a d o p o r : " I ;
p r i n t f ( " \ n I n g . Enrique Zamora P a n o " ) ;
prlntf("\nAsesnrado por: " I ;
p r i n t f ( " \ n M . en C . Alelandco T o r r e s Aldaco"!;
printfI"\n\n\n\n\nPresione c u a l q u i e r t e c l a para c o n t i n u a r
"I;
qetche ( I ;
ClrSCKl);
printf("\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nCualquier
repcoducciCn t o t a l o p a r c i a l d e
e s t e programa"1;
p r i n t f ( " \ n e s una v i o l a c i C n a los d e r e c h o s de a u t o r " ) ;
p r i n t ( " \ n y s e r penado por la Ley.");
printf("\n\n\n\nPresione c u a l q u i e r t e c l a para c o n t i n u a r
I:
getchel);
clrscr(1;
printf("\n\t\t
SimulaciCn d i n mica de t i e n d e e v a p o r a d o r e s " ) ;
p~int("\n\n\t\t\t\t**+
Menu * * + " I ;
printf("\n\n\t\t
11 CimulaciCn d e un evzporador de s i m p l e e f e c t o . " ) ;
printf ("\n\n\t\t
21 S i m u l a c i e n d e dos e f e c t o s en d i r e c t o . " ) ;
printEl"\n\n\t\t
31 SimulaciCn d e dos e f e c t o s en c o n c r a c o r r i e n c e . " l ;
priritf"\n\n\t\t
4 1 SimulaciFn d e t r e s e f e c t o s en p a r a l e l o . " l ;
p r i n t f ( " \ n \ n \ t \ t 5 1 SimulaciCn d e t r e s f ! f e c t o s en c o n t c a f l u j o . " l ;

...

..."

prinF("\n\n\t\t
q l Sslir del simtih3or."I;
printf ("\n\n\n\n\n\riElije la opcicn que desees simular."l;
casos = getche0;
switch (casos1

case ' I ! :
clrscr I I ;
printf ("Ciniiilacitn de un evapcrndor %<esimple efecto.");
printE("\n\n\n\t.\tLas condiciones a siiiiular son las siguientes.. . " l ;
p r t n t f 1 " \ n\ t \ t \ r \ t Da tos de 1 p raob 1ema :" I ;
printf("\nC~i I es el flujo de alimencacicn. (i:g/hl
scarif ( " ? I " , S F I ;
printf ( " C u 1 es el poncentaje en peso de la disolucien. ( % I
s ca n f ( " S 1f" , S:.: P I ;
printE("Cu 1 es el porcentaje en peso del concentrado. ( $ 1
scanf ( " a 1 E", S X C l ;
print("Cu 1 e s la presicn del vapor de agua.(atml
scsnf ("31E", Spvl;
printP("Su 1 e s la presiCn del evaporador. (atml
scan f ( "I1 " ,spe 1 ;
printP("Tii 1 es l a temperatura de alirrientaciCn. 15CI
s c a n f ( " % 1 f" ,6Ta I ;
printf ( " C u 1 es la temperatura de ebulliciCn del agua. ICl
scan ( ''e 1 f",6TeI ;
printf("Cu 1 e5 la temperatura de ebullicien de la disolucien.Cl
scan f 1 " % 1 f" , STd 1 :
printE('<Cu 1 es la temperatura de conrfensaciCn del agua.(SC)
scan ( " % 1 f " ,6Tc I ;
print("Cu 1 es Is entalpia de alimentacicn.lkCal/kgl
scanf ( " % l f " , S H l ;
printf(*'Cu 1 e s la entalpia de la disoluciCr1 concentrada.(kCal/kq)
scanf1 " % 1 f",SHc) ;
printf("Cu 1 es la entalpja del vapor de aqua.(kCal/kgl
scanf ("%lf",SHvl;
printf1"Cu 1 e s el calor d e vapori7acien del. vapor de agua.[kCal/kgl " ) ;
scan E ( " % 1f" , G 1 I ;
"I ;
printf1"Cu 1 es el coeficiente glohal. (W/mZhJCl
scanf ( '"t 1 f",&IJi ;
printf("\n\n\nPresiona cualquier tecla para continuar...");
qetche i ;
datos ( I ;
break;
case ' 2 ' :
clrscrl ) ;
print("SimiilaciCn de dos efectos en direct,>.");
princE("\n\n\n\t\tLas condiciones a simular son las siguientes
I;
print("\n\t\t\t\tDatos del pcoblema:");
"I;
pcintf("\nCii 1 es el flu20 de alimentacien. hq/hl
scanI " 9 1 f " ,6 F I ;
printf("Cu 1 es el poncentaje en peso d e la disolucien.l%l
'*I ;
scanf("%lE",br.Fl;
printf("Cu 1 es el porcentaje en peso del concentrado.l%l
"I;
scanf ('*%lf",GxCl;
printf1"Cu 1 es la presiCn del vapor d e agua.(atml
I;
scan f I " % 1f" , 6pv I :
printf1'Cu 1 es la presiCn del evapcadoc.latm1
"I;
scanf ( " % f",
I &pel i
printf ('Tu 1 es la temperatura de a1imentac:iCn. (SCl
"I;
scanf ("%if",6Tal;
printf1"Cu i es la temperatura de condensacicn del aqua.lCl
"I;
scan( " % 1f ",6Tc I ;
printE("Cu 1 es la temperatura de ebullici<:n a l a salida.lCl
"I;
scanf ("%lf",sTe2l;
printf("Cu 1 es la temperatura de vapoprzaciCn del agua.(Cl
");
scan f i " % 1f" ,6Tv I ;
prntf("r3u 1 es la entalpja d e alimentacitn. (kCal/kgl
"I;
scan E ( "'3 1 f " ,s H fI ;
printf["Cu 1 es la entalpia del vapor del <efecto Z.(kCal/kgI
,' I ;
scanf ( "%If",ShCl I ;
printf1"Cu 1 es la entalpja del vapnr de aqua.lkCal/kql
"I:

..."

I*

12

:Iriexo ,.l. Chfigo (le sintulncibn de tren de evoporticlores

scan E ( " 9 1 E", 6 Hv I ;

printE("Cii 1 es la entalpia de la disolucitn concentrada, (kCal/kg) " I ;
scan f I " 9 1 f " , 6 H ElI ;
printf1"Cu 1 es el calor .de vaporiLciCn del vapor de agua.(kCal/kgl " I ;
sco ri E ( " 1 E", 6 1 I ;
printE("Cu 1 es el coeficientp q l o b i l -del efecto 1 . (W/m2hC)
"1;
scan f ( " ? 1 f " ,6
I
.
I 1;
printf("Cu 1 es el coeficiente q l ~ ~ t . 3,del
1
efecto 2. (W/m2hC)
-1;
scinf("9lf",cU2);
printf("\nPreniona cualqbiier tecla p i r . 3 continuar.. ."I;
getche(1;
datosds ( ) ;
break;
case ' 3 ' :
clrscr1 I ;
printE("5irni~l~~iCn
de dos efectos en c,,ntracorriente."l;
printf("\n\n?ri\t\tl.as condiciones a simbular son las siguientes.. . " I ;
pr i n t f ( " \ n\ t \ r\ t \ t D.3 t?s d i 1 p mt.1ema :" I ;
printEf"\nOi 1 es el Elu]o de aliinentaiiCn. (kg/hl
I;
scan f ( " % 1 E " ,L F I ;
princf("Cu 1 es el poncentaje en peso de la disolucitn.(%l
"I;
scan f ( '"4 1E", sxF1 ;
printf1"Cu 1 es el porcentaje en peso del concentrado. 1 % )
,' I ;
scan( ''%1 f'* , & x C ) ;
printf("Cu 1 es la presiCn del vapor de aqua.(atml
I;
scanf ("%lf",~ p v l ;
");
printE("Cu 1 es la premien del evaporarlor.(;itml
scan f ( " t 1 f" , spe 1 ;
printf1"Cu 1 es la temperatura de alinientac-.en. (SC)
"I;
scanf ("%1E",6Ta);
"1;
p~intE1"Cu1 es la temperatura de condensactcn del agua.(C)
scanf ( " % 1 E", 6TcI ;
printf("Cii 1 es la temperatura de ebulliciCn a la salida.(Cl
** I ;
s c 3 n f I " Ir 1 f " , 6 Te? I ;
printE("Cu 1 es la temperatura de vapoprizacitn del aqua.(Cl
'* I ;
scan f ( " 9 1 i" , 6Tv I ;
printE1"Cu 1 es la entalpia d e alimentaciCn.(kCal/kg)
-1;
scan f 1 " O 1 E" , 6 H f I ;
printf1"Cci 1 es 1 3 entalpia del vapor del efecto Z.IkCal/kgl
"1;
scan1 "'f 1 E", 6hC1) ;
prinrrE("Cti 1 es la entalpia del vapor de agua.lkCal/kqt
"1;
scan f ( " 9 1 f " ,6 Hv I ;
printE("Cu 1 es la entalpia de la disolucicn concentrada.(kCal/kgl
"I;
scanf("Bif",6HEl);
printf1"Cu 1 es el calor de vaporizaciCn del vapor de agua.lkCal/kgl " I ;
scan( "9 1 E", 6 1 I ;
printf1"Cu 1 es el coefieiente global del efecto 1.lW/m2hCI
-1;
scan f 1 " ? 1E " , 6U 1 ;
printf("Cti 1 es el coeficiente global del eCecto 2.(W/mZhC)
"1;
scan f ( " % 1 f'*, 6 U 2 1 ;
printff<'\nPresionacualquier tecla para continuar.. . " I ;
getrhe ( I ;
datoslp( I ;
break;
case ' 4 ' :
ClKSCK( 1;
printE("Simu1aci~nde tres efectos en contraflujo."l;
prlntf["\nPara resolver este problema tomemos dos consideraciones");
printf("\nl.- que el calor de vaporizaciCn es el mismo para cada etapa");
printf("\n2.- que podemos conocer el cP de la soluciCn a concentrar");
printf("\n\n\n\t\tLas condiciones a simular son las siguientes...");
printfl"\n\t\t\t\tDatos del problema:");
printf("\ncu 1 es el flujo de alimentacicn. (kg/ht
"I;
scanf I '*% 1 f", 6 Fl ;
printE1"Cu 1 es el poncentaje en peso d e la disoluciCn. ( % I
-1;
scanf 1 " % 1 " , 6 : ~ F I ;
printf("Cu 1 es el porcentaje en peso del concentrado.(%)
"I;
scan E I " % 1 f" , 6 xC3 I ;
printf("?u 1 em la temperatura de alimentacien. ISCI
-1;
scan f ( " % 1 E", 6Ta I ;
O'

13

p r i n t f 1 " T i i 1 e s l a t m p r a t u r a ,del v i - . p ~ rd e aqua a 1s e n t r i d a . 1 5 C )

"I;
scan L ( " i. 1 f " , CTv 1 ;
p r i n t f 1 " S i i I e s l a t-mpe:atura
de e b u l l i c i e r a ?a s a l i d a . ( S C )
I;
scan f (
1 f " , GTeb3 1 ;
printf1"rl.i 1 P S P I c31nr e s p e r i f ; m ,-io la s c , l u c i C n . lkC.3l/%q5c)
"I;
s c o n f I '*Q I f", Cp);
p r i r t t f ( " r i i 1 c s e l a ~ i 1 c . r <i- v . 3 p 3 r i i 3 c i C n 1 ~ 1vapor de aqua.(kCal/kgl " I ;
scan f ( " p 1 f", 6 I ) ;
p r i n t f ( " c 4 ~1 e s - 1 c o e f i r i ~ n t er j l n t i + l del n f e r t o I. iW/nilhSC)
"I ;
s c a n I " u L C" , I I ) ;
p r i n t f i " r i i 1 e.5 e l c n e f i c i e n r i qI,:.tml d e l e f e c t o 2. (W/m?hSC)
-1;
scan f I "'+ I f " , 61J2 I ;
p r i r i t f ( " c u 1 e s e l c o * f i c i s = n t e q l * - . b 3 l d e l e l e c t o 3. (W/m:hSCI
-1;
scan f ( " 1 i f " ,hI!3 I ;
p r i n t f l"\ri\ri\n\nPrrsicna c u a l q i i i e r t e r l a para c o n t i n u a r . . . " I ;
gecche( I ;
datosCcf0;
break;
case ' 5 ' :
clrscrl t ;
printCi'<.'irniili;iCn
,de t r e s e f e c t o s en . i i r e c t . o . " ) ;
p r i n t ( " \ n P a r a r e s o l v e r e s t e problema tomemos dos c o n s i d e r a c i o n e s " ) ;
p r i n t E ( " \ , r i l . - que e l c a l o r de v a W r i z a c i C n t?s e l m i s m o para cada e t a p a " t ;
p r i n t f ( " \ n Z . - que podemos conoecer e l C P de :.a s o l u c i C n a c o n c e n t r a r " ) ;
p r i n t f ( " \ i i \ n \ n \ i \ t L a sc o n d i r i o n e s u s i i n u l a r son l a s s i g u i e n t e s . . .'O;
p r i n t f 1 " \ n \ t \ t \t \t Da tos d e l p m b l emu : " 1 ;
p r i r ~ t f l ' * \ n C u1 e s e l f l u j o de s l i m e n t a c i t n . : t q / h )
*- I ;
s c a n f ( " % 1 f " ,6 E' I ;
p r i n t f ( * ' C u 1 e s e l p n c e n t a j e en peso de l a d i s o l u c i t n . 1 % )
"I ;
scanf l " ~ l f " , 6 x F l ;
p r i n t f ( " T i i 1 e s e l p o r c e n t a j e en peso d e l c o n c e n t r a d o . ( $ 1
** I ;
s c a n f ( "2 1 E", 6 x C 1 I ;
p r i n t f i " c u ! e s 1 3 temperatura de alimentacrCn.lSC1
"):
scanf("%lf",STaI;
p r i n t f 1 " C u 1 e s l a temperatura d e l vapor de aqua a l a encrada.(SCt
"I;
s c a n f ( "P 1 f '*, 6Tv t ;
p r i n t f ( > ' S i i 1 e s l a temperatura de e b u l l i c i C n a l a s a l i d a . ( C l
"1;
s c s n f l " r 1f", 6Teb3) ;
p r i n t l ' q r i i 1 e s e1 c a l o r e s p e c i f i c o de l a s o l u c i C n . (kCal/kgCl
"1;
, sc~ntI"~lt",sCpl;
, p r i i i t f i " c u 1 e s e l c a l o r de v a p o r i z a c i c n de! vapor de a g u a . l k C a l / k q l " t ;
scan f ( " d 1 E , 6 I 1 ;
p r i n t ( " ' i ~ i 1 e s e l c o e f i c i e n t e g l c , b a l d e l e f e c t o 1 . iW/mZhCl
I' I i
scanf 1"3if".lJ);
p r i n t f ( " 0 i 1 e s el c o e f i c i e n t e g l c b a l d e l e f e c t o 2. (W/m?hC)
"ti
s c a n f I "b 1 f " ,& U ? I ;
p r i n t f ( " C u 1 e s e l c o e f i c i e n t e g l o b a l d e l e f e c t o 3.(W/mZhSC)
,*I ;
scan f i "'i1 f" , U3 I ;
p r i n t f ( " \ n \ n \ n \ n P r e s i n n ac u a l q u i e r t e c l a para c o n t i n u a c . . . " I ;
gecche I 1 ;
datostrfl);
break;
c a s e '3-1': b r e a k ;
default:
clrscrl);
p~intf("\n\n\n\n\n\n\n\n\n\n\n\tTc~tnaste una mala o p c i h , t a l vez q u i e r a s e l e g i r
nuevamente. ** I ;
qetche ( I i
/*Menu( I ; * /
break;

return O ;

ALiHexo 3%
Propiedade, termodinmicas del vapor de agua.

Anmo B. Prooiedndes termuriinniicas del vapor de ogtto.

Diagrama de Mollier

O.nt16228 . 1,0002 ' 206.3

O.OO%XYO 1 . W O O
147.2
0.01 251 3 1.0004 106.4
0.017376 1,0010
77,953
0.02383
1.0018
5734

5
10
15
20

(1.n 1229

25
IO
35

45

0.04 325
0.05733
0.0,520
0.09771

50

0,1257x

55
6o

0.16USi
0.2031
0.2555
0.31;;

10

65

;o

75

80
x5
90
Y5

1 O0
105
1 10
11.5

1LO

125

115
1.10
1-1s

4.854
155
S,540
160 . 6,302
7,14fi
165
170
8,076

t T .i ' 9.101

180 i 10.:25
I 11.456
18.;
190 j 12.XO11
19

4.1 36
3,4 LO
2.830
2.ifil
1.981

0.2425

1,673
1.119
1,210
1 ,O 36
0.8914

0.3121
0.3255

114.265

I.O'IY6

1.04i5
I .01;.1
. 1.15IS

. 1.0.<58
! 1.0603

! 1.0650
1.0697
1.074 6

1.O798

1,0850

'

608,2

0.1679
0.1833 j
0.1984
0.2131
0,22t(O

25.02
30,OO i
34.99 t

12,oi

1.0321

1 so

603,5

1.0121
1.0135

. 1.0159

3.192
3.68s
4.237

1 s,o4
20.03

O. 1 524

1,0171
1,0199
1.0228

'

2,0970
2.0697 .

'

2,0436 f
2,0187
1,9917 '
1,9718
1,9498

0.5B94
0.7149
(1.86 1 Y

2.3666

10.04

O.US76
0.1042

1.0258
1 .O290

1 ?U ; 2.7544

0.0161
0.0536
I 0,0708
,

O
' 597,2
5.03
599.4

2.1865
2.1 55 I
2,1251

4 3,4I
32.93
25.25
1Y.55
15.28

0.3931
O.iB29

1.03323
1.231 8
1,4609
1.72IY
2.0247

. 0.0182

&In10
i ,0044
1.0061
1.0079

I .O099

1 .O906
1,096 3
1,1021
1.1082
1.1144

1.1210
1,1275
1.1345
1.1415
1.1490

9.584
7.682
6.206

5.049

1.7701
0,6680
0.5817

o.5ow
0,4459

a,3924

0.1205
0.1 366

0,2567
0,2708
0.2848

0.2985

0.33R?
0,3519
0.3647

0.3775
0.3901
0.4026
0.4 150
0.4272

6k0,4
612.5

3,9287 ,
1,9085
1.8891
1.8702
1,8522

49,95

54,94

614.0
621.0

59,94

623.2

64,93
69,9 3

625.2

1.8349
i ,8178
1,8015
1.7858
J,7708

74,94

79,95 .
tit,%
89,98
95.01

1.756 1
100.04
1,7419
105.08
1,7282
110012
1,7150
115,IR
1,7018 , 120,3

1 .ti895
1.6772
1,6652 '

1,6f39

'

1,6428

'

125,3
130,4
135.5
140.6
145.8

0,1756
0.4874

1.5914

171.7

0,2166
0,1939
0.1739
O. 156-1
0.1410

0.4991
0.5107
0.5222
0,5336

1,5818
1.5721
1,5629

176,9
182.2
187.5

0.5449

6060

614,7
6163

8,2724
0,2426

0.4395
0.3516
O , 4 637

601.6

39,98 ,
44,96 ;

1.6320
1.6214
1,6112
1,601 2

1.3464
0.3068

1.5538
1,5448

150.9
156.1
161.3

r~i.5

627.3
629.3
631.3
633.2
635.1
637.0

63%9
690,7
642,5
644,3
646,0

6472
'

649,f

650.8
,

.
*

652.5
6 5-1.O
655.5
656.9
658,)
659.6
660.9

I
*

;Ei

662.1
663.2
664.3

665.3
666.2

.-inem . Propierlodes erriiodinhitiicns del iwpor <IC osun.

Concliisin )
I'ic6:on dc
sniiitdclon

Kx

Ern:

15,857
17,SRS
19.456
21,177
23,659

1.1565
I . 16?S
1.1726

26.007
28.531

1.1991 f

1,1812

0,1273
O. I IS 1
U.1043
(1,09472

(1.oa614

1,1900

0,5562

0.60io

1,5358
I,5270
1.5184
1.5099
1.5012

0,6120
0.6229
0.6339
0,644R
0.6558

1.4926
I ,4R40
1.4755
1 :lh69
1.4584

0,6667
0.6776
0,6886
0,6994
0,7103

1.4499
1,4413
1.4327
1,4240

0,5675
0,5788

' 0,5899

214.3
219.8
225,3

668.3
468.8
669.3
I

0.07845
1.2fiR8 i 0.07151
1,2186 ! 0,06530
1,2291
0.05970
1.2400
0.05165

31,239

34.140
3 7.244

I,2S12

40.56
34,lO
47,R7
56.14

0.05006
0.04591

1.%29
1.2755
1.2888
1.302 1

Si,H8

'

O.0421.3
0,03870
0,03557

1.4066
I . 7978
t. m n
1.379;
r,.:W2
( ~ 6 5 1.3706
~

65.46
70.54
75.92
8 1.60

1.3 169
1, 3.I 2 I
1.1484
1.35i
1. I837

0.03272
0.0 30 1 D
O.O?;iI
0.02552
0.02350

li.7212
0.7 72 I
f1.7.1 I1

87.61
93.95
100.64
107.69
11 5.1 1

1.4036
1.425
1.448
1,472
J .499

0.021 6 3
0.0 1991
0.01830
0.01682

0.7767
0.7880
0.3994
0.81 I O

0.01 5 1 4

0.8229

122.95
131.18
139.85
118.96
1.58.54

1.529

0,01415

0.8351
0.8476

60.66

1,642
1,692

0.01076
0,009759

1.747
1.814
1.YCii
2.03
2,23

168.63
179.24
190,42
202.2 I
214.68

217.3
2 19.9
222.5
225.2
22 5.4

;;;;; j

0.01295
0.01183

I 2.30
: 2,38
! 2.50
j

2.79
?,O7

o.oaaso3i 0.9015
0.006%3!

0.00606

1.3613
1.3516
1.3415
1,3312
1,3206

1.3097
I .2982
1,2860
1.2728
1.3586

0.8604
0.8735
0.8871

! 0.007h75 1 0.9173
!

1.II 5'1

0.00500

0.00476
0.00450
0.00418 '
. 0,00365 .
. 0.00314 ,
i

'

0,9353 f
0.9553
0.9842 '

'

0.992
1.002
1.011
1.04

1.2433
1.2263
1,2072
1.1833
1,1506

. 1.142

1.058

1.132
1.116
1.08

259,2
265,o

271.0
277.0
283.0

669,O
66a,4
667.8
666.9
665.9

Propiedades termodinmicas del vapor de agiia saturado.

( entrada por presiones )
Prowl 6
lUiad&m

XC/cm

0.01

28.64

0.05

32,55

0,06

0.10

35.82
41.16
45.45

O,] 2

49,06

0,15
0.20
0.2s
0,30
0,35
0.40

53,60
59,67

64 ,S 6

OJO

110.86

0,60

81.45

0,70

g9.45

1.4
1.5

108.74

1.8

2?0

22
2.4

2.6

58
3,O
312
3.4

3.6
3.8
4 .O
4.5
5.0
5,s
6.0
6.5
7.0

24.19

18.4 5
14,95
1260
10,21
7,795
6.322
5.328
4.614
4.069

3.301
2.78 3

2,409
2.1 2s
1.904
1.725

92.W
96.1 8
99.09
101.76
104.25

1.6

W,64
68,27
5 5,28
46.53
35,46
28.73

68.68
72,24
75.42

0,W
1.o
1.1
1.2

L,3

-?
---

0.04

0.80

I
-

0.02
0,025
0:33

0.08

de-

6.70
12.74
17,20
20,7a
23.f7

0.01 5

T-rrr.

1,S78

1.455
1,350
1,259
1,180
1,111
0,9952
0,9016
0.8246
0.7601

irrb.s6

f 10.79
112.7 3
116.33
119.62
122.65

125.46

0,@24J
0.04>7
0.06 12
0,07 35
0.08 36
0.0 99 I
0.1 I26

O,L232

a, t 4Q2

O,LS38
0.f 650
O, 790
O, 974
0,2120
0.2241
8.2345

132,813

l35.08

0.5804

04Qa

0,5413

0,408 L
i1.4130
0.4 i76
0.4 22 I
0.4 32
0.4422
D.45 1o
0.4591
13,4666

137,18
139,18
141.09
142.92
147.20
151.1 1
154.71
158.08
161.15

IM,i7

0,5196

0.4939
0.4706
0,4213
:

I
a

0.3816
0 . 3 4 ~

0.3213
0.2980
0.2778

2,0655

2.0499

2.0253
2,0064

1.3973

(3.4737

21.65

1.ipO8
1,9664
1,9478
1.9326
1.9140

1.8903

59.61

1,8718

M,49
68,61
7517
f5.36
80.81

1,8267

1,8436
1,83 34
1,%ISO
1,4001
3 m
2
a i ~ m 1,7767
( 1 1 ~ 0 1 8 1,7673
1,7587
CI, 3096
1.7510
I.3L68
1.744-0
0,3235
1,2375
0,3291
(1,3354 1,7395
1.7260
i.3408
1,7209
0.3459
1.711s
435%
1,7029
0.3638
1,6952
a3715
ii.37U6
1.6884
1,6819
0.38!53
1,6759
1.3914

0.7052
0.6578
0.6 166

128.08

Wr=.r

2,1447
2,1096
2,01147

0 3 43 7
Om92
41,2721
(112832

130.55

1,6703
3.665c

I ,660 I
1,6557
1,6514
1,6474
1.638E
1.629;
1,6195
1.6151
1.608j
1,6025

ar;ri

623.1
625.1
626.8

628.2
629,5

631.6
633:4
634,9
636;2

89.43
9239
96.19
99,12

637.4

101,81

638.5
639.4
620.3

11092

642.8
43,5

io.t.32

112.89
116.54
11987
122.9
125.8 .

644.7

13i,Q

645.6
646.8
647.8
648.7
49.5

135,6

6M,9

128.5

650,3

133.4
137.8

139.8
141.8
143.6
148.0
152.1
155,s
159.3

$62.6
165.6

651.6
6522
652.8
653.4
659.7
655.8
656,g
6S7.8
658.7
659.4

( Conclusin )

KP

cm7

... P

. .. .

7.5
8 ,O

I66.96
169.61

8.5

9.0
9.7

10
I1
12
13
14
15

lh

!
!

.
'

17

18

19

20
32

24

26

30
32
3.1

36

3s
,

42
43
36
48

SO
5s

60

65
70
75
80
,

85

l'4,.!3

197.36
2110.43
203.3.5
106. I 4
Z08,8 I

211.38
216,23
220.75
223.99
236.35
239.7i
%.I't
,a5
2.16.1 7
249.18
252.07
254.RT
257.56
260.1 7
262.70
268.69
274.29
27934
284,48
289,17
293.62
297.86

301.92

90
95
1a0
I10

120
1 10
140
\ ?O
1 hrl
1x0
_'o0
224
225.4

lY(l.SI

225.98
232.76

2R

40

172.11
t 74.5 1
176.82
Ii9.01
183.20
I8i.W

30 5.80
109,53
116,.58
123.15

329.30

'

'

315.09
3-1056
345.7-1
?55.35
364.08
37 1.6
17 1.1

--.-

)el liquido

..

0.2602
0.2448
0.2 11 I

n.4803
O.486S

0.21~19

0.4980

0,2080

O.198 1
0,1808

0.1664
0,1541
0,1435
0.1 343
O. 1262
0.1 190
o.1 126

Dcl

0.4923
0.503 3

o, 4 5 s

varxir

I"

168.5

1.sxo

181.2
1115.6
t R9.7

1,5592
1.5526
1.5464
1,5406
1.53.51

0.5503
0,5572

A,

1,5974
3.5922
I.SKI
1.5827
1.5782
1.5661

0.IXU
0,SZiY
0,5352
0.5930

<-

Dri Ilquido

0.06375

0.05995
O.OS65X
0.0535 1

o.o so7 8
0,0482I
0.0 460 1
0.04393

I).U4201
0.04 o ::1

173.9
176.4
178.9

661.0

66 1.9

664.7
665.4
666.0
666.6

191,s

197,l
200,6

667.1

203.9

667.5

667,9
661.2

0,5928
0.6026
.6120
0.6290
0.6168
0.644 3
0.65 15
0.6584
O .6 649
0.671 2
0.6773
0.6832
0,6889
! 0,6944
j 0.7075
0.7 196
0.7311
0.7420

,4793
,1732

1 239,5
243.6

668.5
668.9
669.3
669.5
649.6
. 6Y.7
I 669.;

,167I
,4617
,45164
,4513
.4463
.4415
,4 369
,432.7
1,4286

247.5

1,5160 ! 21'5.8
1.5078
221,2
1.5060
220.1
1.4926
230.8

0,5820

251,2
254.8
258.2
261.6
264,9
268.0
271,2

1,4176
1,4078
1.3985
1,3897
m i 3
1.3i3i

66V.h

669.5
669,)
!

274.2
281.4
288.4
294.8
300.9
0.7524
307,O
. 0.7623
312.6
1;3654 , 31 8.2
0.78 O
1.3576 : 32?,6
0,02096
0,019-l
1.3500
328.8
0,7898
1.3424
314,O
0.01 845
0.7983
0.R I47
113279
344.0
0.01637
1 .3 1 3 8
35 1.9
0.01 462
363.0
0.01 3 12
1.299H
I.ZX58
372.4
0.01 181
3R1.7
1.2713
1.2544 ' 390.1
! !?;I
'
410.2
431.5
I
i,ii83
!
1.01
. 47A
505
0.00314
1.058

0.03616
0.03 310
0.03033
0.02795
0,02Sli

*"

660,2
660.8
66 I .?
662.0
662.5

171,3

o. 1068
0.1016
U.0925 I
0.08492
0.07816
O .O728 8
0,06802

...

Dei vapor

669.0
668.8
668.4
668.0

667.7
667.3
666.2
665.0
661.6
662.1
660.5
658.9
657.0
655.1

653.2
651.1
646.7

64 I .9
636.6

63 I .O

624.9
61 8.3
60L.S
582.3
572