APLICACION DE LA COMPUTACION AL CALCULO DE LAS PROPIEDADES

TERMODINMICAS
Saavedra L., Francisco D., M Sc. Ing. Mecnico Universidad de Carabobo. Facultad de Ingeniera.
Escuela de Ingeniera Mecnica. Departamento de Trmica y Energtica. Brbula, Carabobo 2008,
Venezuela. TeleFax 0241- 8672814.
Resumen
Desde el ao 1982 se han venido desarrollando en la Escuela de Ingeniera Mecnica de la Universidad de
Carabobo algoritmos de clculo para determinar las Propiedades Termodinmicas de diferentes Sustancias
Puras. Hasta ahora se han presentado nueve (9) Trabajos (Cristiano, Saavedra y otros 1984; Agudelo y Durn
1985; Coln 1985; Saavedra y Cristiano 1987; Cristiano y Saavedra 1988; Saavedra y Cristiano 1991; Ramis y
Vera 1992; Saavedra y Cristiano 1993; Saavedra 1993) los cuales han centrado su atencin en el desarrollo de
Programas de Computacin para calcular las Propiedades Termodinmicas de esas Sustancias Puras. Este
Trabajo presenta un Programa de Computacin, desarrollado y compilado en BASIC PDS ver. 7.1, capaz de
calcular las Propiedades Termodinmicas; Temperatura o Presin, Volumen Especfico, Energa Interna
Especfica, Entalpa Especfica y Entropa Especfica, de las Sustancias Puras siguientes: Aire Seco, Alcohol
Proplico, Agua, Amonaco, Argn, Butano, Dixido de Carbono, Etano, Etileno, Heptano, Hexano, Hidrgeno,
Iso-Butano, Iso-Pentano, Mercurio, Metano, Nitrgeno, Octano, Oxgeno, Pentano, Propano, Propileno y
Refrigerantes 11, 12, 13, 14, 22, 23, 114, C-318, 500, 502 y 503. Este Programa es aplicable en cualquier
Micro-Computador que acepte el MS-DOS y tabula las Propiedades Termodinmicas en los intervalos de
Temperatura y Presin que el usuario desee. El error que se comete en el clculo de las Propiedades es menor al
2%.
1. INTRODUCCION
El objetivo principal de este Artculo es presentar un procedimiento para determinar las propiedades
termodinmicas en los estados de saturacin y sobrecalentamiento y un programa de computacin capaz de
calcular y tabular esas Propiedades para una Sustancia Pura seleccionada en los intervalos e incrementos de
presin y temperatura que el usuario seleccione y a una temperatura de referencia preestablecida.
2. PROCEDIMIENTO GENERAL DE CALCULO PARA DETERMINAR LAS PROPIEDADES
TERMODINMICAS
Para obtener las propiedades termodinmicas en los estados de saturacin y sobrecalentamiento es necesario
conocer un conjunto de ecuaciones empricas o modelos matemticos obtenidas a partir de experiencias de
laboratorio. Ese conjunto de ecuaciones es particular de cada sustancia pura.
2.1.- Procedimiento de clculo para determinar las propiedades termodinmicas en el estado de
saturacin.
Para determinar las propiedades termodinmicas en el estado de saturacin (Reynolds 1979; Saavedra 1993) se
deben conocer cuatro ecuaciones empricas; la ecuacin que determina el calor especfico a volumen especfico

constante para densidades o presiones muy bajas, la ecuacin de estado Virial que representa el comportamiento
de la sustancia pura en la fase gaseosa, la ecuacin que representa la variacin de la presin de vapor con la
temperatura, presin de saturacin, y la ecuacin que determina el valor de la densidad del lquido saturado para
diferentes temperaturas de saturacin, es decir:
Cv(T),P = f(,T), Psat = g(T) y f .= (T)
respectivamente.
Con el propsito de ilustrar este procedimiento de clculo, se muestra en la Figura 1 un diagrama de flujo del
mismo.

Figura l.- Diagrama de Flujo para Determinar las Propiedades Termodinmicas en el Estado de
Saturacin.
Los pasos que deben seguirse para cumplir este procedimiento se exponen a continuacin:
1.- Se selecciona un valor de la temperatura de saturacin, Ts.
2.- Se calcula la presin de saturacin.
3.- Se determina la densidad del lquido saturado, f
4.- Se sustituye f en la expresin vf = 1 /f y se determina el volumen especfico del lquido saturado.
5.- Se calcula el volumen especfico del vapor saturado a partir de la ecuacin de estado Virial. Como el
volumen especfico del vapor saturado (vg = 1/f) est implcito en la ecuacin de estado Virial, se necesita de
un mtodo de ensayo y error (Saavedra 1993) para determinar ese volumen especfico (vg). Para ello se procede
de la forma siguiente:
5.a.- Con la temperatura de saturacin (Ts) y la presin de saturacin (Psat), se calcula el volumen especfico del
vapor saturado (vg) a partir de la ecuacin de estado para un gas ideal.

5.b.- Se sustituyen los valores de la temperatura de saturacin (Ts) y el volumen especfico del vapor saturado
(vg) en la ecuacin de estado Virial y se determina una presin de tanteo (Ptan).
5.c.- Se determina el error o la convergencia a partir de:

5.d.- Se verifica si el error es menor que 0.05%. Si es as, se concluye el tanteo, sino, se vara el valor del
volumen especfico (vg) hasta que se consiga un error menor al 0.05%.
6.- Se determina el valor del volumen especfico de vaporizacin (vfg) a partir de:

7.- Se calcula el valor de la entalpa especfica de vaporizacin a partir de la ecuacin de Clausius-Clapeyron.

8.- Se determina la entropa de vaporizacin (sfg) a partir de:

9.- Se obtiene el valor de la energa interna especfica del vapor saturado (ug), mediante la ecuacin siguiente:

donde Cv(Ts) es la ecuacin que determina el valor del calor especfico a volumen especfico constante y uo es
la energa interna especfica de referencia.
10.- Se calcula el valor de la entropa especfica del vapor saturado (sg) a partir de la ecuacin siguiente:

donde Cv(Ts) es la ecuacin que determina el valor del calor especfico a volumen especfico constante y s0 es la
entropa especfica de referencia. La ecuacin anterior tambin se puede escribir as:

11.- Se determina la entalpa especfica del vapor saturado a partir de:

12.- Se calculan la entalpa especfica y la entropa especfica del lquido saturado, usando las relaciones
siguientes:
hf = hg - hfg
y
Sf = Sg Sfg
Siguiendo este procedimiento de clculo se pueden determinar las propiedades termodinmicas para el estado
de saturacin de cualquier sustancia pura.

2.2.- Procedimiento de clculo para determinar las propiedades termodinmicas en el estado de

sobrecalentamiento.
Para obtener las propiedades termodinmicas en el estado de sobrecalentamiento (Reynolds 1979; Saavedra
1993) se deben conocer dos de las cuatro ecuaciones empricas, es decir:
Cv (T) y P = f(,T)
Este procedimiento de clculo es esencialmente el mismo que el presentado en el apartado 2.1, eliminando los
pasos 2,3,4,6, 7,8 y 12 e iniciando el clculo para unos valores de temperatura (T) y presin (P).
3. PROGRAMA DE COMPUTACION TATERMO
La palabra TATERMO es una conjuncin de las palabras Tablas Termodinmicas, nombre que usualmente se
da a las tablas de las propiedades termodinmicas de cualquier sustancia pura.
TATERMO es un programa de computacin concebido para elaborar las tablas de las propiedades
termodinmicas de treinta y tres sustancias puras; treinta y dos de esas sustancias pueden desarrollar sus tablas
termodinmicas para los estados de saturacin a temperatura de saturacin, cuatro sustancias puras pueden
construir sus tablas termodinmicas tambin en los estados de saturacin, pero para presiones de saturacin y
treinta y tres sustancias puras pueden conformar sus tablas termodinmicas en la fase gaseosa o estado de
sobrecalentamiento (Saavedra y Cristiano 1991; Saavedra y Cristiano 1993; Saavedra 1993).
TATERMO es un programa desarrollado para elaborar tablas termodinmicas en intervalos especficos de
temperatura y presin a incrementos deseados y a temperaturas de referencia preestablecidas. Este programa
puede ser empleado en cualquier microcomputador, PC, XT, AT, que trabaje con cualquier versin del Sistema
Operativo MS-DOS.
3.1.- Estructura del programa TATERMO.
El programa TATERMO lo constituyen cinco mdulos y un archivo de datos (Saavedra 1993). Ver Figura 2.
Estos mdulos son:

TATERMO.EXE, COMBUST.EXE, MISCELA.EXE, REFRIGE.EXE y BRT7IEFR.EXE, y el archivo de

datos es TABLAS.DAT.
Cada mdulo tiene una funcin especfica de importancia para el programa TATERMO. Las funciones de los
mdulos son las siguientes:
TATERMO.EXE. Este mdulo, diseado con el compilador BASIC PDS ver 7.1, tiene como objetivo
principal hacer la presentacin en pantalla del programa TATERMO exponiendo las diferentes opciones del
programa y direccionando las ejecuciones necesarias hacia los mdulos que procesan los datos que darn como
resultado las tablas termodinmicas.
COMBUST.EXE. Este mdulo se encarga de calcular las propiedades termodinmicas de las sustancias puras
siguientes: Butano, Etano, Etileno, Heptano, Hexano, Isobutano, Metano, Octano, Pentano, Propano y
Propileno, de acuerdo a las instrucciones y opciones seleccionadas en el mdulo TATERMO.EXE.

Figura 2.- Diagrama de Bloques del Programa TATERMO.

MISCELA.EXE. Este mdulo de clculo se encarga de determinar las propiedades termodinmicas de las
sustancias puras siguientes: Agua, Aire Seco, Alcohol Proplico, Argn, Dixido de Carbono, Hidrgeno,
Mercurio, Nitrgeno y Oxgeno, tambin dadas las instrucciones y opciones del mdulo de presentacin
TATERMO.EXE.
REFRIGE.EXE. Este mdulo ejecuta los clculo necesarios para obtener las propiedades termodinmicas de
las sustancias puras siguientes: Amonaco y Refrigerantes 11, 114, 12, 13, 14, 22, 23, 500, 502, 503 y C318, de
acuerdo a las opciones sealadas en el mdulo TATERMO.EXE.
BRT7IEFR.EXE. Este mdulo se encarga de enlazar el mdulo de presentacin TATERMO.EXE con
cualquier mdulo de clculo, COMBUST.EXE, MISCELA.EXE, y REFRIGE.EXE, y viceversa.
El archivo de datos TABLAS.DAT contiene todos los intervalos de temperatura y presin que por defecto
puede ejecutar, si as se desea, el mdulo TATERMO.EXE para elaborar las tablas de las propiedades
termodinmicas de las treinta y tres sustancias puras. No se recomienda el clculo de las propiedades para
valores fuera de estos intervalos.

Los mdulos de clculo, COMBUST.EXE, MISCELA.EXE, y REFRIGE.EXE estn constituidos por un

conjunto de subrutinas de nombres iguales a las sustancias puras que contienen, las cuales son las que
determinan realmente las propiedades termodinmicas en los estados de referencia, saturacin y
sobrecalentamiento, segn el caso que seleccione el mdulo de presentacin TATERMO.EXE. Adicionalmente
estos mdulos de clculo contienen una subrutina de nombre DIALOGO, la cual se encarga de informar al
usuario acerca de los errores que se comenten en la entrada de datos y/o ejecucin de los mdulos de clculo.

3.1.1.- Manejo bsico del programa TATERMO.

TATERMO es un programa de manejo fcil. Inclusive puede ser operado con la accin de un ratn (mouse). Si
el usuario no posee ratn, las teclas que debe tomar en consideracin para manejar este programa son las
siguientes (Saavedra 1993):
[BARRA ESPACIADORA]: Para confirmar la seleccin <OK> que aparece en algunas ventanas.
[ENTER]: Para activar la opcin seleccionada.
[ALT]: Para activar el men de barras.
[TAB]: Para moverse en una ventana activa.
[SHIFT] + [TAB]: Para moverse en sentido contrario en una ventana activa.
[ESC]: Para cancelar la opcin seleccionada o alguna ventana activa.
Y las teclas direccionadoras del cursor.
En el caso de la pantalla que se presenta en la Figura 3, que aparece cuando se instala el programa TATERMO,
puede cambiarse el monitor pulsando la tecla [TAB], y una vez seleccionado se pulsa [ENTER] o [BARRA
ESPACIADORA] y aparecer la pantalla de presentacin del programa TATERMO que se muestra en la Figura
4, la cual desaparecer en pocos segundos o, si se desea cancelarla, se pulsa [BARRA ESPACIADORA] y
aparecer el men de barras que se encuentra en la Figura 5, la cual puede activarse con la tecla [ALT].

Figura 3.- Seleccin del Monitor.

3.2.- Presentacin del programa TATERMO
El men principal del programa TATERMO, compuesto por un men de barras, el cual se muestra en la
pantalla de la Figura 5, contiene 5 opciones, Tablas, Impresin, Opciones, Presentacin y Salida (Saavedra
1993), y las funciones de stas son:

a. Tablas. Esta opcin presenta un submen que da al usuario la posibilidad de elaborar tablas de
propiedades termodinmicas para temperatura de saturacin, presin de saturacin o estado de
sobrecalentamiento. La Figura 6 muestra el submen de esta opcin.
b. Impresin. La opcin impresin brinda al usuario la posibilidad de obtener la impresin por impresora o
pantalla de los valores de las propiedades termodinmicas de una tabla elaborada y grabada con
anterioridad.

Figura 4.- Presentacin del Programa TATERMO.

Figura 6.- Submen Tablas.

c. Opciones. Da al usuario la posibilidad de cambiar de tipo de monitor y las lectoras donde se encuentran
los mdulos de presentacin, clculo y enlace y los archivos de datos.
d. Presentacin. Facilita al usuario el aprendizaje para manejar con cierta facilidad el programa
TATERMO.
e. Salida. Con esta opcin se finaliza la tarea que estuvo realizando TATERMO. Tiene a su vez tres
opciones para salir hacia el sistema operativo.

4. CONCLUSIONES Y RECOMENDACIONES.
a) Se desarroll un mtodo de ensayo y error (Saavedra 1993) necesario para calcular el volumen
especfico de una sustancia pura a una presin y temperatura dadas. Tomando como variable la
presin para cerrar el proceso de clculo por tanteo. El error relativo que se comete en el mismo
es menor a 0.05%.
b) El tiempo que toma el programa en calcular el valor de las propiedades termodinmicas para un
estado definido es menor a 5 segundos (Saavedra 1993), incluyendo presentaciones en pantalla y
grabacin en disco. Siempre y cuando se emplee un microcomputador AT-286 con un
coprocesador matemtico.
c) Se desarrollaron y codificaron en BASIC los modelos matemticos que calculan los valores de
las propiedades termodinmicas de las sustancias puras siguientes (Saavedra 1993): Aire Seco,
Alcohol Proplico, Agua, Amonaco, Argn, Butano, Dixido de Carbono, Etano, Etileno,
Heptano, Hexano, Hidrgeno, Iso-Butano, Iso-Pentano, Mercurio, Metano, Nitrgeno, Octano,
Oxgeno, Pentano, Propano, Propileno y Refrigerantes 11, 12, 13, 14, 22, 23, 114, C-318, 500,
502 y 503.
d) Los resultados de las propiedades termodinmicas calculadas en este trabajo fueron comparados
con aquellos obtenidos por Reynolds (Reynolds 1979). Los errores relativos que se cometieron
en los clculos realizados en todos los casos son menores al 2%, presentando el mayor error para
la entropa especfica. En cualquier caso el error relativo que se comete en el proceso de clculo
de tanteo para la presin es menor a 0.05%.
e) Para este ao deben incluirse en TATERMO los refrigerantes sustitutos desarrollados por
DuPont Chemicals. SUVA MP39, SUVA MP66, SUVA 134a, SUVA HP80 y SUVA HP81.
5. REFERENCIAS BIBLIOGRAFICAS
1.- AGUDELO, ALONSO y DURAN, JOSE. 1985. "Tablas Termodinmicas - Parte II". Escuela de
Ingeniera Mecnica, Universidad de Carabobo, Venezuela. Profesor Asesor: Francisco Saavedra.
2.- COLAN, LUIS. 1985. "Tablas Termodinmicas - Parte I". Escuela de Ingeniera Mecnica, Venezuela.
Profesor Asesor: Francisco Saavedra.
3.- CRISTIANO, FERNANDO; SAAVEDRA, FRANCISCO; ROMAN, HECTOR; ALCALA, RICHARD;
y DE POOL, CLARA. 1984. "Tablas de Propiedades Termodinmicas: Agua, Amonaco y Refrigerantes
11, 12 y 22". 1era y 2da Edicin, Facultad de Ingeniera, Universidad de Carabobo, Venezuela.
4.- CRISTIANO, FERNANDO y SAAVEDRA FRANCISCO. 1988. "Tablas Termodinmicas I". Escuela
de Ingeniera Mecnica, Universidad de Carabobo, Venezuela.
5.- RAMIS, RAFAEL y VERA, AQUILES. 1992. "Elaboracin de Diagramas Termodinmicos a travs del
Uso de Computadoras". Escuela de Ingeniera Mecnica, Universidad de Carabobo, Venezuela. Profesor
Asesor: Francisco Saavedra.
6.- REYNOLDS, WILLIAM C. 1979. "Thermodynamics Properties in SI. Graphs, Tables and
Computational Equations for Forty Sustances". Departament of Mechanical Engineering, Stanfort
University, California, E.E.U.U.
7.- SAAVEDRA, FRANCISCO y CRISTIANO, FERNANDO. 1987. "Tablas Termodinmicas II". Escuela
de Ingeniera Mecnica, Universidad de Carabobo, Venezuela.
8.- SAAVEDRA, FRANCISCO y CRISTIANO, FERNANDO. 1991. "Tablas Termodinmicas 29
Sustancias". Escuela de Ingeniera Mecnica, Universidad de Carabobo, Venezuela.
9.- SAAVEDRA, FRANCISCO y CRISTIANO, FERNANDO. 1993. "Tablas Termodinmicas de 33
Sustancias Puras". Escuela de Ingeniera Mecnica, Universidad de Carabobo, Venezuela.

10.- SAAVEDRA, FRANCISCO. 1993. "Aplicacin de la Computacin al Clculo de las Propiedades

Termodinmicas Mediante el Uso de los Compiladores Quick BASIC ver 4.5 y BASIC PDS ver 7.1".
Escuela de Ingeniera Mecnica, Universidad de Carabobo, Venezuela.

6.- Nomenclatura

NOTA: Para solicitar un disquete con el programa TATERMO y el manual de Usuario puede
comunicarse con el autor a la direccin y/o telfono que aparece en el encabezado.