Analisis Termodinamico de Un Motor de Ciclo Stirling PDF
Analisis Termodinamico de Un Motor de Ciclo Stirling PDF
Analisis Termodinamico de Un Motor de Ciclo Stirling PDF
STITU
UTO POLIT
P TCNIICO N
NACIO
ONAL
ESCUELA SUPERIO
S OR DE INGENIE
I ERA
MECNICA Y ELCTTRICA
SECCIN DE
E ESTUDIIOS DE P
POSGRAD
DO E INV
VESTIGA
ACIN
UNIDA
AD PROFESIONAL A
ADOLFO LPEZ
L MA
ATEOS
AN
LISIS
S TERM
MODINMICOO DE
UN MOTOOR DE CICLO
C STIRL
LING
TESIS
PARA
A OBTENER
R EL GRAD
DO DE
MAES
STRO E
EN CIEN
NCIAS
CON LA ESPE
ECIALIDAD
D EN
INGE
ENIERA
A MECN
NICA
PRESENTA
I
ING. JOS ARNULFO
A MIRANDA
A MORALE
ES
B
BAJO LA DIR
RECCIN DE:
D
DR. SA
AMUEL ALCCNTARA MONTES
M
Mayoo 2008
AGRADECIMIENTOS
Doy gracias a Dios porque me ha permitido realizar una etapa importante en mi vida.
A mi madre Irma Morales Acosta y a mi padre Fabin Miranda Guzmn, a quienes dedico este trabajo
como una pequea muestra de las cosas tan valiosas que en mi han fomentado
A mi hermano Jos del Carmen Miranda Morales por su desinteresado apoyo tanto econmico como moral.
A mis tos Rudy, Ignacio y ngel Gabriel Morales Acosta por haber confiado en m al momento de iniciar
Un agradecimiento especial al Dr. Samuel Alcntara Montes por haber dirigido este trabajo de tesis y por su
Finalmente se agradece a todas aquellas personas que me apoyaron incondicionalmente para la terminacin
de este trabajo.
INDICE GENERAL
INDICE GENERAL
DESCRIPCIN PAG.
RESUMEN I
ABSTRACT III
INTRODUCCIN V
CAPITULO I
ANTECEDENTES HISTRICOS 3
1.1 Nacimiento del Motor Stirling 3
1.1.1 Patente de1816 4
1.1.2 Mejoras Realizadas por Stirling al diseo original 7
1.1.3 Aplicaciones al motor Stirling a mediados del siglo
9
XVIII
1.1.4 Motor Robinson 10
1.1.5 Motor Lehman 10
1.1.6 Motor Rider 11
1.1.7 Controversia sobre el primer diseo del Motor Rider 11
1.2 La Philips y el Motor Stirling 12
1.2.1 Bautizando al Motor Stirling 13
1.2.2 Mejoras realizadas a los diseos anteriores 13
1.2.3 El Stirling como refrigerador 14
1.2.4 Contribucin de la Philips en la Inundacin de Holanda
15
en 1953
1.2.5 Mecanismo Rmbico 15
1.2.6 Visita de Henry Ford II a los laboratorios Philips 16
1.3 El Motor Stirling en la actualidad 16
1.3.1 El Motor Stirling de pistn libre tipo Beale 17
1.3.2 El motor Stirling en la industria automotriz 17
1.4 Tipos de configuraciones 18
1.4.1 Motor tipo ALFA 19
1.4.2 Motor tipo BETA 19
1.4.3 Motor tipo GAMMA 20
INDICE GENERAL
CAPITULO II
PRIMERA LEY DE LA TERMODINMICA:
23
VOLMENES DE CONTROL
2.1 Anlisis Termodinmico de volmenes de control 23
2.2 Principio de la conservacin de la masa 24
2.2.1 Relaciones de flujo de masa y volumen 25
2.3 Principio de la conservacin de la energa 27
2.3.1 Trabajo de flujo 27
2.3.2 Energa total de un fluido que fluye 28
2.4 Proceso de flujo permanente 29
2.4.1 Conservacin de la masa 30
2.4.2 Conservacin de la energa 31
2.5 Proceso de flujo no permanente 32
2.5.1 Conservacin de la masa 33
2.5.2 Conservacin de la energa 34
2.6 Caso Especial: proceso de flujo uniforme 36
2.7 Ciclo Stirling
37
CAPITULO III
ANLISIS ISOTRMICO 42
3.1 Anlisis Isotrmico 42
3.2 Solucin de las ecuaciones 46
3.3 Anlisis de energa del modelo isotrmico 51
CAPITULO IV
ANLISIS ADIABTICO 58
4.1 Anlisis Adiabtico 58
4.2 Desarrollo del sistema de ecuaciones 62
CAPITULO V
ANLISIS GENERAL 74
5.1 Anlisis General 74
5.2 Desarrollo del sistema de ecuaciones 77
INDICE GENERAL
CONCLUSIONES Y RECOMENDACIONES
Conclusiones 91
Recomendaciones 93
APNDICE A 96
BIBLIOGRAFIA 103
Resumen
RESUMEN
El estudio del modelo isotrmico, se llevo a cabo con la finalidad de proponer una
expresin que permita calcular la presin del fluido de trabajo en funcin de las variaciones de
volumen en los espacios de compresin y expansin, obteniendo as, el trabajo realizado por el
motor; posteriormente se aplic la ley de la conservacin de la energa al modelo isotrmico
para estudiar la transferencia de calor en los diferentes componentes del motor y analizar de
esta manera el comportamiento terico de la mquina.
I
Resumen
manera, se proponen las expresiones para las variaciones cclicas de presin y temperatura del
fluido de trabajo.
Las ecuaciones obtenidas de los dos primeros estudios realizados (modelo isotrmico y
adiabtico) se evalan con el programa matemtica, los resultados que se obtuvieron
demostraron que el estudio isotrmico no cumpli con el objetivo trazado debido a que por la
naturaleza misma del modelo isotrmico, los espacios de compresin y expansin se
mantienen a las temperaturas respectivas del enfriador y del calentador, llegando de esta
manera a una evidente contradiccin.
Aunque los resultados tericos del modelo adiabtico fueron satisfactorios es necesario
compararlos con datos experimentales, por lo tanto, es necesario construir una mquina
Stirling, y con datos reales hacer mejoras a los anlisis propuestos.
II
Abstract
ABSTRACT
Presently work was carried out a thermodynamic analysis, to a motor Stirling of simple
action, with the purpose of deriving a system of equations that describe the main design
parameters like they are the pressure and the temperature
This analysis was developed, using a configuration type alpha and a widespread
diagram of the mechanism of the machine Stirling, in the configuration type alpha was carried
out two studies one isothermic and another adiabatic, in the widespread diagram you carries
out a general analysis whose results can be applied to any configuration type (configuration
alpha, configuration beta and configuration gamma) of motors Stirling of simple action.
The isothermic study you carries out with the purpose of to propose an expression that
allows to calculate the pressure of the work fluid in function of the variations of volume in the
compression spaces and expansion, and to obtain this way, the work carried out by the motor,
later on, the law of the conservation was applied from the energy to the isothermic pattern to
study the transfer of heat in the different components of the motor and to analyze this way the
behavior of theoretical the machine.
The study adiabatic was carried out with the purpose of to obtain a more complete
system of equations and to obtain improvements to the results obtained with the isothermic
pattern, in the development of the system of equations it was applied the state equation and the
conservation equation from the energy to each component of the motor, the obtained results
are related applying the equation of continuity through the whole system.
Finally a general analysis is made, in which, the transference is heat as well as the
variation of the superficial area of the cylinders of compression and expansion, of equal way,
pressure variations of and temperature of the fluid of the work set out the expressions for the
cyclical
The obtained equations of both first made studies (isothermal and adiabatic model) are
evaluated with the mathematical program, the results that were obtained demonstrated that the
III
Abstract
isothermal study did not fulfill the objective drawn up because by the same nature of the
isothermal model, the spaces of compression and expansion stay to the respective
temperatures of the cooler and the heater, arriving this way to an evident contradiction.
Although the theoretical results of the model adiabatic was satisfactory, is necessary to
evaluate them with experimental data, therefore, is necessary to build a machine Stirling, and
with real data to do improvements to the analysis proposed.
IV
INTRODUCCION
INTRODUCCION
Desde su invencin, el motor de ciclo Stirling como ms tarde fue llamado, tuvo un
gran xito, aunque si bien es cierto no ayud a disminuir en gran medida el derroche
energtico en su tiempo, en la actualidad se ha colocado como uno de los candidatos ms
fuertes.
V
INTRODUCCION
La gran aceptacin que tuvo el motor inventado por Robert Stirling por sus
contemporneos fue opacada con la aparicin del motor de combustin interna, ya que este es
mucho ms potente. Con los nuevos avances tcnicos, y ante la creciente contaminacin
ambiental as como la cada vez mayor escasez de energticos, centros de investigacin,
universidades y empresas importantes como la Philips Co., la General Motors., entre otras, han
retomado en las ltimas dcadas el estudio del motor Stirling y sus aplicaciones.
El nacimiento del motor Stirling no tuvo una fundamentacin terica que sustentara su
principio de funcionamiento, mucho menos indicadores que marcaran sus parmetros de
diseo, la casi nula actividad que existe en nuestro pas con respecto a este tipo de motores se
realiza de igual manera a prueba y error.
Objetivo General
Objetivos Particulares
VI
INTRODUCCION
Justificacin
Alcance
Desarrollo de tres anlisis trmicos (isotrmico, adiabtico y general) que nos permite
estudiar a detalle el comportamiento trmico de un motor de ciclo Stirling
Aportaciones
En el capitulo dos se hace una revisin de los principios termodinmicos que se aplican
a sistemas que incluyen flujo de masa a travs de sus fronteras, es decir, volmenes de control,
de igual manera, se analizan los procesos de flujo permanente (modelo para dispositivos de
VII
INTRODUCCION
VIII
CAPTULO
I
ANTECEDENTES HISTORICOS
CAPITULO I ANTECEDENTES HISTORICOS
ANTECEDENTES HISTORICOS
Las mquinas trmicas regenerativas del tipo Stirling, se pueden utilizar tanto para
producir energa mecnica como enfriamiento, se pueden utilizar en bombas de calor para
calefaccin elctrica pasando el calor de una temperatura inferior a otra superior. En Holanda
la compaa Philips ha investigado intensivamente los motores Stirling pequeos con
rendimientos superiores al 30 por ciento, los aspectos mas relevantes de la historia de dicha
mquina se presentan a continuacin.
El motor de ciclo Stirling, o motor de aire caliente* como tambin fue llamado tiene
una interesante y larga historia. Fue inventado en 1816 por el reverendo Robert Stirling,
ministro escocs, a los de 26 aos de edad.
Las capacidades mecnicas de Robert Stirling eran naturales ya que proviene de una
familia de prominentes ingenieros, su abuelo Michael Stirling, invent la primer mquina de
paleta rotatoria en 1756, y su hermano James era un ingeniero civil muy conocido. Cuatro de
los hijos de Robert Stirling tambin fueron ingenieros muy notables. (1)
*
El motor de aire caliente es un trmino amplio y probablemente obsoleto que incluye diferentes tipos de
motores de ciclo abierto, tales como el Ericcson, y de ciclo cerrado como el Stirling.
3
CAPITULO I ANTECEDENTES HISTORICOS
Otros motores de aire se haban ideado antes del diseo de Stirling, como el realizado
por Sir George Caley en 1807, tambin se tiene conocimiento de dispositivos llamados
motores de aire desde 1699, lo ms importante del diseo de Stirling es el empleo del
regenerador, el cual adquiri mucha importancia ya que precedi al nacimiento de la
termodinmica y a los escritos de Sadi Carnot por 40 aos.
La patente que registr Robert Stirling en 1816 se titulaba Mejoras para disminuir
el consumo de combustible, y en particular, una mquina capaz de ser aplicada al
movimiento de la maquinaria basada en un principio enteramente nuevo. La patente se
ocupa en gran parte del principio del regenerador, el cual es un medio para reutilizar el calor,
que de otra manera sera perdido en cualquier tipo de proceso. Se dice que Stirling propuso el
trmino de economizador, pero el trmino regenerador que introdujo John Ericsson fue
ms popular. (1)
Stirling no slo describi la construccin y el uso del regenerador, tambin previ sus
primeras aplicaciones: fabricacin de vidrio, procesos de elaboracin de cerveza, destilacin,
entre otros. Adems, incluy la descripcin del primer motor de ciclo cerrado que funcionaba
con aire caliente, as como el principio de operacin, en el cual el volumen del gas es
4
CAPITULO I ANTECEDENTES HISTORICOS
5
CAPITULO I ANTECEDENTES HISTORICOS
En unas de las primeras conferencias cientficas sobre motores de aire caliente dada
por Michael Faraday, el famoso cientfico tuvo que confesar que no poda explicar como
funcionaban estos motores en su totalidad. 30 aos tuvieron que pasar despus de la invencin
original para que llegaran las primeras explicaciones tericas.
A mediados del siglo XIX gran parte de la comunidad cientfica sostena la errnea
teora calrica del calor. Bajo esta teora, el motor de Stirling era considerado como una
posible mquina de movimiento perpetuo, y gran parte de las primeras investigaciones fueron
hechas con la intencin de perfeccionar tal mquina. Incluso, de acuerdo a la termodinmica
moderna, el ciclo Stirling es tericamente tan eficiente como podra ser la mquina de Carnot,
en cualquier rango de temperatura dado.
Con la patente de Stirling se construy un motor para el bombeo de agua en una mina
de Ayrshire, Escocia, en 1818. Dicho motor tena las mismas dimensiones de la patente
original: una altura del cilindro de casi 10 pies de altura, un dimetro de 2 pies
aproximadamente y una potencia de salida estimada en 2 caballos de fuerza.
6
CAPITULO I ANTECEDENTES HISTORICOS
Aunque el primer motor tuvo aparentemente cierto xito, Stirling reconoci que no
tena la suficiente potencia para satisfacer las demandas energticas de su poca. En 1824
James su hermano menor le sugiri elevar la presin en la mquina como una forma de
obtener mayor potencia. (1)
De esta manera, entre 1824 y 1840, los hermanos Stirling trabajaron en forma conjunta
llevando a cabo diversas innovaciones al diseo original. Entre las principales modificaciones
se encuentra el uso de un nivel de presin, el cual increment considerablemente la eficiencia.
Otra innovacin de la cual se duda, es la construccin de una mquina gemela, es decir, una
maquina de pistn de doble accin y cilindros desplazadores separados, como el mostrado en
la Fig. 1.2. (3)
7
CAPITULO I ANTECEDENTES HISTORICOS
En este modelo el cilindro era ms largo para evitar que se calentaran las tapas y haba
un dispositivo donde se coloc el regenerador y el enfriador fuera del cilindro. (3)
Todos estos modelos no tuvieron gran xito econmico a pesar de que la potencia y
eficiencia eran mayores que en las mquinas de vapor de aquella poca. La razn principal fue
el poco conocimiento en materiales resistentes al calor, lo que trajo como consecuencia que
fcilmente se quemaran los cilindros. Con el descubrimiento de Bessemer de un proceso de
produccin de acero de calidad en serie, las mquinas de vapor llegaron a ser de mayor
8
CAPITULO I ANTECEDENTES HISTORICOS
potencia y mucho ms seguras de funcionar, por tal motivo, la mayora de las mquinas fueron
rediseadas a mquinas de vapor. (3)
A partir de 1860 y hasta la primera guerra mundial, una gran cantidad de motores que
operaban con el principio de Stirling, adems del diseo de Ericcson, fueron construidos en los
Estados Unidos, Inglaterra y gran parte de Europa. La potencia de salida se extendi de una
fraccin de caballo de fuerza a varios caballos de fuerza. Estos motores fueron utilizados en el
bombeo de agua, irrigacin, uso domstico, bombeo de aire a los rganos de las iglesias,
imprentas, mquinas de costura. En algunos talleres eran la nica fuente de energa. Para esta
poca, estos motores haban mejorado y eran ms confiables al ser utilizados en la industria
metalrgica. No se requera de un ingeniero para que operara estas mquinas, por lo tanto las
mquinas de vapor se fueron haciendo obsoletas. Al construir los motores omitieron el
regenerador y no los presurizaban, por consiguiente la potencia y eficiencia eran bajas con
respecto al tamao, sin embargo fueron muy tiles, tanto por su seguridad como por su
potencia, que aunque baja satisfaca muy bien la demanda energtica de la poca,
constituyendo una vez ms una alternativa para las mquinas de vapor.
Duplex Vacuum, construy motores que fueron empleados en los tpicos vagones de
los carritos de palomitas, sus diseos fueron quizs los ms estticos de todos los fabricantes
americanos de motores Stirling de ese periodo.
9
CAPITULO I ANTECEDENTES HISTORICOS
Inglaterra, tambin fabric muchos motores Stirling, entre sus fabricantes mas
importantes figuran Bailey, Gadner y Hayward Tyler, entre otros.
Uno de los motores Stirling ms populares y hermosos de los 60s fue el construido por
Lehmann (fue a este motor al que se le hizo el primer anlisis termodinmico). Esta mquina
era muy similar a la version patentada por Stirling en 1816, y de ella algunas personalidades
de la poca escribieron:
10
CAPITULO I ANTECEDENTES HISTORICOS
Los catlogos de manufactura de aquella poca, algunos de los cuales se han vuelto a
reimprimir, muestran que estos motores eran simples, confiables y podan funcionar todo el
da sin supervisin.
Un tipo de motor Stirling muy usado fue el Rider, desarrollado por A. K. Rider, de
Philadelphia, en 1876. Este motor usaba el ciclo Stirling, pero utilizaba dos pistones de trabajo
para mover el fluido en vez de emplear un pistn de trabajo y un desplazador. Uno de los
pistones fue colocado en el cilindro caliente y el otro en el cilindro de enfriamiento, ambos
cilindros estaban conectados por medio de otro cilindro hueco en el cual haba un regenerador
construido con finas placas de hierro espaciadas y cerradas. Los dos pistones estaban en un
ngulo de 90; es decir, cuando un pistn estaba en la parte superior del cilindro, el otro estaba
a la mitad del recorrido.
Muchos de los motores Rider, fueron fabricados y utilizados para el bombeo de agua.
Se cree que el motor Rider fue la base para la invencin del motor de doble pistn propuesto
por Rinia en 1947, lo cual fue un gran acontecimiento en su momento, pero como sucede con
frecuencia en la historia de la creatividad, aunque fueron muy novedosos cada uno en su
poca, ni el motor Rider ni el motor Rinia eran nuevos, ya que ambos fueron reinventos del
motor de aire caliente patentado por Robert Stirling.
La primera persona que us esta ingeniosa y sutil versin del ciclo Stirling con dos
pistones de trabajo, al parecer, fue Charles Franchot en 1853. Su motor de dos pistones de
trabajo se describi en un fascinante artculo de Babcock, titulado sustituto para el vapor,
impreso en 1885.
11
CAPITULO I ANTECEDENTES HISTORICOS
ninguna otra forma de motor de aire ofrece tantas ventajas, pero tambin tiene sus
dificultades muy peculiares. Si estas dificultades llegan a ser superadas, es muy probable que
pueda convertirse en el motor de aire del futuro
A finales del siglo XIX, los motores elctricos y de gasolina fueron desarrollados; y
debido a su mayor potencia para un tamao dado, comenzaron a sustituir a los motores
Stirling, posteriormente la mquina Stirling slo se utilizaba para determinados trabajos como
el bombeo de agua en las casas o accionar ventiladores en las comunidades rurales. (1)
12
CAPITULO I ANTECEDENTES HISTORICOS
Cuando los ingenieros de la Philips, comenzaron sus investigaciones sobre este tipo de
motor, observaron que el nombre era muy particular y a la vez muy general, el motor
simplemente era conocido como motor de aire, en ese entonces muchos eran los motores
conocidos bajo este nombre, pero que no seguan el principio de Stirling, ah radica su
generalidad y su particularidad en el hecho que el motor funcionaba con aire caliente, de ah el
nombre, pero si se quita el aire como fluido de trabajo y se emplea otro tipo de gas ( por
ejemplo helio o hidrgeno), dejara de ser un motor de aire caliente, as que se vieron en la
necesidad de buscarle un nombre mas adecuado a tal dispositivo. El Dr. Rolf J. Meijer
encargado en ese momento del proyecto propuso darle al motor el nombre de su inventor. Al
no haber una mejor propuesta, la mquina se llam motor Stirling, en honor a su inventor. (2)
13
CAPITULO I ANTECEDENTES HISTORICOS
Philips describi dos tipos de motores Stirling en sus primeras publicaciones. Uno fue
el de pistn-desplazador, que era una copia del diseo original de la patente de 1816, y el otro
fue el de doble pistn diseado por H. Rinia, quien fu unos de los directores del programa
Stirling de la Philips; los dos tipos de motores usaban aire como fluido de trabajo.
Muchas pruebas fueron realizadas en este periodo, pero una prueba en particular que
entusiasm mucho a los ingenieros de la Philips, fue la que hicieron en los canales holandeses
al colocar un motor Stirling de dos y medio caballo de fuerza en un bote de remo; el motor no
era ms grande que una botella de medio galn, por lo cual lograron fcilmente cubrirla con
una caja de cartn. El bote motorizado recorri 50 millas alrededor de los canales de Holanda,
a una velocidad de 80 millas por hora aproximadamente. Las personas que caminaban a esa
hora por la orilla de los canales quedaron totalmente sorprendidas al ver el bote deslizarse
silenciosamente a travs del agua, sin ningn medio de propulsin visible. (1)
Pero los ingenieros de la Philips no slo desarrollaron el Stirling como motor, tambien
hicieron investigaciones del Stirling como refrigerador. Tericamente, el ciclo Stirling es un
ciclo reversible, es decir, si se sustrae calor de un espacio refrigerado (cuerpo fro), la maquina
opera en un ciclo de refrigeracin. Basndose en este principio, la Philips comenz a realizar
pruebas con un prototipo conectado a un motor elctrico. Al encender el motor, el Stirling
giraba en sentido de las manecillas del reloj. En ese momento la mquina actuaba como un
refrigerador capaz de enfriar a -320 0F (-195.56 0C) y licuar el aire sin la necesidad de
precomprensin. Si el motor elctrico era detenido, inmediatamente el Stirling cambiaba su
sentido al contrario de las manecillas del reloj y entonces se converta en una bomba de calor
elevando la temperatura de los cilindros hasta llegar al rojo vivo.
14
CAPITULO I ANTECEDENTES HISTORICOS
El diseo de doble pistn propuesto por Rinia, que era la principal esperanza para los
motores de alto poder, presentaba serios problemas en el sellado, por lo tanto se vieron en la
necesidad de inventar un nuevo tipo de mecanismo. Fue entonces cuando se ide el
15
CAPITULO I ANTECEDENTES HISTORICOS
mecanismo rmbico que permite el arreglo pistn-desplazador sin la necesidad de una caja
de cigeal presurizada. Por otro lado, las dimensiones y los pesos de las partes reciprocantes
del mecanismo pueden seleccionarse para un buen balance aun en motores monocilndricos.
Con la introduccin del mecanismo rmbico, la philips cambi el aire como fluido de
trabajo, por el helio y el hidrgeno. Estos gases proporcionan ms energa y una mejor
eficiencia, reduciendo prdidas del flujo y mejorando la transferencia de calor, aunque al usar
estos gases tenan el problema de sellado. (3)
En 1948 Henry Ford II visit los laboratorios Philips y qued muy impresionado ante
las posibilidades del nuevo motor. Para 1950 la Ford Motors Co. ya tena varios prototipos
cuyo funcionamiento se mostr a especialistas europeos y americanos; dicha presentacin
consisti en un motor Stirling de 40 caballos y de un cilindro, una eficiencia trmica del 38%,
superior a cualquier motor de gasolina y tan alta como la del mejor motor Diesel (en el motor
de automvil comn la eficiencia trmica es de 20%).
Al igual que sucedi con los laboratorios Philips y sin razn aparente, cuando se
trabajaba en un motor de 150 caballos diseado por computadora, la gerencia de la General
16
CAPITULO I ANTECEDENTES HISTORICOS
Motors, en forma sorpresiva, decidi suspender el programa Stirling, dejando este ltimo
motor a medio construir.
Para entonces ya haba varios grupos en el mundo que desarrollaban nuevos prototipos
de uso industrial y la Philips haba otorgado licencias a empresas como la MAN de Alemania,
la United Stirling de Suecia, entre otras, que prosiguieron las investigaciones.
En 1972, 24 aos despus de la visita de Henry Ford II a Holanda, la Ford Motors Co.
anuncia un convenio con la Philips para producir un motor de automvil, el cual, por una
parte, resuelve el problema de las emisiones indeseables de escape y, por otra, los cada vez
ms altos costos de la gasolina, ya que este motor puede trabajar con otros tipos de
combustibles. En 1976 se revel al pblico el primer motor Stirling instalado en un Ford
torino. Mientras la Ford perfecciona el motor Stirling de automvil, la United Stirling de
Suecia ofrece un motor V-4, ideal para casas-remolque y sistemas de clima artificial. Por otras
parte, William Beale funda la Sun Power Inc., que se prepara para producir electricidad con el
calor del sol y su motor de pistn libre, en tanto que la NASA (National Aeronautics and
Space Administration) ensaya con los motores desarrollados por la General Motors y, en
17
CAPITULO I ANTECEDENTES HISTORICOS
Harwell, Inglaterra demuestran cmo funciona una boya martima basada en el Stirling de
pistn libre atmico. (4)
18
CAPITULO I ANTECEDENTES HISTORICOS
Est compuesto por dos pistones en cilindros separados, los cuales son conectados en
serie por una cmara de enfriamiento, un regenerador y una cmara de calentamiento; un
cilindro comprende la zona de compresin a temperatura baja y el otro la zona de expansin a
temperatura alta.
Consiste de un solo cilindro que contiene dos pistones colocados uno frente al otro, un
pistn es conocido como pistn de trabajo y el otro como pistn desplazador, ambos
pistones estn separados por el regenerador. La funcin del pistn desplazador es pasar el
fluido de trabajo de la zona de expansin a la zona de compresin. El motor original de Robert
Stirling en su patente de 1816 muestra una configuracin beta.
19
CAPITULO I ANTECEDENTES HISTORICOS
20
CAPTULO
II
PRIMERA LEY DE LA TERMODINMICA:
VOLMENES DE CONTROL
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
En el presente capitulo se estudian sistemas que incluyen flujo de masa a travs de sus
fronteras. La ecuacin de la conservacin de la energa para un volumen de control general es
bastante compleja. Por esta razn, el anlisis de energa de los volmenes de control se har en
dos etapas.
Un volumen de control tiene tamao y forma fija, pero tambin incluye una frontera
mvil, sin embargo la mayor parte de los volmenes de control, tienen fronteras fijas y, por
ello, no implican ningn trabajo de la frontera mvil. Un volumen de control tambin implica
interacciones de calor y trabajo igual que en un sistema cerrado, adems de la interaccin de
masa.
23
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
E = mc 2 2.1
Esta ecuacin indica que la masa de un sistema cambiar cuando cambie su energa .
No obstante, en todas las interacciones de energa, con excepcin de las reacciones nucleares,
el cambio en la masa es extremadamente pequeo e incluso no puede detectarse.
Donde:
24
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
La cantidad de masa que fluye en una seccin transversal por unidad de tiempo se
denomina relacin de flujo de masa y se denota por m , el punto sobre el smbolo se emplea
para indicar derivada con respecto al tiempo.
d m = Vn dA 2.3
Donde
Vn = Componente de velocidad normal a dA
m = Vn dA 2.4
A
25
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
m = V pro A 2.5
Donde
= Densidad (kg m 3 ) = 1
V pro = Velocidad promedio del fluido normal a A (m s )
( )
A = rea de la seccin transversal normal a la direccin del flujo m 2
El volumen del fluido que circula a travs de una seccin transversal por unidad de
tiempo se llama relacin de flujo de volumen V y esta dada por
(m s )
V = Vn dA = V pro A 3
2.6
A
V
m = V = 2.7
26
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
Q W = E 2.8
Es claro que si no hay masa que entre o salga del volumen de control, el segundo y
tercer trmino se anulan y la ecuacin se reduce a la expresin de la primera ley de la
termodinmica para un sistema cerrado. La ecuacin 2.9 se aplica a cualquier volumen de
control sujeto a cualquier proceso. La transferencia de calor hacia un volumen de control o de
un volumen de control, no debe confundirse con la energa transportada con la masa hacia
dentro o hacia fuera del volumen de control. La energa necesaria para hacer pasar al fluido
dentro o fuera de un volumen de control recibe el nombre de trabajo de flujo, o energa de
flujo; se considera como parte de la energa transportada con el fluido y se analiza a
continuacin.
27
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
A diferencia de los sistemas cerrados, los volmenes de control incluyen flujo de masa
a travs de sus fronteras, por lo que se necesita cierto trabajo para empujar la masa dentro o
fuera del volumen de control. Este trabajo se conoce como trabajo de flujo o energa de flujo,
y es necesario para mantener un flujo continuo a travs de un volumen de control.
El trabajo de flujo por unidad de masa se obtiene al dividir ambos lados de la ecuacin
entre la masa del elemento del fluido:
w flujo = P 2.11
La relacin del trabajo de flujo es la misma sin importar que el fluido se empuje hacia
dentro o hacia fuera del volumen de control.
La energa total de un sistema compresible simple esta formada por 3 partes: energa
interna, energa cintica y energa potencial, con base en una masa unitaria, la energa total se
expresa como
e = u + ec + ep 2.12
El fluido que entra o sale de un volumen de control posee una forma de energa
adicional: la energa de flujo P , como se indic antes. As, la energa total de un fluido que
fluye con base en una masa unitaria (denotada ) se convierte en
= P + e = P + (u + ec + ep ) 2.13
28
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
= h + ec + ep 2.14
Muchos dispositivos de ingeniera operan durante largos periodos bajo las mismas
condiciones y se clasifican como dispositivos de flujo permanente. Los procesos que implican
dispositivos de flujo permanente son representados por medio de un proceso idealizado,
denominado proceso de flujo permanente. Un proceso de este tipo se define como un proceso
durante el cual un fluido fluye permanentemente por un volumen de control, es decir, las
propiedades del fluido cambian de un punto a otro dentro del volumen de control, pero en
cualquier punto fijo permanecern iguales durante todo el proceso. Un proceso de este tipo se
caracteriza por lo siguiente:
1.- Ninguna propiedad (intensiva o extensiva) dentro del volumen de control cambia
con el tiempo.
2.- Ninguna propiedad cambia en las fronteras del volumen de control con el tiempo.
3.- Las interacciones de calor y trabajo entre un sistema de flujo permanente y sus
alrededores no cambia con el tiempo.
29
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
de flujo permanente mediante valores promediados en el tiempo para las propiedades y las
relaciones de flujo de calor a travs de las fronteras.
m en = m sal 2.15
La mayor parte de los dispositivos de ingeniera, son de una sola corriente (nicamente
una entrada y una salida). Para estos casos, el subndice 1 denota el estado de la entrada y el
subndice 2 el estado de la salida, por lo tanto, la ecuacin 2.15 se reduce, en sistemas de flujo
permanente de una sola corriente, a:
m1 = m 2 2.16
Donde
= Densidad
V = Velocidad de flujo promedio en la direccin del flujo
30
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
Q W = m sal sal m en en 2.18
Donde es la energa total del fluido que circula, incluido el trabajo de flujo, por
unidad de masa.
En sistemas de una sola corriente (una entrada, una salida) se eliminan las sumatorias
sobre las entradas y las salidas y los estados de entrada y salida, se denotan por medio de los
subndices 1 y 2, respectivamente. La relacin de flujo de masa por todo el volumen de control
permanece constante m1 = m 2 . La ecuacin de la conservacin de la energa para sistemas
de flujo permanente de una sola corriente se convierte en
V 2 V12
Q W = m h2 h1 + 2 + g ( z1 z 2 ) 2.19
2
Si se divide la ecuacin anterior entre m , se obtiene la relacin de la primera ley en
una base de masa unitaria como
31
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
V22 V12
q w = h2 h1 + + g ( z 2 z1 ) 2.20
2
Donde
Q
q=
(Transferencia de calor por unidad de masa)
m
W
w=
(Trabajo efectuado por unidad de masa)
m
q w = h 2.21
Muchos procesos de inters implican cambios dentro del volumen de control con el
tiempo. Dichos procesos se llaman procesos de flujo no permanente o de flujo transitorio.
Cuando se analiza un proceso de flujo no permanente, es importante seguir de cerca los
contenidos de masa y de energa del volumen de control, as como las interacciones de energa
a travs de la frontera.
32
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
sistema de flujo no permanente es similar a un sistema cerrado, excepto en que la masa dentro
de las fronteras del sistema no permanece constante durante un proceso.
Con frecuencia algunos trminos en la ecuacin 2.23 son cero. Por ejemplo, men = 0 si
no entra masa al VC durante el proceso, m sal = 0 si no sale masa del VC durante el proceso y
m1 = 0 si el VC est vaci inicialmente.
33
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
dmVC
m en m sal = dt
2.24
d
V n dA V n dA = dV 2.25
A en A sal dt V
Q W + en sal = E vc 2.26
34
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
Ven2
en = m
men
en = hen +
men
2
+ gz en men
2.28
o en forma de relacin
V2
en = m en hen + en + gz en 2.29
2
Al repetir esto en cada entrada y salida y sustituir en las ecuaciones 2.26 y 2.27, se
obtiene
35
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
Vsal2 Ven2
Q W = sal 2 + gz sal msal
h + en 2 + gz en men + EVC
h + 2.30
msal
y en forma de relacin
Vsal2 Ven2 dE
Q W = sal sal 2 + gz sal
m h + en en 2 + gz en + dtVC
m h + 2.31
En general los procesos de flujo no permanente son difciles de analizar debido a que
las integraciones en la ecuacin 2.30 no son fciles de efectuar. Sin embargo, algunos
procesos de flujo no permanente puede representarse razonablemente bien mediante otro
modelo simplificado: el proceso de flujo uniforme. Un proceso de flujo uniforme involucra las
siguientes idealizaciones que simplifican el anlisis:
2.- Las propiedades del fluido pueden diferir de una entrada o salida a otra aunque el
flujo del fluido en una entrada o salida sea uniforme y permanente, es decir, las propiedades
no cambian con el tiempo o la posicin sobre la seccin transversal de una entrada o salida. Si
cambian, son promediadas y tratadas como constantes para todo proceso.
36
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
Vsal2 Ven2
Q W = sal sal 2 + gz sal
m h + en en 2 + gz en + (m2 e2 m1e1 )VC
m h + 2.31
Se advierte que si no entra o sale masa al volumen de control o del volumen de control
(men = m sal = 0) , los primeros dos trminos del lado derecho de la relacin anterior se
eliminan y esta ecuacin se reduce a la relacin de la primera ley para sistemas cerrados.
Consideremos una mquina trmica que opera entre un depsito de alta temperatura a
TH y un depsito a baja temperatura TL . Para que el ciclo de la mquina de calor sea
totalmente reversible, la diferencia de temperatura entre el fluido de trabajo y la fuente (o
sumidero) de energa trmica nunca debe exceder una cantidad diferencial dT durante
cualquier proceso de transferencia de calor. Es decir, los procesos de adicin y rechazo de
calor durante el ciclo deben de suceder de modo isotrmico, uno a cierta temperatura TH y el
otro a una temperatura TL . Esto es precisamente lo que sucede en un ciclo de Carnot.
37
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
La figura 2.2 muestra el diagrama P-v del ciclo Stirling, el cual est integrado por
cuatro procesos totalmente reversible:
38
Captulo II Primera Ley de la Termodinmica: Volmenes de Control
39
CAPTULO
III
ANALISIS ISOTERMICO
CAPITULO III ANALISIS ISOTERMICO
ANALISIS ISOTERMICO
Conocer la presin del fluido de trabajo, es muy importante, ya que permite estudiar
cmo diferentes mecanismos de transmisin afectan la potencia de salida, por lo tanto, el
objetivo fundamental del presente capitulo es proponer una expresin que permita calcular la
presin del fluido en funcin de las variaciones de volumen en los espacios de compresin y
expansin.
El anlisis descrito por Schmidt fue el mtodo mas simple y til durante el desarrollo
de la mquina Stirling, se basa en la expansin y compresin isotrmica de un gas ideal, este
mtodo se ha convertido en el anlisis clsico y punto de partida para el estudio de los motores
de ciclo Stirling. Desafortunadamente en el anlisis de Schmidt se obtiene una solucin
cerrada lo cual imposibilita enormemente predecir el ciclo real, no obstante lo utilizamos
como punto de partida para dirigirnos en ltima instancia a un acercamiento ms realista.
Existe una gran variedad de diseos del motor Stirling, los cuales se han agrupado en 3
configuraciones, para realizar el anlisis se tomara como modelo la configuracin tipo alfa, la
cual se muestra en la Fig. 3.1
42
CAPITULO III ANALISIS ISOTERMICO
Estas suposiciones permiten generar una expresin simple para calcular la presin del
fluido de trabajo en funcin de las variaciones del volumen. Esta expresin se puede utilizar
para investigar cmo diversos mecanismos de transmisin afectan la potencia de salida.
Ahora bien, consideremos la Fig. 3.1 como un modelo en serie de cinco componentes,
cada componente estar representado de acuerdo a la siguiente nomenclatura:
1. Zona de Compresin c
2. Enfriador k
3. Regenerador r
4. Calentador h
43
CAPITULO III ANALISIS ISOTERMICO
5. Zona de Expansin e
1. Zona de Compresin (m c , Tc , Vc y p )
2. Enfriador (m k , Tk , V k y p )
3. Regenerador (mr , Tr , Vr y p )
4. Calentador (m h , Th , V h y p )
5. Zona de Expansin (m e , Te , V e y p )
El punto de partida del anlisis es que la masa total del fluido en la mquina es
constante, as:
M = mc + m k + m r + m h + me 3.1
44
CAPITULO III ANALISIS ISOTERMICO
m = pV RT 3.2
Factorizando se obtiene:
p Vc Vk Vr Vh Ve
M= + + + +
R Tk Tk Tr Th Th
3.4
MR
p =
Vc V k V r V h Ve 3.5
+ + + +
Tk Tk Tr Th Th
El trabajo hecho por el sistema sobre un ciclo completo esta dado respectivamente por
la integral cclica de p dV
dVe dVc
W = We + Wc = p dVc + p dVe = p + d 3.6
d d
Donde:
= Angulo de la manivela
45
CAPITULO III ANALISIS ISOTERMICO
MR
p=
( )
Presin
Th
Vc Vk Vr ln Tk Vh Ve
+
Tk Tk + + +
(Th Tk ) Th Th
dVe
Qe = We = p d Calor
d Transferido
dVc
Qc = Wc = p d
d
W = Wc + We Trabajo Efectuado
w
= Eficiencia
Qe
Tabla 3.1 Sistema de Ecuaciones
La solucin a las ecuaciones obtenidas del anlisis isotrmico para el caso especial de
las variaciones de volumen senoidal de los espacios de trabajo con respecto al ngulo de la
manivela ( ) , esta basada en el anlisis de Schmidt..
Vbe
Ve = Vme + (1 cos ) 3.7
2
Vc = Vmc +
Vbc
[1 cos( )] 3.8
2
46
CAPITULO III ANALISIS ISOTERMICO
Donde:
MR
p=
Vmc + 2 [1 cos( )] Vk Vr Vh Vme + 2 (1 cos )
Vbc Vbe
3.9
+ + + +
Tk Tk Tr Th Th
Simplificando obtenemos:
MR
p=
Vbc V V 3.10
a cos + be cos + bc sen sen
2Tk 2Th 2Tk
Donde:
47
CAPITULO III ANALISIS ISOTERMICO
Vbc
c sen = sen 3.11
2Tk
Vbe V
c cos = + bc cos 3.12
2Th 2Tk
Vbc
sen
1 2Tk
= tg 3.13
Vbe Vbc
+ cos
2Th 2Tk
c=
1 (Vbe )2 +
2(Vbe )(Vbc ) (Vbc )
+
2
3.14
2 Th Th Tk Tk
MR
p=
a (c cos cos + c sen sen )
3.15
Simplificando
MR
p=
a (1 bcos )
3.16
Donde:
48
CAPITULO III ANALISIS ISOTERMICO
c
a=
s
MR
p min =
a (1 b )
3.17
MR
p max =
a (1 + b )
3.18
2
1
p prom =
2 p d
0
3.19
2
MR d
p prom =
2 a 0 (1 bcos )
3.20
MR
p prom = 3.21
a 1 b2
49
CAPITULO III ANALISIS ISOTERMICO
2
dV
We = p e d 3.22
0
d
dVe V sen
= be 3.23
d 2
2
MR Vbe sen
We = a (1 bcos )
0
2
d
3.24
2
Vbe MR sen d
We =
2a 0 1 bcos( )
3.25
MR 1 b2 1
We = Vbe sen 3.26
a 1 b 2 b
1 b2 1
We = Vbe Pprom sen 3.27
b
sen( )
2
V MR
Wc = bc
2a 0 1 bcos( )
d 3.28
50
CAPITULO III ANALISIS ISOTERMICO
1 b2 1
Wc = Vbc Pprom sen( ) 3.29
b
Las ecuaciones 3.27 y 3.29 describen el trabajo realizado por el motor, en funcin de
las variaciones de volumen en las zonas de compresin y expansin respectivamente.
Ahora consideremos el modelo isotrmico desde el punto de vista del flujo de energa.
La transferencia de calor hacia y desde los alrededores ocurre a las temperaturas Tk y Th
respectivamente. Normalmente los motores Stirling tendrn intercambiadores de calor
separados del enfriador y el calentador, como se indica en la Fig. 3.1. Para estudiar la
transferencia de calor en estos espacios, es necesario considerar la ley de conservacin de la
energa para el fluido del trabajo.
Una clula representativa se muestra en la Fig. 3.2 la cual puede representar un espacio
de trabajo o un intercambiador de calor.
dW
dQ m, p, T , V
51
CAPITULO III ANALISIS ISOTERMICO
de energa a travs de las fronteras del sistema en forma de calor y trabajo, as como de la
cantidad de energa transportada hacia dentro y fuera del volumen de control mediante la
masa, durante el proceso. Por esta razn es importante seguir de cerca el contenido de energa
del volumen de control, as como las energas de las corrientes entrante y saliente.
Q W + en sal = E vc 3.30
O en forma diferencial
d Q d W + d en d sal = dE vc 3.31
Donde representa la energa total transportada por la masa dentro o fuera del
volumen de control a travs de la entrada o la salida durante el proceso. No es fcil determinar
la energa total transportada por la masa hacia dentro o hacia fuera del volumen de control
( en , sal ) , puesto que las propiedades de la masa en cada entrada o salida pueden cambiar
con el tiempo, as como en la seccin transversal. De modo que para determinar el transporte
de energa a travs de una abertura como un resultado del flujo de masa se deben considerar
masas diferenciales m suficientemente pequeas con propiedades uniformes y sumar sus
energas totales.
V2
= dm = h + + g z dm 3.32
2
O en forma diferencial
52
CAPITULO III ANALISIS ISOTERMICO
v2
d = dm h + + gz 3.33
2
Donde:
v2
La energa cintica por unidad de masa e c = y la energa potencial por unidad de
2
masa e p = g z se consideran despreciables
Obteniendo finalmente
d = dm(h ) 3.35
La energa total de un sistema compresible simple esta formado por 3 partes: energa
interna, cintica y potencial, los cambios en la energa cintica y potencial son despreciables
por lo tanto la energa total del sistema seria igual a la energa interna es decir:
E = U 3.35
En su forma diferencial
dE = dU 3.36
Donde:
h = c pT 3.38
U = mc v T 3.39
Derivando
53
CAPITULO III ANALISIS ISOTERMICO
dU = c v d (mT ) 3.40
Por lo tanto
dQ dW + c p T dm = c v Tdm 3.45
Factorizando
dQ dW + Tdm (c p c v ) = 0 3.46
dQ = dW RTdm 3.47
54
CAPITULO III ANALISIS ISOTERMICO
La transferencia de calor neta Q al fluido de trabajo sobre el ciclo esta dado por la
integral cclica de dQ . Sin embargo, la consideracin implcita de estado esttico en el ciclo
implica que el cambio en la masa (m ) del fluido de trabajo es cero para cada una de las clulas.
Por lo tanto aplicando la ecuacin de la conservacin de energa para cada una de las clulas
isotrmicas e integrando sobre el ciclo, obtenemos para los espacios de compresin y
expansin.
Qc = Wc
Qe = We
Qk = 0
Qh = 0
Qr = 0
Este resultado implica que todos los intercambiadores de calor en el motor ideal
Stirling son redundantes puesto que toda la trasferencia de calor externa ocurre a travs de los
lmites de los espacios de compresin y de expansin. Esta evidente paradoja es resultado
directo de la definicin del modelo isotrmico ideal en el cual los espacios de compresin y de
expansin se mantienen a las temperaturas respectivas del enfriador y del calentador.
Obviamente esto no puede ser correcto, puesto que las paredes del cilindro no se disean para
la transferencia de calor.
55
CAPTULO
IV
ANALISIS ADIABATICO
CAPITULO IV ANALISIS ADIABATICO
ANALISIS ADIABATICO
Debido a esta situacin es necesario considerar un anlisis alternativo, que nos lleve a
un acercamiento mas realista del funcionamiento del motor, para realizar tal anlisis nos
apoyaremos en el proceso adiabtico, ya que de acuerdo a las conclusiones del capitulo
anterior los espacios de trabajo (zonas de compresin y expansin) tienden a ser adiabticos
ms que isotrmicos, al igual que se hizo con la seccin anterior, en este capitulo se tomara la
configuracin tipo alfa mostrada en la Fig. 3.1 como base para realizar el anlisis adiabtico.
58
CAPITULO IV ANALISIS ADIABATICO
1. Zona de Compresin c
2. Enfriador k
3. Regenerador r
4. Calentador h
5. Zona de Expansin e
59
CAPITULO IV ANALISIS ADIABATICO
distribucin de temperaturas de la Fig. 4.1, la cual se llevara a cabo desde el punto de vista
adiabtico.
1. Zona de Compresin (m c , Tc , Vc y p )
2. Enfriador (m k , Tk , V k y p )
3. Regenerador (m r , Tr , Vr y p )
4. Calentador (m h , Th , V h y p )
5. Zona de Expansin (m e , Te , V e y p )
En base a la Fig. 4.1 se define la nomenclatura de las 4 interfaces entre cada uno de los
componentes o clulas de la configuracin tipo alfa:
60
CAPITULO IV ANALISIS ADIABATICO
condicionadas por la direccin del flujo, de esta manera, la temperatura Tck se define de
acuerdo a la siguiente expresin:
Tck = Tc H m ck + Tk 1 H m ck 4.1
H ( x a ) = 1 si x > a
4.2
H ( x a ) = 0 si x < a
H m ck = 1 si m ck > 0
4.3
H m ck = 0 si m ck < 0
The = Th H m he + Te 1 H m he 4.4
61
CAPITULO IV ANALISIS ADIABATICO
dW
dQ m, p, T , V
Donde cp y cv son las capacidades calorficas especficas del gas a presin y volumen
constante respectivamente, los trminos de energa cintica y potencial son despreciados.
62
CAPITULO IV ANALISIS ADIABATICO
Suponemos que el fluido de trabajo es ideal. Esto es una suposicin razonable para los
motores Stirling puesto que los procesos del fluido de trabajo se trasladan lejos del punto
crtico del gas. La ecuacin de estado para cada clula se presenta en su forma estndar y
diferencial como sigue
p V = mR T 4.6
En su forma diferencial
d p dV d m d T
+ = + 4.7
p V m T
El punto de partida del anlisis es que la masa total del gas en la mquina es constante,
as:
M = mc + mk + m r + m h + me 4.8
m = pV RT 4.9
Factorizando se obtiene
p Vc Vk Vr Vh Ve
M= + + + +
R Tc Tk Tr Th Te
4.11
63
CAPITULO IV ANALISIS ADIABATICO
MR
p=
Vc Vk Vr Vh Ve 4.12
+ + + +
Tc Tk Tr Th Te
Los volmenes y las temperaturas respectivas son constantes en cada uno de los
intercambiadores de calor (calentador, regenerador y enfriador), por lo tanto, la ecuacin 4.7
se reduce a:
d p dm
= 4.14
p m
Despejando dm
dm =
dp
(m) 4.15
p
d p V
dm = 4.16
R T
dp V dp Vr dp Vh
dmc + dme + k + + = 0 4.17
R Tk R Tr R Th
Factorizando
64
CAPITULO IV ANALISIS ADIABATICO
dp V V V
dmc + dme + k + r + h = 0 4.18
R Tk Tr Th
Es necesario eliminar dmc y dme en la ecuacin anterior para obtener una ecuacin
explcita para dp .
Por lo tanto
65
CAPITULO IV ANALISIS ADIABATICO
pV
dmc (c p Tck ) = p dVc + cv d c 4.24
R
pV
Derivando cv d c
R
Factorizando p dVc
R + cv cv
dmc (c pTck ) = p dVc + (Vc dp ) 4.27
R R
cp c
dmc (c pTck ) = p dVc + v (Vc dp ) 4.28
R R
Factorizando
66
CAPITULO IV ANALISIS ADIABATICO
c p p dVc + cv Vc dp
dmc (c pTck ) = 4.29
R
p dVc cv dp Vc
dmc = +
c p R Tck
4.30
R Tck
cp
Sustituyendo la relacin = se obtiene
cv
p dVc 1 dp Vc
dmc = + 4.31
R Tck R Tck
dp Vc 1
dmc = p dVc + 4.32
R Tck
dp Ve 1
dme = p dVe +
R The
4.33
dp Vc 1 dp Ve 1 dp Vk Vr Vh
p dVc + + p dVe + + + + =0 4.34
R Tck he R Tk Tr Th
R T
67
CAPITULO IV ANALISIS ADIABATICO
dV dVe
p c +
Tck The
dp = 4.35
Vc V V V V
+ k + r + h + e
Tck Tk Tr Vh The
dp dVc dmc
dTc = Tc +
mc
4.36
p Vc
dp dVe dme
dTe = Te + 4.37
p Ve me
dp
dQ + c p (dmen Ten dmsal Tsal ) = cv V 4.39
R
Para el enfriador
dp
dQk = cv Vk c p (dmck Tck dmkr Tkr ) 4.40
R
68
CAPITULO IV ANALISIS ADIABATICO
Para el regenerador
dp
dQr = cv Vr c p (dmkr Tkr dmrh Trh ) 4.41
R
Para el calentador
dp
dQh = cv Vh c p (dmrh Trh dmhe The ) 4.42
R
W = Wc + We 4.43
mc = pVc RT c
mk = pVk RT k
mr = pVr RT r Masas
mh = pVh RT h
me = pVe RT e
Tabla 4.1 Sistema de Ecuaciones Diferenciales y Algebraicas
69
CAPITULO IV ANALISIS ADIABATICO
MR
p=
Vc Vk Vr Vh Ve
+ + + +
Tc Tk Tr Th Te
dV dVe Presin
p c +
Tck The
dp =
Vc V V V V
+ k + r + h + e
Tck Tk Tr Vh The
dp Vc 1
dmc = p dVc +
R Tck
dp Ve 1
dme = p dVe +
R The
dmk =
dp
(mk ) Acumulacin de masa
p
dmr =
dp
(mr )
p
dmh =
dp
(mh )
p
dmck = dmc
H m ck = 1 si m ck 0 ; H m ck = 0 si m ck 0
Tck = Tc H m ck + Tk 1 H m ck
Condiciones de temperatura
H m he = 1 si m he 0 ; H m he = 0 si m he 0
The = Th H m he + Te 1 H m he
Continuacin Tabla 4.1 Sistemas de Ecuaciones Diferenciales y Algebraicas
70
CAPITULO IV ANALISIS ADIABATICO
dp dVc dmc
dTc = Tc +
p Vc mc
Condiciones de temperatura
dp dVe dme
dTe = Te +
p Ve me
dp
dQk = cv Vk c p (dmck Tck dmkr Tkr )
R
dp
dQr = cv Vr c p (dmkr Tkr dmrh Trh )
R
dp
dQh = cv Vh c p (dmrh Trh dmhe The )
R Energa
dWc = p dVc
dWe = p dVe
dW = dWc + dWe
W = Wc + We
Continuacin Tabla 4.1 Sistema de Ecuaciones Diferenciales y Algebraicas
La solucin del conjunto de ecuaciones mostradas en la tabla 4.1 solo ser posible de
forma numrica, esto debido a la naturaleza no lineal de las ecuaciones.
71
CAPTULO
V
ANLISIS GENERAL
CAPITULO V ANALISIS GENERAL
Los mtodos usados para calcular los principales parmetros de funcionamiento de las
primeras mquinas Stirling fueron muy precarios, operaciones a volumen constante fueron
supuestas, as mismo se tom cualquier ley conveniente (isotrmico, politrpico o adiabtico),
para las fases de compresin y expansin, incluso relaciones lineales de presin-volumen
fueron utilizadas de vez en cuando.
En 1871 Gustav Schmidt public un anlisis en el cual calcul el trabajo de salida por
medio de integracin, el anlisis termodinmico descrito por Schmidt se basa en un ciclo
regenerativo suponiendo fases isotrmicas, de un gas ideal, el cual implica transferencia de
calor infinito (operacin isotrmica) a travs de las paredes del cilindro a velocidad cero.
El mtodo propuesto por Schmidt fue el ms utilizado por los ingenieros de la poca,
convirtindose as en un mtodo clsico en el clculo de los parmetros de funcionamiento de
las mquinas Stirling.
74
CAPITULO V ANALISIS GENERAL
75
CAPITULO V ANALISIS GENERAL
Ahora bien, consideremos la Fig. 5.1 como un modelo en serie de cinco componentes,
cada componente estar representado de acuerdo a la siguiente nomenclatura:
1. Zona de Compresin c
2. Enfriador k
3. Regenerador r
4. Calentador h
76
CAPITULO V ANALISIS GENERAL
5. Zona de Expansin e
1. Zona de Compresin (m c , Tc , Vc y p )
2. Enfriador (m k , Tk , V k y p )
3. Regenerador (m r , Tr , Vr y p )
4. Calentador (m h , Th , V h y p )
5. Zona de Expansin (m e , Te , V e y p )
El primer paso del anlisis es derivar una expresin que relacione la mayor parte de la
temperatura instantnea t del fluido en alguno de los dos espacios de trabajo, con el calor
transferido a travs de las paredes del cilindro durante la operacin regular del motor.
Para derivar esta expresin se debe establecer un balance energtico de la cantidad (m)
del fluido de trabajo que esta contenido en un instante cualquiera en uno de los dos espacios
de trabajo, por lo tanto es necesario igualar la razn de flujo de energa de entrada en uno de
los dos espacios de trabajo con la razn de flujo de energa de salida en este mismo espacio de
trabajo.
Q = a (Ts T ) 5.1
77
CAPITULO V ANALISIS GENERAL
Donde:
Es necesario relacionar los dos trminos de entalpa de los casos anteriores en una sola
expresin, para esto utilizamos la funcin escaln:
H (x ) = 1 cuando x > 0
5.2
H (x ) = 0 cuando x < 0
78
CAPITULO V ANALISIS GENERAL
Ahora bien si dm < 0 habr perdida de energa = dm hsal , por otra parte, si dm > 0
tendremos ganacia de energa = dm hent , por lo tanto:
[ ]
dm hsal + mc p dTsal = a (Ts T ) + dm hsal + H (dm )c p (Tent Tsal ) 5.9
79
CAPITULO V ANALISIS GENERAL
(Tent Tsal ) a
dTsal = + dm H (dm ) 5.10
m c p
(Th Te ) e a e
dTe = + dme H (dme ) 5.11
me c p
(Tk Tc ) c a c
dTc = + dmc H (dmc ) 5.12
mc c p
Donde:
= ngulo de desfase entre el pistn de trabajo y el pistn desplazador
80
CAPITULO V ANALISIS GENERAL
pVe
me = 5.15
RTe
Ve =
Vbarrido
(1 cos ) 5.16
2
me =
pVbarrido
(1 cos ) 5.17
2 RTe
mc =
pVbarrido
[1 cos( )] 5.18
2 RTc
Para el regenerador
pVr
mr = 5.19
RTr
Para el calentador
81
CAPITULO V ANALISIS GENERAL
pVh
mh = 5.20
RTh
Para el enfriador
pVk
mk = 5.21
RTk
M = mc + m k + m r + m h + me 5.22
log me = log
pVbarrido
(1 cos ) 5.23
2 RTe
82
CAPITULO V ANALISIS GENERAL
Despejando dTe
dp sen dme
dTe = + Te 5.26
p 1 - cos m e
dp sen dme (T Te ) e ae
+ Te = h + dme H (dme ) 5.27
p 1 - cos me me c p
Dividiendo entre Te
dp sen dme (T Te ) e ae
+ = h + dme H (dme ) 5.28
p 1 - cos me me Te c p
sen dme T 1 e ae
+ dme H (dme )
dp
+ = h 5.29
p 1 - cos me me Te me c p
Factorizando
sen T a dme Th a
1 + H (dme )
dp
+ h e e = 1 e e 5.31
p 1 - cos me Te c p me Te me c p
83
CAPITULO V ANALISIS GENERAL
sen Th e ae dme
Th
dp
+ = 1 + H (dm ) 1
p 1 - cos pV (1 cos ) me
e
barrido
pV (1 cos )
me barrido c p 5.32
2Rme
2Rme
e ae
me c p
Simplificando
Definimos los siguientes valores adimensionales de masa y presin tanto para la zona
de expansin como de compresin, as mismo, se define la presin para los intercambiadores
de calor
84
CAPITULO V ANALISIS GENERAL
me mc
e = c = 5.36
M M
mk + mr + mh
= 5.37
M
e r L c rL
e = c = 5.38
M cp M cp
Con las notaciones de masa y presin podemos establecer las siguientes ecuaciones
e + c + = 1 5.39
d sen 2 e d e 2 e e (1 cos )
+ = 1 + H (d e ) 1 5.40
1 cos e (1 cos ) e
d sen( ) 2c d c 2 c c (1 cos( ))
+ = 1 + H (d c ) 1 5.41
1 cos( ) c (1 cos( )) c
sen d e 2 e e (1 cos )
d + 2 e = 1 + H (d e ) 1 5.42
1 cos e (1 cos ) e
Simplificando
85
CAPITULO V ANALISIS GENERAL
d e e (1 cos ) sen
d e H (d e )
2
d = 2 e 5.43
e 1 cos e e 1 cos
= 1 e c 5.44
Derivando 5.44
d = d e d c 5.45
(1 e c ) d e (1 e c )
d e H (d e )
2
d e d c = 2 e
e 1 cos e
5.46
(1 e c ) e (1 cos ) (1 e c )sen
e (1 cos )
Factorizando
(1 e c ) (1 e c )
+ H (d e )
2
d c d e 1 + =
e 1 cos e
5.47
(1 e c )(1 cos ) sen (1 e c )
e 2
e 1 cos
86
CAPITULO V ANALISIS GENERAL
(1 e c ) (1 e c )
+ H (d c )
2
d e d c 1 + =
c 1 cos( ) c
5.48
(1 e c )(1 cos( )) sen( )(1 e c )
c 2
c 1 cos( )
Para obtener una solucin de las ecuaciones diferenciales 5.47 y 5.48, estas deben ser
reducidas algebraicamente. Las derivadas de las variables desconocidas e y c con respecto
al tiempo pueden ser escritas explcitamente en trminos de funciones compuestas de e y
c , que no implican otras derivadas excepto el valor conocido de d y las funciones
H (d e ) o H (d c ) . Estas frmulas finales son:
c c e
d e = 5.49
e c 1
e e c
d c = 5.50
e c 1
Donde:
(1 e c ) (1 e c )
+ H (d e )
2
e = 1 + 5.53
e 1 cos e
(1 e c ) (1 e c )
+ H (d c )
2
c = 1 + 5.54
c 1 cos( ) c
87
CAPITULO V ANALISIS GENERAL
Las ecuaciones 5.49 y 5.50 describen completamente las variaciones cclicas de masa,
temperatura y presin.
88
CONCLUSIONES
Y
RECOMENDACIONES
CONCLUSIONES Y RECOMENDACIONES
CONCLUSIONES
Por ultimo el tercer anlisis se efectu a un diagrama que se denomino generalizado del
mecanismo de un motor de ciclo Stirling, debido a que los resultados obtenidos de este anlisis
pueden ser aplicados a cualquier tipo de configuracin de motores Stirling de accin simple.
3. Las ecuaciones 3.40, 3.41, y 3.42, muestran que el calor neto transferido es
realizado por los intercambiadores de calor (calentador y enfriador), por lo
91
CONCLUSIONES Y RECOMENDACIONES
tanto, el trabajo realizado por la mquina es la suma de los trabajos en las zonas
de compresin y expansin, por otra parte, se obtuvo un sistema de ecuaciones
diferenciales no lineales, la naturaleza no lineal de estas ecuaciones, se debe al
hecho de que la distribucin de temperaturas en las zonas de trabajo no es
constante y estn condicionadas por la direccin del flujo.
Se considera que el presente trabajo puede contribuir en gran medida al estudio terico
del motor de ciclo Striling, tambin puede servir como base para desarrollar programas de
computo y poder analizar las variaciones del fluido de trabajo as como la presin y la
temperatura, de una manera ms didctica.
92
CONCLUSIONES Y RECOMENDACIONES
RECOMENDACIONES
93
APENDICE
A
APNDICE A
APNDICE A
MR
La ecuacin p = es una funcin par y puede ser representada en una
a (1 b cos )
serie de fourier en un intervalo completo p < x < p la cual corresponde a la serie de los
cosenos.
a0 n
P ( ) = + anCos A.1
2 n =1 q
Donde:
a0 = p( ) d
1
q0
n
p
an = p( )cos d
1
q0 q
dVe V Sen
= be A.2
d 2
2
dV
We = p e d A.3
0
d
2
a
n Vbe sen
We = 0 +
2
a n cos
q 2
d
A.4
0 n =1
96
APNDICE A
2 2
Vbe a0 n
We =
2 0 2
sen d + a0 n =1
n cos
q
sen d
A.5
2 2
Vbe a0 n
a cos ( ) sen d
2 0 2
We = sen d + A.6
n
0 n =1
a0 2 2
n cos n ( ) sen d
Vbe
We =
2
sen d + a A.7
2 0 n =1 0
2 2 2
Vbe a0
We =
2 2 sen d + a1 cos( ) sen d +
0 0
n=2
a n cos n( ) sen d
0
A.8
2 2
La integral sen d = 0 y la an cos n( ) sen d = 0 , por lo tanto la
n=2
0 0
2
a1 cos( ) sen d
Vbe
2 0
We = A.9
2 2
cos( ) sen d =
0
(cos cos + sen sen ) sen d
0
A.10
97
APNDICE A
2 2 2
2 2 2
cos( ) sen d
0
= cos cos sen d + sen sen sen d
0 0
A.12
2
Integrando cos cos sen d obtenemos:
0
2
cos cos sen d = cos sen 2 ( ) 2
0
A.13
0
2
( )
cos cos sen d = cos sen 2 2 cos sen 2 0 ( ) A.14
0
2
cos cos sen d = 0 A.15
0
1 cos 2
Utilizando la identidad trigonometrica sen 2 = , para simplificar la integral
2 2
2
sen sen sen d obtenemos:
0
2 2 2
1 cos 2 d cos 2
sen = sen 0 2 sen 0 2 d A.16
0
2 2
2 2 2
1 cos 2 sen sen
sen = (sen2 ) A.17
0
2 2 2 0 4 0
98
APNDICE A
2
1 cos 2 sen
sen = (2 ) sen (0) sen sen2(2 ) sen sen2(0) A.18
0
2 2 2 2 4 4
2
sen sen 2 d = sen A.19
0
2
Por lo tanto la integral de cos( ) sen d es:
0
Vbe
We = a1 sen A.21
2
2
n
p( ) cos
1
a1 = d A.22
q 0
q
p ( ) =
MR
A.23
a (1 b cos )
q = n =1 A.24
99
APNDICE A
2
1 MR cos
a1 =
a (1 b cos ) d
0
A.25
2
MR cos
a1 =
a (1 b cos ) d
0
A.26
2
cos
Integrando (1 b cos ) d
0
1 z2 2 dz
cos = ; d = A.27
1+ z 2
1 + z2
1 z 2 2 dz
2
cos
2
1 + z 2 1 + z 2
(1 b cos )
d =
1 z2
A.28
0 0
1 b
2
1+ z
2
cos
2
(
2 1 z 2 dz )
0
(1 b cos )
d = (
0
)(
1 + z 2 1 + z 2 b + bz 2 ) A.29
2
2
(1 z )dz 2
1+b
0
(1 + z ) z + 11 + bb
2 2 A.30
2
2
A B
1 + b 0 z 2 + 1 b
+ dz A.31
z 2 + 1
1+ b
A=
1
; B=
(1 + b )
A.32
b b
100
APNDICE A
2
2
1 (1 + b ) dz
1 + b 0 2 1 b b z 2 + 1
( ) A.33
b z +
1 + b
2 1 2
dz (1 + b ) dz
2
b 0 z 2 + 1
A.34
1 + b b 0 z 2 + 1 b
1+ b
2
2 1 dz (1 + b)
2
dz
1 + b b 0 b 0 z 2 + 12
2
A.35
1 b
z +
2
1+ b
2
2 1 1
tg 1
z
(1 + b ) tg 1 z
( ) 2
1 + b b 1 b 0
A.36
1 b b
1 + b 1 + b
0
2
2
2 1 1 + b 1 1+ b
(1 + b ) tg 1 tg
tg tg
b 2 0
1 b A.37
1 + b b 1 b 2
0
101
APNDICE A
1 1+ b 1+ b
tg 1 1 + b tg 2 tg 1 tg 0
1 b 2 1 b
2 b 1 b
A.38
1 + b (1 + b ) 2
1
b tg tg 2 tg (tg 0)
1
2 1 1+ b
(1 + b )
A.39
1 + b b 1 b b
2 1 2
A.40
b 1 b 2
b
2 2
A.41
b 1 b2 b
2 1
1 A.42
b 1 b 2
Finalmente obtenemos
2
cos 2 1
1
(1 b cos ) d =
b 2
1 b
A.43
0
MR 2
1 1
a1 = A.44
a b 2
1 b
2MR 1
a1 = 1 A.45
a b 2
1 b
102
APNDICE A
Vbe 2 MR 1
sen
We = 1 A.46
2 a b 2
1 b
MR 1 b2 1
We = Vbe sen A.47
a 1 b2 b
MR
Como Pprom = , obtenemos finalmente:
a 1 b2
1 b2 1
We = Vbe Pprom sen A.48
b
103
BIBLIOGRAFIA
BIBLIOGRAFIA
BIBLIOGRAFIA
1. Andy Ross, 1977. Stirling Cycle Engines. Second Edition. Published by Solar Engines
Phoenix.
2. www.tamarugo.cec.uchile.cl
7. www.strilingengines.org.uk
8. www.ent.ohiov.edu
9. www.robotiker.com
10. Finkelstein, Ted Generalized Thermodynamic Analysis of Stirling Engines, New York,
N.Y. Society of Automotive Engineeres, 1960
106