GENERACIN DE ENERGA ELCTRICA CON TURBINAS WELLS
Alfredo Santa* Marco Lozano** Alvaro E. Pinilla***
Abstract
The present article, briefly, describes the theoretical design procedure of Wells turbines, typically utilized in Pneumatic Wave Energy Conversion Systems (PWECS). The principle behind the PWECS consists of extracting the pneumatic energy caused by an oscillating water column by the entrance of water waves in enclosed spaces. The paper establishes the design and experimental testing of a small permanent magnet generator suited for the wells turbine under development. The adequate matching between the Wells turbine and the generator has been conceived as a portable 50w d.c. electricity generation device; to be used, either, to powering telecommunications systems in military operations or to power navigation aids in rivers or at sea
Palabras Claves
Turbina Wells, generacin de energa, aerodinmica, turbo-maquinaria, Energa de las Olas.
NOTACIN
Vio"
A rea proyectada por el disco de la turbina B Nmero de aspas Longitud de cuerda del perfil C Coeficiente de sustentacin C^ Coeficiente de arrastre CD D Fuerza de arrastre Dt Dimetro del disco de la turbina h Relacin dimetros del cubo y del disco L Fuerza de sustentacin P Presin esttica Q Flujo r radio de la turbina, elemento variable R Radio de la turbina T Torque T* Nmero adimensional de torque T/ p .122.R5 U* U/52R Ut Velocidad rotacional inducida U Velocidad de la corriente de aire Velocidad efectiva VEFEC
APo
APo*
a
No
rl p
6
S2
Velocidad geomtrica Potencia Cada de presin a travs de la turbina (P1-P2) Nmero adimensional de presin APo/p .S2z.R2 ngulo de ataque ngulo de velocidad efectiva Eficiencia aerodinmica de la turbina Densidad del aire bajo condiciones estndar Solidez de la turbina Velocidad rotacional angular (SQR/U) Velocidad especfica de rotacin
SUBNDICES
1,2
Condiciones inicial y final
1. INTRODUCCIN
La turbina Wells es el dispositivo aerodinmico utilizado en la actualidad para las instalaciones
Ingeniero Mecnico Universidad EAFIT, Colombia. Ingeniero Mecnico Uniandes, Colombia. Profesor Titular Departamento de Ingeniera Mecnica Uniandes, Bogot-Colombia.
56
UNIVERSIDAD DE LOS ANDES
de extraccin de Energa de las Olas, o ms bien la energa cintica neumtica generada por una columna oscilante de agua. Esta turbina utiliza el flujo bi-direccional del aire, comprimido por una columna oscilante de agua; y en consecuencia su diseo implica que opere en un solo sentido de giro independiente del flujo de aire que incida sobre ella. Para conservar el sentido de giro de la turbina, sin importar la incidencia del aire, el rotor se disea con perfiles aerodinmicos simtricos dispuestos radialmente. [Figura No. 1]. La turbina Wells al igual que otro tipo de turbinas, extrae la energa cintica disponible del aire, sin cambios bruscos en su velocidad de rotacin que impliquen grandes aceleraciones y fuerzas aerodinmicas. Este dispositivo Wells se acopla a un pequeo generador elctrico de imanes permanentes para la generacin de energa elctrica a 12 Voltios y una capacidad nominal de 50 vatios. Las aplicaciones de este pequeo dispositivo pueden ir desde sistemas porttiles para operaciones militares de telefona o para la iluminacin de boyas marinas, entre otras.
respecto al plano de giro del rotor vara desde el radio del cubo hasta el radio exterior de la turbina. Para un elemento diferencial de la turbina (dr), la tasa de cambio del momentum angular por la velocidad de rotacin (Q), puede igualarse al diferencial de la potencia extrada por el elemento de aspa. Esta expresin conduce al principio de conservacin de la energa por unidad de volumen, expresada como:
(PP,}^ pLr ' = p'r'S7'U, (2.1)
Y la potencia extrada por el elemento de aspa ser entonces:
W
-
4Y',= U (27rrdr
U'`
(2.2)
El trmino , representa (P 1 -P2)2nr drU el trabajo entregado por el fluido en la unidad de tiempo. ('hpUtz) 27cr drU, representa la energa cintica rotacional por unidad de tiempo, que se deja en la estela y la cual no se aprovecha en la extraccin de energa. (Pinilla, 2001)
2.1. Fuerzas aerodinmicas sobre el perfil
Las fuerzas aerodinmicas de sustentacin y arrastre pueden descomponerse en una componente axial (F X) y una componente tangencial (FT) (Ver figura No 2). Esta ltima es la responsable del movimiento rotacional del disco de la turbina. = L..senT Dcos9
Ex
2.
DISEO AERODINMICO
La extraccin de energa en el rotor de la turbina, implica un cambio en el momentum angular de la corriente de aire. Este cambio ocurre con la interaccin entre el flujo de aire y los alabes simtricos que componen la turbina. El adecuado diseo de esta consiste en que la componente de energa cintica rotacional aguas arriba de la turbina sea lo ms reducida posible. Dicha componente rotacional o velocidad inducida (lit) tiene un efecto directo en el desempeo aerodinmico de la turbina. (Pinilla,
= LcosT + Dsemp
(2.3)
1985)
La teora utilizada para el diseo aerodinmico de la turbina Wells se basa en la adecuada combinacin del anlisis del elemento de aspa (Blade Element Theory) y la teora de momentum. En la figura No. 2 se ilustra el tringulo de velocidades incidentes y las fuerzas aerodinmicas sobre el elemento de aspa. Para la turbina Wells, el ngulo que describe la velocidad efectiva (cp0)
En un perfil simtrico la direccin de la componente tangencial (F T) debe ser la misma para ngulos de ataque positivos y negativos, produciendo un giro en el rotor de la turbina en una nica direccin. La componente axial es la fuerza que servir para el diseo estructural de la turbina. (Raghunathan y Beattie, 1996). La condicin de diseo que optimiza la fuerza aerodinmica sobre el elemento de aspa y expresado como fuerza de sustentacin y arrastre ser entonces:
Facultad de Ingeniera
57
(2.4)
n
tipos de rotores Wells. El procedimiento se basa en rotores de 8 aspas radiales con un perfil NACA 0015, con la siguiente geometra: Tabla No. 1 Dimensiones (Curran y Gato, 1997)
Dimetro exterior
Combinando las componentes tangencial de la fuerza aerodinmica y la ecuacin de energa, expresada como:
(YP2 )= pUU,cotcp ,, (2.5)
del tringulo de velocidades [Figura No. 2] se puede determinar el ngulo de la velocidad efectiva ((p.) respecto al plano de rotacin de la turbina. (Pinilla, 2001). Este es el primer parmetro para determinar la geometra ptima del aspa de la turbina.
cot cp =
Dimetro del cubo Cuerda Solidez (a) Velocidad de rotacin
590 mm 400 mm 125 mm 0.640 2000 rpm
La solidez (a) del rotor es la relacin entre el rea ocupada por las aspas y el rea total descrita por el mismo. Para el caso de la turbina Wells este valor puede calcularse a travs de:
a=
Qr + U,
U
(2.6)
n D, (l +la)
(2.9)
El diseo ptimo de la turbina considera que la cada de presin a travs de la turbina sea constante para cada elemento de aspa. Asegurando esta condicin se obtiene la ecuacin 2.7, que determina la velocidad inducida (U) para cada posicin radial:
PP, 2pS2 r
(2.7)
La solidez del rotor varia con el reciproco del cuadrado de la velocidad especfica (U*) de diseo. (Pinilla, 1985) Es decir, cuando se requiere una alta velocidad para extraer energa, se requiere valores pequeos de solidez. Otra caracterstica es que el torque de arranque aumenta con el incremento de la solidez, de all el compromiso de balancear estas caractersticas para obtener un diseo optimo de turbina. Realizando un anlisis adimensional, se pueden determinar 4 grupos que tipifican el funcionamiento de la turbina. Estos grupos son funcin de la densidad del fluido, velocidad del aire, dimetro del rotor, cada de presin, velocidad angular, potencia extrada y torque. Estos nmeros se convierten en los parmetros conocidos inicialmente en el diseo del sistema. La eficiencia de la turbina puede ser determinada a travs de la relacin entre la potencia extrada y la energa neumtica disponible.
q
T
S2
Combinando esta expresiones se obtiene la ecuacin de diseo que determina la longitud y distribucin de la cuerda optima a lo largo del aspa. (Pinilla, 2001).
B CL;,
475r . ........
* Y` ,
, /1 + ((} )^
\,
; 1+ dP ,
' (2 +(rY
`1+(i%
(2.8)
3.
DIMENSIONAMIENTO PRELIMINAR DE LA
TURBINA WELLS
Q AP
(2.10)
3.1. Anlisis Dimensional Para determinar unas dimensiones preliminares del rotor de la turbina se recurre a los resultados reportados en Curran y Gato, 1997, en el cual se obtiene la eficiencia para varios
Los valores obtenidos en el experimento de la referencia para los parmetros adimensionales en estudio, son determinados en el punto mximo de eficiencia obtenida para la turbina y mostrados en la tabla No 2.
58
UNIVERSIDAD DE LOS ANDES
Para el valor obtenido de eficiencia, se produjo una velocidad en la corriente de aire de 9.30 m/s, un flujo de L37 m 3/s, una cada de presin de 2362 Pa, un torque en el disco de 7.35 N.m para una potencia til de 1924 W. (Curran
radica en la concentracin del flujo en el dimetro exterior, lo cual se refleja en una rpido arranque. [Figura No. 5] Entre ms delgado sea el perfil aerodinmico, su relacin sustentacin-arrastre ptimos aumenta para un determinado nmero de Reynolds. Se elige un perfil NACA 0018 con una relacin CL/CD de 37.2 que se presenta a un ngulo de ataque de 8 y su nmero de Reynolds de 160000. [Figura No. 4].
y Gato, 1997). Tabla No. 2 Resultados (Curran y Gato, 1997)
OPo* T* U*
11*
0.330 0.04 0.120 59,5%
3.3. Adecuacin geometra preliminar
Con los resultados mostrados en los parmetros adimensionales, se optimiza la geometra del rotor de la turbina con ayuda de las ecuaciones tericas de la seccin 2. Se conserva la consideracin inicial de un rotor de 8 aspas, con una velocidad angular de 3000 rpm. Se obtienen los valores de la velocidad inducida (Ut), ngulo de ataque (a=cpo), factor de diseo (BC . CLopt.m /4pr), nmero de flujo (U*) y velocidad efectiva de la corriente de aire en cualquier elemento de aspa para determinada posicin radial. [Tabla No.4]. La velocidad inducida (Ut) al igual que el ngulo (cpo), son mayores en el cubo y varan hasta un mnimo en la punta, este comportamiento produce una velocidad efectiva mayor en la punta. De igual manera la geometra del perfil representada en la longitud de cuerda, muestra un valor menor en el cubo hasta un mximo en la punta. Este tipo de aspas es conocido como sectorial, y optimiza la extraccin de la potencia disponible en la corriente de aire. [Figura No. 5]. La desventaja de las aspas sectoriales radica en su fabricacin. Otra opcin de construccin consistira en un diseo simplificado conformado por una aspa rectangular con una longitud de cuerda media. [Figura No. 6].
4.
DISEO DEL GENERADOR
3.2. Diseo de la geometra preliminar
Con los parmetros adimensionales del estudio anterior, se realiz un procedimiento de similitud geomtrica representado en el aseguramiento de la solidez (s) y la relacin de dimetros. Como primera aproximacin se elige un dimetro del rotor de 300 mm para una velocidad de rotacin de 3000 rpm. Para las condiciones atmosfricas de Bogot con una densidad de 0.950 Kg/m 3 a 2600 m, los resultados de diseo preliminares para la turbina son los siguientes:
Tabla No. 3 Resultados propuestos
Dimetro exterior Dimetro del Cubo Cuerda Torque Potencia Cada de presin Velocidad aire Caudal de aire Energa neumtica Eficiencia Global 300 mm 204 mm 63 mm 0.285 Nm
89.47 W
708.83 Pa 5.65 m/s 0.215 m3/s 152.32 W 59.8%
La relacin que se presenta entre la velocidad de rotacin y potencia generada, para el rotor propuesto se presenta en la figura No. 3. Se conserva el nmero de aspas para mantener el torque de arranque necesario. La disposicin de las aspas puede ser de manera rectangular o por sectores, la ventaja de las aspas sectoriales
La energa mecnica rotacional suministrada por la turbina debe ser transformada a energa elctrica a travs de un generador elctrico. Se opt por disear un generador de imanes permanentes el cual tiene un menor volumen y peso
Facultad de Ingeniera
59
en comparacin a una mquina tradicional equivalente. (Patrick y Stephen, 1996). Este generador se debe ajustar a los requerimientos de la turbina, tales como: bajo torque de arranque y generacin de potencia entre 50 -100 W. [Figura No. 7]. 4.1. Diseo del embobinado El generador consta de una parte esttica conformada por un ncleo magntico y 8 bobinas dispuestas simtricamente alrededor del ncleo. La parte giratoria es constituida por un volante de motocicleta NIPON DENSO, el cual contiene (4) imanes dispuestos radialmente. El voltaje inducido en el devanado depende del nmero de vueltas de cada una de las bobinas, el material del ncleo y la velocidad relativa entre el devanado y el campo. Para obtener un voltaje suficiente en el devanado, el generador debe tener la mxima cantidad posible de boTabla No. 4
r (mm) VT =S2r (m/s) 32,044 102 t 107 33,552 112 35,060 116 36,568 121 38,076 126 39,584 131 41,092 136 42,600 140 44,108 47,124 150 2
binas y vueltas en cada una de ellas. Considerando las condiciones geomtricas del volante, se disponen 8 bobinas con 16 vueltas cada una; dos bobinas por cada polo del volante. [Foto No. 1] Cada bobina se fabrica con alambre de cobre esmaltado AWG No. 19, con un dimetro de 0.912 mm, y una capacidad mxima de corriente de aproximadamente 12 A. Estas se disponen en un ncleo fabricado con lminas de acero al silicio con alta permeabilidad magntica. Debe tenerse en cuenta que las principales prdidas en una mquina elctrica se deben a las producidas en el ncleo y al cobre del devanado. (Wallace y Rodrguez, 2000) La eficiencia del generador se ve influenciada en primera medida, por el entrehierro que existe entre el ncleo y el campo magntico. Esta distancia en el generador diseado tiene aproximadamente 1.20 mm. (radiales)
Resultados ptimos (Pinilla, 2001)
(Po (grados)
ut (m/s) 11,642 11,119 10,641 10,202 9,798 9,425 9,079 8,758 8,458 7,917
7,375 7,215 7,054 6,894 6,737 6,582 6,431 6,283 6,140 5,866
Vete. (m/s) 44,051 45,028 46,050 47,111 48,207 49,334 50,489 51,668 52,869 55,330
U* C_DISEO c(mm) 0,176 0,2486 54,03 0,2489 56,63 0,169 0,161 0,2491 59,22 61,82 0,155 0,2493 64,41 0,149 0,2494 0,143 0,2496 67,00 0,138 0,2497 69,59 0,133 0,2499 72,18 0,128 0,2500 74,77 0,120 0,2502 79,94
4.2. Generacin de Corriente Cuando una carga externa es conectada al devanado del generador un circuito cerrado se completa, permitiendo la circulacin de una corriente. Es importante recordar que el disco de turbina Wells diseado estar sometido a un flujo de aire oscilante que implicar en el generador un voltaje a una frecuencia variable. Esta consideracin obliga a la utilizacin de un rectificador a la salida del generador para lograr un voltaje de naturaleza continua. El voltaje en el generador diseado es rectificado inicialmente con una conexin de diodos tipo puente de onda completa con una capacidad mxima de 25 A. La onda de voltaje obtenida en las terminales (3 y 4) del puente es directa, la cual debe ser corregida a una onda lo
ms uniforme posible. Para obtener un voltaje continuo, se colocan dos condensadores de 6800 mF y 4700 mF conectados en paralelo a la salida del puente rectificador. 4.3. Montaje experimental Para caracterizar el funcionamiento del generador se disea una carga elctrica conformada por 3 bombillos de alta de automvil dispuestos en paralelo, con el fin de obtener una corriente elctrica por cada una de las lneas. El generador es montado en un torno mecnico con valores fijos de velocidad entre 60 a 2245 RPM, verificando este valor con un tacmetro de luz estroboscopica. El objetivo de este montaje es determinar la potencia extrada del generador, obtenida por la corriente total que cir-
60
UNIVERSIDAD DE LOS ANDES
cula por el circuito elctrico y el voltaje a la salida del puente de rectificacin. [Foto No. 2] Finalmente se mide el voltaje en el dominio del tiempo con ayuda de un osciloscopio, para determinar el comportamiento de sus ondas dentro del sistema.
4.4. Resultados obtenidos
5.
INVESTIGACIN FUTURA
Los conocimientos tericos y experimentales hasta el momento enmarcan el adecuado acople y diseo del disco preliminar de la turbina y el generador requerido para extraer la potencia elctrica necesaria de 60 W.
5.1 Optimizacin del generador y de la turbina
En un generador de imanes permanentes la velocidad de giro del campo magntico es conocida como la velocidad sincrnica de rotacin y es determinada fcilmente por:
n, = l20 * f
(2.11)
Donde:
n s : Velocidad sincrnica de rotacin (RPM) f: Frecuencia del sistema (Hz) p: Nmero de polos (4)
Con los resultados preliminares del generador se debe optimizar la geometra del rotor de la turbina con el fin de definir la geometra del modelo. En esta etapa se define la construccin del rotor y los experimentos que permitirn el refinamiento del modelo terico y la obtencin de parmetros adimensionales. Esta etapa igualmente caracteriza el desempeo mecnico del generador, determinado por curvas de torque y el comportamiento frente a variaciones en la velocidad del rotor. Variaciones que se presentan por las oscilaciones en el flujo de aire a travs de la turbina.
5.2 Diseo final y adecuacin
La frecuencia medida de la onda de voltaje es de 84.7 Hz; considerando la ecuacin anterior se obtiene un valor de velocidad sincrnica de 2541 RPM, valor que es cercano a la velocidad de giro del volante del generador. El comportamiento de la onda final de voltaje se acerca a la forma deseada con pequeas fluctuaciones que debern ser corregidas posteriormente. El voltaje se reduce a medida que la carga en el sistema aumenta, considerando slo la rectificacin en el puente de onda completa el valor mximo obtenido es de 10.58 VDC. Este valor puede aumentarse con la adicin de los condensadores a la salida del puente rectificador. [Figura No. 8] La potencia mxima extrable es de 60.51 W para un valor de 9.10 VCC y una corriente total en el sistema de 6.65 A. [Figura No. 9] Para producir una velocidad de giro de 2500 RPM en el disco preliminar de la turbina Wells, se requiere una potencia mecnica de 51.78 W. Comparando con los resultados de la potencia elctrica del generador, sera posible utilizar el generador para alcanzar el objetivo de 50 W elctricos. El exceso de potencia elctrica podra ser utilizado para optimizar el dimetro exterior del disco o su velocidad de giro.
La etapa final consiste en la construccin del dispositivo sugerido con la adecuacin del generador elctrico y el dispositivo de suministro de aire por medio de un fuelle. Este ltimo involucra los limites ergonmicos, las condiciones de los equipos a los cuales se suministra energa, y el desarrollo de un modelo econmico para su inmediata o futura implementacin. La construccin del sistema propuesto, genera una base preliminar de conocimiento en nuestro pas para implementar y disear futuros sistemas de conversin de energa a pequea escala bajo el principio de turbinas Wells, principalmente en zonas remotas del pas.
6. CONCLUSIN
Los resultados experimentales obtenidos del generador elctrico muestran la posibilidad de extraer una potencia de 60 W a 2545 rpm. Los clculos tericos preliminares para dimensionar la turbina Wells muestran una potencia requerida de 51.78 W para producir un giro de 2500 rpm. Los resultados tericos y experimentales
Facultad de Ingeniera
61
se encuentran bien encaminados para realizar el acople ms adecuado de los elementos para la construccin del dispositivo final. Vale la pena mencionar que este trabajo es el resultado de la primera parte de una tesis de Maestra en Ingeniera Mecnica (Santa) con el apoyo de un estudiante en su proyecto de grado (Lozano). El desarrollo del dispositvo completamente probado en el laboratorio se culminar a mediados del mes de Mayo de 2002.
7.
REFERENCIAS
^
GENERADOR ELECTRICO EJE DE LA TURBINA
PERFIL NACA00XX SIMETRICO
DIRECCIN DE ROTACIN
I
V /
CUBO bE LA TURBINA
1. Curran, R., Gato, L.M.C. The Energy conversion of several types of Wells turbine designs. Proc. Instn. Mech. Engrs, Part A, 1997, 211, 133-145. 2. Dunn, P.D. Renewable Energies: Sources, Conversion and Application. IEE energy series Vol. 2, Peter Peregrinus Ltd., 1986, 290-336. 3. Pinilla, Alvaro E., Wind Powered Pumping Systems for Colombia. Ph. D. Thesis, Reading University, 1985, 4366. 4. Pinilla, Alvaro E., Diseo de turbina Wells. Documento Interno de Trabajo, Universidad de los Andes, 2001, Bogota. 5. Patrick, R., Dale. Stephen, W, Fardo. Rotating electrical machines and power systems. Prentice Hall, second edition, 1996. 6. Raghunathan, S., Beattie, W.,C. Aerodynamic performance of contra-rotating Wells turbine for wave energy generation. Proc. Instn. Mech. Engrs, Part A, 210, 1996. 7. Sharpe, D., J. Layman's Guide to the aerodynamics of Wind Turbines. Wind energy conversion, 1988, 229-242. 8. Wallace, C., Rogel. Rodrguez, E., Alejandro. Diseo de un generador de flujo axial con imanes permanentes para aplicaciones elicas. Revista facultad de Ingeniera, U.T.A. Chile. Vol.7, 2000.
FLUJO DE AIRE OSCILANTE
Figura No. 1 Diagrama esquemtico de una turbina Wells
Ft
Fx D
VcfaV1
^V^C
Figura No. 2 Velocidades y fuerzas aerodinmicas sobre perfil. (finilla, 1985)
DESEMPEO PARA LA TURBINA WELLSO(300_mm)
1000
1500
2000
2500
3000
3500
4000
Velocidad de rotacin fi(rpm)
Figura No. 3 Relacin entre la velocidad angular y potencia generada
62
UNIVERSIDAD DE
LOS
ANDES
:m:w (o +^ y^ i
z mx
^ .m9
xw xw
:4692
a .g
^
(D^ asaJ^ ^ 1ort^ e^ m3o gc lyda^e/em^^ ddw(Camp^l w Imi, geqaeeatr Estructura ptydpdC} 4. 3. ^e26 qtocr ma^ 6 6, Rodamiento G21 2RS e 2 ^ ^ adn^ , 2 #: ;u o m
2.
E2 2^
2^ nde m^ acin \deoo^ CrG erre y 30^ oS y ^I
2 3
r::!:: :::r;r:<:
;a.
:
^ ^ ^
Bombillos
)*BB ; Bombillos
Comportamiento del vottajz aau>w
70.00
bombillo
^ .
2 bom bilio b :
500
1000 RPM
1500
2000
m::
immG potencia 7amd
6}
ERO DAD DE LDS ANDES