Maestría en Ingeniería Mecánica · 2024
Optimización paramétrica de una turbina eólica de eje vertical de álabes rectos para el Cañón del Chicamocha
La microgeneración de energía a través de turbinas eólicas de eje vertical para bajas velocidades de viento en el departamento de Santander, Colombia, específicamente en el cañón del Chicamocha es viable. Sin embargo, existen diversos factores para el desarrollo de una turbina en la región, como lo son, la tipología de la turbina, solidez, relaciones de aspecto, perfil aerodinámico, longitud de cuerda, entre otros. En primer lugar, se desarrolla un algoritmo en el software MATLAB R2021a, en el cual se ingresan las variables físicas de la región tales como densidad, velocidad de viento y las variables de la turbina como el número de álabes, potencia deseada y ángulo de ataque deseado. El algoritmo tiene una serie de restricciones donde se tiene en cuenta la solidez de la turbina, la longitud de la cuerda y la relación de aspecto altura/diámetro. Una vez finalizado el algoritmo se procede a probar la optimización mediante la dinámica de fluidos computacional (CFD), con el fin de predecir el desempeño que tendrá la turbina en operación. Para esto, es necesario conocer las condiciones de contorno adecuadas para llevar a cabo la simulación, la discretización de la malla y modelos de turbulencia a utilizar. Con base en lo mencionado anteriormente, se realiza un análisis en ANSYS FLUENT de las optimizaciones generadas por el algoritmo con el cual se obtuvo el coeficiente de sustentación (Cl) y el coeficiente de arrastre (Cd), con los cuales fue posible calcular el coeficiente de potencia de cada turbina, con el cual se eligió la opción que tenga mejor rendimiento en la región de interés.
Texto completo 153 páginas con texto
Leer la tesis completa Ficha en el repositorio
Contenido
- introducciónp. 19
- 1. Objetivosp. 21
- 1.1 Objetivo Generalp. 21
- 1.2 Objetivos Específicosp. 21
- 2. Marco Teóricop. 22
- 2.1 Turbina Darrieus Tipo Rectop. 22
- 2.2 Turbina Savoniusp. 25
- 2.3 Turbina Helicoidalp. 27
- 2.4 Parámetros Geométricosp. 29
- 2.5 Dinámica De Fluidos Computacional (Cfd)p. 38
- 2.6 Proceso De Mallado O Dominio Computacionalp. 42
- 2.7 Condiciones De Contornop. 45
- 2.8 Modelos De Turbulenciap. 46
- 2.8.1 Modelo De Turbulencia Kω-Sstp. 47
- 2.8.2 Modelo Spalart-Allmarasp. 48
- 2.8.3 Modelo De Turbulencia K- Ɛp. 48
- 2.9 Solucionadoresp. 49
- 3. Metodologíap. 50
- 3.1 Desarrollo Del Algoritmo Desarrollado Para Generar Las Dimensiones De La Turbinap. 52
- 3.2 Geometría Del Perfil Aerodinámicop. 54
- 3.3 Discretización Espacial Y Condiciones De Contornop. 54
- 3.3.1 Zona Rotacional (“Rotate”)p. 55
- 3.3.2 Zona Fija O Ambientep. 55
- 3.3.3 “Inlet” O Entradap. 56
- 3.3.4 “Outlet” O Salidap. 56
- 3.3.5 “Up & Down”p. 56
- 3.3.6 “Aifoil”p. 56
- 3.3.7 “Interface Out” E “Interface In”p. 57
- 3.4 Modelo De Referencia Experimentalp. 58
- 3.5 Propiedades Del Fluidop. 60
- 3.6 Control Y Modelos De Soluciónp. 60
- 3.7 Modelo De Turbulenciap. 61
- 3.8 Residuos, Criterio De Convergencia, Time Steps Y Factores De Relajaciónp. 62
- 4. Resultadosp. 63
- 4.1 Algoritmo De Optimizaciónp. 63
- 4.2 Independencia De Mallap. 65
- 4.3 Contornos De Velocidadp. 73
- 4.4 Contornos De Presiónp. 74
- 4.5 Análisis De La Primera Propuesta De Optimización Ad=0.25 Y C=0.117mp. 75
- 4.6 Análisis De La Segunda Propuesta De Optimización Ad=0.25 Y C=0.173mp. 80
- 4.7 Análisis De La Tercera Propuesta De Optimización Ad=0.4 Y C=0.0929mp. 86
- 4.8 Análisis De La Cuarta Propuesta De Optimización Ad=0.4 Y C=0.137mp. 91
- 4.9 Análisis De La Quinta Propuesta De Optimización Ad=1 Y C=0.0587mp. 97
- 4.10 Análisis De La Sexta Propuesta De Optimización Ad=1 Y C=0.0864mp. 102
- 5. Análisis De Resultadosp. 107
- 5.1 Contornos De Velocidad De La Opción 2 (Ad=0.25; C=0.173)p. 114
- 5.2 Contornos De Presión De La Opción 2 (Ad=0.25; C=0.173)p. 120
- 5.3 Análisis De Y+ De La Opción 2 (Ad=0.25; C=0.173)p. 122
- 6. Conclusionesp. 123
- 7. Recomendacionesp. 124
- Referencias Bibliográficasp. 125
- Apéndicesp. 133
- Tabla 1. Velocidad de arranque automáticop. 34
- Tabla 2. Características de los diseños VAWT propuestosp. 35
- Tabla 3. Cp a varias velocidades de viento para diferentes modelos de turbinasp. 35
- Tabla 4. Situaciones simplificadas de las ecuaciones de gobiernop. 41
- Tabla 5. Características de los parámetros de calidadp. 44
- Tabla 6. Potencial eólico en el Cañón del Chicamochap. 52
- Tabla 7. Restricciones establecidas para desarrollo del algoritmop. 53
- Tabla 8. Modelo geométrico por utilizarp. 59
- Tabla 9. Dimensiones del modelo geométrico a utilizar en función del radio de la turbinap. 59
- Tabla 10. Propiedades del fluidop. 60
- Tabla 11. Controles de soluciónp. 61
- Tabla 12. Modelos de soluciónp. 61
- Tabla 13. Factores de relajación utilizadosp. 62
- Tabla 14. Validación del algoritmop. 63
- Tabla 15. Resultados obtenidos por el algoritmop. 64
- mallap. 65
- Tabla 17. Parámetros de calidad de la malla 1p. 65
- Tabla 18. Parámetros de calidad de la malla 2 y 3p. 66
- Tabla 19. % de Error entre los datos experimentales y el modelo 1 planteadop. 67
- Tabla 20. % de Error entre los datos experimentales y el modelo 2 planteadop. 68
- Tabla 21. % de Error entre los datos experimentales y el modelo 3 planteadop. 70
- Tabla 22. Parámetros de calidad de la malla opción 1 (AD=0.25 y c=0.117m)p. 76
- c=0.117m)p. 77
- Tabla 24. Parámetros de calidad de la malla opción 2 (AD=0.25 y c=0.173m)p. 81
- c=0.173m)p. 82
- Tabla 26. Parámetros de calidad de la malla opción 3 (AD=0.4 y c=0.0929m)p. 86
- c=0.0929m)p. 88
- Tabla 28. Parámetros de calidad de la malla opción 4 (AD=0.4 y c=0.137m)p. 91
- Tabla 30. Parámetros de calidad de la malla opción 5 (AD=1 y c=0.0587m)p. 97
- Tabla 32. Parámetros de calidad de la malla opción 6 (AD=1 y c=0.0864m)p. 102
- Tabla 34. Cp promedio de cada opción a V=4.5 m/s y α=0°p. 107
- Tabla 35. Cp promedio de cada opción a V=4.5 m/s y α=10°p. 109
- Tabla 36. Cp promedio de cada opción a V=7 m/s y α=0°p. 111
- Tabla 37. Cp promedio de cada opción a V=7 m/s y α=10°p. 112
- Figura 1. Turbina eólica Darrieusp. 23
- Figura 2. Cp en varias configuraciones de turbinasp. 24
- Figura 3. Forma final del perfil DU-06-W-200 comparado con el perfil NACA 0018p. 30
- Figura 4. Vista esquemática de la arquitectura de las turbinas Darrieusp. 30
- álabesp. 34
- Figura 6. Distribución de fuerzas y velocidades en el perfil del rotor Darrieusp. 37
- velocidad de punta de 0,5 y 1,5 a θp=0°p. 38
- Figura 8. Discretización unidimensional por volúmenes finitosp. 40
- Figura 9. Tipos de malla, malla no estructurada y malla estructuradap. 43
- Figura 10. Zona estacionaria y rotacionalp. 45
- Figura 11. Diagrama de la metodología utilizada en el presente trabajop. 51
- Figura 12. Diagrama de flujo del algoritmop. 53
- Figura 13. Modelo geométrico planteado. del estudio de independencia de mallap. 54
- Figura 14. Discretización espacial del estudio de independencia de mallap. 57
- Figura 15. Cp vs λ opción 1p. 67
- Figura 16. Cp vs λ opción 2p. 69
- Figura 17. Cp vs λ opción 3p. 70
- Figura 18. Cp vs λ de los modelos planteadosp. 71
- Figura 19. Coeficientes de sustentación (lift) y arrastre (drag) del modelo 2 planteadop. 72
- Figura 20. Contorno de velocidad del estudio de independencia de mallap. 73
- Figura 21. Contorno de presión del estudio de independencia de mallap. 75
- Figura 22. Mallado en la zona rotacional (AD=0.25 y c=0.117m)p. 76
- Figura 23. Mallado detallado del ataque del álabe (AD=0.25 y c=0.117m)p. 77
- Figura 24. Cp vs λ a 4.5 m/s (AD=0.25 y c=0.117m)p. 78
- Figura 25. Cp vs λ a 7 m/s (AD=0.25 y c=0.117m)p. 79
- c=0.117m)p. 77
- Figura 27. Mallado detallado del ataque del álabe (AD=0.25 y c=0.173m)p. 81
- Figura 28. Mallado detallado de la salida del álabe (AD=0.25 y c=0.173m)p. 82
- Figura 29. Cp vs λ a 4.5 m/s (AD=0.25 y c=0.173m)p. 83
- Figura 30. Cp vs λ a 7 m/s (AD=0.25 y c=0.173m)p. 84
- c=0.173m)p. 82
- Figura 32. Mallado detallado zona fija y rotacional (AD=0.4 y c=0.0929m)p. 86
- Figura 33. Mallado detallado del álabe (AD=0.4 y c=0.0929m)p. 87
- Figura 34. Cp vs λ a 4.5 m/s (AD=0.4 y c=0.0929m)p. 88
- Figura 35. Cp vs λ a 7 m/s (AD=0.4 y c=0.0929m)p. 89
- c=0.0929m)p. 88
- Figura 37. Mallado detallado zona fija y rotacional (AD=0.4 y c=0.137m)p. 91
- Figura 38. Mallado detallado del álabe (AD=0.4 y c=0.137m)p. 92
- Figura 39. Cp vs λ a 4.5 m/s (AD=0.4 y c=0.137m)p. 93
- Figura 40. Cp vs λ a 7 m/s (AD=0.4 y c=0.137m)p. 95
- c=0.137m)p. 96
- Figura 42. Mallado detallado zona fija y rotacional (AD=1 y c=0.0587m)p. 97
- Figura 43. Mallado detallado del álabe (AD=1 y c=0.0587m)p. 98
- Figura 44. Cp vs λ a 4.5 m/s (AD=1 y c=0.0587m)p. 99
- Figura 45. Cp vs λ a 7 m/s (AD=1 y c=0.0587m)p. 101
- c=0.0587m)p. 101
- Figura 47. Mallado detallado zona fija y rotacional (AD=1 y c=0.0864m)p. 103
- Figura 48. Mallado detallado del álabe (AD=1 y c=0.0864m)p. 103
- Figura 49. Cp vs λ a 4.5 m/s (AD=1 y c=0.0864m)p. 105
- Figura 50. Cp vs λ a 7 m/s (AD=1 y c=0.0864m)p. 105
- c=0.0864m)p. 106
- Figura 52. Comparativo de las opciones a V=4.5 m/s y α=0°p. 108
- Figura 53. Comparativo de las opciones a V=4.5 m/s y α=10°p. 109
- Figura 54. Comparativo de las opciones a V=7 m/s y α=0°p. 110
- Figura 55. Comparativo de las opciones a V=7 m/s y α=10°p. 113
- Figura 56. Contorno de velocidad de la optimización elegida V=4.5 m/s; α=0°p. 114
- Figura 57. Contorno de velocidad en el álabe V=4.5 m/s; α=0°p. 115
- Figura 58. Contorno de velocidad de la optimización elegida V=4.5 m/s; α=10°p. 116
- Figura 59. Contorno de velocidad en el álabe V=4.5 m/s; α=10°p. 116
- Figura 60. Contorno de velocidad de la optimización elegida V=7 m/s; α=0°p. 117
- Figura 61. Contorno de velocidad en el álabe V=7 m/s; α=0°p. 117
- Figura 62. Contorno de velocidad de la optimización elegida V=7 m/s; α=10°p. 118
- Figura 63. Contorno de velocidad en el álabe V=7 m/s; α=10°p. 119
- Figura 64. Contorno de presión de la optimización elegida V=4.5 m/s; α=0°p. 120
- Figura 65. Contorno de presión de la optimización elegida V=4.5 m/s; α=10°p. 121
- Figura 66. Contorno de presión de la optimización elegida V=7 m/s; α=0°p. 121
- Figura 67. Contorno de presión de la optimización elegida V=7 m/s; α=10°p. 122
- Figura 69. Línea de ruta de la optimización elegida V=4.5 m/sp. 123
- Figura 68. Línea de ruta de la optimización elegida V=7 m/sp. 123
- APÉNDICE A. COORDENADAS DE CADA PROPUESTA DE OPTIMIZACIÓNp. 133