Trabajo de Fin de Grado

Grado en Ingeniería de la Energía

Simulación de instalaciones frigoríficas usando

ASPEN-HYSYS

Autor: Lucía de la Rosa Pineda

Tutor: David Velázquez Alonso

Equation Chapter 1 Section 1

Dpto. Ingeniería Energética

Escuela Técnica Superior de Ingeniería
Universidad de Sevilla
Sevilla, 2018
 Trabajo de Fin de Grado
Grado en Ingeniería de la Energía

Simulación de instalaciones frigoríficas usando

ASPEN-HYSYS

Autor:
Lucía de la Rosa Pineda

Tutor:
David Velázquez Alonso
Profesor titular

Dpto. de Teoría de la Señal y Comunicaciones

Escuela Técnica Superior de Ingeniería
Universidad de Sevilla
Sevilla, 2018

iii
 Proyecto Fin de Carrera: Simulación de instalaciones frigoríficas usando ASPEN-HYSYS

Autor: Lucía de la Rosa Pineda

Tutor: David Velázquez Alonso

El tribunal nombrado para juzgar el Proyecto arriba indicado, compuesto por los siguientes miembros:

Presidente:

Vocales:

Secretario:

Acuerdan otorgarle la calificación de:

Sevilla, 2018

v
El Secretario del Tribunal
 A mi familia y en especial a mi
madre.

vii
 Agradecimientos

Me gustaría agradecer a mi familia por toda su ayuda y apoyo incondicional, a Ismael por sus consejos a la
hora de redactar esta memoria, a mis amigos del grado por estos 4 años, sin ellos no podría haber llegado
donde estoy y a mi tutor, David Velázquez.

Lucía de la Rosa Pineda

ix
 Resumen

La aplicación de simuladores en la industria es cada vez más frecuente. Programas como HYSYS son de gran
utilidad para ingenieros, resulta imprescindible conocerlos. Para poder emplear todos los recursos del
programa se requería mucho tiempo, pero es posible aplicar los conocimientos básicos en gran cantidad de
procesos. En este trabajo se pueden adquirir dichos conocimientos básicos, definición de los componentes, los
métodos termodinámicos empleados, los equipos o la creación de estudios de sensibilidad.
Durante la elaboración de 3 casos distintos, HYSYS ha permitido calcular sus propiedades de forma muy
precisa y analizarlos para distintas variaciones. El primero de los casos ha sido un sistema simple con un único
compresor y evaporador. El segundo y el tercero contaban con dos evaporadores y dos compresores, su
diferencia recaía en los tanques. El segundo contaba con un único tanque del que partían todos los flujos y el
tercero con dos tanques a dos presiones distintas.

xi
 Abstract

The application of simulators in the industry is more and more common. Programs like HYSYS are useful for
engineers and they really need to understand them. A lot of time would be needed in order to be able to use all
the resources of the program; however, it is possible to apply the basic concepts in many processes. In this
project these basic concepts are explained: components are defined, as well as the thermodynamic methods
that are employed and the development of sensibility studies.
For three different scenarios, HYSYS has allowed us to calculate their properties in a very precise way. It also
let us analyze these scenarios when some conditions varied. The first scenario was a simple system consisting
of one single compressor and evaporator. Second and third scenario had two evaporators and two compressors,
and their difference was the number of tanks that they had. The second scenario had just one single tank while
the third scenario had two tanks at different pressures.

xiii
 Índice

Agradecimientos ix
Resumen xi
Abstract xiii
Índice xv
Índice de Tablas xvii
Índice de Figuras xix
Notación xxi
1 Introducción 1
1.1 Simuladores de procesos 1
1.2 HYSYS 1
2 instalaciones frigoríficas 3
2.1 Componentes 3
2.1.1 Evaporador 3
2.1.2 Condensador 3
2.1.3 Compresor 3
2.1.4 Válvula 3
2.1.5 Refrigerantes 3
2.1.6 Tanques 4
3 Métodos termodinámicos 5
3.1 Introducción 5
3.2 Tipos de métodos termodinámicos 5
3.2.1 Ecuaciones de estado (EOS) 5
3.2.2 Métodos coeficientes de actividad (LACM) 7
3.3 Elección del método termodinámico 7
3.4 Cálculo de entalpía y entropía 9
4 Hysys 11
4.1 Propiedades 11
4.1.1 Componentes 11
4.1.2 Paquete termodinámico 12
4.2 Simulación 13
4.2.1 Caso principal 13
4.2.2 Casos de estudio 17
5 Casos Simulados 19
5.1 Caso 1 19
5.1.1 Enunciado 19
5.1.2 Solución manual 20
5.1.3 Solución HYSYS 23
5.1.4 Comparación de los resultados 25

xv
 5.1.5 Análisis de sensibilidad 26
5.2 Caso 2 28
5.2.1 Enunciado 28
5.2.2 Solución manual 29
5.2.3 Solución HYSYS 31
5.2.4 Comparación de los resultados 33
5.2.5 Análisis de sensibilidad 34
5.3 Caso 3 40
5.3.1 Enunciado 40
5.3.2 Solución manual 41
5.3.3 Solución HYSYS 43
5.3.4 Comparación 45
5.3.5 Análisis de sensibilidad 46
6 Conclusiones 49
Bibliografía 61
 ÍNDICE DE TABLAS

Tabla 2-1 Tipos de refrigerantes.3 4

Tabla 3-1 Métodos para compuestos Polares-No polares 8
Tabla 3-2 Métodos para compuestos Polares/No electrolitos 8
Tabla 3-3 Ecuaciones recomendadas para cada proceso, Introducción a la simulación en Ingienería Química29
Tabla 5-1 Propiedades de los componentes del sistema, caso 1. 22
Tabla 5-2 Desplazamiento volumétrico, potencia del evaporador, del compresor y COP, caso 1. 22
Tabla 5-3 Propiedades de los componentes del sistema, caso 1 22
Tabla 5-4 Potencia del evaporador, del compresor y COP, caso 1. 22
Tabla 5-5 Eficiencia y UA del condensador 22
Tabla 5-6 Propiedades de los componentes del sistema (HYSYS), caso 1. 24
Tabla 5-7 Potencia del evaporador, del compresor y COP, caso 1. 24
Tabla 5-8 Propiedades de los componentes del sistema, condensador de agua (HYSYS), caso 1. 25
Tabla 5-9 Potencia del evaporadot, del compresor y COP, evaporador de agua (HYSYS), caso 1 25
Tabla 5-10 Eficiencia y UA del condensador, condensador de agua, caso 1. 25
Tabla 5-11 Comparación de resultados, caso 1 26
Tabla 5-12 Propiedades de los componentes del sistema, caso 2. 30
Tabla 5-13 Desplazamiento volumétrico de los compresores, caso 2. 31
Tabla 5-14 Potencia de los compresores y COP, caso 2. 31
Tabla 5-15 Propiedades de los flujos de materia de HYSYS, caso 2 33
Tabla 5-16 Desplazamiento volumétrico de los compresores (HYSYS), caso 2. 33
Tabla 5-17 Potencia de los compresores y COP, caso 2 33
Tabla 5-18 Resultados Manual y HYSYS, caso 2 34
Tabla 5-19 Propiedades de los componentes del sistema, caso 3 42
Tabla 5-20 Caudales, caso 3. 42
Tabla 5-21 Desplazamientod volumétricos del compresor y COP, caso 3. 43
Tabla 5-22 Propiedades de los flujos de materia de (HYSYS), caso 3 45
Tabla 5-23 Potencia del intercambiador, de los compresores y COP (HYSYS), caso 3 45
Tabla 5-24 Comparativa de resultados, caso 3 46

xvii
 ÍNDICE DE FIGURAS

Figura 4-1. Components Lists 11

Figura 4-2. Componenetes 12
Figura 4-3. Métodos Termodinámicos 12
Figura 4-4. Opciones del método Peng Robinson 13
Figura 4-5 Palette 13
Figura 4-6 Elementos 14
Figura 4-7 Intercambiador 14
Figura 4-8 Esquema Simple 15
Figura 4-9. Incógnitas del sistema 15
Figura 4-10 Flujos de materia, Worksheet 16
Figura 4-11 Esquema simple 2 16
Figura 4-12 Model Summary Grid 17
Figura 4-13. Selección de variables. 18
Figura 4-14 Cotas de las variables independientes. 18
Figura 4-15 Resultados del caso estudio. 18
Figura 5-1 Ciclo simple, condensado por aire (izq) y condensado por agua (der). 20
Figura 5-2 Diagrama p-h, caso 1 23
Figura 5-3 Ciclo simpe (HYSYS), condensado por aire (izq), condensador por agua (der). 23
Figura 5-4 Datos de entrada al sistema (HYSYS), caso 1. 24
Figura 5-5 Presión de condensación. 26
Figura 5-6 Evolución de la temperatura de salida del compresor frente a su presión. 27
Figura 5-7 Evolución de la potencia del compresor frente su presión de salida. 27
Figura 5-8 UA del condensador frente su presión. 27
Figura 5-9 Presión de evaporación. 27
Figura 5-10 Temperatura de salida del agua frente a la presión del evaporador. 28
Figura 5-11 Potencia del compresor frente a la presión del evaporador. 28
Figura 5-12 Esquema 2 29
Figura 5-13 Diagrama P-h, caso 2 31
Figura 5-14 Esquema caso 2 (HYSYS). 32
Figura 5-15 Datos de entrada en de los flujos de materia, caso 2 32
Figura 5-16 Datos de entrada de los flujos de energía, caso 2 33
Figura 5-17 Línea de alta presión. 34
Figura 5-18 Potencia del condensador frente a su presión. 35
Figura 5-19 Potencia del compresor de alta frente a la presión de condensación. 35

xix
Figura 5-20 Evolución del título de vapor 4 frente a la temperatura de salida del condensador. 36
Figura 5-21 Caudal del condensador frente su temperatura de salida. 36
Figura 5-22 Potencia del compresor de alta presión frente a la temperatura de salida del condensador. 36
Figura 5-23 Potencia del condensador frente a su temperatura de salida. 36
Figura 5-24 Linea de media presión, caso 2 37
Figura 5-25 Evolución del caudal del condensador y Evaporador A frente a la potencia del evaporador A.
37
Figura 5-26 Evolución de la potencia del compresor frente a la potenci del evaporador A. 37
Figura 5-27 Evolución de la potencia del condensador frente a la potencia del evaporador A. 38
Figura 5-28 Línea de baja presión, caso 2 38
Figura 5-29 Evolución de la presión del Evaporador B frente a la temperatura a la entrada del evaporador.
39
Figura 5-30 Evolución del título del vapor a la entrada del Evaporador B en función de su temperatura. 39
Figura 5-31 Potencia del compresor de baja frente a la presión de entrada. 39
Figura 5-32 Temperatura de salida del compresor de baja frente a la presión de la entrada 39
Figura 5-33 Evolución del caudal del condensador frente a la temperatura a la salida del compresor de baja.
40
Figura 5-34 Esquema 3 41
Figura 5-35 Diagrama p-h, caso 3 43
Figura 5-36 Esquema Caso 3 44
Figura 5-37 Datos de entrada en los flujos de materia, caso 3. 44
Figura 5-38 Datos de entrada en los flujos de energía, caso 3. 45
Figura 5-39 Estudio alta presión, caso 3 46
Figura 5-40 Evolución del caudal a la salida del Evaporador A frente a su temperatura. 47
Figura 5-41 Evolución del caudal 6 frente a la temperatura de salida del evaporador A. 47
Figura 5-42 Estudio baja presión, caso 3 47
Figura 5-43 Evolución del título del vapor a la salida del intercambiador frente caudal etileno-glicol. 48
Figura 5-44 Evolución del caudal compresor de baja frente al título de vapor de salida del intercambiador.
48
Figura 5-45 Evolución de la potencia del compresor de baja presión frente a su caudal. 48
Figura 5-46 Evolución del caudal del condensador frente al del compresor de baja presión. 48
Figura 5-47 Potencia del compresor de alta frente a su caudal 48
 Notación

AU (W/K)
cpw Calor específico del agua (kJ/kg)
COP Coeficiente de rendimiento
ε Efectividad del intercambiador
hi Entalpía específica (KJ/Kg)
i Línea del circuito.
mR Caudal del refrigerante. (Kg/s)
PCD Presión de condensación (kPa)
PEV Presión de evaporación (kPa)
Q f Calor cedido en el evaporador (kW)

Q MAX Calor máximo que puede cederse en el evaporador (kW)

s Entropía específica (kJ/kgK)
TEV Temperatura de evaporación (ºC)
TWE Temperatura del agua a la entrada. (ºC)
TWS Temperatura de agua a la salida. (ºC)
v Volumen específico (m3/kg)

V Desplazamiento volumétrico del condensador (m3/h)
t

WC Potencia del compresor (kW)

ηV Rendimiento volumétrico
ηS Rendimiento isentrópico

xxi
22 Índice de Figuras
 1 INTRODUCCIÓN

n los últimos años,e el uso de simulaciones en los procesos de ingeniería ha ido en aumento. Han sido

E utilizados tanto en la industria como en las aulas de ingeniería. El uso de simuladores ha permitido en las
aulas la familarización con los equipos por parte de los estudiantes. Nos permiten conocer el
funcionamiento de los equipos y las distintas combinaciones entre ellos. La realización de este trabajo se ha
basado en la simulación de problemas académicos, pues su objetivo principal ha sido demostrar las cualidades
del programa y servir como ayuda en el futuro para futuros usuarios.

1.1 Simuladores de procesos

Mediante los simuladores se pueden conocer los resultados de los procesos antes de implementarlos en la
realidad y de este modo poder corregir errores y evitar accidentes o situaciones no deseadas. Hay una gran
variedad dentro de los simuladores de procesos de ingeniería, una de sus mayores vetajas reside en sus bancos
de datos de las propiedades termodinámicas de los compuestos. También se deben destacar los diferentes
métodos termodinámicos que poseen de esta forma el usuario puede escoger el que más se ajuste a su
necesidad, en un capítulo posterior hablaremos más sobre ellos y los distintos tipos que existen.
Las simulaciones llevan poco tiempo en el mercado, gracias al avance de la tecnología su desarrollo en los
últimos años ha avanzado mucho. No solo permite la optimización y mejora de los sistemas también se
utilizan para el diseño y cálculo de las instalaciones. En las escuelas de ingeniería se enseñan los métodos
manuales de cálculo a partir de ecuaciones y tablas, estos resultan muy útiles para la compresión del proceso,
pero en las aulas también deben considerarse las posibles aplicaciones académicas que poseen estos softwares.
El cálculo de las propiedades de los procesos se basa en modelos matemáticos, estos modelos pueden
resolverse para sistemas sencillos. Sin embargo, cuando se añaden equipos y se vuelven más complejos resulta
más fácil cometer errores. Mediante las simulaciones se pueden evitar esos errores y ahorrar el tiempo que
llevaría una resolución manual.

1.2 HYSYS
El simulador Aspen HYSYS es un software perteneciente a la empresa AspenTech, la Escuela Técnica
Superior de Ingeniería de la Universidad de Sevilla posee licencia de uso para su versión V 8.6, disponible
para todos los estudiantes en los ordenadores del Centro de Cálculo.
HYSYS es un programa de simulación para ingeniería de procesos. Permite la simulación de: procesos
químicos, refinación petrolífera, instalaciones frigoríficas y otras muchas funciones. HYSYS cuenta, además,
con un gran paquete termodinámico con las propiedades de multitud de componenetes permitiendo una
simulación muy realista de los procesos.
El simulador cuenta con una interfaz muy sencilla y va guiando al usuario en todo momento sobre los pasos
que debe seguir. A partir de unos datos de entrada puede calcular todas las condiciones de la corriente de
manera bidireccional, es decir tanto las propiedades anteriores a la corriente definida como las posteriores.
Quizá el mayor inconveniente que puede tener sea el idioma, solo está disponible en inglés, puede llegar a
retrasar a algunos usuarios, aunque el lenguaje utilizado es sencillo y fácil de entender.
El entorno de simulación utiliza bloques para representar los distintos equipos, de esta forma nos podemos
hacer una idea muy visual de como sería la instalación en la realidad. En este proyecto se ha dedicado un
capitulo a la introducción de HYSYS como guía para aquellos que deseen empezar a utilizarlo. De forma clara

1
 2 Introducción

se detallan como deben ser los primeros pasos para realizar una simulación, en este trabajo no pueda abordarse
todas las posibilidades de HYSYS. No obstante, el programa cuenta con una gran ayuda que resolverá la
mayor parte de las dudas que nos puedan surgir durante su utilización.
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 3

2 INSTALACIONES FRIGORÍFICAS

as instalaciones frigoríficas han permitido la obtención de frío de manera artificial ayudando así a la

L conservación de alimentos y a la climatización de locales. Tiene por tanto aplicaciones tanto en la

industria como en el hogar.

2.1 Componentes
2.1.1 Evaporador
Es uno de los elementos más importantes de una instalación frigorífica, es el encargado de producir el efecto
frigorífico. El fluido llega al evaporador en la zona bifásica tras haber pasado por la válvula de laminación. Al
entrar en el evaporador se pone en contacto con otro fluido, la corriente a enfriar, que se encuentra a una
temperatura superior. Se produce la ebullición del refrigerante, para lograr la ebullición el proceso requiere de
gran cantidad de energía. Esta energía es extraída de la segunda corriente disminuyendo su temperatura.

2.1.2 Condensador
El condensador transfiere el calor del refrigerante a otro medio donde no moleste. Al condensador el
refrigerante este pierde la energía que había ganado en el evaporador y puede seguir extrayendo energía al
fluido a enfriar. Los condensadores más utilizados suelen ser de agua, aire o evaporativos.

2.1.3 Compresor
La función del compresor es por una parte mantener el circuito en funcionamiento impulsando al fluido y por
otra elevar la presión del refrigerante para lograr su condensación a temperaturas más altas. Es el único
elemento junto con las bombas que consume electricidad, las bombas suelen despreciarse.

2.1.4 Válvula
Tras condensar el fluido se debe expandir hasta la temperatura de evaporación. Al expandirse el fluido entra en
la zona bifásica, parte del refrigerante se evapora. Con la válvula se completa el sistema del ciclo simple de
refrigeración.

2.1.5 Refrigerantes
Son compuestos químicos que permiten el uso de sus cambios de estado para la producción de frío. Las
características que se deben tener en cuenta a la hora de buscar un fluido son: color, olor, inflamabilidad,
seguridad, corrosión de los metales, solubilidad con el agua, detección de fugas y aplicación. En la Tabla 2-1
se muestran algunos de los refrigerantes más utilizados, normalmente no se utilizan a la presión atmosférica, se
utilizan los cambios de presión para lograr las temperaturas de ebullición deseadas. Es fundamental conocer
las propiedades de los refrigerantes a la hora de utilizarlos.

3
 4 instalaciones frigoríficas

Refrigerante R-114 R-12 R-500 R-22 R-502 R-116 R-13 R-503 R-14

Punto de
Ebullición (ºC) a 4 -30 -33 -41 -46 -78 -82 -89 -128
Presión Atm.
Tabla 2-1 Tipos de refrigerantes.3

2.1.6 Tanques
La cantidad de refrigerante que circula por el sistema varía en función de la carga que requiera el sistema.
Permite absorber las fluctuaciones del sistema y conservar el refrigerante en equilibrio líquido y gaseoso.
También se puede utilizar para recoger y almacenar el refrigerante.
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 5

3 MÉTODOS TERMODINÁMICOS

os métodos termodinámicos permiten el cálculo de las distintas propiedades termodinámicas de un

L componente a partir de su tempertura y/o presión.

3.1 Introducción
Los métodos termodinámicos van acompañados de bases de datos experimentales a partir de los cuales se
pueden hallar las propiedades deseadas. Existen bases de datos de componentes puros con parámentros
dependientes de la presión y la temperatura o bien con datos sobre los valores críticos, normal de
ebullición…También cuentan con bases de datos procedentes de la mezcla de componentes dónde se han
registrado los datos referidos a la interacción de pares de componentes.
La elección del método termodinámico adecuado es de vital importancia cuando se realiza una simulación. Si
se escogiera de manera errónea induciría a error a la simulación invalidando los resultados obtenidos. La
elección dependerá principalmente de los compuestos a simular y del rango de temperatura y presión del
sistema. Debemos tener esto muy en cuenta y evaluar si los resultados obtenidos tras la simulación son los
esperados. Para escoger el mejor sistema se recomienda hacer comprobaciones con varios hasta hallar el que
más se ajuste especialmente cuando se requiera el cálculo de entalpía y exista equilibrio de fases.

3.2 Tipos de métodos termodinámicos

Existen muchos métodos termodinámicos distintos, cada uno de ello permite el cálculo de las propiedades para
compuestos y rangos diferentes. La variedad de métodos permite ampliar las posibilidades pues debido a la
gran cantidad de componentes que existen y todas las combinaciones que podrían realizarse sería imposible
hacer un único método que abarcara a todos. A pesar de ello, los métodos si pueden ser agrupados en 4 tipos:
Ecuaciones de estado (EOS), modelos de coeficientes de actividad (LACM), modelos combinados (EOS
Predictivas) y LACM especiales para sistemas específicos.

3.2.1 Ecuaciones de estado (EOS)

Son válidas para sustancias en estado líquido o gaseoso, relaciona la temperatura, la presión y el volumen
molar de un fluido puro. Su rango de aplicación es muy amplio, para presión y temperatura superior al punto
crítico, lo que permite su uso en multitud de aplicaciones, además, permite el cálculo de gran cantidad de
propiedades. El caso más básico de aplicación de una EOS es mediante el uso de la ecuación de los gases
ideales (PV=nRT).
Aunque se basen en el cálculo de propiedades para sustancias puras pueden utilizarse en aplicaciones con
mezclas de sustancias. Siguiendo las Reglas de mezcla realiza los promedios de los distintos compuestos que
la forman.
Los métodos Peng Robinson (PR) y Soave-Redlich-Kwong (SRK) utilizan tres parámetros para los
componentes puros y uno de interacción binaria para las parejas de componentes en el cálculo de sus
propiedades. Son las más utilizadas para las mezclas líquido-gaseoso

5
6 Métodos termodinámicos

- Ecuación de Peng-Robinson: son útiles para la estimación de las propiedades en equilibrio líquido-
vapor. Se pueden utilizar tanto para presiones y temperaturas moderadas como altas (T>-271ºC y
P<100000 kPa). Sus resultados son similares a los de SRK, pero mejora en la predicción de refiones
cercanas a la crítica y en las densidades de líquidos, en especial con los apolares.

RT α ai (3-1)
=P −
V − bi v(v + bi ) + bi (v − bi )

2
  T 
α=
1 + (0.37464 + 1.5722ω − 0.2699ω ) 1 −
2
 (3-2)
  TC  

R 2TC (3-3)
ai = 0.4572
PC

RTC (3-4)
bi = 0.07780
PC

PC: Presión crítica

TC: Temperatura crítica
R: Constante de los gases
V: Volumen molar
C1 y C2: Funciones empíricas
α: Función del factor acéntrico de Pitzer (ω). HYSYS corrige la ecuación cuando el factor acéntrico es
mayor a 0.49, (3-5).
2
  T 
α=
1 + (0.379642 + 1.48503ω − 0.164423ω + 0.016666ω ) 1 −
2 3
 (3-5)
  TC  

- Ecuación de Soave-Redilich-Kwong (SRK): se utiliza para mezclas débilmente no ideales en

hidrocarburos con presiones y temperaturas moderadas (T>-143ºC, P<35000 kPa), ya sea en fase
líquida o vapor y también para el nitrógeno, monóxido de carbono, dióxido de carbono y sulfuro de
hidrógeno. Para presiones y temperaturas cerca de la crítica no se recomienda por no ser demasiado
exacta, en cambio para presiones y temperaturas superiores sí se aconseja su uso.

RT αa
=P − (3-6)
V − b v (v + b )

2
  T 
α=
1 + (0.48508 + 1.55171ω − 0.1561ω ) 1 −
2
 (3-7)
  TC  

R 2TC
a = 0.42747 (3-8)
PC

RTC
b = 0.08664 (3-9)
PC
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 7

3.2.2 Métodos coeficientes de actividad (LACM)

Para mezclas con comportamientos no ideales en fase líquida se desarrollaron los Métodos coeficientes de
actividad. Estos métodos permiten suplir algunos de los que no pueden ser abordados por las EOS, que en fase
vapor se asumen como gases ideales o equilibrios líquido-líquido. Ofrece ecuaciones algebraicas para los
coeficientes de actividad a partir de la composición y la temperatura. Debido a que los coeficientes de
actividad son parámetros para sistemas ideales no se recomienda su uso en situaciones supercríticas ni para
incondensables. Requieren de parámetros de interacción binaria que en ocasiones son díficiles de conocer o
bien están límitados a unos rangos de presiones y temperaturas.
Estos métodos se basan en la energía libre de Gibbs de exceso, esta es la diferencia entre la energía libre de
Gibbs ideal y real. La existencia de distintos métodos se debe a las diferentes formas en las que se puede
modelar la energía libre de Gibbs.
Algunos de los métodos son:
- Wilson: se utiliza para mezclas de varios componentes en equilibrio líquido-vapor. Es la primera
ecuación de actividad que utiliza el modelo de composición local para derivar la expresión del exceso
de energía de Gibbs. Es bueno prediciendo el comportamiento de mezcla de multiples componentes
en equilibrio gas-líquido, por el contrario este método falla para sistemas con equilibrio de dos
líquidos. Su paquete de propiedades se utiliza con compuestos fuertemente no ideales.
- Non-Random-Liquids (NRTL):es una extensión de la ecuación de Wilson. Puede representar el
comportamiento de equilibrios liquido-vapor (VLE), líquido-líquido (LLE) y líquido-líquido-vapor
(VLL). Combina las ventajas de Wilson y Laar.
- UNIversal Functional Activity Coefficient (UNIQUAC): al igual que con NRTL se utilizan para
mezclas de varios componentes en equilibrio líquido-vapor, líquido-líquido o líquido-líquido-vapor.
Este modelo es más detalladado y sofisticados que los demás LACM. Tiene una buena representación
de VLE y LLE, puede obtener las propiedades para un amplio rango de mezclas no electrolitas
utilizando solo dos parámetros. Puede ser aplicada a mezclas que contengan agua, alcohol,
hidrocarburos y varios más.
- Margules y Van Laar: se utiliza para mezclas binarias en líquido-vapor o líquido-líquido.
- Chien Null: método empírico, permite escoger el método termodinámico en función de la sustancia,
buscando el que más se ajuste. Puede llegar a aplicar dos o más métodos.

3.3 Elección del método termodinámico

Como ya se ha mencionado anteriormente en este mismo capítulo la existencia de tantos métodos se
fundamenta en los distintos comportamientos de los fluidos puros y sus mezclas. Los dos comportamientos
principales se basan en la composición de las moléculas y la existencia de las fuerzas intermoleculares. Las
móleculas donde sus fuerzas intermoleculares sean despreciables tendrá muy poca diferencia de
electronegatividad entre sus átomos y se considerarán móleculas no polares (hidrocarburos). En cambio,
cuando las fuerzas no puedan ser despreciadas, la diferencia de electronegatividad entre los átomos será alta y
se considerarán moléculas polares. Dentro de las móleculas polares se dividen en moléculas débilmente
polares (aldehídos, éteres) y moléculas altamente polares (alcoholes, agua).
- Mezclas de sustancias apolares (hidrocarburos): se utilizan EOS tradicionales en todo el rango de
presión (PR o SRK).
- Mezclas de sustancias polares:
o Baja presión: para sustancias en estado líquido se utilizan métodos LACM. En cambio, para
sustancias en estado gaseoso se utiliza la ley de gases ideales.
o Media presión: en función de la precisión que se requiera y de los parámetros que vayan a
necesitarse se puede elegir entre EOS polares o LACM.

7
 8 Métodos termodinámicos

o Alta presión: se recomienda utilizar EOS polares o bien en caso de no contar con los
parámetros requeridos se utilizan EOS predictivas.
Para escoger el método más adecuado para el proceso se pueden utilizar como árbol de decisiones la Tabla 3-1
y la Tabla 3-2. La primera diferencia las sustancias polares de las apolares, igual que como se mencionaba
anteriormente. La segunda divide las sustancias polares no electrolíticas en función de su presión. También
puede consultarse Tabla 3-3, identifica el tipo de proceso o industria con el método más recomendado.

Polar No Polar
No Pseudocomponente
Electrolito Real
Electrolito P>1 atm Vacío

Peng-Robinson (PR)
Electrolito Grayson, Ideal,
Tabla 3-2 Soave-Redlich-Kwong
NRTL BK10 BK10
(SRK)

Tabla 3-1 Métodos para compuestos Polares-No polares

Polar – No electrolito
P>10 bar P<10 bar
Parámetros de interacción
Sí No Sí No

EOS polares EOS predictivas Líquido-Líquido Líquido-Líquido

Sí No Sí No
SR-Polar SRK-HV2
SRK-BM PSRK Wilson, UNIFAC,
NRTL
NRTL, UNIFAC-LL UNIFAC-Lyngby,
UNIQUAC
UNIQUAC UNIFAC Dortmund

Tabla 3-2 Métodos para compuestos Polares/No electrolitos

 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 9

Tabla 3-3 Ecuaciones recomendadas para cada proceso, Introducción a la simulación en Ingienería Química2

3.4 Cálculo de entalpía y entropía

La entalpía de gas ideal utilizada por HYSYS es igual a la entalpía de formación a 25ºC. Puesto que el valor de
la entalpía dependerá del valor de referencia escogido en este caso los valores de entalpía para presiones
pequeñas serán negativos. Para la entropía se considera 25ºC y 1atm.
- Ecuación de Peng Robinson:

H − H ID 1  da   V + (20.5 + 1)b 
= Z − 1 − 1.5  a − T dt  ln  V − (20.5 − 1)b  (3-10)
RT 2 bRT  

S − So ID  T da   V + (2 + 1)b 
0.5
P A
= ln( Z − B ) − ln o − 1.5  a dt   V − (20.5 − 1)b 
ln (3-11)
R P 2 bRT  

∑∑ x x (a a ) (1 − k )
n n
=a i j i j
0.5
ij (3-12)
=i 1 =j 1

- Ecuación de SRK:

H − H ID 1  da   b 
= Z −1−  a − T dt  ln 1 + V  (3-13)
RT bRT

S − So ID P A  T da   B 
= ln( Z − b) − ln o −  ln 1 + 
B  a dt   V 
(3-14)
R P

∑∑ x x (a a ) (1 − k )
n n
(3-15)
=a i j i j
0.5
ij
=i 1 =j 1

9
10 Métodos termodinámicos
 4 HYSYS

omo ya se ha comentado programa de simulación HYSYS es de gran utilidad en procesos de ingeniería.

C A través de una interfaz sencilla permite la simulación de diferentes procesos. Una de sus grandes
cualidades es la base de datos con la que cuenta, ofrece la posibilidad de utilizar diferentes métodos
termodinámicos facilitando el ajuste al tipo de componente que vaya a utilizarse. Otra gran ventaja que ofrece
el programa es la clara visualización de los equipos que forman el proceso.

4.1 Propiedades
4.1.1 Componentes
La primera vez que se vaya a utilizar el programa debe seleccionarse New para la elaboración de un nuevo
caso. En ocasiones poseteriores podrá abrirse el archivo directamente desde la ventana emergente donde
aparecen los últimos trabajos que han sido abiertos. El programa indicará que información requiere:
componentes, paquetes de propiedades termodinámicas. Primero se añaden los componentes que formarán el
sistema tal como aparece en la Figura 4-1. Al seleccionar Add aparece el listado de componentes que posee en
su base de datos. Pueden filtrarse por familias, pureza o nombre, Figura 4-2. HYSYS permite trabajar con
listados de componentes independientes mediante el uso de distintas listas, a estas listas se les podrá asociar
métodos termodinámicos distintos para el cálculo de sus propiedades.

Figura 4-1. Components Lists

11
12 Hysys

Figura 4-2. Componenetes

4.1.2 Paquete termodinámico

Seleccionados los componentes se escoge el paquete de propiedades termodinámicas. Es muy importante
escoger de forma correcta el método termodinámico a partir del cual se calcularán las propiedades. Si no se
escogiese el método correcto toda la simulación podría quedar invalidada al obtener resultados erróneos. El
programa Aspen HYSYS posee una gran variedad de métodos termodinámicos lo que permite escoger el que
más se ajuste a las necesidades de la simulación. En la Figura 4-3 se encuentran todos los métodos que posee
el programa. Además, pueden importarse métodos propios del usuario en fomato fpk.

Figura 4-3. Métodos Termodinámicos

 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 13

Al seleccionar el método termodinámico deseado este debe asociarse a uno de los listados de componentes,
por defecto el programa lo asocia al primer listado. Para cada método se podrá seleccionar diferentes opciones
para ajustarlo así lo máximo posible al uso que vaya a utilizarse. Puede verse un ejemplo para el método Peng-
Robison en la Figura 4-4.

Figura 4-4. Opciones del método Peng Robinson

4.2 Simulación
4.2.1 Caso principal
Una vez escogidas las propiedades se procede a añadir los elementos que posee el sistema. Aparece una
ventana emergente Palette, Figura 4-5, dónde pueden seleccionarse los elementos que se necesiten. Posee
distintas pestañas donde clasifica los elementos para que el usuario pueda encontrarlos más fácilmente. Para
utilizar cualquiera de los elementos en la simulación basta con arrastrarlos a la pantalla. También deben
añadirse los flujos de materia, flecha azul y energía, flecha roja.

Figura 4-5 Palette

14 Hysys

Al añadir los equipos estos aparecerán en rojo, Figura 4-6, indicando que requieren la conexión de los flujos de
materia y energía. Para conectarlos basta con hacer doble click en el elemento que se quiera conectar e
indicarle cuales son sus flujos de entrada y salida, puede verse en el intercambiador de la Figura 4-7. En el
apartado de Parameters de la pestaña de Design pueden modificarse algunas propiedades del equipo, como la
pérdida de carga sufrida en su interior o el tipo de equipo, para intercambiadores puede escogerse por ejemplo
Simple End Point o Simple Weighted (intercambiador a contracorriente puro). El programa indica en todo
momento que datos necesita, esto supone una gran ayuda tanto para el nuevo usuario como para uno ya
experimentado. Al definir los flujos de entrada y salida se le puede asociar la lista de comoponente “Basis” a la
que debe hacer referencia. La lista incluirá tanto los componetes como el método termodinámico asociado.

Figura 4-6 Elementos

Figura 4-7 Intercambiador

Conectados todos las entradas y salidas de los sistemas estos cambiarán de color, celeste para los flujos de
materia y grisáceo con líneas amarillas para los equipos, Figura 4-8. Estos colores permanecerán hasta que el
usuario defina sus propiedades o bien el programa sea capaz de calcularlas. No deben cumplimentarse todas
las variables existentes en el sistema pues se produciría una sobreespecificación y el programa no tendría
margen para el cálculo. En la Figura 4-9 aparecen algunas de los datos que el usuario podría introducir.
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 15

Figura 4-8 Esquema Simple

Figura 4-9. Incógnitas del sistema

Los datos principales se introducirán en la pestañan Conditions de los flujos de materia, Figura 4-10, los
principales datos a introducir suelen ser: presión, temperatura, caudal másico y fracción de vapor. No siempre
se introducen todos, pero al menos uno de ellos suele ser necesario. En la pestaña Composition se introducirá
la composición del flujo, puede escogerse la fracción en base de la masa, moles… Los componentes
disponibles dependerán de la lista de componentes a la que se asocie el flujo. Si el Solver está activo irá
resolviendo el sistema a medida que se le introduzcan los datos.
16 Hysys

Figura 4-10 Flujos de materia, Worksheet

Definidas o calculadas todas las condiciones del sistema los colores se vuelven mas oscuros, Figura 4-11.
Pueden sacarse tablas donde vengan resumidas las propiedades principales de cada flujo o bien la que se
deseen evaluar. Para obtener las tablas generadas automáticamente por el programa basta con pulsar Model
Summary, esta opción se encuentra en el apartado Summaries de la pestaña Home. En la Figura 4-12 se
observa la tabla referida a los flujos de materia, también pueden visualizarse tablas referidas a los equipos del
sistema como el intercambiador o el compresor. El programa permite guardarlas, copiarlas, resetearlas o bien
enviarlas a Excel. Desde la pestaña de Flowsheet Case también pueden obtenerse tablas junto con el esquema
del circuito, con el botón derecho se despliega una lista de opciones, la última de estas es Add Workbook
Table. Se pueden implementar tres tablas: materia, energía y composición.

Figura 4-11 Esquema simple 2

 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 17

Figura 4-12 Model Summary Grid

4.2.2 Casos de estudio

Identificado el caso principal el programa permite estudiar como varian una serie de variables a partir de la
variación de otras variables independientes. Se pueden crear tantos casos de estudio como se deseen.
Desde la pestaña de Home en el apartado de Analysis pueden añadirse los casos, estos casos se simularán de
manera independiente al caso principal. Los resultados podrán visualizarse sin realizar cambios en el caso
principal. Para realizar el estudio el primer paso es incluir las variables que se desean analizar, Figura 4-13.
Las variables independientes serán las variables de entrada que se han definido en el caso principal. Las
dependientes serán cualquier otra que no haya sido definida, HYSYS permite analizar una gran cantidad de
propiedades tanto de los componentes como de los equipos.
Una vez seleccionadas, se acotan los casos que desean estudiarse, Figura 4-14, en función del tipo de estudio
los estudios se carcaterizarán de una manera u otra:
- Sensibilidad: se debe definir el valor mínimo y máximo de las variables, así como el Step entre uno y
otro o el número de casos a analizar. Al introducir varias variables independientes se analizan los
casos para la modificación de una de ellas cada vez sin mezclarlas. Para dos variables con 10 casos
cada una tendríamos un total de 20.
- Anidado: al igual que con el anterior se debe definir el valor mínimo y máximo de las variables, así
como el Step entre uno y otro o el número de casos a analizar. Este estudio muestra la combinación de
cada una de las variables independietes introducidas para sus distintos valores. Para dos variables
independientes con 10 casos cada unas tendríamos un total de 100 combinaciones.
- Discreto: permite definir los valores de las variables para cada caso que desee estudiarse. Para 5 casos
con 2 variables independientes habrá que definir los valores de cada una para cada caso.
- Base & Shift: permite definir unos casos principales (Base) junto con otros donde se modifica una de
sus variables. Para incluir las modificaciones se debe definir el incremento (Shift) de las variables.
Cada una de las las bases sería independiente de las demás. Si se define una base con dos variables los
resultados mostrarán 3 casos. El primero serían los valores de las variables dependientes para la base
definida. El segundo caso modificaría la primera variable sumando la definida en la base con el Shift,
la segunda variable mantendría el valor de la base constante. El tercer caso por el contrario mantendría
constante la primera y modificaría la segunda. Si se considerasen mas bases cada una tendría sus
resultados de manera independiente, también sería necesario definir la Base y el Shift.
Definido el studio se pulsa Run, en la pestaña de resultados aparecerá la tabla con los resultados obtenidos,
Figura 4-15. Si se introducen varias variables independientes se analizan todos los casos posibles, por ejemplo,
si las variables son presión y temperatura los casos serían la asociación de cada temperatura con cada presión.
18 Hysys

Figura 4-13. Selección de variables.

Figura 4-14 Cotas de las variables independientes.

Figura 4-15 Resultados del caso estudio.

 5 CASOS SIMULADOS

al como se explicaba en el capitulo 4 de este documento, para resolver el sistema con HYSYS basta con

T introducir los componenetes que lo forman y los datos conocidos. El programa cálcula todas las
variables restantes. En el cálculo manual la resolución se ciñe estricatamente al enunciado sin embargo
HYSYS muestra muchos más. Esto permite evaluar parámetros que en principio no se hubiesen considerado.
En este capitulo se podrá comprobar la resolución mediante ambos métodos y su comparación.
Una vez realizada la simulación principal se pueden realizar modificaciones que nos permitan evaluar el
comportamiento de la planta a las variaciones que introduzcamos. Antes de realizar las variaciones es preciso
conocer que parámetros son dependientes y de que variables dependen. Si mantenemos la simulación activa,
pero eliminamos alguna de las variables que habíamos impuesto podremos conocer de manera rápida que
líneas se ven afectadas esto se realizará para cada uno de los casos.

5.1 Caso 1

5.1.1 Enunciado

Una instalación para enfriamiento de agua desarrolla un ciclo simple de compresión mecánica con amoniaco,
utilizando un condensador de aire y una válvula de expansión termostática. Los siguientes datos fueron
medidos:
• Presión manométrica de condensación: 1454 KPa
• Presión manométrica de evaporación: 328.3 KPa
• Temperatura del gas de aspiración: 5ºC
• Temperatura del amoniaco a la salida del condensador: 35ºC
• UA del evaporador: 4032 W/K
• Caudal de agua sobre el evaporador: 3300 l/h
• Temperatura del agua a la entrada al evaporador: 20ºC
• Rendimiento volumétrico del compresor:
• Rendimiento isentrópico del compresor:

En las expresiones de los rendimientos pc y pe son las expresiones absolutas de condensación y evaporación
respectivamente. En estas condiciones calcular:
A. Temperatura del agua a la salida del evaporador. [ºC]
B. Caudal de refrigerante [Kg/s] y desplazamiento volumétrico del compresor [m3/h].
C. Potencia frigorífica [kW], potencia de compresión [kW] y COP de la instalación.

Para mejorar el comportamiento de la instalación se desea cambiar el condensador de aire por un condensador
de agua que consiga una nueva temperatura de condensación de 30ºC, suponiendo que la temperatura de
evaporación permanece constante. Se utiliza como fluido de condensación 8000l/h de agua a 20ºC. Calcular en
esta nueva situación:
D. Nueva potencia frigorífica [kW], potencia de compresión [kW] y COP de la instalación.
E. Efectividad [%] y UA [W/K] que debe de tener el nuevo condensador de agua.

19
20 Casos Simulados

Nota: Suponer que no existe pérdida de presión por rozamiento. Y que para el caso de la condensación por el
agua el subenfriamiento del líquido es nulo.

5.1.2 Solución manual

3 2
3 2

4 1
4 1

Figura 5-1 Ciclo simple, condensado por aire (izq) y condensado por agua (der).

5.1.2.1 Planteamiento

Para obtener la temperatura de agua a la salida del evaporador se utiliza la eficiencia del intercambiador y el
método E-NTU. Este método permite el cálculo de la eficiencia a partir del AU del intercambiador y su cmín
(5-1), En este caso el cmin será el perteneciente al agua. Una vez conocida la eficiencia se puede calcular el
calor transmitido entre los fluidos. Para ello se utiliza la ecuación que define la eficiencia, calor cedido entre el
máximo que podría haberse transmitido. En esta ecuación para simplificarla se puede despejar el mcp del agua
con el cmin y calcular directamente la temperatura de salida del agua (5-2).

− AU (5-1)
m f c pw
ε = 1− e

Q f m wc p (Twe − Tws ) (5-2)

=ε = w

Q m wc p (Twe − TEv )
MÁX w

Hacer el balance de energía en el evaporador, (5-3) y (5-4), permitiría conocer el caudal de refrigerante, no
obstante, primero debe hallarse la entalpía en los flujos de entrada y salida del refrigerante al evaporador. Para
la línea de entrada al compresor el cálculo es sencillo pues se conoce su temperatura y presión, basta con
utilizar estos datos para entrar en la tabla (Anexo D). Para la línea de entrada al evaaporador en principio, no
puede mirarse de forma directa en la tabla pues al encontrarase en la zona bifásica sería necesario conocer la
fracción de vapor. Sin embargo, si se conoce la temperatura y presión a la salida del condensador, el paso del
refrigerante a través de la válvula no produce cambios en la entalpía, por tanto, tendrán la misma entalpía.
- Agua:

Q f m wc p w (TWE − =
= TWS ) vw ρ wc p w (TWE − TWS ) (5-3)
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 21

- Refrigerante:

Q f m R (h4 − h3 )
= (5-4)

Obtenido el caudal de refrigerante que circula por el circuiro puede calcularse el desplazamiento volumétrico
del compresor. Para ello se necesita el rendimiento volumétrico (5-5), obtenido a partir de las presiones a la
entrada y salida del compresor. Utilizando el rendimiento y el caudal se puede determinar el desplazamiento
volumétrico.

PCD (5-5)
η=
V
0.8 − 0.04
PEV

m R v1
v t = (5-6)
ηV

La potencia del compresor puede determinarse mediante el balance de energía en el mismo. Para calcular la
entalpía a la salida del compresor se utiliza en rendimiento isentrópico, (5-7) y (5-8). Para hallar la COP
bastará con dividir la potencia del evaporador entre la potencia del compresor (5-9).
- Compresor:

m R (h2 s − h1 ) (5-7)
W= mR (h2 − h1=
)
C
ηS

PCD (5-8)
η S 0.75 − 0.04
=
PEV

Q f (5-9)
COP =
W C

El cambio del condensdor de aire por el condesandor de agua requiere de la repetición de los cálculos. El
desplazamiento volumétrico se mantiene constante, pero el rendimiento cambia al hacerlo la presión de
condensación. Mediante el uso de (5-5) y (5-6) se determina el nuevo caudal de refrigerante- La nueva
potencia frigorífica, potencia del compresor y COP se hallaría de la misma forma detallada anteriormente, en
este caso se supone un subenfriamiento nulo. Para hallar la efectividad y UA del nuevo condensador bastaría
con hacer balance de energía y aplicar las ecuaciones.

5.1.2.2 Resultados

- Condensador de aire:
Gracias a las ecuaciones descritas anteriormente se obtienen todoss los parámetros que caracterizan el sistema,
Tabla 5-1. Utilizando estos datos junto con el resto de las ecuaciones se obtiene el desplazaiento volmétrico, la
potencia del evaporador, la potencia del compresor y el COP, Tabla 5-2.
22 Casos Simulados

P T R
m h s
(kPa) (ºC) (kg/s) (kJ/kg) (kJ/kgK)
1 429.6 5 0.045 1475.3 5.668
2s 1555.3 99.73 0.045 1664.56 5.668
2 1555.3 0.045 1788.05
3 1555.3 35 0.045 366.25
4 429.6 0 0.045 366.25
5 101.3 20 0,917
6 101.3 7 0,917
Tabla 5-1 Propiedades de los componentes del sistema, caso 1.

Vt (m3/h)  (kW)
Q 𝑾𝑾̇ 𝑪𝑪 (kW) COP
f

73.2 49.87 14.07 3.55

Tabla 5-2 Desplazamiento volumétrico, potencia del evaporador, del compresor y COP, caso 1.
- Condensador de agua:
Al cambiar el tipo de condensador y disminuir la temperatura de condensación también lo hace su
presión. Al disminuir la presión la entalpía a la entrada y salida del condensador también disminuye,
Tabla 5-3. Al disminuir la presión la potencia del compresor necesaria será menor, puesto que la
entalpía a la salida del condensador es menor la potencia del evaporador aumenta. El aumento de la
potencia del evaporador y la disminuición de la potencia del evaporador aumentará el COP, Tabla 5-4.
El nuevo condensador tiene una mayor eficiencia y por tanto mayor AU, Tabla 5-5.

P T m R h s
(kg/s) (kJ/kg (kJ/kgK
(kPa) (ºC)
) )
1 429.6 5 0,0475 1475.3 5.668
2s 1167,20 99.73 0,0475 1617,7 5.668
2 1167,20 0,0475 1697,4
3 1167,20 35 0,0475 341,8
4 429.6 0 0,0475 341,8
5 101.3 20 0,917
6 101.3 7 0,917
Tabla 5-3 Propiedades de los componentes del sistema, caso 1

Q f (KW) 𝑾𝑾̇ 𝑪𝑪 (KW) COP

53,80 10,54 5,01

Tabla 5-4 Potencia del evaporador, del compresor y COP, caso 1.

εc (%) UAC (W/K)

69,2 10965

Tabla 5-5 Eficiencia y UA del condensador

 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 23

El circuito con el condensador de aire se representa en la Figura 5-2. Mediante la representación se puede
comprobar mejor el comportamiento mencionado anteriormente.

3
2

4 1

Figura 5-2 Diagrama p-h, caso 1

5.1.3 Solución HYSYS

Figura 5-3 Ciclo simpe (HYSYS), condensado por aire (izq), condensador por agua (der).

5.1.3.1 Planteamiento

Para la resolución de este caso mediante HYSYS solo se han utilizado dos componentes el amoniaco y el
agua, a pesar de tener un aerocondensador no es necesario definir el aire en la lista. El método termodinámico
escogido ha sido Wilson que como podrá comprobarse da muy buenos resultados tanto para el amoniaco como
para el agua.
El evaporador utilizado ha sido un intercambiador a contracorriente ideal, se ha escogido este debido a que los
cálculos teóricos realizados también se basan en ecuaciones simplificadas. La línea de entrada del agua al
evaporador solo requiere de sí misma para definirse por completo, Figura 5-4.
24 Casos Simulados

Figura 5-4 Datos de entrada al sistema (HYSYS), caso 1.

5.1.3.2 Resultados

Una vez el Solver ha calculado todos los datos que faltaban podemos conocer el estado del refrigerante en cada
punto del sistema.
- Condensador de aire: Se puede comprobar a simple vista que los resultados obtenidos por HYSYS,
Tabla 5-6(5-9), son muy similares a los obtenidos anteriormente de forma manual, Tabla 5-3. En este
caso la temperatura escogida por el sistema no es cero pero si muy cercana.
A partir de las propiedades y con las ecuaciones descritas anteriormente se obtienen las potencias del
evaporador, compresor y el COP, Tabla 5-7.
Nº Presión Temperatura Caudal Vapor / Phase
[kPa] [C] [kg/s] Fraction
1 429,6 5 0,0439 1
2 1555,3 158,04 0,0439 1
3 1555,3 35 0,0439 0
4 429,6 0,32 0,0439 0,132
5 101,3 20 0,9167 0
6 101,3 7,43 0,9167 0
Tabla 5-6 Propiedades de los componentes del sistema (HYSYS), caso 1.

Q f (KW) 𝑾𝑾̇ 𝑪𝑪 (KW) COP

50,01 15 3,33

Tabla 5-7 Potencia del evaporador, del compresor y COP, caso 1.

- Condensador de agua: al utilizar cambiar el condensador a pesar de no haber subenfriamiento a la

salida, el título de vapor a la entrada del evaporador disminuye, Tabla 5-8. La potencia del evaporador
se mantiene prácticamente constante, sin embargo, la potencia del compresor si disminuye
notablemente, Tabla 5-9.
La nueva eficiencia y UA del condensador calculada por HYSYS se refleja en Tabla 5-10
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 25

Nº Presión Temperatura Caudal Vapor / Phase

[kPa] [C] [kg/s] Fraction
1 429,6 5 0,0429 1
2 1157,7 121,07 0,0429 1
3 1157,7 30 0,0429 0
4 429,6 0,32 0,0429 0,112
5 101,3 20 0,9167 0
6 101,3 7,43 0,9167 0
7 101,3 20 2,22 0
8 101,3 26,26 2,22 0
Tabla 5-8 Propiedades de los componentes del sistema, condensador de agua (HYSYS), caso 1.

Q f (KW) 𝑾𝑾̇ 𝑪𝑪 (KW) COP

49.73 10.31 3,54

Tabla 5-9 Potencia del evaporadot, del compresor y COP, evaporador de agua (HYSYS), caso 1

εc (%) UAC (W/K)

62,62 7562

Tabla 5-10 Eficiencia y UA del condensador, condensador de agua, caso 1.

5.1.4 Comparación de los resultados

Para comprobar la fiabilidad del programa se comparan los resultados pedidos iniciamente en el enunciado y
obtenidos manualmente y mediante HYSYS.
Utilizando el condensador de aire el error cometido en el calculo de la temperatura de salida del agua es el
mayor. En cambio, para el cálculo de la potencaia del evaporador potencia del compresor y el COP, el error es
muy pequeño, apenas supone diferencia. A pesar de haber considerado distintos rendimientos en cada uno de
los métodos. En el primero se considera el rendimiento iséntropico del compresor. HYSYS no permite
considerar el rendimiento iséntropico por tanto el que se considera es el politrópico, este se ajusta más a a la
realidad, pero utilizarlo en el modo manual complicaría los cálculos.
Al cambiar el condensador por uno de agua el error aumenta, especialmente en el UA del condensador.
26 Casos Simulados

Manual Hysys %ERROR

TWS (ºC) 7 7,43 6,2%

mR (kg/s) 0.0449 0,0438 2,3%

QF (kW) 49,87 49,73 0,3%

WC (kW) 14,05 14,05 0,0%

COP 3,55 3,54 0,3%

QF’ 53,80 10,31 2,2%

WC’ 10.54 10,31 2,2%

COP’ 5.1 4,83 5,4%

ε’ 69.2 62,62 9,5%

UAC’ 10965 7562 31,0%

Tabla 5-11 Comparación de resultados, caso 1

5.1.5 Análisis de sensibilidad

Tal como hemos visto en este caso solo hay un flujo de materia de refrigerante. Al principio del capitulo se
menciono que para comprobar la dependencia de unas líneas con otras podían eliminarse variables y ver donde
afectaba. Para este caso se han realizado 2 estudios distintos.

5.1.4.1 Presión de condensación

Al modificar la presión a la salida del compresor todo el sistema cambiará.

Figura 5-5 Presión de condensación.

Al variar la presión a la salida del compresor la temperatura y la potencia varía, una mayor presión supondría
mayor temperatura y un mayor consumo de potencia en el compresor. En cambio, el UA del condensador
disminuiría con el aumento de la presión.
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 27

170 15
Temperatura (ºC)
165
14,5

Potencia (kW)
160
14
155
13,5
150

145 13
1400 1450 1500 1550 1600 1650 1700 1400 1450 1500 1550 1600 1650 1700
Presión (kPa) Presión (kPa)

Figura 5-6 Evolución de la temperatura de salida Figura 5-7 Evolución de la potencia del compresor
del compresor frente a su presión. frente su presión de salida.

1420

1400

1380
UA (W/C)

1360

1340

1320

1300
145 150 155 160 165 170
Presión (kPa)

Figura 5-8 UA del condensador frente su presión.

5.1.4.2 Presión de evaporación.

Igual que ocurría con la la presión a la salida del compresor, modificar la presión de entrada también alteraría
todo el sitema.

Figura 5-9 Presión de evaporación.

Aumentar la presión provocaría la elevación de la temperatura de evaporación del refrigerante y por tanto, se
produciría un menor enfriamiento en el agua. Por otra parte, la potencia necesaria del compresor sería menor al
existir un menor salto de presiones entre la entrada y la salida.
28 Casos Simulados

9 22
20
Temperatura (ºC)

Potencia (kW)
18
7
16
6
14
5
12
4 10
350 400 450 500 360 410 460
Presión (kPa) Presión (kPa)

Figura 5-10 Temperatura de salida del agua frente Figura 5-11 Potencia del compresor frente a la presión
a la presión del evaporador. del evaporador.

5.2 Caso 2
5.2.1 Enunciado
La instalación de refrigeración por amoniaco de la figura consta de dos evaporadores, uno (A) alimentado por
una bomba, que produce una potencia frigorífica de 75KW y otro (B) en expansión directa, con una válvula de
expansión termostática que produce un recalentamiento del vapor a la salida del evaporador 5ºC, y que
desarrolla 125kW de potencia frigorífica. La presión a la entrada al evaporador A marca 4.28 bar, la
temperatura de evaporación del evaporador B es de -20ºC y la de condensación 40ºC. Se pide:
A) Calcular el caudal que debe mover la bomba para que el título de vapor del punto 6 sea el 75%.
B) Si el rendimiento volumétrico e isentrópico de los dos compresores puede estimarse en 0.7 y 0.8
respectivamente, calcular el desplazamiento volumétrico de los compresores y el COP de la
instalación comparándolo con el de Carnot.
C) Dibujar sobre el diagrama p-h adjunto todos los puntos de la instalación de la figura, indicando su
temperatura y entalpía.

NOTA: Suponer que las pérdidas de carga por rozamiento en todos los elementos, el consumo de energía de la
bomba y el incremento de presión en la bomba del líquido son despreciables. El subenfriamiento del líquido de
condensador es de 4ºC.
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 29

5.2.2 Solución manual

Compresor
de Alta

3
1

7
6
5
10

Compresor Evaporador A
de Baja

9 8

Figura 5-12 Esquema 2

5.2.2.1 Planteamiento

El caudal del refrigerante que circula por el Evaporador A se obtiene mediante el balance de energía en
él,(5-10). Se conoce la potencia del evaporador, el título de vapor de la salida (permitirá el cálculo de su
entalpía) y la entalpía a la entrada, al porceder de un tanque tendrá la entalpía de saturación para 4.29 bar. A
partir del título de vapor y las entalpías de saturación de gas y líquido se saca el valor de la entalpía a la salida,
(5-11), pudiendo así despejar el caudal.

Q fA m RA (h6 − h5 )
= (5-10)

h6 = 0.75hg + 0.25hl (5-11)

El cálculo del desplazamiento volumétrico de los compresores, así como el COP del sistema requiere de un
previo cálculo de los caudales. Utilizando de nuevo el balance de energía se puede obtener el caudal del
evaporador B, (5-12). La entalpía de entrada, h8, será la misma que la de salida del botellón, h7, entalpía de
saturación líquido. El Evaporador B está alimentado por una válvula de expansión termoestática, esta se
encarga de que se produzca un recalentamiento de 5 grados a la salida del mismo. Conociendo su temperatura
y presión puede acceder a las tablas (Anexo D), y estimar la entalpía a la salida del evaporador, h9. Despejados
los caudales de ambos evaporadores se puede hallar el caudal del condensador mediante el balance de energía
entre los flujos de entrada y salidad del botellón, (5-13).
30 Casos Simulados

Q fB m RB (h9 − h8 )
= (5-12)

m RB h7 + m RAh5 + m RC h1 = m RC h4 + m RB h10 + m RAh6 (5-13)

Para determinar el desplazamiento volumétrico de los condensadores basta con despejar de la ecuación, (5-14),
que representa el rendimiento volumétrico. La potencia aportada al compresor se puede obtener analizando el
salto de entalpía que se produce en el refrigerante. En un compresor ideal la entropía se mantendría constante a
la entrada y a la salida del compresor, sin embargo, en este caso se debe tener en cuenta el rendimiento
isentrópico que definirá la potencia que realmente debe aportar el compresor, (5-15). El calculo de las
potencias de ambos compresores permite conocer el COP del sistema, (5-16).

m R v (5-14)
ηV =
Vt

W h −h (5-15)
η s = CBS = 10 s 9
WCB h10 − h9

Q fA + Q fB (5-16)
COP = 
WCA + WCB

5.2.2.2 Resultados

Al igual que en el caso anterior mediante el uso de las ecuaciones y tablas se obtienen las propiedades del
refrigerante en todo el sistema, Tabla 5-12. Mediante el uso de las propiedades y las ecuaciones se obtiene el
desplazamiento volumétrico de los compresores Tabla 5-13 y las potencias de los compresores y el COP,
Tabla 5-14.
P T Caudal h s
(kPa) (ºc) (kg/s) (kJ/kg) (kJ/kgK)
1 429 0 0,1958 1462,2 5,621
2s 1555,3 93,3 1647,4 5,621
2 1555,3 0,1958 1693,8
3 1555,3 36 0,1958 371,12
4 429 0 0,1958 371,12
5 429 0 0,0792 200
6 429 0 0,0792 1146,7
7 429 0 0,1000 200
8 190,1 -20 0,1000 200
9 190,1 -15 0,1000 1449,6 5,95
10s 429 38,6 1558,28571 5,95
10 429 0,1000 1585,45714
Tabla 5-12 Propiedades de los componentes del sistema, caso 2.
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 31

VTA (m3/h) VTB (m3/h)

328 291

Tabla 5-13 Desplazamiento volumétrico de los compresores, caso 2.

WCA (KW) WCB (KW) COP

45.4 13.6 3.39

Tabla 5-14 Potencia de los compresores y COP, caso 2.

3
2

5,7 4 6 1 10

8 9

Figura 5-13 Diagrama P-h, caso 2

5.2.3 Solución HYSYS

5.2.3.1 Plantemiento

Para la resolución de este caso mediante HYSYS solo se ha necesitado utilizar un componente, amoniaco. Para
el amoniaco el método termodinámico escogido ha sido Peng-Robinson que como podrá compararse en los
resultados da muy buenos resultados.
En la Figura 5-14 se puede observar el esquema del sistema una vez ya se han definido y calculado todas las
incógnitas. Para la simulación de los evaporadores se han escogido dos Heater y para el condensador un
Cooler. Cada uno de ellos lleva asociado un flujo de energía al igual que los compresores. Este sistema
también lleva incluido un bloque Tank donde el fluido se encuentra en equilibrio. El programa permite
multiples entradas en el tanque, pero una única salida de gas o líquido. Para solventar este inconveniente basta
con introducir un bloque Tee que divida el flujo en dos manteniendo las propiedades del líquido saturado, los
puntos 5,7 y 11 tienen las mismas propiedades a diferencia únicamente de su caudal.
32 Casos Simulados

Figura 5-14 Esquema caso 2 (HYSYS).

Colocados todos los elementos y unidos los flujos de materia y energía se introducen los datos de entrada,
Figura 5-15 y Figura 5-16. Además de estos datos es necesario indicar la composición de los flujos de materia,
para este caso en todos los flujos tendremos amoniaco. Una vez introducido todos los datos basta con activar el
Solver.

Figura 5-15 Datos de entrada en de los flujos de materia, caso 2

 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 33

Figura 5-16 Datos de entrada de los flujos de energía, caso 2

5.2.3.2 Resultados

De nuevo una vez el Solver ha calculado todos los datos que faltaban podemos conocer el estado del
refrigerante en cada punto del sistema. Se puede comprobar a simple vista que los resultados obtenidos por
HYSYS, Tabla 5-15, son muy similares a los obtenidos anteriormente de forma manual, Tabla 5-12.
N Pressure Temperature Mass flow Vapor / Phase
[kPa] [C] [kg/s] Fraction
1 429,3 0,30 0,1915 1
2 1555,3 123,46 0,1915 1
3 1555,3 36 0,1915 0
4 429,3 0,30 0,1915 0,14
5 429,3 0,30 0,0772 0
6 429,3 0,30 0,0772 0,75
7 429,3 0,30 0,0982 0
8 187,38 -20 0,0982 0,07
9 187,38 -15 0,0982 1
10 429,3 56,26 0,0982 1
11 429,3 0,30 0,1754 0
Tabla 5-15 Propiedades de los flujos de materia de HYSYS, caso 2

VTA (m3/h) VTB (m3/h)

332 292

Tabla 5-16 Desplazamiento volumétrico de los compresores (HYSYS), caso 2.

WCA (KW) WCB (KW) COP

47,27 14,17 3,26

Tabla 5-17 Potencia de los compresores y COP, caso 2

5.2.4 Comparación de los resultados

De nuevo para comprobar la fiabilidad del programa se comparan los resultados pedidos iniciamente en el
enunciado y obtenidos manualmente y mediante HYSYS. Para este caso el error competido en el calculo del
caudal es pequeño, apenas supone diferencia. En cambio, para el cálculo del desplazamiento volumétrico y de
la potencia del compresor y el COP si hay diferencias bastante considerables. Las primeras pueden deberse al
método termodinámico o la base de datos para el volumen específico del refrigerante. Las diferencias en la
potencia se deben de nuevo rendimiento considerado en cada uno de los métodos. El COP de Carnot al ser un
34 Casos Simulados

cálculo meramente teórico se mantiene constante tanto para uno como el otro.

Manual HYSYS %ERROR

mRA (kg/s) 0,0792 0,0772 2,5%

VTA (m3/h) 291 292 0,3%

VTB (m3/h) 328 332 1,2%

WCB (kW) 13.59 14,17 4,2%

WCA (kW) 45.37 47,27 4,2%

COP 3.39 3,26 4,0%

COPCARNOT 4.66 4,66 0

Tabla 5-18 Resultados Manual y HYSYS, caso 2

5.2.5 Análisis de sensibilidad

Tal como hemos visto en este caso se tienen tres flujos de materia que se unen y divididen en el tanque. Al
principio del capitulo se mencionó que para comprobar la dependencia de unas líneas con otras podían
eliminarse variables y ver donde afectaba. Para este caso se han realizado 4 estudios distintos.

5.2.5.1 Línea de alta presión

Las líneas de baja y media presión son independientes de la línea de alta presión, por tanto, como puede
apreciarse en la Figura 5-17 al eliminar uno de los datos de esta línea las demás permanecen intactas y no se
ven afectadas. Para ella se han realizado dos de los estudios mencionados.

Figura 5-17 Línea de alta presión.

 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 35

Para el primer estudio se ha modificado la presión a la salida del compresor de alta presión y se ha
representado la variación de la potencia del compresor y del condensador. En ambos casos al aumentar la
presión han aumentado las potencias. Estos resultados son coherentes puesto que si queremos elevar la presión
a la salida del condensador el compresor deberá realizar un esfuerzo mayor o en palabras de energía consumir
más energía y requerir más potencia. En el caso del condensador también se comporta como era de esperar, al
aumentar la presión aumenta la temperatura de condensación y requerirá más energía para lograr el
subenfriamiento deseado.

263 49
262,5 48,5
262 48
Potencia (kW)

Potencia (kW)
261,5 47,5
261 47
260,5 46,5
260 46
259,5 45,5
1475 1525 1575 1625 1475 1525 1575 1625
Presión (kPa)
Presión (kPa)

Figura 5-18 Potencia del condensador frente a su Figura 5-19 Potencia del compresor de alta frente a
presión. la presión de condensación.

En el segundo caso estudiado se ha variado la temperatura del refrigerante a la salida del condensador, se han
analizado 4 parámetros distintos: el título del vapor a la salida de la válvula de expansión, VLV-100; el caudal
del refrigerante del condensador, la potencia del intercambiador de alta presión y la potencia del condensador.
Al variar la temperatura a la salida del condensador la entalpía en este punto cambia ya que depende
directamente de la temperatura del fluido. Una disminuición de la temperatura disminuiría la entalpía, del
mismo modo el aumento de entalpía vendrá asociado con el aumento de la temperatura siempre que no se
halle en la zona bifáscia, para nuestro caso T<40C. Puesto que la válvula se ha considerado ideal la entalpía
antes y después de la válvula será la misma, es decir la variación de la temperatura a la salida del condensador
cambirá la entalpía a la salida de la válvula. Al expandir el refrigerante parte de este se evapora entrando en la
zona bifásica, el cambio de entalpía provocará un mayor o menor evaporación, o en otras palabras un mayor o
menor título del vapor.
Tal como se veía en Figura 5-17 las zonas de baja y media presión no se ven afectada, al no existir cambios la
energía intercambiada en el tanque entre las líneas de alta y media presión debe permanecer constante. La
energía intercambiada puede traducirse como el caudal multiplicado a la diferencia de entalpía que sufre la
corriente. Puesto que al variar la temperatura a la salida del condensador se ha variado la entalpía que posee el
refrigerante al entrar en el tanque desde el condensador, el caudal de refrigerante de esta línea deberá aumentar
o disminuir de tal modo que compense el aumento o disminuición del Δh.
En el caso de producirse una disminuición de la temperatura disminuiría la entalpía y aumentaría el salto
entálpico disminuyendo también así el caudal. Al disminuir el caudal el compresor tendría que realizar un
esfuerzo menor, requeriría de menor potencia. La disminuición de caudal también afectaría del condensador
que a pesar de tener que llevar el refrigerante a una temperatura inferior consumiriía menos potencia. Puede
concluirse por tanto que la potencia del condensador es más dependiente del caudal del refrigerante que del
salto entálpico que se produce en él.
36 Casos Simulados

0,16 0,195
0,194
0,15 0,193
Fracción de vapor

0,14 0,192

Caudal (kg/s)
0,191
0,13 0,19
0,189
0,12 0,188
0,11 0,187
0,186
0,1 0,185
29 31 33 35 37 39 41 29 31 33 35 37 39 41
Temperatura (ºC) Temperatura (ºC)

Figura 5-20 Evolución del título de vapor 4 frente a la Figura 5-21 Caudal del condensador frente su
temperatura de salida del condensador. temperatura de salida.

48,5 262,5

48 262
Potencia (kW)

Potencia (kW)

47,5
261,5
47
261
46,5

46 260,5

45,5 260
29 31 33 35 37 39 41 29 31 33 35 37 39 41
Temperatura (ºC) Temperatura (ºC)

Figura 5-22 Potencia del compresor de alta presión Figura 5-23 Potencia del condensador frente a su
frente a la temperatura de salida del condensador. temperatura de salida.

5.2.5.2 Línea de media presión

Al contrario de como sucedia en los anteriores estudios la modificación de un parámetro en la línea de media
presión si provoca cambios en la línea de alta presión, aunque la línea de baja presión si se mantien
independiente, Figura 5-24.
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 37

Figura 5-24 Linea de media presión, caso 2

Para este estudio se ha modificado la potencia en el evaporador y se ha evaluado el caudal en la corriente del
condensador y del evaporador A la potencia en el compresor de alta presión y la potencia en el condensador.
La línea de media presión tiene fijado el salto entálpico pues la entrada de este corresponde con la salida del
líquido del tanque que al encontrarse en equilibrio líquido-vapor corresponderá con líquido saturado, la salida
también está fijada al conocer el título del vapor. Puesto que el salto entálpico está fijado la variación de la
potencia en el evaporador solo podría modificar el caudal que circula por el mismo.
La energía abosrbida por el refrigerante en el evaporador es cedida en el tanque, por tanto, para mantener el
equilibrio la varación de energía en la línea de media presión debe ser la misma que la absorbidaa por la línea
de alta presión. Corriente del condensador al igual que la de media también tiene fijado el salto entálpico, la
salida del tanque es vapor saturado y la entrada viene determinada por la temperatura de salida del
condensador, de nuevo los cambios en la potencia se verán reflejados en la variación del caudal.
Como ya se había visto en el estudio dos de este caso la variación del caudal en la la corriente del condensador
provocará la modificación de la potencia del compresor y de la potencia del condensador. Un aumento del
caudal provocará un aumento en ambas potencias.

0,25 48,5

0,20 48
Potencia (kW)
Caudal (kg/s)

0,15 47,5

0,10 47

0,05 46,5

0,00 46
68 70 72 74 76 78 80 68 70 72 74 76 78 80
Potencia (kW) Potencia (kW)

Figura 5-25 Evolución del caudal del condensador y Figura 5-26 Evolución de la potencia del compresor
Evaporador A frente a la potencia del evaporador A. frente a la potenci del evaporador A.
38 Casos Simulados

268
266
264

Potencia (kW)
262
260
258
256
254
68 70 72 74 76 78 80
Potencia (kW)

Figura 5-27 Evolución de la potencia del condensador frente a la potencia del evaporador A.

5.2.5.3 Línea de baja presión

En el cuarto estudio realizado para este caso se ha modificado el parámetro en la línea de baja presión, la
variación de este parámetro no aferctará a la línea de media presión que será independiente, pero si a la línea
de alta presión, Figura 5-28.

Figura 5-28 Línea de baja presión, caso 2

Se ha analizado como afecta a distintos parámetros la variación de la temperatura a la entrada del evaporador
B. A la salida del tanque el refrigerante se encuentra como líquido saturado, al pasar por la válvula se expande
entrando así en la zona bifásica. Puesto que en la zona bifásica la presión y la temperatura permanecen
constantes al fijar una la otra también queda determinada.
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 39

El título del vapor a la salida de la válvula dependerá de la temperatura y presión del fluido a la entrada y de la
presión a la salida. Al disminuir la temperatura de evaporación también disminuirá la presión, esta
disminuición de la presión irá acompañada con una mayor evaporación del refrigerante, es decir un mayor
título de vapor.
La potencia del compresor dependerá de la diferencia de presión entre la entrada y la salida. Disminuir la
presión, a la entrada, aumentaría la diferencia entre ambas pues la presión de salida ya viene determinada por
el tanque. El aumento de la diferencia de presiones supondría el aumento de la potencia del compresor.
La salida del compresor dependerá de la entrada en un compresor ideal la entropía a la entrada y a la salida se
mantiene constante. Para una mayor entropía a la entrada mayor será la temperatura a la salida, para el caso
real la temperatura a la salida será aún mayor a la ideal. La temperatura a la salida aumentará al disminuir la
presión de entrada, debido a que la entropía crece con la bajada de la presión.
Como hemos dicho la temperatura a la salida del compresor aumenta con la bajada de la presión, con el
aumento de la temperatura también aumentará la entalpía. Este aumento de entalpía supone un mayor aporte
de energía por parte de la línea de baja presión. Sabemos que todas las demás entalpías permanecen constantes
las líneas de salida están fijadas por el tanque, el líquido y gas se encuentran saturados, y las entradas por sus
títulos de vapor que permanecen constantes. Los caudales de las líneas de media y baja presión también
permanecerán constantes. La única variable que se puede ver afectada en el balance es el caudal del
condensador, un aumento de la entalpía supondrá un aumento del caudal, (5-17).

m RB h7 + m RAh5 + m RC h1 = m RC h4 + m RB h10 + m RAh6 (5-17)

250 0,120

200 0,100
Título de vapor
Presión (kPa)

150 0,080

100 0,060

50 0,040
-38 -33 -28 -23 -18 -13 -38 -33 -28 -23 -18 -13
Temperatura (ºC) Temperatura (ºC)

Figura 5-29 Evolución de la presión del Evaporador B frente Figura 5-30 Evolución del título del vapor a la entrada del
a la temperatura a la entrada del evaporador. Evaporador B en función de su temperatura.

30 140
Potencia del compresor (kW)

25 120
Temperatura (ºC)

100
20
80
15
60
10 40
5 20
0
0
80 130 180 230
80 130 180 230
Presión (kPa)
Presión (kPa)

Figura 5-31 Potencia del compresor de baja frente a la Figura 5-32 Temperatura de salida del compresor de baja
presión de entrada. frente a la presión de la entrada
40 Casos Simulados

0,206
0,204
0,202
0,2

Caudal (kg/s)
0,198
0,196
0,194
0,192
0,19
0,188
0,186
30 50 70 90 110 130
Temperatura (ºC)

Figura 5-33 Evolución del caudal del condensador frente a la temperatura a la salida del compresor de baja.

5.3 Caso 3
5.3.1 Enunciado
La instalación de refrigeración por amoniaco de la figura consta de dos evaporadores, uno (A) alimentado por una válvula
de expansión termostática de 5ºC, que produce potencia frigorífica de 50kW y otro (B) alimentado mediante una bomba
que mueve un caudal de refrigerante de m  RB = 0.15 kg/s.
En el manómetro colocado a la entrada al evaporador A marca 4.15 bar. El evaporador B se utiliza para enfriar
un caudal 6.25kg/s de etilen-glicol de -5ºC a -15ºC, utilizando una temperatura de evaporación de -20ºC. Si la
temperatura de condensación es 40ºC, se pide:
A. Calcular la efectividad [%] y el UA [W/K] del evaporador B.
B. Calcular el título del vapor [%] del refrigerante en el punto 10 y los caudales de refrigerante que
circulan por el evaporador A (, por el condensador y por el compresor de baja (.
C. Si el rendimiento volumétrico e iséntropico de los dos compresores puede estimarse en 0.7 y 0.8
respectivamente, calcular el desplazamiento volumétrico de los compresores y el COP de la
instalación.
D. Dibujar sobre el diagrama p-h adjunto todos los puntos, indicando su temperatura y entalpía.

NOTA: Suponer que el consumo de energía y el incremento de presión asociados a la bomba de líquido son
despreciables, el subenfriamiento del líquido en el condensador es nulo y no existen pérdidas de presión por
rozamiento. Calor específico del etilen-glicol 2400J/(Kg·K)
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 41

5.3.2 Solución manual

5
4

Compresor de
Alta

3
6
2

11
Compresor de
12 Baja
Evaporador A
7
1

10

Figura 5-34 Esquema 3

5.3.2.1 Planteamiento.

La eficiencia del evaporador puede obtenerse directamente mediante el salto de temperatuas en el etilen-glicol
y el salto máximo de temperatura, (5-18). A partir de la eficiencia y el método E-NTU puede obtenerse el UAB
del evaporador, (5-19).

Q fB TE − TS (5-18)
=εB =
Q fBMAX T − TEB

− AU B (5-19)
ε = 1− e m EG c pw

Haciendo balance de energía en el evaporador de baja presión se obtienela entalpía, (5-20) y (5-21), y el título
de vapor a la salida del evaporador B, (5-22). El caudal de de amoniaco que circula por el compresor de baja
se obtiene mediante el balance en el tanque de baja presión, (5-23). Puesto que la potencia del evaporador de
alta es conocida, que su entalpía de entrada es igual a la entalpía a la salida del condensador y que la válvula de
alimentación del evaporador A es una válvula de expansión termoestática de 5ºC, se puede calcular su caudal,
(5-24).

=Q fB m EG c p (TE − TS ) (5-20)

Q fB m RB (h10 − h9 )
= (5-21)
42 Casos Simulados

h10 − h9
x10 = (5-22)
h1 − h9

m RD h8 + m RB h10 =m RB h9 + m RD h1 (5-23)

Q fA m RA (h12 − h11 )
= (5-24)

Por último, para el cálculo del caudal del condesador se realiza balance de energía en el tanque intermedio,
(5-25). Los caudales del evaporador A y del compresor de baja ya se han sido determinados, para hallar el
caudal del condensadoro solo faltaría obtener la entalpía a la salida del compresor de baja. La entalpía se
obtiene mediante el rendimiento isentrópico y la entalpía a la salida del compresor si la entropía a la entrada y
salida se mantuviese constante, (5-26).

m RD h2 + m RAh12 + (m RC − m RA )h6 = m RD h7 + m RC h3 (5-25)

h2 s − h1
η SB = (5-26)
h2 − h1

5.3.2.2 Resultados

Una vez más mediante el uso de las ecuaciones y las tablas se obtienen todas las propiedades para todos los
puntos del sistema, Tabla 5-19. Con el uso de las propiedades y las ecuaciones se obtienen los caudales, Tabla
5-20 y el desplazamiento volumétrico y COP de Tabla 5-21.
P T Caudal Entalpía
[kPa] [C] [kg/s] [kJ/kg]
1 190,95 -20 0,1235 1
2 516,3 59,7 0,1235 1
3 516,3 5 0,2048 1
4 1545 99.3 0,2048 1
5 1554,81 40 0,2048 0
6 516,3 5 0,1590
7 516,3 5 0,1235 0
8 190,95 -20 0,1235
9 190,95 -20 0,15 0
10 190,95 -20 0,15
11 516,3 5 0,04585
12 516,3 10 0,04585 1
Tabla 5-19 Propiedades de los componentes del sistema, caso 3

Evaporador A Compresor Baja Condensador

Caudal (kg/s) 0.04585 0.1235 0.2048

Tabla 5-20 Caudales, caso 3.

 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 43

VTA (M3/H) VTB (M3/H) COP

256.05 396.1 3.29

Tabla 5-21 Desplazamientod volumétricos del compresor y COP, caso 3.

El circuito con el condensador de aire se representa en la Figura 5-35. Mediante la representación se puede
comprobar mejor el comportamiento mencionado anteriormente. La forma de los puntos sobre la gráfica es
similar al del caso anterior.

5 4

7 6,11 3 2

10 1
9
8

Figura 5-35 Diagrama p-h, caso 3

5.3.3 Solución HYSYS

5.3.3.1 Planteamiento

Para este tercer caso los componentes utilizados han sido, amoniaco y el etano-glycol. El amoniaco será el
componente principal el fluido refrigerante mientras que el etano-glycol será el fluido a enfriar del Evaporador
B.
Tal como se ve en la Figura 5-36, para el Evaporador A se ha utilizado Heater, en cambio en el caso del
Evaporador B al conocer el fluido que es enfriado se utiliza un intercambiador carcasa tubo. El bloque del
intercambiador a diferencia de l Heater tiene como entrada únicamente flujos de materia. El condensador al
igual que en el caso anterior se ha tomado un bloque Cooler. En este sistema también se cuenta con dos
bloques Tank que mantienen el fluido en la zona bifásica para dos presiones distintas, media y baja presión.
44 Casos Simulados

Figura 5-36 Esquema Caso 3

Figura 5-37 Datos de entrada en los flujos de materia, caso 3.

 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 45

Figura 5-38 Datos de entrada en los flujos de energía, caso 3.

5.3.5.2 Resultados

De nuevo una vez el Solver ha calculado todos los datos que faltaban podemos conocer el estado del
refrigerante en cada punto del sistema. Se puede comprobar a simple vista que los resultados obtenidos por
HYSYS, Tabla 5-22, son muy similares a los obtenidos anteriormente de forma manual, Tabla 5-19.
Name Pressure Temperature Mass flow Vapor / Phase
[kPa] [C] [kg/s] Fraction
1 187,38 -20 0,1231 1
2 516,3 67,74 0,1231 1
3 516,3 5,33 0,2030 1
4 1545 109,61 0,2030 1
5 1545 40 0,2030 0
6 516,3 5,34 0,1582 0,1358
7 516,3 5,33 0,1231 0
8 187,380 -20 0,1231 0,0869
9 187,380 -20 0,15 0
10 187,380 -20 0,15 0,7496
11 516,3 5,34 0,0448 0,1358
12 516,3 9,97 0,0448 1
13 1545 40 0,1582 0
14 1545 40 0,0448 0
15 101,3 -5 6,25 0
16 101,3 -15 6,25 0
Tabla 5-22 Propiedades de los flujos de materia de (HYSYS), caso 3
QB (kW) WA (KW) WB (KW) COP
152,49 21,86 41,93 3,17
Tabla 5-23 Potencia del intercambiador, de los compresores y COP (HYSYS), caso 3

5.3.4 Comparación
Los resultados obtenidos para el modelo teórico y del programa son muy similares, el error cometidos es de
alrededor del 1-2%. El mayor error se encuentra en el COP tal como se ha mencionado anteriormente depende
del rendimiento utilizado en cada caso.
46 Casos Simulados

Manual Hysys %ERROR

ε (%) 66.67 66.67 0

UAB (W/K) 16480,7 16750,3 1,6%

x10 (%) 75,24 74,96 0,4%

mRA (kg/s) 0,04585 0,04484 2,2%

mRC (kg/s) 0,2048 0,20303 0,9%

mRD (kg/s) 0,1235 0,12314 0,3%

COP 3,29 3,17 3,5%

Tabla 5-24 Comparativa de resultados, caso 3

5.3.5 Análisis de sensibilidad

Este caso posee dos tanques y tres presiones distintas para el refrigreante. De nuevo para comprobar la
dependencia del sistema de las variables a estudiar se eliminan. Se han elaborado dos estudios distintos.

5.3.5.1 Línea de alta y media presión

En este estudio se ha la variable de independiente utilizada ha sido la temperatura a la salida del evaporador A.
El cambio de este variable afecta a toda la corriente que circula por el condensador en todo su recorrido.

Figura 5-39 Estudio alta presión, caso 3

Al variar la temperatura se modifican los caudales pertenecientes al evaporador A y a la línea 6, sin embargo,
el caudal que circula por el condensador se mantiene constante. Para un aumento de temperatura el caudal del
evaporador siminuira mientras el caudal 6 aumenta, ambos aumentarán o disminuirán lo mismo. También es
cierto que para el incremento o disminución de un grado los cambios producidos en el caudal pueden llegar a
ser casi despreciables. La disminución del caudal del evaporador se debe al aumento de entalpía que conlleva
el aumento de temperatura. Mediante los cambios en el caudal 6 se compensan las variaciones de energía de
forma que el balance en el tanque se mantenga constante.
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 47

0,04525 0,15875

0,04500 0,1585
Caudal (kg/s)

Caudal (kg/s)
0,04475 0,15825

0,04450 0,158

0,04425 0,15775
5 8 11 14 17 5 8 11 14 17
Temperatura (ºC) Temperatura (ºC)

Figura 5-40 Evolución del caudal a la salida del Figura 5-41 Evolución del caudal 6 frente a la
Evaporador A frente a su temperatura. temperatura de salida del evaporador A.

5.3.5.1 Línea de baja presión

Para este estudio se ha analizado como afectaría al sistema los cambios en el caudal del etileno-glicol. Tal y
como puede verse en Figura 5-42, las únicas líneas del sistema que son independientes del caudal son la
entrada al Evaporador B y la línea del Evaporador A.

Figura 5-42 Estudio baja presión, caso 3

Para analizar este estudio debe empezarse por el Evaporador B, al definir el caudal, tres de las líneas de
materia que entran y salen de este intercambiador quedan completamente determinadas. La única incógnita
sería el título del vapor del refrigerante a la salida del mismo. Puesto que el interés en este estudio del
intercambiador reside en enfriar la corriente de etileno-glicol de -5ºC a -15ºC sea cual sea su caudal, la energía
extraída por el refrigerante deberá adaptarse en función del caudal que circule. Un mayor caudal supondría una
mayor energía para absorber por parte del refrigerante. Si el refrigerante absorbe más energía incrementará su
entalpía y por tanto su título de vapor aumentará y podría llegarse al caso de que se produjera un
recalentamiento.
Los cambios de entalpía en la entrada del botellón provocarán cambios en el caudal de media presión. Esto se
48 Casos Simulados

justifica mediante el balance de energía en el botellón. Para mantener el equilibrio la energía que entra debe ser
igual a la que sale. El único parámetro que puede variar ante el cambio de la energía a la entrada es el caudal
que circula por el compresor de media presión. Un aumento de energía en la entrada supondrá un aumento del
caudal. La variación de este caudal afectará a la potencia del compresor que deberá aumentar para un mayor
caudal o disminuir para uno menor.
El botellón de alta presión se comporta del mismo modo que lo hacía el de baja. Una variación de la energía
pocedente del botellón anterior modificará el caudal del compresor de alta presión y la potencia del compresor
se verá afectada por los cambios del caudal.
0,9 0,15000
0,14000
Título de vapor

0,8

Caudal (kg/s)
0,13000
0,7 0,12000
0,11000
0,6
0,10000
0,5 0,09000
4,5 5,5 6,5 7,5 0,55 0,65 0,75 0,85 0,95
Caudal (kg/s) Título de vapor

Figura 5-43 Evolución del título del vapor a la salida Figura 5-44 Evolución del caudal compresor de baja
del intercambiador frente caudal etileno-glicol. frente al título de vapor de salida del intercambiador.
25
0,23
23
Potencia (kW)

0,21
Caudal (kg/s)

21
0,19
19
0,17
17
0,15
15
0,09000 0,11000 0,13000 0,15000
0,09000 0,11000 0,13000 0,15000
Caudal (kg/s)
Caudal (kg/s)

Figura 5-46 Evolución del caudal del condensador

Figura 5-45 Evolución de la potencia del compresor de
frente al del compresor de baja presión.
baja presión frente a su caudal.
50

45
Potencia (kW)

40

35

30
0,17 0,18 0,19 0,2 0,21 0,22 0,23
Caudal (kg/s)

Figura 5-47 Potencia del compresor de alta frente a su caudal

 6 CONCLUSIONES

L uso de sistemas de producción de frío resultan fundamentales hoy en día tanto para usos industriales

E como en los hogares, cada vez encontramos más utilidades a esta tecnología. Como ingenieros debemos
procurar que los sistemas sean cada vez más eficientes y lograr así un ahorro de energía y recursos.
Para lograr estos objetivos los simuladores, en nuestro caso HYSYS juegan un papel fundamental. En este
trabajo a partir de la simulación de problemas académicos, se ha podido comprobar la gran fiabilidad del
programa. Ha permitido obtener de forma rápida y precisa los datos que requeríamos, los cálculos manuales a
pesar de ser fiables se encuentran siempre sujetos a los errores humanos lo que en ocasiones implica la
repetición de los mismos.
Los análisis de sensibilidad realizados permiten conocer aquellas variables más críticas en cada proceso y
como afectan a las demás. Para un proceso de frío las variables más importantes son las presiones y
temperaturas que determinan la condensación y evaporación del refrigerante. En cambio, una pequeña
diferencia en la presión o temperatura de fluidos como el agua utilizados también en la condensación y
evaporación apenas producen cambios.
A pesar de ser un programa relativamente sencillo es importante conocer bien el proceso y prácticar con
HYSYS. Esto es de especial importancia con la elección de los métodos termodinámicos, este ha sido uno de
los mayores problemas en la realización de los casos analizados. Este trabajo pretende servir como guía para
aquellos estudiantes, docentes o ingenieros que deseen adentrarse en el uso del programa para la realización de
simulaciones.

49
50
Conclusiones

50
 Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 51

7 DESARROLLO FUTURO
ste trabajo es una base del uso de HYSYS para la resolución de instalaciones de frío. Se ha podido

E comprobar que el programa puede resolver los casos teóricos. Para futuros proyectos sería recomendable
la simulación de sistemas reales que incluyeran pérdidas de carga en redes y equipos, factores de
ensuciamiento, diferentes compresores o bombas. Como sabemos en la realidad, en la industria, los equipos no
se comportan igual que los modelos ideales. Para empezar las condiciones de contorno, como la temperatura
del aire de un aerocondensador, no suelen mantenerse constantes. Por tanto, la simulación de problemas reales
resultaría más compleja pero también más interesante.
Para lograr resultados fiables sería necesario contar con una gran base de datos procedentes de la instalación a
simular, solo de esta forma sería posible validar completamente los casos que se realizaran. Una vez que los
casos hubiesen sido validados sería posible la realización de estudios de sensibilidad, tal como se ha visto en
los casos teóricos. Estos análisis podrían incluirse en nuevos trabajos más enfocados a la optimización del
proceso.
52
Desarrollo futuro

52
Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 53

ANEXOS
A. Anexo: Enunciados de los casos
54
Anexos

B. Anexo

54
Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 55

C. Anexo
56
Anexos

D. Anexo Tablas de propiedades R717

56
57
58
Anexos

58
Simulación de instalaciones frigoríficas usando ASPEN-HYSYS 59
60
Anexos

E. Anexo: Diagrama P-h del refrigerante R717

60
 BIBLIOGRAFÍA

[1] León Cohen, Diseño y simulación de procesos químicos, 2001.

[2] M. Cuevas Aranda, Introducción a la simulación en Ingeniería Química, Universidad de Jaén, 2014

[3] J. Doria, G. Gondra, Pedro F. Hernández, L. Ruiz de Guana, S. Makazaga, Instalaciones frigoríficas, 1995

[4] Samira Khellaf Modelado y simulación de una planta de cogeneración mediante el simulador de procesos
HYSYS, Escuela Técnica Superior de Ingeniería (Universidad de Sevilla), 2016.

[5] Tillner-Roth, Harms-Watzenberg, and Baehr, Eine neue Fundamentalgleichung fur Ammoniak, DKV-
Taugungsbericht 20:167-181,1993.

[6] R. González Flacón, D. Velázquez Alonso, Tablas y gráficos de tecnolgía energética, Versión 1, 2014

[7] Aspen HYSYS Help, 28 de Abril 2016

[8] Ángel L. Villanueva Perales, Métodos termodinámucos en simuladores de procesos, Simulación y

optimización de procesos químicos, Dpto Ingeniería Química y Ambiental, Universidad de Sevilla,
2013.

[9] Relación de problemas de producción de frío, Instalaciones y Máquinas Hidráulicas y Térmicas.

[10]Sandra, N. Rodríguez, Introducción al entorno de HYSYS, 2013.

[11]Stewart, Jacobsen, Renoncello, Cap 2.2 Ammonia, Manual: PROPATH, A program package for
thermophysical properties of fluids, Version 3, PROPATH GROUP, 2001

61
62
Bibliografía

62