Conceptos Bsicos para Diseo de motor Stirling con baja diferencia de temperatura

Prof. Roberto Romn L. Departamento de Ingeniera Mecnica Universidad de Chile

1 Introduccin y Objetivos:
En este documento se entregan algunos conceptos bsicos de diseo y metodologa elemental de clculo de un motor Stirling con bajo gradiente trmico. Adicionalmente entregamos algunas indicaciones bsicas referentes al diseo del motor. Para efectos de simplicidad, hemos dividido el documento en las siguientes partes: Ciclo Stirling terico y real. Elementos bsicos del motor. Clculo termodinmico elemental. Aspectos constructivo s.

2 Ciclo Stirling Terico y Real:

El ciclo Stirling terico est compuesto por dos isotermas y dos evoluciones a volumen constante (ver figura 1). En las cuatro evoluciones hay intercambio de calor. Pero la magnitud del calor absorbido en el calentamiento a volumen constante es igual (pero de signo contrario) a la magnitud del calor cedido en el enfriamiento a volumen constante.

Si este calor se almacena en forma intermedia en un regenerador, la teora demuestra que es posible alcanzar el rendimiento de Carnot. La funcin de este regenerador es almacenar (en forma reversible) el calor cedido al exterior en el enfriamiento a volumen constante (evolucin 4-1) y volverlo a entregar al ciclo en el calentamiento a volumen constante (evolucin 3-4). Esto es posible si el regenerador est compuesto de un material poroso capaz de almacenar y ceder calor, que permita el paso del fluido de trabajo y que tenga conductividad trmica axial despreciable. Realmente es posible dimensionar regeneradores bastante buenos. En el caso del ciclo real obviamente el rendimiento es inferior al terico. Adems las evoluciones que se siguen no son exactamente las tericas. En la figura 2 vemos un ciclo real de un motor Stirling de aire caliente con y sin regenerador de calor. Este diagrama p-V ha sido obtenido de un pequeo motor real (ver referencias al final del documento).

3 Elementos Bsicos del Motor:

El motor Stirling es un motor de ciclo cerrado y fuente de calor externa. Se caracteriza fundamentalmente por tener un lado fro y un lado caliente. El fluido de trabajo se trasvasija desde el lado fro al lado caliente (y viceversa) por medio de un desplazador. Cuando el fluido se calienta su presin aumenta, al enfriarse su presin disminuye. Esta fuerza se transmite al exterior por medio de un pistn o membrana de trabajo. Debe existir un mecanismo tal que: Mueva el desplazador de tal forma que el fluido de trabajo pase de la zona fra a la caliente (y al revs). Transmita el esfuerzo del pistn o membrana al exterior. Sincronice el movimiento de pistn y desplazador. 2

El movimiento y sincronismo se logra por medio de bielas y manivelas adecuadas. En particular, la biela que mueve el desplazador est desfasada en aproximadamente 90 con respecto a la biela que es accionada por el pistn (o membrana). Recordemos: El desplazador sirve para mover el fluido del lado fro al caliente y viceversa. Por lo tanto no entrega ningn trabajo til al exterior. El pistn o membrana es el elemento que entrega trabajo til. Este tambin acciona (indirectamente) al desplazador. El movimiento global se asegura por bielas y manivelas que sincronizan pistn con desplazador. Adems suele haber un volante que sirve para almacenar energa cintica y hacer todo el movimiento ms suave y armnico.

4 Clculo Termodinmico Elemental:

4.1 Hiptesis Bsicas: En los prximos puntos plantearemos un clculo simplificado del motor Stirling. Se aplicarn los siguientes conceptos e hiptesis bsicas: El gas en el interior del motor se comporta como gas perfecto. Los volmenes muertos son cero. De existir regenerador, este se considera perfecto. La irreversibilidad termodinmica asociada a la transferencia de calor se tomar en cuenta por un coeficiente emprico.

Para los motores que se van a construir en el Seminario se trata de que operen con el menor gradiente de temperatura posible. Adems su presin media interna va a ser igual a la presin ambiente. Definiremos las siguientes variables fundamentales: V: el volumen medio interno. Vmax: volumen mximo del motor. Vmin: volumen mnimo del motor Tmax: temperatura mxima del fluido detrabajo. Tmin: temperatura mnima del fluido de trabajo.

El objetivo del clculo termodinmico es determinar la presin mxima, mnima y media en el motor a partir de las hiptesis de temperaturas. 3

4.2 Secuencia de Clculo: A partir de la geometra del motor, debemos determinar el volumen medio, mximo y mnimo. Recuerden que el volumen mximo ocurre con mximo desplazamiento del pistn de trabajo y el mnimo en la situacin opuesta. Durante el ciclo de trabajo, se pasa por el mximo y mnimo una vez y por los volmenes medios dos veces. Hay algunos aspectos que hay que tomar en cuenta: El volumen mnimo tiende a ocurrir hacia la temperatura mnima. Mientras ms chica sea la diferencia de volumen, ms grande tiende a ser la diferencia de presiones. La mxima diferencia de presiones se produce para un V de 0. La mxima potencia del motor se tiende a producir a m/m la mitad de la mxima velocidad. Para estos motores a aire, la mxima velocidad est entre los 200 y 300 rpm.

Vamos a ejemplificar con un caso concreto: El motor experimental ST tendr un dimetro de 15 cm, una altura de 2 cm y efectuaremos el clculo con aire e hidrgeno como fluido de trabajo. El desplazador ocupa la mitad del volumen interno. Dejaremos como variable el volumen desplazado. La secuencia de clculo a seguir ser la siguiente: Primero calcularemos la masa de fluido que encierra el motor. Luego, para un diferencial de temperatura fijado, calcularemos la presin mxima y mnima para un V de 0. Conociendo la condicin anterior, probaremos diferentes V hasta escoger uno adecuado. Con el dato anterior, definiremos dimetro y carrera del pistn (o membrana). Definida la carrera, comenzamos a definir la geometra del conjunto bielamanivela.

4.2.1 Masa de fluido: El volumen de fluido dentro del cilindro es de: V = r2h/2 4

Se considera h/2 para tomar en cuenta el volumen del desplazador. En nuestro caso r=7,5 cm y h/2 = 1 cm, de donde: V = 7,521/2 = 88,4 cm3 Ahora debemos calcular la masa del fluido de trabajo. Supongamos una temperatura media de 25C = 298K. Para ello aplicamos la ecuacin de estado de gases perfectos:

p V = m R 'T p V m= R'T 1 10 5 88,4 10 6 m= 287 298 m = 0,1033 [ gra mos ]

4.2.2 Presin mxima y mnima para diferencia de volumen cero: Vamos a suponer que la temperatura mxima ser de 50C y la mnima de 0C (para efectos del clculo). Entonces, si la variacin de volumen es de 0, entonces podemos calcular las presiones mximas y mnimas:

p max = p min p max p max

m R 'Tmax V m R 'Tmin = V 1,03 10 4 287 323 = 88,4 10 6 = 108.389 [ Pa ] 1,03 10 4 287 273 88,4 10 6 = 91.611 [ Pa ]

p min = p min

Es decir, tenemos una presin mxima de 1,08 bar y una mnima de 0,92 bar. El diferencial mximo disponible sera de 0,16 bar, lo que corresponde a 16.778 Pa, es decir 1678 mm de columna de agua. Esto equivale a casi 168 g/cm2 de diferencia. 4.2.3 Con diferencia de volumen:

Ahora supongamos que el pistn de trabajo tenga 2 cm de dimetro y 1 cm de carrera. El nuevo volumen medio ocurrir a la mitad de la carrera. La presin mxima tender a ocurrir a volumen mnimo y temperatura mxima y la presin mnima a volumen mximo y temperatura mnima. Replanteemos las ecuaciones: El volumen medio de fluido dentro del cilindro es de: V = r2h/2 + 2L/2

Con dimetro interno del cilindro y L la carrera.

En nuestro caso r=7,5 cm y h/2 = 1 cm; = 2 cm y L = 1 cm de donde: VM = 7,521/2 + 121/2= 89,93 cm3 VMax = 7,521/2 + 121= 91,50 cm3 VMin = 7,521/2 + 120= 88,4 cm3

Ahora debemos calcular la masa del fluido de trabajo. Supongamos una temperatura media de 25C = 298K. Para ello aplicamos la ecuacin de estado de gases perfectos:

p VM = m R 'T p VM R 'T 1 10 5 89,93 10 6 m= 287 298 m = 0,1051 [ gra mos ] m=

4.2.4 Presiones mximas y mnimas para diferencia de volumen real: Vamos a suponer que la temperatura mxima ser de 50C y la mnima de 0C (para efectos del clculo). Entonces, con la hiptesis que la mxima presin ocurre a volumen mnimo y temperatura mxima y la presin mnima a mximo volumen y temperatura mnima se puede calcular las presiones mximas y mnimas:

p max = p min = p max = p max p min p min

m R 'Tmax Vmin m R 'Tmin Vmax

1,05 10 4 287 323 88,4 10 6 = 110.316 [ Pa ]

1,05 10 4 287 273 = 91,5 10 6 = 90.038 [ Pa ]

El diferencial mximo disponible sera de 20.278 Pa, es decir 2028 mm de columna de agua. Esto equivale a casi algo ms de 202 g/cm2 de diferencia. Naturalmente los valores reales sern bastante menores. Esto se debe a que el ciclo real no sigue exactamente el terico y adems el roce entre cilindro y pistn. Pero los valores obtenidos nos sirven de referencia de diseo. Traduzcamos estos valores a fuerzas a esperar en el pistn de trabajo. La fuerza mxima ser hacia afuera y estar dada por:

Fmax = S ( p max p amb ) Fmax = ( p max p amb ) 2 Fmax Fmax 2 = (110.316 100.000) 1 10 4 2 = 3,24 [ N ]
2 2

La fuerza mnima ser hacia adentro y estar dada por:

Fmin = S ( p amb p min ) Fmax = ( p amb p min ) 2

2 2

2 Fmax = (90.038 100.000) 1 10 4 2 Fmax = 3,13 [ N ]

4.2.5 Hiptesis adicionales: Los anteriores nos dan los clculos bsicos. Para este motor hipottico es razonable suponer que dispondremos de un diferencial de hasta 3 N con respecto a la presin ambiente para trabajar con l (positivo y negativo). En un ciclo, es razonable suponer 2 N positivo (como media) y 2 N negativo. Como la carrera es de 1 cm, el trabajo realizado en un ciclo sera: 2 Ncm en carrera ascendente 2 Ncm en carrera descendente Esto nos da 4 Ncm por ciclo. Si suponemos que el motor funciona a 100 rpm, implica una potencia de:

100 0,01 [W ] 60 Pot = 66,7 [ mW ] Pot = 4

Esta potencia es susceptible de ser aumentada si el gradiente de temperatura aumenta, si la presin media aumenta o si las revoluciones aumentan. Tambin puede mejorar si las propiedades termodinmicas del fluido mejoran. El mismo clculo lo podemos realizar para diferencias de temperatura menores, en vista de nuestro objetivo de diseo. Es importante tomar en cuenta de que si uno utiliza una fuente de calor a temperatura Tc, el fluido en el interior del motor se calentar a T max = Tc - T1. Anlogamente, al enfriar el fluido de trabajo, realmente en el interior este se enfriar a: Tmin = Tf - T2 Esta prdida de temperatura toma en cuenta el hecho de que existe una resistencia trmica a superar tanto en el calentamiento como enfriamiento del fluido. Es decir, realmente se calienta menos de lo terico y se enfra menos. 4.2.6 Geometra bsica del motor: Con los clculos elementales recin enunciados, uno debe pasar ahora a definir la geometra bsica del motor. Esto implica:

Definir excentricidad de bielas para permitir movimiento correcto del desplazador y pistn. Distancia entre ejes y cara superior del motor, de forma de que no exista interferencia entre componentes. Ubicacin de componentes principales.

5 Aspectos Constructivos:
Se podran escribir muchas pginas sobre aspectos constructivos a considerar. Hay sin embargo, algunos puntos bsicos: Evitar que el desplazador toque cara superior o inferior del motor. Reducir volmenes muertos todo lo posible. Mantener el roce al mnimo posible. Equilibrar componentes mviles: no sera raro que fuerza de gravedad sobre desplazador sea similar a la fuerza que origina el pistn, por lo tanto si no est equilibrado, costar que el motor parta.

Lo mejor que tienen estos motores es que son muy tolerantes y sencillos. Es decir, aunque no estn fabricados de forma ptima, igual tienden a funcionar.

BUENA SUERTE!