Resumen: La resistencia a la cavitación es un requerimiento importante en el diseño de turbinas hidrocinéticas para aplicaciones marinas o fluviales dado que se ha demostrado que la cavitación puede contribuir al desgaste, corrosión, vibración y fatiga de los álabes de la turbina. La presencia de cavitación en los álabes puede conducir a la disminución del rendimiento de la turbina y la reducción de su vida útil. Por lo tanto, es fundamental incluir un estudio de cavitación en el análisis y desarrollo de los sistemas hidrocinéticos. En este trabajo, se presentan los elementos a tener en cuenta en un estudio de cavitación de los álabes de las turbinas hidrocinéticas. Como criterio para determinar la ocurrencia de cavitación, se presentó la comparación entre la distribución del coeficiente de presión (Cpre) sobre el perfil hidrodinámico Eppler 420 y el número de cavitación (σ). El Cpre fue calculado mediante simulación numérica usando el software Ansys Fluent. Los resultados mostraron que el perfil hidrodinámico Eppler 420 podría ser empleado para el diseño de los álabes de turbinas hidrocinéticas.
Palabras clave: cavitacióncavitación,número de cavitaciónnúmero de cavitación,inicio de la cavitacióninicio de la cavitación,turbina hidrocinéticaturbina hidrocinética,perfil hidrodinámicoperfil hidrodinámico,coeficiente de presióncoeficiente de presión,CFDCFD,Eppler 420Eppler 420,modelo de turbulencia k - ω SSTmodelo de turbulencia k - ω SST,coeficiente de sustentacióncoeficiente de sustentación.
Abstract: Resistance to cavitation is an important requirement in the design of hydrokinetic turbines for marine or river applications due to cavitation has been found to contribute to the turbine blade wear, corrosion, vibration, and fatigue. The presence of cavitation in the blades can lead to a decrease in the turbine performance and the reduction of its useful life. Therefore, it is crucial to include a cavitation study in the analysis and development of hydrokinetic systems. In this work, the elements to be considered in a cavitation study of the blades of a hydrokinetic turbine are presented. The comparison between the distribution of the pressure coefficient (Cpre) on the Eppler 420 hydrofoil and the number of cavitation (σ) was presented as a criterion to determine the cavitation occurrence. Cpre was calculated by numerical simulation by using the Ansys Fluent software. The results showed that the Eppler 420 hydrofoil could be used for the design of the blades of hydrokinetic turbines.
Keywords: cavitation, cavitation number, cavitation inception, hydrokinetic turbine, hydrofoil, pressure coefficient, CFD, Eppler 420, k -ω SST turbulence model, lift coefficient.
Artículos
Cavitación en perfiles hidrodinámicos para turbinas hidrocinéticas
Cavitation in hydrofoils for hydrokinetic turbines
Ana Barona-Mejía
Sara Gómez-Díaz
Jonathan Aguilar-Bedoya
Ainhoa Rubio-Clemente ainhoarubioclem@gmail.com
Edwin Chica-Arrieta
Recepción: 06 Agosto 2020
Aprobación: 11 Noviembre 2020
El plan energético de Colombia está alineado con los objetivos de desarrollo sostenible al contemplar la diversificación de la matriz energética con recursos energéticos de carácter renovable que garanticen el suministro oportuno, confiable, seguro y eficiente a toda la sociedad, y a su vez contribuyan a la reducción de emisiones de efecto invernadero y a la mitigación del cambio climático [1,2].
En Colombia, a pesar de que existen capacidades instaladas en tecnologías para el aprovechamiento de la energía hidráulica a filo de agua, eólica y solar, como fuentes de origen renovable [3-6], es importante emprender acciones para la apropiación y desarrollo de nuevas tecnologías no convencionales de generación de energía tales como el uso de equipos de generación a pequeña escala , incluyendo las turbinas hidráulicas, y/o el empleo de equipos no convencionales de generación como son las turbinas hidrocinéticas [7].
Hoy en día, la energía hidrocinética es reconocida como una fuente alterna de energía renovable no contaminante y que no requiere de grandes obras de infraestructura comparables con las que sí necesitan las convencionales centrales hidráulicas [8, 9]. Con el fin de aprovechar eficientemente la energía que se encuentra en la velocidad del flujo de una corriente de agua, es necesario desarrollar las tecnologías de transformación energética para las condiciones hidráulicas de los saltos y cauces colombianos [7, 8].
Durante el diseño es importante definir cuidadosamente los requisitos funcionales de la turbina hidrocinética, con miras a maximizar el rendimiento del sistema aumentando la captura de energía durante la vida útil, reduciendo las inestabilidades hidrodinámicas y fallas del sistema [8, 9]. Dentro de los requisitos se encuentran la potencia a generar, durabilidad, facilidad de fabricación, instalación, operación y mantenimiento del sistema entre otros. Un aspecto importante relacionado con la durabilidad es la resistencia de los álabes al fenómeno de la cavitación que puede provocar un deterioro del rendimiento de los álabes, vibraciones y aumento de la turbulencia [10, 11]. Este fenómeno cavitación puede cambiar la respuesta hidrodinámica del rotor ante el flujo debido a que disminuye la sustentación y aumenta el arrastre, por lo tanto, se debe evitar su aparición en los álabes de las turbinas hidrocinéticas [10,11].
El fenómeno de la cavitación se origina en la sección del álabe donde la presión disminuye por debajo de la presión de vapor del fluido. Ante este evento, el líquido se vaporiza instantáneamente, formando una cavidad de vapor que altera el flujo alrededor del álabe. La forma y el tamaño de la burbuja también varía en función del campo de presiones y velocidad alrededor del fluido. Cuando la cavidad de vapor implosiona, la presión sobre la superficie del álabe aumenta, lo que promueve la erosión de su superficie [12].
En este trabajo se presenta, un análisis del fenómeno de cavitación sobre un perfil hidrodinámico que puede ser empleado para el diseño de turbinas hidrocinéticas. El criterio utilizado para garantizar que no existe cavitación fue el del mínimo coeficiente de presión (Cpre). La distribución de Cpre fue determinada usando simulación numérica en el software Ansys Fluent.
El criterio más utilizado para minimizar o evitar la cavitación en turbinas hidráulicas consiste en comparar el número de cavitación (σ) o el factor de cavitación de Thoma's con el mínimo Cpre mediante la Ecuación (1) [13, 14].
donde el número de cavitación puede ser calculado empleando la Ecuación (2).
siendo patm la presión atmosférica; p, la densidad del agua; h, la distancia entre la superficie libre y la posición radial del rotor; pv , la presión de vapor; y V, la velocidad relativa a la sección del álabe. Es importante mencionar que el inicio de la cavitación ocurre en el perfil hidrodinámico cuando la presión local (p) en el perfil cae por debajo de pv , lo cual se puede predecir a partir de la distribución de presión [13, 14]. Con el fin de diseñar un rotor con un buen desempeño a la cavitación, a debería ser más alto que el Cpre , definido mediante la Ecuación (3). Adicionalmente, para garantizar un buen desempeño del perfil, los valores correspondientes al coeficiente de sustentación (CL ) y a la relación entre el coeficiente de sustentación (CL) y el de arrastre (CD) deben ser lo más altos posible. La cavitación se presenta cuando el mínimo Cpre es menor que σ [13, 14].
donde p es la presión local.
La Figura 1 muestra la presión estática en una sección del álabe de una turbina hidrocinética.
Molland y colaboradores en 2004 [14] describieron que las características de cavitación para un perfil hidrodinámico particular se pueden analizar mediante una envolvente de presión mínima (Figura 2a), o un área sin cavitación (Figura 2b) en función de σ. Dado que CL es función de la distribución de presión, los autores determinaron que, para un perfil hidrodinámico particular del álabe, el área libre de cavitación se puede representar como una envoltura CL de limitante a una base de σ . En la Figura 2b se muestra un esquema del área sin cavitación, junto con los tipos probables de cavitación. El ancho (vertical) del área es una medida de la tolerancia del perfil al funcionamiento sin cavitación; es decir, con un área más ancha, el perfil podrá tolerar una variación mucho mayor en el ángulo de ataque (α ) sin cavitación. El ancho y la forma del área libre de cavitación dependerán de las características del perfil, tales como el grosor, la curvatura, la forma general y la forma del borde de ataque. En este sentido, un aumento en el espesor del perfil hidrodinámico tiende a ensanchar el área libre de cavitación, mientras que un aumento en la inclinación del perfil conlleva a mover el mismo ancho y forma vertical del área hacia valores más altos de CL [14].
Los coeficientes CL y CD se pueden determinar empleando las ecuaciones (4) y (5), respectivamente [15].
donde S es la superficie del perfil hidrodinámico, Fx y Fy son las componentes de la fuerza hidrodinámica paralela y ortogonal a la dirección del flujo [15].
En la literatura, existen estudios sobre la cavitación de álabes de turbinas hidrocinéticas. A modo de ejemplo, Molland y colaboradores [14] analizaron cuatro perfiles hidrodinámicos de la serie NACA (4415, 6615, 63-215 y 63-815), con los cuales se mostró mediante el uso del software XFOIL que el contorno libre de cavitación depende de las características del perfil como la forma, el espesor y la cuerda [14].
Por su parte, Bajah y colaboradores [16] analizaron la variación del coeficiente de potencia (CP) y CL con respecto a la inmersión de la turbina. Para ello utilizaron un rotor de diámetro (D) igual a 800 mm, con 3 álabes con variaciones del perfil hidrodinámico de la referencia NACA 63-8xx (812, 815, 818, 821 y 824) y un túnel de cavitación en el que el caudal máximo era de 8 m3/s y la presión se encontraba en el rango de 0.2 a 1.2 atm. Los autores llegaron a la conclusión de que los valores de Cp obtenidos con inmersión poco profunda de la punta del álabe (la altura de la superficie libre del agua y la punta del álabe fue 0.19 D) fueron 10-15% más bajos que los resultados obtenidos con la inmersión profunda del álabe (la altura de la superficie libre del agua y la punta del álabe fue 0.55 D). Los resultados obtenidos de CL para el caso de inmersión poco profunda de la punta de los álabes fueron aproximadamente 5% más bajos que los encontrados en el caso de inmersión profunda de la punta del álabe. Los investigadores también concluyeron que la cavitación no aparece hasta que a se ha reducido aproximadamente a 0.9 con un ángulo en el cubo o hub de 25° [16].
Guo y colaboradores [17] utilizaron el programa Ansys CFX para simular una turbina con similares especificaciones a la utilizada en la investigación de Bajah con un ángulo del cubo de 20° y utilizando un modelo homogéneo de dos fases junto con el modelo de cavitación Rayleigh-Plesset, el cual valida que la cavitación ocurre en las puntas de los álabes cuando se encuentra a poca profundidad de inmersión en aguas poco profundas (la altura de la superficie libre del agua y la punta del álabe fue 0.19 D) [17].
Por su parte, Batten y colaboradores [18] presentaron un modelo basado en la teoría BEM (Blade Element Momentum) para el diseño hidrodinámico de una turbina hidrocinética. En el modelo, la predicción de la cavitación fue llevada a cabo mediante la comparación de a y Cpre. Los autores indicaron que la cavitación en álabes de turbinas hidrocinéticas se puede evitar con el uso de perfiles hidrodinámicos adecuados para el diseño de los álabes; por lo tanto, es importante la correcta selección del perfil durante el diseño de los álabes [18].
Goundar y colaboradores [19] también utilizaron la teoría BEM y modificaron el perfil hidrodinámico S1210 con un aumento de la cuerda y del grosor de un 20%, para obtener un buen rendimiento en las relaciones de velocidad de punta (TSR) sin tener cavitación. Los autores concluyeron que al incrementar el TSR aumentan las posibilidades de cavitación en las secciones de los álabes, por lo que es mejor mantener el TSR lo más bajo posible con el fin de obtener un mayor rendimiento sin que se produzca cavitación [19].
Por otro lado, en la literatura se muestra el diseño de una turbina hidrocinética tipo propela a pequeña escala desarrollado por Schleicher et al. Los autores emplearon el software Ansys Fluent para verificar la presencia de cavitación en los álabes. Los resultados del estudio mostraron que se presenta cavitación cuando hay corrientes de 5 m/s. Sin embargo, en el caso de la turbina objeto de estudio no se observó la existencia de corrientes con esta velocidad puesto que su diseño fue para corrientes inferiores a 2.2 m/s [20].
En 2015, Silva y colaboradores [12] discutieron el diseño del álabe de una turbina hidrocinética considerando la cavitación. En este estudio, se propuso una aproximación matemática para la prevención del fenómeno de cavitación. Cabe destacar que este enfoque usa el criterio del mínimo Cpre como el límite de la cavitación para el flujo sobre los álabes. Los autores referidos modificaron la velocidad relativa local con el objetivo de prevenir la ocurrencia de cavitación en cada sección del álabe. La modificación de la velocidad se logró a través del cambio de la longitud de la cuerda local del álabe sin influir de manera significativa en el Cp de la turbina objeto de estudio [12].
La optimización de la longitud de la cuerda y el ángulo de giro de cada sección trasversal del álabe de una turbina hidrocinética de eje horizontal fue desarrollada por Nigam y colaboradores [21] empleando la teoría BEM. Para lo lograr este propósito, los autores emplearon para el diseño del álabe dos perfiles hidrodinámicos diferentes, el perfil S832 y el perfil Eppler 817, logrando mantener la rigidez en el extremo del álabe y eliminando la cavitación. Los resultados de los autores mostraron que la combinación de perfiles propuesta fue adecuada para evitar la cavitación. Los autores establecieron que mediante el perfil Eppler 817, el cual fue empleado para generar el 25% de la longitud total del álabe desde la punta del álabe hacia la raíz o el hub (cubo), se puede reducir los esfuerzos y la cavitación. El 75% restante de la longitud del perfil fue diseñado con el perfil hidrodinámico S832 [21].
Silva y colaboradores [22] realizaron la optimización del diseño de una turbina hidrocinética de eje horizontal considerando la posibilidad de existencia de la cavitación mediante la teoría BEM. Los autores, como criterio de inspección de la cavitación, también utilizaron la comparación entre el mínimo Cpre y a. El modelo desarrollado por los autores mostró un buen comportamiento ante el fenómeno de la cavitación [22].
Por su parte, Muratoglu y colaboradores [23] compararon distintos perfiles hidrodinámicos (NACA, NREL y RISØ). Para los perfiles NACA, se analizaron las referencias NACA 4418, NACA 63-818 y NACA 63218. Los autores presentaron la distribución de los valores de Cpre de los perfiles estudiados. Los resultados mostraron que el perfil NACA 63-818 tiene un Cpre promedio mayor y una mejor relación CL/CD en comparación con los valores obtenidos con los otros dos perfiles NACA estudiados. Por lo que el perfil NACA 63818, a pesar de mostrar mejor rendimiento en términos de Cp , es más susceptible a tener cavitación. De igual manera, se compararon tres perfiles NREL (S822, S833 y S825) y cuatro perfiles RIS0 (RISØ A-12, RISØ A-15, RISØ A-18 y RISØ A-21). El perfil S833 mostró una mejor relación CL/CD en comparación con los perfiles S822 y S825. Por otro lado, el Cpre más alto se logró con el perfil S825. El valor correspondiente a CL/CD se extendió hasta 12 para un a igual a 12°. El valor de Cpre del perfil hidrodinámico S825 fue más alto que el obtenido con los otros dos perfiles. Los perfiles RIS0 usados generalmente en aerogeneradores tienen un borde de ataque pronunciado y bajos Cpre , lo cual provee un buen desempeño, así como baja cavitación en las turbinas hidrocinéticas. Finalmente, los autores concluyeron que los perfiles NACA 63 y RISØ-A proveen un buen desempeño para ser implementados en turbinas hidrocinéticas.
No obstante, según el análisis de los autores la mayoría de los perfiles NACA no son óptimos para ser usados en las turbinas debido a sus pobres características hidrodinámicas, las cuales conducen a la entrada en pérdida del perfil, comúnmente conocida como stall. Adicionalmente, este tipo de perfiles presenta baja eficiencia estructural cerca del centro, y un desempeño inestable con distintos números de Reynolds. Por ello, la mejor opción entre los perfiles estudiados por los autores son los perfiles RISØ-A [23].
Recientemente, Vaz y colaboradores [24]. propusieron un procedimiento de optimización de una turbina hidrocinética con difusor, dado que se ha demostrado que el difusor puede intensificar el fenómeno de la cavitación en la turbina [24]. La metodología presentada por los autores es una extensión de la metodología de optimización empleada por Vaz y Wood [25], en la cual el criterio del mínimo Cpre es incorporado en el coeficiente de empuje de la turbina para mantener la presión de la punta del álabe por encima de la presión de vapor del agua. La longitud de la cuerda de cada sección del álabe se determina mediante un análisis unidimensional donde se emplea una corrección para evitar la cavitación. Los resultados de la metodología de optimización propuesta por los autores fueron comprados con los estudios desarrollados por Silva y colaboradores [22, 29], en los cuales se evidenció que una mayor distribución de la cuerda puede evitar la cavitación. El modelo desarrollado por los autores referidos tiene un bajo costo computacional y una implementación numérica sencilla. Esta aproximación sugiere que la potencia de salida puede ser incrementada 42% para una velocidad del agua de 2.5 m/s, comparado con el modelo de optimización clásica de los autores Glauert y Durand detalla en la referencia [26].
El análisis de los estudios presentados muestra que al diseñar una turbina hidrocinética es importante seleccionar el perfil hidrodinámico más adecuado para evitar el inicio de la cavitación y maximizar la producción de energía. El criterio usado para determinar la ocurrencia de cavitación consiste en comparar el mínimo Cpre con a. Para este propósito, la herramienta más usada reportada en la literatura para determinar el valor de Cpre consiste en la realización de un análisis numérico usando un software de dinámica de fluido computacional (CFD).
En este trabajo, se realizaron simulaciones fluidodinámicas bidimensionales en el software Ansys Fluent con el fin de describir el comportamiento del flujo alrededor del perfil y determinar el valor de CPre. Las ecuaciones de Navier-Stoke promediadas de Reynolds (Reynolds-averaged Navier-Stoke-RANS equations) se consideraron como las ecuaciones gobernantes con el modelo de turbulencia k -ω SST. Cabe resaltar que el modelo de turbulencia referido se utiliza para simular el comportamiento de turbinas hidrocinéticas [27-31], debido a que ha demostrado tener un buen rendimiento para flujos complejos, incluyendo gradientes de presión adversos y separaciones de flujo, como ocurren en este tipo de turbinas. Este modelo ofrece una predicción mejorada de los gradientes de presión adversos en las regiones cercanas a la pared en comparación con los modelos k - ω y k - ε [32]. Las simulaciones CFD en este trabajo se realizaron utilizando el solucionador de estado estable (steady state solver) del software Ansys Fluent. Debe mencionarse que este enfoque ha sido utilizado por otros autores para el análisis de perfiles aerodinámicos e hidrodinámicos. En la presente investigación, se utilizó este esquema con el fin de reducir el tiempo computacional durante la simulación [27, 28, 33-35].
En la Figura 3 se presenta el dominio computacional empleado para el análisis fluidodinámico del perfil hidrodinámico para determinar los valores de Cpre, CL y CD. El dominio computacional tiene una topología en C y se extiende 10 veces la longitud de la cuerda (radio) aguas arriba y 20 veces la longitud de la cuerda aguas abajo.
En este estudio, se seleccionó el perfil Eppler 420 para conocer su resistencia a la cavitación. El perfil Eppler 420 fue seleccionado debido a que estudios previos mostraron que dicho perfil tiene una mejor relación Cl/Cd (47.77) entre 14 perfiles hidrodinámicos analizados [36, 37].
Por lo tanto, para el estudio se consideró un número de Reynolds (Re ) promedio de 789749.677 a lo largo del álabe y una velocidad relativa (Vrei ) promedio de 5.517 m/s. La determinación de Vrei puede ser consultada en detalle en la referencia [37]. A partir de la Ecuación (6) se determinó una longitud promedio de la cuerda (C), la cual fue igual a 0.1773 m para el perfil hidrodinámico utilizado para el análisis [37].
El dominio computacional fue mallado utilizando una malla no estructurada, la cual se generó para asegurar un y+≤ 1 colocando al menos 30 capas en la región de la capa límite. Se realizó un estudio de la independencia de la malla para asegurar la convergencia de la solución. Para ello, se seleccionaron 5 tamaños de la malla. En el Tabla 1 se muestra el estudio de independencia de la malla para un a del perfil igual a 0°. El CL y CD fueron elegidos parámetros de interés para el estudio de independencia de la malla y los resultados se compararon con los obtenidos con la malla 5. A la vista de estos resultados, la malla 3 fue seleccionada para los estudios numéricos.
Adicional al estudio de convergencia de la malla presentado, se calculó el índice de convergencia de la malla (GCI, Grid Convergence Index). A este procedimiento se le llama extrapolación de Richardson [38, 39, 40].
El GCI se puede obtener mediante la Ecuación (7) cuando se usa la solución de una malla fina o mediante la Ecuación (8) cuando se usa la solución de una malla gruesa [38, 39, 40].
Donde r es la relación entre el tamaño de la malla gruesa y el tamaño característico de la malla fina, es la estimación del error relativo de la solución (f) de una variable de control, que se obtiene como la diferencia entre la solución fina y la solución gruesa dividida entre la solución fina. Se multiplica por un factor de seguridad F que suele tomarse como 3 (a veces como 1.25, cuando se usan muchos tamaños de malla para estimar el orden de convergencia ). El orden de convergencia se estima con la Ecuación (9).
El GCI se suele expresar en tanto por ciento, y permite también estimar que tamaño de malla sería necesario para obtener una solución de cierta calidad, a partir de las mallas que ya se han probado [38, 39, 40]. El GCI para las variables CL y CD fueron iguales a 1.0011 y 0.9987. Un valor de GCI cercano a 1 es ideal para garantizar la independencia de la malla.
Las condiciones de frontera que fueron impuestas al dominio computacional son mostradas en la Tabla 2. Durante las simulaciones se varió a para obtener la mejor relación CL/CD.
A partir del análisis numérico en CFD, se obtuvo una relación CL/CD máxima de 39.05 (siendo CL y CD iguales a 1.42 y 0.036, respectivamente) para un a igual a 3°. En la Figura 5 se muestra la relación CL/CDvs. α con el fin de conocer el desempeño del perfil en un rango más amplio de α.
La Figura 6 muestra la distribución de Cpre para el perfil hidrodinámico Eppler 420. Con el fin de comparar los valores obtenidos de Cpre, se determina el número de cavitación empleando la Ecuación (2) considerando que p es igual a 997 kg/m3, la distancia entre la superficie del agua libre y el centro del rotor de la turbina, h = H y se estima en 1.5 m. Por su parte, g es igual a 9.81 m/s y Pv se corresponde con la presión de vapor a la temperatura de flujo; por lo tanto, se emplearon valores de Pv entre 872.579 Pa a 5 °C y 2339.256 Pa a 20 °C. Adicionalmente, el valor de V se asume igual a Vrei equivalente a 5.517 m/s. Al reemplazar los valores de las variables referidas anteriormente en la Ecuación (2), se obtiene que el valor más bajo de a es 4.235; sin embargo, como criterio de comparación se establece un valor de a igual a 4, valor reportado experimentalmente en la literatura [15].
Al analizar la distribución de |Cpre|, se puede concluir que ningún valor de este coeficiente sobre la superficie del perfil supera el valor de σ; por lo tanto, se puede concluir que el perfil no presentará cavitación, de modo que se podría emplear para el diseño de las alabes de la turbina hidrocinética. El valor más bajo de Cpre se presenta en la superficie superior del perfil cerca del borde de ataque. Si h es igual a 0 en la Ecuación 2 la punta de los álabes se encuentra sobre la superficie del agua, en este caso a es igual a 3.268. Sin embargo, en este escenario cuando el alabe se encuentre en la parte superior tampoco se presenta cavitación debido a que los valores de | Cpre | tampoco supera el valor de α.
En la Figura 7a y 7b se presentan las gráficas del contorno de presión y de velocidad alrededor del perfil Eppler 420 a un valor de α igual a 3°. Debido al flujo alrededor del perfil, se crea una región de alta presión y otra de baja presión en la parte inferior y superior del perfil, respectivamente. Esta diferencia de presión conduce a una diferencia de las velocidades en las regiones a través del principio de Bernoulli. Adicionalmente, la diferencia de presión referida hace que el agua en la región de alta presión genere una fuerza de elevación hacia arriba en el perfil hidrodinámico que contribuye a la generación del torque alrededor del eje del álabe.
En la Figura 7b se visualiza una buena adherencia de la capa límite en el extradós del perfil. Se muestra un desprendimiento de la capa límite del fluido muy cerca del borde de salida del perfil. Este desprendimiento produce alteraciones del campo de presiones que afectan ligeramente la sustentación y la resistencia hidrodinámica.
Estudios reportados en la literatura muestran que la cavitación en los álabes de turbinas hidrocinéticas puede ocasionar una disminución en su rendimiento y eventual falla de los alabes. Generalmente, la cavitación se puede presentar en las partes del álabe con altas velocidades, como ocurre en sus puntas. Por lo tanto, durante el diseño de los álabes de las turbinas hidrocinéticas es fundamental elegir un perfil hidrodinámico resistente a la cavitación. En este estudio, se analizó la resistencia a la cavitación del perfil hidrodinámico Eppler 420 mediante simulaciones numéricas utilizando el software Ansys Fluent. El criterio empleado para el análisis de cavitación del perfil consistió en comparar el mínimo Cpre con a. Los resultados indicaron que el perfil Eppler muestra una buena resistencia a la cavitación; mientras se mantengan las condiciones del análisis presentado para el álabe estudiado, por lo tanto, se recomienda su uso para el diseño de los álabes de una turbina hidrocinética.
Los autores agradecen al Programa Colombia Científica como fuente de financiación, en el marco de las convocatorias Ecosistema Científico (Contrato No. FP44842- 218-2018).
a Correos electrónicos: ana.barona@udea.edu.cobsara.gomez2@udea.edu.co,cjonathan.aguilar@udea.edu.coeedwin.chica@udea.edu.co
