Manual EES (Termodinámica Aplicada)
Manual EES (Termodinámica Aplicada)
Manual EES (Termodinámica Aplicada)
PRÁCTICAS DE TERMODINÁMICA APLICADA
ÁREA DE MAQUINAS Y MOTORES TÉRMICOS
DEPARTAMENTO DE INGENIERÍA TÉRMICA Y FLUIDOS
1
APLICACIÓN DE EES A LA RESOLUCIÓN DE PROBLEMAS DE TERMODINÁMICA EN AULA
INFORMÁTICA
1. Objetivo
El objetivo de la práctica es la presentación del programa de simulación EES® al
alumno como una herramienta válida para resolver problemas de Termodinámica Aplicada.
2. Introducción
EES es el acrónimo de “Equation Engineering Solver”, es decir un programa para
resolver ecuaciones en ingeniería. Además de permitir la resolución de cualquier tipo de
ecuaciones.
El programa dispone de una gran base de datos de las propiedades termofísicas y
termodinámicas de muchas sustancias de uso en Ingeniería (gases ideales, vapor de agua,
refrigerantes, etc).
Hay dos grandes diferencias entre EES y otros programas similares:
1. EES permite a las ecuaciones con variables desconocidas colocarlas en cualquier orden.
EES reordena automáticamente las ecuaciones para obtener una solución eficiente.
2. EES proporciona muchas funciones integradas de propiedades matemáticas y
termofísicas útiles para los cálculos de ingeniería. Las propiedades de transporte
también están disponibles para una gran cantidad de sustancias.
3. Menús y barras de herramientas
El interfaz de usuario del programa es similar al resto de programas desarrollados en
entorno de Windows, con una barra de menús desplegables en la superior, una barra de
herramientas de acceso rápido, y una zona de trabajo.
Figura 3.1. Menús y barras de herramientas.
2
En la ventana “Equations Window” es donde se deberán de introducir las ecuaciones
necesarias para el problema en cuestión.
4. Librería de funciones termodinámicas
Para la determinación de las propiedades termodinámicas de un fluido se pueden
emplear funciones específicas para ello que el solver EES tiene integradas. Para acceder a ellas,
debemos de abrir el menú “Function Information”, donde seleccionamos el tipo de variable a
determinar y el fluido correspondiente. Una vez seleccionada, pulsamos en el botón “PASTE”, y
la introduce en la ventana “Equations Windows”.
Figura 3.2. Librería de funciones termodinámicas.
En el argumento de la función se determina el tipo de fluido de trabajo, así como las
propiedades de entrada necesarias para su determinación.
4.1. Argumentos de una función termodinámica
4.1.1. Nombre de la sustancia
El primer argumento de una función termodinámica es el nombre de la sustancia. La
siguiente tabla muestra todas las sustancias que EES tiene incluidas en su librería:
3
Figura 3.3. Sustancias de EES®.
Las sustancias incluidas en la 1ª columna están modelas como gases ideales, mientras
que el resto están modeladas como fluidos reales. Por tanto, mientras “N2” es
nitrógeno modelado como gas ideal, “Nitrogen” es el mismo gas modelado como
fluido real.
4.1.2. Funciones termodinámica
A continuación se presentan las funciones principales de termodinámica incluidas en
EES®.
AcentricFactor: factor acéntrico.
Conductivity: conductividad térmica. Además del nombre del fluido, para gases ideales
requiere la temperatura como argumento. Para fluidos reales requiere 2 argumentos
que pueden ser cualquier par de propiedades independientes. Para aire húmedo
(AirH2O) se deben proporcionar 3 argumentos, temperatura, presión y humedad
específica (o relativa).
CP, CV: calores específicos a presión y volumen constante. Para sustancias puras que
se comporten como gases ideales sólo es necesario especificar la Tª. Para fluidos reales
es necesario proporcionar 2 propiedades independientes. En condiciones de
saturación, presión y Tª no son independientes, por lo que habrá que proporcionar
además el título de vapor o alguna otra propiedad dependiente del título.
Density: densidad. Para sustancias puras se necesitan 2 argumentos, y 3 para aire
húmedo.
DewPoint: temperatura de rocío. Sólo se puede usar con aire húmedo (AIRH2O) como
sustancia. Requiere 3 argumentos, temperatura, presión y humedad específica (o
relativa).
Enthalpy: entalpía específica. Para gases ideales solo requiere un argumento (Tª o
energía interna), para fluidos reales requiere 2 argumentos y para aire húmedo 3.
Entropy: entropía específica. Requiere 2 argumentos para todas las sustancias puras y
3 para el aire húmedo.
4
Fugacity: fugacidad.
HumRat: humedad específica. Sólo se puede emplear esta función con aire húmedo y
requiere 3 argumentos, uno de los cuales debe ser obligatoriamente la presión.
IntEnergy: energía interna específica. Para gases ideales solo requiere un argumento
(Tª o energía interna), para fluidos reales requiere 2 argumentos y para aire húmedo 3.
IsIdealGas: se trata de una función con un solo argumento (el nombre del fluido) que
devuelve el valor 1 si el fluido se comporta como gas ideal o 0 en caso contrario.
MolarMass: peso molecular.
Pressure: presión. Esta función requiere 2 argumentos y no es aplicable al aire
húmedo (aunque la presión para el aire húmedo puede obtenerse empleando alguna
de las funciones aplicables al aire húmedo que tiene la presión como argumento).
P_Crit: presión crítica.
Prandtl: número de Prandtl. Requiere la Tª como argumento para gases ideales y
presión y Tª para sustancias reales.
Quality: Calidad o título de vapor. Requiere 2 propiedades independientes como
argumento (Tª y presión no son propiedades independientes para sustancias
saturadas).
Relhum: humedad relativa. Esta función requiere 3 argumentos además del nombre
de la sustancia. Estos argumentos son la temperatura, la presión y otra propiedad
independiente de las 2 anteriores.
Specheat: calor específico a presión constante. Además del nombre del fluido, para
gases ideales requiere la temperatura como argumento. Para fluidos reales requiere
presión y Tª.
SoundSpeed: velocidad del sonido a través de un fluido.
SurfaceTension: tensión superficial en la interface líquido vapor de un fluido saturado.
Como argumentos sólo requiere, además del nombre del fluido, la Tª del mismo.
Temperature: temperatura. Para gases ideales puede requerir 1 o 2 argumentos
mientras que para sustancias reales siempre requiere 2 argumentos.
T_Crit: temperatura crítica.
Volume: volumen específico. Para las sustancias puras es necesario introducir 2
argumentos y para el aire húmedo se precisan 3.
V_Crit: volumen específico crítico.
Wetbulb: temperatura de bulbo húmedo. Esta función requiere 3 argumentos además
del nombre de la sustancia. Estos argumentos son la temperatura (o entalpía), la
presión y la humedad relativa (o humedad específica o Tª de rocío).
5
Viscosity: viscosidad. Además del nombre del fluido, para gases ideales requiere la
temperatura como argumento. Para fluidos reales requiere 2 argumentos que pueden
ser cualquier par de propiedades independientes. Para aire húmedo (AirH2O) se deben
proporcionar 3 argumentos, temperatura, presión y humedad específica (o relativa).
4.1.3. Identificación de argumentos
Cualquier argumentos se identifica por una única letra (en mayúsculas) seguida del
signo “=”. Las letras reconocidas como argumento son:
B: Tª de bulbo húmedo.
D: Tª de rocío.
H: entalpía específica.
P: Presión.
R: humedad relativa.
S: entropía específica.
T: temperatura.
U: energía interna específica.
V: volumen específico.
W: humedad específica.
X: título.
5. Definición de unidades
Al iniciar la aplicación, se debe de definir el sistema de unidades empleado, para ello se
tiene que acceder al menú OPTIONS ‐> UNIT SYSTEM
Figura 3.4. Definición de unidades.
Para activar la función de que EES verifique que no existen inconsistencias en las
unidades, debemos marcarlo en la solapa “Options”.
6
Figura 3.5. Verificación de inconsistencias.
6. Diagrama interactivo
El EES permite trabajar con una herramienta interesante: el diagrama interactivo del
proceso (Diagram Window). Este diagrama no es generado por el programa sino que es el
alumno quién tiene la posibilidad de poder llevarlo a cabo. A la hora de desarrollar el
diagrama existe la posibilidad de que en éste aparezcan las variables que intervienen en el
proceso. Así cuando el programa resuelva el problema introducido en la “Equations Window”
se tiene la posibilidad de que los resultados aparezcan en la “Diagram Window“.
Además se pueden introducir en la “Diagram Window” algunas variables necesarias
para la resolución del problema. Si se hace esto, los valores de las variables pueden ser
modificados a través de la “Diagram Window“ y cuando el programa resuelva el problema,
identifica como nuevas entradas los valores especificados en el diagrama interactivo.
Figura 3.6. Diagram Window.
Así, el alumno una vez que haya generado el diagrama interactivo puede optar por
modificar las variables en la ventana del diagrama, y observar cómo se ven afectados los
resultados con la modificación de dichos parámetros. Así que, aunque en un principio esta
7
opción de resolución pueda parecer más compleja (ya que hay que desarrollar el diagrama), a
posteriori permite poder observar con facilidad como varían determinadas incógnitas al variar
determinadas variables. Y por tanto poder sacar conclusiones de la observación de los cambios
que se lleven a cabo.
7. Obtención de curvas paramétricas
El programa también permite optar por dejar una variable sin definir de manera que el
número de incógnitas se vea aumentado en la unidad. El usuario podrá decidir qué variable
quiere que sea paramétrica, y observar cómo varían otras variables en función de ésta. La
variación de variables en función de la paramétrica se obtiene a partir de tablas paramétricas.
7.1. Tabla Paramétrica
El programa permite generar una tabla paramétrica bien pulsando el botón “New
Parametric Table” de la barra de herramientas o bien entrando dentro del Menú Tables. La
ventana que aparece permite elegir la variable paramétrica y las que el programa tiene que
calcular a partir de ésta.
En el caso de la siguiente figura, la variable de entrada es T1 (aparece en color azul), y las
variables calculadas COP y Qevap (aparecen en color negro). Ordenando al programa que
genere la tabla paramétrica pulsando el botón “Solve Table“ de la barra de herramientas se
obtiene la siguiente ventana:
Figura 3.7. Tabla paramétrica.
7.2. Curvas Paramétricas
Como ya se ha dicho, el EES permite observar esta variación también en forma de
gráfica. Por tanto si el alumno quiere representar los valores de la tabla paramétrica en una
curva le basta con elegir la opción “New Plot Window” de la barra de herramientas e indicarle
cuales son las variables que quiere representar en la curva.
8
Siguiendo con el mismo ejemplo si se decide representar las curvas paramétricas de las
variables que aparecen en la tabla el resultado será:
Figura 3.8. Curva paramétrica.
9
8. Problemas propuestos
Resolver los siguientes problemas utilizando EES y presentar una memoria con las ecuaciones
introducidas, variables definidas, gráficos y comentar los resultados obtenidos.
Problema 1: Evaluación de propiedades utilizando las funciones termodinámicas de EES.
“Unidades: SI.”
“Temperatura: ºC.”
“Presión: bar.”
“Energía: kJ.”
Calcular las siguientes funciones termodinámicas.
Conductividad Térmica
K1=CONDUCTIVITY(Air;T=100)
K2=CONDUCTIVITY(Ammonia;T=100;P=200)
K3=CONDUCTIVITY(R134a;T=100;P=200)
K4=CONDUCTIVITY(Water;T=100;P=200)
Densidad
d1=DENSITY(Air;T=100;P=100)
d2=DENSITY(Steam;T=100;P=100)
d3=DENSITY(AirH2O;T=100;P=100;w=10)
Viscosidad
vis1=VISCOSITY(Air;T=300)
vis2=VISCOSITY(R134a;T=40;P=10)
vis3=VISCOSITY(R134a;T=40;x=1)
vis4=VISCOSITY(Steam;T=100;P=20)
vis5=VISCOSITY(AirH2O;T=70;P=14;R=0,7)
Peso Molecular
M_CO2=MOLARMASS(CO2)
M_N2=MOLARMASS(N2)
M_Aire=MOLARMASS(Air)
M_Amoniaco=MOLARMASS(Ammonia)
¿La sustancia es un Gas ideal? Si=1, No=0
GI1=ISIDEALGAS(Air)
GI2=ISIDEALGAS(O2)
GI3=ISIDEALGAS(N2)
GI4=ISIDEALGAS(R134a)
Calores específicos a presión y volumen constante
Cp1=CP(Air;T=100)
Cv1=CV(Air;T=100)
Cp2=CP(Ammonia;T=100;P=200)
Cv2=CV(Ammonia;T=100;P=200)
Cp3=CP(AirH2O;T=70;P=14;R=0,7)
Cv3=CV(AirH2O;T=70;P=14;R=0,7)
10
Entalpía específica
h1=ENTHALPY(Air;T=100)
h2=ENTHALPY(Steam;T=100;P=100)
h3=ENTHALPY(AirH2O;T=70;P=14;w=1)
Entropía específica
s1=ENTROPY(AirH2O;T=70;P=14;R=0,7)
s2=ENTROPY(O2;T=400;P=100)
s3=ENTROPY(Air;T=526,8;P=100)
Energía Interna específica
u1=INTENERGY(Air;T=200)
u2=INTENERGY(Steam;T=200;P=50)
u3=INTENERGY(AirH2O;T=70;P=14; w=0,7)
Presión y presión crítica
P1=PRESSURE(Steam;T=400;H=3263,9)
P2=PRESSURE(R134a;T=100;S=1,1)
P3=PRESSURE(Air;T=526,8;s=5,399)
Pc_agua=P_CRIT(Water)
Pc_R134a=P_CRIT(R134a)
Temperatura y temperatura crítica
T1=TEMPERATURE(Air;H=821,95)
T2=TEMPERATURE(Air;U=592,30)
T3=TEMPERATURE(Air;s=5,399;P=100)
Tc_agua=T_CRIT(Water)
Tc_R134a=T_CRIT(R134a)
Temperatura del punto triple
T_triple_agua=T_TRIPLE(Water)
T_triple_R134a=T_TRIPLE(R134a)
Volumen específico y volumen crítico
v1=VOLUME(Air;T=100;P=100)
v2=VOLUME(Steam;T=400;P=10)
v3=VOLUME(AirH2O;T=70;P=14;R=0,7)
vc_agua=V_CRIT(Water)
vc_R134a=V_CRIT(R134a)
Título o calidad de vapor
x1_agua=QUALITY(Water; H=1500;P=100)
x2_R134a=QUALITY(R134a; H=260;P=10)
x3_R22=QUALITY(R22;H=320;P=10)
Temperatura de Rocío
TR1=DEWPOINT(AirH2O;T=21;P=1;w=0,011)
TR2=DEWPOINT(AirH2O;T=21;P=1;R=0,7)
TR3=DEWPOINT(AirH2O;T=21;P=1;B=17,33)
Humedad específica solo a la sustancia AirH2O ‐‐>aire húmedo
w1=HUMRAT(AirH2O;T=70; P=14; R=1)
11
w2=HUMRAT(AirH2O;T=70; P=14; H=2697)
Humedad relativa solo a la sustancia AirH2O ‐‐>aire húmedo
R1=RELHUM(AirH2O;T=70;P=14;w=0,01418)
R2=RELHUM(AirH2O;T=70;P=14;H=108)
R3=RELHUM(AirH2O;T=70;P=14;B=70)
Temperatura del bulbo húmedo
TB1=WETBULB(AirH2O;T=21;P=1; R=0,7)
TB2=WETBULB(AirH2O;T=70;P=14;R=0,7)
TB3=WETBULB(AirH2O;T=70;P=14;D=30)
Conversión a otras unidades
C_p=4,18 [kJ/kg‐C]
C_p_1=C_p*Convert(kJ/kg‐C; kJ/kg‐K)
C_p_2=C_p*Convert(kJ/kg‐C; Btu/lbm‐F)
C_p_3=C_p*Convert(kJ/kg‐C; Btu/lbm‐R)
C_p_4=C_p*Convert(kJ/kg‐C; kCal/kg‐C)
Problema 2. En un dispositivo de cilindro – émbolo con carga variable y con una rueda de
paletas integrada al cilindro, hay aire. Al principio está a 500 kPa y 27 ºC.
Se hace girar la rueda de paletas meditante un motor eléctrico externo, hasta que se ha
transferido al aire la energía de 50 kJ/kg en forma de trabajo. Durante el proceso se transfiere
calor para mantener la temperatura del aire, y al mismo tiempo se triplica el volumen del gas.
Calcular el calor intercambiado. Suponer el aire gas ideal.
Problema 3. Un dispositivo cilindro‐émbolo contiene 1.5 kg de N2, inicialmente a 100 kPa y
17ºC. Entonces se comprime lentamente en un proceso politrópico durante el cual pV1.3 =
constante, hasta que el volumen se reduce a la mitad. Determinar el trabajo efectuado y la
transferencia de calor para este proceso. Dibujar el proceso en un diagrama p‐V e investigar el
efecto de n en el trabajo y el calor transferido. Hacer variar el exponente politrópico entre 1.1
y 1.6 y trazar las gráficas correspondientes al trabajo y al calor transferido en función del
exponente politrópico y describir los resultados.
Problema 4. Se va a calentar un recinto rígido de 4 mx5mx6m con un calentador eléctrico. Se
desea que el calentador pueda elevar la temperatura de 5 a 25 en 11 minutos. Suponiendo que
no hay pérdidas de calor en el recinto (adiabático) y que la presión atmosférica es de 100 kPa,
calcular la potencia requerida por el calentador. Suponer que los calores específicos son
constantes a la temperatura ambiente.
Problema 5. Expansión de R134a a través de una válvula (I). Un flujo de refrigerante R‐134a,
que se encuentra a 60°C y 750 kPa circula a través de una tubería de 0,015 m2, a una velocidad
de 15 m/s. El refrigerante atraviesa una válvula de expansión de forma que, a la salida de la
misma, tenemos una presión de 300 kPa. Determinar la temperatura del fluido en las
condiciones de salida, el flujo másico y la velocidad a la salida de la válvula si la sección de
salida tiene también 0,015 m2.
12
Problema 6: Expansión de R134a a través de una válvula (II). Repetir el caso anterior,
estableciendo la sección de salida de 0,010 m2, y variando la sección de entrada en los valores:
0,010 m2, 0,015 m2; 0,020 m2, 0,025 m2, 0,03m2
Representar en una gráfica la variación de la temperatura de salida en función de la sección de
entrada.
Problema 7: Flujo volumétrico de un secador (I). Un secador de pelo es básicamente un
conducto en el que se colocan unas capas de resistencias eléctricas y un pequeño ventilador
empuja el aire hacia dentro y pasándolo por las resistencias donde es calentado. El aire entra
en un secador de pelo de 1200 W a 100 kPa y 22 °C y sale a 47 °C, tal como se muestra en la
figura siguiente. El área de sección transversal del secador de pelo a la entrada y a la salida es
de 60 cm2. Despreciando la potencia consumida por el ventilador y las pérdidas de calor a
través de las paredes del secador de pelo, determinar (a) el flujo de volumétrico de aire en la
entrada y (b) la velocidad del aire en la salida.
“Unidades: SI.”
“Temperatura: ºC.”
“Presión: kPa.” A1=60 cm2
“Energía: kJ
Figura 3.9. Flujo volumétrico del secador de pelo.
Suposiciones: 1. Este es un proceso de flujo estable ya que no hay cambio con el tiempo. 2. El
aire es un gas ideal con calores específicos constantes a temperatura ambiente. 3. Los cambios
de energía cinética y potencial son despreciables. 4. La potencia consumida por el ventilador y
las pérdidas de calor son insignificantes.
Considere: La constante R del aire es igual 0,287 kPa.m3/(kg∙K) y el calor específico a presión
constante de aire a temperatura ambiente es cp = 1.005 kJ/(kg∙K).
Problema 8: Flujo volumétrico de un secador (II). Repetir el caso anterior, estableciendo la
sección de entrada de 60 cm2 y variando la sección de salida en valores de 25 cm2 a 75 cm2,
con incrementos de secciones de 5 cm2.
Representar en una gráfica la variación de la velocidad de salida del aire en función de la
sección de salida del secador.
Problema 9: Trabajo de una turbina. Una turbina opera en situación estacionaria tal y como
muestra la figura. El vapor de agua entra a 3 Mpa y 400 ºC con un flujo volumétrico de 85
m3/min. Parte del vapor se extrae de la turbina a una presión de 5 bar y una temperatura de
180 ºC. El resto se expande hasta una presión de 0.06 bar y abandona la turbina con un flujo
másico de 40000 kg/h y un título del 90 %. El calor disipado por la turbina es de 200kW. La
temperatura que rodea a la turbina es la ambiental a 25ºC. Las variaciones de energía cinética
y potencial pueden despreciarse. Calcular: a) El diámetro, en m, del conducto por el que se
extrae el vapor si su velocidad es de 20 m/s. b) La potencia desarrollada, en kW.
13
“Unidades: SI.”
ܳሶ=200kW “Temperatura: ºC.”
Tamb=25 ºC “Presión: bar.”
“Energía: kJ
Figura 3.10. Turbina de vapor.
Problema 10: Ciclo de refrigeración
Disponemos de una cámara frigorífica que opera en las siguientes condiciones, siguiendo un
ciclo frigorífico:
- Temperatura del aire ambiente exterior: 25ºC
- Temperatura del aire en la cámara frigorífica: ‐10ºC
- Diferencia mínima de temperatura en el condensador: 8ºC
- Diferencia mínima de temperatura en el evaporador: 5ºC
El fluido refrigerante empleado es el R‐134a. Las condiciones del fluido a la entrada del
compresor son vapor saturado, y en la salida del condensador se encuentra en condiciones de
líquido saturado. Determinar el COP y COPcarnot. Repetir el cálculo considerando que existe un
rendimiento isoentrópico en el compresor de 0,80.
Presentar los resultados en un esquema similar al de la Figura 3.6.
14