Ley de Betz

Esquema del flujo de fluido a través de un actuador en forma de disco . Para un fluido de densidad constante, el área de la sección transversal varía inversamente con la velocidad.

En aerodinámica , la ley de Betz indica la potencia máxima que se puede extraer del viento, independientemente del diseño de una turbina eólica en flujo abierto. Fue publicada en 1919 por el físico alemán Albert Betz . ^[1] La ley se deriva de los principios de conservación de la masa y el momento de la corriente de aire que fluye a través de un " disco actuador " idealizado que extrae energía de la corriente de viento. Según la ley de Betz, ninguna turbina eólica de ningún mecanismo puede capturar más de 16/27 (59,3 %) de la energía cinética del viento. El factor 16/27 (0,593) se conoce como coeficiente de Betz . Las turbinas eólicas prácticas a escala de servicios públicos alcanzan en el pico el 75-80 % del límite de Betz. ^[2]^[3]

El límite de Betz se basa en un actuador de disco abierto . Si se utiliza un difusor para recoger el flujo de viento adicional y dirigirlo a través de la turbina, se puede extraer más energía, pero el límite sigue aplicándose a la sección transversal de toda la estructura.

Conceptos

La ley de Betz se aplica a todos los fluidos newtonianos , incluido el viento. Si toda la energía proveniente del movimiento del viento a través de una turbina se extrajera como energía útil, la velocidad del viento después caería a cero. Si el viento dejara de moverse a la salida de la turbina, entonces no podría entrar más viento fresco; estaría bloqueado. Para que el viento siga moviéndose a través de la turbina, tiene que haber algún movimiento de viento, por pequeño que sea, en el otro lado con una velocidad del viento mayor que cero. La ley de Betz muestra que a medida que el aire fluye a través de un área determinada, y a medida que la velocidad del viento disminuye debido a la pérdida de energía para la extracción de una turbina, el flujo de aire debe distribuirse a un área más amplia. Como resultado, la geometría limita la eficiencia máxima de cualquier turbina.

Descubrimientos independientes

El científico británico Frederick W. Lanchester derivó el mismo máximo en 1915. El líder de la escuela aerodinámica rusa, Nikolay Zhukowsky , también publicó el mismo resultado para una turbina eólica ideal en 1920, el mismo año que Betz. ^[4] Es, por tanto, un ejemplo de la ley de Stigler , que postula que ningún descubrimiento científico lleva el nombre de su descubridor real.

Prueba

El límite de Betz es la máxima energía posible que puede extraerse mediante un rotor infinitamente delgado de un fluido que fluye a una determinada velocidad. ^[5]

Para calcular la máxima eficiencia teórica de un rotor delgado (por ejemplo, de una turbina eólica ), se supone que se lo sustituye por un disco que absorbe energía del fluido que pasa a través de él. A cierta distancia detrás de este disco, el fluido que ha pasado a través del disco tiene una velocidad reducida, pero no nula. ^[5]

Suposiciones

El disco es ideal, sin buje (de modo que si se lo conceptualizara como un rotor, podría tener una cantidad infinita de aspas que no tienen resistencia). Se supone que cualquier no idealidad reduce la eficiencia.
Como modelo unidimensional, el flujo que entra y sale del disco es axial y todas las velocidades son uniformes transversalmente. Se trata de un análisis de volumen de control; el volumen de control debe contener todo el flujo entrante y saliente para poder utilizar las ecuaciones de conservación.
El flujo no es compresible, la densidad es constante y no hay transferencia de calor.
Se aplica una presión uniforme al disco. (En este modelo 1-D no existe dependencia radial de la presión).

Aplicación de la conservación de la masa (ecuación de continuidad)

Aplicando la conservación de masa al volumen de control, el caudal másico (la masa de fluido que fluye por unidad de tiempo) es

${\dot {m}}=\rho A_{1}v_{1}=\rho Sv=\rho A_{2}v_{2},$

donde v ₁ es la velocidad en la parte delantera del rotor, v ₂ es la velocidad aguas abajo del rotor, v es la velocidad en el dispositivo de potencia de fluido, ρ es la densidad del fluido, es el área de la turbina, y y son las áreas del fluido antes y después de la turbina (la entrada y la salida del volumen de control). ${\estilo de visualización S}$ $Estilo de visualización A_{1}$ $Estilo de visualización A_{2}$

La densidad multiplicada por el área y la velocidad deben ser iguales en cada una de las tres regiones: antes de la turbina, al pasar por la turbina y más allá de la turbina.

La fuerza que ejerce el rotor sobre el viento es la masa de aire multiplicada por su aceleración: ${\begin{aligned}F&=ma\\&=m{\frac {dv}{dt}}\\&={\dot {m}}\,\Delta v\\&=\rho Sv(v_{1}-v_{2}).\end{aligned}}$

Poder y trabajo

El trabajo incremental realizado por la fuerza puede escribirse

$dE=F\,dx,$

y la potencia (tasa de trabajo realizado) del viento es

$P={\frac {dE}{dt}}=F{\frac {dx}{dt}}=Fv.$

Sustituyendo la fuerza F calculada anteriormente en la ecuación de potencia se obtiene la potencia extraída del viento,

$P=\rho Sv^{2}(v_{1}-v_{2}).$

Sin embargo, la potencia se puede calcular de otra manera, utilizando la energía cinética. Aplicando la ecuación de conservación de la energía al volumen de control se obtiene

$P={\frac {\Delta E}{\Delta t}}={\tfrac {1}{2}}{\dot {m}}(v_{1}^{2}-v_{2}^{2}).$

Sustituyendo el caudal másico de la ecuación de continuidad se obtiene

$P={\tfrac {1}{2}}\rho Sv(v_{1}^{2}-v_{2}^{2}).$

Ambas expresiones para la potencia son válidas; una se derivó examinando el trabajo incremental y la otra mediante la conservación de la energía. Al igualar estas dos expresiones se obtiene

$P={\tfrac {1}{2}}\rho Sv(v_{1}^{2}-v_{2}^{2})=\rho Sv^{2}(v_{1} -v_ {2}).$

La densidad no puede ser cero para ningún v y S, por lo que

${\tfrac {1}{2}}(v_{1}^{2}-v_{2}^{2})={\tfrac {1}{2}}(v_{1}-v_ {2})(v_{1}+v_{2})=v(v_{1}-v_{2}),$

$v={\tfrac {1}{2}}(v_{1}+v_{2}).$

La velocidad constante del viento a través del rotor puede tomarse como el promedio de las velocidades ascendentes y descendentes. Podría decirse que esta es la etapa más contraintuitiva de la derivación de la ley de Betz. Es una consecuencia directa de la suposición de "flujo axial", que no permite ningún flujo de masa radial en la región del disco actuador. ^[6] Sin escape de masa y un diámetro constante en la región del actuador, la velocidad del aire no puede cambiar en la región de interacción. Por lo tanto, no se puede extraer energía más que en la parte delantera y trasera de la región de interacción, lo que fija la velocidad del aire del disco actuador como el promedio. (Eliminar esa restricción puede permitir un rendimiento superior al que permite la ley de Betz, pero también deben considerarse otros efectos radiales. ^[6] Este efecto de velocidad constante es distinto de la pérdida de energía cinética radial que también se ignora. ^[7] )

Ley de Betz y coeficiente de rendimiento

Volviendo a la expresión anterior para la potencia basada en la energía cinética:

${\begin{aligned}P&={\tfrac {1}{2}}{\dot {m}}(v_{1}^{2}-v_{2}^{2})\\&={\tfrac {1}{2}}\rho Sv(v_{1}^{2}-v_{2}^{2})\\&={\tfrac {1}{4}}\rho S(v_{1}+v_{2})(v_{1}^{2}-v_{2}^{2})\\&={\tfrac {1}{4}}\rho Sv_{1}^{3}\left(1+\left({\frac {v_{2}}{v_{1}}}\right)-\left({\frac {v_{2}}{v_{1}}}\right)^{2}-\left({\frac {v_{2}}{v_{1}}}\right)^{3}\right).\end{aligned}}$

El eje horizontal refleja la relación $v 2 / v 1$ , el eje vertical es el coeficiente de potencia ^[8] $C P$ .

Al derivar con respecto a para una velocidad de fluido dada $v$ $1$ y un área dada $S$ , se encuentra el valor máximo o mínimo para . El resultado es que alcanza el valor máximo cuando . $P$ ${\tfrac {v_{2}}{v_{1}}}$ $P$ $P$ ${\tfrac {v_{2}}{v_{1}}}={\tfrac {1}{3}}$

Sustituyendo este valor obtenemos:

$P_{\text{max}}={\tfrac {16}{27}}\cdot {\tfrac {1}{2}}\rho Sv_{1}^{3}.$

La potencia que se puede obtener de un cilindro de fluido con área de sección transversal $S$ y velocidad $v 1$ es

$P=C_{\text{P}}\cdot {\tfrac {1}{2}}\rho Sv_{1}^{3}.$

La potencia de referencia para el cálculo de la eficiencia de Betz es la potencia en un fluido en movimiento en un cilindro con área de sección transversal $S$ y velocidad $v 1$ :

$P_{\text{wind}}={\tfrac {1}{2}}\rho Sv_{1}^{3}.$

El coeficiente de potencia ^[9] $C P$ (= $P / P viento$ ) es la relación adimensional entre la potencia extraíble $P$ y la potencia cinética $P viento$ disponible en la corriente no distribuida. ^{[ cita requerida ]} Tiene un valor máximo $C P máx = 16/27 = 0,593$ (o 59,3%; sin embargo, los coeficientes de rendimiento se expresan normalmente como un decimal, no como un porcentaje). La expresión resultante es:

$C_{P}\left({\frac {v_{2}}{v_{1}}}\right)={\tfrac {1}{2}}\left(1+\left({\frac {v_{2}}{v_{1}}}\right)-\left({\frac {v_{2}}{v_{1}}}\right)^{2}-\left({\frac {v_{2}}{v_{1}}}\right)^{3}\right)$

Las grandes turbinas eólicas modernas alcanzan valores pico de $CP$ $en$ el rango de 0,45 a 0,50, ^[2]^{[ cita completa requerida ]} alrededor del 75-85% del máximo teóricamente posible. En alta velocidad del viento, cuando la turbina está operando a su potencia nominal, la turbina gira (inclina) sus aspas para reducir $CP$ y protegerse de daños. La potencia del viento aumenta en un factor de 8 de 12,5 a 25 m/s, por lo que $CP$ $debe$ caer en consecuencia, llegando a un mínimo de 0,06 para vientos de 25 m/s $.$

Entendiendo los resultados de Betz

La relación de velocidad entre el viento saliente y el entrante implica que el aire saliente tiene solo la energía cinética del aire entrante y que se extrajo la energía del aire entrante. Este es un cálculo correcto, pero solo considera el aire entrante que finalmente pasa por el rotor. ${\tfrac {v_{2}}{v_{1}}}={\tfrac {1}{3}}$ $({\tfrac {1}{3}})^{2}={\tfrac {1}{9}}$ ${\tfrac {8}{9}}$

El último paso para calcular la eficiencia de Betz $C p$ es dividir la potencia calculada extraída del flujo por una potencia de referencia. Como potencia de referencia, el análisis de Betz utiliza la potencia del aire que se mueve aguas arriba a $V 1$ a través del área de la sección transversal $S$ del rotor. Dado que en el límite de Betz, el rotor extrae de , o de la energía cinética entrante. $A_{1}={\tfrac {2}{3}}S$ ${\tfrac {8}{9}}$ ${\tfrac {2}{3}}$ ${\tfrac {16}{27}},$

Debido a que el área de la sección transversal del viento que fluye a través del rotor cambia, debe haber algún flujo de aire en las direcciones perpendiculares al eje del rotor. Cualquier energía cinética asociada con este flujo radial no tiene efecto en el cálculo porque este solo considera los estados inicial y final del aire en el sistema.

Límites superiores de las turbinas eólicas

Aunque a menudo se promociona (por ejemplo, ^[10]^[11] ) como el límite superior definitivo de la extracción de energía por cualquier turbina eólica posible, no lo es. A pesar del título engañoso de su artículo, ^[12] Betz (ni Lanchester) nunca hicieron una afirmación tan incondicional . ^[1] Cabe destacar que una turbina eólica que funciona con la máxima eficiencia de Betz tiene una estela de velocidad del viento distinta de cero. ^[2]^[1] Cualquier disco actuador colocado aguas abajo del primero extraerá energía adicional y, por lo tanto, el complejo de actuadores duales combinado excede el límite de Betz. ^[13]^[14] El segundo disco actuador podría estar, pero no necesariamente, en la zona de viento de campo lejano (línea de corriente paralela) para que se cumpla esta consideración. ^[6]

La razón de esta sorprendente excepción a una ley basada únicamente en leyes de conservación de energía y flujo se esconde en la suposición aparentemente modesta de uniformidad transversal del perfil axial del viento dentro de las líneas de corriente. Por ejemplo, la turbina eólica de doble actuador antes mencionada tiene, aguas abajo, un perfil transversal del viento que tiene dos velocidades distintas y, por lo tanto, no está limitado por los límites del disco actuador único. ^[6]

Matemáticamente, la derivación para un solo disco actuador incorpora implícitamente la suposición de que el viento no cambia de velocidad a medida que transita por el actuador "infinitamente delgado"; por el contrario, en el híbrido de actuador dual, el viento sí cambia de velocidad a medida que transita, invalidando el paso clave de la derivación que requiere una velocidad constante. Un solo actuador infinitamente delgado no puede cambiar la velocidad porque de lo contrario no conservaría el flujo, pero en el par híbrido, el flujo puede desprenderse (fuera de la sección transversal) entre los actuadores, lo que permite una velocidad de salida final diferente a la velocidad de entrada. ^[6]^[13]^[14]

Se han analizado turbinas eólicas físicas de rotor multicoaxial. ^[10] Aunque estas no superan el límite de Betz en la práctica, esto puede atribuirse al hecho de que los rotores no solo tienen pérdidas sino que también deben obedecer al momento angular ^[15] y a la teoría del momento del elemento de la pala que limita su eficiencia por debajo del límite de Betz. ^[7]

Las investigaciones modernas han sugerido que se puede lograr un límite superior más relajado cuando se eliminan las "suposiciones innecesarias" en la derivación de la ley de Betz. ^[16] ${\tfrac {2}{3}}$

Relevancia económica

La mayoría de las turbinas eólicas reales son aerodinámicamente "delgadas", lo que las hace aproximarse a los supuestos de la ley de Betz. En la medida en que una turbina eólica típica se aproxima a los supuestos de la ley de Betz, entonces el límite de Betz establece un límite superior aproximado a la energía anual que se puede extraer en un sitio. Incluso si un viento hipotético soplara de manera constante durante un año completo, cualquier turbina eólica que se aproxime bien mediante el modelo de disco del actuador no puede extraer más que el límite de Betz de la energía contenida en el viento de ese año.

En esencia, el aumento de la eficiencia económica del sistema se debe a una mayor producción por unidad, medida por metro cuadrado de exposición de las aspas. Se requiere un aumento de la eficiencia del sistema para reducir el costo de la producción de energía eléctrica. Los aumentos de la eficiencia pueden ser el resultado de la ingeniería de los dispositivos de captura de viento, como la configuración y la dinámica de las turbinas eólicas, que pueden aumentar la generación de energía de estos sistemas dentro del límite de Betz. Los aumentos de la eficiencia del sistema en la aplicación, transmisión o almacenamiento de energía también pueden contribuir a un menor costo de la energía por unidad.

Puntos de interés

Las suposiciones de la derivación de Betz ^[17] imponen algunas restricciones físicas sobre la naturaleza de las turbinas eólicas a las que se aplica (velocidad de entrada/salida idéntica, por ejemplo). Pero más allá de esas suposiciones, el límite de Betz no depende de la mecánica interna del sistema de extracción de viento, por lo tanto, $S$ puede tomar cualquier forma siempre que el flujo viaje desde la entrada al volumen de control hasta la salida, y el volumen de control tenga velocidades de entrada y salida uniformes. Cualquier efecto extraño solo puede disminuir el rendimiento del sistema (generalmente una turbina) ya que este análisis se idealizó para ignorar la fricción. Cualquier efecto no ideal restaría valor a la energía disponible en el fluido entrante, lo que reduciría la eficiencia general.

Algunos fabricantes e inventores han afirmado haber excedido el límite mediante el uso de toberas y otros dispositivos de desviación del viento, generalmente tergiversando el límite de Betz y calculando solo el área del rotor y no la entrada total de aire que contribuye a la energía eólica extraída del sistema.

El límite de Betz no tiene relevancia al calcular la eficiencia de una turbina en una aplicación móvil, como un vehículo eólico , ya que aquí la eficiencia podría acercarse teóricamente al 100% menos las pérdidas de las palas si el flujo de fluido a través del disco de la turbina (o equivalente) solo se retardara de manera imperceptible. Como esto requeriría una estructura infinitamente grande, los dispositivos prácticos rara vez alcanzan el 90% o más. La cantidad de energía extraída del flujo de fluido a altas eficiencias de la turbina es menor que el límite de Betz, que no es el mismo tipo de eficiencia. ^{[ cita requerida ]}

Desarrollo moderno

En 1934, H. Glauert derivó la expresión para la eficiencia de la turbina, cuando se tiene en cuenta el componente angular de la velocidad, aplicando un balance de energía a lo largo del plano del rotor. ^[15] Debido al modelo de Glauert, la eficiencia está por debajo del límite de Betz y se aproxima asintóticamente a este límite cuando la relación de velocidad de la punta tiende al infinito.

En 2001, Gorban , Gorlov y Silantyev introdujeron un modelo exactamente solucionable (GGS), que considera la distribución de presión no uniforme y el flujo curvilíneo a través del plano de la turbina (cuestiones no incluidas en el enfoque de Betz). ^[7] Utilizaron y modificaron el modelo de Kirchhoff , ^[18] que describe la estela turbulenta detrás del actuador como el flujo "degenerado" y utiliza la ecuación de Euler fuera del área degenerada. El modelo GGS predice que la eficiencia máxima se logra cuando el flujo a través de la turbina es aproximadamente el 61% del flujo total, lo que es muy similar al resultado de Betz de 2 ⁄ 3 para un flujo que resulta en eficiencia máxima, pero el GGS predijo que la eficiencia máxima en sí es mucho menor: 30,1%.

En 2008, se aplicaron cálculos viscosos basados en dinámica de fluidos computacional (CFD) al modelado de turbinas eólicas y demostraron una concordancia satisfactoria con el experimento. ^[19] La eficiencia óptima calculada está, típicamente, entre el límite de Betz y la solución GGS.

Referencias

^ abc Betz, A. (1966) Introducción a la teoría de las máquinas de flujo . (DG Randall, Trad.) Oxford: Pergamon Press.
^ abc "Familia E de Enercon, 330 kW a 7,5 MW, especificación de turbina eólica".
^ Tony Burton et al., (ed), Wind Energy Handbook , John Wiley and Sons 2001, ISBN 0471489972 , página 65.
^ Gijs A. M. van Kuik, El límite Lanchester–Betz–Joukowsky. Archivado el 9 de junio de 2011 en Wayback Machine , Wind Energ. 2007; 10:289–291.
^ ab Manwell, JF; McGowan, JG; Rogers, AL (febrero de 2012). La energía eólica explicada: teoría, diseño y aplicación . Chichester, West Sussex, Reino Unido: John Wiley & Sons Ltd., págs. 92-96. ISBN 9780470015001.
^ abcde Strauss, C.È.M. (octubre de 2021). "Superando la ley de Betz: un límite superior fundamental más grande para la recolección de energía eólica". arXiv : 2110.14612 [physics.flu-dyn].
^ abc Gorban', AN; Gorlov, AM; Silantyev, VM (diciembre de 2001). "Límites de la eficiencia de la turbina para el flujo de fluido libre" (PDF) . Revista de tecnología de recursos energéticos . 123 (4): 311–317. doi :10.1115/1.1414137.
^ "Coeficiente de potencia". Archivado desde el original el 31 de octubre de 2009.
^ "Asociación de la Industria Eólica Danesa". Archivado el 31 de octubre de 2009 en Wayback Machine.
^ ab Shinde, Sachin Manohar; Chaudhari, Mohit; Jeurkar, Tejas; Kadam, Sanket; Salunkhe, Kiran B. (2019). "Diseño y análisis de turbinas eólicas de rotor coaxial". Actas de la Conferencia internacional sobre fabricación inteligente y automatización . Apuntes de la conferencia sobre ingeniería mecánica. págs. 69-80. doi :10.1007/978-981-13-2490-1_7. ISBN 978-981-13-2489-5.
^ El límite de Betz y la máxima eficiencia de los aerogeneradores de eje horizontal
^ Betz, A. (2013). "El máximo de la explotación teóricamente posible del viento por medio de un motor eólico". Ingeniería eólica . 37 (4): 441–446. Bibcode :2013WiEng..37..441B. doi :10.1260/0309-524X.37.4.441. JSTOR 43857254.
^ ab Newman, BG (1983). "Teoría de discos actuadores para turbinas eólicas de eje vertical". Revista de ingeniería eólica y aerodinámica industrial . 15 (1–3): 347–355. Código Bibliográfico :1983JWEIA..15..347N. doi :10.1016/0167-6105(83)90204-0.
^ ab Newman, BG (1986). "Teoría de múltiples discos actuadores para turbinas eólicas". Revista de ingeniería eólica y aerodinámica industrial . 24 (3): 215–225. Código Bibliográfico :1986JWEIA..24..215N. doi :10.1016/0167-6105(86)90023-1.
^ ab White, FM, Mecánica de fluidos , 2.ª edición, 1988, McGraw-Hill, Singapur
^ Strauss, Charlie EM (20 de febrero de 2023). "Superando la ley de Betz: un límite superior fundamental más grande para la recolección de energía eólica". arXiv : 2110.14612 [physics.flu-dyn].
^ Rauh, A.; Seelert, W. (1984-01-01). "La eficiencia óptima de Betz para molinos de viento". Applied Energy . 17 (1): 15–23. Bibcode :1984ApEn...17...15R. doi :10.1016/0306-2619(84)90037-0. ISSN 0306-2619.
^ Milne-Thomson, LM (1960). Hidrodinámica teórica (4.ª ed.). Nueva York: Macmillan. pág. 632.
^ Hartwanger, D.; Horvat, A. (10 y 11 de junio de 2008). Modelado 3D de una turbina eólica mediante CFD (PDF) . Conferencia NAFEMS del Reino Unido de 2008: Simulación de ingeniería: uso eficaz y mejores prácticas. Cheltenham, Reino Unido. Archivado desde el original (PDF) el 7 de agosto de 2009.

Enlaces externos

Pierre Lecanu, Joel Breard, Dominique Mouazé. Límite de Betz aplicado a la teoría de turbinas eólicas de eje vertical