Se formula un modelo para la predicción de los parámetros característicos de una turbina eólica de eje horizontal basado en la aplicación de la teoría de la cantidad de movimiento, con la incorporación de la rotación del disco actuador. Se determinan los coeficientes de potencia, tracción y torque, ahora afectados por la rotación del disco que simula a la turbina, determinando las pérdidas de energía en un modelo invíscido en donde se simula la presencia de una estela rotacional corriente abajo del actuador.A simplified model has been developed to predict the characteristic coefficients of horizontal-axis wind turbines. The model is based on the momentum theory and the inclusion of a rotating actuator disc. The power, axial force and torque coefficients are evaluated under the condition of disc rotation, in order to determine the energy losses in a model where a rotational inviscid wake is developed downstream of the rotating actuator disc that simulates the wind turbine.Asociación Argentina de Energías Renovables y Medio Ambiente (ASADES

Prado, Ricardo Adolfo

Servicio de Difusión de la Creación Intelectual

       TEORÍA DEL DISCO ACTUADOR ROTANTE – FUENTE PUNTUAL PARA TURBINAS EÓLICAS DE EJE HORIZONTAL   Ricardo A. Prado Departamento de Mecánica Aplicada - Facultad de Ingeniería - Universidad Nacional del Comahue calle Buenos Aires N°1400, Q8300BCX Neuquén - Argentina Tel. +54 - (0)299 – 449·0338, email: prado@uncoma.edu.ar  RESUMEN: Se formula un modelo para la predicción de los parámetros característicos de una turbina eólica de eje horizontal basado en la aplicación de la teoría de la cantidad de movimiento, con la incorporación de la rotación del disco actuador. Se determinan los coeficientes de potencia, tracción y torque, ahora afectados por la rotación del disco que simula a la turbina, determinando las pérdidas de energía en un modelo invíscido en donde se simula la presencia de una estela rotacional corriente abajo del. actuador.  Palabras clave: modelo invíscido, disco actuador, rotación, turbina eólica, eje horizontal.  INTRODUCCIÓN  Entre los más simples modelos desarrollados para predecir los coeficientes característicos de las turbinas eólicas de eje horizontal, se encuentra el de Betz, originado por la adaptación de la teoría de la cantidad de movimiento de Rankine-Froude, y descrito, entre otros, por Le Gourières (1982). Este modelo no incluye la rotación del denominado disco actuador que simula a la turbina ni pérdidas de energía en el aire que cruza al disco. Entre los modelos que incorporan la rotación del disco se encuentra la adaptación al modelo de Glauert, el cual representa una extensión del mencionado modelo de Betz, con lo cual comparte las falencias de la mencionada representación (Eggleston y Stoddard, 1987).  Tomando como base la propuesta de teoría presentada por Dyment (1989), el modelo descrito inicialmente por Prado (1995) simula al tubo de corriente -el cual encierra al aire que fluye a través del disco actuador que representa a la turbina cólica de eje horizontal-, mediante la superposición potencial de una fuente puntual y una corriente uniforme. Esta representación conduce a una serie de parámetros característicos de la turbina eólica que incluyen la presencia de pérdidas de energía, con lo cual el modelo difiere sustancialmente del de Betz, si bien no se incluyó el efecto de la rotación del aire que fluye detrás del disco actuador.   En un modelo anteriormente propuesto por Prado (1996), el cual consideraba una representación semejante a la presentada en Eggleston y Stoddard (1987), se incorporaba un disco actuador rotante que imprimía dicha rotación solamente a la porción del fluido que se encuentra corriente abajo del mismo, pero se asumía que el aire respondía fundamentalmente a un modelo irrotacional de flujo, con la consecuente distribución de presiones en el tubo de corriente generado aguas abajo de la posición del actuador. Con la finalidad de mejorar a dicho modelo, la presente reformulación incorpora la presencia de un flujo rotacional, el cual simula la presencia de una estela helicoidal de radio no constante a medida que se desarrolla corriente abajo del disco actuador, totalmente contenida dentro del tubo de corriente τ, cuyo interior de la porción posterior rota, si bien no lo hace como un cuerpo rígido.  DEFINICIONES DEL MODELO  En el presente trabajo se considera un disco actuador de radio R, infinitamente delgado, perfectamente permeable al flujo, ubicado en coincidencia con el origen de coordenadas del eje x, donde se ubica asimismo la fuente puntual de acuerdo al modelo propuesto por Prado (1995). El disco se encuentra inmerso en el flujo coaxial de un fluido idealmente no viscoso. Adicionalmente, dicho disco rota a velocidad angular uniforme Ω, dejando tras de sí un flujo rotacional que simula la presencia de una estela rotacional no viscosa.   El tubo de corriente que envuelve al aire que pasa a través del disco actuador, se identifica mediante τ. Por definición de tubo de corriente, no existe flujo alguno a través de sus paredes. Por fuera de dicho tubo de corriente, el flujo resulta potencial (irrotacional y asumido no viscoso) por lo cual no presenta pérdidas y admite la aplicación de la ecuación de Bernoulli. La componente axial del flujo de aire dentro de τ es considera unidimensional, es decir V = V(x), asumiéndose entonces la presencia de deslizamiento entre los flujos interno y externo a τ. Como fuera indicado, solo el flujo de aire interno a la porción posterior de τ es el que rota.   El flujo dentro de τ  en x = -∞  está representado por la corriente de aire no perturbada de velocidad V∞. Sobre la posición del disco actuador (x = 0), la velocidad axial resulta reducida a V0 = V∞ (1-a), siendo a el coeficiente de reducción de velocidad en el disco. Suficientemente lejos corriente abajo, la velocidad axial dentro de τ resulta aún más reducida al valor                 ASADES Avances en Energías Renovables y Medio Ambiente Vol. 8, Nº 2, 2004.  Impreso en la Argentina.  ISSN  0329-5184  06.07   V∞ (1-a)/(1+a) (Prado, 1995). La Figura 1 muestra un segmento del tubo de corriente τ, representado alrededor de la posición del disco actuador, para un coeficiente de reducción de velocidad en el disco, a, igual a 0,3.  Para aplicar la metodología del volumen de control, interesa particularmente determinar de la geometría del tubo de corriente τ, los radios infinitamente lejos por delante (x = -∞) y por detrás (x = +∞) de la posición del disco actuador, resultando (Prado, 1995 a1RR −=∞−                  a1RR +=∞+                  1a0 ≤≤  (1) La introducción de la rotación del disco actuador conduce a la introducción de una nueva variable adimensional, conocida como velocidad específica o relación de velocidad en el borde del disco, definida por la razón entre la velocidad tangencial en el borde del disco y la velocidad del viento no perturbada, es decir, ∞∞Ω=λVR  (2) Dado que la velocidad del viento se reduce hasta el valor V0 = V∞ (1-a) en el disco actuador, resulta en la estación x= 0 la siguiente definición de la razón de velocidades sobre el disco )a1(VR00−λ=Ω=λ ∞  (3)  Figura 1: Sistema coordenado, porción del  tubo de corrienteτ  y disco actuador, siendo a = 0,3.  ECUACIONES GENERALES DEL MODELO  Ecuaciones integrales del método del volumen de control Las ecuaciones de gobierno del movimiento de fluidos (Shames, 1992), escritas bajo la consideración de la metodología del volumen de control, considerando a las mismas bajos condiciones estacionarias, resultan las siguientes ecuaciones, evaluadas sobre la superficie (SC) que envuelve al volumen de control adecuadamente definido: ecuación de continuidad, 0)dAnˆ·V(CS=ρ∫∫ ρ , (4) ecuación de la cantidad de movimiento (segunda ley de Newton) ∫∫∑ ρCS)dAnˆ·V(V=Fρρρ, (5) ecuación del momento de la cantidad de movimiento lineal ∫∫∑ ρ×CS0 )dAnˆ·V()Vr(=Mρρρρ (6) 06.08   y la ecuación para la conservación de la energía (primer principio de la termodinámica) despreciando los efectos gravitatorios ∫∫ ρρ=−CS221s )dAnˆ·V(p+V+u~WQρ&&  (7) donde,  ρ densidad del aire, ρV vector velocidad local del flujo (incluyendo al viento y a la rotación de la estela), con módulo V, ∃n   versor normal a la superficie de control, ρF  fuerzas actuantes sobre el volumen de control, tanto superficiales como másicas, activas o reactivas, tanto inerciales como las provenientes de posibles efectos no inerciales, 0Mρ momentos referidos al origen de coordenadas de las fuerzas actuantes sobre el volumen de control, ρr  vector posición al origen de coordenadas, &Q  calor suministrado, por unidad de tiempo, &Ws  trabajo cedido por el aire, por unidad de tiempo, como potencia al eje, ~u  energía interna por unidad de masa y  p  presión estática absoluta  Asimismo, se ha considerado que los volúmenes de control, y por ende sus superficies, son fijos e indeformable, el flujo es estacionario y el sistema de referencia inercial se encuentra fijo, con origen en la fuente puntual ubicada en el centro del disco. El sistema coordenado utilizado es el cilíndrico, aprovechando la simetría axial del problema., denotando a la coordenada radial mediante r y a la axial mediante x.   Efectos de la rotación de la corriente Previo al planteo de las ecuaciones integrales de gobierno, se deben plantear las condiciones de contorno del problema, por lo cual debe determinar la distribución de presiones lejos corriente abajo originada por el efecto de rotación de la corriente dentro de τ. Según Schlichting (1972), resulta que [ ]2)x(r)r,x(rpωρ=∂∂  (8) para compensar el efecto de la aceleración centrípeta presente en el flujo. Entonces, [ ] )x(C)x(r)r,x(p 2221 +ωρ=  (9) Dado que en x = +∞  y por fuera de τ, la presión estática resulta ser la atmosférica ∞=++∞ p)a1R,(p  (10) se determina entonces la constante de integración C(∞), con lo cual [ ] [ ]22221 r)a1(R)(rp)r,(p −++∞ωρ−=−+∞ ∞             a1Rr0 +≤≤  (11) Promediando dicha distribución sobre la sección de τ en x = +∞ , resulta ( ))a1(1q'''p)(p 221manman+λ−==+∞∞∞ (12) donde el subíndice man indica que la medición es manométrica y ( )221 Vq ∞∞ ρ=  (13) denota la presión dinámica de la corriente libre.  Siguiendo un proceso semejante al detallado y planteando el equilibrio de fuerzas, incluida la variación de la cantidad de movimiento axial, resultan las presiones a ambos lados del disco actuador ( )  λ−+−−==+ ∞∞ 2)a1()a1(a24aq''p)0(p222manman [ ])a78(aq41'p)0(p manman −==− ∞   (14) observándose que la reducción de la presión detrás del disco rotante (con respecto a la condición sin rotación, λ ∞ = 0), trae como consecuencia el incremento de la tracción sobre el rotor. 06.09   DOMINIO I DOMINIO II  ω(0) = Ω p’’(r) p∞p∞p∞                  vx = V∞ V∞p∞                  vx = V∞V∞ V∞(1-a)/(1+a)  V∞ p’’’(r) ω(+∞) (r) x rR -∞ R +∞ R  τ  V∞(1-a) vx (0,r) = V∞ p(0,r)  p’ ω(x<0) = 0 Obtención de los coeficientes característicos La Figura 2 muestra, esquemáticamente, los dominios en que subdivide el espacio para el planteo de las ecuaciones integrales de gobierno, con sus correspondientes condiciones de contorno.                               Figura 2: Dominios para la aplicación del método del volumen de control. Condiciones de contorno.  Del planteo de la ecuación de la cantidad de movimiento axial (5) resulta el coeficiente de tracción ( )2)a1()a1(a4AqTC2T∞∞λ++−== (15) donde A = π R² denota el área del disco. La tracción T tiene sentido de x(+), representando una resistencia a la corriente. Planteando la ecuación de la conservación de la cantidad de movimiento angular (6) tanto en la totalidad el dominio II, como  solamente en el interior del tubo de corriente inmerso en dicho dominio, resultan )a1(RR2+Ω≡Ω=ω∞+∞+ (16) y el coeficiente de torque )a1(ARqMC 0Q −λ== ∞∞ (17)  Figura 3: variaciones a lo largo del eje del tubo de corrienteτ de las magnitudes adimensionales V(x)/V∞  y ω(x)/Ω.  06.10   La variación axial de la velocidad angular de la estela, la cual se asume iniciarse discontinuamente una vez que el aire atraviesa el disco, es representada en la Figura 3, conjuntamente con la variación de la velocidad dentro del tubo de corriente, considerando un coeficiente de reducción de velocidad en el disco a = 0,3. En general, estas resultan  2)x(RR)x( Ω=ωτ          x > 0                              2)x(RR)a1(V)x(V −=τ∞          x > -∞     (18) donde Rτ (x) es el radio del tubo de corriente τ en la estación x, el cual verifica la representación funcional presentada en (Prado, 1995), y resulta asimismo una función del parámetro a. En este modelo se asume que no existe deslizamiento entre las velocidades tangenciales del rotor y las del aire en contacto con la cara posterior del disco.  Como lo indican Eggleston y Stoddard (1987), la energía rotacional impartida al aire no resulta económicamente recuperable, por lo cual representa una energía desperdiciada en hacer girar la estela, limitándose así aún más la eficiencia del rotor. Dicha energía es proporcional al producto entre la cupla o torque y la velocidad angular de rotación, M0Ω, definiendo el correspondiente coeficiente de potencia ( ) )a1(AVqMC 20Pot −λ=Ω=∞∞∞ (19)  Vale aclarar que los únicos efectos no inerciales que aparecen en el presente análisis de deben a la presencia de aceleraciones centrípetas ( )rρρρ ×ω×ω  en la estela, en donde aparece la componente tangencial de la velocidad del aire, vθ = ωr. Por otra parte, la aceleración de Coriolis es nula pues V//ρρω y no existen aceleraciones lineales o angulares del tubo de corriente.  De la aplicación de la ecuación de la energía (7) en el dominio I, en el dominio II y en la totalidad del espacio, resultan 22II a432/)V(u~=∆=ζ∞ (20) para el coeficiente de pérdida de energía en la porción I de τ, +−−=∆=ζ∞2222IIII )a1()a1(41a2/)V(u~  (21) para el coeficiente de pérdida de energía en la porción II de τ, +−−−=ζ+ζ−=ζ 222III )a1()a1(1a)a1())(a1(C  (22) para el coeficiente de pérdida de potencia total, y )a1()a1(a4AVqWC2sWs+−==∞∞& (23) para el coeficiente potencia al eje de la turbina.  Dado que las pérdidas representadas por la ecuación (22) y graficadas en la Figura 4 no son consideradas por el modelo de Betz, resulta un coeficiente de potencia al eje menor que el propuesto por dicha teoría. En la Figura 5 se grafica la variación del coeficiente de potencia al eje, ecuación (23), con el factor de reducción de velocidad en el disco. Se observa que el máximo valor del coeficiente apenas supera el valor 0,45 frente al (reconocidamente elevado) 0,60 que predice la teoría de Betz. Por otra parte, el presente modelo admite soluciones para a > 0,5 mientras que el modelo de Betz deja de ser válido en el rango 0,5 < a ≤ 1. Por su parte, la teoría de Dyment (1989) subestima el valor del máximo coeficiente de potencia obtenible (resultando éste aproximadamente 0,30). En cuanto a la comparación con resultados experimentales, Heijdra y Hendriks (2001) muestran que los máximos coeficientes de potencia observados en su análisis comparativo se encuentran entre 0,37 y 0,44 (Table 2, p.28).   Por otra parte, para definir el adecuado valor del parámetro λ∞ del modelo, se considera la relación presentada por Le Gourières (1982) proveniente de un estudio llevado a cabo por Hutter (p.54 op.cit.), donde se representa la variación funcional entre la velocidad específica λ∞ y la solidez del rotor, σ, que es la razón entre la superficie de las palas y el área barrida por las mismas. Dicha relación funcional (ver figura 54, p.55 op.cit) es aproximada como  [ ]( ) ( )∞λ−≈σ 1010 log523,12%log   (24) En el caso un disco, donde la razón de solidez es unitaria (σ = 100%), le corresponde por la relación (24) el valor de velocidad específica λ∞ = 1. Por lo cual, resulta que Ω = V∞/R es el valor de rotación más representativo de un disco actuador.06.11     Figura 4: Coeficiente de pérdida de potencia dentro del tubo de corriente τ, como función del parámetro a.    Figura 5:Coeficiente de potencia en función del factor de disminución de velocidad en el disco, a. Modelo presente y de Betz  CONCLUSIONES  Se propone un modelo de disco actuador rotante que complementa a la representación del tubo de corriente que envuelve al aire que fluye a través del rotor de una turbina eólica de eje horizontal mediante la superposición potencial de una corriente uniforme de dirección del eje del disco y de una fuente puntual ubicada en coincidencia con el centro del disco.   Se desarrollan las expresiones de los coeficientes característicos de un rotor eólico considerando el modelo descrito. Si bien el presente modelo es todavía muy simplificado, representa un avance frente al modelo de Betz y a aquellos modelos que, incorporando la rotación del disco, mantienen las hipótesis de la teoría de Betz en su desarrollo.  REFERENCIAS  Eggleston D.M. and Stoddard F.S. (1987) Wind Turbine Engineering Design. Van Nostrand Reinhold Co. Inc., New York. Dyment A. (1989) A modified form of the Betz’ wind turbine theory including losses. J. Fluids Engineering 111, 356-358. Heijdra J.J. and Hendriks H.B. (ed.) (2001) DOWEC Concept Study – Task 11: Selection of feasible concepts. ECN Windenergie, ECN-C--01-079. Le Gourières D. (1982) Energie éolienne. Deuxième édition, Eyrolles, Paris. Prado R.A. (1994) Reformulación de la teoría de la cantidad de movimiento aplicada a turbinas eólicas. Anales del VI Congreso Nacional de Ingeniería Mecánica, Santiago de Chile, pp.723-728. Prado R.A. (1995) Reformulation of the momentum theory applied to wind turbines. Journal of Wind Engineering and Industrial Aerodynamics 58, 277-292. Prado R.A. (1996) Modelo de disco actuador rotante para la simulación de una turbina eólica de eje horizontal. ASADES’96, Mar del Plata. Schlichting H. (1972) Teoría de la capa límite, Ediciones URMO, Bilbao, España. Shames I.H. (1992) Mecánica de Fluidos. McGraw-Hill, Colombia.  ABSTRACT: A simplified model has been developed to predict the characteristic coefficients of horizontal-axis wind turbines. The model is based on the momentum theory and the inclusion of a rotating actuator disc. The power, axial force and torque coefficients are evaluated under the condition of disc rotation, in order to determine the energy losses in a model where a rotational inviscid wake is developed downstream of the rotating actuator disc that simulates the wind turbine.  Keywords: inviscid model, actuator disc, rotation, horizontal-axis wind turbine. 06.12 

Teoría del disco actuador rotante : Fuente puntual para turbinas eólicas de eje horizontal

SEDICI - Repositorio de la UNLP

  
 
 
 
 
 
TEORÍA DEL DISCO ACTUADOR ROTANTE – FUENTE PUNTUAL 
PARA TURBINAS EÓLICAS DE EJE HORIZONTAL 
 
 
Ricardo A. Prado 
Departamento de Mecánica Aplicada - Facultad de Ingeniería - Universidad Nacional del Comahue 
calle Buenos Aires N°1400, Q8300BCX Neuquén - Argentina 
Tel. +54 - (0)299 – 449·0338, email: prado@uncoma.edu.ar 
 
RESUMEN: Se formula un modelo para la predicción de los parámetros característicos de una turbina eólica de eje 
horizontal basado en la aplicación de la teoría de la cantidad de movimiento, con la incorporación de la rotación del disco 
actuador. Se determinan los coeficientes de potencia, tracción y torque, ahora afectados por la rotación del disco que simula 
a la turbina, determinando las pérdidas de energía en un modelo invíscido en donde se simula la presencia de una estela 
rotacional corriente abajo del. actuador. 
 
Palabras clave: modelo invíscido, disco actuador, rotación, turbina eólica, eje horizontal. 
 
INTRODUCCIÓN 
 
Entre los más simples modelos desarrollados para predecir los coeficientes característicos de las turbinas eólicas de eje 
horizontal, se encuentra el de Betz, originado por la adaptación de la teoría de la cantidad de movimiento de Rankine-Froude, 
y descrito, entre otros, por Le Gourières (1982). Este modelo no incluye la rotación del denominado disco actuador que 
simula a la turbina ni pérdidas de energía en el aire que cruza al disco. Entre los modelos que incorporan la rotación del disco 
se encuentra la adaptación al modelo de Glauert, el cual representa una extensión del mencionado modelo de Betz, con lo 
cual comparte las falencias de la mencionada representación (Eggleston y Stoddard, 1987). 
 
Tomando como base la propuesta de teoría presentada por Dyment (1989), el modelo descrito inicialmente por Prado (1995) 
simula al tubo de corriente -el cual encierra al aire que fluye a través del disco actuador que representa a la turbina cólica de 
eje horizontal-, mediante la superposición potencial de una fuente puntual y una corriente uniforme. Esta representación 
conduce a una serie de parámetros característicos de la turbina eólica que incluyen la presencia de pérdidas de energía, con lo 
cual el modelo difiere sustancialmente del de Betz, si bien no se incluyó el efecto de la rotación del aire que fluye detrás del 
disco actuador.  
 
En un modelo anteriormente propuesto por Prado (1996), el cual consideraba una representación semejante a la presentada en 
Eggleston y Stoddard (1987), se incorporaba un disco actuador rotante que imprimía dicha rotación solamente a la porción 
del fluido que se encuentra corriente abajo del mismo, pero se asumía que el aire respondía fundamentalmente a un modelo 
irrotacional de flujo, con la consecuente distribución de presiones en el tubo de corriente generado aguas abajo de la posición 
del actuador. Con la finalidad de mejorar a dicho modelo, la presente reformulación incorpora la presencia de un flujo 
rotacional, el cual simula la presencia de una estela helicoidal de radio no constante a medida que se desarrolla corriente 
abajo del disco actuador, totalmente contenida dentro del tubo de corriente τ, cuyo interior de la porción posterior rota, si bien 
no lo hace como un cuerpo rígido. 
 
DEFINICIONES DEL MODELO 
 
En el presente trabajo se considera un disco actuador de radio R, infinitamente delgado, perfectamente permeable al flujo, 
ubicado en coincidencia con el origen de coordenadas del eje x, donde se ubica asimismo la fuente puntual de acuerdo al 
modelo propuesto por Prado (1995). El disco se encuentra inmerso en el flujo coaxial de un fluido idealmente no viscoso. 
Adicionalmente, dicho disco rota a velocidad angular uniforme Ω, dejando tras de sí un flujo rotacional que simula la 
presencia de una estela rotacional no viscosa.  
 
El tubo de corriente que envuelve al aire que pasa a través del disco actuador, se identifica mediante τ. Por definición de tubo 
de corriente, no existe flujo alguno a través de sus paredes. Por fuera de dicho tubo de corriente, el flujo resulta potencial 
(irrotacional y asumido no viscoso) por lo cual no presenta pérdidas y admite la aplicación de la ecuación de Bernoulli. La 
componente axial del flujo de aire dentro de τ es considera unidimensional, es decir V = V(x), asumiéndose entonces la 
presencia de deslizamiento entre los flujos interno y externo a τ. Como fuera indicado, solo el flujo de aire interno a la 
porción posterior de τ es el que rota.  
 
El flujo dentro de τ  en x = -∞  está representado por la corriente de aire no perturbada de velocidad V∞. Sobre la posición del 
disco actuador (x = 0), la velocidad axial resulta reducida a V0 = V∞ (1-a), siendo a el coeficiente de reducción de velocidad 
en el disco. Suficientemente lejos corriente abajo, la velocidad axial dentro de τ resulta aún más reducida al valor   
              ASADES 
Avances en Energías Renovables y Medio Ambiente 
Vol. 8, Nº 2, 2004.  Impreso en la Argentina.  ISSN  0329-5184 
 
06.07 
  
V∞ (1-a)/(1+a) (Prado, 1995). La Figura 1 muestra un segmento del tubo de corriente τ, representado alrededor de la posición 
del disco actuador, para un coeficiente de reducción de velocidad en el disco, a, igual a 0,3. 
 
Para aplicar la metodología del volumen de control, interesa particularmente determinar de la geometría del tubo de corriente 
τ, los radios infinitamente lejos por delante (x = -∞) y por detrás (x = +∞) de la posición del disco actuador, resultando 
(Prado, 1995 
a1RR −=
∞−
                  a1RR +=
∞+
                  1a0 ≤≤  (1) 
La introducción de la rotación del disco actuador conduce a la introducción de una nueva variable adimensional, conocida 
como velocidad específica o relación de velocidad en el borde del disco, definida por la razón entre la velocidad tangencial 
en el borde del disco y la velocidad del viento no perturbada, es decir, 
∞
∞
Ω
=λ
V
R  (2) 
Dado que la velocidad del viento se reduce hasta el valor V0 = V∞ (1-a) en el disco actuador, resulta en la estación x= 0 la 
siguiente definición de la razón de velocidades sobre el disco 
)a1(V
R
0
0
−
λ
=
Ω
=λ ∞  (3) 
 
Figura 1: Sistema coordenado, porción del  tubo de corrienteτ  y disco actuador, siendo a = 0,3. 
 
ECUACIONES GENERALES DEL MODELO 
 
Ecuaciones integrales del método del volumen de control 
Las ecuaciones de gobierno del movimiento de fluidos (Shames, 1992), escritas bajo la consideración de la metodología del 
volumen de control, considerando a las mismas bajos condiciones estacionarias, resultan las siguientes ecuaciones, evaluadas 
sobre la superficie (SC) que envuelve al volumen de control adecuadamente definido: 
ecuación de continuidad, 
0)dAnˆ·V(
CS
=ρ∫∫ ρ , (4) 
ecuación de la cantidad de movimiento (segunda ley de Newton) 
∫∫∑ ρ
CS
)dAnˆ·V(V=F
ρρρ
, (5) 
ecuación del momento de la cantidad de movimiento lineal 
∫∫∑ ρ×
CS
0 )dAnˆ·V()Vr(=M
ρρρρ
 (6) 
06.08 
  
y la ecuación para la conservación de la energía (primer principio de la termodinámica) despreciando los efectos gravitatorios 
∫∫ ρ



ρ
=−
CS
2
2
1
s )dAnˆ·V(
p+V+u~WQ
ρ&&  
(7) 
donde, 
 ρ densidad del aire, ρ
V vector velocidad local del flujo (incluyendo al viento y a la rotación de la estela), con módulo V, 
∃n   versor normal a la superficie de control, ρ
F  fuerzas actuantes sobre el volumen de control, tanto superficiales como másicas, activas o reactivas, tanto inerciales 
como las provenientes de posibles efectos no inerciales, 
0M
ρ
 momentos referidos al origen de coordenadas de las fuerzas actuantes sobre el volumen de control, 
ρ
r  vector posición al origen de coordenadas, 
&Q  calor suministrado, por unidad de tiempo, 
&Ws  trabajo cedido por el aire, por unidad de tiempo, como potencia al eje, 
~u  energía interna por unidad de masa y 
 p  presión estática absoluta 
 
Asimismo, se ha considerado que los volúmenes de control, y por ende sus superficies, son fijos e indeformable, el flujo es 
estacionario y el sistema de referencia inercial se encuentra fijo, con origen en la fuente puntual ubicada en el centro del 
disco. El sistema coordenado utilizado es el cilíndrico, aprovechando la simetría axial del problema., denotando a la 
coordenada radial mediante r y a la axial mediante x.  
 
Efectos de la rotación de la corriente 
Previo al planteo de las ecuaciones integrales de gobierno, se deben plantear las condiciones de contorno del problema, por lo 
cual debe determinar la distribución de presiones lejos corriente abajo originada por el efecto de rotación de la corriente 
dentro de τ. Según Schlichting (1972), resulta que 
[ ]2)x(r)r,x(
r
p
ωρ=
∂
∂  (8) 
para compensar el efecto de la aceleración centrípeta presente en el flujo. Entonces, 
[ ] )x(C)x(r)r,x(p 2221 +ωρ=  (9) 
Dado que en x = +∞  y por fuera de τ, la presión estática resulta ser la atmosférica 
∞
=++∞ p)a1R,(p  (10) 
se determina entonces la constante de integración C(∞), con lo cual 
[ ] [ ]22221 r)a1(R)(rp)r,(p −++∞ωρ−=−+∞ ∞             a1Rr0 +≤≤  (11) 
Promediando dicha distribución sobre la sección de τ en x = +∞ , resulta 
( )
)a1(
1q'''p)(p 22
1
manman
+
λ−==+∞
∞∞
 (12) 
donde el subíndice man indica que la medición es manométrica y 
( )221 Vq ∞∞ ρ=  (13) 
denota la presión dinámica de la corriente libre. 
 
Siguiendo un proceso semejante al detallado y planteando el equilibrio de fuerzas, incluida la variación de la cantidad de 
movimiento axial, resultan las presiones a ambos lados del disco actuador 
( ) 

 λ
−
+
−
−==+ ∞
∞ 2)a1(
)a1(a2
4
aq''p)0(p
222
manman
 
[ ])a78(aq
4
1'p)0(p manman −==− ∞  
 
(14) 
observándose que la reducción de la presión detrás del disco rotante (con respecto a la condición sin rotación, λ ∞ = 0), trae 
como consecuencia el incremento de la tracción sobre el rotor. 
06.09 
  
DOMINIO I DOMINIO II 
 ω(0) = Ω 
p’’(r) 
p∞
p∞
p∞                  vx = V∞
 V∞
p∞                  vx = V∞
V∞
 V∞(1-a)/(1+a) 
 V∞ 
p’’’(r) 
ω(+∞) (r) 
x 
r
R -∞ 
R +∞ 
R  
τ 
 V∞(1-a) 
vx (0,r) = V∞ 
p(0,r)  
p’ 
ω(x<0) = 0 
Obtención de los coeficientes característicos 
La Figura 2 muestra, esquemáticamente, los dominios en que subdivide el espacio para el planteo de las ecuaciones 
integrales de gobierno, con sus correspondientes condiciones de contorno. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Figura 2: Dominios para la aplicación del método del volumen de control. Condiciones de contorno. 
 
Del planteo de la ecuación de la cantidad de movimiento axial (5) resulta el coeficiente de tracción 
( )
2)a1(
)a1(a4
Aq
TC
2
T
∞
∞
λ
+
+
−
==
 
(15) 
donde A = π R² denota el área del disco. La tracción T tiene sentido de x(+), representando una resistencia a la corriente. 
Planteando la ecuación de la conservación de la cantidad de movimiento angular (6) tanto en la totalidad el dominio II, como  
solamente en el interior del tubo de corriente inmerso en dicho dominio, resultan 
)a1(R
R
2
+
Ω
≡



Ω=ω
∞+
∞+
 
(16) 
y el coeficiente de torque 
)a1(
ARq
MC 0Q −λ== ∞
∞
 (17) 
 
Figura 3: variaciones a lo largo del eje del tubo de corrienteτ de las magnitudes adimensionales V(x)/V∞  y ω(x)/Ω.  
06.10 
  
La variación axial de la velocidad angular de la estela, la cual se asume iniciarse discontinuamente una vez que el aire 
atraviesa el disco, es representada en la Figura 3, conjuntamente con la variación de la velocidad dentro del tubo de corriente, 
considerando un coeficiente de reducción de velocidad en el disco a = 0,3. En general, estas resultan  
2
)x(R
R)x( 



Ω=ω
τ
          x > 0                              
2
)x(R
R)a1(V)x(V 



−=
τ
∞
          x > -∞     
(18) 
donde Rτ (x) es el radio del tubo de corriente τ en la estación x, el cual verifica la representación funcional presentada en 
(Prado, 1995), y resulta asimismo una función del parámetro a. En este modelo se asume que no existe deslizamiento entre 
las velocidades tangenciales del rotor y las del aire en contacto con la cara posterior del disco. 
 
Como lo indican Eggleston y Stoddard (1987), la energía rotacional impartida al aire no resulta económicamente recuperable, 
por lo cual representa una energía desperdiciada en hacer girar la estela, limitándose así aún más la eficiencia del rotor. Dicha 
energía es proporcional al producto entre la cupla o torque y la velocidad angular de rotación, M0Ω, definiendo el 
correspondiente coeficiente de potencia 
( ) )a1(
AVq
MC 20Pot −λ=
Ω
=
∞
∞∞
 (19) 
 
Vale aclarar que los únicos efectos no inerciales que aparecen en el presente análisis de deben a la presencia de aceleraciones 
centrípetas ( )rρρρ ×ω×ω  en la estela, en donde aparece la componente tangencial de la velocidad del aire, vθ = ωr. Por otra 
parte, la aceleración de Coriolis es nula pues V//
ρρ
ω y no existen aceleraciones lineales o angulares del tubo de corriente. 
 
De la aplicación de la ecuación de la energía (7) en el dominio I, en el dominio II y en la totalidad del espacio, resultan 
2
2
I
I a4
3
2/)V(
u~
=
∆
=ζ
∞
 (20) 
para el coeficiente de pérdida de energía en la porción I de τ, 



+
−
−=
∆
=ζ
∞
2
2
2
2
II
II )a1(
)a1(
4
1a
2/)V(
u~  
(21) 
para el coeficiente de pérdida de energía en la porción II de τ, 



+
−
−−=ζ+ζ−=ζ 2
2
2
III )a1(
)a1(1a)a1())(a1(C  
(22) 
para el coeficiente de pérdida de potencia total, y 
)a1(
)a1(a4
AVq
WC
2
s
Ws
+
−
==
∞∞
&
 
(23) 
para el coeficiente potencia al eje de la turbina. 
 
Dado que las pérdidas representadas por la ecuación (22) y graficadas en la Figura 4 no son consideradas por el modelo de 
Betz, resulta un coeficiente de potencia al eje menor que el propuesto por dicha teoría. En la Figura 5 se grafica la variación 
del coeficiente de potencia al eje, ecuación (23), con el factor de reducción de velocidad en el disco. Se observa que el 
máximo valor del coeficiente apenas supera el valor 0,45 frente al (reconocidamente elevado) 0,60 que predice la teoría de 
Betz. Por otra parte, el presente modelo admite soluciones para a > 0,5 mientras que el modelo de Betz deja de ser válido en 
el rango 0,5 < a ≤ 1. Por su parte, la teoría de Dyment (1989) subestima el valor del máximo coeficiente de potencia obtenible 
(resultando éste aproximadamente 0,30). En cuanto a la comparación con resultados experimentales, Heijdra y Hendriks 
(2001) muestran que los máximos coeficientes de potencia observados en su análisis comparativo se encuentran entre 0,37 y 
0,44 (Table 2, p.28).  
 
Por otra parte, para definir el adecuado valor del parámetro λ∞ del modelo, se considera la relación presentada por Le 
Gourières (1982) proveniente de un estudio llevado a cabo por Hutter (p.54 op.cit.), donde se representa la variación 
funcional entre la velocidad específica λ∞ y la solidez del rotor, σ, que es la razón entre la superficie de las palas y el área 
barrida por las mismas. Dicha relación funcional (ver figura 54, p.55 op.cit) es aproximada como 
 
[ ]( ) ( )
∞
λ−≈σ 1010 log523,12%log  
 
(24) 
En el caso un disco, donde la razón de solidez es unitaria (σ = 100%), le corresponde por la relación (24) el valor de 
velocidad específica λ∞ = 1. Por lo cual, resulta que Ω = V∞/R es el valor de rotación más representativo de un disco 
actuador.
06.11 
  
 
 
Figura 4: Coeficiente de pérdida de potencia dentro del tubo de corriente τ, como función del parámetro a. 
 
 
 
Figura 5:Coeficiente de potencia en función del factor de disminución de velocidad en el disco, a. Modelo presente y de Betz 
 
CONCLUSIONES 
 
Se propone un modelo de disco actuador rotante que complementa a la representación del tubo de corriente que envuelve al 
aire que fluye a través del rotor de una turbina eólica de eje horizontal mediante la superposición potencial de una corriente 
uniforme de dirección del eje del disco y de una fuente puntual ubicada en coincidencia con el centro del disco.  
 
Se desarrollan las expresiones de los coeficientes característicos de un rotor eólico considerando el modelo descrito. Si bien 
el presente modelo es todavía muy simplificado, representa un avance frente al modelo de Betz y a aquellos modelos que, 
incorporando la rotación del disco, mantienen las hipótesis de la teoría de Betz en su desarrollo. 
 
REFERENCIAS 
 
Eggleston D.M. and Stoddard F.S. (1987) Wind Turbine Engineering Design. Van Nostrand Reinhold Co. Inc., New York. 
Dyment A. (1989) A modified form of the Betz’ wind turbine theory including losses. J. Fluids Engineering 111, 356-358. 
Heijdra J.J. and Hendriks H.B. (ed.) (2001) DOWEC Concept Study – Task 11: Selection of feasible concepts. ECN 
Windenergie, ECN-C--01-079. 
Le Gourières D. (1982) Energie éolienne. Deuxième édition, Eyrolles, Paris. 
Prado R.A. (1994) Reformulación de la teoría de la cantidad de movimiento aplicada a turbinas eólicas. Anales del VI 
Congreso Nacional de Ingeniería Mecánica, Santiago de Chile, pp.723-728. 
Prado R.A. (1995) Reformulation of the momentum theory applied to wind turbines. Journal of Wind Engineering and 
Industrial Aerodynamics 58, 277-292. 
Prado R.A. (1996) Modelo de disco actuador rotante para la simulación de una turbina eólica de eje horizontal. ASADES’96, 
Mar del Plata. 
Schlichting H. (1972) Teoría de la capa límite, Ediciones URMO, Bilbao, España. 
Shames I.H. (1992) Mecánica de Fluidos. McGraw-Hill, Colombia. 
 
ABSTRACT: A simplified model has been developed to predict the characteristic coefficients of horizontal-axis wind 
turbines. The model is based on the momentum theory and the inclusion of a rotating actuator disc. The power, axial force 
and torque coefficients are evaluated under the condition of disc rotation, in order to determine the energy losses in a model 
where a rotational inviscid wake is developed downstream of the rotating actuator disc that simulates the wind turbine. 
 
Keywords: inviscid model, actuator disc, rotation, horizontal-axis wind turbine. 
06.12 


Spanish

http://sedici.unlp.edu.ar/bitstream/handle/10915/81758/Documento_completo.pdf?sequence=1

Teoría del disco actuador rotante : Fuente puntual para turbinas eólicas de eje horizontal

Abstract

Similar works

Full text

Available Versions

Servicio de Difusión de la Creación Intelectual

SEDICI - Repositorio de la UNLP